---
title: Two-Site Bose-Hubbard Model
url: https://www.emergentmind.com/topics/two-site-bose-hubbard-model
type: topic
---

# Two-Site Bose-Hubbard Model

The two-site Bose-Hubbard model is the minimal interacting bosonic lattice model with two modes, usually interpreted as a bosonic double well, a bosonic Josephson junction, or a two-mode condensate. In its standard form it combines intersite tunnelling with onsite interaction, while conserving the total particle number. Because the fixed-\(N\) sector is finite dimensional and admits both exact and controlled approximate treatments, the model has become a canonical setting for analyzing Josephson oscillations, self-trapping, fragmentation, Schrödinger-cat formation, current algebras, Bethe-ansatz integrability, and semiclassical collapse-and-revival dynamics [1101.4762] [1505.06793] [1409.5115].

## 1. Hamiltonian structure and Fock-space formulation

A standard parameterization for \(N\) bosons in two weakly coupled modes is
\[
\hat{H} = -\hbar J\left(\hat{a}_1^\dagger \hat{a}_2+\hat{a}_2^\dagger \hat{a}_1\right) + \frac{\hbar U}{2} \left( \hat{a}_1^{\dagger 2}\hat{a}_1^2 + \hat{a}_2^{\dagger 2}\hat{a}_2^2 \right),
\]
where \(J>0\) is the hopping or tunnelling amplitude and \(U\) is the onsite interaction strength, with \(U>0\) corresponding to repulsive interactions. A biased form,
\[
H=-\varepsilon \left(a_L^\dagger a_L - a_R^\dagger a_R\right) - J \left(a_L^\dagger a_R + a_L a_R^\dagger\right) + \frac{U}{2}\left( a_L^\dagger a_L^\dagger a_L a_L + a_R^\dagger a_R^\dagger a_R a_R \right),
\]
is used to select one branch of the attractive problem, while the canonical Josephson Hamiltonian is often written as
\[
\hat H=\frac{K}{8}(\hat N_1-\hat N_2)^2-\frac{\Delta\mu}{2}(\hat N_1-\hat N_2) -\frac{\mathcal E_J}{2}\left(\hat a_1^\dagger \hat a_2+\hat a_2^\dagger \hat a_1\right).
\]
More general two-site Hamiltonians allow \(K_{ij}\hat N_i\hat N_j\), on-well energies, chemical potentials, and asymmetric tunnelling amplitudes \(\Omega_{ij}\) [1101.4762] [1108.4851] [1505.06793] [1511.05026] [1610.07111].

The conserved total number \(\hat N=\hat N_1+\hat N_2\) decomposes the Hilbert space into fixed-particle-number sectors. In the occupation basis \(|l,N-l\rangle\), with \(l\) bosons in the left well and \(N-l\) in the right well, the state amplitudes \(c_l(t)\) obey a tridiagonal system,
\[
i\frac{dc_l}{dt} = -\left(\kappa_l c_{l+1}+\kappa_{l-1}c_{l-1}\right) +V_l c_l,
\qquad l=0,1,\dots,N,
\]
with
\[
\kappa_l = J\sqrt{(l+1)(N-l)}, \qquad
V_l=\frac{U}{2}\left[l^2+(N-l)^2-N\right].
\]
This representation is structurally decisive: the two-site problem becomes a one-dimensional tight-binding chain in Fock space, with nonuniform nearest-neighbour couplings and site-dependent diagonal energies. In the fixed-\(N\) sector the model therefore interpolates between a tunnelling-dominated delocalized regime and interaction-dominated occupation-number localization [1101.4762].

For attractive interactions, one frequently introduces the dimensionless control parameter
\[
\Lambda \equiv \frac{NU}{J},
\]
which organizes the crossover from a nearly coherent condensate to fragmented and cat-like regimes. In the large-\(N\) classical description the bifurcation occurs near \(|\Lambda|=2\), but for finite \(N\) and weak bias the crossover is smooth and shifted [1108.4851].

## 2. Algebraic structure, spin mappings, and current variables

