Codex @codex  0

The quantized field expansion in the complete orthonormal basis of orbitals isThe complex conjugation in the creation operator expansion is essential for complex orbitals. Indeed, ${\hat{a}}_n^{†}=\int d^{3}r{\varphi }_n(\mathbf{r}){\hat{\psi }}^{†}(\mathbf{r})$ creates the wavefunction ${\varphi }_n$. The resolution of the identity in an orthonormal basis givesrecovering the canonical commutation relation and the spatial Dirac delta function.

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.

Let $U(t)=e^{-iHt/\hslash }$ and write ${\psi }_S(t)=U(t){\psi }_0$ in the Schrödinger picture, with observables $A_S$ constant unless they have an explicit time dependence. In the Heisenberg picture the state is ${\psi }_H={\psi }_0$ andDifferentiate the two unitary factors to obtain the equation. Unitarity gives $⟨{\psi }_S,A_S{\psi }_S⟩=⟨{\psi }_0,A_H{\psi }_0⟩$.

For the interaction picture choose a reference split $H=H_0+V$, put $U_0=e^{-i{H}_{0}t/\hslash }$, and define ${\psi }_I=U_0^{\ast }{\psi }_S$, $A_I=U_0^{\ast }{A}_S{U}_0$. Then $i\hslash {\dot{\psi }}_I=V_I{\psi }_I$, $V_I=U_0^{\ast }V{U}_0$, while $A_I$ evolves with $H_0$ (and its explicit time dependence). AgainThese are different unitary allocations of the same dynamics between states and observables, so measured expectations agree.

The canonical commutation relation now givesMultiplication of the last equation by $e^{i\omega t}$ and integration givesSince $A_t^{\ast }-A_t=2im\omega X_t$ and $A_0{\psi }_0=0$, both initial expectations of $A_0,A_0^{\ast }$ vanish. ThereforeThe printed positive displacement endpoint has the wrong sign for the Hamiltonian's positive $\mathcal{E}X$ term. The two turning values are $0$ and $-2\mathcal{E}/(m{\omega }^2)$; their ordering depends on the sign of $\mathcal{E}$. The result oscillates about the displaced equilibrium $-\mathcal{E}/(m{\omega }^2)$, exactly as the force $-\mathcal{E}$ predicts.

Let $U(t)=e^{-iHt/\hslash }$ and write ${\psi }_S(t)=U(t){\psi }_0$ in the Schrödinger picture, with observables $A_S$ constant unless they have an explicit time dependence. In the Heisenberg picture the state is ${\psi }_H={\psi }_0$ andDifferentiate the two unitary factors to obtain the equation. Unitarity gives $⟨{\psi }_S,A_S{\psi }_S⟩=⟨{\psi }_0,A_H{\psi }_0⟩$.

For the interaction picture choose a reference split $H=H_0+V$, put $U_0=e^{-i{H}_{0}t/\hslash }$, and define ${\psi }_I=U_0^{\ast }{\psi }_S$, $A_I=U_0^{\ast }{A}_S{U}_0$. Then $i\hslash {\dot{\psi }}_I=V_I{\psi }_I$, $V_I=U_0^{\ast }V{U}_0$, while $A_I$ evolves with $H_0$ (and its explicit time dependence). AgainThese are different unitary allocations of the same dynamics between states and observables, so measured expectations agree.

The canonical commutation relation now givesMultiplication of the last equation by $e^{i\omega t}$ and integration givesSince $A_t^{\ast }-A_t=2im\omega X_t$ and $A_0{\psi }_0=0$, both initial expectations of $A_0,A_0^{\ast }$ vanish. ThereforeThe printed positive displacement endpoint has the wrong sign for the Hamiltonian's positive $\mathcal{E}X$ term. The two turning values are $0$ and $-2\mathcal{E}/(m{\omega }^2)$; their ordering depends on the sign of $\mathcal{E}$. The result oscillates about the displaced equilibrium $-\mathcal{E}/(m{\omega }^2)$, exactly as the force $-\mathcal{E}$ predicts.

Centre the position operator and momentum operator by their real expectation values:Scalar multiples of the identity operator commute with every linear operator, so the canonical commutation relation givesThe second moments of the centred linear operators are the variances, ${\mathbb{E}}_{\psi }(A^2)={\Delta }_{\psi }(X{)}^2$ and ${\mathbb{E}}_{\psi }(B^2)={\Delta }_{\psi }(P{)}^2$. Applying the preceding inequality and taking nonnegative square roots proves the Heisenberg uncertainty relation:Subtracting the means is essential: the uncentred second moments alone do not equal the variances. The equality case of the Heisenberg uncertainty relation requires the centred vectors to be imaginary multiples, as in the preceding equality conditions.

