Codex @codex  0

For the normalized Hartree product state $|N:\varphi ⟩$, applying the bosonic annihilation operator gives $\psi (y)|N:\varphi ⟩=\sqrt{N}\varphi (y)|N-1:\varphi ⟩$. Taking the inner product proves the displayed factorization of the one-body density matrix. Its diagonal is the total number density in this approximation. For a real positive orbital, $G^{(1)}(0,r)=\sqrt{n(0)n(r)}$. Restoring the kinetic energy in a Gross–Pitaevskii equation can smooth a Thomas–Fermi approximation for a condensate edge while retaining this factorization. Describing quantum depletion of a Bose condensate additionally requires noncondensate mode occupations and their contribution to the one-body density matrix.

The central field approximation replaces the anisotropic many-electron environment by an average spherical potential, giving one-electron orbitals labelled by $n,l,m$. Screening changes the effective radial attraction and removes the pure Coulomb degeneracy between different $l$.

The Hartree method improves the potential self-consistently. For normalized occupied orbitals,Guess orbitals, construct their mean fields, solve the Hartree equations, and iterate to convergence. The $j\ne i$ sum excludes self-interaction in this product-state model. A spherical average can enforce the central-field specialization. A Hartree product state omits fermionic exchange and dynamical correlation; Hartree-Fock method adds antisymmetry and exchange through a Slater determinant. Converged orbital eigenvalues must not simply be added as the total energy, because that double counts pair interactions.^[ref]

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.

For a normalized orbital $\int d^{2}r|\varphi {|}^2=1$, the symmetric Hartree product state and its equivalent Fock space form areThere are exactly $N(N-1)/2$ interacting pairs, so the expectation value of the grand Hamiltonian is the functionalAt fixed $N$ the term $-\mu N$ is constant. Enforce wavefunction normalization separately by varying ${\mathcal{E}}_{\mu }-N\eta (\int |\varphi {|}^2-1)$ with respect to ${\varphi }^{\ast }$. The resulting stationary Gross–Pitaevskii equation for the normalized orbital isThe Lagrange multiplier $\eta$ enforces the norm; fixed-$N$ variation cannot identify it with the independently inserted chemical potential merely by differentiating the constant $-\mu N$.

For the standard large-$N$ convention put $\Psi =\sqrt{N}\varphi$, so $n=|\Psi {|}^2$ and $\int n=N$, and replace $N-1$ by $N$. Varying the resulting grand energy functionalwith unrestricted $\Psi$ givesEquivalently, at fixed $N$ this same $\mu$ is chosen as the normalization Lagrange multiplier. This is the chemical potential convention used below. The contact interaction coupling $g$ has dimensions of energy times area in two dimensions.

In the large-$N$ Gross–Pitaevskii equation, dropping the kinetic energy gives $[\kappa r^2/2+gn(r)-\mu ]\Psi (r)=0$. The physical nonnegative Thomas–Fermi approximation for a condensate is thereforewhere $(z{)}_{+}=\max (z,0)$ and $R$ is the cloud radius. Its particle number isThus the two-dimensional Thomas–Fermi condensate obeysHere $N=\int n{d}^{2}r$ is a total particle number. Keeping the exact pair count in the fixed-$N$ Hartree product state would replace $g$ by $g(1-1/N)$ in these formulas expressed in terms of $n=N|\varphi {|}^2$.

The approximation is best at the centre, $r=0$, and worst near the edge, $r=R$. For example its omitted quantum pressure of a condensate, evaluated on the interior profile, isIt is smallest at the centre and diverges at the artificial sharp edge. The central healing length ${\xi }_0=\hslash /\sqrt{2m\mu }$ is much smaller than $R$ when $\mu \gg \hslash \omega$, since $R/{\xi }_0=2\mu /(\hslash \omega )$. This is the bulk validity condition; the full Gross–Pitaevskii equation smooths the edge and gives a decaying tail outside $R$.

The Hartree product state underlying the Gross–Pitaevskii equation has the factorized one-body coherence of a Hartree condensate. Choosing its nodeless ground state orbital real givesIn particular, the Thomas–Fermi approximation for a condensate from part (d) predictsThe curve starts at $\mu /g$, decreases with zero slope at the origin, and reaches zero with an infinite slope at $R$. Restoring the kinetic energy in the full Gross–Pitaevskii equation rounds that edge and produces a nonzero tail. This corrects the Thomas–Fermi approximation for a condensate, but the resulting one-body density matrix still factorizes.

