Codex @codex  0

For the catenoid profile $y(x)=E\cosh (x/E)$ spanning fixed-endpoint rings, the second variation of $F[y]=2\pi \int y\sqrt{1+y^{\prime 2}}dx$ isThis refers to axisymmetric graph variations, not all possible changes in surface topology. Strict positivity for every nonzero admissible variation gives a strict local minimum in the $C^1$ topology: the positive lowest Dirichlet eigenvalue of the regular Sturm-Liouville problem gives coercivity, and sufficiently small coefficient perturbations preserve it. The zero variation must be excluded from a strict-positivity hypothesis.

For a regular self-adjoint Sturm-Liouville problem with positive weight $w$ and self-adjoint boundary conditions, eigenfunctions with distinct eigenvalues are orthogonal in the weighted inner product $\int w\overline{u}v$. Integration by parts proves this by equating $\int \overline{u}Lv$ and $\int \overline{Lu}v$.

Let primes denote radial derivatives of the equilibrium and put $\Delta u=u^{\prime }+2u/r$. The linearized continuity, adiabatic pressure and Poisson equation areTake the time derivative of the Poisson equation and substitute continuity. Integration in radius givesCentre regularity excludes a changing point mass at the origin, so $C(t)=0$ andThe perturbed self-gravity must be retained; this calculation does not make the Cowling approximation.

The radial linear momentum equation is $\rho {\partial }_{t}u=-{\partial }_r\delta p-\delta \rho {\Phi }^{\prime }-\rho {\partial }_r\delta \Phi$. Differentiate it in time and use the preceding equations:The equilibrium has $p^{\prime }=-\rho {\Phi }^{\prime }$ and ${\Phi }^{\prime \prime }+2{\Phi }^{\prime }/r=4\pi G\rho$. The remaining terms simplify asConsequently the radial stellar oscillation equation isThe derivative acts on the full variable coefficient $\gamma (r)p(r)$; discarding ${\gamma }^{\prime }$ would change the result.

For $u=r\xi (r)e^{-i\omega t}$, $\Delta u=(3\xi +r{\xi }^{\prime })e^{-i\omega t}$. Substitute and multiply by $-r^3$. Combining the two terms proportional to ${\xi }^{\prime }$ into a total derivative gives the radial stellar pulsation equation,

Set $P=\gamma p{r}^4$, $W=\rho r^4$ and $Q=r^3[(4-3\gamma )p{]}^{\prime }$. The Sturm-Liouville operator isFor regular physical perturbations, $\xi$ is finite at the centre and $P{\xi }^{\prime }\to 0$ at the free surface. For a nonzero-frequency mode the fluid displacement is $u/(-i\omega )$, so vanishing Lagrangian pressure perturbation is equivalent to $\gamma p(3\xi +r{\xi }^{\prime })\to 0$. With $p(R_s)=0$ and finite $\xi$, this gives $P{\xi }^{\prime }\to 0$. Zero-frequency modes use the same displacement boundary condition directly. Assume bounded $\gamma$, positive $\rho ,p,\gamma$ in the interior and the usual finite-energy endpoint domain. Integration by parts givesThus this physical self-adjoint differential operator has real eigenvalues ${\omega }^2$. The endpoint domain is essential: the fact that $p$ vanishes does not by itself allow arbitrary singular trial functions. The regular free-surface realization of the Sturm-Liouville problem is the one used here.

The Rayleigh-Ritz variational principle gives the fundamental squared frequency as the infimum of the weighted stellar pulsation Rayleigh quotient,The infimum is over admissible finite-energy functions satisfying the physical endpoint conditions. If this quotient is nonnegative for every such function, all radial frequencies are real and there is no exponentially growing radial mode. A negative value for even one trial function proves a negative eigenvalue and an exponentially growing solution, because ${\omega }^2<0$. A zero lowest value is marginal and needs separate treatment of neutral motion.

