Order parameter

An order parameter vanishes in a symmetric phase and becomes nonzero when that symmetry is spontaneously broken.

A field conjugate to a scalar order parameter $m$ contributes $-hm$ to its free-energy density. Differentiating the equilibrium free energy with respect to $h$ gives $-m$, with density or extensive normalization stated consistently. In a ferromagnet the conjugate field is magnetic; for a fluid-density order parameter it is related to a chemical-potential difference. It can explicitly break the symmetry whose spontaneous breaking defines the ordered phase.

A vector order parameter has several components transforming together as a vector under spatial rotations. A polar order parameter is an example. Its functional derivative $\delta F/\delta \mathbf{p}$ has one component for each component of the order parameter.

A conserved order parameter is a field whose spatial integral remains fixed under the dynamics. A local composition difference in a closed mixture is a standard example.

A topological defect is a localized obstruction to extending an ordered field smoothly and nonvanishingly through space. Its charge is determined by the homotopy class of the order parameter around the defect.

A spontaneous symmetry breaking chooses order-parameter directions independently in regions that have not been in causal contact. Mismatches around disconnected components, loops or enclosing spheres can trap domain walls, cosmic strings or magnetic monopoles. The correlation length cannot exceed the causal horizon, so one monopole per horizon volume estimates the smallest causal initial abundance up to geometrical and efficiency factors.

// Target: classical-field-theory-soliton.bigb

A domain wall is a codimension-one configuration interpolating between disconnected components of a vacuum manifold. In a one-dimensional system with broken discrete symmetry it is a point defect of finite energy, so positional entropy produces a nonzero density at every positive temperature.

For a Landau-Ginzburg theory with positive gradient stiffness $\kappa$, $r<0$ and quartic coupling $u>0$, the degenerate minima are $\pm M_0$, $M_0^2=|r|/u$. A planar domain wall solves $\kappa {\varphi }^{\prime \prime }=r\varphi +u{\varphi }^3$. Multiplying by ${\varphi }^{\prime }$ and using the limiting minima gives $\kappa ({\varphi }^{\prime }{)}^2/2=u({\varphi }^2-M_0^2{)}^2/4$. Separating this equation yields the displayed profile, and integrating its excess energy gives the positive surface tension $2\sqrt{2\kappa }|r{|}^{3/2}/(3u)$. Its area cost is subextensive compared with bulk volume, permitting phase-separated common-tangent construction states when the average order parameter is constrained.

A vortex is a topological defect around which the phase of a complex order parameter has nonzero integer winding. In two dimensions an isolated vortex in a phase-only model has energy proportional to the logarithm of the system size.

A quantum vortex has integer winding number of a condensate's complex argument. For atomic mass $m$, its superfluid velocity gives circulation $2\pi \hslash q/m$ for integer $q$. A straight vortex has a density-depleted core and, outside it, tangential velocity $q\hslash /(mr)$. Its long-range flow energy is proportional to $q^2\log (R/a)$, while its angular momentum has the sign of $q$.

In healing-length units a straight quantum vortex obeys $f^{\prime \prime }+f^{\prime }/r-{\mathcal{N}}^{2}f/r^2+(1-f^2)f=0$ and $f\to 1$. Inserting $f=1-a/r^2+\cdots$ makes its order-$r^{-2}$ balance $-{\mathcal{N}}^2+2a=0$, proving the tail. The physical healing length convention must be stated; the deficit of density $1-f^2$ has twice the leading coefficient of the amplitude deficit $1-f$.

For a centered quantum vortex in a two-dimensional Thomas–Fermi condensate, use an intermediate radius ${\xi }_{\mathrm{core}}\ll L\ll R$. The approximately uniform core contributes $\pi n_c{\hslash }^2{\mathcal{N}}^2[\log (L/\xi )+L_{0|\mathcal{N}|}]/m$. The outer superfluid velocity is $\hslash \mathcal{N}/(mr)$, so its kinetic contribution is $(\pi {\hslash }^2{\mathcal{N}}^2/m){\int }_L^R{n}_c(1-r^2/R^2)dr/r$. Integration cancels $L$ and produces the constant $-1/2$. The displayed result is the matched small-core approximation at fixed particle number, not a formula for arbitrarily large winding at fixed cloud size. The angular momentum is exactly $L_z=N\hslash \mathcal{N}$ for a normalized axisymmetric single-winding state.

