Codex @codex  0

In the Bogoliubov quasiparticle vacuum, ${\alpha }_{\mathbf{k}}|0_{\alpha }⟩=0$, but the inverse Bogoliubov transformation gives $⟨b_{\mathbf{k}}^{†}{b}_{\mathbf{k}}⟩=v_k^2$. Consequently the quantum depletion of a Bose condensate isThis expression assumes zero temperature. At nonzero temperature, thermal Bogoliubov quasiparticles add to the depletion:Thus the printed formula is the quantum contribution, or a low-temperature approximation when the thermal contribution is negligible. It is nonnegative for a stable quadratic Hamiltonian and vanishes in the noninteracting zero-temperature limit.

For $H={\hat{p}}^2/(2m)+V(\hat{x})$, divide elapsed time $t$ into $J$ slices $\Delta t=t/J$ and apply the Trotter product formula. Inserting complete position and momentum bases gives the regulated phase-space path integralThe endpoint values are $x_0=x_i$, $x_J=x_f$. Each Gaussian momentum integration in a phase-space path integral can be performed explicitly, leaving the configuration-space path integralIts measure is the time-sliced limit of $(m/(2\pi i\hslash \Delta t){)}^{J/2}{\prod }_{j=1}^{J-1}d{x}_j$, with the usual infinitesimal convergence prescription. The integral sums amplitudes, not positive classical probabilities. Operator ordering remains part of the time-slice definition. Expanding a short-time kernel recovers the Schrödinger equation; composing kernels reproduces the time-evolution operator.

In the classical limit, rapid phase cancellation suppresses paths away from stationary action. Varying $S$ with fixed endpoints gives $m\ddot{x}=-V^{\prime }(x)$. The semiclassical result sums the classical-path phases multiplied by fluctuation functional determinants; for a quadratic Hamiltonian the fluctuation integral is Gaussian and this expansion is exact.

For quantum statistical mechanics, continue $t=-i\tau$ by a Wick rotation. The Euclidean action and thermal trace giveThe trace closes position-space paths. In a many-particle trace one also sums over endpoint permutations with the appropriate exchange signs. In a coherent state field formulation, bosonic fields are periodic and fermionic coherent states instead give antiperiodic Grassmann fields, hence bosonic and fermionic Matsubara frequencies $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. ^[ref]

The quantum-to-classical imaginary-time mapping views a local $d$-dimensional quantum problem as a classical-looking field theory in $d+1$ dimensions, with imaginary-time extent $\beta \hslash$. A genuine classical probability interpretation additionally requires a nonnegative weight; fermion signs and Berry terms can obstruct it.

For a smooth potential in the high-temperature or small-$\hslash$ limit, imaginary-time paths fluctuate only over the thermal de Broglie wavelength ${\lambda }_T=\sqrt{2\pi \beta {\hslash }^2/m}$. If this is short compared with the potential's variation scale, the diagonal free kernel givesThis classical phase-space limit of the thermal kernel is the classical phase space partition function for one degree of freedom. A dilute identical-particle gas also needs its $1/N!$ factor; when exchange is appreciable, classical statistics fails. As a check, a single quantum harmonic oscillator gives $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}\to 1/(\beta \hslash \omega )$, precisely the classical phase-space integral.

The spin commutation relations give $[S_{\mathrm{tot}}^z,S_n^{\pm }]=\pm S_n^{\pm }$. Using the product rule for a commutator,and the conjugate term also vanishes. The longitudinal field commutes with every $S_n^z$, so $[H,S_{\mathrm{tot}}^z]=0$. The Hamiltonian is invariant under the global spin rotation $U(\theta )=e^{-i\theta S_{\mathrm{tot}}^z}$: its continuous symmetry group is $U(1)$. Simultaneous energy and total-magnetization eigenstates can be chosen, and evolution preserves each fixed-magnetization sector. Local magnetization can still flow between sites; generally $[H,S_n^z]\ne 0$. In the fermion representation the conserved number will be $N_f=N/2-S_{\mathrm{tot}}^z$.

