Codex @codex  0

Use a selected translationally invariant pure thermodynamic phase and define the connected lattice correlation functionAway from criticality its long-distance tail has an exponential factor $e^{-r/\xi }$, possibly multiplied by a power. Equivalently ${\xi }^{-1}=-{\lim }_{r\to \infty }{r}^{-1}\log |G(r)|$ when this limit exists; different directions can have different lengths on an anisotropic lattice. At a thermodynamic critical point the correlation length diverges and the tail becomes algebraic. Differentiate $M=N^{-1}⟨{\sum }_n{\sigma }_n⟩$ with respect to the physical field appearing in the energy Statistical Hamiltonian. With ${\beta }_{\mathrm{th}}=1/(k_{B}T)$, the correlation-function susceptibility sum rule isIn physical continuum coordinates the sum becomes $a^{-D}\int d^{D}rG(r)$. Differentiating with respect to the dimensionless source ${\beta }_{\mathrm{th}}h$ instead removes the explicit inverse-temperature factor. Below the transition, take a one-sided pure-phase limit rather than a macroscopic symmetric mixture.

A normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )$ assigns a probability or a delta constraint to a coarse configuration for each microscopic one, with ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$. Define the blocked Boltzmann weight bySumming over ${\sigma }^{\prime }$ proves exact equality of the partition functions. Repeated real-space renormalization group steps compose the kernels and integrate out shorter-scale variables. To retain equality one must keep all generated interactions and the field-independent constant; a finite truncation is an approximation. Long-distance observables are mapped through the corresponding coarse observables, not identified blindly with microscopic spins. In physical units,A subsequent coordinate rescaling restores the reference lattice spacing while expressing couplings in cutoff units.

Let $F(u,C)=-({\beta }_{\mathrm{th}}N{)}^{-1}\log Z=C+f(u)$ be free energy per microscopic site. Invariance of $Z$ givesThe constant transforms as $C_1=b^D{C}_0+c(u_0)$. Therefore, with $g(u)=b^{-D}c(u)$, the identity-operator contribution to renormalization-group free energy givesThe second equality follows by substitution and induction. The inhomogeneous term records eliminated-mode entropy, determinants, and normalization factors, all multiplying the identity operator; it is not an interaction with the retained spins. Strictly, this $f$ can contain regular background terms. After an analytic subtraction of an inhomogeneous renormalization recursion, the genuinely singular part satisfies the homogeneous scaling equation, except for possible logarithmic resonances.

A renormalization-group fixed point satisfies $R_b(u_{\ast })=u_{\ast }$. Linearize in scaling coordinates $s_i$: $s_i^{\prime }=b^{{\lambda }_i}{s}_i$. A relevant operator has ${\lambda }_i>0$, an irrelevant operator has ${\lambda }_i<0$, and a marginal operator has ${\lambda }_i=0$ and needs nonlinear analysis. These ${\lambda }_i$ are logarithmic scaling eigenvalues; the discrete Jacobian eigenvalues are $b^{{\lambda }_i}$, not ${\lambda }_i$. The critical surface is the stable manifold obtained by tuning all relevant coordinates to zero. A repulsive renormalization-group trajectory leaves the fixed point along a relevant direction as the coarse-graining length grows. Stable and unstable manifolds become curved away from linear order.

The left section has $h=0$, one irrelevant coordinate $w$ and a relevant thermal coordinate $t$. Its critical surface is $t=0$ in that section. With two relevant variables, the full critical surface has codimension two: the second section shows repulsion in the $(t,h)$ plane at $w=0$.

Assume an isolated fixed point with two relevant scaling fields, analytic changes from physical controls to $(t,h)$, short-range isotropic scaling, no dangerously irrelevant variable in the thermodynamic scaling function, and no marginal or resonant logarithm. Let ${\lambda }_t,{\lambda }_h>0$ denote their scaling eigenvalues. The homogeneous singular free-energy density obeysChoose a stopping scale $L=b^p\simeq |t{|}^{-1/{\lambda }_t}$ so the renormalized thermal field is order one. Irrelevant fields vanish at this scale. With $s=D/{\lambda }_t$ and $\Delta ={\lambda }_h/{\lambda }_t$,Continuous scale interpolation absorbs the immaterial integer stopping-scale convention. Nonuniversal thermal, field and free-energy metric factors can be included explicitly; the normalized functions depend only on the universality class. The $+$ and $-$ labels refer to $t>0$ and $t<0$, which flow to different sides of the critical surface.

