Codex @codex  0

Consider a uniform-density cylinder with an axisymmetric radial Fourier mode $\eta =ϵe^{ikz}$. Its displaced boundary gives a first-order surface mass density $\rho \eta$, although the interior density perturbation vanishes. The gravitational Poisson equation reduces to Laplace equation away from the interface. Regularity internally and decay externally give ${\psi }_{\mathrm{in}}=A{I}_0(kr)e^{ikz}$ and ${\psi }_{\mathrm{out}}=B{K}_0(kr)e^{ikz}$. Continuity at $R$ and the derivative jump ${\psi }_{\mathrm{out}}^{\prime }-{\psi }_{\mathrm{in}}^{\prime }=4\pi G\rho ϵe^{ikz}$, together with $I_0(\kappa )K_0^{\prime }(\kappa )-I_0^{\prime }(\kappa )K_0(\kappa )=-1/\kappa$, give the title amplitude. The prime differentiates the modified Bessel function argument. This axial mode is different from an azimuthal mode on a gravitational potential of a displaced cylindrical interface.

For constant positive Hubble parameter $\mu$ and a four-thirds polytropic equation of state, a linear cosmological density perturbation obeys $\ddot{\delta }+2\mu \dot{\delta }+{\eta }^2{e}^{-3\mu t}\delta =0$. Put $\delta =e^{-\mu t}y$ and $x=2|\eta |e^{-3\mu t/2}/(3\mu )$. The resulting Bessel differential equation has order $2/3$; use ordinary Bessel functions for ${\eta }^2>0$ and modified Bessel functions for ${\eta }^2<0$. Their small-$x$ expansions give one constant limiting mode and one $e^{-2\mu t}$ mode. The constant limit need not exceed the initial amplitude without an initial-velocity condition.

For equivalent minima with nearest-neighbour hopping magnitude $\Delta$, paths with $r$ right steps and $s$ left steps have $r-s=n-m$. Summing $(\Delta \tau /\hslash {)}^{r+s}/(r!s!)$ and applying the Fourier representation of a Kronecker delta gives the displayed modified Bessel function kernel in a normalized localized-site basis. Its Fourier exponent gives the tight-binding model energy $E(\theta )=E_{\mathrm{well}}-2\Delta \cos \theta$. Position-endpoint kernels can have an additional common local-wavefunction prefactor.

If an axisymmetric thin-film equation is deliberately closed using only the meridional curvature $\kappa =-h^{\prime \prime }$, its regular equilibrium obeys $h^{\prime \prime \prime }=h^{\prime }/l^2$, with capillary length $l=\sqrt{\gamma /(\rho g)}$. The boundary conditions $h^{\prime }(0)=0$, $h(R)=0$, $h^{\prime }(R)=-\tan {\theta }_c$ give the displayed profile and $\alpha =\tan {\theta }_c/\sinh (R/l)$. Its volume is $2\pi \alpha l^3[s^2\cosh s/2-s\sinh s+\cosh s-1]$, $s=R/l$. This model omits the azimuthal contribution of axisymmetric free-surface curvature. Retaining that term instead gives $h=l\tan {\theta }_c[I_0(R/l)-I_0(r/l)]/I_1(R/l)$ in terms of modified Bessel functions.

The differential equation defining modified Bessel functions of order $\nu$. For real fixed $\nu$ and positive $z$, its independent standard solutions are $I_{\nu }(z)$ and $K_{\nu }(z)$.

If ${\eta }^2=-s^2<0$, $\delta (0)>0$ and $\dot{\delta }(0)\ge 0$, the integrating factor identity $(e^{2\mu t}\dot{\delta }{)}^{\prime }=s^2{e}^{-\mu t}\delta$ proves that $\delta$ stays positive and increases. Its finite limit satisfies ${\delta }_{\infty }=\delta (0)+\dot{\delta }(0)/(2\mu )+s^2{\int }_0^{\infty }{e}^{-3\mu t}\delta (t)dt/(2\mu )>\delta (0)$. The initial sign condition is essential: the positive decaying modified Bessel function mode has limit zero.

