---
title: Three-Mode Soliton Josephson Junction
url: https://www.emergentmind.com/topics/three-mode-soliton-josephson-junction-tmsjj
type: topic
---

# Three-Mode Soliton Josephson Junction

The three-mode soliton Josephson junction (TMSJJ) is a bosonic Josephson system introduced as a three-mode extension of the previously proposed soliton Josephson junction bosonic model, with the specific aim of generating tripartite \(N00N\)-like states for quantum metrology in sensor networks [2507.15348]. In the formulation presently available, it consists of three weakly coupled bright solitons with attractive interactions, implemented most directly as atomic Bose–Einstein condensates in a symmetric three-well atomtronic circuit, and by analogy in Kerr nonlinear waveguides or microring resonators with weak inter-site coupling [2507.15348]. Its central result is an interaction-driven change in the ground state from a Gaussian-like coherent state to a superposition of entangled Fock states that rapidly approaches the three-mode \(N00N\) state, thereby linking nonlinear Josephson dynamics, multipartite entanglement, and multiparameter quantum metrology [2507.15348].

## 1. Conceptual placement and terminology

In the TMSJJ, “Josephson junction” is used in the bosonic sense: weakly coupled nonlinear modes exchange population and phase in close analogy to the bosonic Josephson junction. A closely related earlier optical construction is the weakly coupled nonlinear dielectric waveguide–surface-plasmon system, where the two-mode dynamics can be written in terms of a population imbalance \(Z\), a phase difference \(\phi\), and an amplitude-dependent coupling \(q(Z)\), and where self-trapped oscillatory states appear at nonzero fractional populations with zero and \(\pi\) time-averaged phase difference [1105.5937]. The TMSJJ extends that general bosonic-Josephson logic from two coupled nonlinear modes to three.

This usage is distinct from the superconducting multi-terminal Josephson-junction literature, where an \(N\)-terminal junction is a coherent conductor described by an \(N\)-terminal scattering matrix, with \(N-1\) independent superconducting phase differences and a critical-current manifold that forms the boundary of a \((N-1)\)-dimensional region of simultaneously allowed supercurrents [1812.06017]. A plausible source of confusion is that both literatures employ the language of phases, modes, and Josephson couplings, but the TMSJJ introduced in current work is a bright-soliton three-mode bosonic model rather than a three-terminal superconducting weak link.

## 2. Physical realization and microscopic model

The TMSJJ is formulated for three weakly coupled bright solitons with attractive interactions. The main implementation discussed is a symmetric three-well atomtronic circuit populated by attractive atoms such as \(^{7}\mathrm{Li}\), with one bright soliton per well; the same model is stated to apply to bright optical solitons in Kerr nonlinear waveguides or microring resonators [2507.15348]. The geometry has full cyclic symmetry and equal, each-to-each tunnel coupling.