For a cubic-quintic Gross–Pitaevskii equation with normalized bulk amplitude one, a straight quantum vortex of integer winding number $\mathcal{N}$ has radial equation $R^{\prime \prime }+R^{\prime }/r-{\mathcal{N}}^{2}R/r^2+(1-\alpha R^2-\beta R^4)R=0$. With $\alpha +\beta =1$, the derivative of the nonlinear term at $R=1$ is $-2(1+\beta )$. Balancing the order-$r^{-2}$ centrifugal term with the bulk amplitude restoring term proves the displayed tail for a stable nondegenerate bulk state, $1+\beta >0$. A vortex's amplitude tends to one; its full complex field tends to $e^{i\mathcal{N}\theta }$, not a direction-independent constant. At $\beta =-1$ the stated inverse-square expansion cannot cancel the centrifugal term.

If a continuous nonzero complex field tends uniformly to one on large circles, its values on a sufficiently large circle lie in a disk about one that excludes zero. That disk is contractible, so the field's winding number on that circle is zero. A lone unit-charge quantum vortex therefore cannot satisfy the literal boundary $\psi \to 1$. For a vortex it is the wave amplitude, rather than the entire complex field, that can approach one; $f(r)e^{i\theta }$ gives a simple illustration of direction-dependent far-field complex argument.

A quantum vortex becomes energetically favourable in a rotating reference frame when $\Delta E-\Omega \Delta L_z<0$. For a singly quantized vortex in a long radially Thomas–Fermi condensate, the cutoff-shell estimate isIntegrating the trapped number density against $u_{\theta }^2$ produces the logarithm and constant; integrating it against angular momentum per particle gives $\Delta L_z\simeq \pi H{n}_0\hslash R_r^2/2$. Resolved-core corrections can alter the nonlogarithmic constant. This criterion compares equilibrium rotating-frame energies and does not remove a dynamical vortex-entry barrier.

Single-valuedness of the condensate field gives total complex argument change $2\pi q$ on a closed loop avoiding zeros, where $q$ is an integer winding number. Integrating the superfluid velocity around that loop gives the displayed circulation. It is unchanged by smooth deformations of the loop that do not cross a quantum vortex.

A quench is a change of control parameters faster than the system can equilibrate. The subsequent relaxation can trap domains and topological defects.

A quench abruptly changes temperature or another control parameter and then studies the subsequent nonequilibrium relaxation of the order parameter.

A liquid crystal has orientational or partial translational order while retaining some fluid mobility.

An isotropic phase has no preferred spatial direction. In a liquid crystal, its orientational nematic order parameter vanishes, although the molecules themselves may be anisotropic.

A polar liquid crystal distinguishes head from tail and is described at leading order by a vector polar order parameter $\mathbf{p}$.

A nematic liquid crystal identifies the director $\mathbf{n}$ with $-\mathbf{n}$. Its rotational order is described by a symmetric traceless tensor such as $Q_{ij}=S(n_i{n}_j-{\delta }_{ij}/3)$.

The twisted nematic field effect switches a nematic liquid crystal between a twisted optically rotating state and a largely field-aligned state. With adiabatic optical following in a twisted nematic and crossed polarizers aligned to the two planar surface directors, a quarter-turn cell passes light at low voltage and suppresses it at high voltage. A rear reflector permits reflective operation with little power expended on light generation.

When the birefringence separation of optical propagation constants $\Delta k=2\pi \Delta n/\lambda$ is large compared with the spatial rate of rotation of the optical axes, light with linear polarization in a local eigenpolarization follows the twisting nematic director. Thickness alone is insufficient if birefringence vanishes. Strong field alignment along the propagation direction suppresses the transverse optical anisotropy and hence the quarter-turn polarization rotation.

The Fréedericksz transition is a field-driven orientational instability of a confined nematic liquid crystal. Elasticity and surface anchoring resist the electric field or magnetic torque; a lowest allowed distortion mode becomes unstable at a finite threshold. Positive dielectric anisotropy of a nematic permits a normal electric field to tilt a planar anchored director.