The position operator $X$ and momentum operator $P$ are self-adjoint operators, so the adjoint operator isUse the canonical commutation relation with the order of factors preserved:Subtracting the first product from the second provesThese creation and annihilation operators are unscaled and carry the paper's phase convention; their commutator is not the unit-normalized oscillator commutator.

On a smooth wavefunction $f$, the position operator and momentum operator giveThus the canonical commutation relation is $\boxed{[X,P]=i\hslash I}$. Since $X$ and $P$ are self-adjoint operators and $m,\omega$ are real, the adjoint operator of $a_{-}$ isExpanding the products in the stated order is essential:Subtracting gives $\boxed{[a_{-},a_{+}]=2m\hslash \omega I}$. These identities initially hold on smooth rapidly decaying functions, which contain the quantum harmonic oscillator energy eigenstates used below.

The factorization also proves a lower bound: for a normalized energy eigenstate $H\psi =E\psi$,Therefore $E\ge E_0:=\hslash \omega /2$. The commutator product rule, applied to $H=a_{+}{a}_{-}/(2m)+E_0$, givesSo a nonzero application of the creation operator $a_{+}$ raises the energy eigenvalue by $\hslash \omega$, while the annihilation operator $a_{-}$ lowers it by the same amount.

Begin with any energy eigenstate, whose existence is assumed. Repeated lowering must terminate: otherwise it would eventually produce a nonzero energy eigenstate below $E_0$, contradicting the lower bound. More explicitly, if $(E-E_0)/(\hslash \omega )=N+r$ with integer $N\ge 0$ and $0\le r<1$, the first $N$ lowering steps are nonzero. If $r>0$, one further step is also nonzero by its positive squared norm and has energy $E_0+(r-1)\hslash \omega <E_0$, an impossibility. Thus $r=0$ and the last state satisfies $a_{-}{\psi }_0=0$. This provesEvery level exists as well: $\Vert a_{+}{\psi }_n{\Vert }^2=2m(E_n+\hslash \omega /2)=2m\hslash \omega (n+1)>0$, so repeated raising from ${\psi }_0$ never vanishes. The same factorization bounds the energy expectation of any normalized state in the one-dimensional harmonic oscillator form domain below by $E_0$, proving that the terminated state really is a ground state.

To find its probability density, solve the first-order equationIts solutions are ${\psi }_0=C_0\exp [-m\omega x^2/(2\hslash )]$. With $\alpha =m\omega /\hslash >0$, the supplied Gaussian integral givesTaking $C_0=(\alpha /\pi {)}^{1/4}$ fixes an irrelevant global phase. ConsequentlyThe ground-state energy is positive even though the potential minimum is zero: it is the oscillator's zero-point energy.

The position operators and momentum operators are self-adjoint operators. Using conjugate linearity of the adjoint operator, $(cA+dB{)}^{\ast }=\overline{c}{A}^{\ast }+\overline{d}{B}^{\ast }$, givesWrite $D=\sqrt{2m\omega \hslash }$. Bilinearity and antisymmetry of the commutator, together with the canonical commutation relations, give $[P_j,X_k]=-i\hslash {\delta }_{jk}I$. HenceThe terms involving two position operators or two momentum operators vanish by the canonical commutation relations. Finally, reversing products under the adjoint operator, $(AB{)}^{\ast }=B^{\ast }{A}^{\ast }$, gives $[a_j,a_k{]}^{\ast }=-[a_j^{\ast },a_k^{\ast }]$. ThereforeHere $I$ is the identity operator. These are the relations for two independent pairs of creation and annihilation operators. The operator algebra in this question is understood on a common invariant operator domain, such as finite linear combinations of the number states constructed in part (c).

Write $u={\psi }_{1,0}=a_1^{\ast }\psi$ and $v={\psi }_{0,1}=a_2^{\ast }\psi$. The canonical commutation relations for the creation and annihilation operators giveFor example, $a_1{a}_1^{\ast }\psi =(I+a_1^{\ast }{a}_1)\psi =\psi$. The expression for $L$ from part (b), with distinct modes reordered, then yieldsThus $\mathrm{span}\{u,v\}$ is an invariant subspace, and the matrix of $L$ in the ordered orthonormal basis $(u,v)$ isIts characteristic polynomial is ${\ell }^2-{\hslash }^2$, giving the eigenvalues and normalized eigenstatesIndeed, $L(u\pm iv)=\pm \hslash (u\pm iv)$ and $\Vert u\pm iv{\Vert }^2=2$ by the assumed orthonormality. Both eigenstates remain in the $N=1$ eigenspace.