For a bulk bond, $K_{n+1}=K_n(1-2n_n)$ and $c_n(1-2n_n)=-c_n$. Reordering the remaining different-site operators then givesThe periodic spin boundary condition requires separate attention. Set $P={\prod }_{n=1}^N(1-2n_n)=(-1{)}^{N_f}$. The end bond becomes $S_N^{+}{S}_1^{-}=-P{c}_1^{†}{c}_N$, with the conjugate relation for its adjoint. Consequently the exact finite-ring Hamiltonian isfermion parity is conserved, so diagonalize in a sector $P=\pm 1$. There the end bond can be written like every other bond provided $c_{N+1}=-P{c}_1$. This parity-dependent boundary condition of the Jordan-Wigner chain is antiperiodic for even fermion number and periodic for odd fermion number. Simply assuming periodic fermions in both sectors changes the finite-size spectrum.

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]

For one coordinate with Hamiltonian $H={\hat{p}}^2/(2m)+V(\hat{x})$, the propagator is $K(x_f,t_f;x_i,t_i)=⟨x_f|e^{-iHT/\hslash }|x_i⟩$, where $T=t_f-t_i>0$. Divide $T$ into $M$ intervals $\delta t$, use the Trotter product formula to separate kinetic and potential evolution, and insert complete position and momentum bases at each slice. The momentum integral is Gaussian:Multiplying the short-time kernels gives the time-sliced configuration-space path integral,with $x_0=x_i$, $x_M=x_f$. Continuum notation abbreviates this limit as ${\int }_{x_i}^{x_f}\mathcal{D}x{e}^{iS[x]/\hslash }$, where $S=\int dt[m{\dot{x}}^2/2-V(x)]$. The measure includes the displayed normalization; it is not an unweighted count of paths. A small damping or analytic-continuation prescription defines the oscillatory integrals. The paths contribute complex amplitudes, so alternative histories interfere before probabilities are calculated. This is the central difference from a classical statistical mixture of trajectories.^[ref]

When actions vary rapidly relative to $\hslash$, stationary phase selects classical paths. Put $x=x_{\mathrm{cl}}+\eta$, with $\eta$ zero at the endpoints. Vanishing first variation gives $m{\ddot{x}}_{\mathrm{cl}}+V^{\prime }(x_{\mathrm{cl}})=0$. The quadratic term isThe resulting Gaussian fluctuation determinant is the semiclassical propagator prefactor. In one dimension, away from caustics,The branch convention makes $\nu =0$ for a direct free-particle path before any caustic. The Maslov index tracks the phase changes as fluctuations pass through a focal singularity. In several dimensions replace the endpoint derivative by the Van Vleck determinant. The approximation sums all relevant classical paths; real-time action need only be stationary, not a minimum. It fails in this simple form at a caustic or coalescing saddles, where a uniform approximation is needed. Quadratic actions have no higher fluctuation terms, making their Gaussian evaluation exact. For example $S_{\mathrm{cl}}=m(x_f-x_i{)}^2/(2T)$ gives the exact free-particle kernel $\sqrt{m/(2\pi i\hslash T)}{e}^{i{S}_{\mathrm{cl}}/\hslash }$. Tunnelling involves classically forbidden regions and can be described by imaginary-time saddles rather than real classical trajectories.

Quantum statistical mechanics follows by Wick rotation $t=-i\tau$. The imaginary-time evolution operator is $e^{-\tau H/\hslash }$, with Euclidean actionTaking a thermal trace identifies and integrates the endpoint positions:This thermal trace as a closed Euclidean path integral explains the periodic imaginary-time direction. A density-matrix element instead retains open endpoints. For many identical particles the trace also sums over endpoint permutations, with a plus sign for bosons and permutation parity for fermions. In coherent-state language this becomes periodic bosonic fields and antiperiodic Grassmann fields, with corresponding Matsubara frequencies.

