N. A. Lopez
Abstract:The propagation of electromagnetic waves in a linearly-varying index of refraction is a fundamental problem in wave physics, being relevant in fusion science for describing certain wave-based heating and diagnostic schemes. Here, an exact solution is obtained for a given incoming wavefield specified on the boundary transverse to the direction of inhomogeneity by performing a spectral, rather than asymptotic, matching. Two case studies are then presented: a Gaussian beam at oblique incidence and a speckled wavefield at normal incidence. For the Gaussian beam, it is shown that when the waist $W$ is sufficiently large, oblique incidence manifests simply as rigid translation and focal shift of the corresponding diffraction pattern at normal incidence. The destruction of the hyperbolic umbilic caustic (corresponding to a critically focused beam) as $W$ is reduced is then demonstrated. The caustic disappears once $W \lesssim \delta_a \sqrt{L}$ ($L$ being the medium lengthscale normalized by the Airy skin depth $\delta_a$), at which point the wave behavior is increasingly described by Airy functions, but experiences less focusing as a result. To maximize the intensity of a launched Gaussian beam at a turning point, one should therefore minimize the imaginary part of the launched complex beam parameter while having the real part satisfy critical focusing. For the speckled wavefield, it is shown that the transverse speckle pattern only couples to the Airy longitudinal pattern when the coupling parameter $\eta = \sqrt{L}/f_\#$ is large, with $f_\#$ being the f-number of the launching aperture. When $\eta \ll 1$, a reduced description of the total wavefield can be obtained by simply multiplying the incoming speckle pattern with the Airy swelling.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the classical problem of electromagnetic wave propagation in a linearly varying refractive index medium, namely the so - called "linear layer problem". Specifically, the paper focuses on how the wave field is distributed throughout the domain when the incident field on the given boundary is specified. Previous studies usually could only obtain the asymptotic solution to this problem, while this paper obtains the exact solution by using the spectral - matching scheme. This exact solution (Equation 3.11) is expressed as an integral, which contains the product of the Fourier transform of the predetermined incident boundary field and a kernel function. The kernel function contains not only the familiar Airy function but also the related Skorokhod function.
### Main research contents and results:
1. **Derivation of the exact solution**:
- The paper derives an exact solution by the spectral - matching method. This solution is applicable to any setting of the incident wave field on the boundary and is not affected by the relative magnitude of the incident wavelength and the medium length scale.
- The form of the solution is an integral, and the integral kernel contains the Airy function and the Skorokhod function, ensuring the non - singularity and universal validity of the solution.
2. **Special case analysis of Gaussian beams**:
- When the incident wave field is a Gaussian beam, the paper gives the corresponding exact solution (Equation 4.5).
- The analysis shows that when the beam is wide enough (i.e., Gaussian - focused plane wave), the oblique incidence angle will lead to a simple translation of the focal length and a rigid displacement of the diffraction pattern.
- Explicit expressions (Equations 4.8 and 4.9) are provided, which correct the similar expressions in Lopez et al. (2023) and remove the Goos - Hanchen displacement that does not occur when there is no finite beam waist.
- It is further shown that as the beam waist decreases, the hyperbolic umbilic singularity corresponding to critical focusing will gradually soften and eventually disappear, resulting in a decrease in the maximum intensity of the Gaussian beam at the turning point. This suggests that in applications that depend on the wave field intensity near the turning point, the imaginary part of the complex beam parameter of the emitted beam should be reduced as much as possible.
3. **Special case analysis of speckle plane waves**:
- The normally incident speckle plane wave field is studied.
- A coupling parameter \(\eta=\sqrt{L}/f\#\) is derived, where \(L\) is the medium length scale and \(f\#\) is the f - number of the emitting aperture.
- The results show that when \(\eta\ll1\), the transverse speckle pattern is decoupled from the longitudinal Airy pattern, and at this time, the complete solution can be obtained by simply multiplying the given incident speckle pattern by a single Airy profile.
- Conversely, when \(\eta\gg1\), the speckle pattern has a significant impact on the longitudinal profile of the total wave field, showing a complex interference pattern, and its longitudinal expansion is significantly different from the standard Airy expansion.
### Conclusion:
In this paper, the exact solution of the linear layer problem is successfully solved by the spectral - matching method, which is applicable to the setting of any incident wave field. For the special case analysis of Gaussian beam and speckle plane wave fields, in - depth physical insights are provided, especially regarding critical focusing and the influence of speckle patterns. These results are of great significance for wave heating and diagnostic applications in plasma physics and fusion research.