Codex @codex  0

The Holstein–Primakoff transformation converts a spin algebra into bosonic operators while retaining the finite physical spin Hilbert space. With $[a_n,a_m^{†}]={\delta }_{nm}$, writeThe Holstein–Primakoff occupation constraint is $0\le n_n\le 2S$. Acting on $|n⟩$, the raising matrix element is $\sqrt{n(2S-n+1)}$ and the lowering element is $\sqrt{(n+1)(2S-n)}$. Their squared difference is $2(S-n)$, reproducing $[S^{+},S^{-}]=2S^z$. Operator ordering inside the square roots is essential. For low boson occupation relative to $S$, expand them to obtain $S^{+}\simeq \sqrt{2S}a$ and $S^{-}\simeq \sqrt{2S}{a}^{†}$. Keeping quadratic terms gives the linear spin-wave approximation; higher terms describe magnon interactions. Large $S$ helps control the expansion, but a low-dimensional infrared divergence can invalidate it even at large finite $S$.

For the ferromagnet choose $H_F=-J{\sum }_n{\mathbf{S}}_n\cdot {\mathbf{S}}_{n+1}$, $J>0$, with unit lattice spacing and $\hslash =1$. Expand around the all-up state. A longitudinal bond contributes $-J{S}^2+JS(n_n+n_{n+1})$, and the transverse exchange contributes $-JS(a_n^{†}{a}_{n+1}+a_{n+1}^{†}{a}_n)$. After a Fourier transform,These are ferromagnetic magnons. The polarized vacuum is an exact ground state of the full ferromagnetic Hamiltonian: every bond is in its maximum total-spin sector, minimizing that bond energy. It therefore has no quantum zero-point reduction of its magnetization. The one-magnon dispersion is also exact, although interacting many-magnon states need not be free bosons. Rotating the polarization generates the maximum-total-spin ground-state multiplet. Selecting a definite direction in a thermodynamic symmetry-breaking description breaks $SU(2)$ to $U(1)$; the two transverse spin components form a canonical pair, producing a single quadratic type-B Goldstone boson. A small aligning field gaps it by the Zeeman energy.

At nonzero temperature the one-dimensional ferromagnet is different: its thermal magnon number contains $\int dk[e^{\beta JS{k}^2}-1{]}^{-1}$, with infrared behaviour $\int dkT/k^2$. The divergent magnetization reduction signals the absence of finite-temperature long-range order in an infinite short-range isotropic chain, consistent with the Mermin-Wagner theorem. It does not disallow the polarized zero-temperature ground states just constructed.^[ref]

For the antiferromagnet use $H_{AF}=+J{\sum }_n{\mathbf{S}}_n\cdot {\mathbf{S}}_{n+1}$, still $J>0$, on an even bipartite chain. The classical reference state alternates up and down. Rotate one sublattice through $\pi$ so that both local quantization axes point along their reference spins. In unrotated physical components, to leading order $S_A^z=S-a^{†}a$, $S_B^z=-S+b^{†}b$, $S_A^{+}\simeq \sqrt{2S}a$ and $S_B^{+}\simeq \sqrt{2S}{b}^{†}$. Choosing Fourier phases using the actual site positions givesThe reduced Brillouin zone has $M=N/2$ momenta and length $\pi$. The pair creation and annihilation terms express the fact that an antiparallel reference product state is not an eigenstate. A bosonic Bogoliubov transformation with $u_k=\cosh {\theta }_k$, $v_k=\sinh {\theta }_k$, $\tanh 2{\theta }_k=\cos k$ diagonalizes each two-mode block. With the sign convention $a_k=u_k{\alpha }_k-v_k{\beta }_{-k}^{†}$ and $b_{-k}=u_k{\beta }_{-k}-v_k{\alpha }_k^{†}$, the result isThus the harmonic antiferromagnetic modes have linear, rather than quadratic, long-wave dispersion. In an unfolded Brillouin-zone description the same result is represented by zeros at both zero and the antiferromagnetic wavevector; these are not extra independent degrees of freedom beyond the two-sublattice accounting.

