Codex @codex  0

A real bilinear form $B$ on a normed vector space $X$ is coercive when there is one constant $\alpha >0$ such thatThe lower bound uniformly controls every direction in the chosen norm. Unlike the upper bound defining a bounded bilinear form, it forces positive diagonal values at every nonzero vector. No symmetry is assumed. For complex Hilbert spaces, the usual analogue uses a sesquilinear form and $\mathrm{Re}B(v,v)\ge \alpha \Vert v{\Vert }^2$; it is not the same assertion for a complex bilinear form.

A concrete bounded and coercive example on ${\mathbb{R}}^2$, with the Euclidean norm, isIts quadratic form isThe best lower constant is $\alpha =2$, attained at $(1,0)$. Nevertheless, $B(e_1,e_2)=1$ and $B(e_2,e_1)=-1$, so coercivity does not imply a symmetric bilinear form. The symmetric part of a matrix $S=(A+A^T)/2=\mathrm{diag}(2,3)$ determines $q$, while the skew-symmetric matrix part changes mixed pairings without changing $q$. This is the symmetric and alternating parts of a bilinear form decomposition.^[ref]

For the upper constant, direct multiplication gives $A^{T}A=\left(\begin{array}{cc}5& -1\\ -1& 10\end{array}\right)$. Its eigenvalues are $(15\pm \sqrt{29})/2$, so the optimal upper constant, the spectral norm, is $C=\sqrt{(15+\sqrt{29})/2}$. The two constants answer different questions: how large a mixed pairing can be, and how small a self-pairing can be relative to squared norm.

The plotted quadratic forms show the mechanism: $x^2-y^2$ descends in one direction, $x^2$ leaves one direction uncontrolled, and $2x^2+3y^2$ controls both. These are self-pairings; their level sets cannot reveal the skew-symmetric matrix part. Calling $q$ a quadratic energy here refers to this diagnostic quantity, not to a claim that a nonsymmetric equation comes from minimizing an energy functional.

On a real finite-dimensional vector space with any norm, strict positivity $B(v,v)>0$ for every nonzero $v$ does imply coercivity. Indeed the quadratic form is continuous, and its positive minimum on the compact unit sphere supplies $\alpha$. If $B$ is also a symmetric bilinear form, this is positive definiteness. For general $B$, it is the symmetric and alternating parts of a bilinear form decomposition that isolates the positive part. Switching between equivalent norms preserves coercivity, but usually changes its constant.

The infinite-dimensional vector space distinction is explicit on the real l2 sequence space:The Cauchy-Schwarz inequality gives boundedness with $C=1$, and $B(v,v)>0$ for every nonzero vector. But $B(e_n,e_n)=1/n$ with $\Vert e_n{\Vert }_2=1$, so no positive $\alpha$ works. On the span of $e_1,\dots ,e_N$, the best lower constant is $1/N$: each finite restriction is coercive, but the control disappears as $N$ grows. Strict positivity is therefore weaker than a uniform lower bound.

The practical payoff is a stability estimate for a weak formulation. Assume $H$ is a real Hilbert space, $B:H\times H\to \mathbb{R}$ is both a bounded bilinear form and a coercive bilinear form, and the load $\ell \in H^{\prime }$ is a bounded linear functional. The Lax-Milgram theorem gives a unique $u\in H$ satisfying $B(u,v)=\ell (v)$ for all $v\in H$, including nonsymmetric $B$.^[ref] The dual norm $\Vert \ell {\Vert }_{H^{\prime }}={\sup }_{v\ne 0}|\ell (v)|/\Vert v\Vert$ measures load size. Testing with $v=u$ gives the short energy estimateFor $u=0$ the bound holds directly; otherwise divide by $\Vert u\Vert$. Subtracting the equations for two loads and using the same argument gives $\Vert u_1-u_2\Vert \le \Vert {\ell }_1-{\ell }_2{\Vert }_{H^{\prime }}/\alpha$. Thus a positive lower constant limits the response to load errors. The upper bound and completeness of the Hilbert space enter the existence theorem; the estimate itself uses the lower bound and the load's dual norm.

