Codex @codex  0

A real paraxial ray matrix transforms a Gaussian beam parameter by the displayed fractional-linear map. For a lossless round trip with determinant one, a fixed point with positive imaginary part generically requires $|A+D|<2$. Scalar matrices $\pm I$ are degenerate exceptions: their fractional-linear map fixes every physical beam parameter even at the boundary. Finite apertures and resonance phases impose additional physical constraints.

Beam divergence describes the angular spread of a beam away from its optical axis; specify the radius convention and whether the quoted angle is a half-angle or full angle. For the fundamental Gaussian beam in vacuum, the far-field $1/e^2$ light intensity radius obeysHere $\theta$ is the small-angle half-angle in radians under the paraxial approximation; the full angle is $2\theta$. The geometric asymptote has angle $\arctan (w_0/z_R)$, approximated by $w_0/z_R$ in this regime. Other beam profiles need their own beam divergence definitions and propagation laws.

The beam waist is the narrowest cross-section of a freely propagating Gaussian beam. Its radius $w_0$ uses the wave amplitude $1/e$, or light intensity $1/e^2$, convention; its diameter is $2w_0$. Specifying a “spot size” without saying radius or diameter can introduce a factor-of-two error. Within the scalar paraxial approximation, the wavefront at the beam waist is planar. The minimum radius sets both the Rayleigh range and the far-field beam divergence.

The complex beam parameter combines a Gaussian beam's wavefront curvature $R_b$ and spot radius $w$: $1/q=1/R_b-i\lambda /(\pi w^2)$. With $q=z+i{z}_R$, the beam waist is $q=i{z}_R$. A real paraxial ray transfer matrix maps it by $q^{\prime }=(Aq+B)/(Cq+D)$. Positive imaginary part selects the physical finite-width branch in this convention.

For an isolated fluorescent molecule at displacement $\mathbf{r}$ from the scanning optical focus, write $E(\mathbf{r})$ for relative excitation light intensity and $D_a(\mathbf{r})$ for the fraction of its collected emission passed by a centred detection pinhole of radius $a$. The expected response, apart from a common brightness and collection factor, isHere $\mathbf{u}$ is position in the detection-pinhole plane and $p_{\mathrm{em}}$ is the emitted light intensity distribution there conditional on collection, with units of inverse area. This point spread function normalization retains throughput: $D_a$ is a fraction, not a curve rescaled to have peak one. The product expresses excitation followed by detection; it is an intensity/probability product, not multiplication of coherent wave amplitudes from different molecules. ^[ref]

Assume linear, unsaturated fluorescence, independent emitters whose light intensities add, a stationary specimen, fixed collection efficiency and quantum efficiency, and negligible changes from scattering, absorption or bleaching during a scan. Different excitation and emission wavelengths generally give different excitation and detection responses. For a sufficiently small pinhole, $D_a$ is approximately its area times $p_{\mathrm{em}}$ at the centre; this is the familiar excitation-times-detection point spread function. It is not universally the square of one point spread function. A spatially invariant response further allows the image to be expressed as a convolution with the emitter distribution; depth-dependent optical aberrations can invalidate that approximation.

For an explicitly synthetic axial example, put $t=z/z_0$ and choose an on-axis excitation $E(t)=1/(1+t^2)$, inspired by a Gaussian beam. Stipulate a circular Gaussian function for the detection-plane light intensity, with standard deviation $s(t)=s_0\sqrt{1+t^2}$. In dimensionless coordinates $\mathbf{v}=\mathbf{u}/s_0$, put $\alpha =a/s_0$. ThenThe expression for $D_{\alpha }$ follows by integrating $p$ over a disk using polar coordinates. The assumptions specify a toy model with matched axial scales; its parameters are not a calibration for a real objective optical lens. It omits Airy disk rings, vector diffraction, and unequal excitation and emission wavelengths.