The Bogoliubov vacuum has no $\alpha$ or $\beta$ excitations but contains original spin bosons, with $v_k^2=\frac{1}{2}(|\sin k{|}^{-1}-1)$. Its staggered magnetization is $m_s=S-M^{-1}{\sum }_{k\in \mathrm{RBZ}}{v}_k^2$. This quantum depletion of Néel order diverges logarithmically at small $k$ in one dimension. Linear spin-wave theory therefore does not establish a stable ordered antiferromagnetic ground state for this chain. Its nominal gapless dispersion is a result of expanding around an invalid ordered reference, not a general proof of the exact excitation spectrum. For the nearest-neighbour isotropic chain, integer-spin systems have a Haldane gap, whereas half-integer-spin systems are gapless; the spin-one-half chain in particular has fractional spin excitations rather than a simple dilute gas of ordered-state magnons. This distinction requires physics beyond the harmonic expansion, including the topological term in the continuum theory.^[ref]

Consider a conventional spin-singlet, weak-coupling, three-dimensional superconductor. Set $\hslash =k_B=1$ in the microscopic action and restore electromagnetic units below. An attraction $\lambda >0$ in an energy shell about the Fermi surface gives a Grassmann actionThe attractive sign permits pairing even when the single-particle Fermi gas is otherwise stable. A complex Hubbard–Stratonovich transformation introduces a pair field $\Delta$ proportional at the saddle to $\lambda ⟨{\psi }_{\downarrow }{\psi }_{\uparrow }⟩$. The action is then quadratic in fermions. Integrating them out gives the BCS pairing-field effective actionHere the trace uses one independent Nambu spinor pairing block; a fully doubled particle-hole representation instead carries its corresponding one-half factor. At a uniform saddle the Bogoliubov--de Gennes Hamiltonian can be written asThe negative branch describes the paired vacuum; excitations are Bogoliubov quasiparticles mixing particles and holes. Extremizing the effective action givesThe physical cutoff or a renormalized interaction is part of the microscopic model; an unregulated contact integral must not be used indiscriminately. With constant density of states $N(0)$ per spin and shell half-width ${\Omega }_c$, the zero-temperature weak-coupling solution is ${\Delta }_0\simeq 2{\Omega }_c{e}^{-1/[N(0)\lambda ]}$. The exponentially small scale follows from the logarithmic pairing susceptibility of states near the Fermi surface. For this ideal weak-coupling model, $2{\Delta }_0/(k_B{T}_c)\simeq 3.53$. The microscopic theory relates the gap, pairing and thermodynamic transition.^[ref]

Close to a continuous transition, for fields varying slowly compared with microscopic scales, expand the effective action in powers and gradients of $\Delta$. Its static part is the Ginzburg-Landau theory of superconductivity:The normalization of $\Delta$ fixes the coefficients; here $a(T)$ changes sign at $T_c$. The uniform minimum has $|{\Delta }_0{|}^2=-a/b$ for $a<0$, condensation free-energy density $-a^2/(2b)$, and characteristic coherence-length scale $\xi =\sqrt{K/|a|}$. Above the transition the minimum is at zero. The quartic stabilization and gradient cost are precisely the ingredients of Landau theory for a complex order parameter. This controlled near-transition expansion supplies a bridge between microscopic BCS theory and macroscopic magnetic response; its coefficients and validity depend on the material and fluctuation regime.^[ref]

Write $\Delta =(|{\Delta }_0|+\eta )e^{i\varphi }$. The radial fluctuation has positive local restoring curvature, whereas the phase has no local potential in a neutral theory with global $U(1)$ symmetry. Freezing the amplitude givesThe neutral case supports a gapless collective phase mode. In the charged case the electromagnetic field couples to this phase rigidity. Under $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$ and $\varphi \mapsto \varphi +q\chi /\hslash$, the displayed combination is invariant. Locally choosing a gauge with constant phase leaves a quadratic term in $\mathbf{A}$. Varying the free energy gives the transverse London current and magnetic screening scale,Combining this relation with static Maxwell equations gives exponential field penetration, the Meissner effect. Single-valuedness of the pair field around a loop with negligible current gives flux in units ${\Phi }_0=2\pi \hslash /|q|=\pi \hslash /e$. These are gauge-invariant physical consequences of the condensate stiffness.