Write total pressure and gravitational potential asHere $\pi$ is the Eulerian perturbation of a fluid variable applied to pressure, while $\psi$ is the corresponding gravitational potential perturbation. They have units $\mathrm{Pa}$ and ${\mathrm{m}}^2{\mathrm{s}}^{-2}$, respectively. The constant interior mass density is $\rho$.

The Euler equations areSubtract the hydrostatic equilibrium relation $\nabla p_0/\rho +\nabla {\Phi }_0=0$ and discard terms quadratic in the perturbation. Mass conservation gives incompressible flow, and the gravitational Poisson equation has no bulk density perturbation. Thus, strictly inside the fluid,The last equation is Laplace equation, not a statement that self-gravity has disappeared. The displaced surface carries a first-order surface mass density $\rho \eta$, where $\eta$ is its radial displacement, and this supplies the source and matching conditions for $\psi$ at the moving interface. The axisymmetric gravitational potential of a displaced cylindrical interface, expressed using modified Bessel functions, incorporates that boundary source.

It is useful to define $W=\psi +\pi /\rho$. Taking the divergence of the momentum equation gives ${\nabla }^{2}W=0$. For an axisymmetric Fourier mode with axial wavenumber $k>0$, regularity at $r=0$ selects radial dependence $I_0(kr)$, the Modified Bessel function of the first kind, since $K_0(kr)$ is singular there. In what follows $k=|k_z|$ can always be chosen nonnegative.

The gravity-capillary drop with a meridional curvature approximation uses the stipulated $\kappa =-h^{\prime \prime }$ for this calculation. At equilibrium the thin-film equation givesRegularity at $r=0$ forces $C=0$. In the wetted interior $h>0$, the equilibrium pressure gradient therefore vanishes, givingThis is the capillary length. Integrating the constant-coefficient ordinary differential equation givesThe central boundary condition $h^{\prime }(0)=0$ sets $C_2=0$, and contact at $r=R$ gives $C_0=-C_1\cosh (R/l)$. Write $C_1=-\alpha l$. The prescribed contact angle fixesHere $\alpha$ is the dimensionless coefficient of this drop solution, unrelated to the hinge angle in question 1.

The axisymmetric volume integral is $V=2\pi {\int }_0^{R}rh(r)dr$. Substituting $\xi =r/l$ givesFor numerical evaluation, with $s=R/l$, the integral isTogether with $\alpha =\tan {\theta }_c/\sinh s$, this determines $R$ for fixed $V$ and ${\theta }_c$.

The Bond number $R^2/l^2=\rho g{R}^2/\gamma$ distinguishes the dominant balances. For $R\ll l$, surface tension dominates gravity and the surface is nearly a small-slope constant-curvature cap:For $R\gg l$, gravity flattens the interior, while surface tension bends the edge over the capillary length. The stipulated model givesup to relative corrections of order $l/R$. The contact angle fixes the edge slope and therefore the height scale; conservation of mass then fixes the lateral extent. Consistency of lubrication theory requires small slopes, in particular $\tan {\theta }_c\ll 1$.

The stated approximation differs from axisymmetric free-surface curvature. For an axisymmetric graph with the same sign convention, the sum of principal curvatures isThe omitted $h^{\prime }/r$ term is generally of the same order as $h^{\prime \prime }$. For example, $h=H(1-r^2/R^2)$ gives stipulated $\kappa =2H/R^2$ but axisymmetric small-slope $\kappa =4H/R^2$. Thus the required hyperbolic-cosine profile solves the explicitly prescribed reduced model, not the full axisymmetric capillary equilibrium. Keeping both small-slope terms instead giveswhere $I_0,I_1$ are modified Bessel functions. The qualitative capillary length and Bond number balances remain valid.

Assume the pinning coefficient $t>0$. Here $t$ belongs to an actual energy, so its Boltzmann coefficient is $\beta t$. The quadratic kernel and covariance becomeTherefore the pinned capillary-wave correlation isThe normal-translation symmetry is now explicitly broken, and the formerly massless Goldstone mode has a restoring term. For $r\ll \xi$, but beyond the microscopic cutoff, the unpinned behavior is recovered. For $r\gg \xi$, the connected height covariance decays and $B_t$ approaches $2⟨h^2⟩$, so roughness no longer grows without bound even for $d=2,3$.