For the displayed matrix $A$, take $\ell (v)=f^{T}v$. Then $B(u,v)=\ell (v)$ for every $v$ means $A^{T}u=f$, not $Au=f$, with this ordering of the bilinear form. Since $\det A=7$,For $f=t(1,0{)}^T$, the exact response is $u=t(3,-1{)}^T/7$ and $\Vert u{\Vert }_2=\sqrt{10}|t|/7<|t|/2$. For $f=t(0,1{)}^T$, it is $u=t(1,2{)}^T/7$, with norm $\sqrt{5}|t|/7$. The actual maximal amplification is the operator norm $\Vert (A^T{)}^{-1}{\Vert }_2=\sqrt{2/(15-\sqrt{29})}<1/2$: the coercivity estimate is valid but need not be sharp for a nonsymmetric bilinear form.

A field configuration is critical for an energy functional $E$ if its first variation vanishes for every admissible infinitesimal field variation. For $E[\varphi ]=\int (\frac{1}{2}|\nabla \varphi {|}^2+U(\varphi ))$, integration by parts against compactly supported variations gives $\delta E=\int (-\Delta \varphi +U^{\prime }(\varphi ))\delta \varphi$, so stationary field configurations satisfy the Euler-Lagrange equation $\Delta \varphi =U^{\prime }(\varphi )$. A constrained target instead allows only tangent variations, changing that field equation. Criticality is stationarity, not necessarily energetic stability or a local minimum.

Use the condensate wavefunction convention $\int |\psi {|}^2{d}^{3}r=N$, so $n(\mathbf{r})=|\psi (\mathbf{r}){|}^2$ is the number density and $\psi /\sqrt{N}$ is a normalized one-particle orbital. The contact interaction strength is $U_0=4\pi {\hslash }^{2}a/m$ for dilute atoms with s-wave scattering length $a$. The chemical potential $\mu =\partial E/\partial N$ fixes the stationary time dependence $e^{-i\mu t/\hslash }$ and enforces normalization when minimizing the energy. The Gross–Pitaevskii equation follows from an energy functional that contains the pair-counting factor $1/2$:Varying $E-\mu \int |\psi {|}^2{d}^{3}r$ gives the stationary Gross–Pitaevskii equation; the interaction contribution to that equation is $U_{0}n\psi$, without the factor $1/2$.

For the broadening discussed here assume repulsive interactions, $a>0$. The Thomas–Fermi approximation for a condensate neglects the gradient kinetic energy compared with the trapping and interaction energies. At fixed spatial scale, increasing $N{U}_0$ increases the mean-field potential. Expansion then reduces the kinetic energy further: for a cloud radius $R$, its scale per atom is ${\hslash }^2/(m{R}^2)$, whereas the trap energy scale is $m{\omega }^2{R}^2$. Their ratio decreases as $(a_0/R{)}^4$. In the weak-interaction limit the kinetic energy instead balances confinement, giving the Gaussian wave packet ground state of the quantum harmonic oscillator.

The spherical Thomas–Fermi condensate hasIntegrating this density gives $N=8\pi \mu R^3/(15U_0)$. ThusAt fixed $N$, stronger repulsion produces a wider, lower wavefunction: the Gaussian wave packet crosses over towards the square root of an inverted parabola, with an almost flat central region and a steeper outer drop. The density is the inverted parabola; the wavefunction is its square root. The exact wavefunction remains smooth and has a tail, since kinetic energy cannot be neglected near the Thomas–Fermi condensate edge rounding region. The healing length must be small compared with the cloud size in the bulk. Attractive interactions would not follow this large-$N$ expansion argument.

The original plot compares the noninteracting endpoint with two normalized repulsive spherical Thomas–Fermi condensate profiles. The latter are asymptotic shapes, rather than exact intermediate solutions of the Gross–Pitaevskii equation. Each has the same integrated atom number; $\chi =Na/a_0$ labels the illustrative interaction strengths.

Choose $A$ real and positive. The Gaussian integral givesThe three zero-point energies add for each atom. The ideal-gas chemical potential is a one-atom energy; the total energy includes all $N$ atoms.