For planar strong anchoring and a quarter-turn base profile, ${\varphi }^{\prime }=\pi /(2L)$ and $\theta =0$. Expanding the one-dimensional twisted nematic energy through quadratic order gives $\Delta F=\frac{1}{2}\int [K_1{\theta }^{\prime 2}+\{(K_3-2K_2)(\pi /(2L){)}^2-\Delta ϵE^2/(4\pi )\}{\theta }^2]dz$. The lowest zero-endpoint mode is $\sin (\pi z/L)$. Its coefficient vanishes at the displayed threshold, with $V=EL$, provided the bracket and anisotropy are positive. Halving the electric coupling doubles $V_c^2$.

The dielectric response differs along and perpendicular to the nematic director. In Gaussian units the orientational electric contribution at prescribed field is $-\Delta ϵ(\mathbf{E}\cdot \mathbf{n}{)}^2/(8\pi )$; positive anisotropy favors alignment with the electric field. The electric term does not receive the additional one-half convention used in writing the Frank elastic terms.

A nematic director is a unit vector specifying the local unoriented axis of a nematic liquid crystal. Head-tail symmetry identifies $\mathbf{n}$ with $-\mathbf{n}$, so the order-parameter space is a Real projective space rather than a sphere.

For a unit nematic director, the bulk elastic density is $[K_1(\nabla \cdot \mathbf{n}{)}^2+K_2(\mathbf{n}\cdot \nabla \times \mathbf{n}{)}^2+K_3|\mathbf{n}\times \nabla \times \mathbf{n}{|}^2]/2$. The three positive moduli penalize nematic splay, nematic twist, and nematic bend. Boundary elastic terms, chirality, or varying nematic order require additional contributions.

For a laterally uniform nematic director $(\cos \theta \cos \varphi ,\cos \theta \sin \varphi ,\sin \theta )$, its divergence is $\cos \theta {\theta }^{\prime }$ and its twist invariant is $-{\cos }^2\theta {\varphi }^{\prime }$. Thus $f=K_1{\cos }^2\theta +K_3{\sin }^2\theta$ and $g=K_2{\cos }^4\theta +K_3{\sin }^2\theta {\cos }^2\theta$. A normal electric field contributes $-c{\sin }^2\theta$. The Euler-Lagrange equations for two fields are $(g{\varphi }^{\prime }{)}^{\prime }=0$ and $f{\theta }^{\prime \prime }+f_{\theta }{\theta }^{\prime 2}/2-g_{\theta }{\varphi }^{\prime 2}/2+2c\sin \theta \cos \theta =0$.

Nematic bend curves the integral lines of the nematic director. For a unit field, $(\mathbf{n}\cdot \nabla )\mathbf{n}=-\mathbf{n}\times (\nabla \times \mathbf{n})$, so the bend energy measures curvature along those lines.

Nematic twist rotates the nematic director around an axis perpendicular to it. The helix $\mathbf{n}=(\cos qz,\sin qz,0)$ has zero nematic splay and nematic bend, and twist density $K_2{q}^2/2$.

Nematic splay describes local spreading or convergence of the nematic director field. Its leading elastic cost is quadratic in the divergence and invariant under the head-tail identification $\mathbf{n}\sim -\mathbf{n}$.

A nematic disclination is a topological line defect of the unoriented director field. In a two-dimensional sample its point cross-section can have half-integer winding because $\mathbf{n}$ and $-\mathbf{n}$ represent the same state.

A neutral pair of opposite nematic disclinations can annihilate under dissipative relaxation, because the total topological charge of a two-dimensional nematic disclination is zero. Its far-field elastic energy increases as $2\pi K{\lambda }_0^2{q}^2\log (R/r_0)$ with separation $R$, so the pair has an attractive elastic force proportional to $-1/R$. As opposite cores meet, the nematic order parameter may vanish locally and the two windings disappear. Integer defects may first undergo energetic fission of an integer nematic disclination, leading to encounters of elementary half-charge defects. The actual paths and annihilation time require a dynamical law; the free energy alone does not specify core mobility. Pinning, imposed boundary winding or thermal defect production can obstruct a defect-free final state.

With constant amplitude outside a core, the planar nematic divergence-elasticity identity reduces the elastic density to $K{\lambda }_0^2|\nabla \chi {|}^2/2$. A topological charge of a two-dimensional nematic disclination $q$ has $\chi =q\theta +{\chi }_0$ and $|\nabla \chi {|}^2=q^2/r^2$. Integrating over an annulus gives the displayed energy, plus finite core and boundary terms. The outer cutoff is the sample size, compensating-defect separation, or another screening length. An isolated defect in an infinite ordered plane has divergent energy. A neutral pair $q,-q$ separated by $R$ instead has $E_{\mathrm{pair}}=2\pi K{\lambda }_0^2{q}^2\log (R/r_0)$ plus core terms, since its far-field gradients cancel.