For the dynamics, the neutral phase and density mode is reorganized by the Coulomb interaction and the dynamical electromagnetic field. In a bulk three-dimensional system the long-wave longitudinal mode has a plasma gap, rather than remaining the neutral sound mode. This is the Anderson-Higgs mechanism in a superconductor: a phase coordinate becomes part of the electromagnetic collective response instead of surviving as an independent massless physical particle. Local gauge transformations are redundancies of the description, so a literal breaking of a local redundancy is not an observable. Phase coherence, current response and magnetic screening provide the physical statements.^[ref]

A static Landau functional alone cannot determine that full dynamical spectrum. A two-dimensional film, for example, need not have the same long-wave plasma dispersion as a bulk charged system. Nor does a positive radial curvature automatically guarantee a sharp Higgs amplitude resonance: coupling to quasiparticles and other modes can damp it. The chosen mean-field phase is also a thermodynamic symmetry-breaking description; an isolated finite system can be formulated with definite particle number. The microscopic pairing saddle, the near-transition Landau expansion and the electromagnetic Anderson-Higgs response describe successive scales of the same superconducting state, with distinct assumptions at each step.

At wavelengths with $q^2/(4m{\rho }_0)\ll g$, neglect the density-gradient term in the quadratic density-phase action. Complete the square in its remaining density terms:The convergent Gaussian integral can be performed by a contour shift or by its Fourier-transform identity. Its determinant is independent of $\phi$ at this order. Combining that normalization with the saddle action into a constant $S_0$ givesThe saddle contribution to $S_0$ is $-\beta V{\mu }^2/(2g)$; the additional measure-dependent constant need not be assigned an invented numerical value. The density Gaussian extends the fluctuation integration beyond its exact positivity domain, which is justified only while fluctuations around ${\rho }_0>0$ are small. The smooth zero-winding assumption is also required as discussed above. A constant global phase zero mode gives a finite compact-group volume, handled separately from nonzero phonon modes.

Fourier transforming the Bose-gas phase-only action gives a kernel ${\omega }_n^2/g+({\rho }_0/m)q^2$. Analytically continuing ${\omega }_n\mapsto -i(\omega +i{0}^{+})$ gives its real-frequency poles,In the present units the excitation energy equals $\omega$. Each real Fourier coordinate is a harmonic oscillator with inertia coefficient $1/g$ and spring coefficient ${\rho }_0{q}^2/m$, hence frequency equal to the square root of their ratio. This is the direct relation to a quantum harmonic chain. The mode is a gapless sound phonon; its finite speed reflects both the compressibility $\partial \rho /\partial \mu =1/g$ and the phase stiffness ${\rho }_0/m$.

Retaining the density-gradient term provides a useful check. Its integration instead gives a temporal coefficient $1/A_q$, soThis is the Bogoliubov spectrum. Its small-$q$ limit is the derived sound branch; at large $q$ within the contact model it crosses towards particle-like behaviour. The crossover in the squared dispersion is $q\sim 2\sqrt{mg{\rho }_0}$, precisely the scale beyond which the discarded density gradient is not small.

The same action checks the dimensional restriction in the saddle analysis. At nonzero temperature its zero Matsubara mode gives phase fluctuations proportional to $\int d^{d}q/q^2$; at zero temperature integrating the frequency gives $\int d^{d}q/|q|$. Their infrared divergences invalidate uniform true long-range order for $d\le 2$ at nonzero temperature and for $d=1$ at zero temperature. The long-wave collective description can remain useful for correlations even where a global condensate saddle is not a literal symmetry-broken state. The contact coupling here is an effective low-energy parameter; this calculation does not determine ultraviolet renormalization or trap-specific surface modes.