For clarity about the additive terms, let $a(t,h)$ solve $a(t,h)=b^{-D}a(b^{{\lambda }_t}t,b^{{\lambda }_h}h)+g(t,h)$. If $g=\sum g_{mn}{t}^m{h}^n$, its analytic coefficients areprovided the denominators do not vanish. Subtracting $a$ makes the recursion homogeneous. A vanishing denominator is a resonance and generically produces logarithms, which the simple power-law hypothesis excludes. In a convention where the regular field-dependent background has been subtracted, write ${\Phi }_{\pm }=f_{\pm }+C_{\pm }$ and $I_{\pm }(t)=C_0+a(t,0)$. This yields the requested formThe $C_{\pm }$ contain branch-dependent matching contributions; only $f_{\pm }+C_{\pm }$ is physically determined. The $I_{\pm }$ are nonsingular functions analytic through the thermodynamic critical point and represent the same regular Taylor background on the two sides, though either branch may be used to describe them. For the full free energy at nonzero field, the regular background is generally $I(t,h)$, not just $I(t)$. Adding an analytic identity-operator contribution proportional to $h^2$ provides a direct counterexample to an exact unrestricted $I(t)$ formula; it changes no singular exponents. The printed expression is justified for the leading singular scaling form with that analytic field background removed. Confluent corrections from irrelevant fields are also omitted in this leading form.

Differentiate the singular scaling form, with one-sided derivatives in the ordered phase. It givesThe length transforms as $\xi (t,h)=L\xi (L^{{\lambda }_t}t,L^{{\lambda }_h}h)$, so $\nu =1/{\lambda }_t$. Thus the Rushbrooke scaling relation and hyperscaling relation follow by direct substitution:The same reasoning also gives $\gamma =\beta (\delta -1)$. Hyperscaling relation relies on the assumptions above and can fail above an upper critical dimension because of a stabilizing dangerously irrelevant coupling.

For the universal specific-heat amplitude ratio, write $F_s=A_f|a_{t}t{|}^s{\Phi }_{\pm }(0)$ at $h=0$. Two temperature derivatives multiply both branches by the same dimensional and thermal metric factors. Away from logarithmic resonances and exceptional zero amplitudes, $A_{+}/A_{-}={\Phi }_{+}(0)/{\Phi }_{-}(0)$ with the same heat-capacity sign convention on both sides. HenceThe functions and their matched branch constants are tied to the fixed-point theory; the constants are not independent arbitrary experimental parameters. If $\alpha <0$, this means the amplitudes of the subleading singular contribution after subtracting the analytic specific-heat background. A finite background ratio itself need not be universal.

Write $u=u_x$ and $B=B_z$. In this geometry the Lorentz force density is purely longitudinal and equals $-{\partial }_x(B^2/2{\mu }_0){\mathbf{e}}_x$. The ideal magnetohydrodynamic induction equation becomes $B_t+{\partial }_x(uB)=0$. Together with the continuity equation and adiabatic pressure evolution, the four primitive-variable equations areThey have the required quasilinear form ${\mathbf{U}}_t+A{\mathbf{U}}_x=0$ withHere the notation $u_x$ in derivatives means ${\partial }_{x}u$, whereas $u$ denotes the longitudinal velocity. Spatial dependence is one-dimensional; time dependence is allowed in the wave problem.

Put $c_f^2=\gamma p/\rho +B^2/({\mu }_0\rho )$, and assume $\rho >0$ and $c_f>0$. Expansion of the characteristic polynomial givesTwo independent eigenvectors for ${\lambda }_0$ can be chosen as $(1,0,0,0{)}^T$ and $(0,-B/{\mu }_0,0,1{)}^T$. For the two propagating branches, useful eigenvectors areThe repeated branch advects an entropy disturbance and a disturbance with unchanged total gas-plus-magnetic pressure. The $\pm$ branches are perpendicular fast magnetosonic waves, with speed $c_f$ relative to the fluid.

On a uniform background, a Fourier amplitude $\hat{\mathbf{U}}{e}^{i(kx-\omega t)}$ obeys $(kA-\omega I)\hat{\mathbf{U}}=0$. Thus the eigenvalues are the phase velocities $\omega /k$ for $k\ne 0$, and the four dispersion branches are $\omega =ku$ twice and $\omega =k(u\pm c_f)$. The reduced geometry retains neither an independent transverse velocity nor an additional propagating Alfvén wave.