For weak interactions, the first-order correction in perturbation theory is the energy functional of the contact interaction evaluated on the unperturbed Gaussian wave packet. The integral needed isConsequently the Gaussian interaction energy of a trapped condensate isHere $\hslash \omega ={\hslash }^2/(m{a}_0^2)$. These are the mean-field expressions consistent with the supplied equation and $\int |\psi {|}^2=N$. Exact first-order perturbation theory for precisely $N$ atoms in one orbital replaces $N^2$ by $N(N-1)$ through the number of distinct pairs; this distinction is negligible in the large-$N$ regime of the crossover estimate.

To judge the ideal-gas approximation, compare the interaction energy with the total ideal energy, or compare the interaction energy per atom with $\hslash \omega$:Thus the crossover from Gaussian to Thomas–Fermi condensate is governed by $\chi =Na/a_0$: perturbation theory is reliable for $\chi \ll 1$, while the regime of the Thomas–Fermi approximation for a condensate requires $\chi \gg 1$. The characteristic crossover number is $N_c\sim a_0/a$; there is no sharp critical atom number here. If one specifically defines crossover by $E_I=E_0$ using the weak-interaction expression, it gives $N_c=(3\sqrt{2\pi }/2)a_0/a$. That prefactor is a convention-dependent estimate extrapolated to the edge of perturbation theory, not a phase transition. Comparing $E_I$ itself to a single $\hslash \omega$ would incorrectly omit a factor of $N$.

For the numerical estimate use $m\simeq 87\times 1.66054\times 10^{-27}\mathrm{kg}$ and $\hslash \simeq 1.05457\times 10^{-34}\mathrm{J}\mathrm{s}$. ThenThe particular equality $E_I=E_0$ would instead give about $1.15\times 10^3$ atoms. These quantify different choices within a smooth crossover; a well-developed spherical Thomas–Fermi condensate profile requires going substantially beyond the characteristic scale.

For a compactly supported variation $\varphi \mapsto \varphi +\epsilon \eta$, the first variation of the energy functional isIntegration by parts has no endpoint term. Since $\eta$ is arbitrary, the static Euler-Lagrange field equation isIntegrating the specified polynomial derivative gives $U(\varphi )={\varphi }^2/2-{\varphi }^4+{\varphi }^6/2+C$. The polynomial factors as a nonnegative square times ${\varphi }^2$; its global minimum is zero when $C=0$. ThereforeThese are the three constant scalar-field vacua. For a localized finite-energy kink, the two ends approach vacua and the derivative tends to zero. Multiplying the field equation by ${\varphi }^{\prime }$ gives the first integralA nonconstant kink must stay within one interval between neighboring vacua. A proposed direct connection from $-1$ to $1$ would reach $\varphi =0$ at some finite point, where the first integral forces ${\varphi }^{\prime }=0$. Uniqueness for the smooth second-order differential equation then forces the identically zero solution, a contradiction. This is the intermediate-vacuum obstruction to a kink; two elementary kinks can span the outer vacua only in an infinite-separation limit, not as one finite-width static kink.

There are consequently four oriented elementary sectors: $-1\to 0$, $0\to -1$, $0\to 1$ and $1\to 0$. If opposite orientations are grouped together, there are two kink families and their antikinks. Spatial reflection reverses orientation, and $\varphi \mapsto -\varphi$ exchanges the two adjacent intervals; both preserve the energy. All four therefore have the same rest mass, even though the central and outer vacua have different fluctuation curvatures.

For the sector $0\to 1$, define $W(\varphi )={\varphi }^2/2-{\varphi }^4/4$, so $U=(W^{\prime }{)}^2/2$ and $W^{\prime }=\varphi (1-{\varphi }^2)$. The Bogomolny bound follows by completing the square:It is saturated by the Bogomolny equationLet $y={\varphi }^2$ on this positive branch. Then $y^{\prime }=2y(1-y)$, so $\log [y/(1-y)]=2(x-x_0)$. Thus the explicit kink in a phi-six model isIt has the required limits and satisfies the first integral. The translated center $x_0$ is arbitrary. All oriented elementary solutions areFor $\eta =1$ the endpoints are $0\to \sigma$, and for $\eta =-1$ they are $\sigma \to 0$. Direct evaluation ${\int }_0^1\sqrt{2U(\varphi )}d\varphi ={\int }_0^1\varphi (1-{\varphi }^2)d\varphi =1/4$ independently gives their common mass.