The normal ordering of both linear maps puts every annihilation operator to the right, soIn particular the stated Hamiltonian has already subtracted the free vacuum energy. For a one-particle state, evaluate the commutators using $[AB,C]=A[B,C]+[A,C]B$:Acting on the Fock vacuum therefore givesThe extra factor $\sqrt{2E_{\mathbf{k}}}$ in the chosen relativistic normalization of a one-particle state is a scalar multiplier and does not change either eigenvalue. Applying the same argument to a multiparticle quantum state gives additive total energy and momentum.

The Fock vacuum condition and the creation operator commutator giveEquality of the square-root product and $2E_{\mathbf{p}}$ holds on the support of the Dirac delta function.

The positive-energy mass shell has the Lorentz invariant measureHere both $p^2$ and the four-dimensional volume element are invariant, and positive energy is preserved by a proper orthochronous Lorentz transformation. Consequently the three-dimensional Jacobian determinant is $J=E_{\mathbf{\Lambda }p}/E_{\mathbf{p}}$. Transforming the Dirac delta function yieldsso its factor exactly cancels that of $2E_{\mathbf{\Lambda }p}$. Thus the displayed inner product is invariant under simultaneous transformation of both momenta. This is also why the one-particle completeness measure is $d^{3}p/((2\pi {)}^{3}2E_{\mathbf{p}})$.^[ref]

For a Lagrangian density depending on $\varphi$, ${\partial }_t\varphi$, ${\partial }_i\varphi$ and ${\partial }_i{\partial }_j\varphi$, direct variation givesTake $\delta \varphi$ to have compact support, so every boundary term vanishes. One integration by parts for the first-derivative terms and two for the second-derivative term yield the Euler-Lagrange field equation with second derivativesThe plus sign on the last term is essential. For the action in a, it is particularly transparent to vary the Laplacian directly:Thus the classical equation of motion isHere ${\nabla }^4=({\nabla }^2{)}^2$ is the biharmonic operator. Substitution of the plane wave $\varphi =A{e}^{-i\omega t+i\mathbf{p}\cdot \mathbf{x}}$ givesA real classical solution combines complex-conjugate plane waves. Quantization assigns the positive branch to the excitation energy, giving the required dispersion relation. The negative-frequency branch supplies the creation operator part of the field, rather than negative excitation energy. Independently, the Hamiltonian in a gives $\dot{\varphi }=\pi$, $\dot{\pi }=-{\kappa }^2{\nabla }^4\varphi -m^2\varphi$, confirming the same sign.

Define $E_{\mathbf{p}}=\sqrt{{\kappa }^2|\mathbf{p}{|}^4+m^2}$ and ${\int }_{\mathbf{p}}=\int d^{3}p/(2\pi {)}^3$. The Hermitian quantum field has the mode expansion of a free fieldwhere the annihilation operators and creation operators obey the bosonic canonical commutation relationsBecause $\pi =\dot{\varphi }$, its mode expansion of a free field isThe two nonzero cross commutators then giveand the equal-time field-field and momentum-momentum commutators vanish. This verifies the factor $1/\sqrt{2E_{\mathbf{p}}}$ without requiring Lorentz invariance.

Substitution into the Hamiltonian, followed by normal ordering, yieldsConsequently $[:H:,a_{\mathbf{p}}^{†}]=E_{\mathbf{p}}{a}_{\mathbf{p}}^{†}$: one application of a creation operator adds precisely the required quasiparticle energy. Before normal ordering there is the usual constant zero-point vacuum energy. If $m=0$, the exactly zero-momentum mode has $E_0=0$ and is a free canonical coordinate, not an oscillator with a normalizable oscillator vacuum state; in a finite box handle that mode separately or introduce a small positive mass. The formula applies directly to all modes with $E_{\mathbf{p}}>0$.