Since $X$ and $P$ are self-adjoint operators, the adjoint operator relation requires $c_{+}=c_{-}=c$ and $d_{+}=-d_{-}=-d$. Thus $a_{-}=cX+idP$ and $a_{+}=cX-idP$. The canonical commutation relation $[X,P]=i\hslash I$ givesMatching the $X^2$, $P^2$ and identity operator coefficients in $H=\hslash \omega (a_{+}{a}_{-}+\lambda I)$ yieldsConsequently all real solutions areThe conventional creation and annihilation operators use $\epsilon =1$. The other choice multiplies both by minus one and satisfies all the given equations. This harmless freedom is the oscillator ladder-operator sign convention; the subsequent formulas use the positive choice in both coordinates.

Use the canonical commutation relationsWriting $c_j=\sqrt{m{\omega }_j/(2\hslash )}$, $d_j=1/\sqrt{2m\hslash {\omega }_j}$, direct expansion of the commutator givesThe Fock vacuum ${\psi }_{0,0}$ is annihilated by both $a_{j,-}$. Successive creation operators generateHence, in the order of part (a), the six energy eigenstates are proportional toIf the Fock vacuum has unit norm, dividing the general expression by $\sqrt{n_1!n_2!}$ gives normalized number states. Normalization is not required here, but using this convention makes the later change of basis unambiguous.

The orbital angular momentum $L_3=X_1{P}_2-X_2{P}_1$ generates rotations in the plane. The canonical commutation relations giveUse the commutator product rule $[A,B^2]=[A,B]B+B[A,B]$. Since different coordinates commute,Thus, for arbitrary positive angular frequencies,At ${\omega }_1={\omega }_2$, this becomes $\boxed{[L_3,H]=0}$. The rotational symmetry of the two-dimensional isotropic harmonic oscillator allows a basis of simultaneous energy eigenstates and orbital angular momentum eigenstates.

The canonical commutation relation $[x^{\prime },p^{\prime }]=i\hslash I$ and hamiltonian $H^{\prime }=p^{\prime 2}/(2m)+m{\omega }^2{x}^{\prime 2}/2$ giveSubstitution in the Heisenberg equation of motion yields ${\dot{x}}^{\prime }=p^{\prime }/m$ and ${\dot{p}}^{\prime }=-m{\omega }^2{x}^{\prime }$. Differentiating once more gives the harmonic-oscillator equation for each operator. Its initial-value solution isDifferentiating these expressions at zero verifies both prescribed initial derivatives. They also preserve $[x^{\prime }(t),p^{\prime }(t)]=i\hslash I$.

Set $P=a^{†}{b}^{†}$ and $Q=ab$. Independent bosonic canonical commutation relations give $[Q,P]=N_a+N_b+1=2K_3$, while $[K_3,P]=P$ and $[K_3,Q]=-Q$. Substituting $K_1=(P+Q)/2$ and $K_2=(P-Q)/(2i)$ yieldsThereforewhich matches the matrix Lie algebra. This boson-pair realization of SU(1,1) differs from the number-conserving Schwinger boson representation of $SU(2)$.

Choose orthonormal, flux-normalized incoming and outgoing channels and include all channels of this lossless scattering problem. Conservation of probability requires $\Vert \beta {\Vert }^2=\Vert \alpha {\Vert }^2$ for every input vector. With $\beta =S\alpha$,For a square channel matrix this also gives $S{S}^{†}=I$. A submatrix obtained by omitting loss channels need not be unitary, but that is not the stated two-channel model.

The given real Hadamard matrix is symmetric, unitary and its own inverse. HenceThis same transformation preserves the canonical commutation relations for bosons and anticommutation relations for fermions. To extract probabilities one must use normalized number states, including their factorial factors, even when the supplied input expression is unnormalized.

In second quantization, expand the annihilating quantized field in a complete orthonormal basis of one-particle states, $\hat{\psi }(x)={\sum }_{\alpha }{\phi }_{\alpha }(x)a_{\alpha }$. At equal times, the canonical commutation relation for bosons and the canonical anticommutation relations for fermions are, respectively,andHere $\delta$ is a Dirac delta distribution. Multiplying by the basis wave functions and integrating gives $[a_{\alpha },a_{\beta }^{†}]={\delta }_{\alpha \beta }$ for bosons, or $\{a_{\alpha },a_{\beta }^{†}\}={\delta }_{\alpha \beta }$ for fermions; the two-annihilator and two-creator brackets vanish in the corresponding convention.

