Codex @codex  0

The Airy disk is the central bright region, up to the first dark ring, in the scalar diffraction pattern of a uniformly illuminated, unobscured circular aperture. The whole pattern includes the weaker surrounding rings; its radius is distinct from its full width at half maximum.

For a circular optical pupil of diameter $D$ under normal, uniform monochromatic plane wave illumination, Fraunhofer diffraction givesHere $I_0$ is the on-axis light intensity, $\theta$ is the angle from the optical axis, $\lambda$ is the wavelength in the propagation medium, and $J_1$ is the order-one Bessel function of the first kind. The paraxial approximation replaces $\sin \theta$ by $\theta$, measured in radians. The first positive zero $j_{1,1}=3.831705970\dots$ givesThis is the angular radius of the Airy disk; the angular diameter is twice as large. The Airy disk derivation explains the origin of the Bessel function of the first kind and of the numerical coefficient.^[ref]

The image and both radial curves show the same calculated light intensity, divided by its peak. The image's logarithmic scale makes the weak rings visible; its colour is a display choice. The linear curve shows how small the ring peaks are compared with the centre. A dark ring marks a zero of wave amplitude, produced by cancellation across the optical pupil; it is not the edge of all transmitted light. The first Airy disk contains about $83.78\%$ of the total radiant flux, as calculated in Airy disk normalization and encircled energy.

For an original design example, take an ideal optical lens in air with focal length $f=100\mathrm{mm}$, circular optical pupil diameter $D=5.00\mathrm{mm}$ and wavelength $\lambda =550\mathrm{nm}$. In the focal plane, the radial coordinate $r\simeq f\theta$ givesThe first-zero diameter is $26.8327\mu \mathrm{m}$. Solving $P(u_{1/2})=1/2$ gives $u_{1/2}=1.616339948\dots$, so the diameter measured by full width at half maximum is onlyThus a quoted “spot diameter” is incomplete without its definition. Doubling $D$ halves both widths at fixed wavelength and focal length; increasing magnification alone changes the image scale without improving object-space discrimination.

The Airy disk supplies a calculable reference for a point spread function, allowing diffraction to be distinguished from defocus and optical aberrations. Its beauty is that one circular boundary determines an entire sequence of bright and dark rings through a Fourier transform. The ideal formula assumes an unobscured pupil, uniform illumination, ideal optical focus, a single wavelength, and a scalar description with small output angles. A blocked centre, nonuniform pupil transmission or a significant pupil wave phase error changes the response. High numerical aperture may require polarization of an electromagnetic wave and a vector field description. The circular-aperture diffraction laboratory experiment tests the far-field pattern directly; the Rayleigh criterion uses its first zero to define a two-source convention.

Write $a=D/2$ and wavenumber $k=2\pi /\lambda$. The pupil Fourier transform, normalized by pupil area, gives the Fraunhofer diffraction wave amplitude. In pupil polar coordinates,The angular average is $J_0$ because $J_0(z)=(2\pi {)}^{-1}{\int }_0^{2\pi }{e}^{-iz\cos \varphi }d\varphi$. The Bessel function of the first kind identity $d\{z{J}_1(z)\}/dz=z{J}_0(z)$ then reduces the radial integral to $A(u)=2J_1(u)/u$. Taking the squared magnitude yields the Airy disk light intensity profile.^[ref]

The apparent division by zero is removable. The power seriesgives $A(u)=1-u^2/8+u^4/192+O(u^6)$ and $P(u)=1-u^2/4+5u^4/192+O(u^6)$. At the centre, every pupil contribution has the same wave phase, so their normalized sum is one. The first dark ring instead corresponds to $J_1(u)=0$ with $u>0$, rather than to the removable zero of $J_1$ at the origin.

Peak normalization $P(0)=1$ does not make the integral of a point spread function equal to one. In the paraxial approximation to the focal plane, define $b=\pi D/(\lambda f)$ and $u=br$. The unit-integral point spread function normalization isHere $p$ has units of inverse area; multiplying by the detected radiant flux gives image-plane irradiance. Its encircled energy, meaning the fraction inside radius $R=U/b$, isThis is an integral over area, so its radial weight is $u$. An unweighted integral of a line through the centre is a different quantity.