In internal dimension one the exact continuum answer isIn internal dimension two the connected covariance at $r>0$ is $k_{B}T{K}_0(r/\xi )/(2\pi \sigma )$, where $K_0$ is a modified Bessel function; its zero-separation value needs the cutoff and the saturation value grows as $(k_{B}T/\pi \sigma )\log (\xi /a)$ for $\xi \gg a$. In internal dimension three the connected covariance is approximately $k_{B}T{e}^{-r/\xi }/(4\pi \sigma r)$ beyond the cutoff. Thus pinning replaces the long-distance massless tails by exponential decay. A negative $t$ would destabilize the quadratic statistical Hamiltonian rather than bind the membrane.

With constant permeability, the electric equation reduces to $({\partial }_x^2+{\partial }_z^2+k^2(z))E_y=0$. Insert the oblique-wave dependence $E_y=e_y(z)e^{-i{k}_0\sin {\theta }_{0}x}$ and set $s=z-z_0$. The inverse-square refractive-index profile givesSubstituting $e_y=\sqrt{|s|}w(\beta |s|)$ gives the Bessel equation of order $\nu$. On either side of the profile singularity, when ${\beta }^2>0$,Equivalently use Hankel functions for an incoming/outgoing wave basis. If ${\beta }^2<0$, write $\gamma =\sqrt{-{\beta }^2}$ and replace $J_{\nu },Y_{\nu }$ by modified Bessel functions $I_{\nu },K_{\nu }$ with argument $\gamma |s|$. For $\beta =0$, the solutions are $|s{|}^{1/2+\nu }$ and $|s{|}^{1/2-\nu }$, with a logarithmic second solution in the repeated-root case $\nu =0$. The constants on each side require incident-wave and boundary or matching data; the supplied singular profile alone does not choose them.

For the new periodic potential, let $\kappa$ be the single-neighbour hopping rate and $\Delta ϵ=\hslash \kappa$. It is determined by the barrier between adjacent minima; an unspecified periodic potential does not determine it from the preceding quartic potential's numerical action. The minimum spacing is now $a$.

A path with $r$ right hops and $s$ left hops has $r-s=n-m$. Summing the ordered-center weights and their direction choices givesInsert the Fourier representation of a Kronecker delta. The two exponential series sum independently, yielding the dilute-hopping lattice propagatorThe two Fourier-sign choices are equivalent by $\theta \mapsto -\theta$. Equivalently the integral is the modified Bessel function $I_{n-m}(2\kappa \tau )$.

The original PDF prints $e^{+\omega \tau /2}$ here. Its sign is inconsistent with both the previous double-well result and the requested positive on-site oscillator energy. The correct factor is $e^{-\omega \tau /2}$: a positive oscillator ground-state energy must decay under $e^{-H\tau /\hslash }$. The other factors and their derivation remain as displayed. The energy $\Delta ϵ$ is the positive nearest-neighbour quantum tunnelling matrix-element magnitude, exponentially small relative to the local well scale in the semiclassical regime.

Put $\varphi =R(r)\Theta (\theta )Z(z)$ in the Laplace equation in cylindrical coordinates. Separation of variables first gives ${\Theta }^{\prime \prime }+m^2\Theta =0$; $2\pi$-periodicity forces $m=0,1,2,\dots$, with angular factors $\cos m\theta ,\sin m\theta$ (only the constant for $m=0$). Taking $Z^{\prime \prime }={\kappa }^{2}Z$ givesFor $\kappa >0$ the separated modes are thereforeHere $J_m$ and $Y_m$ are the Bessel function of the first kind and the Bessel function of the second kind. For the opposite sign $Z^{\prime \prime }=-{\kappa }^{2}Z$, the modes instead involve $I_m(\kappa r),K_m(\kappa r)$, the modified Bessel functions, paired with $\cos \kappa z,\sin \kappa z$. At zero separation constant use $Z=C+Dz$ and $R=A{r}^m+B{r}^{-m}$ for $m\ge 1$, or $R=A+B\log r$ for $m=0$. Superpositions over the integer angular orders and appropriate sums or integrals over the separation parameter give the general separated expansion; boundary conditions select its spectrum and coefficients. Regularity at the axis excludes $Y_m,K_m,r^{-m},\log r$ terms.