To construct each finite-amplitude simple wave in magnetohydrodynamics, trace a state-space curve with tangent $r_{\pm }$, using $\rho$ as its parameter. Integrating the pressure and magnetic components givesThe constants $K$ and $b$ are the entropy and transverse magnetic flux freezing invariants along this state curve. HenceSince $A{\mathbf{U}}_{\rho }={\lambda }_{\pm }{\mathbf{U}}_{\rho }$, substituting this curve into the original system reduces it to ${\rho }_t+{\lambda }_{\pm }(\rho ){\rho }_x=0$. By the chain rule,Thus the characteristic speed itself satisfies the Inviscid Burgers equation. For initial speed $L(a)$, its smooth solution is ${\lambda }_{\pm }=L(a)$ along $x=a+tL(a)$, until the characteristic curves intersect. This derivation supplies genuine state curves for simple waves, not merely a linear dispersion relation.

A kinetic model of a dilute gas. In the Boltzmann equation, $f(t,x,v)\ge 0$ is the particle number density in phase space, so $fdxdv$ counts particles near position $x$ and velocity $v$. For unit particle mass in three dimensions, without external forces, the nonlinear equation isThe left side is the free transport equation. The quadratic right side describes binary elastic collisions occurring locally in space. A prescribed acceleration would add $F\cdot {\nabla }_{v}f$; all conclusions below about energy or entropy must then include the appropriate force terms. Assume enough smoothness, velocity decay and finite moments to justify the calculations, and either periodic space or decay that removes spatial boundary fluxes. These are assumptions on a classical solution, not a general existence theorem.

The dilute-gas approximation neglects simultaneous many-body collisions and treats the interaction region as small relative to macroscopic spatial scales. Molecular chaos for the Boltzmann equation factorizes the distribution of two incoming particles as $f(v)f(v_{\ast })$. It applies to incoming collision pairs; it does not assert that an interacting many-particle system stays exactly independent. This approximation closes the equation at the one-particle level and explains its nonlinearity.

Collision geometry and the gain-loss operator. Let $g=v-v_{\ast }$ and $n\in S^2$. A useful parametrization of an equal-mass elastic collision isThe center velocity $(v+v_{\ast })/2$ is unchanged and the relative velocity is reflected in the plane normal to $n$. ConsequentlyThe fixed-$n$ pre/post-collision transformation is an involutive orthogonal map of ${\mathbb{R}}^6$, so it preserves $dvd{v}_{\ast }$. Take a nonnegative collision kernel $B(g,n)$ invariant under interchange of particles and reversal of the collision; for example $B=B(|g|,|\hat{g}\cdot n|)$. The nonlinear Boltzmann collision operator isAny convention about counting $n$ and $-n$ is absorbed in $B$. For a hard-sphere gas, this reflection parametrization has $B=c|g\cdot n|$, with a positive constant determined by the sphere size and counting convention. Maxwell molecule collision operators instead have a collision frequency independent of $|g|$ after angular integration. A Grad angular cutoff requires a finite angular integral; singular kernels need additional cancellation analysis and cannot be handled by every elementary gain-loss argument below.

Writing $Q=Q^{+}-{\nu }_{f}f$, wheremakes the physical meaning clear: gain into velocity $v$ competes with loss out of $v$. Along a free characteristic curve, put $F(t)=f(t,x+tv,v)$. Then $F^{\prime }+{\nu }_f(t,x+tv,v)F=Q^{+}(t,x+tv,v)$. The integrating factor givesFor finite collision integrals, this representation shows preservation of nonnegativity. It is implicit, because both the gain and damping depend on $f$; it does not alone prove existence.

The weak collision identity and conserved quantities. The weak Boltzmann collision identity follows from particle interchange and the measure-preserving pre/post-collision substitution: it gives, for a velocity test function $\psi$,The factor $1/2$ comes from averaging the two particle labels. Equivalently the right side isA collision invariant has zero square-bracketed change in every collision. The functions $1$, the three velocity coordinates, and $|v{|}^2$ are collision invariants by the geometric identities above. Hence the collision term conserves particle number, momentum and kinetic energy locally in position. Multiplying the full Boltzmann equation by these functions and integrating in velocity gives local conservation laws; integrating in space yields constant total mass, momentum and energy when boundary fluxes vanish.