Continue with $m=0$. Along a characteristic $x=x_0+t$, the chain rule givesThus the initial wave packet ${\psi }_R(0,x)=f(x)$ evolves as ${\psi }_R(t,x)=f(x-t)$. Similarly, along $x=x_0-t$, the left component is constant, givingA feature initially at $x_0$ therefore appears at $x_0+t$ for the right component and $x_0-t$ for the left component. The profiles translate without changing shape at the speed of light, which is 1 in these units. For positive-energy plane waves, the right component has $E=p>0$ and the left component has $E=-p>0$, another way to identify their propagation directions.

The conclusion concerns the massless Dirac equation. For $m\ne 0$, the components couple as in c, and a positive-energy packet has group velocity $p/\sqrt{p^2+m^2}$, rather than a universal speed $\pm 1$. For example, a massive positive-energy rest solution has $p=0$ and ${\psi }_R={\psi }_L$, so neither component alone describes an independently moving light-speed packet.

Use the combinations $a_{\pm }=a\pm b$ at each momentum. Part a gives $[Q,a_{+}]=0$, $[Q,a_{-}]=-2a_{-}$, while their creation operators obey $[Q,a_{+}^{†}]=0$, $[Q,a_{-}^{†}]=2a_{-}^{†}$. Applying b givesThe creation operators transform by the Hermitian adjoint formulas with $e^{+2i\lambda }$. At $\lambda =\pi /2$, both exponentials are $-1$, soSubstituting into the fixed-time mode expansion of a free field of the complex scalar field givesThus $C$ exchanges particles and antiparticles, as required for charge conjugation. To see the charge reversal explicitly, the conserved charge $Q_{\mathrm{em}}={\int }_{\mathbf{p}}(a_{\mathbf{p}}^{†}{a}_{\mathbf{p}}-b_{\mathbf{p}}^{†}{b}_{\mathbf{p}})$ obeys $C{Q}_{\mathrm{em}}{C}^{†}=-Q_{\mathrm{em}}$. The generator $Q$ in the question is distinct from this charge.

The two occupation terms in the Hamiltonian are exchanged, with the same relativistic energy $E_{\mathbf{p}}$, and thereforeIn fact $[Q,H]=0$ because its mixing terms exchange equal-energy modes. With normal ordering, $:Q:=2{\int }_{\mathbf{p}}{c}_{-,\mathbf{p}}^{†}{c}_{-,\mathbf{p}}$ for $c_{-}=(a-b)/\sqrt{2}$. Then $C=e^{i\pi N_{-}}$ fixes the vacuum state and satisfies $C^2=1$. The construction is a bosonic species-swap generator.

Keep $\lambda$ for the cubic coupling and $j$ for the external source. Expanding the interaction part of the scalar field path integral yieldsA functional derivative of the source exponential inserts $i\varphi (y)$. In particular, ${\delta }_j(y{)}^3{e}^{i\int j\varphi }=-i\varphi (y{)}^3{e}^{i\int j\varphi }$, so $-i\lambda \varphi (y{)}^3$ can be replaced by $\lambda {\delta }_j(y{)}^3$. The functional derivatives commute and can be taken outside the free path integral. Resumming the exponential series givesThis is a formal perturbation series; a real cubic scalar potential is unbounded and the expression does not assert a stable nonperturbative vacuum state.

To evaluate the free Gaussian path integral, use integration by parts to writeThe inverse kernel isCompleting the square with $\varphi =\chi -K^{-1}j$ changes the exponent to $i\int \chi K\chi /2-i\int j{K}^{-1}j/2$. The translation-invariant functional measure makes the remaining Gaussian path integral precisely $Z[0;0]$, henceFor $ϵ>0$, the exponent includes $-ϵ\int {\varphi }^2/2$, which damps the regulated finite-dimensional integrals. The Feynman i-epsilon prescription displaces the positive-energy pole below the real $p^0$ axis and the negative-energy pole above it, selecting the time-ordered vacuum expectation value. Indeed the Feynman propagator is ${\Delta }_F=iD$. Thus $iϵ$ specifies both convergence and the vacuum state boundary conditions, with $ϵ\to 0^{+}$ taken after integration.

