Condensate wavefunction Created 2026-10-20 Updated 2026-10-25
A complex amplitude describing a coherent condensate in mean-field theory. In the convention , its squared modulus is the condensate density. Dividing by gives a normalized orbital; formulas for interaction strengths and nonlinear terms must use the same normalization convention.
For a repulsive contact-interacting Bose-Einstein condensate, the classical Hamiltonian density generating evolution in a frame rotating at the chemical potential is . Its integral is the grand canonical ensemble generator , not the physical energy alone. For the uniform mean-field theory minima satisfy , with a freely chosen constant phase. In the physical frame, the corresponding field has time factor .
Mean-field approximation Created 2026-09-24 Updated 2026-10-18
The mean-field approximation replaces interactions with a self-consistent average field, neglecting correlated fluctuations around that average.
This is a use of mean-field theory, whose general framework also applies beyond a Landau free-energy model.
Mean-field theory Created 2026-10-13 Updated 2026-10-18
Mean-field theory replaces fluctuating interactions with a self-consistent average field. For example, an Ising spin feels a field determined by the average neighboring spin magnetization; a self-consistency equation approximates the interacting Ising model. Neglecting correlated fluctuations can fail near a continuous phase transition. The mean-field approximation in Landau theory is a particular use of this general framework.
A strong negative surface potential creates exponentially enhanced positive counterion density and exponentially suppressed negative co-ion density. The co-ion population remains positive, unlike an extrapolated linear approximation. Finite ion size and correlations eventually limit the applicability of point-ion mean-field theory.
Use . The given state is a two-species uniform boson number state: only the mode of each species is occupied. Applying an annihilation operator to an unoccupied mode gives zero, while has eigenvalue . Thus
The two species have independent fixed occupation numbers, so mixed expectation values factorize. The inverse discrete Fourier transforms yield
In the momentum-space Hamiltonian operator, only contributes: each annihilation operator must remove an occupied boson, and each creation operator must restore it to the same mode. Hence
There is no replacement of by in this energy: in the real-space calculation the relevant neighboring-site commutator is zero. Both expressions assume distinct neighbouring sites, as in ii.
To match the prescribed spin length in the mean, impose . At fixed integer , the energy is minimized by maximizing . Consequently the best choices within this relaxed family are
When is even, and . When is odd, either assignment of gives . Without fixing , this auxiliary boson Hamiltonian operator has no lower bound: taking sends its expectation value to .
Next, . Each term changes the two species' occupation numbers, producing a state orthogonal to . Therefore
For distinct sites modulo , expand the product of the two spin operators. The terms and again have zero expectation value. The remaining terms give
Here and the Kronecker delta vanishes at distinct sites. Consequently
The same holds for the components. At the balanced choice and as , each transverse two-point correlation function tends to independently of separation, whereas the one-point transverse magnetization is zero. This is transverse spin order in a number eigenstate: the state has a definite total and retains the global U(1) symmetry of rotations about , yet has the transverse correlations of an ordered state with no selected common azimuth. The sum of the and correlations tends to , describing ferromagnetic long-range order within this approximation. The thermodynamic limit is understood with both the separation and growing, without identifying the two sites periodically.
Fixed-spin qualification. The proposed state is not actually in the physical fixed- Hilbert space for . It satisfies the local spin constraint in the Schwinger boson representation only in expectation value. For example, with ,
The nonzero component has local total occupation numbers and hence local spins , instead of . More generally the local total number operator has variance
Each species contributes a binomial distribution at a given site, with success probability , and the two contributions are independent. In the physical subspace that variance must be zero. The proposed state even has
by the Casimir operator identity from a. Thus the calculated energy is a relaxed mean-field theory estimate, not a Rayleigh-Ritz variational principle bound for the stated fixed- quantum XY chain. Projecting onto at every site would give a physical trial state and generally change its expectation values; in the two-site example the projected state has energy , whereas the unprojected value is . The calculated long-range order characterizes the supplied approximation and does not establish long-range order in the exact one-dimensional ground state.
Full dissociation gives bulk concentrations in of for sodium, for potassium, for calcium, for chloride, and for hydrogen phosphate. The positive and negative charge concentrations are both , so the bulk is electrically neutral. At ,
The Boltzmann distribution then gives
For completeness the coion concentrations are and . They are depleted, not zero. Using and multiplying molar concentrations by the Avogadro constant in the Grahame equation gives
This calculation retains every species. Since , a linear Debye–Hückel approximation is not justified. The large accumulation, especially of divalent counterions, follows from the valence in the Boltzmann factor. However, multi-molar local concentrations also warn that point-ion mean-field theory neglects finite ion size, hydration, specific adsorption, correlations, and changes in local permittivity. These are predictions of the stated idealized model, not precise measured interfacial concentrations.
An ionic atmosphere. A negatively charged molecule attracts positive counterions and repels negative co-ions; signs reverse for a positive molecule. The resulting electric double layer partially screens the bare electric field. Far from an isolated molecule, the electrolyte recovers its neutral bulk composition. At thermodynamic equilibrium, spatial uniformity of each species' electrochemical potential requires
The logarithm implicitly uses a common reference density. This result for Boltzmann ion densities in a symmetric electrolyte neglects ion correlations, finite size and nonideal activities, but shows immediately why equal concentrations of the two signs do not persist close to a charged surface.
Self-consistent screening. Insert these densities into Poisson's equation, with solvent permittivity and fixed charge density :
This is the Poisson-Boltzmann equation. For a symmetric monovalent salt with bulk density of each sign, outside the fixed charge it becomes
For , linearizing and using bulk neutrality gives the Debye–Hückel approximation
Thus added salt shortens the Debye length: electrostatic screening becomes stronger over a shorter distance. The factor of two in a salt comes from two ionic species, not from double-counting the density of salt pairs.
For example, a charged sphere of radius and total charge has the exterior linearized solution
assuming fixed uniform surface charge and no ions inside. It follows from the decaying spherical solution and . At zero salt it reduces to the Coulomb potential; at finite salt it decays exponentially as well as geometrically. The excess mobile charge integrates to in an infinite neutralizing reservoir, although a finite sphere surrounding the molecule generally contains only part of that ionic atmosphere.
Biophysical scales and limitations. In water at room temperature, the monovalent formula is approximately , with the ionic strength. It gives about at and at . Consequently the same charged protein can exert very different long-range interactions mediated by electrostatics at different salt concentrations. Close to strongly charged DNA, the linear approximation fails; nonlinear screening, counterion condensation, ion-specific adsorption and changes of surface ionization through charge regulation may matter. Mean-field theory also misses strong multivalent-ion correlations and finite-size saturation: the exponential law cannot be extrapolated to arbitrarily large near-surface concentrations.
Define . Each site has exchange-connected neighbours. In mean-field theory, replacing neighbours by their mean spin produces the effective field that yields
Thus
for Mn, Fe, Co and Ni. Negative means antiferromagnetic exchange in this convention; a Hamiltonian written with instead of would have twice these numerical values. The increasingly negative values indicate a growing effective exchange energy, although different orbital occupancy and anisotropy mean that they are not a simple universal function of -electron number.
A three-dimensional six-connected framework is more favourable to mean-field reasoning than isolated chains, but the theory neglects short-range correlations, anisotropic exchange and fluctuations. The ratios show substantial departures from the simplest unfrustrated isotropic mean-field relation. The values are therefore effective high-temperature exchange estimates, not a complete microscopic fit to the ordered phases.
In the Weiss molecular field approximation, the local response obeys
Linearize around in the stable phase of paramagnetism. Substitution of a gives
Hence the Curie–Weiss law and Weiss temperature are
For positive , is the mean-field theory Curie temperature: the zero-magnetization solution becomes unstable there, and below it a spontaneous solution appears. For negative coupling, a negative fitted Weiss temperature indicates predominantly opposing interactions, but this one-uniform-magnetization treatment does not determine an antiferromagnetic ordering pattern or its transition temperature. More generally the fitted Weiss temperature need not equal the observed ordering temperature.
For , the paramagnetic curve is positive for , diverges as , and tends to at high temperature; its inverse is the straight line . The expression is not the stable susceptibility below . For example substituting would predict a negative response around the unstable branch. The nonlinear self-consistency equation must instead be expanded about the appropriate spontaneous magnetization.
Figure 1.
Paramagnetic Curie-Weiss susceptibility and its inverse
. The positive-Weiss-temperature mean-field law is shown only in its stable paramagnetic range . Its inverse extrapolates to zero at . The curve below the transition is not obtained by continuing the displayed rational expression. Details.
Ferromagnetic interactions enhance the response above the ordering transition. The Curie–Weiss law improves the independent-spin approximation:
with positive Weiss temperature and the same high-temperature coefficient . In a simple mean-field theory, take . More generally is fitted and need not equal , and the simple pole does not give exact critical behaviour arbitrarily close to the transition.
At zero temperature, a selected ferromagnetic orientation has all magnetic moments aligned. Its saturation magnetization is
Equivalently, there are per conventional cell. The moment points along either allowed easy-axis direction at zero magnetic field. This is the local or single-domain spontaneous magnetization; a macroscopic multidomain sample can have smaller, even zero, net magnetization in the absence of a selecting field.
Choose the positive spontaneous branch for . Its magnetization starts at , decreases with temperature, and reaches zero continuously at . It remains zero for . There is an equally valid negative branch below ; domain averaging can instead give zero net moment.
A fixed small positive magnetic field intensity selects the positive branch. At zero temperature it still saturates at ; at intermediate temperature it increases the magnetization relative to the zero-field branch. It rounds the zero-field transition and leaves a positive magnetization above , tending in the linear-response region to . The magnetization therefore does not vanish at in the nonzero-field curve.
For a concrete labelled sketch, the figure uses the illustrative mean-field theory equation
It shows and a fixed small . The curve is an illustrative equilibrium equation of state, not a determination of the precise critical exponents of GdN. At no symmetry-breaking transition remains along this temperature sweep.
Figure 1.
Spontaneous and field-selected magnetization in the two-direction model
. Details.
Let be uniform, with dimensionless spin operators . The factor in the source Hamiltonian removes double counting of pairs. Decoupling about a uniform polarization gives the one-site mean-field theory Hamiltonian
A body-centred-cubic site has eight nearest neighbors and six next-nearest neighbors, so
For a spin one half in this effective field, the mean-field self-consistency equation is
Linearizing at the continuous onset gives , hence
Equivalently, the general spin- result is . The antiferromagnetic next-nearest-neighbor exchange reduces the uniform ordering scale rather than contributing its absolute magnitude. The negative nuclear g-factor does not change this exchange instability; aligning the spins also aligns their moments with one another.
For these particular competing couplings, the uniform mode is indeed the largest exchange eigenvalue. Writing , and gives
with equality at modulo reciprocal lattice vectors. To see the bound, replace a negative by its absolute value for an upper bound, let , and use . The expression is at most , whose maximum on is five. Thus another modulation does not preempt the stated ferromagnetic onset at this level of mean-field theory. Fluctuations can change the actual transition temperature.
The Weiss temperature is the extrapolated zero of in a Curie–Weiss law fit. It summarizes the leading high-temperature effect of interactions between magnetic moments. For the sign convention
where are dimensionless angular momentum operators, mean-field theory gives
Thus positive net exchange interaction in this convention favours parallel magnetic moments, giving and tendencies towards ferromagnetism; negative net exchange interaction favours antiparallel magnetic moments, giving and antiferromagnetic tendencies. Competing interactions can partially cancel in this sum. The Weiss temperature need not equal the actual ordering temperature, and its sign alone is not a complete determination of the low-temperature ordered structure.
For competing exchanges on a body-centred cubic antiferromagnet, use the Heisenberg model convention , with each bond counted once, dimensionless spins, and positive antiferromagnetic . There are eight nearest and six next-nearest neighbours. Bond counting in the two classical collinear structures gives
For type II the eight contributions cancel, while all six bonds are antiparallel. Their crossing yields . One can also exclude a lower spiral within the classical bilinear model: writing gives
For a favourable sign of the product, . Hence when its minimum is at , the type-II ordering vectors. When , compare with the type-I value and use to obtain the complementary bound. Thus
This is a classical/mean-field energy comparison; an exact quantum ground state is not generally a collinear product state.
For type I a magnetic cell is the simple cube of side , with opposite moments at and . Its magnetic structure factor is proportional to . Allowed magnetic reflections have integer with odd. The first two families are and , with symmetry-equivalent members understood.
For type II choose magnetic primitive translations , and , enclosing four spins in volume . Equivalently, a side- cube is a convenient nonprimitive magnetic cell. On each original simple-cubic sublattice the moment at is proportional to . Accordingly,
The translated body-centre basis gives , so there is no further extinction. The first two families are and . The neutron polarization factor may additionally suppress a reflection for a particular moment direction; the stated rules are the positional selection rules.
Each sublattice has number density and saturation magnetization . Replacing neighbouring spins by their averages makes their exchange energy equivalent to a molecular magnetic field. Its minus sign expresses antiferromagnetic exchange, not an externally applied field. Since ,
For , . Linearizing the Brillouin magnetization formula gives
The g-factor cancels because exchange is specified in energy units. A convention with doubles these exchange energies and .
Let , and let denote the cubic coefficient of the Brillouin magnetization formula. The argument is . Thus implies
Consequently in mean-field theory. The positive branch falls from saturation at zero temperature to zero at and remains zero above it. An elastic magnetic neutron-reflection intensity scales as ; a log–log fit against has slope . Background subtraction and a suitable critical interval are needed. Fluctuations can change the measured exponent from its mean-field value; uniform bulk magnetization vanishes for a compensated antiferromagnet and is not itself this order parameter.
Figure 1.
Sublattice magnetization and its square near a mean-field antiferromagnetic transition
. Details.
To obtain the Boltzmann equation, take a dilute gas with a short-range pair potential. Its interaction range is small compared with the mean free path and the spatial scale on which the phase-space distribution function changes. The binary collision approximation treats encounters as isolated elastic collisions; simultaneous encounters involving three or more particles can be neglected. The encounter duration is short compared with the collision relaxation time and macroscopic evolution time, and the external force changes negligibly during an encounter. This permits local, effectively instantaneous scattering at one position .
The essential statistical closure is molecular chaos for the Boltzmann equation: the distribution of incoming pairs factorizes as . Outgoing pairs need not be independent. Follow each incoming pair through its two-body scattering, and use conservation of momentum and kinetic energy, together with time reversal and the invariant scattering measure, to express both gain and loss in the same variables. This gives
Here ; when the measure is written as , its Jacobian is included in the Boltzmann scattering kernel. In a scattering-angle convention, that kernel contains relative speed times differential scattering cross-section. The hemisphere prevents double counting equivalent outcomes for indistinguishable particles; the precise factor is included in . The primed velocities are connected to the unprimed ones by the reversible elastic collision map. This closure is appropriate for a dilute short-range gas, rather than a generic long-range interacting system with an important collective mean-field theory contribution.
The further replacement by the BGK model is a phenomenological relaxation closure. Replace the full nonlinear Boltzmann collision operator by relaxation towards a local Maxwellian, with a single, velocity-independent collision relaxation time :
It models relaxation towards a local Maxwellian without resolving the detailed scattering. Treating as constant is an additional simplification; it is of the order of a mean time between effective collisions. Near a local Maxwellian, departures in the nonconserved velocity moments decay on this timescale.
Choose to match the collision invariants , and of . For unit particle mass, this means
Here is the mass density, is the mean velocity, and is the temperature in squared-velocity units, before setting . In particular,
Thus relaxation preserves local mass, momentum and kinetic energy rather than relaxing them to externally prescribed values.
The Landau-Ginzburg theory for a scalar order parameter supplements a local symmetry-allowed potential with a gradient penalty. For a symmetry at zero field, a stable quartic model is
with near a mean-field transition. Equilibrium fluctuations have weight . In the Landau approximation, minimize the functional. A uniform stationary state obeys . At zero field the minimum is for and for . The latter pair is spontaneous symmetry breaking of the discrete symmetry.
The uniform free-energy density at its minimum is zero above the transition and below, up to regular backgrounds. Differentiating the equation of state gives and . At , . Gaussian fluctuations around the minimum have inverse covariance above and below. Their correlation lengths are and . These results give the usual mean-field critical exponents
The heat capacity has a finite jump rather than a power-law divergence. Negative quartic coefficient requires higher stabilizing terms and can produce a first-order transition; a vanishing quartic coefficient requires a tricritical analysis. Thus the assumed positivity of matters.
The Ginzburg criterion tests whether this saddle description is self-consistent. Below the transition compare fluctuations of the order parameter averaged over a correlation volume with its squared mean-field value . Modes with dominate this coarse-grained fluctuation estimate:
where the positive dimensionless depends on the averaging window and the factor of two used in . Dividing by gives
It is the fluctuation of a correlation-volume average, not the full cutoff-dominated microscopic variance, that should be compared with the mean-field order parameter.
For , diverges as , so the asymptotic critical region is fluctuation dominated. Away from it, mean-field theory can still work when
For , the long-wavelength ratio tends to zero near criticality, consistently with Gaussian control and mean-field thermodynamic exponents. At , quartic coupling is marginal; scale-dependent logarithms require a renormalization-group analysis rather than interpreting a constant bare ratio as a proof of no corrections. This identifies four as the upper critical dimension of the short-range scalar quartic model.
The Landau-Ginzburg theory is valuable because symmetries restrict a manageable list of local operators, the gradient term describes correlations and interfaces, and the Ginzburg criterion identifies when neglected fluctuations become decisive. Below the upper critical dimension an interacting renormalization-group fixed point, rather than the bare saddle, determines universal exponents. For continuous order parameters the manifold of ordered minima also has Goldstone modes, and lower-critical-dimension fluctuations require separate attention. Mean-field minimization supplies an approximation; the Ginzburg ratio decides whether its neglect of fluctuations is justified.
Reading the lowest point of the order-disorder boundary in the supplied experimental diagram gives approximately
(about –, allowing for the graphical resolution). The ordered region just above that minimum is lamellar.
This is a reading of the supplied polyisoprene-polystyrene phase diagram. It should not be replaced by the often quoted symmetric-chain mean-field theory value : that refers to a different idealization and is not the minimum drawn here.
A coarse-grained scalar order parameter describes a short-range Ising-like transition with sign symmetry. Write the dimensionless Landau-Ginzburg functional as
Here , with the relevant critical-mass shift understood when fluctuations are included. Higher local powers and gradients are negligible only within a justified long-distance expansion. The partition function is with a microscopic cutoff. The local functional is not already the fluctuation-integrated thermodynamic free energy.
Mean-field theory replaces this integral by its stable uniform saddle, . At ,
The branch curvature is : it is above and below. Hence
Together with the finite heat-capacity jump from , this yields
The equal- susceptibility amplitude ratio is two. Uniform saddles also miss domain walls and spatial texture. For an -component order parameter, symmetry allows a vector generalization; ordered transverse directions additionally carry Goldstone fluctuations.
To include small fluctuations set . The Gaussian action is , giving
The displayed is a spectral covariance per Fourier delta, avoiding an extra volume factor. As , a growing correlation length lets increasingly large correlated regions fluctuate. The singular Gaussian heat correction contains , which is infrared singular for and logarithmic at four. This signals failure of a simple saddle expansion; the Gaussian heat exponent obtained by extrapolation is not the exact interacting critical exponent.
The Ginzburg criterion compares the fluctuation of a correlation-volume average with , rather than using all ultraviolet fluctuations at a point. Keeping below the transition gives, up to finite convention-dependent constants,
Dividing by yields
For , writing gives a nonuniversal Ginzburg reduced-temperature window
Inside , interacting fluctuations must be treated. For example, in three dimensions in this normalization. Above four dimensions the long-wavelength ratio vanishes at criticality; at four it is marginal and interaction-generated logarithms require further analysis.
Figure 1.
Mean-field wells and the Ginzburg test
. Original scalar quartic potentials and dimensionless illustrative Ginzburg ratios. The three-dimensional ratio grows towards criticality, the four-dimensional power estimate is constant and the five-dimensional one falls. A unit-ratio line marks a chosen crossover convention, not measured material data. Details.
The spatial gradient term supplies correlations absent from uniform Landau theory. The Ginzburg test bounds where its mean-field saddle remains self-consistent; it does not compute nonclassical exponents or a universal material-independent window.[ref]
Use the reduced attractive reduced BCS pairing Hamiltonian with volume V, and attraction :
A specified pairing shell or an ultraviolet renormalization is part of the model. The trace is grand canonical. Insert fermionic coherent states into its time-sliced trace. The resulting Grassmann variables are antiperiodic, and with inverse-energy imaginary time the dimensionless action is
This boundary condition is essential for Fermi-Dirac statistics; replacing it by a bosonic periodic trace changes the answer. A complex Hubbard–Stratonovich transformation uses the local identity
with time-slice factors understood. The pairing field renders the fermion action quadratic. For a uniform static field the independent Nambu spinor block has
Its determinant is . The phase-independent determinant sign cancels under consistent reference normalization. Integrating the fermions gives the BCS pairing-field effective action, , with the normal-ordering/reference term retained. The fermionic pairing-block determinant normalization uses one independent pairing block for every k; a fully doubled spin-particle-hole trace instead needs one half. The subtracted grand potential of the reduced BCS model is particularly unambiguous:
Performing the Matsubara frequency sum yields
up to a common constant. At zero gap this reduces exactly to the two-spin free-gas result. The Matsubara derivation of the BCS saddle follows by differentiating with respect to and using gives
The second equation fixes the chemical potential at specified number. This stationary action is BCS mean-field theory, not an exact integration of pairing-field fluctuations. With constant single-spin density of states per volume and a symmetric energy cutoff , weak coupling gives and . The zero-temperature condensation grand-potential density is to leading weak-coupling order; using a spin-summed density would change these displayed factors. The BCS coherence factors are and , and the associated BCS ground state is the normalized product of empty and occupied pair states. [ref]
Allowing a slowly varying complex pairing field generates phase stiffness and, near the transition, the Ginzburg-Landau theory of superconductivity. A neutral paired system has a gapless phase sound mode; in a charged bulk superconductor electromagnetic coupling produces the Anderson-Higgs mechanism in a superconductor and the Meissner effect. These collective fluctuations are distinct from the gapped fermionic excitations. A microscopic dynamical calculation is needed to fix collective-mode dispersion and damping; the static saddle alone does not determine them. [ref]