Codex @codex  0

The functional derivative is $\delta \mathcal{E}/\delta h=A{h}^{(4)}$, so the overdamped filament bending equation is $\zeta h_t=-A{h}^{(4)}$ with the free-end conditions from (a). Expand the initial shape in the complete orthonormal basis:The overdamped relaxation of a free filament isEach bending mode relaxes on time ${\tau }_n=\zeta /(A{k}_n^4)$; the two zero modes remain constant. Consequently the long-time shape is the affine projection of $H$,The total displacement and its first moment are conserved: integrate $h_t$ and $(x-L/2)h_t$, using $h^{\prime \prime }=h^{(3)}=0$ at the ends. This also checks why the final affine part cannot generally be set to zero. The series defines the $L^2$ solution for arbitrary square-integrable initial data; a classical solution at time zero additionally requires compatible boundary values.

Normalize the vacuum path integral by $Z[0]$. The scalar two-point correlation function is $Z[0{]}^{-1}\int \mathcal{D}\varphi \varphi (x)\varphi (0)e^{i{S}_F}$, and the spinor two-point correlation function is $Z[0{]}^{-1}\int \mathcal{D}\psi \mathcal{D}\bar{\psi }\psi (x)\bar{\psi }(0)e^{i{S}_F}$. Vacuum boundary conditions make these time-ordered Feynman propagators. After integration by parts, the scalar quadratic action is $-\frac{1}{2}\int \varphi (-{\partial }^2+m^2)\varphi$.

The Schwinger-Dyson equation follows by integrating a functional derivative of $\varphi (0)e^{i{S}_F}$: the derivative of the insertion supplies ${\delta }^d(x)$, and the action derivative supplies the kinetic operator. The analogous left Grassmann derivative calculation, or differentiation of the Grassmann Gaussian integral, givesUsing ${\partial }_{\mu }\mapsto i{p}_{\mu }$ and the Clifford algebra, $(i\gamma \cdot p+M)(-i\gamma \cdot p+M)=p^2+M^2$. ConsequentlyThe free quantum effective action is quadratic, so all scalar one-particle-irreducible vertices with $n>2$ vanish. Its two-point vertex is the stated inverse kinetic form $-p^2-m^2$.

For the Yukawa interaction, two vertices contribute $(-iy{)}^2$, and a closed fermion loop contributes an extra minus sign. Tracing the two spinor numerators gives $4[M^2-k\cdot (k-p)]$, since the one-gamma traces vanish. Removing the overall $i$ from the amplitude gives the displayed loop integral. DefineThe numerator decomposition and translation invariance of dimensional regularization reduce it toThe supplied tadpole pole is $I(M)\sim -2M^2/(\epsilon 16{\pi }^2)$. A Feynman parameter combines the two bubble denominators. Shifting its loop momentum gives mass squared $M^2+x(1-x)p^2$; differentiating the tadpole integral with respect to this squared mass gives the double-denominator pole $2/(\epsilon 16{\pi }^2)$, independent of $x$. Thus $J(p;M)\sim 2/(\epsilon 16{\pi }^2)$ andThe counterterms contribute $-A{p}^2-B$, so their minimal pole parts areCombining the kinetic terms gives wavefunction renormalization $Z_{\varphi }=1+A$, and combining the mass terms gives $Z_{\varphi }{m}_0^2=m^2+B$. ThereforeThese are the Yukawa scalar self-energy pole coefficients; finite parts depend on the chosen renormalization condition.

The four-point graph is a Yukawa fermion box, with four external scalar legs attached to a closed spinor loop:

Each high-momentum dirac propagator is $O(k^{-1})$. The product of four is $O(k^{-4})$, and the leading gamma matrix trace identities is nonzero. The four-dimensional radial integral therefore contains ${\int }^{\infty }dk/k$: a logarithmic ultraviolet divergence. This local four-scalar divergence cannot be absorbed by scalar mass or field normalization. Add $-\lambda {\varphi }^4/4!$, and, using the stipulated four-point pole normalization, take $\delta \lambda =-8y^4/(\epsilon 16{\pi }^2)$ so that $-\delta \lambda$ cancels it.