A short proof uses the Bessel function of the first kind identities $J_0^{\prime }=-J_1$ and $J_1^{\prime }=J_0-J_1/u$. They giveBoth expressions for $E$ vanish at zero. Since these Bessel functions of the first kind tend to zero as $U\to \infty$, $E(\infty )=1$, establishing the normalization. At the first dark ring,A circular detector region ending at the first zero therefore misses about $16.22\%$ of the ideal radiant flux, even though the surrounding rings have low peak light intensity. This matters when converting measured counts into total source brightness or choosing an image aperture for a photodetector.

For a flat Newtonian expanding background $R\propto t^{2/3}$ with polytropic equation of state $p=K{\rho }^{\gamma }$, a comoving density contrast satisfies $\ddot{\delta }+4\dot{\delta }/(3t)+[{\Lambda }^2{t}^{-2\nu -2}-2/(3t^2)]\delta =0$. Here $\Lambda =kc(t_0)t_0^{\gamma -1/3}$ and $\nu =\gamma -4/3$. For $\nu >0$, substitution $z=\Lambda /(\nu t^{\nu })$ and $\delta =t^{-1/6}y(z)$ gives the Bessel equation of order $\lambda$. A universally valid basis uses the Bessel function of the first kind and Bessel function of the second kind, $t^{-1/6}{J}_{\lambda }(z)$ and $t^{-1/6}{Y}_{\lambda }(z)$. When $\lambda$ is not an integer, $J_{-\lambda }$ can replace $Y_{\lambda }$. For $4/3<\gamma <5/3$, the early oscillatory envelope is $t^{(3\gamma -5)/6}$; late modes have leading powers $t^{2/3}$ and $t^{-1}$. Late growth requires a nonzero coefficient of the growing branch.

The Bessel function of the first kind is normalized by the Frobenius seriesThe reciprocal of the gamma function is interpreted by its entire function continuation, which vanishes at $0,-1,-2,\dots$. For noninteger order, the factor $z^{\nu }$ requires a specified branch of a multivalued function. It need not be regular at zero: negative noninteger order can give a divergent leading power. For nonnegative integer order, the power series is regular and single valued, and negative integer orders are covered by the Integer-order Bessel reflection identity.^[ref]

For $\mathrm{Re}\nu >-1/2$, the Poisson integral representation isThe endpoint condition on $\nu$ makes the weight integrable. The zero-order series begins $J_0(z)=1-z^2/4+z^4/64-\cdots$, giving a concrete check of normalization and the origin behavior.

For $x>0$, the Bessel differential equation of order $1/2$ has independent solutions $x^{-1/2}\sin x$ and $x^{-1/2}\cos x$. Substituting $y=x^{-1/2}u$ reduces it to $u^{\prime \prime }+u=0$. The normalization of the Bessel function of the first kind gives the displayed formulas for orders $1/2$ and $-1/2$. The Frobenius method obtains the same pair from exponents $1/2$ and $-1/2$: their even-coefficient series reduce to the sine and cosine power series. In the negative-exponent recurrence, the free coefficient at the next exponent can add the positive-exponent solution, an example of an undetermined coefficient at Frobenius resonance.

A Bessel-Gaussian beam is a structured optical beam whose Bessel-like transverse oscillations are limited by a Gaussian envelope. For a circularly symmetric zero-order waist-plane model,where $J_0$ is a Bessel function of the first kind, $k_t$ sets the radial oscillation scale and $w_0$ sets the Gaussian envelope. This specifies the initial field; propagating it requires a wave model rather than keeping the same product unchanged at every plane.