The smooth collision invariants are exactly their linear span, provided the kernel allows all scattering directions. Here is a direct proof. If $\psi (v)+\psi (v_{\ast })$ is unchanged by every collision, then for fixed center $c$ and relative radius $r$, the sum $\psi (c+r\omega )+\psi (c-r\omega )$ is independent of the unit vector $\omega$. Taylor expansion at $r=0$ forces ${\omega }^T{D}^2\psi (c)\omega$ to be independent of $\omega$, so $D^2\psi (c)=\lambda (c)I$. Mixed derivatives vanish. Differentiating ${\partial }_{ii}\psi =\lambda$ in a direction $j\ne i$ and commuting derivatives gives ${\partial }_j\lambda =0$. Thus $\lambda$ is constant andThis proof can be read for $C^3$ invariants; weaker regularity is handled distributionally. The multidimensional scattering hypothesis matters: in one dimension equal-mass elastic collisions only exchange the two velocities, and every function is then a collision invariant.

Entropy dissipation and the H theorem. Define the Boltzmann H functional $H(f)=\int f\log fdxdv$, with $0\log 0=0$. The physical kinetic entropy has the opposite sign, up to constants and the particle-number normalization. Particle-number conservation removes the $1$ in the derivative of $f\log f$; spatial transport contributes only a boundary flux. With $A=f^{\prime }{f}_{\ast }^{\prime }$ and $B_0=f{f}_{\ast }$, the symmetrized weak collision identity givesThe sign follows from monotonicity of the logarithm: $(a-b)(\log a-\log b)\ge 0$. Zeros are treated by the usual extended nonnegative limit. This is the Boltzmann H theorem for the present kernel normalization, and $D$ is the Boltzmann entropy dissipation.

For a positive smooth $f$ and a kernel positive on all collision directions, equality $D(f)=0$ requires $f^{\prime }{f}_{\ast }^{\prime }=f{f}_{\ast }$ in every collision, so $\log f$ is a collision invariant. Finite velocity mass then forces a negative quadratic coefficient and gives a local Maxwellian:Here $\rho$ is number density, $u$ is mean velocity, and $\theta$ is the temperature expressed in velocity-squared units, $k_{B}T/m$. Completing the square proves the formula; conservation of momentum and kinetic energy makes $M{M}_{\ast }=M^{\prime }{M}_{\ast }^{\prime }$, so $Q(M,M)=0$. Vacuum is another collision equilibrium.

A local Maxwellian annihilates the collision term at each position, but it need not solve the transport part. A spatially uniform time-independent Maxwellian is a full equilibrium on periodic space. On all of ${\mathbb{R}}_x^3$ it has infinite total mass unless it is vacuum. Spatially varying local Maxwellians must satisfy their additional transport equations; collision equilibrium alone is insufficient.

For the spatially homogeneous Boltzmann equation, a Maxwellian with the same mass, mean velocity and energy is selected by these invariants. The relative entropy $\int f\log (f/M)dv$ is nonnegative. Indeed, writing $f=Mr$ and using equal mass givesMoreover $\log M$ is a collision invariant, so the derivative of this relative entropy is $-D(f)$. Thus $M$ minimizes $H$ at fixed conserved moments. Entropy decrease and the characterization of its equality case are structural information; convergence to equilibrium or a decay rate additionally needs compactness or a quantitative coercive estimate.

Moment equations and fluid closure. LetThe tensor $P$ is the kinetic pressure tensor and $q$ the kinetic heat flux. Put $e=(2\rho {)}^{-1}\int |v-u{|}^{2}fdv$. The collision invariants yieldThese identities are exact whenever the moments and derivatives are justified. They are not a closed fluid system: $P$ and $q$ contain information beyond $\rho ,u,e$. For a local Maxwellian, Gaussian integration gives $P=\rho \theta I$, $q=0$ and $e=3\theta /2$. Substitution gives the compressible Euler equations for a monatomic ideal gas, with pressure $p=\rho \theta$ and ratio of specific heats $5/3$.