In a strictly planar nematic liquid crystal, a charge-one nematic disclination is a winding-two map into the circle ${\mathbb{RP}}^1$. Its nonzero fundamental group class prevents complete unwinding while the surrounding loop stays fixed, but permits splitting into two charge-$1/2$ defects. For separation $r_0\ll d\ll L$, the two near fields contribute $\pi K{\lambda }_0^2(1/4+1/4)\log (d/r_0)$ while their common charge-one far field contributes $\pi K{\lambda }_0^2\log (L/d)$. Compared with the unsplit defect, the change isIt is negative at sufficiently large scale separation. Thus integer-defect fission is energetically favorable despite conservation of total topological charge. This leading-log argument proves failure of global energetic stability, not absence of kinetic barriers or pinning. If escape into the third dimension is allowed, an integer winding can also become topologically trivial in ${\mathbb{RP}}^2$.

For a planar tangential charge-one nematic disclination, $\nabla \cdot Q=-({\lambda }^{\prime }/2+\lambda /r){\mathbf{e}}_r$. This follows by differentiating the tensor dyad, using $\nabla \cdot {\mathbf{e}}_{\theta }=0$ and $({\mathbf{e}}_{\theta }\cdot \nabla ){\mathbf{e}}_{\theta }=-{\mathbf{e}}_r/r$. With bulk density $a{\lambda }^2/2+b{\lambda }^4/4$ and divergence elastic density $K({\lambda }^{\prime }/2+\lambda /r{)}^2/2$, radial variation of $2\pi \int rfdr$ givesFor $a<0$, a regular melted core has $\lambda =O(r^2)$ and an ordered far field has ${\lambda }_0=\sqrt{-a/b}$. Linearizing the bulk derivative gives $\lambda /{\lambda }_0=1-K/(2|a|r^2)+O(r^{-4})$. The core scale is of order $\sqrt{K/|a|}$; the amplitude disturbance has an integrable algebraic tail, rather than exactly compact support.

For a continuous director angle around a closed counterclockwise circuit, the enclosed topological charge is $q=\Delta \theta /(2\pi )$. Since $\theta$ is defined modulo $\pi$, elementary nematic disclinations can have $q=\pm 1/2$.

A two-dimensional director winding can sometimes unwind when the director may tilt out of the plane. For a three-dimensional nematic the fundamental group of the order-parameter space is ${\pi }_1({\mathbb{RP}}^2)=\mathbb{Z}/2\mathbb{Z}$, so two disclination lines are topologically equivalent to none.

In the one-elastic-constant approximation, slow director distortions have a quadratic gradient energy with one modulus. For a planar angle field it reduces to $F=(K/2)\int |\nabla \theta {|}^2$ up to the convention-dependent effective modulus.

For a planar nematic director $\mathbf{n}=(\cos \chi ,\sin \chi )$ and constant amplitude ${\lambda }_0$, the nematic order parameter isDifferentiation gives $\nabla \cdot Q={\lambda }_0(-\sin 2\chi {\chi }_x+\cos 2\chi {\chi }_y,\cos 2\chi {\chi }_x+\sin 2\chi {\chi }_y)$. The coefficient matrix is orthogonal, proving the identity. Hence $K|\nabla \cdot Q{|}^2/2$ is the angle-gradient energy with effective modulus $K{\lambda }_0^2$. The full tensor-gradient norm instead obeys ${\sum }_{ijk}({\partial }_k{Q}_{ij}{)}^2=2{\lambda }_0^2|\nabla \chi {|}^2$ in this normalization; the two local densities must not be confused.

A uniaxial nematic order parameter has one distinguished nematic director and equal eigenvalues in its perpendicular plane. In the displayed normalization the eigenvalues are $\lambda ,-\lambda /2,-\lambda /2$, giving $\mathrm{Tr}{Q}^2=3{\lambda }^2/2$ and $\mathrm{Tr}{Q}^3=3{\lambda }^3/4$. Positive and negative $\lambda$ describe prolate and oblate orientational anisotropy, respectively; reversing the nematic director leaves the tensor unchanged.