The two-site model admits an \(SU(2)\)-type reformulation in terms of collective currents. A particularly useful choice is
\[
\hat{\mathcal I}=\frac12(\hat N_1-\hat N_2), \qquad
\hat{\mathcal J}=\frac{1}{2i}\left(\hat a_1^\dagger \hat a_2-\hat a_2^\dagger \hat a_1\right), \qquad
\hat{\mathcal T}=\frac12\left(\hat a_1^\dagger \hat a_2+\hat a_2^\dagger \hat a_1\right),
\]
respectively the imbalance current, the Josephson tunnelling current, and the coherent correlation tunnelling current. They satisfy the closed algebra
\[
[\hat{\mathcal T},\hat{\mathcal J}]=+i\hat{\mathcal I},\qquad
[\hat{\mathcal T},\hat{\mathcal I}]=-i\hat{\mathcal J},\qquad
[\hat{\mathcal J},\hat{\mathcal I}]=+i\hat{\mathcal T}.
\]
With \(\hat L_x=\hbar\hat{\mathcal T}\), \(\hat L_y=\hbar\hat{\mathcal J}\), and \(\hat L_z=\hbar\hat{\mathcal I}\), this is the standard angular-momentum algebra [1505.06793] [1511.05026].

The corresponding Casimirs are the conserved total number and the quadratic invariant
\[
\hat{\mathfrak C}_2=\hat{\mathcal T}^2+\hat{\mathcal I}^2+\hat{\mathcal J}^2
=\frac{\hat N}{2}\left(\frac{\hat N}{2}+1\right).
\]
This gives the model a Bloch-sphere geometry: the dynamics in a fixed-\(N\) sector is constrained to an effective spin of length \(N/2\). In semiclassical treatments this becomes explicit because the rescaled spin operators satisfy commutators proportional to \(1/N\), so \(\hbar_{\mathrm{eff}}=1/N\) and the limit \(N\to\infty\) is the classical limit [1505.06793] [1409.5115].

A closely related identification appears in the hopping-only dimer. Restricting the Hamiltonian
\[
\mathcal H=\hat a_2^\dagger \hat a_3+\hat a_3^\dagger \hat a_2
\]
to the \(k\)-particle sector \(\mathbb H^{\odot k}\), one obtains exactly the spin projection operator \(S_x\) for spin \(s=k/2\). In the symmetric and antisymmetric modes
\[
\hat c=\frac{\hat a_2+\hat a_3}{\sqrt 2},\qquad
\hat d=\frac{\hat a_2-\hat a_3}{\sqrt 2},
\]
the hopping Hamiltonian becomes
\[
\mathcal H=\hat c^\dagger \hat c-\hat d^\dagger \hat d,
\]
making the equally spaced spectrum in each fixed-\(k\) sector manifest [2605.04598].

## 3. Integrability, Bethe ansatz, and explicit eigenstate constructions

The two-site Bose-Hubbard model is one of the standard bosonic systems accessible to the algebraic Bethe ansatz. In one integrable parameterization,
\[
\hat H=(\hat N_1-\hat N_2)^2-\mu(\hat N_1-\hat N_2)-E_J(\hat a_1^\dagger \hat a_2+\hat a_2^\dagger \hat a_1),
\]
the construction starts from the rational \(gl(2)\)-invariant \(R\)-matrix, a bosonic Lax operator, and the monodromy matrix \(T(u)\). The Bethe vector state is
\[
|\{v_i\}\rangle=\prod_{i=1}^N C(v_i)\,|0\rangle,
\]
with Bethe roots satisfying
\[
\eta^2(v_i^2-\omega^2)=\prod_{j\neq i}^{N}\frac{v_i-v_j-\eta}{v_i-v_j+\eta}.
\]
A binomial-operator method rewrites \(C(v_j)\) as \(f_j\hat a_1^\dagger + D\), expands the Bethe vector in powers of \(D\), and yields a fully explicit combinatorial Fock-basis expression. The same construction gives scalar products, norms, and imbalance-current form factors as finite sums in symmetric polynomials of the Bethe roots and recursively generated coefficients \(C_{n,j}\) [1503.07885].

