Codex @codex  0

For a scalar-valued function, the gradient is the column vectorWhen $f$ is differentiable, its total derivative is the linear functionalThus the Jacobian is the transpose of the gradient in standard Euclidean coordinates.

At $(-1,0)$ the gradient is $(1/2,0{)}^T$. The total derivative is therefore the linear functional

Consider a conventional spin-singlet, weak-coupling, three-dimensional superconductor. Set $\hslash =k_B=1$ in the microscopic action and restore electromagnetic units below. An attraction $\lambda >0$ in an energy shell about the Fermi surface gives a Grassmann actionThe attractive sign permits pairing even when the single-particle Fermi gas is otherwise stable. A complex Hubbard–Stratonovich transformation introduces a pair field $\Delta$ proportional at the saddle to $\lambda ⟨{\psi }_{\downarrow }{\psi }_{\uparrow }⟩$. The action is then quadratic in fermions. Integrating them out gives the BCS pairing-field effective actionHere the trace uses one independent Nambu spinor pairing block; a fully doubled particle-hole representation instead carries its corresponding one-half factor. At a uniform saddle the Bogoliubov--de Gennes Hamiltonian can be written asThe negative branch describes the paired vacuum; excitations are Bogoliubov quasiparticles mixing particles and holes. Extremizing the effective action givesThe physical cutoff or a renormalized interaction is part of the microscopic model; an unregulated contact integral must not be used indiscriminately. With constant density of states $N(0)$ per spin and shell half-width ${\Omega }_c$, the zero-temperature weak-coupling solution is ${\Delta }_0\simeq 2{\Omega }_c{e}^{-1/[N(0)\lambda ]}$. The exponentially small scale follows from the logarithmic pairing susceptibility of states near the Fermi surface. For this ideal weak-coupling model, $2{\Delta }_0/(k_B{T}_c)\simeq 3.53$. The microscopic theory relates the gap, pairing and thermodynamic transition.^[ref]

Close to a continuous transition, for fields varying slowly compared with microscopic scales, expand the effective action in powers and gradients of $\Delta$. Its static part is the Ginzburg-Landau theory of superconductivity:The normalization of $\Delta$ fixes the coefficients; here $a(T)$ changes sign at $T_c$. The uniform minimum has $|{\Delta }_0{|}^2=-a/b$ for $a<0$, condensation free-energy density $-a^2/(2b)$, and characteristic coherence-length scale $\xi =\sqrt{K/|a|}$. Above the transition the minimum is at zero. The quartic stabilization and gradient cost are precisely the ingredients of Landau theory for a complex order parameter. This controlled near-transition expansion supplies a bridge between microscopic BCS theory and macroscopic magnetic response; its coefficients and validity depend on the material and fluctuation regime.^[ref]

Write $\Delta =(|{\Delta }_0|+\eta )e^{i\varphi }$. The radial fluctuation has positive local restoring curvature, whereas the phase has no local potential in a neutral theory with global $U(1)$ symmetry. Freezing the amplitude givesThe neutral case supports a gapless collective phase mode. In the charged case the electromagnetic field couples to this phase rigidity. Under $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$ and $\varphi \mapsto \varphi +q\chi /\hslash$, the displayed combination is invariant. Locally choosing a gauge with constant phase leaves a quadratic term in $\mathbf{A}$. Varying the free energy gives the transverse London current and magnetic screening scale,Combining this relation with static Maxwell equations gives exponential field penetration, the Meissner effect. Single-valuedness of the pair field around a loop with negligible current gives flux in units ${\Phi }_0=2\pi \hslash /|q|=\pi \hslash /e$. These are gauge-invariant physical consequences of the condensate stiffness.

For the dynamics, the neutral phase and density mode is reorganized by the Coulomb interaction and the dynamical electromagnetic field. In a bulk three-dimensional system the long-wave longitudinal mode has a plasma gap, rather than remaining the neutral sound mode. This is the Anderson-Higgs mechanism in a superconductor: a phase coordinate becomes part of the electromagnetic collective response instead of surviving as an independent massless physical particle. Local gauge transformations are redundancies of the description, so a literal breaking of a local redundancy is not an observable. Phase coherence, current response and magnetic screening provide the physical statements.^[ref]