For a vanishingly small pinhole, the peak-normalized synthetic response tends to $(1+t^2{)}^{-2}$, whose full width at half maximum is $2\sqrt{\sqrt{2}-1}{z}_0\simeq 1.2872z_0$. With all collected emission accepted it tends to $(1+t^2{)}^{-1}$, of width $2z_0$. Meanwhile $H_{\alpha }(0)=1-e^{-{\alpha }^2/2}\to 0$ as the pinhole closes. The narrower normalized curve therefore accompanies a vanishing detected signal; peak normalization must not conceal that loss.

A symmetric confocal optical resonator has two concave mirrors of equal curvature $R$, separated by $R$. Its symmetric fundamental Gaussian beam has a central beam waist and Rayleigh range $R/2$. Its ideal transverse spectrum has degeneracies, so transmission at one resonance need not certify exclusive excitation of the fundamental normal mode.

Suppose an extended source fills every position of the entrance aperture with a cone of half-angle ${\theta }_{\mathrm{in}}$, and all these light rays must reach a one-sided planar receiver. If the receiver accepts at most a cone of half-angle ${\theta }_{\mathrm{out}}$, its available etendue bounds the admitted etendue:ConsequentlyThis original derivation combines conservation of etendue with the exact cone integral; the second inequality allows the complete output hemisphere. It bounds capacity and does not prove that a chosen shape attains the bound.

For the $30^{\circ }$ input cone with $n_{\mathrm{in}}=n_{\mathrm{out}}=1$, $C\le 4$: a $10\mathrm{mm}$ entrance cannot deliver the entire filled bundle to a circular receiver smaller than $5\mathrm{mm}$ in diameter. If the output half-angle is restricted to $60^{\circ }$, $C\le 3$ and the receiver diameter must be at least $10/\sqrt{3}\approx 5.7735\mathrm{mm}$. For a receiver medium with model refractive index $n_{\mathrm{out}}=1.5$, the full-hemisphere capacity bound instead gives $C\le 9$; that larger value also depends on transmitting the bundle into that medium.

Attempting $C=10$ in the air-to-air example forces rejection of some input. For the uniform radiance already specified, even filling the entire output hemisphere givesThe maximum mean irradiance gain remains $\eta C\le 4$. With nonuniform radiance, rejecting a given fraction of geometric etendue need not reject the same fraction of radiant flux.

These bounds use a filled extended-source bundle in passive reversible geometrical optics, comparing the same wavelength and the same radiation channels. An ideal optical point source does not fill that four-dimensional set. Diffraction, optical coherence, or coupling a Gaussian beam into an optical fiber require a physical optics calculation; the cone bound is not a universal guarantee for a single wave mode. Likewise, optical gain or fluorescence changes the propagation assumptions.

A flying focus is an optical pulse whose position of maximum focusing varies with time because different frequencies focus at different axial locations and arrive at different times. For a broad incident Gaussian beam, $f(\omega )=f_0+\beta (\omega -{\omega }_0)$ and spectral phase $\psi =\eta (\omega -{\omega }_0{)}^2$ give $z=f(\omega )$ and $t={\psi }^{\prime }(\omega )+z/c$, hence $v_f/c=[1+2\eta c/\beta {]}^{-1}$. This is the velocity of a focusing pattern, which can exceed $c$ or be negative without individual portions of the field propagating faster than light.

A fundamental Gaussian beam is an ideal beam of light whose transverse wave amplitude has a Gaussian function profile. Here consider a circularly symmetric beam in vacuum, with propagation along $z$, two transverse coordinates $x,y$, and radius $r=\sqrt{x^2+y^2}$ from the optical axis. With the beam waist at $z=0$, total power $P>0$, and vacuum wavelength $\lambda$, its time-averaged light intensity, meaning irradiance in $\mathrm{W}{\mathrm{m}}^{-2}$, isHere $w_0$ is the minimum beam radius, and $z_R$ is the Rayleigh range. Far from the beam waist, $w(z)\sim \theta |z|$ with beam divergence half-angle $\theta =\lambda /(\pi w_0)$. These are the fundamental circular Gaussian beam laws within the scalar paraxial approximation, appropriate when $w_0\gg \lambda$ and the beam divergence is small. The underlying envelope solves the parabolic wave equation exactly; this does not make it an exact vector solution of Maxwell equations. ^[ref].