Choose the zero of energy so that the translation-invariant vacuum state has $P^{\mu }|\Omega ⟩=0$. A one-particle state of rest mass $m_{\psi }$ obeys the relativistic energy-momentum relationThe translation covariance of a quantum field in c remains valid for the full interacting four-momentum operator. Acting its exponential functions on the two momentum eigenstates givesHere $p\cdot x=E_{\mathbf{p}}t-\mathbf{p}\cdot \mathbf{x}$. A different additive vacuum energy must be subtracted consistently from both energy eigenvalues before identifying the physical $p^0$.

Quantum decoherence arises when different system alternatives correlate with distinguishable environmental states. For an initially uncorrelated state, a measurement-like unitary time evolution can produceTaking the environmental partial trace givesThe environmental overlap suppresses locally observable quantum coherence. Scattered photons, surrounding gas and uncontrolled apparatus variables can carry away such records. The interaction often singles out robust states forming a pointer basis. Decoherence can be much faster than thermal energy relaxation because recording a macroscopic distinction need not transfer much energy.

The overall state can remain pure and evolve unitarily even when the system alone looks like a classical mixture. Discarding environmental information is therefore not, by itself, a proof of fundamental wavefunction collapse or a solution to the definite-outcome measurement problem. A small controlled environment can return coherence, as in the single-qubit environment dephasing model; practical irreversibility arises when recovering the many dispersed records becomes inaccessible. In an interferometer, varying gas pressure or an independently controllable path marker tests the dependence of fringe contrast on this environmental coupling. Quantum computation must suppress or correct such unwanted system–environment entanglement to preserve the intended interference.

Write $h_0=-{\hslash }^2{\nabla }^2/(2m)+V(\mathbf{r})$ for the one-particle Hamiltonian operator. The second quantization of the kinetic energy, external potential energy and repulsive contact interaction isThe factor $1/2$ counts each unordered particle pair once. Normal ordering excludes self-interaction: the interaction acts on two distinct bosons. The fixed particle number sector has $\hat{N}=\int d^{3}r{\hat{\psi }}^{†}\hat{\psi }=N$. In three dimensions, this contact interaction is understood as a regulated low-energy interaction with a matched coupling, as in renormalization of a three-dimensional dilute-gas contact coupling; the finite-mode calculation below uses its projected matrix elements.

The two-mode approximation for a Bose gas retains the real localized orbitals:It is valid when the lowest doublet is well separated from higher energy levels, and the temperature and interaction-induced changes of the orbitals are small on that separation. More explicitly, if $\Delta$ is the energy gap to discarded orbitals, sufficient low-energy conditions are $k_{B}T,t,N{U}_{0}I\ll \Delta$, where $I=\int d^{3}r|{\phi }_L{|}^4$. A deep double-well potential suppresses interwell interaction overlaps. Large $N$ alone does not justify the truncation: increasing $N$ at fixed coupling can violate the interaction condition.

Let $ϵ=(E_0+E_1)/2$ and use the even and odd eigenfunctions to evaluate the projected one-particle Hamiltonian operator:ThusThe same quantum tunnelling matrix element is $t=-\int d^{3}r{\phi }_L{h}_0{\phi }_R$. Reflection symmetry makes $\int |{\phi }_L{|}^4=\int |{\phi }_R{|}^4=I$. Neglecting mixed contact interaction integrals leavesThe last step is the normal-ordered boson pair count, following from the canonical commutation relation. Consequently the two-site Bose-Hubbard model has precisely the requested coefficientsandHere $\mu$ is the coefficient obtained by rewriting the fixed-particle number Hamiltonian operator, rather than an independently prescribed thermodynamic chemical potential. Adding a reservoir term $-{\mu }_{\mathrm{res}}\hat{N}$ would replace it by $\mu +{\mu }_{\mathrm{res}}$.

The negligible-overlap instruction applies to the mixed contact interaction integrals while the small quantum tunnelling matrix element is retained. Exactly disjoint smooth orbitals for a local $h_0$ would also give $t=0$; setting every overlap-related matrix element identically to zero would discard the physical splitting one is trying to describe.