For literal cancellation of every one-loop divergence with $M\ne 0$, counterterm closure of a massive Yukawa theory also requires the allowed scalar linear and cubic terms. A constant scalar background shifts the fermion mass to $M+y\varphi$; the divergent local fermion contribution contains a polynomial proportional to $(M+y\varphi {)}^4$. Its linear and cubic terms are not forbidden by a symmetry when the fermion mass is nonzero. A closed renormalizable family therefore haswith field, mass and coupling redefinitions for both scalar and spinor fields. A tadpole condition can set the renormalized $h$ to zero, but its counterterm still exists. The vacuum constant is needed if vacuum energy is retained. If an exact discrete chiral symmetry is imposed with $M=0$, the scalar potential can be even and the odd terms are forbidden; the essential new interaction is then the quartic one.

Using functional derivatives of the interaction functional,Dividing the flow by $e^{-S_1}$ therefore givesIn the first term, remove one leg from each of two interaction vertices and join them with a line weighted by ${\dot{C}}_{\Lambda }(p)$. This produces the tree joining of two vertices. In the second, remove two legs from one vertex and contract them with that line, raising the loop order by one, with tadpole diagrams as the simplest example. The factor one half accounts for interchanging the contracted ends; the displayed signs are the signs in the interaction-action flow.

The varying cutoff replaces an internal propagator by its cutoff derivative. Repeated tree joins and loop closures express how eliminated high-momentum fluctuations generate the vertices of the Wilsonian effective action. The action is not restricted to one-particle-irreducible Feynman diagrams: connected tree joins also occur. Field-independent vacuum contributions can again be absorbed into normalization.

Take an oriented closed string in flat, critical bosonic string theory, with ${\alpha }^{\prime }=1/(2\pi T)$ and a mostly-plus target Minkowski metric. After continuation of the worldsheet to Euclidean signature, the Polyakov path integral sums over embeddings $X$ and worldsheet metrics, divided by worldsheet diffeomorphisms and Weyl transformations. In a flat target its kinetic action isTarget-time continuation or analytic continuation of external momenta defines the Lorentzian scattering amplitude; a naive real Euclidean Gaussian for timelike $X^0$ would not be convergent.

An external tachyon is represented by the tachyon vertex operator $V_p(z,\bar{z})=:e^{ip\cdot X}:$, with conformal weights $h=\bar{h}={\alpha }^{\prime }{p}^2/4$. The physical integrated vertex has $(h,\bar{h})=(1,1)$, so $p^2=4/{\alpha }^{\prime }$. Schematically its tree amplitude isThe Faddeev-Popov determinant from conformal gauge is represented by the worldsheet ghost fields. At tree level the worldsheet is a Riemann sphere, whose unpunctured complex structure has no moduli. Its residual conformal automorphisms are the Möbius transformations, $PSL(2,\mathbb{C})$. Fix three insertion points, accompanying their unintegrated vertices by the required $c\tilde{c}$ ghost factors. The remaining $N-3$ complex insertion positions are integrated over the sphere; equivalently one integrates all positions and divides by the residual conformal group.

The embedding fields are free, withTheir zero-mode integral gives momentum conservation, and their nonzero-mode Gaussian integral gives the Koba-Nielsen factorThis explains the sphere tachyon position integral without needing its evaluation. For four tachyons the resulting Virasoro–Shapiro amplitude displays the exchanged string spectrum directly.

Use all-incoming external momenta and introduce $u=-(p_1+p_4{)}^2/(8\pi T)$ alongside the two printed dimensionless Mandelstam variables. Since $p_a^2=8\pi T$ and ${\sum }_a{p}_a=0$, one obtains $s+t+u=-4$. The physical center-of-mass energy squared in the $s$ channel is $s_{\mathrm{phys}}=8\pi Ts=4s/{\alpha }^{\prime }$. The $t$ channel measures the analogous crossed momentum transfer, with sign determined by the mostly-plus convention. Rewriting the Gamma factors in a symmetric form givesThe Gamma function poles imply, at generic fixed values of the other invariant,The denominator Gamma factors can remove residues at special intersecting channel kinematics; the statement concerns a generic single-channel limit. An $s$-channel pole occurs when the intermediate momentum $p_1+p_2$ satisfies the mass-shell condition for a closed-string state:The pole at $s=-1$ exchanges the ground-state tachyon, the pole at $s=0$ exchanges massless states, and the positive integer poles exchange the infinite massive tower. The $t$-channel interpretation is the crossed version. Factorization means that each residue is a sum of products of couplings to intermediate physical states that couple to the chosen external particles. It need not expose every representation at that mass.