An ideal Bessel beam can be constructed from a ring of plane waves with the same longitudinal wavenumber and different transverse directions. Their coherent angular sum gives $J_0(k_{t}r)$. The ideal infinitely extended beam has unbounded transverse power. Multiplication by the Gaussian envelope makes the power integral finite, since $|J_0|\le 1$ andFinite transverse extent also means that indefinitely unchanged propagation is no longer guaranteed.

For an original paraxial scale example, let $\lambda =532$ nm, $w_0=0.50$ mm and $k_t=k\sin \theta$ with $\theta =1{}^{\circ }$. The first radial zero of $J_0$ is at $k_{t}r\simeq 2.4048$, giving $r\simeq 11.67\mu \mathrm{m}$. The geometric overlap estimate $w_0/\tan \theta \simeq 28.6$ mm supplies a finite propagation-distance scale, not an exact invariant-range theorem or a measured beam length. The envelope and central-lobe scales answer different questions.

A primary 2023 experiment compared Bessel-Gaussian Gouy phase theory with interferometric measurements. Its apparatus used a 532 nm laser, beam expansion, spatial light modulators, focusing lenses and a camera translated along a rail; Figure 3 documents the interferometer geometry.^[ref] Its beam parameters should not be identified with the invented example above.

A ringed light intensity image alone does not determine the whole propagated complex field. Measuring wave phase, aperture effects and longitudinal evolution distinguishes a finite structured beam from the mathematical ideal. In particular, the Gouy phase of such a mode should not automatically be replaced by the fundamental Gaussian beam formula.

At equal argument and large integer order, the Bessel function of the first kind has phase $n(\sin t-t)=-n{t}^3/6+\cdots$ near the endpoint. The cubic stationary endpoint yields the displayed positive constant. It is also the zero-argument value of the Airy function in the uniform Bessel turning-point approximation.

For fixed $0<\alpha <\pi /2$, the Bessel function of the first kind has an interior nondegenerate stationary point in its cosine integral. The stationary phase method gives the displayed oscillatory leading term and an additive $O(n^{-3/2})$ remainder. The approximation loses uniformity near $\alpha =0$, where the stationary point reaches the endpoint and the curvature vanishes. For negative integer-order arguments, parity supplies a factor $(-1{)}^n$.

Diffraction-limited resolution is a scale for distinguishing details with a specified ideal diffraction response and a specified criterion. For two equally bright, mutually incoherent optical point sources imaged through a uniformly transmitting, unobscured circular optical pupil, the Rayleigh criterion places the centre of one Airy disk at the first zero of the other. At small angles in air this separation isHere $D$ is pupil diameter, $\lambda$ is wavelength in air, and $j_{1,1}$ is the first positive zero of the order-one Bessel function of the first kind. This is the radius to the first zero, not the diameter of the central spot. Fraunhofer diffraction gives the radial peak-normalized light intensityThe removable value at zero follows from $J_1(q)\sim q/2$. The pupil Fourier transform reduces to a radial integral by circular symmetry, which produces the Bessel function of the first kind; its first zero supplies the factor $1.22$.^[ref]

In the usual scalar circular-pupil approximation for a fluorescence microscope, the corresponding object-plane Rayleigh criterion isHere ${\lambda }_0$ is the vacuum emission wavelength, and the object-side numerical aperture already includes the medium's refractive index. Do not divide ${\lambda }_0$ by $n$ again. For small cone angles in air, $\mathrm{NA}\simeq D/(2f)$ and image-plane radius $f{\theta }_{\mathrm{R}}$ agrees with this form. The assumptions include a single wavelength, ideal optical focus, and mutually incoherent emitters. High numerical aperture can require a vector description with polarization of an electromagnetic wave and emitter orientation; confocal microscopy has its own excitation-and-detection response. The same $0.61$ formula therefore does not describe every microscope or illumination method.^[ref]

In this original model, ${\lambda }_0=550\mathrm{nm}$ and numerical apertures $0.25$ and $0.50$ give first-zero radii $1.342\mu \mathrm{m}$ and $0.671\mu \mathrm{m}$. Each point spread function is divided by its own peak to compare shape. The pupil axes are transverse wave vector components multiplied by ${\lambda }_0/(2\pi )$, so pupil radius equals numerical aperture. A wider pupil accepts a wider range of transverse wave vectors, producing a narrower image. Increasing magnification alone enlarges the same spot and does not reduce its object-plane radius.