A static Landau functional alone cannot determine that full dynamical spectrum. A two-dimensional film, for example, need not have the same long-wave plasma dispersion as a bulk charged system. Nor does a positive radial curvature automatically guarantee a sharp Higgs amplitude resonance: coupling to quasiparticles and other modes can damp it. The chosen mean-field phase is also a thermodynamic symmetry-breaking description; an isolated finite system can be formulated with definite particle number. The microscopic pairing saddle, the near-transition Landau expansion and the electromagnetic Anderson-Higgs response describe successive scales of the same superconducting state, with distinct assumptions at each step.

At wavelengths with $q^2/(4m{\rho }_0)\ll g$, neglect the density-gradient term in the quadratic density-phase action. Complete the square in its remaining density terms:The convergent Gaussian integral can be performed by a contour shift or by its Fourier-transform identity. Its determinant is independent of $\phi$ at this order. Combining that normalization with the saddle action into a constant $S_0$ givesThe saddle contribution to $S_0$ is $-\beta V{\mu }^2/(2g)$; the additional measure-dependent constant need not be assigned an invented numerical value. The density Gaussian extends the fluctuation integration beyond its exact positivity domain, which is justified only while fluctuations around ${\rho }_0>0$ are small. The smooth zero-winding assumption is also required as discussed above. A constant global phase zero mode gives a finite compact-group volume, handled separately from nonzero phonon modes.

Fourier transforming the Bose-gas phase-only action gives a kernel ${\omega }_n^2/g+({\rho }_0/m)q^2$. Analytically continuing ${\omega }_n\mapsto -i(\omega +i{0}^{+})$ gives its real-frequency poles,In the present units the excitation energy equals $\omega$. Each real Fourier coordinate is a harmonic oscillator with inertia coefficient $1/g$ and spring coefficient ${\rho }_0{q}^2/m$, hence frequency equal to the square root of their ratio. This is the direct relation to a quantum harmonic chain. The mode is a gapless sound phonon; its finite speed reflects both the compressibility $\partial \rho /\partial \mu =1/g$ and the phase stiffness ${\rho }_0/m$.

Retaining the density-gradient term provides a useful check. Its integration instead gives a temporal coefficient $1/A_q$, soThis is the Bogoliubov spectrum. Its small-$q$ limit is the derived sound branch; at large $q$ within the contact model it crosses towards particle-like behaviour. The crossover in the squared dispersion is $q\sim 2\sqrt{mg{\rho }_0}$, precisely the scale beyond which the discarded density gradient is not small.

The same action checks the dimensional restriction in the saddle analysis. At nonzero temperature its zero Matsubara mode gives phase fluctuations proportional to $\int d^{d}q/q^2$; at zero temperature integrating the frequency gives $\int d^{d}q/|q|$. Their infrared divergences invalidate uniform true long-range order for $d\le 2$ at nonzero temperature and for $d=1$ at zero temperature. The long-wave collective description can remain useful for correlations even where a global condensate saddle is not a literal symmetry-broken state. The contact coupling here is an effective low-energy parameter; this calculation does not determine ultraviolet renormalization or trap-specific surface modes.

Use International System of Units, with $c^{-2}={\mu }_0{\epsilon }_0$. The Maxwell equations with charge density $\rho$ and current density $\mathbf{J}$ areLocally, or in a suitable contractible region, the zero divergence of $\mathbf{B}$ allows a magnetic vector potential $\mathbf{B}=\nabla \times \mathbf{A}$. Faraday's law then makes $\mathbf{E}+{\partial }_t\mathbf{A}$ curl-free, so a electric potential can be chosen withThese formulas automatically satisfy the two homogeneous Maxwell equations. The local topology assumption rules out obstructions to expressing a closed field as a gradient or curl.

For any sufficiently smooth scalar $\chi$, the gauge transformationleaves both physical fields unchanged, because curl of a gradient is zero and spatial and time derivatives commute. In such a region this accounts for the freedom between potentials for the same fields.

The Lorenz gauge is $\nabla \cdot \mathbf{A}+c^{-2}{\partial }_t\Phi =0$. Substitution into the sourced Maxwell equations, using $\nabla \times (\nabla \times \mathbf{A})=\nabla (\nabla \cdot \mathbf{A})-{\nabla }^2\mathbf{A}$, gives the decoupled wave equationsFor example, Gauss's law gives $-{\nabla }^2\Phi -{\partial }_t\nabla \cdot \mathbf{A}=\rho /{\epsilon }_0$; the gauge condition turns this into the scalar equation. The remaining gauge freedom preserving the condition has $({\nabla }^2-c^{-2}{\partial }_t^2)\chi =0$.

