## Abstract

Since the development of laser light sources in the early 1960s, laser beams are everywhere. Laser beams are central in many industrial applications and are essential in ample scientific research fields. Prime scientific examples are optical trapping of ultracold atoms, optical levitation of particles, and laser-based detection of gravitational waves. Mathematically, laser beams are well described by Gaussian beam expressions. Rather well covered in the literature to date are basic expressions for scalar Gaussian beams. In the past, however, higher accuracy mathematics of scalar Gaussian beams and certainly high-accuracy mathematics of vectorial Gaussian beams were far less studied. The objective of the present review then is to summarize and *advance* the mathematics of vectorial Gaussian beams. When a weakly diverging Gaussian beam, approximated as a linearly polarized two-component plane wave, say (${E}_{x},{B}_{y}$), is tightly focused by a high-numerical-aperture lens, the wave is “depolarized.” Namely, the prelens (practically) missing electric field ${E}_{y},{E}_{z}$ components suddenly appear. This is similar for the prelens missing ${B}_{x},{B}_{z}$ components. In fact, for any divergence angle (${\theta}_{d}<1$), the ratio of maximum electric field amplitudes of a Gaussian beam ${E}_{x}\text{:}{E}_{z}\text{:}{E}_{y}$ is roughly $1\text{:}{\theta}_{d}^{2}\text{:}{\theta}_{d}^{4}$. It follows that if a research case involves a tightly focused laser beam, then the case analysis calls for the mathematics of *vectorial* Gaussian beams. Gaussian-beam-like distributions of the six electric–magnetic vector field components that nearly exactly satisfy Maxwell’s equations are presented. We show that the near-field distributions with and without evanescent waves are markedly different from each other. The here-presented nearly exact six electric–magnetic Gaussian-beam-like field components are symmetric, meaning that the cross-sectional amplitude distribution of ${E}_{x}(x,y)$ at any distance ($z$) is similar to the ${B}_{y}(x,y)$ distribution, ${E}_{y}(x,y)$ is similar to ${B}_{x}(x,y)$, and a 90° rotated ${E}_{z}(x,y)$ is similar to ${B}_{z}(x,y)$. Components’ symmetry was achieved by executing the steps of an outlined *symmetrization procedure*. Regardless of how tightly a Gaussian beam is focused, its divergence angle is limited. We show that the full-cone angle to full width at half-maximum intensity of the dominant *vector* field component does not exceed 60°. The highest accuracy field distributions to date of the less familiar higher-order Hermite–Gaussian *vector components* are also presented. Hermite–Gaussian $\mathbb{E}\text{-}\mathbb{B}$ vectors only approximately satisfy Maxwell’s equations. We have defined a *Maxwell’s-residual power* measure to *quantify* the approximation quality of different vector sets, and each set approximately (or exactly) satisfies Maxwell’s equations. Several vectorial “applications,” i.e., research fields that involve vector laser beams, are briefly discussed. The mathematics of vectorial Gaussian beams is particularly applicable to the analysis of the physical systems associated with such applications. Two user-friendly “Mathematica” programs, one for computing six *high-accuracy* vector components of a Hermite–Gaussian beam, and the other for computing the six practically *Maxwell’s-equations-satisfying components* of a focused laser beam, supplement this review.

© 2019 Optical Society of America

Full Article | PDF Article**OSA Recommended Articles**

Qiwen Zhan

Adv. Opt. Photon. **1**(1) 1-57 (2009)

Pavel Zemánek, Giorgio Volpe, Alexandr Jonáš, and Oto Brzobohatý

Adv. Opt. Photon. **11**(3) 577-678 (2019)

Ajanta Barh, Peter John Rodrigo, Lichun Meng, Christian Pedersen, and Peter Tidemand-Lichtenberg

Adv. Opt. Photon. **11**(4) 952-1019 (2019)