The Euclidean weight resembles a classical Boltzmann weight. A time-sliced local quantum field in $d$ spatial dimensions becomes a classical-looking field model in $d+1$ dimensions, with anisotropic couplings and extra extent $\beta \hslash$. Thus zero-temperature quantum criticality can be treated using an infinitely extended imaginary-time direction, while nonzero temperature imposes a finite thickness. The dynamical scaling of time need not equal the spatial scaling: the dynamical critical exponent specifies their relationship. This quantum-to-classical imaginary-time mapping requires qualification when weights have signs or imaginary terms; such a functional integral is not automatically an ordinary positive classical probability distribution. Even with a positive weight, the mapping concerns equilibrium statistics, not literal dissipative real-time trajectories. For a free particle the Euclidean kernel is a diffusion kernel with diffusion coefficient $\hslash /(2m)$, illustrating the mathematical correspondence.

The classical high-temperature limit is a further approximation, rather than simply another name for imaginary time. If quantum level spacings and exchange effects are negligible, a one-particle partition function in $d$ dimensions becomesFor indistinguishable particles there is a factor $1/N!$ in the dilute classical gas, and exchange is negligible when the number density times the thermal de Broglie wavelength to power $d$ is small. For a harmonic oscillator, the exact quantum result $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}$ tends to $1/(\beta \hslash \omega )$ as $\beta \hslash \omega \to 0$, recovering the classical phase-space integral. stationary phase explains classical mechanics; the short imaginary-time, weak-exchange limit explains classical statistical mechanics. These are related but distinct limits of the same quantum formulation.

For a normal-ordered bosonic Hamiltonian, insert coherent-state resolutions between short imaginary-time evolution factors. The overlaps generate $\bar{\psi }{\partial }_{\tau }\psi$ and the Hamiltonian symbol is obtained by replacing creation and annihilation operators by $\bar{\psi }$ and $\psi$. Including the chemical potential gives, in the stated $\hslash =1$ convention,The trace imposes periodic imaginary-time boundary conditions,so bosonic Matsubara frequencies are ${\omega }_n=2\pi n/\beta$. The two complex coherent-state integration variables are treated independently before imposing the physical saddle contour. There is no extra zero-point ordering term: the supplied interaction is already normal ordered.

The spatial confinement is not specified beyond the volume. The uniform bulk results below use a homogeneous box with periodic spatial boundaries, or the interior of a large vessel where boundary effects can be neglected. A hard wall instead requires a vanishing field at that wall and an inhomogeneous saddle near it. Temporal periodicity follows from the thermal trace and does not by itself impose periodic spatial boundaries.

For occupations small compared with $S$, use $S^z=S-n$, $S^{-}\simeq \sqrt{2S}{a}^{†}$ and $S^{+}\simeq \sqrt{2S}a$. The longitudinal bond contributes $-J\Delta S^2+JS\Delta (n_m+n_{m+1})+O(S^0)$; the transverse pair terms contribute $JS(a_m{a}_{m+1}+a_m^{†}{a}_{m+1}^{†})+O(S^0)$. Thus the direct real-space expansion isTake $a_m=N^{-1/2}{\sum }_k{e}^{ikm}{a}_k$, with $k=2\pi j/N$ over one full Brillouin zone and the lattice spacing set to one in $k$. Symmetrizing $k$ and $-k$ changes the phase in the pair sum to $\cos k$, givingnormal ordering the second Nambu spinor component then yieldsThe printed matrix has a factor-of-two error on its diagonal. With a full-zone sum its printed $2\Delta$ would produce $4\Delta a_k^{†}{a}_k$, and even its vacuum expectation would differ from the direct expansion by $J\Delta NS$. A half-zone convention cannot repair both coefficients: one must then group each $k,-k$ pair and double its off-diagonal pair coefficient as well. The formula above follows from the specified Hamiltonian and has the expected gapless isotropic limit; subsequent physical results use this consistent expansion. Terms quartic in bosons first enter at order $S^0$.