The Rayleigh criterion is a useful visual convention, not an absolute information barrier. For equal incoherent sources with total detected radiant flux $F$ and separation vector $\mathbf{s}$, the ideal image isBelow the Rayleigh criterion this still differs from $Fp(\mathbf{x})$; for small $\mathbf{s}$, symmetry cancels the linear change and the first change is quadratic in separation. Estimating two unknown source positions is then increasingly sensitive to noise, background, relative brightness and uncertainty in the point spread function. Localizing one isolated source and recovering an arbitrary extended object are different inference tasks. Measurements that separate optical modes can also extract separation information differently from direct image light intensity measurement.^[ref]

There remains a real limit for arbitrary images in this model: the optical transfer function is zero above the transmitted spatial frequency band. Its cutoff for scalar incoherent imaging is ${\nu }_c=2\mathrm{NA}/{\lambda }_0$ in object-plane cycles per unit length. The Nyquist–Shannon sampling theorem therefore suggests object-referred detector pixel pitch $p_{\mathrm{pix}}\le {\lambda }_0/(4\mathrm{NA})$, with finer sampling providing margin. Finite detector pixels also integrate over their area, suppressing contrast; coarse sampling can cause aliasing. Finer sampling cannot alone recover spatial frequencies removed by the optics, and Fourier deconvolution amplifies noise where the optical transfer function is small. Photon shot noise, read noise, optical aberrations and defocus can all worsen practical discrimination. A useful resolution claim therefore reports the object model, criterion, wavelength, numerical aperture, sampling and signal-to-noise ratio in photon counting.

With the Fourier transform convention $\hat{f}(\omega )={\int }_{\mathbb{R}}f(t)e^{-i\omega t}dt$, the regular order-zero Bessel function of the first kind has the tempered distribution transformThe density is integrable despite its endpoint singularities. Fourier inversion gives the equivalent integral representationThe second formula follows by $\omega =\cos \theta$ and cancellation of the odd imaginary part. It is the standard real integral representation of $J_0$.^[ref]

Indeed, transforming the order-zero Bessel differential equation $t{y}^{\prime \prime }+y^{\prime }+ty=0$ gives $(1-{\omega }^2)Y^{\prime }-\omega Y=0$. Inside $(-1,1)$, its solutions are $C(1-{\omega }^2{)}^{-1/2}$. The compactly supported choice also solves the equation as a distribution: rewrite it as $[(1-{\omega }^2)Y{]}^{\prime }+\omega Y=0$, and observe that $(1-{\omega }^2)Y=C\sqrt{1-{\omega }^2}{1}_{(-1,1)}$ is continuous and zero at the two endpoints. Its distributional derivative has no endpoint Dirac delta terms. The inverse is regular at $t=0$, and $y(0)=1$ fixes $C=2$.

The transformed equation alone also allows separate constants on $\omega <-1$ and $\omega >1$. Selecting the compactly supported density amounts to selecting the solution regular at the origin, rather than every possible singular solution of the Bessel differential equation.

For the Bessel differential equation $x^2{y}^{\prime \prime }+x{y}^{\prime }+(x^2-{\nu }^2)y=0$, substitute $y=x^s{\sum }_{n\ge 0}{a}_n{x}^n$. The indicial equation is $s^2-{\nu }^2=0$, and coefficients obey $[(n+s{)}^2-{\nu }^2]a_n=-a_{n-2}$. For $s=\nu \ge 0$, $a_1=0$ and all odd coefficients vanish. Iteration gives $a_{2m}=(-1{)}^m{a}_0\Gamma (\nu +1)/[4^{m}m!\Gamma (m+\nu +1)]$. Choosing $a_0=1/[2^{\nu }\Gamma (\nu +1)]$ gives the Bessel function of the first kind $J_{\nu }$.