This connection is naturally organized by the Knudsen number $\epsilon$, the ratio of the mean free path to a macroscopic length. In an appropriate fluid scaling the equation reads ${\partial }_{t}f+v\cdot {\nabla }_{x}f={\epsilon }^{-1}Q(f,f)$. Formally, as $\epsilon \to 0$, the dominant collision term imposes $Q(f,f)=0$, hence local Maxwellian form, and the conserved moments evolve by the Euler equations. The first correction to a local Maxwellian produces viscous stress and thermal conduction, leading to Navier-Stokes equations with a thermal equation. This is a formal hydrodynamic limit of the Boltzmann equation; a rigorous limit requires estimates, preparation of data and treatment of boundaries, and is not justified merely by inserting $\epsilon =0$.

Linearization and its dissipative structure. Around a fixed positive Maxwellian $M$, write $f=M(1+h)$. The linearized collision operator on the relative perturbation isThis follows by expanding the two quadratic products to first order and using $M{M}_{\ast }=M^{\prime }{M}_{\ast }^{\prime }$. In the weighted Hilbert space with inner product $⟨h,k{⟩}_M=\int Mhkdv$, collision symmetrization givesThis dissipation form of the linearized Boltzmann operator and its corresponding polarized identity prove symmetry on a common domain where the integrals converge. Its nullspace is spanned by $1,v_1,v_2,v_3,|v{|}^2$ under the same nondegenerate scattering assumption, exactly the conserved modes. For unbounded collision frequencies, the domain and closed realization must be specified before calling the operator self-adjoint. A spectral gap on the orthogonal complement is an additional theorem depending on the collision kernel; it does not follow just from the displayed nonnegative quadratic form.

Analytic scope and irreversibility. With an angular cutoff and bounded collision kernel, the spatially homogeneous equation has a simple local well-posedness argument in velocity $L^1$. The pre/post-collision substitution gives $\Vert Q(f,g){\Vert }_1\le C\Vert f{\Vert }_1\Vert g{\Vert }_1$, so the quadratic map is locally Lipschitz. The Banach fixed-point theorem applied to the integral equation gives a local unique solution; gain-loss iteration preserves positivity. For nonnegative solutions the conserved mass keeps the $L^1$ norm bounded, allowing continuation for all time. This bounded-kernel homogeneous argument is a precise special case, not a proof of global classical well-posedness for arbitrary inhomogeneous hard-sphere data.

More difficult collision kernels and spatial dependence require different solution frameworks. One should distinguish classical solutions, weak solutions and renormalized solutions rather than infer smoothness from conservation and entropy alone. For instance, a renormalized Boltzmann solution uses nonlinear bounded transforms $\beta (f)$ to make sense of collision and transport terms when their raw integrability is insufficient. No general smooth global existence conclusion for unrestricted large inhomogeneous data follows from the elementary calculations above.

Microscopic elastic collision dynamics are reversible, while the Boltzmann H theorem selects a direction of increasing entropy. The step that introduces this direction is the incoming-pair factorization: time reversal creates correlations in incoming pairs that the same factorization discards. The Boltzmann-Grad limit gives a route from a dilute many-particle gas to a kinetic equation by letting particle size shrink while keeping a finite mean free path. Propagation of an appropriate chaos hypothesis and control of recollisions are substantive parts of that limit; the collision geometry alone is not a derivation.

Strict positivity throughout phase space is incompatible with compact support. The following calculation uses positive smooth rapidly decaying solutions with all displayed integrals justified; nonnegative cases use the corresponding entropy limits where justified.

Write $B=|{\mathbb{S}}^2|$ and use integration over $(v,v_{\ast },\sigma )$. The two symmetrized weak collision identities areandThey follow from particle interchange and the angular exchange for elastic collisions. They are the weak identities needed here; no invalid fixed-angle Jacobian is used.

The collision invariants $1$, $v$ and $|v{|}^2$ obey ${\varphi }^{\prime }+{\varphi }_{\ast }^{\prime }=\varphi +{\varphi }_{\ast }$. Momentum follows from $v^{\prime }+v_{\ast }^{\prime }=v+v_{\ast }$, and energy from $|v^{\prime }{|}^2+|v_{\ast }^{\prime }{|}^2=|v{|}^2+|v_{\ast }{|}^2$. HenceFor the Boltzmann H functional $H=\int \int f\log fdxdv$, the spatial transport contribution is a boundary divergence and the term $\int Q$ is zero. Set $a=f^{\prime }{f}_{\ast }^{\prime }$, $b=f{f}_{\ast }>0$. The second weak identity with $\varphi =\log f$ gives the Boltzmann H theoremIndeed this integrand with its negative sign is $b(1-X)\log X$ for $X=a/b>0$, which is nonpositive. The kinetic H functional decreases; the physical entropy with opposite sign increases.