The radius convention matters. Writing $U$ for the complex wave amplitude envelope, the same $w(z)$ is the wave amplitude $1/e$ radius and the light intensity $1/e^2$ radius:At $r=w(z)$ these ratios are $e^{-1}\approx 0.367879$ and $e^{-2}\approx 0.135335$, respectively. The light intensity $1/e$ radius is instead $w(z)/\sqrt{2}$. None of these radii is a sharp edge: the ideal Gaussian function has nonzero tails at every finite radius. Kogelnik and Li explicitly use the wave amplitude $1/e$ radius in their fundamental-mode treatment. ^[ref].

For a numerical example, take vacuum wavelength $\lambda =632.8\mathrm{nm}$ and beam waist radius $w_0=0.500\mathrm{mm}$. Substitution gives Rayleigh range $z_R\approx 1.24114754\mathrm{m}$ and far-field beam divergence half-angle $\theta \approx 0.402852992\mathrm{mrad}$. At $z=z_R$, the beam radius is $w=\sqrt{2}{w}_0\approx 0.707107\mathrm{mm}$ and the on-axis irradiance is half its waist value. At the same wavelength, halving the beam waist radius to $0.250\mathrm{mm}$ quarters the Rayleigh range to $0.310286885\mathrm{m}$ and doubles the beam divergence half-angle to $0.805705984\mathrm{mrad}$. At fixed power, that smaller beam waist also quadruples the peak irradiance.

The contours in the left panel are $r=\pm w(z)$ in a plane through the optical axis; they mark an light intensity level, not an outer boundary. The right panel compares the exact model radii with the far-field asymptotes $w=\theta z$. The horizontal and vertical axes use different length units, so the drawn slopes are not physical angles.

The normalization expresses conserved power during free propagation. In transverse polar coordinates, an annulus has area $2\pi rdr$, soA proof outline is the substitution $s=2r^2/w(z{)}^2$ in this integral. Spreading reduces peak irradiance as $w^{-2}$ while increasing transverse area as $w^2$, leaving power unchanged. In particular, the disk of radius $w(z)$ contains $1-e^{-2}\approx 0.8646647168$, or $86.4665\%$, of the total power; the local light intensity at its rim is only $13.5335\%$ of the on-axis light intensity. These are different quantities.

The left panel compares normalized wave amplitude magnitude and irradiance; the right panel integrates the transverse area to show enclosed power. For an ideal absorbing circular aperture of radius $a$ centered on the optical axis, the same fraction gives the power passing the aperture at that plane. An aperture of radius $w$ clips about $13.5335\%$ of the incident power; radius $2w$ clips only $e^{-8}\approx 0.000335463$, or $0.0335463\%$. After clipping, extra diffraction means the transmitted field no longer obeys the unmodified ideal Gaussian beam law. Every finite aperture cuts some ideal tail. Kogelnik and Li distinguish their large-aperture beam analysis from their finite-aperture diffraction analysis. ^[ref].

The useful tradeoff is that a small bright spot cannot remain narrow indefinitely. An optical lens can create a smaller beam waist, but diffraction then makes the beam spread faster. The particularly simple relation $w_0\theta =\lambda /\pi$ ties localization directly to beam divergence. A proof outline for the propagation law starts with $U(x,y,0)=U_0{e}^{-(x^2+y^2)/w_0^2}$. The parabolic wave equation separates into two transverse factors, and a Gaussian integral in each coordinate givesHere $q$ is dimensionless and $k$ is the wavenumber; the carrier convention is $e^{i(kz-\omega t)}$, with angular frequency $\omega$. Taking the squared complex modulus gives the stated $w(z)$ and its $w^{-2}$ light intensity scaling. Each transverse factor contributes $q^{-1/2}$, explaining the $q^{-1}$ prefactor here and the different prefactor in the preserved Gaussian beam with one transverse coordinate.