Choose the homologous trial function $\xi =1$, which corresponds to radial velocity proportional to $r$. Its gradient contribution is zero, andThe denominator $\int \rho r^{4}dr$ is positive. HenceThis pressure-weighted radial instability criterion is sufficient, not necessary: another trial function can detect instability even if this one does not. For a constant stellar adiabatic exponent it recovers instability below $4/3$. At constant $\gamma =4/3$, the homologous mode is neutral. For constant $\gamma >4/3$, $Q=(4-3\gamma )r^3{p}^{\prime }>0$ in a normally stratified hydrostatic equilibrium, so the quotient is positive and the star is radially stable.

Take a static, nonrotating, nonmagnetic spherical equilibrium with $\nabla P=-\rho \nabla \Phi$ and ${\nabla }^2\Phi =4\pi G\rho$. Use a fluid displacement $\mathbf{\xi }(\mathbf{r})e^{-i\omega t}$. Write $p_1,{\rho }_1,{\varphi }_1$ for Eulerian fluid perturbations. The Lagrangian pressure perturbation and corresponding mass density change obey $\Delta P=p_1+\mathbf{\xi }\cdot \nabla P$, $\Delta \rho ={\rho }_1+\mathbf{\xi }\cdot \nabla \rho$. Linearizing mass conservation, the Euler equations for an inviscid fluid and Poisson equation for Newtonian gravity givesThere is no equilibrium acceleration to multiply a perturbed mass density. The adiabatic process condition, with composition carried by the parcel, closes the system:These are the complete linear adiabatic stellar oscillation equations, including the perturbation of self-gravity. Neglect of heat exchange is appropriate when oscillation periods are short compared with relevant thermal relaxation times; it does not determine nonadiabatic excitation or damping.

For the pressure and buoyancy modes, separate angular dependence using spherical harmonics:The horizontal amplitude ${\xi }_h$ has dimensions of length. Angular differentiation gives $\nabla \cdot \mathbf{\xi }=[r^{-2}(r^2{\xi }_r{)}^{\prime }-\ell (\ell +1){\xi }_h/r]Y_{\ell }^m$, and tangential momentum gives ${\omega }^{2}r{\xi }_h=p_{\ell }/\rho +{\varphi }_{\ell }$. Define the adiabatic sound speed, stellar buoyancy frequency and Lamb frequency byThe adiabatic mass density relation becomes ${\rho }_{\ell }/\rho =p_{\ell }/(\rho c_s^2)+(N^2/g){\xi }_r$. Substitution produces a radial form of the full oscillation equations for nonzero $\omega$:The apparent $N^2/g$ factor is evaluated through its defining gradient at the centre rather than by dividing two zeros.

Regularity at the centre excludes singular solutions. At a free surface the Lagrangian pressure perturbation vanishes, $\Delta P=0$; outside the star the gravitational perturbation decays as $r^{-\ell -1}$. For a model with mass density tending to zero at its surface, continuity of ${\varphi }_{\ell }$ and its radial derivative gives ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R$. If the equilibrium mass density jumps to vacuum, include the displaced-surface mass sheet: the outward-minus-inward derivative jump is $4\pi G\rho (R){\xi }_r(R)$, so the interior condition is ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R-4\pi G\rho (R){\xi }_r(R)$. An atmospheric boundary condition can replace the ideal free surface.

These conditions make ${\omega }^2$ an eigenvalue, not an arbitrary local sound frequency. With conservative boundary conditions the adiabatic operator is self-adjoint, giving real ${\omega }^2$; negative values describe instability. The frequencies depend on $P(r),\rho (r),{\Gamma }_1(r)$, self-gravity, stratification and boundaries, as well as angular degree $\ell$ and radial order. In a spherical nonrotating star they are degenerate in $m$. For homologous equilibrium structures, their scale isthus the typical oscillation time measures inverse square root of mean mass density, while individual frequencies probe the interior adiabatic sound speed and stellar buoyancy frequency profiles.