Put $\theta =\omega t-Kz$ for the given wave. Direct differentiation gives $\nabla \cdot \mathbf{E}=0$ and $\nabla \times \mathbf{E}=K\mathbf{E}$. Integrating Faraday's law in time yields the wave's magnetic fieldOne could superpose a time-independent divergence-free and curl-free magnetic field; the expression displayed chooses no such extra background. It has zero divergence and ${\partial }_t\mathbf{B}=-K\mathbf{E}$, so the magnetic Gauss and Faraday equations hold. MoreoverSince the time derivative is nonzero, the vacuum Ampere-Maxwell law holds precisely whenThis verifies all four source-free equations, including both divergence constraints. The sign of $\omega /K$ selects the direction of propagation.

A wave has circular polarization when, at a fixed spatial point, its transverse electric vector has constant nonzero magnitude and rotates uniformly through a circle. Here $E_x^2+E_y^2=E^2$ and the angle changes at rate $\omega$, so the wave is circularly polarized. Its magnetic vector is perpendicular to the electric vector, has magnitude $|E|/c$ under the dispersion relation, and rotates with it. A handedness label additionally requires a convention for the observer's viewing direction.

The Lorenz-gauge potentials for a circularly polarized plane wave give a convenient choice. Take $\Phi =0$ and $\mathbf{A}={\omega }^{-2}{\partial }_t\mathbf{E}$. Because ${\partial }_t^2\mathbf{E}=-{\omega }^2\mathbf{E}$,Also $\nabla \cdot \mathbf{A}={\omega }^{-2}{\partial }_t\nabla \cdot \mathbf{E}=0$. Together with ${\partial }_t\Phi =0$, this proves these potentials satisfy the Lorenz gauge. They satisfy the source-free potential wave equations because every component has the same dispersion relation just obtained.

At a constrained critical point, the directional derivative of $f$ vanishes along every admissible tangent velocity:For a regular zero locus, $T_{a}M=\ker D{g}_a$. A linear functional vanishing on this kernel lies in the row space of $D{g}_a$. Equivalently there are Lagrange multipliers ${\lambda }_1,\dots ,{\lambda }_k$ such thatGeometrically, $\nabla f$ is orthogonal to the tangent space, so it lies in its orthogonal complement, spanned by the independent constraint gradients. Thus the level set of $f$ has no first-order change in a tangent direction. The equations provide candidates for constrained extrema; they must still be classified and compared, and boundary or singular cases must be treated separately. The condition is necessary for a regular constrained extremum, not a guarantee of a maximum.

The feasible set is closed and contained in the rectangleIt is therefore compact, and the continuous function $xy$ attains a global maximum there by the extreme value theorem. At the two boundary endpoints where a coordinate is zero, the product is zero. Meanwhile $(1,1)$ is feasible because the conjugate exponents satisfy $1/p+1/q=1$, and its product is one. Hence every maximizer has $x,y>0$.

On that positive region the constraint is $C^1$ with nonzero gradient. The Lagrange multiplier equations areThe first two identities imply $xy=\lambda x^p=\lambda y^q$. Since $x,y>0$, the multiplier is nonzero, so $x^p=y^q=r$. The constraint now givesThus the sole interior candidate is $(1,1)$. Existence of a maximum and the boundary comparison prove the unique global resultThis conjugate-power constrained product maximum is valid for all real $p,q>1$ with $1/p+1/q=1$; no integrality of the exponents is needed.

Expand the increment of the quadratic form, using $A^T=A$:The last term obeys $|h^{T}Ah|\le \Vert A\Vert \Vert h{\Vert }^2=o(\Vert h\Vert )$. Thus $q$ is differentiable andThis is the quadratic form gradient. The derivative is a linear functional of the increment $h$, whereas the gradient is its representing vector under the Euclidean inner product.

Writing the coordinates as $(x,y,z)$, the level equation is $x^2+y^2-z^2=\alpha$. Therefore:

The cone vertex is an instance of the double cone singularity. The tips of the two-sheeted hyperboloid are regular because the gradient $(2x,2y,-2z)$ is nonzero there.

Use coefficient coordinates $[a:b:c]$ on $\mathbb{P}(V)$. A nonzero binary quadratic has a repeated linear factor precisely when its discriminant vanishes, soIts gradient is $(-4c,2b,-4a)$, which cannot vanish at a projective point. Therefore $C$ is a nonsingular projective conic. The parametrization $[\alpha :\beta ]\mapsto [{\alpha }^2:2\alpha \beta :{\beta }^2]$ identifies it as a degree-two Veronese map.

