Gaussian Beam Theory
The Gaussian Beam node releases rays with a Gaussian distribution of intensity or power. It can only be used when solving for either intensity or power, or both.
-
For ray propagation over relatively short distances, choose Collimated source (large Rayleigh range) from the Beam type list.
-
For ray propagation over relatively large distances, choose Point source (small Rayleigh range) from the Beam type list.
Effect of Space Dimension
The Gaussian Beam node is available in 3D, 2D axisymmetric, and 2D geometries. The interpretation of the beam is different in 2D than in 3D or 2D axisymmetric models.
In 2D axisymmetric models, the irradiance distribution is still assumed to be rotationally symmetric about the beam axis, but the beam axis is restricted to either the positive or negative z direction. If any other direction were chosen as the beam axis, the beam would not actually be Gaussian in the corresponding 3D solid of revolution.
Gaussian Beam Equation in 3D
Without loss of generality we assume the beam is propagating in the positive z direction with a waist centered at the origin, although in a 3D geometry any other waist position and beam axis may be specified.
The electric field amplitude of a paraxial Gaussian beam is most conveniently expressed in cylindrical polar coordinates,
(3-11)
where
ρ (SI unit: m) is the radial coordinate (with the beam centered around ρ = 0),
z (SI unit: m) is the axial coordinate (chosen so that z = 0 is the beam waist),
E(ρz) (SI unit: V/m) is the electric field at an arbitrary location (ρz),
E0 = E(0, 0) is the electric field at the center of the beam waist, and
k (SI unit: rad/m) is the wave number in the medium through which the beam propagates, k = 2π/λ = 2πn/λ0 where n (dimensionless) is the absolute refractive index of the medium and λ0 (SI unit: m) is the vacuum wavelength.
The beam parameter w0 is called the waist radius. Of course, there is no well-defined cutoff distance at which the electric field goes to zero, so there is a certain degree of arbitrariness involved in assigning a value to the beam radius. As will be shown in the following sections, the beam waist radius is defined such that the fraction of the total beam power that passes through a circle or radius w0 centered at the origin is 1 − e2.
The remaining functions in Equation 3-11 are defined as follows:
The function w(z) is the beam radius at any axial coordinate z,
The function R(z) describes how quickly the beam phase varies with radial position in a plane located at any axial coordinate z,
The term η(z) is an additional phase delay known as the Guoy phase,
In these definitions, zR (SI unit: m) is the Rayleigh range,
The Rayleigh range is the axial distance from the beam waist at which the radius has increased by a factor of 2, or equivalently the distance from the beam waist at which the cross sectional area of the beam has doubled.
The exponent in Equation 3-11 can be separated into amplitude and phase terms,
where δ (SI unit: rad) is the phase factor,
(3-12)
Gaussian Beam Equation for Irradiance
It is also possible to write Equation 3-11 in terms of the irradiance rather than electric field amplitude. As the irradiance I (SI unit: W/m2) is proportional to E ⋅ E, the phase factor vanishes,
(3-13)
Asymptotic Limits of the Paraxial Gaussian Beam
In the asymptotic limit as the z-coordinate is much greater than the Rayleigh range, the beam radius approaches a linear function of z,
so a beam with radius w(z) approaches a biconical shape (with both cone vertices meeting at the origin) in the limits z << −zR and z >> zR. These asymptotes each make an angle of α (SI unit: rad) with the beam axis, where
α is called the beam divergence half-angle, to avoid confusion with references such as Ref. 8 that may define a full angle equal to 2α.
Making the substitution w(z) ≈ w0z/zR into Equation 3-13 then yields
At this point it is convenient to convert to spherical polar coordinates,
where r (SI unit: m) is the distance from the origin and θ (SI unit: rad) is the polar angle measured from the positive z-axis. Thus
Note also from the definitions of the Rayleigh range and the divergence half-angle that
Thus
Finally, recall that the Gaussian beam is paraxial (small θ). Taylor series expansions of the cosine and tangent functions reveal that
So that for small θ,
(3-14)
Equation 3-14 shows that the far-field limit of a Gaussian beam (z >> zR) is a cone of half-angle α in which the ray intensity or power is proportional to a Gaussian function of the angle between each ray and the beam axis.
Asymptotic Limit of Phase
The Geometrical Optics interface includes the option to treat rays as fully coherent wavefronts in which the instantaneous phase is solved for, by selecting the Compute phase check box. When the rays are treated as a collimated source, their phase is initialized as if all rays are released at locations along the beam waist. Substituting z = 0 into Equation 3-12 for the phase factor simply yields δ = 0, so the initial phase of all released rays is equal.
Now consider once again the asymptotic limit of z >> zR.
So for large z, Equation 3-12 becomes
Again converting to spherical polar coordinates,
And in fact a Taylor series expansion shows that
So the phase of the released rays can be considered uniform, so long as
where L is the maximum distance the rays propagate in the model.
Total Beam Power and Encircled Power
The total power transmitted by the Gaussian beam Psrc (SI unit: m) is found by integrating the irradiance from Equation 3-13 over a plane perpendicular to the z-axis,
As expected, the total power crossing a plane perpendicular to the z-axis does not depend on the value of z, otherwise energy would be accumulating at some location in the beam.
Compare the total source power Psrc with the encircled power Pcir (SI unit: W/m2) that passes through a circle of radius w(z) centered at the z-axis,
Thus, the radio of the encircled power to the total source power is
In addition, by substituting ρ = 0 and ρ = w(z) into Equation 3-13, it follows that
Thus the beam radius w(z) is sometimes called the 1/e2 radius. If the above integral is repeated, replacing the upper limit integration with cw(z), where c is a positive constant (called the Beam cutoff ratio in the settings window), the result is
Summary of Beam Parameters (3D)
When the wavelength of the beam λ and any one of the three beam parameters (waist radius w0, divergence half-angle α, and Rayleigh range zR) are defined, the other two beam parameters can be derived:
For a beam of given waist radius w0 and any of the three intensity-related parameters (intensity at the beam waist I0, total source power Psrc or encircled power Pcir), the remaining two parameters can be derived:
Gaussian Beam Equation in 2D
The Gaussian beam in 2D has a different functional form because instead of having rotational symmetry about the beam axis, it is uniform in the out-of-plane direction,
Comparing the 2D form with the 3D form from Equation 3-11, the key differences are a square root on the term w0/w(x) and a factor of 1/2 in the Guoy phase. In this expression x has been considered the propagation direction but the beam may actually propagate in any direction in the xy-plane.
The corresponding expression for the beam irradiance is
The total beam power per unit distance in the out-of-plane direction is
The 2D analogue to the encircled power from the previous section, the beam power per unit distance that passes between coordinates cw(x) and cw(x) is
Because the irradiance is an even function of y, this may be simplified to
Thus the ratio of the encircled power to the total source power is
where the error function erf is a built-in function in COMSOL,
For c = 1 the encircled power is about 95% of the total beam power.