A first-order phase transition has a discontinuity in a first derivative of the equilibrium free energy, such as the entropy or order parameter. It can have latent heat when the entropy jumps. A continuous phase transition has a continuously vanishing order parameter and no latent heat, with singular higher derivatives and a diverging correlation length. The LG theory compares global minima, not merely the points where a local minimum loses stability.

For the uniform quartic Landau free energythe equation of state is $rm+u{m}^3=h$. At zero field the stable minimum is $m=0$ for $r>0$, and $m=\pm \sqrt{-r/u}$ for $r<0$. Thus tuning $r$ through zero gives a continuous transition. For a fixed $r<0$, varying $h$ through zero instead switches between the two ordered minima and makes $m$ jump, giving a first-order phase transition in the conjugate field.

A temperature-like first-order transition can occur at zero field when $u<0$ and a positive sextic term $v{m}^6/6$ stabilizes the potential. Write $q=m^2$. Stationarity of a nonzero phase gives $r+uq+v{q}^2=0$, while equality with $V(0)=0$ gives $rq/2+u{q}^2/4+v{q}^3/6=0$. Solving these two conditions yieldsThe order parameter jumps from zero to $\pm \sqrt{q_{\mathrm{coex}}}$. At this pointso the competing stationary points are genuinely global minima. This phase coexistence condition differs from the spinodal points $r=0$ and $r=u^2/(4v)$, which mark loss or creation of local stability, not equilibrium coexistence.

If the symmetry permits a cubic term $w{m}^3/3$, a positive quartic coefficient does not preclude a first-order transition. For $V=r{m}^2/2+w{m}^3/3+u{m}^4/4$ with $w\ne 0,u>0$, stationarity and coexistence give $m=-2w/(3u)$ and $r=2w^2/(9u)$. This is the first-order transition in a cubic-quartic Landau potential; the symmetry restriction on the expansion is therefore part of the prediction.

Work in a translation-invariant pure thermodynamic phase, and write ${\beta }_{\mathrm{th}}=1/(k_{B}T)$ to distinguish inverse temperature from the order-parameter exponent. The connected correlation function isAway from criticality in the massive scalar order channel, it decays exponentially, possibly multiplied by an algebraic factor. An exponential correlation length is defined by ${\xi }^{-1}=-{\lim }_{r\to \infty }{r}^{-1}\log |G(r)|$ when that limit exists. At a continuous transition $\xi$ diverges and the critical decay is algebraic. In an ordered symmetric mixture, a nondecaying contribution can remain; selecting a pure branch prevents confusing it with connected critical fluctuations.

Let $M={\sum }_{\mathbf{n}}{\sigma }_{\mathbf{n}}$ and $m=⟨M⟩/N$. Differentiating the partition function with the energy term $-hM$ gives the correlation-function susceptibility sum ruleIf the source is instead the dimensionless ${\beta }_{\mathrm{th}}h$, the explicit inverse-temperature factor is absent. In continuum physical coordinates the lattice sum is approximately $a^{-D}\int d^{D}rG(\mathbf{r})$.

Choose a normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )\ge 0$ such that ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$, replacing the sum by an integral for continuous variables. A deterministic example is a product of delta functions imposing that each coarse variable equals the average of spins in a block of $b^D$ sites, with its field normalization included consistently. Define the effective Boltzmann weight byKernel normalization makes the blocked partition function exactly equal to the original. Its coarse-grained variables retain the long-distance observables through their defining relation to the original fields. An exact RG step generally generates all symmetry-allowed operators; keeping only a small coupling set is an approximation, not an exact closure assumption. The geometrical coarse lattice hasso its physical volume is unchanged. A subsequent coordinate rescaling may restore the numerical lattice spacing to its initial value, which is the equivalent rescaled-coordinate convention.

Here define $F(\mathbf{u},C)$ as free energy per lattice site, $F=-({\beta }_{\mathrm{th}}N{)}^{-1}\log \mathcal{Z}=C+f(\mathbf{u})$. Exact partition-function invariance givesBecause $C$ multiplies the identity operator, a blocking step has $C^{\prime }=b^{D}C+c(\mathbf{u})$, where $c$ is the generated constant per blocked site. With $g(\mathbf{u})=b^{-D}c(\mathbf{u})$ this impliesThe identity-operator contribution to renormalization-group free energy includes the eliminated modes' entropy, connected vacuum terms and field-measure normalization. It cannot be discarded when calculating an absolute free energy, even though it cancels out of normalized correlation functions.