For noninteger $\nu$, choosing the even-power series with $s=-\nu$ replaces $m+\nu$ by $m-\nu$ and yields $J_{-\nu }$ with its standard normalization. These solutions are independent. For integer $\nu \ge 0$, the first solution has coefficients $(-1{)}^m/[2^{2m+\nu }m!(m+\nu )!]$ and defines an entire function: the term ratio tends to zero at every finite complex argument. For noninteger order, the even analytic series converges everywhere, while the prefactor $x^{\pm \nu }$ requires a branch off the positive real axis.

The Hankel functions are the combinations $H_{\nu }^{(1)}=J_{\nu }+i{Y}_{\nu }$ and $H_{\nu }^{(2)}=J_{\nu }-i{Y}_{\nu }$ of the Bessel function of the first kind and Bessel function of the second kind. They form a standard independent pair for the Bessel differential equation.

For real $k,r>0$ and time dependence $e^{-i\omega t}$, their leading radial factors are proportional to $e^{ikr}$ and $e^{-ikr}$. Constant wave phase $kr-\omega t$ moves outward with increasing time, whereas $-kr-\omega t$ moves inward. Reversing the time-sign convention reverses these labels. This makes the combinations useful when imposing a Sommerfeld radiation condition, rather than regularity at the origin.^[ref]

For complex arguments, the Hankel asymptotic sectors and the chosen branch of a multivalued function are part of the approximation. A formula continued around zero must not be treated as one global exponential on every argument sheet.

Invert the two-dimensional Fourier transform in the long-path projection of a magnetic correlation, then perform the angular integration in polar coordinates. The angular integral is $2\pi J_0(kr)$, a Bessel function of the first kind, giving the displayed Hankel transform with the stated Fourier-Hankel normalization for an isotropic spectrum. The result recovers the scalar tensor coefficient $H$, whose shell-integrated magnetic energy spectrum has an additional factor $k^2$.

For real $p>0$, the Laplace transform $G$ of the order-zero Bessel function of the first kind satisfies $(1+p^2)G^{\prime }+pG=0$. This follows by transforming $x{J}_0^{\prime \prime }+J_0^{\prime }+x{J}_0=0$, using $\mathcal{L}[xg]=-{\partial }_p\mathcal{L}[g]$ and $J_0(0)=1$, $J_0^{\prime }(0)=0$. Thus $G=C/\sqrt{1+p^2}$; the initial value theorem for Laplace transforms gives $pG(p)\to 1$ and fixes $C=1$.

Write $\eta =H{e}^{\sigma t}{J}_0(\alpha r)$ using the Bessel function of the first kind. The wall boundary conditions at $z=\pm h$ selectBecause $J_0^{\prime }(s)=-J_1(s)$, the sidewall no-penetration boundary condition requireswhere $j_{1,n}$ are the positive zeros of $J_1$. These radial eigenfunctions automatically preserve layer volume becauseThe constant radial mode $\alpha =0$ cannot change the interface height in a closed tank: the Laplace equation and the impermeable lids then give zero vertical velocity, and a uniform displacement would change each layer's volume.

The kinematic boundary condition fixesThus the interface values are ${\varphi }_1(0)=\sigma \eta \coth (\alpha h)/\alpha$ and ${\varphi }_2(0)=-\sigma \eta \coth (\alpha h)/\alpha$. Substitution in the normal stress condition gives the finite-depth Rayleigh-Taylor dispersion relation:soHere $\mathcal{A}$ is the Atwood number. Every allowed mode has a positive growth rate and is subject to Rayleigh-Taylor instability. Since $d[\alpha \tanh (\alpha h)]/d\alpha =\tanh (\alpha h)+\alpha h{\mathrm{sech}}^2(\alpha h)>0$, shorter radial scales grow faster. The limits areThere is no fastest finite wavenumber in the ideal sharp-interface model without surface tension, viscosity or diffusion. The zero-initial-velocity condition uses both exponential branches: an initial mode $H{J}_0({\alpha }_{n}r)$ evolves as $H\cosh ({\sigma }_{n,+}t)J_0({\alpha }_{n}r)$. If the mass density ordering were reversed the same calculation would give oscillatory stable modes; equal mass densities give neutral modes.