For the easy-axis XXZ antiferromagnetic spin-wave spectrum, the consistent quadratic Hamiltonian in part (d) has $A=\Delta$ and $B_k=\cos k$. The inverse transformation acts on the Nambu matrix by a congruence. Its diagonal and off-diagonal entries become $A\cosh 2{\theta }_k-B_k\sinh 2{\theta }_k$ and $B_k\cosh 2{\theta }_k-A\sinh 2{\theta }_k$. The pair term vanishes whenFor $\alpha >0$ the transformation is real and finite for every $k$. Summing the two Nambu number operators over the full zone givesHere $E_0$ names the constant in the requested oscillator form; the actual quadratic ground-state energy is $E_{\mathrm{gs}}^{(2)}=E_0+\frac{1}{2}{\sum }_k{ϵ}_k=-J\Delta N{S}^2+JS{\sum }_k({\Omega }_k-\Delta )$. This distinguishes the normal-ordering constant from the total zero-point contribution. bosonic Bogoliubov diagonalization is a canonical, not an ordinary unitary, diagonalization of the Nambu matrix.

At $\alpha =0$ and $k\ne 0,\pi$, the same formula holds with ${\Omega }_k=|\sin k|$. The exact zero modes have $|\tanh 2{\theta }_k|=1$, so no finite squeezing angle diagonalizes them as positive-frequency oscillators. Treat collective rotations separately or retain a small anisotropy regulator before taking the limit. Finally, if one instead diagonalizes the literal erroneous printed matrix, it gives $\tanh 2{\theta }_k=\cos k/(2\Delta )$ and ${ϵ}_k^{\mathrm{print}}=2JS\sqrt{4{\Delta }^2-{\cos }^{2}k}$. That result remains gapped at $\alpha =0$ and is not the spin-wave spectrum of the Hamiltonian stated at the start.

Take a uniform nearest-neighbour harmonic chain with displacement $u_j$, conjugate momentum $p_j$, mass $m_0$, spring stiffness $\kappa$, lattice spacing $a$ and periodic length $\ell =Na$. These specify a standard model for the requested notes, not additional numerical exam data. Its Lagrangian and Hamiltonian areVariation gives $m_0{\ddot{u}}_j=\kappa (u_{j+1}-2u_j+u_{j-1})$. For a mode $e^{i(kja-\omega t)}$, periodicity requires $k=2\pi n/\ell$ andThe lattice makes $k$ periodic modulo $2\pi /a$, so choose one Brillouin zone. At long wavelength ${\omega }_k\simeq v|k|$, where $v=a\sqrt{\kappa /m_0}$ is the sound speed. The continuum limit of a harmonic chain expands $u_{j+1}-u_j\simeq a{\partial }_{x}u$ and replaces ${\sum }_j$ by $a^{-1}\int dx$. Defining mass density $\rho =m_0/a$ and one-dimensional elastic coefficient $K=\kappa a$ givesThe momentum density is $\Pi =\rho {\partial }_{t}u$, and the Hamiltonian density is ${\Pi }^2/(2\rho )+K({\partial }_{x}u{)}^2/2$. Translational invariance of $u$ makes total momentum conserved. The continuum approximation describes acoustic wavelengths large compared with $a$; it is not a justification for extending the lattice dispersion and zero-point sums to arbitrarily large wavevector.^[ref]

canonical quantization promotes the coordinates and momenta to operators with $[u_j,p_l]=i\hslash {\delta }_{jl}$ and the other equal-time commutators zero. In the continuum, $[u(x),\Pi (x^{\prime })]=i\hslash {\delta }_{\mathrm{per}}(x-x^{\prime })$. The nonzero normal modes can be represented by independent phonon operators obeying $[b_k,b_q^{†}]={\delta }_{kq}$. A useful lattice expansion iswhere $X=N^{-1}{\sum }_j{u}_j$ is the centre-of-mass displacement. The corresponding momentum isIncluding the centre-of-mass term restores the full site commutator: the nonzero-mode sum alone gives $i\hslash ({\delta }_{jl}-1/N)$. Substitution diagonalizes the Hamiltonian,A phonon is one quantum of a collective normal mode, carrying energy $\hslash {\omega }_k$ and lattice momentum $\hslash k$ modulo a reciprocal-lattice vector. Reality of $u_j$ relates its Fourier coordinates by $u_k^{†}=u_{-k}$; it does not remove the independent oppositely travelling phonon modes. Using one complete Brillouin zone avoids double-counting real standing-wave coordinates.^[ref]

