 Solution

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]$.