For example, the Gamma function recurrence and Gamma function residue at a nonpositive integer give the dimensionless Virasoro–Shapiro amplitude pole residueThe residue polynomial has degree $2n$, consistent with maximum spin $2n$ in the exchanged level. The amplitude also has the corresponding $u$-channel poles by crossing symmetry.

The massless fields can be treated as target backgrounds rather than separate asymptotic insertions. Write the target metric as $G_{mn}={\eta }_{mn}+h_{mn}$, introduce a Kalb–Ramond field $B_{mn}$, and a dilaton $\Phi$. In conventional Euclidean signs their string nonlinear sigma model action isHere ${ϵ}^{ab}$ is the antisymmetric tensor density. The three backgrounds correspond to the graviton, antisymmetric tensor and scalar states at closed-string level one. Expanding the vacuum functional in $h,B,\Phi$, then Fourier expanding the backgrounds, produces exactly their integrated string vertex operators. Its functional derivatives therefore generate the amplitudes with massless external strings. The connected vacuum functional organizes connected amplitudes; the spacetime effective action organizes the corresponding vertices after treating massless propagation consistently.

At momenta small compared with $1/\sqrt{{\alpha }^{\prime }}$, massive string propagators can be expanded in powers of momenta over their masses. The analytic part of the amplitudes consequently determines local higher-derivative interactions, ordered by powers of ${\alpha }^{\prime }$. Massless exchange poles are retained through propagation of the massless fields, rather than expanded into local contact terms. Up to field redefinitions, the leading massless-sector action in the string-frame metric isThe terms denoted $O({\alpha }^{\prime })$ contain additional derivatives, including curvature-squared terms in the bosonic theory. The expansion concerns the massless sector around the perturbative bosonic background; the tachyon instability remains and is not cured by omitting its field from this displayed action. Thus this is a formal perturbative effective description, not a claim of a stable bosonic vacuum.

Finally split the dilaton into a constant and its variation, $\Phi ={\Phi }_0+\varphi$. By the Gauss-Bonnet theorem, its dilaton Euler-characteristic weighting on a connected closed oriented surface follows fromDefine the string coupling by $g_s=e^{{\Phi }_0}$. A genus-$g$ path integral is weighted byThe sphere carries $g_s^{-2}$, the torus carries $g_s^0$, and each extra handle adds $g_s^2$. At higher genus one integrates over complex-structure moduli as well as insertion points, with the associated antighost insertions supplying the correct moduli measure. With canonically normalized external vertices an $N$-point genus-$g$ amplitude scales as $g_s^{2g-2+N}$.

Summing connected worldsheets of every genus of a surface yields the string-loop effective action expansionEach coefficient has its own low-energy ${\alpha }^{\prime }$ expansion. The two parameters have different roles: ${\alpha }^{\prime }$ resolves finite string size through higher derivatives, while $g_s^2$ counts additional string loops.

Use natural units and the Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$. A real scalar field assigns a real variable $\varphi (t,\mathbf{x})$ to each spatial point. Its Lagrangian density can be taken to beThe principle of stationary action gives the Euler-Lagrange equation $\square \varphi +V^{\prime }(\varphi )=0$. With $V(\varphi )=m^2{\varphi }^2/2$, this is the Klein-Gordon equation. An additional nonlinear part of $V$ describes interactions.

The canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$. The Legendre transform in mechanics gives the canonical Hamiltonian density of a real scalar fieldThe Hamiltonian equations $\dot{\varphi }=\pi$ and $\dot{\pi }={\nabla }^2\varphi -V^{\prime }(\varphi )$ recover the same field equation. In canonical quantization, the fields become operators satisfying the equal-time canonical commutation relationsA spatial lattice makes the analogy with many coupled quantum-mechanical coordinates precise. Each lattice field value is a coordinate, with its own conjugate momentum. The path integral is another representation of the same quantum evolution.