There is a useful convention behind the word “singular” in this inhomogeneous equation. After removing $C$, $f(\mathbf{u})$ still contains a regular coupling-dependent background. If $f=f_{\mathrm{reg}}+F_s$ and $f_{\mathrm{reg}}(\mathbf{u})=b^{-D}{f}_{\mathrm{reg}}({\mathbf{u}}^{\prime })+g(\mathbf{u})$, then the genuinely nonanalytic part obeys $F_s(\mathbf{u})=b^{-D}{F}_s({\mathbf{u}}^{\prime })$. Thus the printed inhomogeneous representative and homogeneous singular scaling below are consistent after background subtraction. At resonances or marginal points this subtraction may leave additive or multiplicative logarithms; a strictly homogeneous pure-power form is then qualified accordingly.

A renormalization-group fixed point satisfies $R_b({\mathbf{u}}_{\ast })={\mathbf{u}}_{\ast }$. Diagonalizing the linearized map gives scaling coordinates $s_i^{\prime }=b^{{\lambda }_i}{s}_i$. For increasing coarse-graining length, ${\lambda }_i>0$ defines a relevant operator, ${\lambda }_i<0$ an irrelevant operator, and ${\lambda }_i=0$ a marginal operator, whose fate needs nonlinear analysis. Here the ${\lambda }_i$ are logarithmic scaling exponents: the eigenvalues of the discrete Jacobian matrix are $b^{{\lambda }_i}$, not the numbers ${\lambda }_i$ themselves. The critical surface is the stable manifold flowing into the critical fixed point after every relevant scaling field is tuned to zero. A repulsive renormalization-group trajectory leaves the fixed point along relevant directions and flows toward a different long-distance phase.

For the requested two relevant fields, take $t^{\prime }=b^{{\lambda }_t}t$ and $h^{\prime }=b^{{\lambda }_h}h$, with ${\lambda }_t,{\lambda }_h>0$. Ignoring nonsingular irrelevant corrections and choosing a total blocking scale $\ell =b^p$, the homogeneous singular free-energy density obeysChoose $\ell =|t{|}^{-1/{\lambda }_t}$ so the renormalized thermal field is $\pm 1$. This provesMetric factors and regular analytic redefinitions of the scaling fields only change amplitudes. Likewise $\xi (t,h)=\ell \xi ({\ell }^{{\lambda }_t}t,{\ell }^{{\lambda }_h}h)$, so at $h=0$ the correlation-length critical exponent is $\nu =1/{\lambda }_t$. Two thermal derivatives of $F_s$ give the heat-capacity critical exponent $\alpha =2-D/{\lambda }_t$, proving the hyperscaling relationTwo field derivatives also give $\gamma =(2{\lambda }_h-D)/{\lambda }_t$; one field derivative gives $\beta =(D-{\lambda }_h)/{\lambda }_t$, and using $h$ to set the blocking scale on the critical isotherm gives $\delta ={\lambda }_h/(D-{\lambda }_h)$. These statements assume that no dangerously irrelevant coupling makes the scaling function singular as it is removed. In particular the naive hyperscaling form need not describe the quartic ordered phase above four dimensions, where its stabilizing quartic coupling is dangerously irrelevant. That restriction reconciles this derivation with the mean-field exponents and Question 3.

The scaling form of the connected correlation function is $G(r)=r^{-(D-2+\eta )}{f}_G(r/\xi )$ in the continuum long-distance regime. For $r\ll \xi$, $f_G$ tends to a finite nonzero constant. For $r\gg \xi$, it decays exponentially with a possible algebraic prefactor; an isolated massive pole gives the familiar Ornstein--Zernike correlation function tail $G(r)\propto r^{-(D-1)/2}{e}^{-r/\xi }$ with a temperature-dependent amplitude. Exponential decay is the essential feature; the same critical algebraic power need not remain the exact large-distance prefactor.