Under rotation by a rotation matrix, a symmetric second-rank tensor transforms by conjugation, so its traces of powers are invariant. For a three-dimensional traceless second-rank tensor, the Cayley-Hamilton theorem implies $\mathrm{Tr}{Q}^4=(\mathrm{Tr}{Q}^2{)}^2/2$. In two dimensions its eigenvalues are $s,-s$, so $\mathrm{Tr}{Q}^3=0$. These identities organize the Landau-de Gennes free energy.

The Landau-de Gennes free energy is a rotational symmetry preserving Landau free energy of the nematic order parameter. One three-dimensional normalization through fourth order isThe cubic invariant generically produces a first-order phase transition. Its absence in two dimensions is an algebraic identity, not a guarantee that a mean-field approximation captures all fluctuations.

The leading scalar coupling of an electric field to the nematic order parameter is quadratic in the field and linear in the tensor. For uniaxial nematic order, it is $-3\chi \lambda E^2({\cos }^2\vartheta -1/3)/2$. When $\chi \lambda >0$, minimization aligns the nematic director with the field axis and adds $-\chi E^2\lambda$ to the scalar Landau free energy.

For a scalar Landau-de Gennes free energy with $\bar{b}>0$, $\bar{c}<0$, and coupling $-\chi E^2\lambda$ with $\chi >0$, the critical endpoint of a quartic Landau free energy is reached at the displayed field strength and ${\bar{a}}_c=3{\bar{c}}^2/(8\bar{b})$, ${\lambda }_c=-\bar{c}/(4\bar{b})$. It terminates the field-biased first-order phase transition when both field and quadratic coefficient can be tuned. This is a quartic small-field mean-field approximation; higher terms may alter its location.

Nonconserved relaxational order-parameter dynamics has ${\partial }_t\psi =-\Gamma \delta F/\delta \psi$ plus thermal noise. Its zero-wavenumber mode may decay.

For quadratic free energy with positive $a,\kappa$ and kinetic coefficient $\Gamma$, a deterministic initial Dirac delta distribution evolves under nonconserved order-parameter dynamics to the meanTaking the expectation removes the centered Gaussian white noise; the resulting diffusion-reaction equation is solved by the heat kernel. The integrated mean is $A{e}^{-\Gamma at}$, so spreading does not conserve the initial excess. The connected Fourier mode covariance grows from zero to $k_{B}T/(a+\kappa q^2)$ as $1-e^{-2\Gamma (a+\kappa q^2)t}$ by the Itô isometry. A realization retains equilibrium fluctuations after its mean has decayed. Pointwise continuum white-noise quantities require a coarse-graining cutoff.

For nonconserved order-parameter dynamics, the positive kinetic coefficient $\Gamma$ converts the functional derivative of the free energy into a relaxation rate through $\dot{p}=-\Gamma \delta F/\delta p$. Microscopic reversibility fixes the accompanying Gaussian white noise variance by the Model A fluctuation-dissipation relation.

For additive Gaussian white noise of variance ${\sigma }^2$, a Model A trajectory has path weight proportional toComparing a path with its time reverse yields the free-energy change and hence the thermal fluctuation-dissipation relation.

For additive Gaussian noise, a stochastic trajectory has probability proportional to the exponential of minus the squared noise history required to generate it. Rewriting that history in terms of the trajectory gives an Onsager--Machlup action; comparing it with the time-reversed action yields a path-probability ratio.

The Onsager–Machlup functional is the action in the exponential of an additive Gaussian stochastic path probability. For $\dot{x}=b(x)+\sigma \Lambda$, its elementary discretization contains $\int |\dot{x}-b(x){|}^2/(2{\sigma }^2)dt$, with convention-dependent Jacobian terms for multiplicative or state-dependent drift.

For additive Gaussian white noise and a time-even order parameter, a midpoint discretization of the Onsager–Machlup functional evaluates drift derivatives at the midpoint configurations. Reversing a trajectory visits the same midpoints in reverse order, so the noise-to-path Jacobian determinant is invariant under reversal. It may depend on the trajectory, and must be distinguished from the constant normalization of the noise measure.

For nonconserved order-parameter dynamics with $\dot{p}=-\Gamma \mu [p]+f$, the factor at one time step, apart from a path-independent power of the step size, isHere the fields have first been restricted to a finite spatial grid, and ${\mathcal{H}}_n$ is the Hessian matrix of the free energy.