To see its origin, first consider one coordinate with $H=p^2/(2M)+V(q)$. Split a time interval into $N$ steps of length $\epsilon$ and insert position and momentum resolutions of the identity. The short-time kernel isMultiplying the kernels and integrating over intermediate positions gives the phase-space path integralThe endpoints of $q$ are fixed. The momentum integrals are Gaussian integrals; completing the square produces the configuration-space path integralAt finite slicing its normalization contains $(M/(2\pi i\epsilon ){)}^{N/2}{\prod }_{j=1}^{N-1}d{q}_j$. This fixes the composition law and the initial delta-function kernel. One sums over all paths, not merely solutions of the classical equation. Restoring $\hslash$ replaces the weight by $e^{iS/\hslash }$; stationary phase explains the emergence of classical trajectories.

For the field, use scalar field configuration eigenstates $|\phi ⟩$, satisfying $\hat{\varphi }(\mathbf{x})|\phi ⟩=\phi (\mathbf{x})|\phi ⟩$. Insert their completeness relations on every time slice. This givesThe endpoint field configurations are fixed. Integrating the Gaussian momentum variables leaves the scalar field path integral $\mathcal{N}\int \mathcal{D}\varphi e^{iS[\varphi ]}$. The functional measure means a regulated product over the field variables. A spacetime lattice or another ultraviolet cutoff makes this product finite before the continuum limit; interacting continuum calculations may require renormalization. The oscillatory Minkowski weight is an amplitude, not a positive probability density.

For vacuum expectation values, the boundaries must select the vacuum rather than arbitrary field configurations. Long imaginary-time evolution suppresses excited states: $e^{-TH}|n⟩=e^{-T{E}_n}|n⟩$, so after normalization only the lowest-energy component remains as $T\to \infty$. This is vacuum projection by imaginary time. The corresponding Feynman i-epsilon prescription in the real-time integral specifies the vacuum boundary conditions and the poles of the propagator. With $t=-i\tau$, the Euclidean path integral has the weight $e^{-S_E}$, whereIt is often a useful regulated starting point; analytic continuation returns the vacuum time-ordered quantities.

Introduce a classical source $J(x)$ and define the normalized vacuum generating functionalwith the same vacuum prescription in numerator and denominator. A functional derivative brings down $i\varphi (x)$. The order of the time slices makes the operator insertion time-ordered. Thus source differentiation inserts time-ordered field operators:The denominator removes vacuum diagrams and gives normalized expectation values. It is essential that these are time-ordered products; differentiating this vacuum functional does not directly give every possible operator ordering.

The free theory illustrates the method. Its quadratic kernel is $K=-\square -m^2$ with the vacuum pole prescription, and completing the square gives the Gaussian evaluation of a free scalar generating functionalTwo source derivatives give ${\Delta }_F=⟨\Omega |T\{\hat{\varphi }(x)\hat{\varphi }(y)\}|\Omega ⟩$. Higher derivatives give all pairings, the content of Wick theorem. For an interaction $S_{\mathrm{int}}$, one may use path-integral perturbation by source derivatives:Expanding this expression generates Feynman diagrams and their Wick contractions. The connected generating functional $W[J]=-i\log Z[J]$ retains connected contributions; in particular $\delta W/\delta J=⟨\varphi {⟩}_J$. These functionals turn the computation of field-operator expectations into source differentiation of an ordinary regulated integral.

Use units $\hslash =c=1$, the Minkowski metric ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$, and the Fourier transform convention $f(x)=\int d^{4}p{e}^{-ip\cdot x}\tilde{f}(p)/(2\pi {)}^4$. Here gamma matrices satisfy $\{{\gamma }^a,{\gamma }^b\}=2{\eta }^{ab}{I}_4$, and $\bar{\psi }={\psi }^{†}{\gamma }^0$ is the Dirac adjoint. A single Dirac field describes both Electrons and Positrons; these are the particle and antiparticle sectors of that field.

