Gaussian Beam Optics
Gaussian Beam Optics
Gaussian Beam Optics
[HechtCh.13.1pages594596NotesfromMellesGriotandNewport]
Readings:
For details on the theory of Gaussian beam optics, refer to the excerpts from the Melles Griot and
Newport catalogs. Melles Griot, anlong with Newport Corporation, is a major manufacturer of
optical components used for research. You will also find tutorials on other subjects of optics and
photonics at their web site:
http://www.mellesgriot.com/products/optics/toc.htm (CVI Melles Griot Optics Guides)
http://www.newport.com/servicesupport/Tutorials/default.aspx?id=110 (Newport Tutorials).
GaussianBeamCalculators:
There are many calculators or programs to determine Gaussian beam propagation in free space or
through a combination of lenses. Below is a short list of program you may use to solve the
homework problems.
Gaussian Beam Calculator for simple lenses
http://www.photonics.byu.edu/Gaussian_Beam_Propagation.phtml
http://www.newport.com/OpticalAssistant/FocusCollimatedBeam.aspx
http://www.originalcode.com/downloads/GBC8p2.zip (need LabView)
http://www.novajo.ca/abcd/ (for Mac only, ABCD Gaussian Beam Propagation)
Gaussian Beam Intensity/Irradiance
http://www.calctool.org/CALC/phys/optics/gauss_power_dist
Gaussian Beam and Spatial Filtering
http://www.newport.com/OpticalAssistant/SpatialFilterPinhole.aspx
Introduction
One usually thinks of a laser beam as a perfectly collimated beam of light rays with be beam
energy uniformly spread across the cross section of the beam. This is not an adequate picture for
discussing the propagation of a laser beam over any appreciable distance because diffraction
causes the light waves to spread transversely as they propagate, Fig. 1.
Additionally, the energy (irradiance) profile of a laser beam is typically not uniform. For the
most commonly used He-Ne lasers (operating in the TEM00 mode) the irradiance (the power
carried by the beam across a unit area perpendicular to the beam = W/m2) is given by a Gaussian
function:
2
2
2 P 2 r 2 / w2
I ( r ) = I 0 e2 r / w =
e
, (1)
w2
1
where w is defined as the distance out from the center axis of the beam where the irradiance
drops to 1/e2of its value on axis. P is the total power in the beam. r is the transverse distance from
the central axis. w depends on the distance z the beam ahs propagated from the beam waist. w0 is
the beam radius at the waist. [The beam waist is defined as the point where the beam wave front
was last flat (as opposed to spherical at other locations).] For a hemispherical laser cavity such as
the one used for the He-Ne laser used in the lab, the waist is located roughly at the output mirror.
w0 is related to w by
Experimentally, one can use a CCD detector array to measure how the irradiance various across
the beam for several values of z>>zR. Then fit the data for each z using Eq. (1), which will yield
values for w(z). Then one can use Eqs. (3) and (4) to determine w0 and calculate zR.
plane (i.e. the two lenses are separated by f0+fe). The rays reaching the eye are again parallel,
but appear to subtend a much larger angle than the original object. From Fig. 2 it is easy to see
that the angular magnification is
Beamexpander:
Because Gaussian beams do not follow the rules of ray optics, we cannot use the lens equation to
design a beam expander. However, as discuss in the Melles Griot Optics Guide, if you consider
the object to be the beam waist of the incoming beam and the image to be the beam waist after
the beam passes through the lens, then you can use a modified lens equation:
Lets now apply this to an inverted astronomical telescope with the focal length of the first lens
being 5 cm and the second 40 cm. In the astronomical telescope the two lenses are separated by
the sum of the focal lengths of the two lenses 45 cm in this case. We assume a red He-Ne laser
(633nm) with beam waist radius of 0.4 mm. We first use Eq. (3) to get zR=0.80 m. For the first
lens, f = 5 cm, and the beam waist for the laser is close to the exit of the laser. We put the lens as
close to the laser as possible and assume s=0. The using Eq. (5) we get
Material Properties
Optical Specifications
Fundamental Optics
2ch_GuassianBeamOptics_f_v2.qxd
6/6/2005
12:51 PM
Page 2.2
Optical Coatings
www.mellesgriot.com
2.2
For virtually all laser cavities, the propagation of an electromagnetic field, E(0), through one round trip in an optical resonator
can be described mathematically by a propagation integral, which
has the general form
E (1) ( x, y ) = e jkp
InputPlane
(2.1)
K ( x, y, x0 , y0 ) E nm x0, y0 dx0dy0
InputPlane
(2.2)
k x2 + y 2
(2.3)
2ch_GuassianBeamOptics_f_v2.qxd
6/6/2005
12:51 PM
Page 2.3
100
60
40
I ( r ) = I 0e 2r
20
13.5
41.5w 4w
Figure 2.1
0
CONTOUR RADIUS
/ w2
1.5w
The invariance of the form of the distribution is a special consequence of the presumed Gaussian distribution at z = 0. If a uniform irradiance distribution had been presumed at z = 0, the pattern
at z = would have been the familiar Airy disc pattern given by a
Bessel function, whereas the pattern at intermediate z values would
have been enormously complicated.
Simultaneously, as R(z) asymptotically approaches z for large
z, w(z) asymptotically approaches the value
w (z) =
lz
p w0
w (z)
l
=
.
z
p w0
(2.8)
(2.4)
(2.7)
Material Properties
(2.6)
Optical Specifications
2P 2r 2 / w 2
,
e
pw 2
where w=w(z) and P is the total power in the beam, is the same at
all cross sections of the beam.
pw2
R ( z ) = z 1 + 0
lz
PERCENT IRRADIANCE
80
Fundamental Optics
the radius of the 1/e2 contour after the wave has propagated a distance z, and R(z) is the wavefront radius of curvature after propagating a distance z. R(z) is infinite at z = 0, passes through a minimum
at some finite z, and rises again toward infinity as z is further
increased, asymptotically approaching the value of z itself. The
plane z=0 marks the location of a Gaussian waist, or a place where
the wavefront is flat, and w0 is called the beam waist radius.
direction
of propagation
and
lz 2
w ( z ) = w0 1 +
p w02
1/ 2
(2.5)
Figure 2.2
Optical Coatings
where z is the distance propagated from the plane where the wavefront is flat, l is the wavelength of light, w0 is the radius of the 1/e2
irradiance contour at the plane where the wavefront is flat, w(z) is
2.3
Fundamental Optics
2ch_GuassianBeamOptics_f_v2.qxd
6/6/2005
12:51 PM
Page 2.4
between near-field divergence and mid-range divergence, is the distance from the waist at which the wavefront curvature is a maximum.
Far-field divergence (the number quoted in laser specifications)
must be measured at a distance much greater than zR (usually
>10#zR will suffice). This is a very important distinction because
calculations for spot size and other parameters in an optical train
will be inaccurate if near- or mid-field divergence values are used.
For a tightly focused beam, the distance from the waist (the focal
point) to the far field can be a few millimeters or less. For beams coming directly from the laser, the far-field distance can be measured
in meters.
w
w0
e2
irradiance surface
ne
tic co
pto
asym
Typically, one has a fixed value for w0 and uses the expression
w0
lz 2
w ( z ) = w0 1 +
p w02
z
w0
Optical Specifications
Figure 2.3 Growth in 1/e2 radius with distance propagated away from Gaussian waist
Material Properties
to calculate w(z) for an input value of z. However, one can also utilize this equation to see how final beam radius varies with starting
beam radius at a fixed distance, z. Figure 2.5 shows the Gaussian
beam propagation equation plotted as a function of w0, with the particular values of l = 632.8 nm and z = 100 m.
1/ 2
The beam radius at 100 m reaches a minimum value for a starting beam radius of about 4.5 mm. Therefore, if we wanted to achieve
the best combination of minimum beam diameter and minimum
beam spread (or best collimation) over a distance of 100 m, our
optimum starting beam radius would be 4.5 mm. Any other starting value would result in a larger beam at z = 100 m.
We can find the general expression for the optimum starting
beam radius for a given distance, z. Doing so yields
(2.9)
At the beam waist (z = 0), the wavefront is planar [R(0) = ]. Likewise, at z=, the wavefront is planar [R()=]. As the beam propagates from the waist, the wavefront curvature, therefore, must
increase to a maximum and then begin to decrease, as shown in
figure 2.4. The Raleigh range, considered to be the dividing line
lz
w0 (optimum ) =
p
.
1/ 2
(2.10)
Using this optimum value of w0 will provide the best combination of minimum starting beam diameter and minimum beam
z=q
planar wavefront
laser
2w0 2
z=0
planar wavefront
Gaussian
profile
2w0
z = zR
Optical Coatings
maximum curvature
Figure 2.4
2.4
Gaussian
intensity
profile
2ch_GuassianBeamOptics_f_v2.qxd
6/6/2005
12:51 PM
Page 2.5
beam waist
2 w0
80
beam expander
60
40
w(zR) = 2w0
20
0
10
w(zR) = 2w0
zR
zR
APPLICATION NOTE
.
w ( z ) = 2 (w0 )
Fundamental Optics
100
p w02
l
with
w ( zR ) = 2w0 .
BEAM EXPANDERS
Melles Griot offers a range of precision beam
expanders for better performance than can be
achieved with the simple lens combinations shown
here. Available in expansion ratios of 3#, 10#,
20#, and 30#, these beam expanders produce
less than l/4 of wavefront distortion. They are
optimized for a 1-mm-diameter input beam, and
mount using a standard 1-inch-32 TPI thread. For
more information,
see page 16.4.
Optical Coatings
This result can now be used in the problem of finding the starting beam radius that yields the minimum beam diameter and beam
spread over 100 m. Using 2(zR) = 100 m, or zR = 50 m, and
l = 632.8 nm, we get a value of w(zR) = (2l/p) = 4.5 mm, and
w0 = 3.2 mm. Thus, the optimum starting beam radius is the same
as previously calculated. However, by focusing the expander we
achieve a final beam radius that is no _larger than our starting beam
radius, while still maintaining the 2 factor in overall variation.
Material Properties
Do you need . . .
2.5
O p t i c s 555
2 r2
I r = I0 exp 2
e-2
0
0.59
0.71
0.83
0.865
0.01
1.0
2
1
1.52
NORMALIZED RADIUS (r/0)
ENCIRCLED POWER
RELATIVE INTENSITY
0.5
e-1
Figure 1
2r2
P r = P 1 exp 2
() ( )
ACCESSORIES
x=0
R(x)
Figure 2
TECHNICAL REFERENCE
0.23
0.5
() ( )
1.0
()
()
POLARIZATION OPTICS
2r 2
IS = ES ES* = E0 E0* exp 2
()
P = 0 I 0
2
2
I 0 =P
BEAMSPLITTERS
r2
ES = E0 exp 2
WINDOWS
2
x
(x) = 1 + 2
0
2 2
R(x) = x 1 + 0
x
POLARIZATION OPTICS
BEAMSPLITTERS
WINDOWS
556 O p t i c s
2
0
ACCESSORIES
xR =
(x)
TECHNICAL REFERENCE
x=0
R(x)
F
Figure 3
qout =
[q A + B]
[q C + D]
in
( ) R(x)
qx
()
w x
in
4
2w0
O p t i c s 557
D
= (f /#)1
F
4 F
2w 0 =
D
8
632,8 nm
10 mm
,
1 mm
TECHNICAL REFERENCE
ACCESSORIES
8 F
DOF =
D
POLARIZATION OPTICS
mm
) 101mm
,
BEAMSPLITTERS
4
632,8 nm
WINDOWS