For radial modes, write ${\xi }_r=r\eta$. Eliminating the pressure and gravitational perturbations gives the radial stellar pulsation equationThis is a Sturm-Liouville problem. Multiplication by $\eta$ and integration, with vanishing boundary terms, gives its Rayleigh quotientFor constant ${\Gamma }_1>4/3$, $P^{\prime }<0$ makes both numerator contributions nonnegative. At ${\Gamma }_1=4/3$, a homologous displacement is neutral; for constant ${\Gamma }_1<4/3$ the same trial displacement makes the quotient negative. For the uniform-density stellar model, $\eta =$ constant is an exact mode and ${\omega }^2=(3{\Gamma }_1-4)GM/R^3$. With varying ${\Gamma }_1$, the integral criterion, rather than a universal pointwise threshold, controls radial stability.

The Cowling approximation neglects ${\varphi }_{\ell }$ while retaining the equilibrium gravitational field. It is useful for short-wavelength modes, but is not needed for the full derivation above. In a locally slowly varying region, take both remaining amplitudes proportional to $e^{i\int k_{r}dr}$ and retain the leading derivative terms. ThenEliminating either amplitude gives the acoustic-gravity propagation relationPositive $k_r^2$ is oscillatory propagation; negative $k_r^2$ means an evanescent wave. The high-frequency branch, ${\omega }^2>S_{\ell }^2,N^2$, describes stellar acoustic modes, restored chiefly by compressibility and pressure. At frequencies well above $N$, this gives ${\omega }^2\simeq c_s^2[k_r^2+\ell (\ell +1)/r^2]$. The low-frequency propagating branch in stable stratification, ${\omega }^2<N^2,S_{\ell }^2$, instead describes stellar gravity modes, restored by buoyancy. At $\ell =0$ there is no such nonradial gravity-wave cavity.

A stellar acoustic mode is trapped between an inner turning point near $\omega =S_{\ell }$ and an outer reflecting region. Low-$\ell$ modes penetrate deeply; radial modes reach the centre. Higher-degree modes turn farther out. Standing waves require the WKB quantization conditionwhere the phase $\alpha$ depends on the central or turning-point behaviour and surface reflection. For high radial order and small degree, the leading acoustic travel-time result is the large frequency separationThe $\ell /2$ term is the leading central angular phase shift; smaller frequency separations depend on detailed interior gradients. The near-surface mass density stratification sets the acoustic cutoff frequency. In a plane-parallel isothermal atmosphere with mass density scale height $H_{\rho }$, ${\omega }_{\mathrm{ac}}\simeq c_s/(2H_{\rho })$. Modes below this cutoff can reflect and form a cavity; waves sufficiently above it escape and need an outgoing-wave boundary condition. Nonadiabatic damping, driving and rotation alter real-star mode properties, but the adiabatic frequency problem isolates their dependence on the equilibrium structure.

Substitute part (b) into the spectral representation. The terms containing $\psi (ik)$ areThe transform $\psi (ik)$ is analytic for $\mathrm{Im}k>0$. Write $k=a+ib$. In the first quadrant,which decays since $x,y>0$. The Cauchy integral theorem in that quadrant equates the two integrals in the first parentheses. In the second quadrant,which decays since $a<0$, $b>0$, and $y<\ell$. The same contour integration equates the two integrals in the second parentheses. Therefore $\psi (ik)$ makes no contribution. This is upper-quadrant cancellation of reflected boundary transforms, and it uses the whole paired expression rather than attempting to discard an individual integral.

The terms involving $c$ cancel by the same two quadrant arguments with the analytic factors $k\pm \gamma$ and $\psi$ omitted. There is a sign defect in the printed reconstruction prefactor. The sides specified in part (a) form a clockwise boundary, so the Cauchy integral theorem gives a negative reconstruction sign. Indeed, integrating each spectral ray first yields $i/(z-\zeta )$, and hence the sum of the three ray integrals is $i{\oint }_{\mathrm{clockwise}}{q}_z(\zeta )/(z-\zeta )d\zeta =-2\pi q_z(z)$. Thus the printed positive prefactor reconstructs $-q_z$. The requested cancellation above holds with either overall sign, but the actual data-dependent derivative isDecay at infinity fixes the integration constant when recovering $q$ from its Wirtinger derivative.