The free Maxwell Lagrangian and Dirac action, with a linear covariant gauge condition for the photon, giveBefore gauge fixing, the Maxwell Lagrangian has zero modes in field theory along $A_a\mapsto A_a+{\partial }_a\lambda$, so its quadratic operator cannot be inverted on all potentials. In the linear Lorenz gauge, the Faddeev-Popov determinant is $\det \square$. It is independent of $A$ and may be absorbed into the normalization; the corresponding Faddeev-Popov ghost fields have no interacting vertices in this Abelian linear gauge. A residual gauge symmetry is removed by the specified boundary conditions.

The free generating functional isThe Dirac field variables and their sources are independent Grassmann fields in this path integral. They anticommute; replacing them by ordinary commuting fields would give the wrong statistics and the wrong functional determinant. The bosonic Gaussian functional integral contributes an inverse square root of a determinant, and the Grassmann Gaussian integral contributes a determinant. If $K_A$ is the photon quadratic operator and $K_D=i{\gamma }^a{\partial }_a-m$, completing the square givesIn the last expression $D_F=i{K}_A^{-1}$ and $S_F=i{K}_D^{-1}$ are the photon propagator and Dirac propagator, with Feynman i-epsilon prescription. The products include the appropriate spacetime integrals and index contractions. Functional derivatives with respect to $J$, and consistently ordered left or right Grassmann derivatives with respect to the fermionic sources, generate the time ordering of the corresponding fields. In Feynman gauge, $\alpha =1$, the momentum-space two-point functions areThe Feynman i-epsilon prescription specifies vacuum boundary conditions rather than an arbitrary inverse of the differential operator.

The covariant photon propagator uses four potential components. Its operator formalism counterpart is Gupta-Bleuler quantization: impose $({\partial }_a{A}^a{)}^{(+)}|\mathrm{phys}⟩=0$ and take the Gupta-Bleuler null-state quotient. This leaves a positive physical state space with two transverse photon polarization vectors. The temporal and longitudinal oscillator components occur in intermediate covariant expressions; the Ward identity removes their dependence from physical amplitudes. The free Dirac field uses the canonical anticommutation relations, producing the same fermionic signs as its Grassmann fields in the path integral.

The operator-path-integral equivalence can be seen directly with a regulator. Divide time into small intervals and insert complete sets of field-coordinate states for bosons, and resolutions in fermionic coherent states for fermions. The bosonic matrix elements produce the phase-space factor $\exp i\sum [\pi \Delta \varphi -H\Delta t]$; integrating out the quadratic canonical momentum produces the bosonic action. The fermionic coherent state overlaps produce the first-order term $i{\psi }^{†}\dot{\psi }$ and the Berezin integral measure. Multiplying the short-time kernels recovers the path integral. Projecting the remote endpoints onto the Fock vacuum with an infinitesimal damping selects the same Feynman propagators as the operator formalism. Field insertions become time-ordered products under this construction. Conversely, their quadratic generating functional obeys the canonical free-field equations and has precisely the oscillator two-point functions, so its higher free correlators agree by the Wick theorem. This establishes the equivalence for the regulated free theory and order by order in the perturbation series.

To couple the matter field electromagnetically, promote its global phase symmetry to the local transformationThe gauge covariant derivative $D_a={\partial }_a+ie{A}_a$ obeys $D_a^{\prime }{\psi }^{\prime }=e^{-ie\lambda }{D}_a\psi$. Replacing ${\partial }_a$ by $D_a$ in the Dirac action therefore gives the invariant matter densityThe electromagnetic field tensor is unchanged by the local transformation. Thus the unfixed quantum electrodynamics action is gauge-invariant. The Dirac current $j^a=\bar{\psi }{\gamma }^a\psi$ is conserved by the matter equations, and the Electron and Positron excitations carry opposite charges. The added gauge fixing density selects a representative and is not itself invariant under arbitrary local transformations; it does not change gauge-invariant observables. At a free fermion vertex,For external on-shell Dirac spinors the corresponding current contraction vanishes. The quantum extension is the Ward identity, which makes physical amplitudes insensitive to adding a multiple of the photon momentum to its polarization vector. A gauge-compatible regularization preserves this vector-current identity.