A generalized two-site Hamiltonian with interaction matrix \(K_{ij}\), local terms \(U_i-\mu_i\), and asymmetric tunnelling \(\Omega_{12}\neq\Omega_{21}\),
\[
\hat{H} = \sum_{i,j=1}^2 K_{ij}\hat{N}_i\hat{N}_j - \sum_{i=1}^2 (U_i - \mu_i)\hat{N}_i - \sum_{i\neq j}^2 \Omega_{ij}\hat{a}_i^\dagger \hat{a}_j,
\]
remains exactly solvable by a new parametrization of a bosonic Lax operator. The transfer matrix generates commuting conserved quantities, and the Hamiltonian can be written as a combination of them. In the no-interaction limit \(K_{ij}\to 0\) with \(\eta\ll 1\), the Bethe equations collapse to a sphere constraint,
\[
\sum_{i=1}^N \left[v_i+\frac{1}{2}(\omega_1-\omega_2)\right]^2 = R_N^2,
\]
which becomes an \((N-1)\)-sphere \(S^{N-1}\subset \mathbb R^N\) if all roots are real [1610.07111].

For the attractive model, Quantum Inverse Scattering Method results can be converted into high-accuracy closed energy formulas. Using the auxiliary Hamiltonian
\[
\hat{H}=a^\dagger b+ab^\dagger + \Delta b^\dagger b+ c^2 a^\dagger a\, b^\dagger b,
\qquad c^2=\frac{U-V}{J},\quad \Delta=\frac{2\epsilon}{J},
\]
the exact energies are written in terms of Bethe roots \(\lambda_j^\sigma\). The ground-state roots are approximately real, negative, and nearly equidistant,
\[
\lambda^{0}_n \approx -\frac{\Delta}{c}-c(n-1),
\]
which yields a simple approximation to the ground-state energy and, after mapping back to the physical Hamiltonian, accurate formulas for both the ground and first excited states over a broad parameter range. The first excited state changes character at \(c^2=\Delta\), reflecting a reorganization of the root pattern [1708.08696].

## 4. Dynamical regimes: Josephson oscillations, self-trapping, pair tunnelling, and revivals

In the tunnelling-dominated regime the basic observable is the population imbalance. For the \(N\)-boson Fock-space chain the normalized imbalance can be written as
\[
P(t)=\sum_{l=0}^N \frac{N-2l}{N}|c_l(t)|^2.
\]
For \(U=0\), the onsite term \(V_l\) vanishes and the effective Fock-space lattice becomes exactly self-imaging because the spectrum is equally spaced. If one initializes the system at one boundary of Fock space, the imbalance oscillates between \(-1\) and \(1\) with period
\[
z_R=\frac{\pi}{J}.
\]
As \(U\) increases, the oscillations are distorted and eventually become self-trapped: the probability distribution remains localized near the initially occupied side, corresponding physically to interaction-induced persistence of imbalance [1101.4762].

The \(N=2\) sector isolates correlated pair transport. In the basis \(|2,0\rangle\), \(|1,1\rangle\), \(|0,2\rangle\), the amplitudes satisfy
\[
i \frac{dc_0}{dt} = -\sqrt{2}J c_1 + Uc_0,\qquad
i \frac{dc_1}{dt} = -\sqrt{2}J(c_0+c_2),\qquad
i \frac{dc_2}{dt} = -\sqrt{2}J c_1 + Uc_2.
\]
As \(U/J\) increases, the intermediate state \(|1,1\rangle\) is suppressed and the bosons tunnel together as a pair. This is the simplest setting in which the two-site model already exhibits second-order pair tunnelling rather than merely independent single-particle transfer [1101.4762].

In the current-algebra formulation, exact second-order Heisenberg equations distinguish parameter regimes particularly clearly. In the extreme Rabi regime \(K/\mathcal E_J\ll N^{-2}\), and practically \(K\to 0\), the symmetric case \(\Delta\mu=0\) gives uncoupled simple-harmonic motion for \(\hat{\mathcal I}\) and \(\hat{\mathcal J}\) with frequency \(\mathcal E_J/\hbar\), while \(\hat{\mathcal T}\) is conserved and the period is \(T=2\pi\hbar/\mathcal E_J\). When \(\Delta\mu\neq 0\), \(\hat{\mathcal J}\) remains harmonic with frequency \(\hbar^{-1}\sqrt{(\Delta\mu)^2+\mathcal E_J^2}\), but \(\hat{\mathcal I}\) and \(\hat{\mathcal T}\) become coupled and interfere [1505.06793].