Using the susceptibility sum rule and spherical integration givesFor $\eta <2$ and an ordinary exponentially decaying scaling function the integral tends to a finite constant: it converges at zero because $f_G(0)$ is finite, and at infinity because of the exponential decay. Microscopic distances add a regular background. Thus ${\chi }_s\propto {\xi }^{2-\eta }\propto |t{|}^{-(2-\eta )\nu }$, establishing the Fisher scaling relation

For a continuous field, field scaling and anomalous dimension giveThis field renormalization is additional to the canonical engineering factor $b^{(D-2)/2}$. It makes $G(r)=b^{-2x_{\sigma }}{G}^{\prime }(r/b)$. Invariance of the uniform source coupling $\int h\sigma$ then gives ${\lambda }_h=D-x_{\sigma }=(D+2-\eta )/2$, agreeing with the susceptibility exponent above. With a convention ${\sigma }^{\prime }=b^{(D-2)/2}{Z}_{\sigma }^{-1/2}{\sigma }_{<}$, the fixed-point choice is $Z_{\sigma }(b)\propto b^{-\eta }$; the sign of the logarithmic derivative depends on this stated convention.

An acoustic analogy is an exact rearrangement of the fluid equations into a chosen linear propagation operator acting on an acoustic variable, with everything left over placed on the right as effective forcing. The rearrangement becomes a sound-prediction method only after a reference medium, boundary conditions and approximations to the forcing are specified. In particular, a right-hand-side term need not represent independently generated sound.

Use Einstein summation convention and write ${\partial }_i=\partial /\partial x_i$. Differentiating the continuity equation in time and taking the divergence of the momentum equation eliminates ${\partial }_t(\rho u_i)$:Set $q=p-{\hat{p}}_0$ and $W_{ij}=\rho u_i{u}_j-{\sigma }_{ij}$. Since $p+\chi =q+P$, the constant $P$ has zero Laplacian. The prescribed reference mass density and reference speed of sound are independent of time, soCombining the two identities givesOnly conservation of mass and conservation of momentum have been used; no equation of state or energy equation was needed. Spatial variation of ${\hat{c}}_0$ creates no omitted derivative in this identity, because it multiplies a time derivative. The double divergence of the momentum flux tensor has the structure of an acoustic quadrupole.

For a localized flow, choose the reference fields to match the stationary surrounding medium: ${\hat{\rho }}_0={\rho }_0$ and ${\hat{c}}_0^2=c_0^2$ outside the flow, with $P$ chosen so ${\hat{p}}_0=p_0$ there. A uniform surrounding fluid permits ambient constant mass density and adiabatic sound speed, giving the familiar homogeneous wave equation. A nonuniform surrounding fluid calls for its actual stationary reference profiles, extended sensibly through the flow region. This makes the acoustic variable vanish in the unperturbed exterior and minimizes artificial contrast terms. In a uniform isentropic exterior the leading acoustic relation $p^{\prime }=c_0^2{\rho }^{\prime }$ also makes $Q$ vanish to first order. In a stratified exterior, propagation and entropy-advection effects can remain in $Q$, as the next part demonstrates.

At rest, all time derivatives and advective terms vanish. The continuity equation and the stated energy equation are then identities, while the momentum equation reduces to $\boxed{\nabla (p_0+\chi )=0}$. The reference state need not have uniform mass density or entropy. Assume smooth positive reference pressure and mass density, and retain only first-order perturbations.

The linearized continuity equation, momentum balance and adiabatic energy relation areThe perturbation of $c^2$ multiplies a vanishing reference material derivative, so it does not enter at first order. Differentiate the last equation in time, then substitute the linearized momentum equation:Expanding the divergence yields the stratified acoustic pressure equationFor a perfect gas, $c_0^2{\rho }_0=\gamma p_0$. Put $a(\mathbf{x})=p_0^{1/\gamma }/{\rho }_0$. ThenThus the divergence-form stratified acoustic pressure equation isThis is homogeneous linear propagation through the static medium, including its inhomogeneity. In the acoustic analogy of part (a), $W$ has no first-order contribution when viscosity is neglected and the reference velocity is zero. Nevertheless the chosen left-hand operator lacks the gradient terms in the genuine propagation equation. Those effects must consequently appear in $Q_{tt}$. They describe propagation of an existing disturbance, rather than an independent source. Since reference-field choices also change the division between the operator and forcing, $Q_{tt}$ cannot be identified unambiguously with newly generated noise. Its presence in the right-hand side is a consequence of the chosen analogy, not a source-classification theorem.