With ${\mathcal{L}}_{\mathrm{int}}=-e\bar{\psi }{\gamma }^a\psi A_a$, expand $\exp (i\int {\mathcal{L}}_{\mathrm{int}})$ in powers of $e$. A term of order $n$ contains $n$ spacetime integrations and $1/n!$. Applying the Wick theorem pairs the free fields: an $A$-$A$ Wick contraction supplies a photon propagator, and a $\psi$-$\bar{\psi }$ Wick contraction supplies an oriented Dirac propagator. Each insertion supplies an interaction vertex. The permutations of contractions cancel the expansion factorials except for the Feynman-diagram symmetry factor. Interchanging Grassmann fields produces the fermionic sign, including a minus sign for every closed fermion loop. This is how Feynman diagrams arise from expectation values, rather than an extra dynamical assumption.

The normalized vacuum generating functional removes components with no external insertions. In the operator formalism, for an interacting-vacuum expectation value of an inserted product $\mathcal{O}$, the same cancellation appears asVacuum projection and the Feynman i-epsilon prescription are implicit. The linked-cluster theorem exponentiates all connected vacuum bubbles into the same factor in numerator and denominator, so it cancels. Normalize by $Z[0]$ to remove every vacuum component. This cancellation of vacuum bubbles still leaves products of disconnected diagrams that each contain external insertions. If only connected correlators are wanted, differentiate the connected generating functional $W=-i\log [Z[J,\bar{\eta },\eta ]/Z[0]]$.

The resulting momentum-space QED Feynman rules for the bare theory can be stated in Feynman gauge as follows:

The illustrated interaction vertex has incoming fermion momentum $p$, incoming photon momentum $k$, and outgoing fermion momentum $p^{\prime }=p+k$. Its solid-line arrows indicate fermion flow. The photon wavy line has no fermion-flow arrow. The propagators, the vertex $-ie{\gamma }^a$, momentum conservation, and the fermionic signs determine the perturbative amplitudes.

The statistical Hamiltonian in this question is already in thermal units, as indicated by $e^{-H}$ without an additional ${\beta }_T$. Define $W[h]=\log Z[h]$ and $F[h]=-W[h]$. The mean field and its linear response areThese are functional derivatives of the connected generating functional; the subtraction defines the connected correlation function. In the Landau approximation, neglect loop corrections and evaluate the field integral at a stable saddle $m_h$. It satisfies the Euler-Lagrange equationThis follows by varying the gradient term and integrating by parts, with periodic, decaying, or otherwise appropriate boundary conditions.

The two requested free energy functionals, following the source and Legendre conventions of the question, areIn the scalar-field source Legendre transform on the chosen stable branch, $h$ is chosen to produce $m$, and $\delta {\Gamma }_{\mathrm{L}}/\delta m=h$. Thus the imposed-source Helmholtz free energy and the fixed-order-parameter Gibbs free energy have the appropriate opposite source derivatives. These names are used in the question's magnetic-ensemble convention; the defining sign relation is what fixes the calculation. At leading Landau approximation there is no fluctuation-determinant term in $F_{\mathrm{L}}$.

Differentiate the saddle equation with respect to $h(y)$. The response obeysThereforeEquivalently, the inverse Hessian relation for a connected two-point function states that $G$ is the inverse kernel of ${\delta }^2{\Gamma }_{\mathrm{L}}/\delta m(x)\delta m(y)$. The factor $u_0/2$ follows from twice differentiating the quartic term $u_0{m}^4/4!$; it is not $u_0/6$. This tree-level connected response is obtained by varying the saddle. It does not require replacing the exact connected correlator by a product of the saddle values, which would incorrectly give zero.

For the requested single-momentum formula, assume a homogeneous source and a translationally invariant equilibrium phase, so $m_h(x)=m_0$ and $G(x,y)=G(x-y)$. SetThe Fourier transform with the printed positive sign sends $-{\nabla }^2$ to $q^2$, while the Dirac delta function transforms to one. Hence the Ornstein--Zernike correlation function hasThe inverse convention is $G(x)=\int d^{D}q(2\pi {)}^{-D}{e}^{-iq\cdot x}\tilde{G}(q)$. For general inhomogeneous $h(x)$, $K$ has nonconstant coefficients and $G$ depends separately on its two positions; the displayed momentum-diagonal formula then does not follow. The preceding differential equation still holds in that case.

