--- title: Far Field Gaussians (FFG) Overview url: https://www.emergentmind.com/topics/far-field-gaussians-ffg type: topic --- # Far Field Gaussians (FFG) Overview Searching arXiv for the cited papers and related usage of “Far Field Gaussians”. Far Field Gaussians (FFG) denotes a family of far-field representations in which the field, kernel, or long-range interaction is written in Gaussian or Gaussian-derived form and then propagated or evaluated analytically or spectrally. The available usage suggests that FFG is better understood as a modeling viewpoint than as a single standardized formalism. In optical beam theory, it appears in the far-field analysis of four-petal Gaussian beams (FPGBs), including hard-edged diffraction by a circular aperture, where Gaussian envelopes, Gaussian aperture expansions, vector angular spectra, and stationary-phase asymptotics yield closed-form far-field expressions [1003.5768][1004.4786]. In quasi-2D electrostatics and periodic Yukawa molecular dynamics, it refers to the long-range component represented as a sum of smooth Gaussians, enabling Fourier–Chebyshev or random-batch Fourier treatments with \(O(N\log N)\) and \(O(N)\) complexity, respectively [2412.04595][2501.10946]. A related time-domain wave-propagation literature addresses the same far-field evaluation problem through compressed sums of exponentials rather than Gaussian bases [1409.5893]. ## 1. Terminological scope and principal usages The literature represented here contains two direct FFG constructions and one closely related analogue. In the optical setting, the far field is obtained from Gaussian-type beams and, in the apertured case, from a hard-edged circular aperture expanded into a finite sum of complex Gaussian functions [1004.4786]. In the electrostatic and Yukawa settings, the far field is the smooth long-range part of the interaction kernel, represented explicitly as a sum of Gaussians and then evaluated in Fourier space or through hybrid spectral methods [2412.04595][2501.10946]. In time-domain wave propagation, the far-field recovery problem is solved by spherical-harmonic decomposition and compressed sums of exponentials; this is not a Gaussian-basis method, but it is directly concerned with efficient and stable far-field evaluation [1409.5893]. | Context | Object represented | Core machinery | |---|---|---| | FPGB optics | Far-field beam structure | Vector angular spectrum, stationary phase, Gaussian aperture expansion | | Quasi-2D electrostatics | Long-range Coulomb part | Sum-of-Gaussians, Fourier spectral solver, Fourier–Chebyshev proxy points | | 3D Yukawa MD | Far-field interaction | Sum-of-Gaussians, random batch approximation, adaptive importance sampling | | Time-domain waves | Far-field signal recovery | Multipoles, teleportation kernels, compressed sum of exponentials | A common misconception is to treat FFG as synonymous with a single beam family or a single numerical algorithm. The cited literature does not support that identification. Instead, the recurring structure is the replacement of a difficult far-field object—diffracted beam, long-range kernel, or asymptotic waveform—by a compact representation built from analytically tractable basis functions. ## 2. Optical FFG for four-petal Gaussian beams For FPGBs, the incident beam at \(z=0\) is linearly polarized along \(x\) and propagates along \(+z\), with transverse field \[ E_x(x,y,0)=G_n\left(\frac{xy}{w_0^2}\right)^{2n}\exp\!\left(-\frac{x^2+y^2}{w_0^2}\right),\qquad E_y(x,y,0)=0, \] where \(n\) is the beam order, \(G_n\) is an amplitude constant, and \(w_0\) is the Gaussian waist radius [1003.5768][1004.4786]. The near-field four-petal structure arises from the factor \((xy)^{2n}\), and the separation of the diagonal petals at \(z=0\) is \[ d=2\sqrt{2n}\,w_0 \] in the apertured treatment [1004.4786]. The optical papers further state that the FPGB can be regarded as a superposition of two-dimensional Hermite–Gaussian modes; in one formulation, the beam with order \(n\) is described as a superposition of \(n^2\) two dimensional Hermite-Gaussian modes [1003.5768]. This modal interpretation is central to the eventual multi-lobe far-field pattern. The hard-edged aperture is a circular aperture of radius \(R\) at \(z=0\), \[ \mathrm{circ}\!\left(\frac{\rho}{R}\right)= \begin{cases} 1,& \rho\le R,\\ 0,& \rho>R, \end{cases} \qquad \rho=\sqrt{x^2+y^2}, \] and the field immediately behind the aperture is the product of the incident field with this aperture function [1004.4786]. The key analytical step is that the aperture function is expanded into a finite sum of complex Gaussian functions, \[ \mathrm{circ}\!\left(\frac{\rho}{R}\right)=\sum_{l=1}^{N}A_l\exp\!\left(-B_l\frac{\rho^2}{R^2}\right), \] with numerical complex coefficients \(A_l\), \(B_l\), and \(N\) taken from Wen & Breazeale’s table [1004.4786]. This converts hard-edge truncation into a weighted sum of Gaussian windows. Two dimensionless parameters organize the analysis. The truncation parameter is \[ \beta=\frac{R}{w_0}, \] which measures aperture size relative to beam waist, and the \(f\)-parameter is \[ f=\frac{1}{k w_0}=\frac{\lambda}{2\pi w_0}, \] which controls nonparaxiality [1004.4786]. The power transmissivity through the aperture is given explicitly by \[ T_n=1-\frac{\Gamma(4n+1,2\beta^2)}{4n\,\Gamma(4n)}, \] with the fundamental Gaussian case \[ T_0=1-e^{-2\beta^2}, \] so the transmitted power is governed directly by \(\beta\) [1004.4786]. ## 3. Vectorial far-field structure in the optical formulation The optical FFG derivation uses the vector angular spectrum method. The field is represented as a superposition of plane waves in normalized spatial frequencies \(p,q\), with \[ \gamma= \begin{cases} \sqrt{1-p^2-q^2},& p^2+q^2\le 1,\\ i\sqrt{p^2+q^2-1},& p^2+q^2>1, \end{cases} \] and the angular spectrum just behind the aperture is obtained by Fourier transforming the truncated field [1004.4786]. In the un-apertured case, the same formalism leads directly to closed-form expressions for \(A_x(p,q)\) involving the confluent hypergeometric function \({}_1F_1\), with \(A_y(p,q)=0\) [1003.5768]. The field is then decomposed into transverse electric and transverse magnetic parts, \[ \vec{E}(\vec r)=\vec{E}_{TE}(\vec r)+\vec{E}_{TM}(\vec r),\qquad \vec{H}(\vec r)=\vec{H}_{TE}(\vec r)+\vec{H}_{TM}(\vec r), \] and the stationary phase method is applied to obtain explicit far-field expressions [1003.5768][1004.4786]. The resulting far-field formulas contain a spherical phase factor \(\exp(ikr)\), geometrical prefactors in \(x,y,z,r,\rho\), and angular dependence encoded by confluent hypergeometric functions evaluated at arguments proportional to \(-x^2/(4f^2r^2)\) and \(-y^2/(4f^2r^2)\) in the apertured case [1004.4786]. This is the precise sense in which the far field is Gaussian-like but not merely Gaussian: the Gaussian envelope is modulated by polynomial or hypergeometric structure inherited from the FPGB order. The TE and TM electric fields are orthogonal in the far field, \[ \vec E_{TE}\cdot \vec E_{TM}=0, \] and the corresponding magnetic fields are also orthogonal, \[ \vec H_{TE}\cdot \vec H_{TM}=0, \] so the far-field energy flux is the sum of separate TE and TM contributions without mixed TE–TM interference terms [1003.5768][1004.4786]. The TE electric field is purely transverse, while the TM electric field contains a longitudinal component proportional to \(-\rho^2\hat e_z\), which is one of the signatures of nonparaxial vectorial structure [1003.5768]. This optical formulation also clarifies a second misconception: the far field of a four-petal Gaussian beam is not, in general, a four-petal Gaussian. The cited papers show that the initial four-petal pattern evolves into multiple small petals in the far field, and that the number of petals increases with beam order \(n\) [1003.5768][1004.4786]. ## 4. Energy flux, nonparaxiality, and limiting cases in FPGB far fields The time-averaged far-field Poynting flux is written as \[ \langle S_z\rangle=\langle S_z\rangle_{TE}+\langle S_z\rangle_{TM}, \] with separate closed forms for the TE term, the TM term, and the total beam [1003.5768][1004.4786]. In the apertured formulation, the total flux depends on a function \(S_n(x,y,z,f,\beta)\) built from the Gaussian-expansion coefficients \(A_l,B_l\) and the two confluent hypergeometric factors in \(x\) and \(y\) [1004.4786]. The flux therefore carries explicit dependence on beam order \(n\), aperture parameter \(\beta\), and nonparaxial parameter \(f\). The numerical behavior reported for the apertured beam is specific. For \(n=1\), \(w_0=10\lambda\), and \(z=2000\lambda\), decreasing \(\beta\) from \(\infty\) to \(\sqrt{2}\) to \(1\) increases the central spot full-width-at-half-maximum from approximately \(41.6\lambda\) to \(61.8\lambda\) to \(83.8\lambda\), and the ratio of first side-lobe peak to central peak increases from approximately \(0.19\) to \(0.32\) to \(0.34\) [1004.4786]. This is the hard-edge diffraction effect: smaller apertures broaden the central lobe and elevate relative side-lobe levels. The paraxial and nonparaxial regimes are separated by the magnitude of \(f\). In the paraxial limit \(f\ll 1\), the far-field radial distance can be approximated by \(r\approx z\), and the total far-field energy flux simplifies to a form with prefactor \(1/r^2\), with asymmetry carried only by \(S_n\) [1004.4786]. Numerical results show that for moderate \(f\le 1\) and \(\beta\sim 2\), paraxial and nonparaxial results are nearly identical, whereas for \(f=1.5\) nonparaxial corrections become significant, the central lobe becomes elliptical, and the paraxial prediction fails particularly in the small-angle region [1004.4786]. The same study states explicitly that the \(f\)-parameter is more critical for nonparaxiality than \(\beta\). The principal limiting case is \(\beta\to\infty\), corresponding to the un-apertured beam. In that limit the aperture expansion is no longer needed, the factor \(\beta^2/(\beta^2+B_l)\) tends to \(1\), and the apertured formulas degenerate into the un-apertured case [1004.4786]. This degeneracy is part of the analytical consistency of the hard-edge model. ## 5. Sum-of-Gaussians far fields in quasi-2D electrostatics In quasi-2D electrostatics, FFG has a different object but the same structural idea: the long-range part of the kernel is represented by smooth Gaussians. The problem consists of point charges in a rectangular box with periodic boundary conditions in \(x,y\) and free boundary conditions in \(z\), and the quasi-2D Coulomb lattice sum is built from the kernel \(K(\mathbf r)=1/|\mathbf r|\) with periodic replication in the periodic directions [2412.04595]. Charge neutrality is imposed. The central step is a sum-of-Gaussians approximation to the Laplace kernel, \[ \mathcal F_b^\sigma(r)=\sum_{\ell=0}^{M}w_\ell e^{-r^2/s_\ell^2}, \] with \[ w_\ell=(\pi/2)^{-1/2}b^{-\ell}\sigma^{-1}\log b,\qquad s_\ell=\sqrt{2}\,b^\ell\sigma, \] and a near/far decomposition \[ \frac{1}{r}=\mathcal N_b^\sigma(r)+\mathcal F_b^\sigma(r), \] where the near field is singular but compactly supported and the far field is smooth [2412.04595]. The paper further refines the far part into near-field, mid-range, and long-range components. The long-range component consists of Gaussians with large variance and is precisely what the source text identifies as Far Field Gaussians. The quasi-2D solver uses the separability of Gaussians. Mid-range Gaussians are handled by a procedure similar to nonuniform fast Fourier transforms in three dimensions, while the long-range component is treated with polynomial interpolation/anterpolation in the free dimension and a Fourier spectral solver in the other two dimensions on proxy points [2412.04595]. The long-range solver is formulated in a Fourier–Chebyshev basis: Fourier in the periodic coordinates and Chebyshev in \(z\). A theorem in the paper gives the Chebyshev approximation error bound \[ \left|e^{-z^2/s_\ell^2}-\sum_{n=0}^{P-1}{}' a_n T_n(2z/L_z)\right| \le \frac{1}{\sqrt{P!}(2\sqrt{2}\eta)^P}, \qquad s_\ell\ge \eta L_z, \] so few Chebyshev terms suffice for large-variance Gaussians [2412.04595]. The method’s main algorithmic consequence is that zero padding in the free direction and upsampling in NUFFT-like steps are not required. The paper attributes this to the smoothness and rapid Fourier decay of Gaussians, in contrast with fast Ewald summation and truncated-kernel approaches [2412.04595]. It reports rigorous error analysis and overall \(O(N\log N)\) complexity with a small prefactor. The broader context given in the paper connects this construction to the u-series work of Predescu et al., to slit-geometry work by Maxian et al., to periodic FMM methods of Yan–Shelley and Pei–Askham–Greengard–Jiang, and to standard NUFFT formulations such as Nestler–Pippig–Potts. ## 6. Random-batch FFG for periodic Yukawa systems For 3D Yukawa systems with periodic boundary conditions, the kernel is \[ Y(r)=\frac{e^{-r/\lambda}}{r}, \] and the FFG construction is again a sum-of-Gaussians decomposition, now adapted to the screened interaction [2501.10946]. The paper starts from an integral representation \[ Y(r)=\int_{-\infty}^{\infty} f(t,r)\,dt, \] then discretizes the integral by a trapezoidal rule to obtain \[ \mathcal F_h^{t_0}(r)=h\sum_{m=-M_1}^{M_2} f(t_m,r), \] which is the smooth far-field Gaussian part [2501.10946]. The near-field residual is \[ \mathcal N_h^{t_0}(r)= \begin{cases} Y(r)-\mathcal F_h^{t_0}(r),& r