The contrast is weak: $(n_1-n_2)/n_2=0.01$. Use a scalar weak-guidance Helmholtz equation for a transverse field component, with longitudinal wave number $b$ and vacuum wave number $k_0=2\pi /{\lambda }_0$. Separating variables as $\psi =F_l(r)e^{il\varphi }{e}^{ibz}$ givesA bound mode requires $n_2{k}_0<b<n_1{k}_0$. DefineThe step-index fiber normalized frequency is $V$; $u$ is a dimensionless core transverse wave number. Regularity at $r=0$ selects the Bessel function of the first kind, while decay at infinity selects the Modified Bessel function of the second kind:In this scalar approximation, continuity of the field and its radial derivative givesAt cutoff the evanescent wave extends arbitrarily far into the cladding, so $b\downarrow n_2{k}_0$, $w\downarrow 0$ and $u\to V$. For the first azimuthally varying family $l=1$, the recurrence relations reduce the last equality toAs $w\to 0$, $K_1(w)\sim 1/w$ and $K_0(w)\sim -\ln w$, making the right side tend to zero. Thus the onset is at the first positive zero of $J_0$:The scalar $l=0$ fundamental family has no positive cutoff. The first additional spatial family is the $l=1$ family just found; the next radial $l=0$ and first $l=2$ families require the larger zero $j_{1,1}\simeq 3.832$. This establishes the first multimode boundary rather than merely finding one possible cutoff. The two polarizations of the fundamental family are not counted as two spatial modes. Exact vector interface conditions split some higher-mode degeneracies at finite contrast; the displayed condition is the weak-guidance result intended by the supplied scalar equations.

The separated disk modes. Away from the plane of a razor-thin disk, there is no volume source, so the Newtonian gravitational potential obeys the Laplace equation in cylindrical coordinates. Separation of variables givesHere $k>0$ is a radial wavenumber, with units of inverse length. Regularity on the axis selects the Bessel function of the first kind $J_0(kR)$ rather than the singular second solution. Decay away from the disk, continuity across it and reflection symmetry select $Z=e^{-k|z|}$. These conditions apply to each nonzero-wavenumber mode; a uniform sheet is a separate zero-mode problem.

Integrate the Poisson equation for Newtonian gravity across the disk to obtain the jump conditionFor the unit-amplitude separated mode, the jump is $-2k{J}_0(kR)$. ThusAn individual basis density is signed; positive physical surface densities are built by superposition.

An arbitrary axisymmetric surface density. Define its order-zero Hankel transform byMatching the coefficients of the mode densities gives the Hankel representation of thin-disk gravity,The boundary conditions exclude added external fields. For sufficiently decaying surface density, this is the isolated-disk potential; for infinite scale-free disks, only potential differences or regularized forces need be finite.

The Mestel disk. Put $C={\Sigma }_0{R}_0$. The Bessel function integral ${\int }_0^{\infty }{J}_0(k{R}^{\prime })d{R}^{\prime }=1/k$, understood with a decaying regulator, gives $\hat{\Sigma }=C/k$. Differentiate the potential before removing the regulator and use $J_0^{\prime }=-J_1$:The last integral is $1/R$, since ${\int }_0^{\infty }{J}_1(u)du=[-J_0(u){]}_0^{\infty }=1$. Therefore the circular speed and enclosed mass areConsequently $v_c^2=GM(R)/R$. This equality follows from the specific Mestel disk calculation; the spherical shell theorem cannot be applied to general disks. Its total mass and absolute potential referenced to infinity diverge. A convenient finite midplane reference is $\Phi (R,0)=v_0^2\log (R/R_0)$.