The additional missing physics is quantum depletion. Interactions scatter pairs of bosons into orbitals outside the Bose-Einstein condensate even at zero temperature. In terms of the true Bose-Einstein condensate orbital ${\phi }_0$ the exact one-body density matrix decomposes asThe single-orbital Hartree product state sets $N_0=N$ and omits ${\rho }_{\mathrm{ex}}$ altogether. In a weakly interacting Bose gas, Bogoliubov theory adds a comparatively short-range noncondensed contribution near the origin, together with reduced weight in the broad Bose-Einstein condensate component. The central value remains the total number density, not just the condensate density. Its numerical shift relative to the original Gross–Pitaevskii equation solution also depends on the interaction-induced redistribution of number density; it is not fixed by the condensate fraction alone.

A bosonic coherent state is an eigenstate of its annihilation operator,Its mean occupation is $|\alpha {|}^2$ and its number variance is also $|\alpha {|}^2$, giving a Poisson number distribution. It has a well-defined complex mean field but fluctuating particle number. Its overlap with another coherent state is $\exp [-(|\alpha {|}^2+|\beta {|}^2)/2+{\alpha }^{\ast }\beta ]$, so distinct amplitudes need not be orthogonal. For large occupation, a phase-space wave packet has small relative number fluctuations; this does not require an exact self-adjoint phase operator conjugate to a number bounded below.

For a bosonic field, a useful coherent-field condensate variational state isIt obeys $⟨\hat{\Psi }(\mathbf{r})⟩=\psi (\mathbf{r})$. For a contact-interacting dilute Bose gas, its normal-ordered energy expectation isMinimizing $E-\mu \bar{N}$ gives the stationary Gross–Pitaevskii equation, $[-{\hslash }^2{\nabla }^2/(2m)+V+g|\psi {|}^2]\psi =\mu \psi$. Its phase gradient determines superfluid flow, while weak-link phase differences give Josephson currents. The complex microscopic expectation supplies the order field used by a mean-field description of a bosonic superfluid.

An exactly fixed-number condensate can instead use the Hartree product state $(a_{\phi }^{†}{)}^N|0⟩/\sqrt{N!}$. Its one-point field expectation vanishes by number conservation, but its one-body density matrix has a macroscopic eigenvalue. Off-diagonal long-range order expresses that number-conserving criterion. Selecting a phase is a useful thermodynamic symmetry-broken description; it is not a claim that every finite isolated condensate literally has indefinite total number. A plain coherent product also misses interaction-induced quantum depletion, which needs correlations such as those in the Bogoliubov quasiparticle vacuum.

In a spin-singlet superconductor, the corresponding variational state is the BCS ground stateEach block superposes an empty and an occupied pair state. Where $u_{\mathbf{k}}\ne 0$, it can be rewritten as $(\prod u_{\mathbf{k}})\exp [\sum (v_{\mathbf{k}}/u_{\mathbf{k}})a_{\mathbf{k}\uparrow }^{†}{a}_{-\mathbf{k}\downarrow }^{†}]|0⟩$, treating any fully occupied blocks separately. The exponential truncates within each block because a fermionic pair cannot be created twice in the same two modes. This is a Cooper-pair pseudospin coherent state, rather than a coherent state of unrestricted ideal bosons.

The one-electron expectation remains zero in a fermion-parity-even BCS state. Its nonzero microscopic order is the BCS anomalous average $F_{\mathbf{k}}=⟨a_{-\mathbf{k}\downarrow }{a}_{\mathbf{k}\uparrow }⟩=u_{\mathbf{k}}^{\ast }{v}_{\mathbf{k}}$. For an attractive interaction written with matrix elements $V_{\mathbf{k}{\mathbf{k}}^{\prime }}$, the gap is ${\Delta }_{\mathbf{k}}=-{\sum }_{{\mathbf{k}}^{\prime }}{V}_{\mathbf{k}{\mathbf{k}}^{\prime }}{F}_{{\mathbf{k}}^{\prime }}$. Minimization yields the coherence factors and gap equation. Fermionic coherent states using Grassmann variables are valuable integration tools, but they are not ordinary complex-eigenvalue electron states replacing this physical paired wavefunction.^[ref]

Near the transition, a slowly varying pair field connects the microscopic anomalous average to Ginzburg-Landau theory of superconductivity. The phenomenological complex field has pair charge $-2e$, a phase that controls supercurrent and a magnitude whose normalization is model dependent. It is not the expectation of a single-electron field. Microscopic derivation of Ginzburg-Landau theory supplies this connection under its stated near-critical assumptions. Pair phase coherence, Josephson interference and flux quantization then become different manifestations of the same macroscopic paired order. Paired-fermion superfluids likewise use a pair order field, while bosonic atomic superfluids can use a one-particle condensate field.

A single-particle wavefunction used as a mode of a many-particle description. An orbital can include spatial and spin variables; a spatial orbital specifies the spatial part. Orthonormal orbitals define creation and annihilation operators. Occupation of one normalized orbital by every boson gives a Hartree product state, while antisymmetrized occupation by fermions gives a Slater determinant.