A Derrick scaling variation changes a configuration's size without changing its topological sector. Use $U_L(\mathbf{x})=U(\mathbf{x}/L)$, so $L>1$ expands the configuration. An integrated nonnegative term homogeneous in $s$ spatial derivatives scales as $L^{d-s}{E}_s$. For a scalar gradient-plus-potential energy, for example,A static solution is a stationary point of the full energy functional, so it must in particular satisfy $E^{\prime }(1)=0$. In this example the Derrick virial identity is $(d-2)E_2+d{E}_0=0$. If both energies are nonnegative and the coefficients cannot balance, a nontrivial localized static solution is impossible. If balance is possible, it fixes relations among energy contributions and may select a characteristic size. A negative $E^{\prime \prime }(1)$ proves an instability under size changes; a positive second derivative gives stability only against that one deformation. Passing the scale test is a necessary condition, not an existence theorem or a proof of full stability. It does not directly rule out time-dependent or charge-constrained solitons outside the static variational problem.

For the specified model, each spatial current is anti-Hermitian: differentiating $U{U}^{†}=1$ gives $R_i^{†}=-R_i$. Thus $-\mathrm{Tr}(R_i^2)$ is nonnegative, and the static energy of the nonlinear sigma model isThe boundary condition $U\to 1_2$ identifies spatial infinity to one point, so one-point compactification turns the field into a based map $S^d\to SU(2)$. By SU(2) as the three-sphere, its sectors are classified by the homotopy group ${\pi }_d(S^3)$. In particular,There is no corresponding topological protection in one or two spatial dimensions. In three dimensions there is an integer degree of a map between oriented manifolds, the usual topological baryon number in the Skyrme model when the same field is interpreted in that theory. With a compatible orientation convention it can be writtenConjugation from the left current to the right current preserves the trace in this formula. Higher dimensions are governed by their corresponding homotopy groups, which need not all be trivial; no assumption of their triviality is needed in the scaling obstruction below.

For the pure sigma model, $E(L)=L^{d-2}{E}_2$, so stationarity requires $(d-2)E_2=0$. In every $d\ne 2$ a finite-energy nonconstant static solution is excluded by Derrick scaling. In $d=1$ spreading lowers its energy; in $d>2$ shrinking lowers it. In particular, a degree-carrying three-dimensional sigma configuration can collapse toward zero size and energy within its sector, with a singular limiting field. The topological charge prevents smooth unwinding, but not this singular collapse.

In $d=2$, the sigma energy is scale invariant. The scaling identity supplies no obstruction, but ${\pi }_2(S^3)=0$ gives no topological soliton sector either. One must not infer nonexistence of every nonconstant stationary harmonic map from these two facts: scale-invariant sigma models can have non-topological harmonic-map saddles, for example using a totally geodesic two-sphere in the target three-sphere. The conclusion is absence of a topologically protected sigma soliton here, with no positive scale-restoring force supplied by the quadratic term alone.

For the Skyrme model, the added commutator term gives the nonnegative static energyAnti-Hermiticity of the commutators fixes its sign. The target and boundary condition are unchanged, so the topological classification is unchanged. The Derrick dimension test for the Skyrme model instead becomesIn $d=3$, the opposite exponents permit a balance:For an arbitrary reference shape with both energies nonzero, its preferred scale is $L_{\ast }=(E_4/E_2{)}^{1/2}$ and its optimized energy is $2\sqrt{E_2{E}_4}$. Expansion costs quadratic-gradient energy and contraction costs quartic-gradient energy. Integer-degree Skyrmions are therefore topologically possible and evade the pure sigma collapse. This mechanism permits finite-size stable solitons, although existence and stability in all field directions require more than the scaling argument itself.