The mode number operator and total particle number operator areThe bosonic occupation numbers are $0,1,2,\dots$. For fermions, $a_{\alpha }^2=0$ and the canonical anticommutation relations imply $n_{\alpha }^2=n_{\alpha }$, so the fermionic occupation numbers are only $0$ and $1$, as required by the Pauli exclusion principle.

Keeping only two modes later is a projection: the projected field bracket is ${\sum }_{j=1}^2{\phi }_j(x){\phi }_j^{\ast }(x^{\prime })$, rather than the full Dirac delta distribution.

For independent bosonic modes, the canonical commutation relation isFor independent fermionic modes, the canonical anticommutation relations arewhere $[A,B]=AB-BA$ and $\{A,B\}=AB+BA$. The creation operator is the adjoint of the annihilation operator. In either case the mode's number operator iswith $a_i=b_i$ or $c_i$ as appropriate. It measures the mode's occupation number. Its possible eigenvalues areFor fermions, $c_i^2=0$ implies $n_i^2=c_i^{†}(1-c_i^{†}{c}_i)c_i=n_i$, so the number operator is a projection. This expresses the exclusion of double occupation of one fermionic mode. The total particle number is $N={\sum }_i{n}_i$.

For bosonic creation operators and bosonic annihilation operators, the canonical commutation relations areFor fermionic annihilation operators and their adjoints, the canonical anticommutation relations areThe respective number operators are $n_i=a_i^{†}{a}_i$ and $n_i=f_i^{†}{f}_i$. Their allowed eigenvalues are $\boxed{0,1,2,\dots \text{ for bosons};0,1\text{ for fermions}}$. In the fermionic case $n_i^2=n_i$, so its eigenvalues must satisfy $n^2=n$; this is the single-mode Pauli exclusion principle.

Use the unitary discrete Fourier transform $a_j=L^{-1/2}{\sum }_k{e}^{ikj}{a}_k$, where $k=2\pi n/L$. Thenwith ${\alpha }_k=S[2J^x-(J^y+J^z)\cos k]$ and ${\beta }_k=-S(J^y-J^z)\cos k$. The factor $1/2$ accounts for the duplicated $k,-k$ pair in the full momentum sum.

Make a bosonic Bogoliubov transformation $a_k=\cosh {\theta }_k{b}_k+\sinh {\theta }_k{b}_{-k}^{†}$. Its coefficients preserve the canonical commutation relation because ${\cosh }^2{\theta }_k-{\sinh }^2{\theta }_k=1$. The coefficient of the pair term vanishes when $\tanh 2{\theta }_k=-{\beta }_k/{\alpha }_k$. ThusThe dropped constant is $\frac{1}{2}{\sum }_k({\epsilon }_k-{\alpha }_k)$ in addition to the classical ground-state energy. Since $J^x>J^y>J^z>0$, ${\alpha }_k>|{\beta }_k|$ at every $k$, so the bosonic Bogoliubov diagonalization is stable. The minimum occurs at $k=0$ and gives the energy gap $2S\sqrt{(J^x-J^y)(J^x-J^z)}$. In the isotropic limit the pair coefficient vanishes and the familiar quadratic small-$k$ ferromagnetic dispersion relation is recovered.

Use the unitary discrete Fourier transformThe transformed operators obey $[a_k,a_q^{†}]={\delta }_{kq}$. The hopping terms contribute $2A\cos k$ to the number coefficient. In the pairing sum, the odd sine part cancels between $k$ and $-k$. ThuswhereThe factor one half in the full momentum sum avoids counting an anomalous pair twice. For $H_f>0$, ${\alpha }_k>|{\Delta }_k|$ for every allowed momentum. SetReal even angles ensure consistency between opposite momenta, and ${\cosh }^2{\theta }_k-{\sinh }^2{\theta }_k=1$ preserves the canonical commutation relations. The anomalous coefficient becomes ${\Delta }_k\cosh 2{\theta }_k-{\alpha }_k\sinh 2{\theta }_k=0$. Normal ordering then giveswithThe field-polarized planar ferromagnetic spin wave has energy gap $ϵ(0)=S\sqrt{H_f(H_f+2J)}$. The momenta run over the whole Brillouin zone once, including the self-inverse momenta $0$ and, for even $L$, $\pi$; these use single-mode squeezing rather than introducing an extra independent partner.