For completeness, the missing condition at infinity and the unrestricted printed $\gamma$ deserve an explicit solvability of a decaying Robin strip check. Let $Y_n$ be orthonormal eigenfunctions from a Sturm-Liouville problem of $-d^2/d{y}^2$ with $Y_n^{\prime }(0)=\gamma Y_n(0)$ and $Y_n^{\prime }(\ell )=-\gamma Y_n(\ell )$, and eigenvalues ${\lambda }_n$. With $g_n={\int }_0^{\ell }g{Y}_{n}dy$, separation of variables givesIndeed, each coefficient obeys $q_n^{\prime \prime }-{\lambda }_n{q}_n=0$ and $q_n^{\prime }(0)=g_n$. A positive eigenvalue has one decaying exponential; a zero eigenvalue has only affine solutions and a negative eigenvalue only oscillatory solutions, neither of which decays unless its coefficient vanishes. This also proves uniqueness in the decay class and justifies reflection symmetry there. For $\gamma >0$, all eigenvalues are positive. For $\gamma =0$, the constant eigenfunction requires ${\int }_0^{\ell }gdy=0$; the apparent zero of $D$ at the origin is then removable. For $\gamma <0$, compatibility with every nonpositive eigenmode is necessary. These restrictions cannot be inferred from smoothness and reflection symmetry alone. The printed boundary value problem without any far-field condition permits additional growing or nondecaying homogeneous solutions, whereas the spectral transforms select the decaying branch.

For the real coefficients of the stated self-adjoint Sturm-Liouville problem, define the weighted inner productSince $w>0$, a nonzero eigenfunction has positive squared norm. Let $Ly=-(p{y}^{\prime }{)}^{\prime }+qy$. Integration by parts and the Dirichlet boundary conditions giveTaking $u=v$ to be an eigenfunction shows that its eigenvalue is real, since $\lambda \int w|u{|}^2=\int (p|u^{\prime }{|}^2+q|u{|}^2)$ is real. For two eigenfunctions with eigenvalues $\lambda ,\mu$, the same identity givesHence distinct eigenvalues imply orthogonality in the weighted inner product. This is the orthogonality of Sturm-Liouville eigenfunctions. The calculation assumes the regularity needed for the stated self-adjoint boundary problem, so all boundary terms and integrals are defined.

Multiply the equation by $x^2$ and combine the derivative terms:Thus the self-adjoint form isThe Dirichlet boundary conditions remain at $1$ and $e$, and the weight $w=x$ is positive on this interval, as required for the Sturm-Liouville problem.

Assume the integrand is sufficiently smooth near the positive admissible profiles. Define the second variation as ${\delta }^{2}F[y,\xi ]={d^{2}F[y+\epsilon \xi ]/d{\epsilon }^2|}_{\epsilon =0}$, so it is twice the quadratic Taylor coefficient. Differentiation under the integral givesIntegrating the middle term by parts, with $\xi (\alpha )=\xi (\beta )=0$, yieldsAll coefficients are evaluated along $y$.

For the surface-area integrand $f=2\pi y\sqrt{1+p^2}$, the Beltrami identity for a stationary profile readsEquivalently, its Euler-Lagrange equation is $y{y}^{\prime \prime }=1+y^{\prime 2}$, solved by the positive catenoid profile $y=E\cosh ((x-x_0)/E)$ with $E>0$. Equal boundary radii at $x=\pm L$, for $L>0$, require $x_0=0$, soDirect substitution also verifies the Euler-Lagrange equation. This is a stationary profile whenever a positive $E$ satisfying the boundary relation exists; such an $E$ does not exist for arbitrary ring separations and radii.

On this profile, with $z=x/E$, the coefficients are $f_{yy}=0$, $f_{pp}=2\pi E{\mathrm{sech}}^{2}z$ and $d{f}_{yp}/dx=(2\pi /E){\mathrm{sech}}^{2}z$. Changing variables therefore givesThe printed strict inequality must mean every nonzero admissible variation: the zero function makes the integral zero, so the literal quantifier including it is false.