The normalized Hartree product state and its equivalent second quantization form areThe product is already symmetric under exchanging bosons; no additional $1/\sqrt{N!}$ belongs in the coordinate wavefunction. The factor in the second quantization expression instead normalizes repeated creation operators acting on the vacuum state. This is a fixed-number two-mode condensate, with every boson in the same normalized orbital.

Set $c=\cos \theta$ and $s=\sin \theta$. Its binomial distribution of occupations givesThe expectation value of the Hamiltonian operator is thereforeKeeping $N(N-1)$ rather than replacing it by $N^2$ retains the exact pair count of the trial quantum state. Equivalently, with $ϵ=w-\mu$, its interaction contribution is $wN(N-1)(c^4+s^4)$ and its one-particle contribution is $ϵN-tN\sin (2\theta )\cos \chi$.

Choose the nonredundant parameter range $0\le \theta \le \pi /2$, with $\chi$ modulo $2\pi$. For $t>0$, minimizing over the relative phase gives $\chi =0$. Putting $x=\sin (2\theta )\in [0,1]$, the remaining nonconstant energy is $-tNx-\frac{1}{2}wN(N-1)x^2$, which decreases with $x$ for the repulsive coupling $w>0$. HenceThe optimal orbital is the symmetric ${\varphi }_0=({\phi }_L+{\phi }_R)/\sqrt{2}$. Both the quantum tunnelling term and the repulsive interaction favor this balanced occupation within the Hartree product state family. This result uses a variational method: for strong repulsion relative to quantum tunnelling, the true many-body ground state can have smaller occupation fluctuations than this fixed-number two-mode condensate.

The printed identity has the wrong sign for the specified relative phase. Differentiating the energy found above and evaluating the operator commutator giveThus the phase derivative of a two-mode condensate energy isA concrete counterexample to the printed sign is $N=100$, $\theta =\pi /4$, $\chi =\pi /2$ and $t>0$: its left side is $100t$, whereas its printed right side is $-100t$, independently of $w$ and $\mu$.

There is also a direct derivation without inserting the Hamiltonian operator. Since $|\theta ,\chi ;N⟩=e^{i\chi {\hat{n}}_R}|\theta ,0;N⟩$, differentiating this unitary operator gives ${\partial }_{\chi }E=i⟨[\hat{H},{\hat{n}}_R]⟩$. Conservation of the total particle number, $[\hat{H},{\hat{n}}_L+{\hat{n}}_R]=0$, then supplies the minus sign for the left well.

By the Heisenberg equation of motion, $\hslash d⟨{\hat{n}}_L⟩/dt=i⟨[\hat{H},{\hat{n}}_L]⟩$. Define the bosonic Josephson junction particle current from left to right as $J_{L\to R}=-d⟨{\hat{n}}_L⟩/dt$. The corrected identity meansThe energy slope with respect to the relative phase drives particle transfer; an extremum with respect to that phase has zero instantaneous particle current. The printed sign would work with $e^{-i\chi }$ in the right-well orbital, or with $[{\hat{n}}_L,\hat{H}]$ in place of $[\hat{H},{\hat{n}}_L]$.

For a uniform continuum gas with pair interactions depending on relative positions, a Galilean boost changes every particle momentum by $p_0\hat{\mathbf{x}}$ and leaves the interaction energy unchanged. Write ${\hat{P}}_x={\sum }_j{\hat{p}}_{j,x}$. The transformed many-body Hamiltonian operator isThe resting ground state has $⟨{\hat{P}}_x⟩=0$, so the Galilean boost energy of an interacting gas isThe mass here is the bare particle mass, because this is an exact transformation of the microscopic kinetic energy. Substituting the Fermi-liquid effective mass $m^{\ast }$ would omit the change in interaction energy encoded by the Landau quasiparticle energy functional. This conclusion requires the Principle of Galilean relativity; a gas in a lattice or a fixed external potential energy landscape need not obey it.