For mixed conserved and nonconserved order-parameter dynamics, the required noises for a joint configuration and flux trajectory arewhere $\mathsf{D},\mathsf{G}$ represent the divergence and gradient, and $\Delta =\mathsf{DG}$ represents the Laplacian. The Jacobian matrix isTaking its Schur complement gives the joint factorFor periodic boundary conditions, discretizing the two spatial operators compatibly gives $\mathsf{D}=-{\mathsf{G}}^{\mathsf{T}}$, so $\Gamma I-M\Delta$ is the positive relaxation operator when $\Gamma ,M>0$. Reversal leaves ${\mathcal{H}}_n$ unchanged and reverses the sign of the flux. The product of the joint factors is consequently the same for both histories, justifying its cancellation from the joint path probability of an order parameter and its flux. The functional chain rule used in the action ratio holds in the continuum limit of this midpoint convention.

For nonconserved order-parameter dynamics, microscopic reversibility with equilibrium weight proportional to $e^{-F/(k_{B}T)}$ requiresIndeed the Onsager--Machlup path probability gives $\log (P_F/P_B)=-2\Gamma \int \dot{p}\delta F/\delta pdt/{\sigma }^2=-2\Gamma (F_2-F_1)/{\sigma }^2$. Equating this to the detailed balance value proves the relation. The functional chain rule is understood with the same midpoint convention as the path action.

After a sudden change from quadratic kernel $K_I(q)$ to $K_F(q)$, each Fourier mode is an Ornstein-Uhlenbeck process with decay rate $r(q)=\Gamma K_F(q)$. Its equal-time covariance interpolates exponentially between $k_{B}T/K_I(q)$ and $k_{B}T/K_F(q)$.

Conserved relaxational order-parameter dynamics has ${\partial }_t\varphi =M{\nabla }^2(\delta F/\delta \varphi )$ plus conserved thermal noise. Its relaxation rate vanishes as $q^2$ at small wavenumber.

The positive order-parameter mobility $M$ relates the diffusive transport flux to the chemical potential gradient: $\mathbf{J}=-M\nabla \mu$. Combining this law with local conservation gives the Cahn--Hilliard equation when $\mu$ is the functional derivative of a composition free energy.

For the flux $W_{ij}=-M{\partial }_i{\mu }_j+N_{ij}$ of a vector order parameter, microscopic reversibility requires independent Gaussian white noise components with covarianceThe current changes sign under time reversal. The resulting path-probability ratio contributes $-(k_{B}T{)}^{-1}\int W_{ij}{\partial }_i{\mu }_j$ to its logarithm; integration by parts converts this to the free-energy loss from the conserved dynamics.

A hydrodynamic mode is a long-wavelength collective mode whose decay rate tends to zero because of a conservation law. Coupling to fast nonconserved variables renormalizes its diffusivity and eigenvector.

With constant positive kinetic coefficient $\Gamma$ and mobility $M$, two independent relaxation channels giveThe deterministic mobility operator is $\Gamma -M{\nabla }^2$. Each channel separately satisfies microscopic reversibility when its noise obeys the corresponding Model A fluctuation-dissipation relation or fluctuation-dissipation relation for a conserved flux.

Put $u_j={\dot{p}}_j+{\partial }_i{W}_{ij}$. The additive independent noises give the forward actionFor a time-even order parameter and time-odd flux, the backward action replaces $u,W$ by $-u,-W$ while keeping $\mu$ fixed on the corresponding configurations. Common normalization and time-reversal invariance of a path Jacobian then giveAt the two thermal noise strengths, periodic boundary conditions cancel the two spatial terms by integration by parts, leaving $-(F_2-F_1)/(k_{B}T)$. Thus the joint dynamics satisfies microscopic reversibility.

A compositional order parameter is a local concentration difference in a mixture. In a closed system its spatial integral is fixed, so its dynamics must take continuity-equation form.

A binary fluid mixture contains two molecular species whose local concentration difference can undergo mixing or phase separation.

If a nonconserved surfactant polarization has local energy $\nu |\mathbf{p}{|}^2/2+\lambda \mathbf{p}\cdot \nabla \varphi$, eliminating $\mathbf{p}$ changes the composition-gradient coefficient from $\kappa$ to $\tilde{\kappa }=\kappa -{\lambda }^2/\nu$. This expresses the surfactant's reduction of interfacial energy.