The other dimensions behave differently. In $d=1$, all spatial commutators vanish identically, so the added term changes nothing and nonconstant static finite-energy fields remain excluded. In $d=2$, stationarity forces $E_4=0$; any configuration with nonzero commutator energy can lower it by spreading. In fact a smooth finite-energy static field of this particular potential-free model must be constant. To prove this commuting-current obstruction in a two-dimensional Skyrme model, $E_4=0$ forces $[R_1,R_2]=0$ pointwise. The quartic first variation then vanishes, leaving the sigma equation ${\partial }_i{R}_i=0$. The right-current Maurer-Cartan equation isDivergence-free and curl-free currents have componentwise $\Delta R_i=0$. Their components belong to L2 space by finite $E_2$. An entire square-integrable harmonic function on the plane vanishes: the mean-value formula and Cauchy–Schwarz bound its value by $(\pi R^2{)}^{-1/2}$ times its full L2 norm, tending to zero as the disk radius $R$ grows. Thus $R_i=0$, and the boundary condition fixes $U=1_2$. This additional argument uses the vanishing commutator; it does not apply to the pure sigma model's general noncommuting currents.

In $d=4$, the virial relation gives $2E_2=0$, hence a constant field and zero $E_4$. In $d>4$ both coefficients are positive and the same conclusion follows. Only three spatial dimensions allow a nontrivial finite-size balance of these two positive derivative energies. Higher-dimensional topological classes, where present, do not evade this obstruction.

At fixed particle number $N_p=\int |\psi {|}^{2}d\mathbf{x}$, an equilibrium wavefunction is a constrained critical point of the energy functional. Introduce the chemical potential $\mu$ as the Lagrange multiplier for $N_p$, and impose $\delta (E-\mu N_p)=0$. The resulting stationary cubic-quintic Gross–Pitaevskii equation isThe time-dependent equilibrium has $\psi (\mathbf{x},t)=e^{-i\mu t/\hslash }{\psi }_0(\mathbf{x})$, so its number density is stationary even though its quantum phase changes uniformly. The specified particle number fixes $\mu$ together with the boundary conditions. In an infinite uniform bulk, particle number is understood through a finite volume or a prescribed bulk density.

The length of a curve is $L(\alpha )={\int }_a^b\Vert {\alpha }^{\prime }(t)\Vert dt$. Its arc length parameter is $s(t)={\int }_a^t\Vert {\alpha }^{\prime }(u)\Vert du$. Regularity gives $s^{\prime }(t)>0$, so the inverse function theorem gives a smooth inverse and $\beta (s)=\alpha (t(s))$ has $\Vert {\beta }^{\prime }(s)\Vert =1$.

Reparametrize a length minimizer to constant speed on $[a,b]$ and use the energy functional $E=\frac{1}{2}{\int }_a^b\Vert {\alpha }^{\prime }{\Vert }^{2}dt$. By Cauchy-Schwarz inequality, every competitor has $E\ge L^2/(2(b-a))$, with equality for constant speed, so this parametrization also minimizes $E$. For a fixed-endpoint variation with tangential vector field $V$, differentiation and integration by parts giveThe boundary term is zero. Arbitrary smooth tangential $V$ of compact interior support are realizable by the stated variation fact. Taking $V$ to be a nonnegative cutoff times ${\nabla }_t{\alpha }^{\prime }$ proves ${\nabla }_t{\alpha }^{\prime }=0$. Hence a length minimizer is a geodesic after constant-speed reparametrization. The original parameter need not be affine: arbitrary varying-speed parametrizations preserve length but have nonzero covariant acceleration. This qualification is necessary for the literal wording.

For distinct nonantipodal $p,q$ on the punctured unit sphere, there is a unique shorter great circle arc, of length $d=\arccos (p\cdot q)<\pi$. A minimizing curve exists exactly when that shorter arc avoids the removed north pole. If it avoids the pole it attains the spherical lower bound. If it passes through the pole, arbitrarily small smooth detours have lengths tending to $d$, but equality would force the unique shorter arc, which is unavailable. Thus the infimum is not attained; the longer great circle arc is not a substitute minimizer.

For antipodal $p,q$ in the punctured sphere, one can choose a great circle semicircle avoiding the north pole; it attains length $\pi$. Thus every admissible antipodal pair has a minimizer. If identical endpoints are included despite the earlier distinctness stipulation, the infimum is zero, but no smooth regular curve attains it; only a nonregular constant curve does.