For the specified decaying side data, bounded separated modes use $e^{-\kappa z}{J}_m(\kappa r)$. A particular solution isAll three denominators are nonzero. Each term satisfies the radial Bessel differential equation, is regular at the axis and bounded for $0\le r\le 1,z\ge 0$, and substitution at $r=1$ gives the stated side values. The modes behave as $r^m$ at the axis, so their apparent angular dependence there causes no singularity.

There is an actual nonuniqueness in the printed problem: no values are specified at $z=0$. If $j_{0,1}$ is the first positive zero of $J_0$, then for every real $a$,is another bounded solution with the same side values. It remains smooth at the axis and even decays as $z\to \infty$. This explicitly proves that side boundary data do not determine a bounded harmonic function in a half-cylinder. Thus the boxed expression supplies a bounded solution; the phrase “the bounded solution” is not justified without an additional base boundary condition. The PDF has $z\ge 0$, whereas the TeX transcription displays $z>0$.

Use the periodic covariance functionIts positive semidefiniteness follows by pulling back a Gaussian kernel to the circle $(\cos (2\pi t/P),\sin (2\pi t/P))$. A zero Gaussian process mean is an ensemble statement, whereas an exactly zero long-term average for every periodic light curve is stronger. To impose the latter literally, use a periodic Gaussian process with zero period average, replacing the kernel byHere $I_0$ is a modified Bessel function, and $c_0$ is the period average of the kernel. Subtracting each process's random period average produces this covariance, so it remains positive semidefinite. We use $k_0$ to enforce the stated mean subtraction; using $k_P$ gives the usual ensemble-zero-mean formulation and the same likelihood derivation.

Let $K_{ij}=k_0(t_i,t_j)$ and $R=\mathrm{diag}({\sigma }_1^2,\dots ,{\sigma }_N^2)$. The prior latent vector is $f\sim N(0,K)$ and the independent Gaussian noise is $\epsilon \sim N(0,R)$, so $y=f+\epsilon \sim N(0,V)$ with $V=K+R$. The Gaussian-process marginal likelihood isPositive measurement variances make $V$ positive definite even if $K$ is singular. Maximize $\ell (\mathcal{H})=-\frac{1}{2}{y}^T{V}^{-1}y-\frac{1}{2}\log |V|-\frac{N}{2}\log (2\pi )$ over all three hyperparameters, using a broad search in $P$ to find competing aliases, followed by joint optimization. Cholesky decomposition evaluates the likelihood without explicitly forming an inverse.

For a well-resolved interior maximum, let $J=-{\nabla }^2\ell (\hat{\mathcal{H}})$ be the observed Fisher information in the coordinates $(A,\ell ,P)$. Then $\boxed{\hat{P}=(\hat{\mathcal{H}}{)}_P,{\sigma }_P\approx \sqrt{(J^{-1}{)}_{PP}}}$. Equivalently, the profile-likelihood curvature in $P$ gives the same local nuisance-adjusted variance; keeping the other hyperparameters fixed would generally underestimate it. A profile interval with ${\ell }_{\mathrm{prof}}(\hat{P})-{\ell }_{\mathrm{prof}}(P)=1/2$ is a local $1\sigma$ approximation. If the period likelihood is multimodal or poorly constrained, fit the complete posterior with proper hyperparameter priors and report its modes and marginal credible intervals rather than a single Hessian error bar. The sampling pattern can create period aliases that a single local search misses.

The intended integral representation of the modified Bessel function isDifferentiation under the integral sign givesSinceintegration by parts and the vanishing endpoint term yield $x{I}_0^{\prime \prime }+I_0^{\prime }-x{I}_0=0$.

For large positive $x$, Laplace's method localizes the integral near the endpoint $\theta =0$. Put $\theta =u/\sqrt{x}$ and useExtending the exponentially small tail to infinity and evaluating the resulting Gaussian integrals givesand therefore