A microemulsion is a thermodynamically stable mixture of immiscible fluids and surfactant with structure on a mesoscopic length scale. In a gradient theory, a negative square-gradient coefficient stabilized by a positive fourth-gradient term favors fluctuations at nonzero wavevector.

The Cahn--Hilliard equation transports a conserved composition field down gradients of its chemical potential: ${\partial }_t\varphi =\nabla \cdot (M\nabla \mu )$ with $\mu =\delta F/\delta \varphi$. Advection adds $\mathbf{v}\cdot \nabla \varphi$ to the left-hand side.

In Ostwald ripening, the higher interfacial chemical potential of small droplets drives diffusion toward larger droplets, increasing the characteristic domain size. The Gibbs--Thomson relation supplies the curvature-dependent surface value, while conserved order-parameter dynamics controls transport through the surrounding phase.

For a three-dimensional droplet with concentration jump $\Delta \varphi$, constant mobility $M$, and quasi-static exterior chemical potential,The first formula solves the exterior Laplace equation. The second follows from conservation: the total outward flux $4\pi MR({\mu }_s-{\mu }_{\infty })$ removes excess composition at the rate $-4\pi R^2\Delta \varphi \dot{R}$. For a symmetric mixture, $\Delta \varphi =2{\varphi }_B$ and the Gibbs--Thomson relation gives ${\mu }_s=\gamma /({\varphi }_{B}R)$.

A quasi-static chemical potential outside a fixed-shape droplet solves the same Laplace equation as an electrostatic potential outside an equipotential conductor. If $\mu ={\mu }_s$ on the droplet and zero on the surrounding reservoirs, its outward diffusive current iswhere $C$ is the corresponding electrostatic capacitance in a medium of permittivity $ϵ$. The arbitrary $ϵ$ cancels against the capacitance's proportionality to $ϵ$. Thus electrostatic geometry directly determines evaporation rates.

For a spherical droplet at fixed center height $h>R$ above a flat equilibrium reservoir, the diffusion-capacitance analogy and sphere-plane capacitance giveAt large separation, $c(R/h)=1+R/(2h)+O((R/h{)}^2)$, soThe nearby reservoir increases the total flux and shortens the lifetime. The spherical-shape assumption controls the geometry; the local current density is not uniform over the surface.

Model H couples a conserved scalar composition field to an incompressible momentum density. It combines the advective Cahn--Hilliard equation with the Navier-Stokes equation and the thermodynamic force density $-\varphi \nabla \mu$.

For constant positive order-parameter mobility $M$, dynamic viscosity $\eta$ and mass density $\rho$, deterministic Model H dynamics dissipates the sum of composition free energy and kinetic energy:Assume smooth fields and periodic or closed, no-work boundaries. Using the chemical potential of a composition field, integration by parts and incompressible flow gives $\dot{F}=\int \varphi \mathbf{v}\cdot \nabla \mu -M\int |\nabla \mu {|}^2$. Dotting the momentum equation with velocity gives ${\dot{E}}_{\mathrm{kin}}=-\int \varphi \mathbf{v}\cdot \nabla \mu -\eta \int |\nabla \mathbf{v}{|}^2$. Pressure and convective work vanish under the boundary assumptions; the two capillary transfer terms cancel. This proves both the sign of the Korteweg force and its reversible exchange of energy with the composition field.

The Korteweg force is the capillary body-force density generated by composition gradients in a diffuse-interface fluid. Up to a pressure gradient it can be written $-\varphi \nabla \mu$, where $\mu$ is the chemical potential.

In bicontinuous phase separation, both fluid phases form interpenetrating connected domains. During late-stage coarsening their morphology is often statistically self-similar after lengths are divided by one characteristic domain size $L(t)$.

The dynamical-scaling hypothesis replaces all macroscopic lengths by one scale $L(t)$ and the characteristic velocity by $\dot{L}$. Then inertial, viscous, and capillary force densities scale respectively as $\rho \ddot{L}+\rho {\dot{L}}^2/L$, $\eta \dot{L}/L^2$, and $\sigma /L^2$.