A semiclassical analysis of the weak-coupling regime \(u=UN/J<1\) and large \(N\) shows that the occupation difference exhibits rapid oscillations together with collapses and revivals. For the initial condition with all particles on one site, the collapse time, revival time, and broader blurring time are
\[
T_c=\frac{\sqrt{2N}}{Ju},\qquad
T_R=\frac{\pi N}{uJ},\qquad
T_B\approx \frac{\pi\sqrt{2}\,N^{3/2}}{3u^2J}.
\]
The first-order \(u\) correction sets the revival period, while the \(u^2\) terms control the detailed width and shape of the revival peaks [1409.5115].

Extensions with nonlinear atom-pair tunnelling add further dynamical structure. In the repulsive, symmetric case of the extended Hamiltonian
\[
H = U_1 \hat n_1^2 + U_2 \hat n_2^2 -\frac{1}{2}\Delta(\hat n_1-\hat n_2) -\frac{J}{2}\left(\hat a_1^\dagger \hat a_2 + \hat a_2^\dagger \hat a_1\right) -\frac{\Omega}{2}\left(\hat a_1^\dagger \hat a_1^\dagger \hat a_2 \hat a_2 +\hat a_2^\dagger \hat a_2^\dagger \hat a_1 \hat a_1\right),
\]
the semiclassical fixed-point analysis identifies Josephson, phase-locking, and self-trapping states, with phase boundaries \(\gamma=-\tfrac12\), \(\gamma=\lambda\), and \(\gamma=-1-\lambda\) in the \((\lambda,\gamma)\) plane [1606.08816].

## 5. Ground states, fragmentation, cat states, and entanglement

For the attractive biased two-site Hamiltonian, exact diagonalization reveals a sequence of qualitatively distinct ground-state regimes as \(\Lambda=NU/J\) becomes more negative. For weak attraction, \(|\Lambda|<2\), the ground state remains close to a single mean-field condensate and the Fock distribution is roughly binomial. For \(N=50\), a fragmented or cat-like window appears roughly in the interval \(-3.2 \lesssim \Lambda \lesssim -2.6\): the ground and first excited states become nearly degenerate, the Fock-space decomposition develops two peaks, and the ground and first excited states are approximately symmetric and antisymmetric combinations of the same bimodal structure. For stronger attraction the external bias selects one well, \(z\to 1\), \(\sigma_z\to 0\), and the cat-like superposition disappears. A central point is that fragmentation begins before a fully developed cat state appears [1108.4851].

To capture this crossover, an improved variational ansatz uses a superposition of two exchanged mean-field states,
\[
|\Psi_{\rm var}\rangle = A|\phi\rangle_1 + B|\phi\rangle_2
= \frac{A}{\sqrt{N!}}\left(\alpha a_L^\dagger+\beta a_R^\dagger\right)^N|0\rangle
+ \frac{B}{\sqrt{N!}}\left(\beta a_L^\dagger+\alpha a_R^\dagger\right)^N|0\rangle.
\]
Because \(A\) and \(B\) are independent, the ansatz interpolates continuously between a single condensate, a nearly symmetric cat-like state, and a bias-selected localized state. In the pre-bifurcation regime it captures the onset of fragmentation missed by a single mean-field state, and in the cat-like window it reproduces the exact Fock distribution, imbalance, fluctuations, and one-body density-matrix eigenvalues substantially better than earlier fixed-amplitude cat ansätze [1108.4851].

A complementary analytical treatment uses a permanent-based interpolation controlled by the single-particle overlap \(\omega\). In the repulsive regime this gives explicit formulas for the energy, visibility, Fisher information, number fluctuations, decay rates, and entanglement entropy in terms of hypergeometric polynomials. In the attractive regime it replaces naive symmetry-broken mean-field states by a symmetric superposition of two imbalanced coherent-like states, yielding a finite-\(N\) Schrödinger-cat description. Within that framework the semiclassical symmetry-breaking point
\[
U_c=-\frac{t}{(N-1)}
\]
is valid only in the limit of infinite number of bosons, while for large but finite \(N\) the relevant crossover scale is
\[
U_c^*\simeq U_c\left(1+\frac{1}{\sqrt{2N}}\right).
\]
This distinction addresses a recurrent misconception: at finite \(N\), the exact ground state remains symmetric even when the semiclassical equations admit broken-symmetry solutions [1108.6188].