Historically, Gaussian beam theory developed alongside early laser resonator theory and lens-guided transmission. In their 1966 review, Kogelnik and Li credit Boyd and Gordon, and Goubau and Schwering, with early beam treatments; they summarize earlier work rather than claim the whole construction as a new discovery. Their paper brings beam propagation, focusing and finite-aperture diffraction into a common practical account. ^[ref].

For a circular Gaussian beam with intensity proportional to $e^{-2r^2/w(z{)}^2}$, the displayed $w(z)$ is the intensity $1/e^2$ radius. Its minimum $w_0$ is the beam waist radius, while the larger radius at a cavity mirror is not itself a waist. The Rayleigh range is $z_R=\pi w_0^2/\lambda$.

For the free parabolic wave equation $E_x=i{E}_{zz}/(2k)$, input $E(0,z)=e^{-A{z}^2}$ with $\mathrm{Re}A>0$ propagates asThe branch is continuous from $q(0{)}^{-1/2}=1$. One-dimensional transverse Fresnel propagation or two Gaussian integrals prove the formula. With $A=D^{-2}+ik/(2F)$, the squared envelope is $|q{|}^{-1}\exp [-2z^2/(D^2|q{|}^2)]$. The amplitude power is $-1/2$ because there is only one transverse coordinate. A two-transverse-coordinate Gaussian beam instead has power $-1$ when both coordinates share the same $q$.

Optical mode matching makes the spatial field of an incident beam equal to a desired optical cavity normal mode. For fundamental circular Gaussian beams, matching the complex beam parameter fixes the beam waist size and location; the optical axis, polarization and frequency must also match. A single converging optical lens connects waist Rayleigh ranges $z_A,z_B$ only when $f^2\ge z_A{z}_B$. With $s=\sqrt{f^2/(z_A{z}_B)-1}$, waist-to-lens distances obey $u=f\pm z_{A}s$, $v=f\pm z_{B}s$ with correlated signs, followed by physical distance constraints.

The incident intensity radius is $w_L=2.5\mathrm{mm}$, half the quoted diameter. For an ideal collimated Gaussian beam focused by a thin lens, the diffraction-limited waist isThe incident Rayleigh range $\pi w_L^2/{\lambda }_0\simeq 24.5\mathrm{m}$ is much larger than $f$, justifying this approximation. The pulse energy is $U=P_{\mathrm{av}}/R=0.25\mu \mathrm{J}$. Taking its duration as an effective energy-to-peak-power width gives $P_{\mathrm{pk}}\simeq U/\tau =7.14\times 10^5\mathrm{W}$. Since $I(r)=I_0{e}^{-2r^2/w_0^2}$ has integrated power $I_0\pi w_0^2/2$,If $350\mathrm{fs}$ denotes the temporal intensity full width at half maximum of a Gaussian pulse, the time integral instead gives $P_{\mathrm{pk}}=2\sqrt{\ln 2}U/(\sqrt{\pi }\tau )$, and $I_{\mathrm{pk}}=1.03\times 10^{13}\mathrm{W}\mathrm{c}{\mathrm{m}}^{-2}$. The pulse shape is unspecified, so the estimate's order-one convention is explicit. Aberration, clipping and propagation through glass were excluded from this calculation as stipulated.

For bulk electric-dipole second-harmonic generation, the crystal must lack inversion symmetry and have a nonzero relevant component of its second-order nonlinear optical susceptibility. It should transmit both pump and harmonic, tolerate the required intensity, and allow phase matching by suitable birefringence or another phase-compensation mechanism. A large nonlinear coefficient, low optical loss and good optical quality are desirable. Noncentrosymmetry permits a response but does not guarantee one in every chosen geometry.