The free centre of mass has no restoring force and must not be treated as an oscillator of zero frequency. Fixing it, weakly pinning it, or specifying a finite centre-of-mass volume defines the corresponding zero-mode factor. With a fixed centre of mass, the internal thermal partition function and occupations areThe internal energy is ${\sum }_{k\ne 0}\hslash {\omega }_k(⟨n_k⟩+1/2)$; its high-temperature heat capacity approaches $(N-1)k_B$. For a large one-dimensional chain at temperatures below the lattice cutoff but above its finite-size mode spacing,This one-dimensional acoustic phonon heat capacity shows that quantization turns classical normal waves into bosonic excitations with a linear low-temperature heat capacity in one dimension. The Euclidean field action $S_E=\frac{1}{2}{\int }_0^{\beta \hslash }d\tau dx[\rho ({\partial }_{\tau }u{)}^2+K({\partial }_{x}u{)}^2]$ is Gaussian, with periodic imaginary-time boundary conditions; its Fourier modes reproduce these oscillators. Absolute displacement correlations require the zero-mode convention, and the lattice regulates short-distance fluctuations.

For one bosonic mode, define the unnormalized bosonic coherent state $|z⟩=e^{z{a}^{†}}|0⟩$. The commutator $[a,a^{†}]=1$ gives $a|z⟩=z|z⟩$ and $⟨z^{\prime }|z⟩=e^{{\bar{z}}^{\prime }z}$. The normalized state is $e^{-|z{|}^2/2}|z⟩$; it is not orthogonal to a different coherent state. The coherent-state resolution of identity follows from Gaussian moments,For many modes take products of these states and measures. The complex amplitudes have the natural scale of square root occupation; coherent states turn normal-ordered matrix elements into classical symbols without replacing the quantum trace by a single classical configuration.

For fermions, $\{c_i,c_j^{†}\}={\delta }_{ij}$ and a nilpotent annihilation operator cannot have an arbitrary nonzero ordinary complex eigenvalue. Introduce independent Grassmann variables ${\eta }_i,{\bar{\eta }}_i$, anticommuting among themselves and with the fermionic operators, and define $|\eta ⟩=\exp (-{\sum }_i{\eta }_i{c}_i^{†})|0⟩$. Then $c_i|\eta ⟩={\eta }_i|\eta ⟩$ and $⟨{\bar{\eta }}^{\prime }|\eta ⟩=e^{{\sum }_i{\bar{\eta }}_i^{\prime }{\eta }_i}$. With an ordered Berezin integral measure satisfying $\int d\bar{\eta }d\eta \bar{\eta }\eta =-1$ for one pair,for an even operator $O$; this is the fermionic coherent-state trace. In one mode the trace integrand is $e^{-\bar{\eta }\eta }(O_{00}-\bar{\eta }\eta O_{11})$, whose integral is $O_{00}+O_{11}$. Without the sign reversal it would be the supertrace, not the thermal trace. These formal Grassmann-valued coherent states do not represent ordinary complex-valued fermion amplitudes.

To construct the many-body thermal path integral, write $\mathcal{Z}=\mathrm{Tr}{e}^{-\beta K}$ with $K=H-\mu \hat{N}$ and divide the imaginary-time interval $\mathcal{T}=\beta \hslash$ into slices of width $\delta \tau$. Insert the appropriate coherent-state resolution of identity between factors. For a normal-ordered Hamiltonian, the adjacent-slice kernel iscoherent-state time slicing makes the overlap and Gaussian measure combine into the first-order time derivative. The grand-canonical coherent-state action in the regulated continuum limit isFor bosons the integration variables are complex, with the conjugate labels on the original integration contour; for fermions the barred and unbarred Grassmann variables are independent. The continuum notation retains the time-sliced ordering prescription. For example, $n(n-1)=a^{†2}{a}^2$ has normal symbol ${\bar{z}}^2{z}^2$, whereas $n^2$ has that symbol plus $\bar{z}z$. Dropping that distinction changes a theory. The coherent-state thermal boundary conditions follow because the bosonic trace closes paths periodically and the fermionic sign in the trace closes them antiperiodically:The barred fields obey the same respective boundary conditions. Hence Matsubara frequencies are $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. The derivative term encodes the canonical structure; omitting it would remove the quantum dynamics.^[ref]