The rotating distribution function. Use specific energy $E=(v_R^2+v_{\varphi }^2)/2+\Phi$ and $L_z=R{v}_{\varphi }$. For the printed positive exponential, take $\epsilon =-E$ and $s:=\sigma >0$; $s$ has units of squared velocity. If instead $\epsilon$ means ordinary total energy and $\sigma >0$, the velocity integral diverges, so that literal convention cannot define the requested galactic distribution function. Equivalently one can retain ordinary energy and choose the printed $\sigma$ negative.

With the binding-energy convention, write the amplitude as $A$ to distinguish it from the function:Integration over both radial velocities and positive azimuthal velocities givesThe half-line integral follows by setting $u=v_{\varphi }^2/(2s)$ in the defining Gamma function integral and requires $q>-1$. Matching its radial power and normalization to the Mestel disk yields the rotating Mestel disk distribution function parametersSince $v_0^2,s>0$, the required integrability condition is automatic. The amplitude uses the stated potential reference; changing the additive potential constant changes $A$ without changing the physical distribution.

For an axisymmetric potential, insert $\Psi (r,\varphi )=R(r)\Phi (\varphi )$ into the time-independent Schrödinger equation and divide by $R\Phi$. The polar Laplacian separates the angular equation ${\Phi }^{\prime \prime }+m^2\Phi =0$ fromSingle-valuedness under $\varphi \mapsto \varphi +2\pi$ gives $m\in \mathbb{Z}$, with an angular eigenbasis $\boxed{\Phi =e^{im\varphi }}$. Linear combinations of the degenerate $m,-m$ modes are also possible.

Inside the circular infinite quantum well, put $\rho =r\sqrt{2\mu E/{\hslash }^2}$. The equation becomesThe physical solution must be regular at the origin, behaving as $r^{|m|}$ for $m\ne 0$, or with finite value and zero radial derivative for $m=0$. At the infinite wall impose $R(r=a)=0$. The printed question's boundary phrase uses capital $R=a$; it means the radial position $r=a$, not a value imposed on the radial wavefunction. The singular Bessel function of the second kind is excluded; mere square-integrability of its logarithmic $m=0$ singularity would not put it in the ordinary finite-energy Hamiltonian domain.

For $m=0$, multiply the equation by ${\rho }^2$ and insert $R={\sum }_{k\ge 0}{A}_k{\rho }^k$. The coefficient of $\rho$ gives $A_1=0$, and subsequent coefficients giveAll odd coefficients vanish, $A_2=-A_0/4$, and $A_4=A_0/64$. The regular solution is a multiple of the Bessel function of the first kind $J_0$. Its requested fourth-order estimate isTherefore the first positive zero of this truncation is $\boxed{\rho \approx \sqrt{8}}$, givingThe ground state is in the $m=0$ sector because nonzero angular modes add nonnegative centrifugal kinetic energy to the radial energy form. This quartic truncation of the Bessel function J0 is only a rough local-series estimate, not an exact eigenvalue or a certified variational bound; its double zero warns that the truncation does not reproduce the simple zero of the full function.

Invert the two-dimensional Fourier transform from part (c):In polar coordinates, the angular integral is $2\pi J_0(ks)$, where $J_0$ is a Bessel function of the first kind. The Fourier-Hankel normalization for an isotropic spectrum therefore yieldsThis Hankel inversion of a rotation-measure correlation recovers the scalar spectral tensor coefficient. To obtain the shell-integrated magnetic energy spectrum, take the trace of the three-dimensional tensor:With magnetic energy density $B^2/(8\pi )$ in Gaussian units, the one-dimensional convention ${\int }_0^{\infty }{E}_B(k)dk=\mathbb{E}|\mathbf{B}{|}^2/(8\pi )$ gives $\boxed{E_B(k)=k^{2}H(k)/(8{\pi }^3)}$. A spectrum normalized instead to $\mathbb{E}|\mathbf{B}{|}^2/2$ has coefficient $k^{2}H(k)/(2{\pi }^2)$; the physical normalization must be stated.