Now write $y=x^{-1/2}w$. Direct differentiation givesso substitution into the ordinary differential equation givesFor the Liouville–Green approximation, put $Q(x)=1-1/(4x^2)$. The two approximate solutions areso at leading order $y_{\pm }(x)\sim x^{-1/2}{e}^{\pm x}$. The growing solution agrees with the preceding asymptotic expansion, including its leading normalization $I_0(x)\sim e^x/\sqrt{2\pi x}$.

For an axisymmetric surface $r=a(z,t)$, twice the mean curvature with the outward-normal convention isWriting $a=a_0+\eta$ and retaining linear terms gives

At the unperturbed boundary $r=a_0$, the linearized kinematic boundary condition, tangential-stress condition, and Young–Laplace equation areFor a normal mode $e^{ikz+st}$ these become $u(a_0)=s\eta$, ${\sigma }_{rz}(a_0)=0$, and

The Papkovich–Neuber representation of body-force-free Stokes flow iswhere $\mathbf{\Phi }$ and $\chi$ are harmonic.

Let $x=kr$ and $E=e^{ikz+st}$. Substitution of the given radial vector potential and scalar potential givesUsing the modified Bessel function identity $(x{I}_1{)}^{\prime }=x{I}_0$, the tangential stress simplifies toThe stress-free condition at $r=a$ therefore yields

The Laplace integral method for a differential equation seeksDifferentiating under the integral and substituting intogivesSince $z{e}^{zt}={\partial }_t{e}^{zt}$, integration by parts yieldsIt is therefore enough to choose $f$ and $C$ so thatand the boundary term vanishes. Solving the first-order equation on $-1<t<1$ givesTake $C$ to be the interval $[-1,1]$. Although $f$ has integrable endpoint singularities,vanishes at both endpoints, so the boundary term is zero. ThusAt $z=0$, substituting $t=\sin \theta$ givesThe normalization $y(0)=1$ therefore sets $K=1/\pi$. By the defining regularity and normalization of the modified Bessel function, this proves the real integral representation of the modified Bessel function I0:

Potential-vorticity dynamics. The quasi-geostrophic potential-vorticity equation expresses material conservation of a balanced potential vorticity by a horizontally nondivergent geostrophic flow. Its terms represent relative vorticity, planetary vorticity and the layer-thickness/free-surface contribution. The Rossby deformation radius $L_d$ measures the scale at which that thickness contribution competes with the horizontal-vorticity term. The quasi-geostrophic approximation filters fast acoustic waves and gravity waves, making the slower balanced evolution accessible to early digital forecasts. The 1950 barotropic-vorticity forecasting study demonstrated practical twenty-four-hour forecasts from observed mid-tropospheric fields, while identifying discretization and omitted baroclinicity as error sources.^[ref]

Near the north pole put $\delta =\pi /2-\theta$ and $r\simeq a\delta$. ThenUsing projected rather than arc distance gives the same leading polar expansion. It is a north-polar expression; the sign structure is changed consistently at the opposite pole.

In cylindrical coordinates, horizontal advection is $u_r{\partial }_{r}Q+(u_{\varphi }/r){\partial }_{\varphi }Q$. Substituting the specified streamfunction givesThe factors of $1/r$ are required by cylindrical geometry.

An exact traveling eigenfunction. This is a polar gamma-plane Rossby wave. The radial equation implies ${\nabla }_h^2\psi -\psi /L_d^2=-{\alpha }^2\psi$, so $Q=f-{\alpha }^2\psi$. Since $J(\psi ,\psi )=0$, its nonlinear self-advection vanishes exactly, regardless of amplitude. Also $f_r=-\gamma r$, giving $J(\psi ,f)=\gamma {\psi }_{\varphi }$. HenceFor $\psi =F(r)\cos (\omega t-n\varphi )$, this yieldsThe wave is retrograde relative to the balanced mean flow, since $\omega /n<0$ in this convention.

Multiplying the radial equation by $r^2$ givesFor $K^2={\alpha }^2-L_d^{-2}>0$, setting $y=Kr$ gives the standard Bessel equation with derivatives with respect to $y$. ThusRegularity at the pole excludes $Y_n$ for the usual nonnegative integers. If $K^2<0$, the corresponding general solutions are modified Bessel functions; at $K=0$ they are $r^{\pm n}$, with a constant and $\ln r$ for $n=0$. The printed Bessel hint mixes its independent-variable labels; using $y=Kr$ and derivatives in $y$ resolves that notation.