As a check, a free diagonal mode with $K=\xi a^{†}a$ gives $Z_B=(1-e^{-\beta \xi }{)}^{-1}$ for $\xi >0$, whereas $K=\xi c^{†}c$ gives $Z_F=1+e^{-\beta \xi }$. Gaussian integration produces an inverse determinant for bosons and a determinant for independent complex fermions, with their respective temporal boundary conditions. Differentiating $\log Z$ recovers Bose-Einstein and Fermi-Dirac occupations. Interactions become higher powers in the normal-ordered symbols, and sources generate time-ordered correlations. The coherent-state construction preserves statistics, ordering and boundary conditions together; a formal continuum expression without its measure and temporal prescription is incomplete.

Keep the electromagnetic units of the supplied Lagrangian, including its factor $c$, and regard $e$ as the signed coupling used there. Set $a_B=e{\Phi }_B/(2\pi c)$. The canonical momentum is $p_{\varphi }=I\dot{\varphi }+a_B$, so the Legendre transform givesSingle-valued periodic wavefunctions have normalized eigenstates $e^{in\varphi }/\sqrt{2\pi }$ with integer $n$. Defining dimensionless flux $f=e{\Phi }_B/(hc)$ gives $E_n={\hslash }^2(n-f{)}^2/(2I)$, and the trace is thereforeThe flux-threaded quantum rotor spectrum is periodic under $f\mapsto f+1$ after relabelling $n$, corresponding to flux period $hc/|e|$ in these units. A local gauge change can remove the magnetic vector potential from the differential operator but then imposes a twisted boundary condition with phase $e^{-2\pi if}$; it cannot remove the flux dependence from both the Hamiltonian and boundary condition simultaneously. In SI conventions the flux coupling has no $c$, so one must not mix the two unit systems.

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.

Let ${\int }_{\mathbf{p}}=\int d^{3}p/(2\pi {)}^3$. The conjugate momentum fields are ${\pi }_{\varphi }={\dot{\varphi }}^{†}$ and ${\pi }_{{\varphi }^{†}}=\dot{\varphi }$. Inserting the given mode expansion of a free field into $i\int d^{3}x({\dot{\varphi }}^{†}\varphi -{\varphi }^{†}\dot{\varphi })$ produces spatial Dirac delta factors. The $b^{†}{c}^{†}$ and $cb$ terms vanish after integration: they contain ${\omega }_{\mathbf{p}}-{\omega }_{\mathbf{q}}$ on ${\delta }^3(\mathbf{p}+\mathbf{q})$, and ${\omega }_{-\mathbf{p}}={\omega }_{\mathbf{p}}$. The number-conserving terms instead add their two frequency factors, cancelling the mode normalization. The un-subtracted ordered expression isThe last term is an infinite vacuum constant, not a finite physical vacuum charge. Choosing the Fock vacuum to have zero Noether charge, equivalently using normal ordering, givesA state containing $n_c$ quanta created by $c^{†}$ and $n_b$ by $b^{†}$ has Noether charge $n_c-n_b$. Thus the two families describe opposite unit Noether charges; with this source's current the $c^{†}$ quantum has positive Noether charge and the $b^{†}$ quantum negative Noether charge. This is an internal Noether charge; identifying it with electric charge requires a specified electromagnetic coupling. The vacuum subtraction of a complex scalar charge leaves all commutators unchanged.^[ref]

Scope of the expansion: a fixed-frequency expansion with time-independent oscillator coefficients is the free-field or asymptotic-field form. For arbitrary interacting $f$, it can be used as canonical oscillator coordinates at a chosen time, rather than as the exact interacting time dependence. For example, $f(s)=m^{2}s+\lambda s^2$ gives $(\square +m^2)\varphi +2\lambda {\varphi }^{\ast }\varphi \varphi =0$, which a generic free plane-wave superposition does not satisfy. The Noether charge and its generator action remain valid in the interacting theory; the symmetry does not depend on using a free Hamiltonian.