The same permanent-based approach yields the relation
\[
F=\frac{4\,\sigma}{N^2}
\]
between Fisher information and on-site number variance for both repulsive and attractive interactions. The entanglement entropy is maximal close to \(U_c^*\) and exceeds its coherent value throughout the interval \(2U_c\lesssim U<0\). In the NOON limit the variance reaches \(N^2/4\) and \(F=1\), while in the coherent limit \(F=1/N\) [1108.6188].

## 6. Optical realization, open-system variants, and generalized dimers

A particularly direct realization of the two-site Bose-Hubbard dynamics is obtained in classical optics. The fixed-\(N\) Fock-space equations are exactly reproduced by light propagation in an engineered finite waveguide array of \(N+1\) sites,
\[
i\frac{d\phi_l}{dz} = -\left(\kappa_l \phi_{l+1}+\kappa_{l-1}\phi_{l-1}\right) +V_l\phi_l,
\]
under the identification \(\phi_l(z)\leftrightarrow c_l(t)\) and \(z\leftrightarrow t\). Each waveguide represents one occupation-number partition, the spacing \(d_l=x_l-x_{l-1}\) controls the coupling \(\kappa_{l-1}\), and the index contrast \(\Delta n_l\) controls the propagation-constant shift \(V_l\). The mapping is exact at the level of Fock-space amplitudes, so the optical power distribution \(|\phi_l(z)|^2\) directly visualizes \(|c_l(t)|^2\) [1101.4762].

This photonic construction makes several many-body phenomena spatially visible: Josephson-like oscillations, self-imaging at \(U=0\), the onset of self-trapping as \(U\) increases, and correlated pair tunnelling in the three-waveguide \(N=2\) realization. Its limitation is equally precise: the analogy is classical and linear, so it does not reproduce measurement backaction, coupling to an environment, or genuinely nonclassical photon statistics. It is exact only for the linear coupled-mode equations associated with a fixed-\(N\) Bose-Hubbard sector [1101.4762].

Generalized two-site models preserve much of the same algebraic structure. Rewriting the Hamiltonian
\[
\hat{H} = \sum_{i,j=1}^2 K_{ij}\hat{N}_i\hat{N}_j - \sum_{i=1}^2 (U_i - \mu_i)\hat{N}_i - \sum_{i\neq j}^2 \Omega_{ij}\hat{a}_i^\dagger \hat{a}_j
\]
in current variables produces coefficients \(\alpha\), \(\beta\), \(\gamma\), \(\rho\), and \(\xi\) that separate interaction nonlinearity, interaction asymmetry, bias, and collective energy shifts. In the special case \(\alpha=\beta=\gamma=0\), \(\hat{\mathcal T}\) is conserved and \(\hat{\mathcal I}\) and \(\hat{\mathcal J}\) are independent harmonic oscillators, recovering pure Rabi/Josephson oscillations. When \(\alpha=\beta=0\) but \(\gamma\neq 0\), \(\hat{\mathcal J}\) remains an independent oscillator while \(\hat{\mathcal I}\) and \(\hat{\mathcal T}\) stay coupled [1511.05026].

Open-system and nonlinear variants extend the model beyond closed Hamiltonian evolution. In a two-mode system with pair damping, thermo field dynamics and Hartree-Fock decoupling reduce the master equation to an \(su(1,1)\times su(1,1)\) squeezing problem. The short-time dynamics generates two-mode squeezing, hence entanglement, while the purity-like quantity \(Tr[\rho^2(t)]\) decays faster as the damping increases [1302.5238]. In the hopping-only dimer, evolution under \(\mathcal H^2\) rather than \(\mathcal H\) generates two-component Schrödinger cat states; at time \(\pi/2\), a coherent state becomes a superposition of coherent states with opposite amplitudes [2605.04598].

Taken together, these constructions show why the two-site Bose-Hubbard model retains a special status. It is small enough to be analytically tractable, but large enough to exhibit interaction-driven localization, macroscopic coherence, fragmentation, near-degenerate cat-like structure, integrable root geometry, nontrivial current dynamics, and exact analog simulation in synthetic lattices.

Source: https://www.emergentmind.com/topics/two-site-bose-hubbard-model