Replacing domain-scale viscous stresses by a local linear drag density $-\bar{\eta }\mathbf{v}$ gives the scales $L_1=(\sigma \rho /{\bar{\eta }}^2{)}^{1/3}$ and $t_1=\rho /\bar{\eta }$. The single-scale force model isFor a regular power-law asymptotic $g\sim C{u}^y$, $g^{\prime }\sim Cy{u}^{y-1}$, $g^{\prime \prime }=Cy(y-1)u^{y-2}+o(u^{y-2})$, $y>0$, the two inertial-to-drag ratios scale as $(y-1)/u$ and $y/u$. Thus drag dominates inertia at late time, and drag-capillary balance integrates to $g^3\sim 3(c_s/c_d)u$. This proves the displayed $1/3$ exponent and its independence of mass density. The scaling assumes a homogenized stationary network and dominant advective transport. Diffusive Ostwald ripening can have the same exponent, so its omission must be justified separately; order-parameter mobility can otherwise enter the prefactor.

Under single-length dynamical scaling of binary-fluid coarsening, mass density $\rho$, dynamic viscosity $\eta$ and surface tension $\sigma$ determine the displayed unique length and time scales by dimensional analysis. Writing $L=L_{0}f(t/t_0)$ reduces the force-density estimate to $\alpha f^{\prime \prime }+\beta f^{\prime 2}/f=\gamma f^{\prime }/f^2+\delta /f^2$. Resistive and driving terms need opposite signed geometric coefficients. Viscous hydrodynamic coarsening has $f\sim u$ and inertial hydrodynamic coarsening has $f\sim u^{2/3}$. These follow respectively from density-independent and viscosity-independent powers in $L\propto {\eta }^{2-3p}{\rho }^{p-1}{\sigma }^{2p-1}{t}^p$. They require thin interfaces, statistically self-similar bicontinuous domains, dominant advection and no additional relevant macroscopic scales. A degeneracy in the signed inertial coefficients can invalidate the pure-power force balance without invalidating the dimensional calculation.

Balancing viscous and capillary forces gives $L(t)\sim (\sigma /\eta )t$. This regime is independent of the mass density and applies below the viscous-to-inertial crossover scale.

Balancing inertia and capillarity gives $L(t)\sim (\sigma /\rho {)}^{1/3}{t}^{2/3}$. This regime is independent of viscosity and applies above the viscous-to-inertial crossover scale.

If stirring imposes a time scale $\tau$, replacing time derivatives by $1/\tau$ predicts a steady domain size. Viscous-capillary balance gives $\Lambda \sim \sigma \tau /\eta$, while inertial-capillary balance gives $\Lambda \sim (\sigma {\tau }^2/\rho {)}^{1/3}$.

A polar order parameter is a vector whose direction distinguishes head from tail. In an isotropic medium without an external polar field, spatial inversion makes the free energy even under $\mathbf{p}\mapsto -\mathbf{p}$.

The polar molecular field is the functional derivative of the free energy with respect to the polar order parameter. In the positive-derivative convention, dissipative dynamics is $D_t\mathbf{p}=-\Gamma \mathbf{h}$ with positive kinetic coefficient. For a quartic local term and an isotropic gradient penalty,Some authors instead define the molecular field with a minus sign, which changes the associated stress and evolution conventions.

The material derivative of a polar order parameter contains rigid rotation with the local vorticity and a material-dependent alignment response to the symmetric strain-rate tensor.

In simple shear flow, a persistent polar order parameter with $|\xi |<1$ has no fixed angle and rotates continuously. Integrating $\dot{\theta }=g(\xi \cos 2\theta -1)/2$ through $2\pi$ gives the displayed period. An unoriented nematic director repeats after half that time. At $\xi =0$, the magnitude relaxes independently while the angular speed is $-g/2$.

For a uniform polar order parameter in simple shear flow $\mathbf{v}=(gy,0)$, the angular equation is $\dot{\theta }=g(\xi \cos 2\theta -1)/2$. A fixed orientation for nonzero shear therefore requires $|\xi |\ge 1$. For $|\xi |>1$, its angular stability condition is $g\xi \sin 2\theta >0$; equality at $|\xi |=1$ requires nonlinear analysis. The ordered magnitude must separately satisfy the radial stationarity equation and be nonzero.

An advected order parameter contributes stress because a velocity gradient changes its free energy. For a polar vector this includes distortion, antisymmetric rotational, and symmetric flow-alignment stresses involving the molecular field.