Chemical potential of a composition field Created 2026-10-05 Updated 2026-10-14
For a composition free energy , the local chemical potential is the functional derivative
At phase coexistence, exchanges between the bulk phases require equal chemical potentials. Mechanical equilibrium additionally requires equal pressure for a flat interface, or a Young–Laplace equation pressure jump for a curved interface.
Generating functional 2026-09-24
A generating functional is a path integral with a source coupled to a field. Its functional derivatives with respect to generate time-ordered correlation functions.
Kinetic coefficient Created 2026-10-05 Updated 2026-10-06
For a vacuum generating functional , normalization gives and . Functional derivatives of give normalized vacuum expectation values. The linked-cluster theorem shows that this ratio cancels Vacuum Feynman diagrams, while gives connected correlators. Components attached to external insertions are not removed merely by the normalization.
One-particle-irreducible vertex Created 2026-10-07 Updated 2026-10-14
A functional derivative of the quantum effective action gives an amputated one-particle-irreducible correlation function. Its perturbative graphs stay connected after cutting any single internal propagator. A free quadratic effective action has no vertices of order greater than two; interactions generate higher vertices.
Order-parameter mobility 2026-10-05
The positive order-parameter mobility relates the diffusive transport flux to the chemical potential gradient: . Combining this law with local conservation gives the Cahn--Hilliard equation when is the functional derivative of a composition free energy.
A linear operator on finite-dimensional vectors acts by . For a field the discrete index becomes a continuous coordinate, and a suitable operator is represented by an integral kernel, possibly a distribution:
The identity operator has kernel . If an inverse exists on the chosen domain and with the chosen boundary prescription, its kernel satisfies
with the analogous reversed identity. A differential operator is a distributional kernel: for example, . Its inverse is a Green function. Boundary conditions, a causal prescription and the treatment of zero modes are part of the inverse's definition, not consequences of formal matrix notation. An operator with an untreated zero mode does not have an inverse on the full field space.
The functional derivative is defined through the first variation: for arbitrary smooth test variation ,
In particular,
For a symmetric kernel, differentiating gives . The two terms from varying the two fields coincide; for a nonsymmetric kernel the symmetric part appears instead. Integration by parts is understood when is differential and the boundary variations vanish. A finite-volume, finite-mode regulator makes these statements ordinary matrix identities before a continuum limit.
For a finite-dimensional integral, . Expanding a polynomial interaction formally and interchanging each finite-order derivative with the regulated integral gives
where is the free Gaussian integral. This is a perturbation-series identity; it does not assert convergence of the full expansion for an arbitrary potential.
For the Minkowski functional integral, each field insertion is generated by
Expand , replace its field factors by these functional derivatives, and take the derivatives outside the free integral. The result is polynomial interactions as source differential operators:
If the field interaction term contains derivatives of the field, the corresponding spacetime derivatives act on the source-insertion operators as well. For the scalar polynomial interactions the operator in the exponent is
The free quadratic source exponential then supplies all Wick contractions, reproducing the Feynman rules. For normalized correlation functions, divide the resulting interacting functional by ; the interaction operator changes this denominator too, so it cannot simply be ignored.
Use a real scalar field and Minkowski action , with the vacuum selected by a Feynman propagator pole prescription or adiabatic vacuum boundary conditions. The normalized correlation functions are
Here is the interacting vacuum and denotes time ordering. These are the canonical vacuum expectation values, rather than unordered products of noncommuting fields. The functional integral formula assumes a regulated measure and the same vacuum prescription on both sides.
Couple an external source linearly to the scalar field and choose the normalized vacuum generating functional
Then functional derivatives insert the fields:
In the convention used below, the connected generating functional is
The logarithm removes products of disconnected components. For example the connected correlation function ; at three points one subtracts all three products and adds . Normalization has already removed vacuum bubbles.
For a bosonic quadratic action , completing the square in the generating functional gives . Two functional derivatives therefore give the Feynman propagator , with boundary conditions and a pole prescription included in the inverse. Zero modes must be removed or fixed before an inverse exists.
Take Minkowski metric , Fourier convention , and the Yang-Mills field strength . The interaction terms do not contribute to the quadratic action. Integration by parts in yields
The color kernel is diagonal. For , use the projectors and . They obey , , , so
Multiplication explicitly gives the identity on both transverse and longitudinal components. Continuing this inverse with the Feynman propagator prescription gives the gluon propagator
The second term has the causal double-pole prescription. Away from the pole this is the usual . For it becomes the Feynman-gauge adjoint propagator; gives a transverse covariant Landau gauge propagator.
The limit removes gauge fixing and has no finite inverse on the full field space. The longitudinal eigenvalue tends to zero and its inverse diverges. A contraction with conserved external currents kills the longitudinal term, but that does not define a gauge-unfixed Gaussian integral. This massless limit is not the massive-vector unitary-gauge limit: here it exposes the unfixed gauge orbit degeneracy.
Because , a periodic primitive exists. Fix its additive constant by . Integrating the conserved-field equation once in space gives : the possible integration constant vanishes on taking the spatial mean. Treat the spatial average as the inner product defining functional derivatives. Integration by parts, with no boundary terms under periodic boundary conditions, gives
Thus the conserved-field real amplitude equation is a gradient flow ascending , or descending , with positive diagonal mobilities. Along a classical solution,
For the mixed term put and . The Cauchy-Schwarz inequality and a completed square give
Consequently, for ,
The negative leading quartic coefficient proves the upper bound. The stronger pointwise square completion
is also useful: a superlevel set bounds the amplitude in , its spatial derivative in , and in , because a quartic polynomial with negative leading coefficient dominates its lower-order terms.
The dynamics dissipates toward the steady-state set rather than sustaining nonconstant periodic motion. Indeed increases to a finite limit, and the time integral of is finite. For global regular trajectories with the usual parabolic precompactness, every limit point is stationary, by the LaSalle invariance principle. A nonconstant periodic or recurrent orbit would force a strictly increasing to return to its original value and is impossible. Boundedness of a monotone functional alone does not prove convergence to one particular equilibrium; isolation of the limiting equilibrium or an additional convergence theorem supplies that stronger claim.
Use a Minkowski external source term , where is the integration variable, and normalize the vacuum functional:
The usual vacuum boundary prescription is implicit. This is the logarithmic convention for the connected generating functional. It makes the given quadratic expression consistent with . Another common convention is , with ; mixing the two conventions would change the factors of below.
Functional derivative generates the full and connected correlation functions:
These are time-ordered products in the vacuum, and by the stated assumption. Differentiating four times partitions the four insertions into connected blocks. All partitions with a singleton vanish, leaving
The shorthand arguments stand for spacetime points; this is a functional cumulant identity, not an assumption of a free field.
To construct the quantum effective action, define the source-dependent classical field
Locally invert this relation to obtain , assuming the response integral kernel is invertible with the chosen regulator and boundary prescription. Take its Legendre transform
The chain-rule terms proportional to cancel, so . A field-independent constant is fixed separately and is immaterial to the vertices. At zero external source the classical field is zero, and there is no linear term.
For the two-point vertex, differentiate the classical field once. If , then
The vertex coefficients in this problem are derivatives of , not of alone. Therefore
For the three-point vertex, the derivative of the connected two-point integral kernel at a general external source is
Combine this with and differentiate the inverse identity, . This gives the logarithmic-source inverse-kernel vertices relation
All kernels here are finally evaluated at zero external source. The connected three-point function is amputated by its three inverse quantum field theory propagators.
For the specified quadratic , the external source response is , hence . Substitution in the transform gives
For the standard free quantum field theory propagator , this is the usual quadratic free action, up to the vacuum prescription. The generating functional is Gaussian, so connected functions beyond second order vanish and there are no interaction vertices or nontrivial scattering amplitudes. Full higher even-point functions can still be nonzero: they are products of two-point functions, as the four-point formula shows.
Differentiate the product of Dirac delta functions defining , regarding and as independent density arguments. The random continuity equation is
For , the causal functional derivatives follow from the explicit strain solution and :
Consequently,
In the Furutsu–Novikov formula, the covariance samples the upper endpoint of the causal time integral. Symmetric regularization of the white noise gives half the mass there. The response vanishes as , while the strain response tends to one, so . Averaging the continuity equation yields the Fokker-Planck equation
There is no direct diffusion term in : its equation contains the time integral of colored noise, whereas the strain equation receives Gaussian white noise. The same result follows from the Itô diffusion with drift and diffusion vector .
With , the logarithm is the connected generating functional. Its derivative is the mean field. The specified Legendre transform satisfies , and the quantum effective action generates amputated one-particle-irreducible vertices by further functional derivatives. In this question is the negative of the commonly used Minkowski effective-action convention; consequently its tree term is . Its diagrammatic content is unchanged if the convention is used consistently.
Expand the action around a stationary background in the presence of the source. For , integration by parts gives
The linear term vanishes at the tree-level saddle, . Keeping the Gaussian fluctuation approximation, and fixing the field-independent normalization, gives
to one-loop order. The shift from the tree saddle to the exact mean field cancels from the tree Legendre transform at this order. Hence the one-loop scalar effective action in the question's convention is
The functional determinant has a Feynman i-epsilon prescription before Wick rotation. Higher loops come from the omitted cubic and higher fluctuation vertices.
To see the one-loop one-particle-irreducible Feynman diagrams, write . Its determinant expansion is
Each trace closes a single scalar loop, with propagators and background insertions . No internal edge is a bridge, so these graphs are one-particle irreducible. Expanding the background factors produces their external legs.
For a constant , set . The quantum effective potential is defined here by . A plane-wave basis diagonalizes , and the supplied Wick rotation converts its correction per unit volume to
Here the power of the renormalization scale keeps the continued integral in four-dimensional units. The Gamma function recurrence and expansion give
Consequently
The divergence is local, so subtract it with a potential counterterm. In the modified minimal subtraction scheme, remove its associated constants as well. The one-loop scalar effective potential becomes
up to an arbitrary field-independent vacuum-energy constant. Another allowed finite subtraction replaces the nonlogarithmic constant. Restoring a loop-counting parameter would multiply the displayed correction by . The logarithm assumes positive curvature ; a background with negative curvature requires an analytic continuation and can have an imaginary part signalling instability.
The quartic case is special because of quartic closure of scalar effective-potential counterterms. If has degree at most four, then has degree at most two and has degree at most four: its pole is removable by the original polynomial couplings and a constant. For example, gives
The last two terms reproduce the one-loop massive phi-fourth counterterms, and the first is a vacuum-energy subtraction. A degree- potential with instead generates degree , forcing new operators. This is the potential-level manifestation of power counting in quantum field theory and four-dimensional renormalizability, even though the finite effective potential itself need not remain polynomial.
Treat and its complex conjugate as independent variables in the functional derivative. Integrating the gradient variation by parts, with vanishing boundary terms, gives
The factors follow from and . Thus the cubic-quintic Gross–Pitaevskii equation is
Use decaying variations with the same complex bulk field one, so that integration by parts produces no variation boundary term. The functional derivatives of energy and renormalized momentum of a condensate are
For the second identity, variation of the two momentum terms and integration of the term each supply , so together they give . It follows that
is exactly the stationary Gross–Pitaevskii solitary wave equation. Thus each solitary wave is a constrained energy critical point with multiplier .
Let parametrize a differentiable family of these waves. Apply the critical-point identity to the variation , keeping the multiplier fixed at its value for that member of the family. Including the complex conjugate variation gives
Hence the energy–momentum slope of a solitary wave is
wherever . At a momentum turning point the parametrized identity remains the correct statement; a globally single-valued energy as a function of momentum is not needed. Differentiating as a composite function would include an extra term, so that composite derivative must not be set to zero.
Replace by , with vanishing on the boundary, and differentiate the action at zero. Integrate the derivative of by parts:
The coefficient of is the functional derivative .
For the electromagnetic action, . Antisymmetry combines these terms and integration by parts gives
Thus the inhomogeneous Maxwell equations are . The definition of also gives the homogeneous equations , equivalently . These are the two four-vector Maxwell systems.
Unrestricted variation of the alternative component-gradient action instead gives . This is equivalent to the first Maxwell system only when supplemented by the Lorenz gauge , because
The alternative action is not by itself gauge-invariant; the additional condition must be supplied consistently, not silently omitted.
Expand the tensor square: . Hence
Integrating this cross term by parts twice expresses its integral as a boundary term plus . Under Lorenz gauge it is therefore purely a boundary term. With fixed or vanishing boundary contributions, the two actions have the same variational equations once the additional gauge condition is included. This is the Lorenz-gauge reduction of the electromagnetic action.
Write and . The proposed gauge-field transformation law becomes . Consequently
Both and its gauge covariant derivative therefore transform as the same multiplet, with no derivative of the local parameter left over. The original algebraic invariance of the Lagrangian density now applies pointwise to , so this constructs the local symmetry.
Also transform by . Antisymmetry gives , so the source coupling is invariant. Let denote the functional derivative at a point. Vary the action functional and integrate the derivative of by parts, taking the local parameter to have compact support:
Since every is arbitrary, its coefficient vanishes. This proves the requested local identity, including the derivative sign and the order of the structure-constant indices.
For the generating functional, the integral of the first term is zero as a change-of-variable identity:
The divergence of this vector field in field space is proportional to . Thus an invariant regulated functional measure, with no gauge anomaly, leaves no Jacobian term. Inserting the classical identity into the integral gives
The same first-order differential operator annihilates the connected generating functional , because . Its source term is .
For the specified Legendre transform, vary . Since , the terms involving cancel, leaving
Antisymmetry also gives . Substitution in the identity for , followed by multiplication by , proves the background gauge invariance of a scalar effective action:
The derivatives of the quantum effective action are proper, amputated one-particle-irreducible vertices. At each loop order they receive contributions from one-particle-irreducible Feynman diagrams with the indicated external scalar legs, including tree vertices and the required counterterms. Diagrams which disconnect after cutting one internal line are not additional proper-vertex contributions. Translation invariance supplies the overall momentum-conserving delta function factored out in the definition. The signs of these vertices must follow the chosen overall convention for .
Differentiate the effective-action identity twice with respect to and and set . In an invariant zero-field vacuum this gives the Ward identity contact terms
The term proportional to vanishes. With the printed Fourier factors , integration by parts turns into . The first delta function combines with , and the second combines with . Using and symmetry of the scalar two-point kernel then yields
This derivation fixes the relative signs independently of any vertex rule.
There is an overall-sign inconsistency in the final printed example. At tree level, stationarity of the source integral gives , hence the specified transform gives . For the free scalar multiplet this implies, after factoring out the delta function,
The mixed result follows by expanding to first order in : its contribution is . The two scalar derivatives yield . Momentum conservation gives
so both sides of the Ward identity equal . Because the scalar theory is Gaussian for prescribed , its fluctuation determinant depends on but not on ; it cannot add further two-scalar vertices. Thus this quadratic verification is sufficient, not merely the first term of an omitted scalar-loop correction.
The printed mixed vertex instead has a minus sign. It satisfies the same Ward identity if one uses the opposite effective-action convention at tree level, because its two-point vertex is then as well. Either overall convention works, but the printed Legendre transform and printed mixed vertex cannot both be used unchanged. For example, with and a nonzero generator, they give opposite sides of the claimed identity.
Write and prescribe a deterministic entrance field . Substitution of the carrier into the Helmholtz equation gives
The usual weak-index paraxial model drops and linearizes , giving
All the mean-field formulas below use this model. If the term in is retained, its extra mean and fluctuations must also be included; Gaussian averaging of the linear random potential alone is then insufficient.
For a colored medium, taking expectations gives the exact propagation identity
Stationarity and Gaussianity alone do not close the last term. To make its dependence precise, let and let be the random paraxial propagator. Assuming a jointly Gaussian random field, rather than only Gaussian one-point marginals, the functional derivative of the solution is
The Furutsu–Novikov formula therefore yields
This is a propagation equation for the mean, with the remaining correlation displayed rather than silently replacing it by a product of means.
The general solution can be written explicitly as a Gaussian-averaged phase-screen product. Put , , and let
be the integral kernel of the Fresnel propagator. Successive free-propagation and random-phase steps, followed by Gaussian averaging, give the colored Gaussian paraxial mean propagator
For smooth colored media this is the time-slicing representation, with oscillatory integrals understood by the usual Fresnel regularization. At finite the exponential follows directly from the characteristic function of the jointly Gaussian phase sum. It exhibits why a general transverse covariance cannot be replaced by its value at zero separation along every diffraction path.
For completeness, the unlinearized paraxial model has and the corresponding mean equation includes . Its colored-medium solution uses the same time-slicing integral, replacing the Gaussian factor by
This follows by completing the square in the finite-dimensional Gaussian integral, with the determinant branch continuous from . It supplies the generic colored-medium mean when the term is retained. A square of ideal white noise is not defined by the covariance in part (b); that limit uses the linear random-potential model and its stated drift.
If transverse diffraction of the fluctuations is neglected, the simpler straight-ray approximation gives
Indeed and the variance of its phase is . If is constant transversely, this formula is exact with replaced by the free diffracted field . For example , a single standard normal variable constant throughout space, satisfies stationarity and unit variance but gives . Thus a constant-rate exponential attenuation is not a consequence of part (a)'s assumptions alone. It becomes the closed paraxial result under the longitudinal white-noise assumption in the next part.
Take a positive gradient coefficient , a positive squared mass, a finite-volume regulator and an ultraviolet cutoff. Use the dimensionless statistical action convention : inverse temperature is absorbed into the coefficients. If denotes dimensional energy instead, the exponent below gains a factor and the covariance gains . These conventions cannot be mixed.
Choose the Fourier transform convention
Integration of the quadratic gradient and source terms gives
Complete the square with . Translation of the regulated real Gaussian variables leaves their measure unchanged, so the Gaussian integral gives
Here is the source-free Gaussian partition function. For independent regulated real coordinates and a positive quadratic matrix , it is under Lebesgue normalization. The continuum determinant needs the stated regulator and a measure convention. Reality pairs the nonzero Fourier modes, so they must not be counted as two independent real fields.
A connected correlation function subtracts the product of the means. The logarithm generates its cumulants: differentiating twice with respect to the source, using the measure in the functional derivative convention, yields
The source changes the mean but not the covariance; the stated zero-source limit therefore follows immediately. Momentum conservation uses , appropriate for a real field without complex conjugation on either displayed factor.
For the continuum massive Gaussian field correlation tail, set and . Fourier inversion gives
The large-argument Modified Bessel function of the second kind has . Hence is proportional to at large distance and
At zero mass the exponential length diverges. For the critical continuum tail is algebraic, ; in lower dimensions infrared regulation needs additional care.
For a uniform source , perform a momentum-shell renormalization group step: separate slow modes and fast modes . In a Gaussian field theory they decouple; the uniform source couples only to the retained zero mode. Integrating fast modes contributes only a field-independent determinant to the free energy. Restore the cutoff by , and
The factors from the measure, gradient and two fields cancel in the gradient term. The mass term gains , and the linear source gains . Therefore
Here is the positive mass; for signed thermal coordinates it is the squared-mass relation that applies. The Gaussian fixed point has thermal scaling eigenvalue , source eigenvalue and zero anomalous dimension.
Under the formal homogeneous Gaussian scaling hypothesis, substitution of those eigenvalues gives
The last two values are independently visible in and , agreeing with Landau-Ginzburg theory. The specific-heat power follows directly from the determinant: after two thermal derivatives its singular contribution is proportional to
The change of variables proves this power because the rescaled integral converges at large exactly when . The two-dimensional free energy itself has a resonance, while its second derivative still has the indicated power.
The printed magnetization exponent is a formal Gaussian scaling index, not a proved ordered-phase law of the purely quadratic model. The ordered-phase obstruction in a Gaussian scalar model is explicit: the uniform field has energy density . At it tends to as , and the partition function diverges. For , the mean is and tends to zero as , so there is no spontaneously ordered branch. Thus no literal for can be established from this Statistical Hamiltonian alone. The value follows from the scaling dimension of the field if one assumes such a branch; for its nonpositive value further warns against that interpretation. Adding a stabilizing quartic interaction defines an ordered phase, but below four dimensions that interaction is relevant and generally leads to the Wilson-Fisher fixed point, not these Gaussian thermodynamic powers.
For the Gaussian specific-heat infrared threshold, at the same determinant integral grows like , so as a power exponent, with a logarithmic divergence. For , it has a finite cutoff-dependent limit; the leading heat capacity is a regular constant and has no divergent power, again described by in the requested convention. After subtraction of this background the Gaussian singular term can still have the negative power index ; it is not correct to erase that distinction. In a stable quartic Landau-Ginzburg theory above four dimensions, the quartic coefficient is a dangerously irrelevant coupling: it is needed for the ordered phase and contributes , giving the mean-field specific-heat jump and while invalidating naive hyperscaling relation. At four dimensions, marginal interactions can add further logarithmic corrections. The agreement with Landau theory is in the leading heat-capacity power; it does not assert the absence of logarithms or supply a missing stable Gaussian ordered phase.