In the undepleted plane wave limit, contributions from successive layers give a harmonic amplitude proportional to ${\int }_0^L{e}^{i\Delta kz}dz$, where $\Delta k=k_{2\omega }-2k_{\omega }$. HenceTune orientation, temperature and pump polarization to maximize the nonlinear coupling and cancel phase mismatch. For continuous operation, an impedance-matched enhancement optical resonator, or intracavity doubling, raises the circulating intensity. Focus the Gaussian beam to balance intensity against diffraction and spatial walk-off; arbitrarily tight focusing shortens the useful overlap. Antireflection surfaces, stable resonance locking and efficient heat removal reduce avoidable loss and drift.

A longitudinal mode makes an integer round-trip phase accumulation. For the same transverse mode and an empty nondispersive cavity, increasing that integer by one changes the optical frequency by the free spectral rangeFrequency-dependent mirror phases or a dispersive intracavity medium would change the round-trip group delay used here.

For a paraxial ray, use the height-slope column $(x,\theta {)}^T$. Free propagation and reflection from a concave mirror are represented, in unfolded coordinates, byThe reflection changes the slope by $-2x/R$. Because the mirrors are identical, the same half-round-trip matrix is applied repeatedly:Its two eigenvalues are reciprocal. They lie on the unit circle without a growing component when $|1-l/R|<1$, whereas $l>2R$ gives a real eigenvalue of modulus greater than one. At $l=2R$, the repeated eigenvalue has a nontrivial Jordan block and generic rays drift out. Thus the nondegenerate stable interval isThe confocal value $l=R$ is included: there $A^2=-I$, so the repeated round trip is bounded. Applying a strict trace inequality to the full-round-trip matrix without checking this exact matrix would incorrectly discard that point. These are the paraxial conditions for a symmetric concave-mirror resonator.

Symmetry places the Gaussian beam waist halfway between the mirrors. Its Rayleigh range is $z_R=\pi w_0^2/\lambda$, and its wavefront radius at distance $z$ from the waist is $R_b(z)=z+z_R^2/z$. Matching that wavefront to the mirror at $z=l/2$ givesHenceA real positive Rayleigh range requires the same stable interval. At its limiting geometries the finite confined Gaussian mode degenerates; the formula is not a way to realize an arbitrarily small physical waist without leaving the paraxial model.

For a paraxial Gaussian beam, choose the complex beam parameter withHere $R_b$ is the wavefront radius and $w$ the intensity $1/e^2$ radius. The ABCD law isFor a real, lossless round-trip ray transfer matrix $M$ in free space, $AD-BC=1$. A reproducing cavity mode requires $q=(Aq+B)/(Cq+D)$, henceExcept for a scalar identity matrix, a finite physical solution requires a complex-conjugate root pair, one with positive imaginary part. The generic strict stability condition is thereforeEquivalently, the ray eigenvalues are distinct unit-modulus numbers $e^{\pm i\mu }$, so repeated round trips do not cause unbounded ray excursions. Adequate mirror apertures and longitudinal resonance are also needed for a low-loss resonant mode.

The strict inequality has a degenerate exception relevant to this paper. If $M=-I$, then $q^{\prime }=q$ for every physical $q$ even though $A+D=-2$. A symmetric confocal resonator has exactly this matrix and supports Gaussian modes; it is not excluded as a physical cavity. Thus the printed strict condition should be understood as the generic, nondegenerate stability criterion, not a universal necessary condition covering the confocal case in c.

Put $s=z-z_0$. In the source's complex beam parameter convention,Substitution separates the transverse wave amplitude and phase:Thus the Gaussian beam is circularly symmetric, has a Gaussian function transverse light intensity, and narrows to a beam waist at $z_0$ before spreading by diffraction. Its wavefront is planar at the beam waist, converging before it and diverging after it. The quadratic transverse phase identifies the signed radius of curvature:Define spot size as the wave amplitude $1/e$ radius, equivalently the light intensity $1/e^2$ radius. ThenThe Rayleigh range is where $w=\sqrt{2}{w}_0$; the far-field beam divergence half-angle is $\lambda /(\pi w_0)$. The factor $1/q$ also produces a longitudinal phase shift relative to a plane wave, changing by $\pi$ across the beam waist; its sign here follows the given $e^{i(kz-\omega t)}$ convention. The on-axis light intensity decreases as $1/(s^2+z_R^2)$ while transverse integration leaves the power constant. These are scalar paraxial approximation results.