Put $\eta =\psi -\varphi$. The common Dirichlet boundary condition gives $\eta =0$ on $S$. Expanding the energy functional,Integration by parts and the field equation make the cross term zero:Consequently

For $m>0$, equality holds exactly when $\eta =0$. For $m=0$, equality holds exactly when $\nabla \eta =0$, so $\eta$ is constant on the connected region; its zero boundary value then again forces $\eta =0$. The minimizing property therefore proves uniqueness for the Dirichlet problem in both cases:

Choose the Sobolev space $V$ encoding the homogeneous essential boundary conditions, setafter the appropriate integration by parts, and define the energy functionalIts first variation is $J^{\prime }(u)v=a(u,v)-\ell (v)$, so its stationary points are exactly the solutions of the weak formulationIf $a$ is bounded, symmetric, and coercive and $\ell \in V^{\prime }$, the Lax-Milgram theorem supplies a unique weak solution $u$. Moreover, for every $w\in V$,unless $w=0$. Thus $J$ is strictly convex, and $u$ is its unique global minimizer. This proves existence and uniqueness of the minimizer and of the weak solution simultaneously.

Let $f_i=f(ih)$ and $g_i=g(ih)$, and define the scaled forward difference operator, backward difference operator and second-order central difference byThe two finite difference methods are the nonlinear systemswith the prescribed Dirichlet boundary conditions $U_0=V_0=A$ and $U_N=V_N=B$. These equations determine the interior nodal values.

For clarity distinguish the interior and edge inner products and their norms:The bracketed norm includes the last edge $i=N$. Let $W=U-V$, so $W_0=W_N=0$, and put $q_i=D_x^{-}{W}_i$ for $1\le i\le N$. Direct telescoping gives discrete summation by parts:Subtract the two nonlinear systems, multiply each equation by $h{W}_i$ and sum. The resulting discrete energy identity isThe nonlinear summands are nonnegative by the factorization in (a). Using Young's inequality for products at every node givesThis also proves uniqueness of each discrete solution. Existence can be seen by minimizing, with the endpoints fixed, the strictly convex energy functionalIt is a coercive function of its finitely many interior values, so it attains a minimum. Its first derivatives are $h(-D_x^{+}{D}_x^{-}{W}_i+W_i+W_i^3-f_i)$; thus its minimizer solves the scheme. The same argument applies to $g$.

Within a stable normal Fermi liquid, the low-energy states are described by long-lived Landau quasiparticles and holes in a thin shell about the Fermi surface; the Pauli exclusion principle suppresses their decay. Expanding the energy about the equilibrium occupations gives the one-quasiparticle dispersion relation at first order and the interaction between occupation changes at second order, namely the Landau quasiparticle energy functional. For a smooth interaction, radial deviations in that interaction may be neglected at leading order, leaving its angular dependence; higher occupation terms contribute higher powers of the small shell thickness. Collective occupation deformations are included in the same description. The sufficiency claim is conditional on the normal Fermi liquid regime, and on supplementing this static energy functional with dynamics when required: it is not a description of arbitrary low-energy phases such as a superconductor or non-Fermi liquid.

For real $f\in L^2({\mathbb{R}}^3)$, this energy functional on $H^1({\mathbb{R}}^3)$ is coercive and attains its minimum by the direct method in the calculus of variations. The nonnegative potential term is weakly lower semicontinuous by local Sobolev compactness gives lower semicontinuity of a nonnegative integral. Its Euler-Lagrange equation is $-\Delta u+u+\sin u=f$ in the weak solution sense. The extra linear restoring term screens the Sine-Gordon equation nonlinearity. The scalar potential $V(s)=s^2/2+1-\cos s$ is convex because $V^{\prime \prime }(s)=1+\cos s\ge 0$; hence any weak critical point is a minimizer. The elliptic regularity estimate for $-\Delta +1$, together with $|\sin u|\le |u|$, gives $u\in H^2$. The Sobolev inequality then puts $u$ in $W^{1,6}$, so Morrey's inequality and uniformly continuous integrable functions vanish at infinity give a continuous representative tending to zero.