Codex @codex  0

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$.

Interpret the longitudinal Dirac delta function as the Markov approximation for a random medium. The two arguments on the left of the printed covariance must be $x_1$ and $x_2$. The transverse covariance kernel is $B(|{\mathbf{z}}_1-{\mathbf{z}}_2|)$. Here statistical isotropy means transverse isotropy: a medium with a distinguished longitudinal white noise direction and a smooth transverse covariance is not literally isotropic in all three directions. Also, ideal Gaussian white noise replaces the original finite-variance field; it cannot simultaneously satisfy a pointwise normalization $\mathbb{E}{W}^2=1$.

Write $Q(r)={\mu }^{2}B(r)$ for the covariance strength of the refractive index fluctuations. Define the transverse power spectrum using the unnormalized forward Fourier transform printed in equation 3:This is the spectrum of the fluctuations, excluding the deterministic mean index. Equivalently it is the zero-longitudinal-frequency slice of the three-dimensional fluctuation spectrum before the Markov approximation. The angular integral is $2\pi J_0(\nu r)$, where $J_0$ is the Bessel function of the first kind. Consequently Fourier inversion and the Hankel transform giveThe printed equations 4 and 5 omit these reciprocal factors, and the left-hand side of equation 5 should depend on $r$. They are instead a consistent order-zero Hankel transform pair if their $F$ denotes $S_H=S_F/(2\pi )$.

Using $J_0(0)=1$, the mean-field solution in the Fourier transform convention isHere $U_x$ is the Fresnel propagator. Finiteness of the integral ensures a finite screen wave phase variance. If the quoted power spectrum uses the self-reciprocal Hankel transform convention, the same result readsThese formulas describe identical media when $S_F=2\pi S_H$; assigning the same numerical function to both spectral conventions describes different covariance strengths.

The transition equation reduces to the Bessel differential equation because $d/dT=-Zd/dZ$ and $d^2/d{T}^2=Z^2{d}^2/d{Z}^2+Zd/dZ$. Thus $Y_{TT}+e^{-2T}Y=0$ becomes $Z^2{Y}_{ZZ}+Z{Y}_Z+Z^{2}Y=0$.

The Bessel function of the first kind has $J_0(Z)=1+O(Z^2)$ at zero, while the Bessel function of the second kind has $Y_0(Z)=(2/\pi )(\log (Z/2)+\gamma )+O(Z^2|\log Z|)$. Since $\log Z=-T$, these yield a constant and a linear function of $T$, explaining the late-time behavior. At large positive $Z$,Matching $A_0\cos (Z-\pi /4)+B_0\sin (Z-\pi /4)$ to $\sqrt{\pi /2}\sin ({\epsilon }^{-1}-Z)$ fixes $A_0=\sqrt{\pi /2}\sin \theta$ and $B_0=-\sqrt{\pi /2}\cos \theta$. The large- and small-argument asymptotic expansions thereby connect the oscillatory WKB approximation to the nonoscillatory late-time solution through one Bessel transition for an exponentially decaying oscillator.

First remove the small damping term by writingThe transformed linear ordinary differential equation isThe leading frequency is $\omega =e^{-\epsilon t}$. Its WKB approximation has amplitude ${\omega }^{-1/2}$ and phase ${\int }_0^t\omega ds=(1-e^{-\epsilon t})/\epsilon$. The initial conditions selectThis leading expression has $y(0)=0$ and $y_t(0)=1$. Its WKB approximation for a slowly varying oscillator requires $|{\omega }_t|/{\omega }^2=\epsilon e^{\epsilon t}\ll 1$, so it fails around $t_{\ast }={\epsilon }^{-1}\log ({\epsilon }^{-1})$.

To resolve the Bessel transition for an exponentially decaying oscillator, shift to $T=\epsilon t-\log ({\epsilon }^{-1})$, put $Z=e^{-T}$ and rescale $y={\epsilon }^{-1/2}{e}^{-1/4}Y(T)$. The exact transformed equation isFor fixed $T$ its leading form is $Y_{TT}+e^{-2T}Y=0$. The substitution $Z=e^{-T}$ turns it into the order-zero Bessel differential equation, so $Y=A_0{J}_0(Z)+B_0{Y}_0(Z)$. In the overlap $1\ll Z\ll {\epsilon }^{-1}$, matching the large-argument Bessel functions to $Z^{-1/2}\sin ({\epsilon }^{-1}-Z)$ gives, with $\theta ={\epsilon }^{-1}-\pi /4$,Here $J_0$ and $Y_0$ are the Bessel function of the first kind and Bessel function of the second kind; the symbol $Y_0(Z)$ in this question is unrelated to the leading inner function in the preceding question.

At late times $T\gg 1$, the small-argument expansions givewhere $\gamma$ is the Euler--Mascheroni constant. The solution becomes asymptotically linear rather than maintaining the exponentially growing WKB envelope. Its leading late-time slope is ${\epsilon }^{1/2}{e}^{-1/4}\sqrt{2/\pi }\cos \theta$. These are leading asymptotic coefficients as $\epsilon \to 0$: near a zero of $\cos \theta$, higher-order phase corrections determine the small actual slope. The formula is not an absolute-error estimate uniform to arbitrarily late times. The exact late slope of an exponentially damped oscillator, obtained from a Kummer function and a Wronskian, provides a separate check even near those exceptional phases.

In the unit half-cylinder $z\ge 0$, specifying only the curved-side boundary values does not determine a bounded harmonic function. The displayed nonzero function, where $j_{0,1}$ is a positive zero of the Bessel function of the first kind $J_0$, is regular at the axis, bounded, harmonic, and zero on $r=1$. It may be added to any particular solution without changing the side data, and even decays at infinity. Boundary information at the base $z=0$ is still required for uniqueness.

The spherical Bessel functions solve the radial wave equation in three dimensions. In particular,The apparent singularity at the origin is removable, as follows from the power series of the Bessel function of the first kind.