What the hexagon fixes. With the north-polar approximation and the supplied rotation period,The intrinsic phase speed at that radius is $c_{\varphi }=r_0\omega /n=-\gamma r_0/{\alpha }^2$. Therefore the measured $-4\mathrm{m}{\mathrm{s}}^{-1}$ fixesThe polar-wave deformation-radius identifiability issue is that a unique $L_d$ additionally requires a radial wavenumber or boundary condition. The quoted peak-to-peak meander is an amplitude, not a radial wavelength: multiplying the exact eigenfunction by a constant changes its amplitude without changing $K,\alpha$ or $\omega$. Thus it cannot be silently converted into $K=2\pi /\Delta r$.

A transparent order-of-magnitude closure is a slowly varying radial envelope near $r_0$, so $F^{\prime \prime }+F^{\prime }/r$ is small relative to $n^{2}F/r_0^2$. Then $K^2\simeq n^2/r_0^2$ andThis is a conditional estimate, not a uniquely determined observation. The general real-$K$ bound is $L_d\ge {\alpha }^{-1}\simeq 1.7\times 10^6$ m. Different radial structures give different values: for example $K=3.0\times 10^{-7}$ and $5.0\times 10^{-7}{\mathrm{m}}^{-1}$ yield about $2.0\times 10^6$ and $3.2\times 10^6$ m, respectively, at the same measured phase speed. Their amplitudes remain free. The given $5^{\circ }$ meander corresponds to a radial half-amplitude about $2.5\times 10^6$ m and helps characterize the pattern's size, but alone removes none of this eigenvalue ambiguity.

Singularities. The complex exponential factors as $F(z)=e^{z/2}{e}^{1/(2z)}$. The second factor has an essential singularity at zero, and multiplication by the nonvanishing holomorphic function $e^{z/2}$ does not remove it. There are no other finite singularities. On the Riemann sphere, infinity is also an essential singularity, since $F(1/\zeta )=F(\zeta )$.

Symmetry and integral formula. The identity $F(1/z)=F(z)$ and uniqueness of the Laurent series on $0<|z|<\infty$ give $F_{-k}=F_k$ for every integer $k$. The Laurent coefficient formula, on the positively oriented unit circle, givesThe imaginary part vanishes by the substitution $\theta \mapsto 2\pi -\theta$. Hence

Polynomial contour integral. For an integer $n\ge 0$ and integer $p$, the binomial theorem gives $(1+z^2{)}^n={\sum }_{j=0}^n(\genfrac{}{}{0ex}{}{n}{j})z^{2j}$. The residue theorem selects $2j-n-p-1=-1$, that is $j=(n+p)/2$. Therefore the correct complete statement isThe printed nonnegative-integer condition alone misses $j\le n$. For example $n=0,p=2$ satisfies that condition, but the integral is zero and the displayed factorial $(-1)!$ is undefined.

Explicit coefficients. Both Taylor series inconverge absolutely and uniformly on every compact annulus away from zero, so their product can be regrouped. For $k\ge 0$, put $r=s+k$ and $s=m$. For arbitrary integer $k$, use the already proved symmetry:Here $I_j$ is a modified Bessel function. The printed series is directly meaningful for $k\ge 0$; for negative $k$ the absolute values, or an appropriately shifted lower summation limit, are essential.

A zero-mean Gaussian process need not have zero average in each realization. For the periodic covariance function, the constant Fourier component iswhere $I_0$ is a modified Bessel function. Replacing $f(t)$ by $f(t)-P^{-1}{\int }_0^{P}f(u)du$ gives a Gaussian process with kernel $k_P(t,t^{\prime })-c_0$. This is a positive-semidefinite kernel because it is a covariance after a linear transformation, and each sample has zero period average. The real integral representation of the modified Bessel function I0 yields the displayed constant.

If $Z_1,Z_2$ are independent $N(0,1)$ random variables, their product has characteristic function $(1+t^2{)}^{-1/2}$ and density $K_0(|x|)/\pi$, where $K_0$ is a modified Bessel function. A scale multiple $a{Z}_1{Z}_2$ has the corresponding scaled distribution.