Place the optical lens at $z=0$. For a sufficiently broad incident Gaussian beam, the beam waist positions approach the focal lengths frequency by frequency:Add the group delay at the optical lens to propagation through vacuum:Here an arbitrary common time origin has been omitted. Differentiating both equations with respect to angular frequency yieldsFor $\beta >0$, define $u=2\eta c/\beta$. The velocity graph in the preceding figure has horizontal asymptote zero and a vertical asymptote at $u=-1$. For $\eta >0$, $0<v_f<c$: the higher-frequency components both enter later and travel farther. At $\eta =0$, $v_f=c$. For $-\beta /(2c)<\eta <0$, $v_f>c$: the farther-focusing components leave earlier, partly compensating their extra travel. At $\eta =-\beta /(2c)$ every component focuses at the same time, so the trajectory is not a finite-velocity function of time. For $\eta <-\beta /(2c)$ the farther positions focus first, giving $v_f<0$; it tends to zero from below as $\eta \to -\infty$.

A flying focus is a moving pattern assembled from different portions of an optical pulse. Each portion propagates causally at $c$ in this model. The superluminal or backward motion of the focus therefore does not describe a material object or information overtaking light. The calculation assumes negligible additional frequency-dependent group delay from the optical lens; any such delay must be included in $\psi$.

Use the complex beam parameter conventionHere $R_b$ is the wavefront curvature, distinguished from the mirror curvature $R$. In an unfolded propagation coordinate, a concave mirror has the ray transfer matrix $\left(\begin{array}{cc}1& 0\\ -2/R& 1\end{array}\right)$: a ray at height $x$ changes direction by $-2x/R$. The return to the original plane is thereforeThe Gaussian beam transformation $q^{\prime }=(Aq+B)/(Cq+D)$ givesFor $R=2L$, $M=\left(\begin{array}{cc}0& L\\ -1/L& 0\end{array}\right)$ and $q^{\prime }=-L^2/q=i{L}^2/z_R$. The returning Gaussian beam consequently has a beam waist at the same plane, withThis alone does not make it the same Gaussian beam; equal Rayleigh ranges require $z_R=L$.

Adding the plane mirror requires a self-consistent optical resonator normal mode, $q^{\prime }=q=i{z}_R$. Since $A=D$, the fixed-point equation reduces to $C{q}^2=B$, givingFor a finite stable Gaussian beam, $0<L<R$. The same formula remains valid at $R=2L$, where self-consistency selects $z_R=L$. At the geometrical hemispherical limit $R=L$ it predicts $z_R=0$; that marginal boundary has no finite-width stable fundamental Gaussian beam. Thus the source's word “hemispherical” is used for the plano-concave family, including the near-hemispherical case below.

The source optical resonator has Rayleigh rangeThe symmetric target optical resonator has its fundamental Gaussian beam beam waist halfway between the mirrors. At distance $\ell =L_2/2$ from that beam waist, wavefront curvature is $R_2=\ell +z_B^2/\ell$. ThusHere $f^2=2.25{\mathrm{m}}^2>z_A{z}_B=0.866025{\mathrm{m}}^2$, so b gives two real branches:Only the first branch has the target beam waist downstream of the optical lens. It is also geometrically accessible: the optical lens is $u-L_1\simeq 0.62642\mathrm{m}$ beyond the source output mirror, and the target input mirror is $v-L_2/2\simeq 7.44786\mathrm{m}$ beyond the optical lens. The total separation between beam waists is $u+v\simeq 14.07428\mathrm{m}$.