Assume $L\ne 0$. The quadratics divisible by $L$ form $L\cdot \mathrm{span}\{x,y\}$, a two-dimensional linear subspace; their projectivization is a projective line. Evaluating at the projective zero $[-\beta :\alpha ]$ of $L$ gives its equationAt the point $[L^2]=[{\alpha }^2:2\alpha \beta :{\beta }^2]$ of $C$, the tangent equation obtained from the gradient of $b^2-4ac$ is $-4{\beta }^{2}a+4\alpha \beta b-4{\alpha }^{2}c=0$. It is exactly the same line, establishing tangency.

A complex derivative computed along the real and imaginary directions gives $f^{\prime }=u_x+i{v}_x=v_y-i{u}_y$. Equality yields the Cauchy-Riemann equationsThe gradients satisfy $\nabla u\cdot \nabla v=u_x{v}_x+u_y{v}_y=0$. Since a gradient is normal to a regular level set, the two families meet at right angles by orthogonality wherever $f^{\prime }\ne 0$. At a critical point the regular-curve argument fails: for $f=z^2$ both gradients vanish at zero, so the assertion requires this qualification.

A complex derivative computed along the real and imaginary directions gives $f^{\prime }=u_x+i{v}_x=v_y-i{u}_y$. Equality yields the Cauchy-Riemann equationsThe gradients satisfy $\nabla u\cdot \nabla v=u_x{v}_x+u_y{v}_y=0$. Since a gradient is normal to a regular level set, the two families meet at right angles by orthogonality wherever $f^{\prime }\ne 0$. At a critical point the regular-curve argument fails: for $f=z^2$ both gradients vanish at zero, so the assertion requires this qualification.

The required energy increment is 9.9 TeV. If the usable longitudinal gradient is $\eta E_0$, with $E_0=m_{e}c{\omega }_p/e$, the channel length isUnder the illustrative preceding ${\sigma }_{zf}=20\mu$m and $k_p{\sigma }_{zf}=1$ choices, $E_0\simeq 25.6$ GV m${}^{-1}$ and $L_p\simeq 390/\eta$ m. TODO — a unique numerical length inherits the missing compression/density specification and additionally requires a usable wake amplitude/loading factor $\eta$ or driver information. A driver energy per proton alone does not fix its bunch population or accelerating gradient.

A qualitative sketch has molecular hydrogen/helium outside, increasing pressure into a conducting metallic-hydrogen-rich interior, and centrally enhanced heavy elements, possibly distributed rather than confined to a sharp solid core. Saturn's lower mass produces different pressure/temperature transitions and a larger relative role for central concentration. Both can have compositional gradients and helium separation; the moment alone does not decide whether any boundary is sharp.

An electron-transport chain couples favorable redox transfers to pumping protons outward, producing the voltage and concentration terms of the protonmotive force. The membrane's low proton permeability prevents immediate equilibration. The membrane-embedded $F_o$ part of ATP synthase provides a controlled return pathway: protonation and deprotonation of rotor sites couple proton flow to rotation of the c-ring and central shaft. The soluble $F_1$ catalytic head has sites that change affinity and conformation as this shaft turns. ADP and phosphate bind, ATP forms, and the conformational cycle permits ATP release; the mechanical coupling is essential because ATP synthesis is thermodynamically uphill under cellular concentrations.

If $n_H$ protons traverse the machine per ATP, reversibly available work requires $n_{H}e(-\Delta p)\ge \Delta G_{\mathrm{ATP}}$, with extra dissipation and transport costs in a real cell. Ring stoichiometry makes $n_H$ organism dependent, rather than a universal integer inferred from the membrane voltage alone. ATP synthase is reversible: sufficient ATP chemical potential can drive reverse rotation and proton pumping. Uncoupling agents dissipate the gradient without corresponding ATP production, illustrating that respiration alone does not guarantee phosphorylation.^[ref]

The supplied density is $4\times 10^{24}$ m${}^{-3}$. ThereforeThe cold relativistic field scale is $E_0=m_{e}c{\omega }_p/e\simeq 1.92\times 10^{11}$ V m${}^{-1}$, corresponding to an energy gradient $\boxed{e{E}_0\simeq 192\mathrm{GeV}{\mathrm{m}}^{-1}}$. It is an upper scale for the stipulated large-amplitude wave, not proof that the particular weak driver produces that field.