At zero source and for $u_0>0$, the homogeneous saddle is $m_0=0$ for $r_0>0$, and $m_0^2=-6r_0/u_0$ on either selected stable ordered branch for $r_0<0$. Thus the Landau scalar correlation length isIn the ordered phase the negative bare quadratic coefficient is compensated by the positive curvature at the nonzero saddle. Retaining $m_0=0$ below the transition would instead give an unstable kernel and is not a physical correlation length. With the usual analytic thermal tuning $r_0=a_0(T-T_c)$, $a_0>0$, both branches diverge as $|T-T_c{|}^{-1/2}$, so the correlation-length critical exponent isThe high-temperature amplitude is $\sqrt{2}$ times the low-temperature amplitude for the same $a_0$. This statement is within Landau theory; fluctuations can change critical behavior outside the mean-field regime.

The $-\mu \psi$ term means that the generator is the grand potential, rather than just the physical energy. With boundary variations vanishing, the Hamiltonian functional isIts functional derivative isso the Gross–Pitaevskii equation is $i\hslash {\psi }_t=\delta H/\delta {\psi }^{\ast }$. Treating $\psi$ and ${\psi }^{\ast }$ as independent variables in this variation gives the factor $g$, not $2g$.

For repulsive interactions $g>0$ and $\mu >0$, complete the square in the local energy:The gradient contribution is nonnegative, so the homogeneous ground state has arbitrary constant complex argument andThe infinite-volume constant may be subtracted to define a relative grand potential; it does not affect the functional derivative. At fixed particle number one instead minimizes the physical energy $H+\mu N$, with $\mu$ as its Lagrange multiplier. For $g>0$, $\mu \le 0$ the grand-potential minimum over $n\ge 0$ is the vacuum $n=0$. For $g<0$ the displayed unconstrained functional is unbounded below; the positive-density stable uniform state assumed in the following parts requires repulsion.

The field $\varphi$ is a compositional order parameter in a binary fluid mixture. Its spatial integral is conserved in a closed sample, so advection and diffusion enter a continuity equation; the material derivative follows the incompressible fluid. The diffusive current $\mathbf{J}=-M\nabla \mu$, with positive order-parameter mobility $M$, moves material down the chemical potential gradient. The chemical potential of a composition field is the functional derivative of the symmetric quartic Landau free energyso $\mu =a\varphi +b{\varphi }^3-\kappa {\nabla }^2\varphi$. Positive $b,\kappa$ stabilize large amplitudes and penalize sharp interfaces; $a<0$ permits phase separation.

The momentum equation is the Navier-Stokes equation with dynamic viscosity $\eta$, constant mass density $\rho$, and the Korteweg force $-\varphi \nabla \mu$ arising from composition stresses. The pressure enforces incompressible flow, $\nabla \cdot \mathbf{v}=0$, and can absorb pure gradient contributions to that thermodynamic force. For example, $-\varphi \nabla \mu =\mu \nabla \varphi -\nabla (\varphi \mu )$ gives an equivalent force after redefining the pressure. These equations constitute deterministic Model H dynamics.

Their dissipative signs can be checked by the energy dissipation identity for Model H. Under periodic or no-flux, no-work boundary conditions, composition advection and capillary work cancel between $F$ and the kinetic energy, givingThus diffusion and viscosity dissipate energy, while the reversible advective coupling transfers it between composition and motion. Thermal noises are omitted at this deterministic hydrodynamic level.

Choose the asymptotic time translation: $N=1+O(r^{-1})$, $N^i=O(r^{-1})$, with their derivatives $O(r^{-2})$, in asymptotically Cartesian coordinates. Use the usual differentiable falloff $h_{ij}={\delta }_{ij}+O(r^{-1})$, ${\partial }_k{h}_{ij}=O(r^{-2})$, ${\pi }^{ij}=O(r^{-2})$, and analogous falloff for allowed variations; impose the standard additional conditions needed for finite bulk integrals. These variations need not fix the Arnowitt-Deser-Misner energy.