For completeness, the strict positivity here is sufficient for a strict local minimum in the fixed-endpoint $C^1$ topology, among positive axisymmetric graph profiles. On this finite interval the weight $p(z)={\mathrm{sech}}^{2}z$ has a positive minimum. The regular Sturm-Liouville problem associated with $Q[\xi ]=\int p({\xi }_z^2-{\xi }^2)$ has an attained lowest Dirichlet eigenvalue ${\lambda }_1$ in the Rayleigh quotient. Its first eigenfunction is a nonzero smooth admissible variation; the assumed strict positivity gives ${\lambda }_1>0$. Hence $Q\ge {\lambda }_1\Vert \xi {\Vert }_{L^2}^2$. Combining this with $Q\ge p_{\min }\Vert {\xi }_z{\Vert }_{L^2}^2-p_{\max }\Vert \xi {\Vert }_{L^2}^2$ proves coercivity in the Sobolev space $H_0^1$.

The coefficients of the second variation depend continuously on $(y,y^{\prime })$ while $y>0$. For a sufficiently small $C^1$ perturbation, the quadratic forms along the segment from the stationary profile to the perturbed one retain a uniform positive coercivity bound. Taylor's formula with integral remainder, and the vanishing first variation, then gives $F[y+v]-F[y]\ge c\Vert v{\Vert }_{H^1}^2>0$ for nonzero sufficiently small admissible $v$. The stationary area is therefore a strict local minimum under the stated nonzero-variation condition. This supplies the uniform bound needed in an infinite-dimensional problem, rather than relying on pointwise positivity alone in an arbitrary function space.

Use the weighted Neumann heat kernel with constant drift. The spatial operator has the self-adjoint formso its homogeneous Neumann boundary conditions give a regular Sturm-Liouville problem in the weighted inner productThe constant eigenfunction is ${\varphi }_0=1$, with eigenvalue ${\lambda }_0=0$ and squared norm $N_0=(e^{\beta L}-1)/\beta$. For $k\ge 1$, put $q_k=k\pi /L$. Substitution of $\varphi =e^{-\beta x/2}v$ in $\mathcal{L}\varphi =-\lambda \varphi$ gives $v^{\prime \prime }+(\lambda -{\beta }^2/4)v=0$ and the Robin boundary conditions $v^{\prime }=\beta v/2$ at both endpoints. HenceIndeed,The Sturm-Liouville eigenfunction expansion is complete and orthogonal. Its heat kernel, relative to the measure $e^{\beta y}dy$, isFor a coefficient $a_k(t)=⟨{\varphi }_k,u(\cdot ,t){⟩}_{\beta }/N_k$, two applications of integration by parts yieldSolving these scalar linear differential equations and summing gives the Neumann boundary-forcing formula:This depends only on the given data. In particular, the zero eigenvalue retains the changing weighted mean:The endpoint derivatives of the boundary integral are understood as interior limits. Differentiating each homogeneous Neumann eigenfunction at the endpoint before summing would incorrectly discard the prescribed boundary forcing: the short-time heat kernel makes that interchange invalid. The stated compatibility conditions give the classical boundary traces.

LetTheir Bessel equations in self-adjoint form areMultiply the first by $v$, the second by $u$, subtract, and integrate. The left side becomes a boundary term:Regularity at zero makes the lower boundary term vanish, whileTherefore

If $\alpha ={\gamma }_k$ and $\beta ={\gamma }_{\ell }$ with $k\ne \ell$, both endpoint values of $J_0$ vanish, so the integral is zero. This is the weighted orthogonality of distinct eigenfunctions in a Sturm-Liouville problem. For the norm, hold $\alpha ={\gamma }_k$ and let $\beta \to {\gamma }_k$. Sincethe quotient tends to

Multiplication by the integrating factor $e^{2x}$ givesHence the Sturm-Liouville problem is

For a radial stellar oscillation, put ${\xi }_r=r\eta$. The self-gravitating adiabatic equation is $({\Gamma }_{1}P{r}^4{\eta }^{\prime }{)}^{\prime }+r^3[(3{\Gamma }_1-4)P{]}^{\prime }\eta +\rho r^4{\omega }^2\eta =0$. It is a Sturm-Liouville problem whose Rayleigh quotient tests radial stability. The familiar constant-exponent threshold is ${\Gamma }_1=4/3$; a variable exponent requires the full integral test.