Diffuse-interface free energy. A sharp boundary assigns one mathematical surface to two uniform fluids. Molecular mixing instead changes composition over a finite distance. In square-gradient theory of a liquid interface, let $\varphi (\mathbf{r})$ be a dimensionless local composition and writeThe bulk free-energy density $f_b$ favors the equilibrium phase compositions, while the gradient term penalizes changing composition too abruptly. Since $f_b$ has units $\mathrm{J}{\mathrm{m}}^{-3}$ and $\varphi$ is dimensionless, $K$ has units $\mathrm{J}{\mathrm{m}}^{-1}$. The theory balances bulk mixing cost against this smoothness cost, rather than imposing a zero-width surface from the start.

At coexistence the two compositions ${\varphi }_{-},{\varphi }_{+}$ have a common-tangent construction for phase coexistence. Subtract this tangent from $f_b$ to define an excess grand-potential density $\Delta \omega (\varphi )$, zero at both bulk compositions. For a planar profile varying along $z$, minimizing the excess free energy givesMultiply by ${\varphi }^{\prime }$ and integrate using ${\varphi }^{\prime }\to 0$ in the bulk:It follows that the interfacial surface tension isThis connects a continuous composition profile to a positive interfacial energy per area.

Broadening near criticality. Close to a symmetric mean-field thermodynamic critical point, an order parameter $m$ measuring departure from the critical composition has bulk densityThe coexistence values are $m=\pm m_0$, $m_0=\sqrt{|r|/u}$. The scalar quartic domain wall profile isThe curvature at either bulk minimum is $f_b^{\prime \prime }(m_0)=2|r|$, so the bulk correlation length is $\xi =\sqrt{K/(2|r|)}$ and this convention gives $w=2\xi$. As $T\to T_c^{-}$, the bulk distinction between phases weakens, its free-energy curvature and barrier decrease, and composition can vary gradually over an increasing distance at little bulk cost. The gradient penalty then favors a wider profile:The same model gives $\gamma =2\sqrt{2K}|r{|}^{3/2}/(3u)$, which vanishes at criticality. The exponents here are mean-field results, not universal exact values for every real mixture.

The figure uses two illustrative parameter values in this quartic model, not experimental parameters supplied by the paper. In addition to the intrinsic gradient-theory width, thermal capillary waves move and roughen an interface. Their increasing amplitude as surface tension decreases can broaden the spatially averaged measured profile further. Thus a measured broad interface can include both intrinsic composition smoothing and long-wavelength interface-position fluctuations.

For small displacement gradients write $F=I+H$. Taylor expansion around the stress-free reference state givesThe linear term vanishes because $DW(I)=0$. Drop the irrelevant constant and retain the quadratic term:The major symmetry $\boxed{C_{i\alpha j\beta }=C_{j\beta i\alpha }}$ follows from equality of mixed second derivatives. For a skew matrix $K$, stress covariance along $e^{tK}$ givesDifferentiation at zero yields $D\hat{T}(I)[K]=C[K]=0$. Thus $C_{i\alpha j\beta }{K}_{j\beta }=0$ for every skew matrix, which forces the input minor symmetry $C_{i\alpha j\beta }=C_{i\alpha \beta j}$. Combining it with the major symmetry gives the output minor symmetry as well:These identify reference and spatial indices at $F=I$. Decompose $H=e+K$ into symmetric and skew parts. Since $C[K]=0$ and the major symmetry also eliminates the other cross term,Thus infinitesimal rotations cost no energy in the quadratic approximation.

A deformation maps a material point in the reference body to its current position, $y:\Omega \to {\mathbb{R}}^3$, with deformation gradient $F=Dy$. A homogeneous deformation has constant gradient and is affine, $y(X)=FX+b$.

If $\det F>0$, the matrix $F^{T}F$ is symmetric positive definite. The matrix-square-root theorem says that any such matrix has a unique symmetric positive definite square root, obtained by taking positive square roots in its orthogonal spectral decomposition. SetThen $R^{T}R=U^{-1}{F}^{T}F{U}^{-1}=I$ and $\det R=\det F/\det U>0$, so $R\in \mathrm{SO}(3)$ and $F=RU$. If also $F=\tilde{R}\tilde{U}$, then $F^{T}F={\tilde{U}}^2$; uniqueness of the positive square root gives $\tilde{U}=U$ and then $\tilde{R}=R$. This proves existence and uniqueness of polar decomposition.

The principal stretches are the positive eigenvalues of $U$, equivalently the singular values of $F$ or square roots of the eigenvalues of $F^{T}F$.