The starting many-body Hamiltonian is
\[
\hat{H} =
\sum\limits_{j=1}^3 \hat{a}_j^\dag \left(
-\frac{1}{2}\frac{\partial^2}{\partial x^2}
- \frac{u}{2}\hat{a}_j^\dag\hat{a}_j
\right)\hat{a}_j
- \kappa \sum\limits_{j=1}^3\sum_{i\neq j}\hat{a}^\dag_{i}\hat{a}_j,
\]
where \(\hat a_j \equiv \hat a_j(x)\) are bosonic field operators, \([\hat a_i(x),\hat a_j^\dag(x')]=\delta_{ij}\delta(x-x')\), \(u=2\pi|a_{sc}|/r_0\) is the nonlinear interaction parameter, \(a_{sc}<0\) is the scattering length, \(r_0=\sqrt{\hbar/(M\omega_0)}\) is the transverse trap size, and \(\kappa\) is the equal tunnel coupling [2507.15348].

A fixed-\(N\) Hartree variational state is taken as
\[
\left|\Psi_N\right\rangle =
\frac{1}{\sqrt{N!}}
\Bigg[
\int_{-\infty}^\infty
\Big(
\psi_1(x)\hat{a}_1^\dag(x)
+\psi_2(x)\hat{a}_2^\dag(x)
+\psi_3(x)\hat{a}_3^\dag(x)
\Big)\,dx
\Bigg]^N
\left|0\right\rangle,
\]
with
\[
\sum_{j=1}^3 \int_{-\infty}^\infty |\psi_j(x)|^2 dx
=
\sum_{j=1}^3 \frac{N_j}{N}
\equiv
\sum_{j=1}^3 n_j
=1.
\]

| Quantity | Definition | Role |
|---|---|---|
| \(u\) | \(u=2\pi|a_{sc}|/r_0\) | nonlinear interaction parameter |
| \(\kappa\) | equal coupling on all links | tunnel coupling |
| \(\Lambda\) | \(\Lambda=\dfrac{u^2(N-1)^2}{16\kappa}\) | interaction–tunneling control parameter |

The presence of equal pairwise tunneling is essential: the TMSJJ is not presented as a chain with nearest-neighbor-only coupling, but as a symmetric three-mode junction.

## 3. Semiclassical reduction and nonlinear Josephson structure

In the uncoupled limit \(\kappa=0\), each mode is a bright soliton,
\[
\psi_j(x)=
n_j \frac{\sqrt{u(N-1)}}{2}\,
\mathrm{sech}\!\left[
\frac{u(N-1)}{2}n_j x
\right] e^{i\theta_j},
\]
with
\[
\theta_j=\frac{u^2(N-1)^2 n_j^2}{8}t.
\]
When weak coupling is restored, the populations \(n_j(t)\) and phases \(\theta_j(t)\) become dynamical variables [2507.15348].

Substitution of the soliton ansatz into the Hamiltonian yields the effective semiclassical energy per particle
\[
H_{\text{eff}}
=
-2\kappa\sum_j\left(
\frac{\Lambda}{3}n_j^3
+
\frac{1}{4}\sum_{i\neq j} I_{ij}\cos(\theta_j-\theta_i)
\right),
\]
where
\[
I_{ij}\equiv n_{ij}(1-z_{ij}^2)\,(1-0.21z_{ij}^2),\qquad
n_{ij}=n_j+n_i,\qquad
z_{ij}=\frac{n_j-n_i}{n_{ij}},
\]
and
\[
\Lambda=\frac{u^2(N-1)^2}{16\kappa}.
\]
This \(\Lambda\) is the principal control parameter of the model. Increasing \(u\) or \(N\), or decreasing \(\kappa\), drives the system toward the strongly nonlinear regime.

The semiclassical stationary-phase analysis identifies families of phase-structured three-mode \(N00N\)-like states. For extremely imbalanced populations \(n_1\approx1\), \(n_2\approx n_3\approx0\), the stationary phases satisfy
\[
\cos[\Theta_{12}]=\cos[\Theta_{31}]
=
\frac{\Lambda - 0.5\cos[\Theta_{23}]}{1.58}.
\]
Two special cases are distinguished. In the out-of-phase case,
\[
\Theta_{12}=\Theta_{31}\equiv \Theta_-,\qquad
\Theta_{23}=-2\Theta_-,
\]
with
\[
\cos[\Theta_-]=\sqrt{1.124+\Lambda}-0.79,\qquad \Lambda\le 2.08.
\]
In the in-phase case,
\[
\Theta_{12}=-\Theta_{31}\equiv \Theta_+,\qquad
\Theta_{23}=0,
\]
with
\[
\cos[\Theta_+]=\frac{\Lambda-0.5}{1.58},\qquad \Lambda\le 2.08.
\]
These stationary solutions underlie the later fully quantum \(N00N\)-like states and show that the TMSJJ possesses a genuinely three-phase nonlinear Josephson structure rather than a trivial extension of a two-mode dimer [2507.15348].

## 4. Quantum Hamiltonian, spectrum, and phase transition

The fully quantum TMSJJ is obtained by quantizing the effective semiclassical model in terms of number operators \(\hat N_i=\hat a_i^\dag \hat a_i\), the phase-number representation \(\hat a_i=\sqrt{\hat N_i}e^{i\hat\theta_i}\), and the imbalance operator
\[
\hat z_{ij}
=
\frac{\hat a_j^\dag\hat a_j-\hat a_i^\dag\hat a_i}
{\hat a_j^\dag\hat a_j+\hat a_i^\dag\hat a_i}.
\]
The tunneling prefactor is dressed by the series
\[
\sqrt{1-\hat z_{ij}^2}
=
\sum_{k=0}^\infty (-1)^k C_{0.5}^k \hat z_{ij}^{2k},
\qquad
C_{0.5}^k
=
\frac{1}{k!}\prod_{l=0}^{k-1}(0.5-l).
\]
The resulting Hamiltonian contains an effective cubic on-site term,
\[
-\frac{\Lambda}{3N^3}\sum_i (\hat a_i^\dag \hat a_i)^3,
\]
and a tunneling term dressed by imbalance corrections \(1-0.21\hat z_{ij}^2\) and the series in \(\hat z_{ij}^{2k}\) [2507.15348].

At fixed total particle number \(N\), the state is expanded as
\[
\left|\Psi(\tau)\right\rangle
=
\sum_{N_1=0}^N \sum_{N_2=0}^{N-N_1}
A_{N_1,N_2}(\tau)\,
|N_1,N_2,N_3\rangle,
\qquad
N_3=N-N_1-N_2,
\]
with \(\tau=2\kappa t\). The amplitudes obey coupled occupation-space hopping equations with diagonal coefficients
\[
\alpha_{N_i,N_j}
=
-\frac{\Lambda}{3}\,
\frac{N_i^3+N_j^3+(N-N_i-N_j)^3}{N^3},
\]
and tunneling coefficients \(\beta_{N_i,N_j}\) determined by the imbalance-dressed hopping matrix elements [2507.15348].

Numerical diagonalization for \(N=20\) shows a sharp quantum phase transition at
\[
\Lambda_{cr}\approx 3.30272.
\]
For \(\Lambda=0\), the ground state is Gaussian-like and centered near equal populations \(N_1\approx N_2\approx N_3\approx N/3\). At \(\Lambda=\Lambda_{cr}\), the ground state exhibits coexistence of a broad Gaussian-like central structure with substantial weight at the edges \((N,0,0)\), \((0,N,0)\), and \((0,0,N)\). For \(\Lambda=3.305>\Lambda_{cr}\), the Gaussian part collapses and the weight is sharply concentrated at those three extreme Fock configurations, so the ground state is essentially the three-mode \(N00N\) state [2507.15348].

The same analysis gives an approximate ground-state energy below the transition,
\[
E/\kappa N \approx -1.911 - 0.008\Lambda,\qquad \Lambda<\Lambda_{cr},
\]
while the symmetric three-mode \(N00N\) state has
\[
E = -\kappa N\,\frac{2\Lambda}{3}.
\]
The crossing of these behaviors is the spectral signature of the coherent-to-\(N00N\) ground-state transition.

## 5. Three-mode \(N00N\) states and metrological function

The symmetric three-mode \(N00N\) state is
\[
|N00N\rangle
=
\frac{1}{\sqrt{3}}
\left(
|N,0,0\rangle
+
|0,N,0\rangle
+
|0,0,N\rangle
\right),
\]
and more general phase-structured variants are
\[
|N00N\rangle_{\pm}
=
\frac{1}{\sqrt{3}}
\Big(
|N,0,0\rangle
+
e^{iN\Theta_\pm}|0,N,0\rangle
+
e^{\pm iN\Theta_\pm}|0,0,N\rangle
\Big),
\]
with \(\Theta_\pm\) given by the stationary-phase conditions of the semiclassical analysis [2507.15348]. These states are highly mode-entangled and serve as probe states for multiparameter sensing.

The metrological analysis introduces the General Heisenberg Limit
\[
\sigma_{GHL}^{(k)}\ge \frac{1}{N^k},
\]
with \(k=1\) for linear metrology and \(k=3\) for nonlinear metrology with solitons. For multiparameter estimation of \(\boldsymbol{\chi}=\{\chi_j\}\), the overall accuracy is
\[
\sigma_{\boldsymbol{\chi}}
\equiv
\left(\sum_{j=1}^d \sigma_{\chi_j}^2\right)^{1/2},
\]
bounded by the multiparameter quantum Cramér–Rao relation
\[
\sigma_{\boldsymbol{\chi}}
\ge
\left(\mathrm{Tr}(\hat F^{-1})\right)^{1/2},
\]
where \(\hat F\) is the quantum Fisher-information matrix [2507.15348].

For the balanced multipartite \(N00N\) input, the bound becomes
\[
\sigma_{\boldsymbol{\chi}}
\ge
\frac{1}{N^k}\sqrt{\frac{d(d+1)}{2}},
\]
and for the optimized unbalanced choice discussed in the same work,
\[
\sigma_{\boldsymbol{\chi}}
\ge
\frac{1}{N^k}\frac{\sqrt d(\sqrt d+1)}{2}
\equiv \sigma_{OS}^{(k)}.
\]
In the three-mode setting used for two-parameter estimation, \(d=2\), so the balanced limit is \(\sqrt{3}\,\sigma_{GHL}^{(k)}\), while the optimized-state limit is approximately \(\sqrt{2.914}\,\sigma_{GHL}^{(k)}\) [2507.15348].

The TMSJJ is used as a source of these three-mode probe states in a sensor network with two sensing arms and one reference arm. Phase encoding is modeled by
\[
\hat U_{PS}
=
\exp\!\left[
i\chi_1(\hat a_1^\dag \hat a_1)^k
+
i\chi_2(\hat a_2^\dag \hat a_2)^k
\right].
\]
For the Gaussian three-mode coherent-like state, the benchmark limits are
\[
\sigma_{SQL}=\sqrt{\frac{3}{N}},\qquad
\sigma_{NQL}\approx \sqrt{\frac{27}{N^5}},
\]
for \(k=1\) and \(k=3\), respectively [2507.15348].

The same work models losses by fictitious beam splitters with transmittance \(\eta\), leading to interferometric limits
\[
\sigma_{SIL}=\sqrt{\frac{3}{\eta N}},\qquad
\sigma_{NIL}\approx \sqrt{\frac{27}{\eta N^5}}.
\]
The reported numerical result is that, for \(\Lambda>\Lambda_{cr}\), the TMSJJ ground state approaches the optimized-state limit and remains close to the GHL in the lossless case, while under weak losses it still surpasses the classical interferometric limits [2507.15348]. For single-parameter estimation of
\[
\chi\equiv \frac{\Lambda}{N^2}=\frac{u^2}{16\kappa},
\]
the phase-structured \(N00N\) states yield
\[
\sigma_{\chi_+}=\frac{1.58}{N^3},\qquad
\sigma_{\chi_-}=\frac{1.22}{N^3},
\]
so both achieve \(1/N^3\) scaling, with the out-of-phase state slightly better because of nonlinear soliton phase counter-accumulation [2507.15348].

## 6. Relation to adjacent Josephson and soliton literatures

The TMSJJ sits at the intersection of several previously separate lines of work. One antecedent is the optical two-mode “new type of Josephson junction” formed by a nonlinear dielectric waveguide and a surface plasmon, where the coupling \(q(Z)=\exp[-kd\sqrt{2\Lambda(1+Z)}]\) depends dynamically on the population imbalance and gives rise to self-trapped oscillatory states [1105.5937]. Another is the long \(0\)-\(\pi\) Josephson-junction literature, in which phase-shift defects support localized modes, double-well mode tunneling, and bright and dark solitons, with multi-mode approximations derived from the underlying sine-Gordon dynamics [1009.5922]. These works are physically distinct from the TMSJJ, but they show that “soliton Josephson junction” has a broader prehistory involving nonlinear localized modes rather than only superconducting weak links.

A separate, and non-equivalent, branch is the superconducting multi-terminal literature. In top-down InAs/Al devices, the multi-terminal Josephson effect is described in terms of phase-dependent Andreev bound states, scattering matrices, and critical-current manifolds in current space [1812.06017]. Gate-defined quantum point contacts have since provided selective control of conductance modes in three-terminal Josephson devices and access to the single-mode regime coexistent with superconducting coupling [2201.01373]. Gate-tunable mesoscopic three-terminal InAs/Al junctions have also been characterized by differential-resistance maps and modeled successfully by a resistively and capacitively shunted-junction network at zero field [1905.11730]. PbTe nanowire three-terminal junctions add evidence for few-mode transport and \(\pi\)-shifted Cooper quartets, encoded phenomenologically by three-phase Josephson-energy terms beyond pairwise couplings [2312.17703]. A different microscopic route to nontrivial three-terminal physics appears in a triple-terminal junction with parallel-coupled double quantum dots, where a \(4\pi\)-periodic Josephson effect and spinful many-particle Majorana modes were reported in theory [1507.07164].

This broader record suggests that “three-mode” presently has at least two distinct meanings in Josephson research. In the TMSJJ proper, the modes are bright-soliton bosonic modes in a symmetric three-well or three-resonator setting [2507.15348]. In superconducting electronic junctions, “few-mode” usually denotes a small number of conductance channels or Andreev modes in a three-terminal weak link [2201.01373]. The terminological overlap is substantial, but the physical models, Hamiltonians, and intended applications differ.

Two further points delimit the present state of the TMSJJ. First, the explicit TMSJJ spectrum and transition are currently presented numerically for \(N=20\), with larger-\(N\) behavior argued qualitatively rather than solved exactly [2507.15348]. Second, the loss analysis treats independent Markovian particle losses through fictitious beam splitters and does not include more elaborate dephasing or correlated-noise models [2507.15348]. Even with those limitations, the TMSJJ is already defined sharply enough to constitute a specific research object: a symmetric, fully coupled three-bright-soliton Josephson trimer whose interaction-driven ground-state transition furnishes multipartite \(N00N\)-like resources for linear and nonlinear quantum metrology.

Source: https://www.emergentmind.com/topics/three-mode-soliton-josephson-junction-tmsjj