Laser Physics
Laser Physics
LASER
PHYSICS
Lecture notes
VILNIUS UNIVERSITY PRESS
Faculty of
Physics
LASER
PHYSICS
Lecture notes
VILNIUS
UNIVERSITY
PRESS
2023
Reviewed by:
assoc. prof. dr. Domas Paipulas (Vilnius University)
Contents
1 Principles of lasers in a nutshell 6
1.1 What is a Laser? . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2 Historical perspective . . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Light Amplification by Stimulated Emission of Radiation . . 13
1.4 Quantifying laser radiation . . . . . . . . . . . . . . . . . . . . 16
1.5 Laser safety classes . . . . . . . . . . . . . . . . . . . . . . . . 25
3 Optical resonators 49
3.1 Light rays and ray (or ABCD) matrices . . . . . . . . . . . . . 49
3.2 Laser resonator stability . . . . . . . . . . . . . . . . . . . . . 53
3.3 Hermite-Gaussian beams and Laguerre Gaussian beams . . . . 57
3.4 Resonance frequencies . . . . . . . . . . . . . . . . . . . . . . 61
3.5 Fabry-Perot resonator . . . . . . . . . . . . . . . . . . . . . . 64
3.6 Losses in resonators . . . . . . . . . . . . . . . . . . . . . . . . 68
9 EN - LT Glossary 165
This page is intentionally left blank
1 Principles of lasers in a nutshell 6
Figure 1: Comparison of light patterns coming from three different light sources.
The left image is adapted from Wikimedia under GFDL 1.2 license.
Stark differences between light coming from a laser and the other two sources
are immediately visible allowing to recognize distinctive features of laser ra
diation:
• Laser light has low divergence (is highly directional): laser
beam is almost aligned to the direction of propagation, with little an
gular spread from the point of origin. In fact, if the authors had chosen
a He-Ne laser instead of a cheap laser pointer, the laser beam in the
photo would have been even smaller, i.e., with even lower divergence.
For He-Ne laser typical full angle beam divergence can be as low as
1 mrad.
• Laser light is nearly monochromatic: in general, laser emits light
with very narrow spectrum. When it passes through a dispersive prism,
1.1 What is a Laser? 7
These properties make laser a unique light source, which offers a wealth of
specific applications not possible with any other light source. To mention a
few: laser pointers make use of low divergence of laser light; laser-induced
fluorescence utilizes the fact that laser radiation is monochromatic, so we can
efficiently excite fluorescence of only the desired molecules; interferometry,
holography and optical sensors rely on coherence of laser light.
Summarizing the above, a following definition of a laser could be given:
Laser is a light source producing low divergence, coherent and monochro
matic light in the optical region of the electromagnetic spectrum.
the beams were overlapped on the screen. Although none of the presented ex
amples can match low divergence, monochromaticity and coherence of lasers,
they prove that the initial definition of a laser is rather quantitative, depend
ing on quantitative description of what we call low divergence, coherence and
monochromaticity.
A truly accurate description of what is a laser comes from the principle
of its operation, i.e. how light with such properties is generated:
ing to Boltzman law and absorption is the dominating effect. The idea of
optical amplification using stimulated emission in gases was first formulated
by Russian physicist V. Fabricant, but it still was not clear how to realize
this idea in practice. After World War II, W. Lamb and R. C. Retherford
noticed that population inversion can be achieved during Nuclear Magnetic
Resonance (NMR) and later stimulated emission of radio waves was observed.
Population inversion at that time was interpreted as ”negative absolute tem
perature” as would follow from the Boltzmann distribution. In 1951 C. H.
Townes conceived an idea that electromagnetic radiation could be ampli
fied via stimulated emission in a resonator, and in 1954 Townes and J.
P. Gordon invented the first device operating on that principle, the maser,
which emitted electromagnetic radiation in the microwave (millimeter wave)
range. The term ”maser” is an acronym for Microwave Amplifier based on
Stimulated Emission of Radiation. In 1964 A. Penzias and R. Wilson using
maser as a sensitive microwave detector, discovered the cosmic microwave
background radiation, which is the relic radiation from the Big Bang. For
this fundamental discovery they received the Nobel Prize in Physics in 1978.
In 1958 C. H. Townes and A. L. Schawlow theoretically predicted the
possibility to amplify electromagnetic radiation in the optical range and put
forward the key principles of laser generation3 . They also foresaw that a
Fabry-Perot etalon could be used as a simple resonator in the optical range.
What was still missing, was to find suitable materials and methods to cre
ate population inversion, which became a subject of very intense research by
many groups at the time. At present, it is widely recognized that the first
laser was invented in 1960 by T. Maiman, when he reported the laser
action in a ruby crystal (where natural chromium impurities exist) with alu
minum coated ends serving as a resonator mirrors. To achieve the population
inversion, Maiman used a white light flashlamp as a pump source. The word
”recognized” is used here by purpose, since Maiman described his results very
vaguely and it might be that he only observed the stimulated emission below
the laser generation threshold. These doubts are also based on the fact that
Physical Review Letters rejected Maiman’s manuscript, but he managed to
publish a short message in Nature 4 . It can be claimed with a confidence that
the first laser, where light amplification and directional radiation in the form
of a coherent beam, was reported a few months later by R. J. Collins, D.
F. Nelson, A. L. Schawlow, W. Bond, C. G. B. Garrett and W. Kaiser5 . At
3
A. L. Schawlow and C. H. Townes, Infrared and optical masers, Physical Review112,
1940 (1958).
4
T. H. Maiman, Stimulated optical radiation in ruby, Nature 187, 493-494 (1960).
5
R. J. Collins, et al., Coherence, narrowing, directionality, and relaxation oscillations
in the light emission from ruby, Physical Review Letters 5, 303 (1960).
1.2 Historical perspective 11
sion of Radiation”, highlighting the fact that, unlike masers which amplify
microwaves, lasers amplify the electromagnetic radiation in the optical range.
Interestingly, the term ”laser” was first used in 1959 in a conference paper8 by
a graduate student G. Gould, who at that time worked on a doctoral thesis
about the energy levels of excited atoms of thallium. Two years before, he
coined some general ideas about how lasers could be built and where they
could be used in a lab notebook which he notarized. In 1959 he filed a patent
application with the high-tech company he worked in. The U. S. Patent Office
denied his application and granted it to Bell Labs, which resulted in a twenty
eight year-long ”patent war”. Eventually, in 1987 G. Gould won the first pat
ent lawsuit and was issued patents for optically pumped and gas-discharged
lasers. However, the historical question about assigning the credits for laser
invention is difficult since patents tell little about how the ideas arose and
spread among researchers. C. H. Townes, N. Basov and A. Prochorov were
awarded the Nobel Prize in Physics in 1964 for their input which led to in
vention of the laser. Another key researcher in the field, A. Schawlow, was
awarded the Nobel Prize in Physics in 1981 for laser applications in spectro
scopy. The beginning of the laser era in 1960 gave birth to a new branch of
optics – nonlinear optics. Since the invention of the laser, a huge variety of
new laser materials was developed, the power of lasers increased enormously
and the duration of emitted pulses decreased significantly, yielding a coher
ent light with unprecedented intensity. The very first pulsed lasers generated
light pulses with a duration of a few microseconds, while modern ultrashort
pulse lasers can provide light pulses with a few femtosecond duration – a
billion of times shorter! The pulse energy and power grew tremendously,
and at present lasers are reliable turn-key devices, without which modern
physics, biology, chemistry, medicine and materials science laboratories and
future progress of these sciences is unimaginable. Interestingly, the laser was
not invented for a specific purpose, the skeptics even claimed that invention
of the laser was ”not a solution of a specific problem, but, on the contrary,
a solution that created many problems”. Despite this, lasers gradually and
imperceptibly became an integral part of many scientific, industrial, techno
logical, military and many other fields, so it can be concluded that laser is
the most important optical device invented in the last six decades.
8
G. Gould, The LASER, Light Amplification by Stimulated Emission of Radiation, in
The Ann Arbor Conference on Optical Pumping (University of Michigan, USA, June 15-18
1959)
1.3 Light Amplification by Stimulated Emission of Radiation 13
and the phase, direction and polarization of the emitted light wave are ran
dom. Apart from the two aforementioned processes, yet another interaction
process is also possible. This process is called stimulated emission. In this
case when an electron is in the higher energy orbital (the atom is excited),
another photon can trigger the electron to jump to the ground orbital state
emitting a photon. The stimulated emission process occurs only when the
triggering photon energy is equal to the energy difference between higher
energy state and ground state orbitals. More importantly, the emitted and
corresponding light wave has identical phase, direction, polarization and
wavelength to the initial triggering photon.
Laser operation is based on light amplification by stimulated emission of
radiation. Let us discuss this in more detail.
Consider that the lasing medium is at ground state (Fig.7a). When a pump
source (e.g., an electric flash) is turned on, some energy is transferred to
the lasing medium and some atoms become excited (denoted by red circles
in Fig.7b). The first key moment is spontaneous emission that starts from
the excited atoms and the emitted photons can then trigger stimulated emis
sion of other excited state atoms (Fig.7c). The second important moment
is that the amount of excited electrons (population of an excited energy
level) must be greater than the amount of electrons in the lower (ground
state) energy level. Such state is called population inversion. The term
”inversion” indicates that this kind of population is impossible to achieve un
der usual conditions: according to Boltzman distribution, the population of
higher energy levels decreases exponentially. When the emitted light (both
1.4 Quantifying laser radiation 16
Figure 8: Time evolution of light electric field emitted from pulsed (black lines)
and continous wave (CW) lasers (red line). τp is the pulse duration, tr is the pulse
repetition period (inversely proportional to pulse repetition rate).
For lasers that operate in CW regime we can define these main paramet
ers:
9
C. Walther, M. Fischer, G. Scalari, R. Terazzi, N. Hoyler, and J. Faist,”Quantum
cascade lasers operating from 1.2 to 1.6 THz”, Applied Physics Letters 91, 131122 (2007)
1.4 Quantifying laser radiation 18
Figure 10: Optical spectrometer (left) and example of measured He-Ne laser
spectrum (right).
Figure 11: Thermal power meter. It consists of two connected parts: power sensor
and power meter device which displays the measured power and is used to control
the power meter.
laser radiation. Such power meters are very fast, broadband (response
is independent of wavelength in a broad spectral range, usually from
200 nm to 20000 nm) and do not require thermal equilibrium with
the environment. However, pyroelectric power meters are fragile, less
accurate, cannot measure CW laser radiation power and are sensitive to
vibrations. Therefore, they are usually used to measure low repetition
rate laser power for each pulse.
Figure 12: Laser beam size measured at different intensity levels: w1/e2 – radius
at 13.6% of maximum intensity level; wF W HM – radius measured as half of beam
diameter at full width at half maximum intensity level.
The beam diameter can differ in the horizontal (x) and vertical (y)
directions; such beams are called elliptical since the beam spot on a
screen looks like an ellipse instead of a circle. Commonly used intens
ity level values: FWHM – full width at half maximum, 1/e2 – 1/e2
(13.6%) of maximum intensity level, 1/e – 1/e (36.8%) of maximum
electric field amplitude level. There is also 4σ notation defined as 4
times the standard deviation of the horizontal or vertical marginal dis
tributions respectively. It is the ISO international standard definition
for beam width. Despite the fact that 4σ notation is the international
ISO standard for beam width definition, other beam width definitions
are also commonly used. Beam size can be measured in several differ
ent ways but the simplest one is to use a charge-coupled device (CCD)
camera which displays laser beam intensity distribution and estimates
its size (Fig.13). Professional beam width measurement software usu
ally provides beam width values at various intensity levels. Although
1.4 Quantifying laser radiation 21
Figure 13: CCD camera (left) and example of measured laser beam intensity pro
file in 2D (top right) and 3D (bottom right). Color coding indicates the intensity.
P ower
I= (1)
Area
Figure 14: Electric field oscillations in an ultrashort light pulse and its duration
estimation at FWHM level (τp ). Thick blue line above the oscillations represents
what is called the pulse envelope.
10
T. Gaumnitz et al, Streaking of 43-attosecond soft-X-ray pulses generated by a pass
ively CEP-stable mid-infrared driver, Optics Express 25, 27506-27518 (2017).
1.4 Quantifying laser radiation 23
Ep
Pp = (3)
τp
Peak power of ultrashort laser pulses can vary in the range from MW
to PW, however, we should have in mind that such extreme power is
achieved only for a very short moment of time (duration of the pulse).
Figure 16: Spectrum of the ultrashort laser pulse and its main characteristics:
central frequency (ωc ) and spectral width at FWHM level (∆ωF W HM ).
• Class 2: Only continuous wave (CW) lasers emitting visible (400 nm–
700 nm) light whose power does not exceed 1 mW fall into this class.
Direct exposure to such laser radiation is considered not dangerous if
the exposure time does not exceed 0.25 s (the duration of eye blinking
reflex). Sometimes a subclass is distinguished: class 2M laser system is
a class 2 laser using magnifying optics. It includes visible wavelength
lasers incapable of causing injury in 0.25 seconds unless collecting optics
are used.
• Class 4: It is the most dangerous class. All the lasers whose power
or pulse energy is greater than 3B subclass lasers fall into this class.
All ultrashort pulse (picosecond and femtosecond) lasers also fall into
this class. Both direct and scattered (reflected from matted or diffuse
surfaces e.g., sheet of paper) radiation as well as reflections from op
tical or mechanical element surfaces (lenses, filters, mounts, etc.) is
dangerous to skin and eyes. The most dangerous in that regard are
lasers emitting ultraviolet and infrared radiation since it is not directly
visible to human eye. Maximum safety precautions are necessary when
working with such lasers.
1.5 Laser safety classes 28
Summary of Chapter 1
• Laser operation is based on stimulated emission of radiation. The
term ’laser’ stands for Light Amplification by Stimulated Emis
sion of Radiation.
Figure 18: Diffraction of the light beam and relevant parameters. Red circles
mark an idealized view of photons in the beam.
in that plane have finite precision defined by the uncertainty relations for
photon position and momentum
h h
∆x∆px ≥ ; ∆y∆py ≥ . (8)
π π
These expressions show that the photon propagation direction in one
plane will be undefined as much as is defined its localization in the x-y
plane (the plane perpendicular to propagation direction). For a collection
of propagating photons (the light beam) with a radius w0 = ∆x (Fig.18) the
propagation direction uncertainty for photons in the beam is
h∆⃗k
∆⃗p ≥ . (9)
2π
For wavevector uncertainty of ∆⃗k = ⃗kθ, where θ is the beam divergence
angle, the uncertainty relation for the beam radius and divergence angle reads
as
λ
w0 θ ≥ . (10)
π
This expression can be interpreted as follows: the smaller is the beam
radius (w0 ), the greater is its divergence angle (θ) and vice versa. For ex
ample, taking the laser wavelength of 500 nm and the beam radius of 100
µm then, according to expression in Eq.(10), the divergence angle of such a
λ
beam cannot be smaller than πw 0
≈ 1.6 mrad. This is a fundamental lim
itation following from the aforementioned quantum mechanical uncertainty
principle.
Two boundary cases can be pointed out. Firstly, if we can accurately
define photon propagation direction, then we cannot define where in the x-
y plane the photon is localized. In this case the light is described as an
infinite plane wave which does not diffract. Secondly, if we can accurately
define localization of the photon, then its propagation direction is completely
undefined. In this case the light is described as a spherical wave emitted from
a point source. The real light beams represent an intermediate case between
these two boundary cases. Eq.(10) describes the minimum beam divergence
angle for a given beam radius, such a beam is called diffraction-limited.
Beam with an ideal Gaussian intensity distribution (Fig. 22) are diffraction-
limited beams.
The time-energy uncertainty is expressed as
h
∆E∆τ ≥ (11)
2π
2.1 Photons and their properties 32
1
∆ν∆τ ≥ . (12)
2π
In the case of emission, ∆τ is directly related to the lifetime of an excited
energy level of an atom or molecule: the shorter is the lifetime of the excited
level, the greater is spread in the energy (frequency uncertainty ∆ν) of the
emitted photons or, in other words, the greater is spectral broadening of
the energy level. For example, energy levels that have long lifetime (the so-
called metastable levels) represent an opposite case and emit radiation with a
narrow spectrum. The above considerations may be applied to light pulses as
well; then parameters in expression ∆ν and ∆τ (Eq.(12)) have meanings of
the spectral bandwidth and pulse duration, respectively. Then the expression
in Eq.(12) could be rewritten in a slightly different form:
∆ν∆τ ≥ C , (13)
where C is a certain constant related to the shape of the pulse envelope (Fig.
14). For example, in the case of a Gaussian-shaped pulse C = 0.441, for
sech2 shaped pulses C = 0.315, for rectangular pulses C = 0.886, etc (Fig.
19).
Figure 19: Examples of different pulse shapes. τ0.5 marks pulse duration at
FWHM level. Note that some pulse shape models (such as rectangular) are purely
theoretical models and in practice only approximations of such shapes are possible
to achieve.
transform-limited. This also means that the shorter is the pulse, the
broader is its frequency spectrum and vice versa. For example, continu
ous wave (infinite pulse duration) lasers are essentially monochromatic (their
emission spectrum corresponds to their intrinsic spectral bandwidth), while
femtosecond lasers are polychromatic, they generate a broadband frequency
spectrum. For example, in the near-IR range (λ=1 µm) the spectral band
width of a ∆τ =1 ps pulse (FWHM level) Gaussian pulse is:
0.44
∆υ = = 4.4 · 1011 (Hz) (14)
∆τ
or
∆υλ2
∆λ = = 1.47(nm), (15)
c
whereas the spectral bandwidth of the same central wavelength ∆τ =10 fs
Gaussian pulse is ∆λ=147 nm at FWHM. In the visible range the spectrum
of a latter pulse would essentially be white light. It is important to mention
that artificial narrowing (extension) of the bandwidth-limited pulse spectrum
results in an increase (decrease) of its duration. If the left hand side of
expression in Eq.(13) is significantly greater than the right hand side, the
light pulse is long, but its frequency spectrum is broad. Such pulses are called
phase modulated (or chirped). Manipulation of the pulse phase and
frequency modulation and its control is the basis of temporal compression and
stretching, which is widely exploited in modern laser amplifiers. Compared
with other light sources, lasers are unique in their ability to generate light
pulses and beams with the smallest possible uncertainty: diffraction-limited
beams and bandwidth-limited pulses.
⃗ = 0,
∇·B (19)
→ → →
where H is the magnetic field strength, j is the electric current density, D
→
is the electric field flux density (electric induction), B is the magnetic field
flux density (magnetic induction), ρ is the electric charge density. × is the
curl operator
p indicating
d vector product, while · indicates scalar product and
∂ ∂ ∂
∇ = ∂x , ∂y , ∂z is the nabla operator (expressed in Cartesian coordinate
system).
These equations, named after the physicist and mathematician James
Clerk Maxwell, who published them in early forms in the 1860s, provide a
mathematical model unifying the phenomena of electricity and magnetism
and put forward the foundation of classical electrodynamics. Eq.(16) essen
tially states the Ampere’s circuital law: electric current and time-varying
electric fields induce magnetic fields curling around them. Eq.(17) represents
the Faraday’s law of electromagnetic induction: time-varying magnetic
fields induce electric fields curling around them. The third equation (Eq.(18))
is essentially the Gauss’ law: electric charges create electric fields. The last
equation (Eq.(19)) can be interpreted as the Gauss’s law for magnetism:
magnetic fields do not diverge, they curl around. In other words, Eq.(19)
indicates that magnetic charges do not exist.
For propagation of electromagnetic waves in free space (vacuum) a more
practical equation can be derived from Maxwell equations. This equation
describes light propagation in free space and is called Helmholtz equation.
In one-dimensional case it is expressed as:
∂ 2E 1 ∂ 2E
− 2 2 =0 (20)
∂z 2 c ∂t
where E is the wave amplitude. Solutions of this equation are certain func
tions E(z, t) = f (z ± ct), which are called fundamental waves. Helmholtz
equation is also valid for a harmonic wave:
formula:
eiα = cos α + i sin α. (22)
Using it and a derived expression:
1
sin α = [exp (iα) − exp (−iα), (23)
2i
it is convenient to express the harmonic wave in the complex form:
1
E(z, t) = E0 exp [i(kz − ωt)] + c.c., (24)
2
where E0 is the complex amplitude of the wave and c. c. marks the complex
conjugate part of the equation: if we analyze wave propagation in free space,
the complex conjugate part of expression in Eq.(24) can be ignored.
There are several approximations of the light waves that are important
when analyzing propagation of a laser beam. Firstly, for the sake of sim
plicity, let us assume that light waves are infinite in time, i.e. they have
only a single frequency component (are monochromatic). Plane waves are
waves whose phase at fixed moments of time is stationary in the entire plane
perpendicular to the propagation direction (Fig. 20 a).
(a) (b)
It is written as
E(z, t) = E0 exp [i(⃗k⃗r ± ωt)], (25)
where ⃗r = (x, y, z), ⃗k = (kx , ky , kz ) and ⃗k⃗r = const. Such wave repeats itself
in the ⃗k direction after a distance equal to λ. This wave is infinite in space
and its k vectors are unidirectional, so such wave does not diffract. In reality
plane waves do not exist, but the concept of plane waves is often applied as
a convenient approximation when diffraction of light can be neglected.
2.2 Light waves 36
Figure 21: Principal view of spherical wave wavefront evolution during propaga
tion in free space: far from the source the wavefront can be approximated by a
paraboloidal function and as it propagates further the wavefront becomes very
similar to plane.
Bearing in mind the above, it is clear that plane and spherical waves
represent the two opposite cases considering angular and spatial energy dis
tributions. A plane wave does not have angular spread and its energy is
evenly distributed through the entire space, whereas a spherical wave comes
from a point source and spreads at every direction with its amplitude rapidly
decreasing at a given direction. The light beams represent an intermediate
case between plane and spherical waves. They possess features of both, plane
and spherical waves. The wave nature of light rules out the possibility for
light to propagate in free space without angular spread, but in the case of
light beams, the energy can be sufficiently well (as much as the uncertainty
principle allows) localized along the propagation direction.
The light waves whose wavefront normal vectors make a small angle with
the propagation axis are called paraxial. Such waves are also solutions of the
Helmholtz equation. One of the most important approximations of paraxial
waves is the Gaussian beam – the light beam whose spatial intensity dis
tribution is expressed by a Gaussian function (Fig.22).
2.3 Gaussian beams 37
2
w0 2r2
I(r, z) = I0 exp − , (27)
w(z) w2 (z)
where I0 = |E02 |, E0 is the amplitude of electric
l field, w0 is the minimum
beam radius at the waist (see Fig.23), r = x + y 2 . For every propagation
2
I0
I= 2 , (28)
1 + zz2
R
πw02
zR = . (29)
λ
√
It defines the distance at which beam radius increases by 2 times. The
minimum Gaussian beam wavefront radius of curvature is at z = ±zR , thus
the Rayleigh length is often called the diffraction
l length. On the other
hand, the beam radius can be expressed as w0 = λzR /π.
The entire power of a Gaussian beam can be determined by integrating
the following expression:
∞
1
P = I(r, z)2πrdr = I0 πw02 . (30)
0 2
The beam radius can be defined at different intensity levels (see Fig.12).
Usually, the beam radius is defined at 1/e2 ≈ 0.135 intensity level, as shown
in Fig.22. Eq.(30) is valid for such beam radius definition. For the electric
field amplitude, this would correspond to 1/e ≈ 0.368 level. If the integration
limit in Eq.(30) is changed from ∞ to w(z), such definition of w(z) corres
ponds to concentration of 86% of the whole beam power, while 99% of the
beam power fits into 1.5w(z) intensity level. In general, an ideal Gaussian
beam is infinite, but only a negligible part of its power is contained in the
far periphery. In practice, Gaussian beam radius can also be defined at a
half maximum intensity level (HWHM – half width at a half maximum) and
the beam diameter is determined at the same intensity level (FWHM – full
width at a half maximum). Beam radii definitions at FWHM level and at
1/e2 level are related as:
wF2 W HM
exp −2 2
= 0.5. (31)
w1/e 2
√
which for the diameter of the beam is 2wF W HM = 2 ln 2w1/e2 , thus a factor
of 2 ln 2 is necessary in the Gaussian beam power expression if the beam
radius is defined at FWHM level.
The radius (w(z)) and wavefront curvature (R(z)) of a Gaussian beam
change according to expressions (Fig.23):
z2 1/2
w(z) = w0 1 + , (33)
zR2
and
zR2
R(z) = z 1 + , (34)
z2
Figure 23: Gaussian beam parameters and their evolution during propagation.
Let us determine how radius of a Gaussian beam changes far from the
beam waist. When z ≫ zR in Eq.(33) ”1” can be omitted and the beam
radius increases linearly with z as depicted in Fig.24. This linear dependence
can be expressed as
w0
w(z) ≈ z = θz (35)
zR
where θ is the divergence angle:
w0 λ
θ= = . (36)
zR πw0
2.3 Gaussian beams 40
kr2
ϕ(r, z) = kz − ξ(z) + . (37)
2R(z)
When considering the phase of the beam we usually use the concept of
wavefront. The wavefront is a surface of equal phase with respect to the
beam propagation direction. The wavefront of a plane wave is flat (the plane
perpendicular to the propagation direction) and the wavefront of a spher
ical wave is spherical. As follows from expression in Eq.(37), the Gaussian
beam wavefront is neither flat, nor spherical. The phase of the beam on the
propagation axis (r = 0) is expressed as ϕ(0, z) = kz − ξ(z). The first term
describes a plane wave. So at the waist (z = 0) the Gaussian beam is a plane
wave, i.e. its wavefront is flat and radius of curvature R(z) is infinite. The
2.3 Gaussian beams 41
second term describes the phase shift of the on-axis beam part compared to
either plane or the spherical wave. When z = ±zR , the phase shift is ±π/4
and approaches ±π/2 when z → ±∞. This means that during propagation
of the Gaussian beam through infinitely long distance the maximum accumu
lated axial phase shift is π. This is called the Gouy effect. The third term in
Eq.(37) describes the phase shift for off-axis points compared with the phase
of the on-axis points. In other words, it shows the degree of wavefront bend
ing. It can be shown that a surface of equal phase is parabolic with a radius
of curvature R. The evolution of the radius of curvature for a Gaussian beam
during propagation is described by Eq.(34), and is graphically illustrated in
Fig.25. It is clear that the minimum radius of curvature of a Gaussian beam
is at z = ±zR : R = ±2zR , i.e., the wavefront here is maximally curved. The
radius of curvature both sides from zR is only increasing with propagation.
Figure 26: Comparison of wavefronts:(a) plane waves, (b) spherical waves and (c)
Gaussian beam.
2.4 Coherence
Coherence is a distinctive feature of laser radiation. In general, the term
coherence is used to describe phase correlation between two monochro
matic waves. Waves with correlated phases are called coherent and waves
with uncorrelated (random) phases are called incoherent. Two types of
coherence can be distinguished:
Figure 27: Principal illustration of temporal and spatial coherence of light waves.
tc and lt mark temporal coherence time and spatial coherence length, respectively.
In the spatial coherence part interference intensity pattern is depicted in the cases
of perfectly coherent, partially coherent and incoherent light source.
Figure 29: A sequence of harmonic waves with random length and phase.
c λ2
lc = ≈ , (38)
∆ν ∆λ
which essentially reflects the boundary case of the uncertainty principle; here
∆λ is the natural spectral linewidth of a given light source. Having this in
mind, we can now explain what is depicted in Fig.28. When optical path
difference between interferometer arms (d) approaches the temporal coher
ence length (lc ), the contrast of interference fringes drops and when d exceeds
temporal coherence length, the interference pattern disappears. This simply
reflects the fact that there are no ideal monochromatic light sources, there
fore any real light source has a finite coherence time and temporal coherence
length.
The temporal coherence of several representative light sources is com
pared in Table 1. For any CW light source (including CW lasers) the tem
poral coherence in this case describes the monochromaticity of such light
2.5 Spatial coherence 45
rλ λ
lt = = , (39)
s θ
where r is the distance from the light source, s is the distance between the
points of the source and θ is the angle at which these points are visible.
According to this definition, an ideal plane waves and point light sources
are fully spatially coherent. Any other incoherent light source can be made
coherent by turning it into a point light source, however, this is possible only
at the cost of its radiation power..
Let us discuss the meaning of spatial coherence in the case of a Gaussian
beam. Comparing expressions in Eq.(36) and Eq.(39), it is clear that they
match with precision of a constant, so the radius of a Gaussian beam waist
corresponds to the length of spatial coherence.
To characterize a real Gaussian beam, another parameter called M 2 is
also introduced. This parameter indicates the difference between the diver
gence angles of real and ideal beams with the same beam waist. Moreover,
this parameter is often used as a beam ”quality” parameter. The term ”qual
ity” is more often defined as a measure for how well the beam can be focused
(Fig.31) and it is quantified as either M 2 parameter:
πw0 θ
M2 = , (40)
λ
or beam parameter product (BPP):
BP P = w0 · θ . (41)
When we calculate diameter of a focused beam (df ) we have to take into
account the M 2 parameter:
4λf
df = M 2 , (42)
πd0
2.5 Spatial coherence 47
Figure 31: Gaussian beam size evolution when it is focused with a lens. The
red contours depict laser beam intensity distribution and d1/e2 denotes laser beam
diameter at 1/e2 of maximum intensity level.
Using the light interference terminology, the Gaussian beam can be viewed
as a result of interfering multiple plane waves propagating at different dir
ections. Since these waves are phased in a certain way, due to interference
they create a specific spatial modulation which is observed as the light beam.
Thus, laser light sources are also unique with regard to spatial coherence as
laser beams are spatially coherent and carry all the power emitted by the
laser source.
Concluding the coherence issues, the concepts of temporal and spatial
coherence could be merged by defining the volume (three-dimensional) co
herence:
lv = lc lt2 , (43)
which describes the spatio-temporal coherence and directly indicates the
volume where the energy of light is localized. The concept of volume (spatio
temporal) coherence applies to pulsed lasers.
2.5 Spatial coherence 48
Summary of Chapter 2
• Light is described either as electromagnetic waves or as particles
(photons). Each concept attributes to different properties of laser
radiation.
3 Optical resonators
An optical resonator (optical cavity) is an extremely important component
of any laser. The simplest optical resonator consists of two reflecting mirrors,
which play several relevant roles in laser operation:
In this chapter we will discuss the role of resonators in lasers: light propaga
tion analysis in optical resonators, resonator stability criteria, transverse and
longitudinal modes and energy loss mechanisms in optical resonators.
optical axis, as depicted in Fig.32. These two parameters define the ray vec
tor (x, θ). An effect of any optical element is equivalent to multiplication of
the ray matrix defining the optical element by ray vector, resulting in new
values of the parameters characterizing the ray:
[ [ [
xout A B xin
= (44)
θout C D θin
where indices in and out mark vector parameters of the initial (incident) and
final (outgoing) rays respectively. We will now discuss the principle how ray
matrices for various optical elements are constructed and how their A, B, C
and D values are determined. Assuming that inclination angle of a ray to the
optical axis is small (which is essentially the condition of paraxial approxim
ation as was already stated), the ray parameter changes are described by a
two equation system:
∂xout ∂xout
xout = xin + θin
∂xin ∂θin
∂θout ∂θout
θout = xin + θin . (45)
∂xin ∂θin
These equations are equivalent to expression in Eq.(44), so matrix ele
ments can be expressed as partial derivatives where
∂xout ∂xout ∂θout ∂θout
A= , B= , C= , D= . (46)
∂xin ∂θin ∂xin ∂θin
Matrix elements A and D can be understood as spatial and angular mag
nification, respectively. With accordance to these expressions, let us write
ABCD matrices for several optical media and elements which are important
for optical resonators. To start with, let us determine the ABCD matrix for
3.1 Light rays and ray (or ABCD) matrices 51
Figure 34: Propagation of a light ray through the intersection of two media.
3.1 Light rays and ray (or ABCD) matrices 52
At the intersection of two media, a distance of the light ray from the
optical axis remains unchanged: xout = xin , whereas the change of the in
clination angle can be calculated from the Snell’s law: θout = nn12 θin . Con
sequently, the ray matrix for the intersection of two media with different
refractive indices is
[
1 0
MI = . (50)
0 n1 /n2
In analogy, the ray matrix for the intersection of two optical media with
different refractive indices and curved (with R radius of curvature) intersec
tion surface can be written as:
[
1 0
MCI = n1 /n2 −1 . (51)
R
n1 /n2
In that case, if we assume a lens as two curved surfaces, the ray matrix
for a lens can be derived as a product of two ray matrices of curved surfaces:
[
1 0
ML = MCI1 MCI2 = . (52)
− f1 1
Here we also assumed that the refractive index of the surrounding medium
is n1 = 1 and the refractive index of a lens material is n2 = n. Then f is the
focal distance of the lens and is expressed as:
1 1 1
= (n − 1) − , (53)
f R2 R1
where R1 and R2 are the radii of curvature of the lens’ surfaces. This expres
sion is also known as the Lens Maker’s formula. Notice that if f > 0 the
lens is convex and if f < 0 it is concave. Based on a similar consideration, we
can determine ABCD matrix for a curved mirror, whose center of curvature
is located on the optical axis:
[
1 0
MCM = , (54)
− R2 1
where R is the mirror’s radius of curvature. Comparing the expressions for
ML and MCM it is clear that a curved mirror focuses rays as a lens with a
focal length of fv = R/2. Therefore a mirror with R > 0 is regarded as
concave and a mirror with R < 0 as convex. In the case of a plane mirror,
the same ray matrix can be used by setting R = ∞, i.e.:
[
1 0
MPM = . (55)
0 1
3.2 Laser resonator stability 53
Concluding the results of this section, we can note that an ABCD matrix
can be written for any optical element no matter how complex it is. This
means that light ray propagation can be modelled in sophisticated optical
systems consisting of various optical elements whose ABCD matrices are
known. Such method is called ray tracing. Moreover, ABCD matrix form
alism can be directly applied to a Gaussian beam assuming that the optical
element transforms its complex propagation parameter in the following way:
Aqin + B
qout = , (56)
Cqin + D
where qin and qout are incident and passing (through the optical element)
complex propagation parameters of the Gaussian beam, respectively. The
complex propagation parameter is expressed via two real functions: beam
radius w(z) and wavefront curvature R(z), which fully ascribe the Gaussian
beam:
1 1 λ
= −i 2 . (57)
q(z) R(z) πw (z)
From this we can determine how relevant parameters (i.e., w(z) and R(z))
of the Gaussian beam change during beam propagation in various optical
systems.
Figure 35: Laser beam in an empty optical resonator consisting of two curved
mirrors.
It is obvious that in order for the laser to work, light should not escape
the resonator: after each reflection the beam must maintain a stable size fit
ting into the resonator mirror aperture which is finite. Light will be confined
in such optical system only if it fully reproduces its spatial distribution after
each round-trip between the resonator mirrors. This is called beam conver
gence condition, which implies that only the laser beams with a specific
spatial distributions of electric field (intensity) are able to do so. These are
called spatial (transverse) modes of the resonator. On the other hand,
beam convergence condition poses strict requirements on the resonator para
meters and geometry: mirrors radii of curvature R1 and R2 and resonator
length l. By applying ABCD law for a Gaussian beam that fully reproduces
itself after each resonator round-trip (qin = qout ):
Aq + B
q= , (58)
Cq + D
where A, B, C, D values are derived from Mresonator = MCM 2 MS MCM 1 MS
and skipping mathematics for the sake of simplicity we can derive laser
resonator stability condition:
l l
0≤ 1− 1− ≤ 1. (59)
R1 R2
Here R1 and R2 are resoantor mirros radii of curvature and l is the length of
a resonator. gi = 1 − Rli are called g parameters of the particular resonator
mirror. This inequality is depicted graphically in Fig.36 which is called res
onator stability diagram. Blue area in the diagram shows the space of g
parameter values, where resonator is stable.
3.2 Laser resonator stability 55
Figure 36: Resonator stability diagram, where stability zone is shown by blue
area, and several resonator types corresponding to different points in the stability
diagram. Image adapted from Wikimedia under CC BY-SA 3.0 license.
resonator geometry.
Therefore, summarizing both facts, we can claim that Gaussian beam is
a spatial mode of a spherical mirror resonator. Generally, any resonator can
be regarded as an approximation of a spherical resonator, so the Gaussian
beam is a spatial mode of any optical resonator, and consequently, all lasers
produce radiation in the form of a Gaussian beam. Furthermore, the
M 2 factor introduced in the previous chapter quantitatively indicates how
much a beam generated in a certain resonator differs from an ideal Gaussian
beam generated in a spherical mirror resonator.
The current discussion may raise a question whether a laser with an un
stable resonator can operate at all? It is possible indeed, however, lasers
operating with unstable resonators typically have high gain which balances
large diffraction losses. In an unstable resonator light will be ejected after a
certain number of round-trips and output coupling is such lasers is also differ
ent: light leakage due to diffraction is often taken as the useful laser output
(Fig,38). Unstable resonators also have transverse modes, but their intensity
distributions are very complex and beyond the scope of this textbook. In
terestingly, unstable resonators are less sensitive to alignment instabilities.
3.3 Hermite-Gaussian beams and Laguerre Gaussian beams 57
where ξ(z) = tan−1 (z/zR ) and R(z) is the wavefront radius of curvature, E0 is
the amplitude of electric field, w0 is the minimum beam radius (see Fig.23).
The first line in this expression describes the amplitude evolution during
propagation in the z direction, whereas the second line in the expression
describes the phase evolution.
Lasers with very stable resonators may produce beams with more complex
electric field (intensity) distributions, which reproduce themselves after each
round-trip, and satisfy the convergence condition. These beams are called
higher order spatial (transverse) resonator modes.
Two cases of higher order modes can be distinguished. In stable optical
resonators with no cylindrical symmetry higher order transverse modes
can be ascribed as Hermite-Gaussian beams, whose their electric field
amplitude is expressed as
3.3 Hermite-Gaussian beams and Laguerre Gaussian beams 58
�√ � �√ � [ 2
2x 2y w0 x + y2
E(x, y, z) = E0 Hm Hn × exp − 2 × (61)
w(z) w(z) w(z) w (z)
k(x2 + y 2 )
[
exp [−i (kz − (1 + m + n)ξ(z))] × exp −i .
2R(z)
where quantities E0 , w0 , w(z), R(z) and ξ(z) are defined in the same manner
as in the case of a Gaussian beam, while functions Hm and Hn are called
Hermite polynomials, whose indices m and n mark the order of the polyno
mial. Hermite polynomials are defined as
H0 (u) = 1 (62)
H1 (u) = u (63)
H2 (u) = 2u2 − 1 etc. (64)
√
wmn = 2m + n + 1 w00
√
θmn = 2m + n + 1 θ00 , (65)
where index 00 stands for the Gaussian beam (TEM00 mode). Higher order
spatial resonator modes, compared to a Gaussian beam, have larger beam
radius and larger angular divergence. This tells us that such laser beam
quality is lower than that of a Gaussian beam and that higher order modes
can only be excited in a very stable optical resonator where diffraction losses
are small. On the other hand, the laser output with such a complex intensity
distribution is rarely used in practice, thus usually higher order mode gener
ation is suppressed: mitigated by specifically inducing losses (e.g., diffraction
losses) to force laser resonator generate the TEM00 mode.
In stable optical resonators with cylindrical symmetry higher order
transverse modes can be ascribed as Laguerre-Gaussian beams and their
3.3 Hermite-Gaussian beams and Laguerre Gaussian beams 60
where ρ, φ and z are the coordinates (in cylindrical coordinate system) and
Llm is the generalized Laguerre polynomial with indices m and l. These
polynomials are the solutions of the Laguerre differential equation. They are
defined via recurrence relation with the first ones being:
Its phase front has a helical trajectory with an undefined phase and in
tensity minimum at the center of the beam (Fig.41).
Figure 41: Intensity distribution (left) and wavefront (right) of an optical vortex.
Dashed black line marks phase front rotation pattern.
∆φ = qπ , (70)
Figure 42: Standing waves in plane parallel resonator and the corresponding
longitudinal mode index q.
z2 z1
kq l − (m + n + 1) arctan − arctan = qπ. (73)
zR zR
Setting m and n constant, the former expression can be rewritten as
π
kq+1 − kq = , (74)
l
3.4 Resonance frequencies 63
and expressing k = 2πνn0 /c, where n0 is the refractive index of the medium
filling the resonator, and transferring into frequency domain we get
c
νq+1 − νq = . (75)
2n0 l
This frequency difference is called the intermodal distance (Fig.43),
which is inversely proportional to resonator length. Assuming that m and
n are equal to 0 (the case of a Gaussian beam), the intermodal distance is
independent of resonator configuration.
Figure 43: Intermodal distance in the case of a TEM00 mode (Gaussian beam).
Figure 44: Example of multimode He-Ne laser spectral structure with 500 MHz
intermodal distance.
the complex amplitude of the incident wave, while B1 , B2 , ... and A1 , A2 , ...
as complex amplitudes of reflected and transmitted waves, respectively, as
shown in Fig.45.
A1 = tt′ Ai
A2 = tt′ r′2 Ai eiδ
A3 = tt′ r′4 Ai e2iδ ... (79)
3.5 Fabry-Perot resonator 66
etc, and the full transmitted wave complex amplitude is At = A1 +A2 +A3 +...:
√
1 − eiδ R
Ar = Ai
1 − Reiδ
T
At = Ai . (81)
1 − Reiδ
Let us rewrite the latter expressions as wave intensity ratio assuming that
I = AA∗ . The ratio of reflected and incident wave intensities then is
Ir Ar A∗r 4R sin2 2δ
= = . (82)
Ii Ai A∗i (1 − R)2 + 4R sin2 2δ
In a similar manner, the ratio of transmitted and incident wave intensities
is:
It At A∗t (1 − R)2
= = . (83)
Ii Ai A∗i (1 − R)2 + 4R sin2 2δ
Since there are no losses in the resonator, it is easy to determine that It +
Ir = Ii . From expression in Eq.(83) it follows that a maximum transmittance
of the Fabry-Perot resonator is when δ = 2πm, where m is any integer
number. Then It /Ii = 1 and does not depend on the reflection coefficient
R. On the other hand, we can notice that maximum Fabry-Perot resonator
transmittance coincides with a standing wave condition (Eq.70) and intrinsic
frequencies can be expressed as
c
νm = m , (84)
2n0 l cos θ
and the intermodal distance is
c
∆νm = νm+1 − νm = . (85)
2n0 l cos θ
3.5 Fabry-Perot resonator 67
In the case of large reflection coefficient, the transmittance maxima are sharp
and minima are deep. When mirror reflectance decreases, the maxima spread
and the amplitudes of minima increase. In that regard it is convenient to
define the spectral width of a longitudinal mode which depends on the res
onator length and mirror reflectance coefficient:
c 1−R
δν = √ . (86)
2n0 l cos θ π R
Fabry-Perot etalon can be also used as a spectral device: the intermodal
distance can be varied by adjusting l and the spectral linewidth can be varied
by changing R. In practice, Fabry-Perot etalon is inserted into a resonator of
a CW laser when single longitudinal mode operation and very narrow spec
tral bandwidth is desired. Fig.47 shows an example of He-Ne laser spectral
structure and a Fabry-Perot etalon modal structure. The infinite set of res
onator longitudinal modes is limited by the width of the lasing medium gain
curve and the Fabry-Perot etalon transmittance: without the Fabry-Perot
etalon there would be a number of longitudinal modes in the emitted He-Ne
3.6 Losses in resonators 68
laser spectrum, but the Fabry-Perot etalon filters most of them. Leaving just
a single longitudinal mode that falls under the gain curve.
Figure 47: He-Ne laser spectral structure and Fabry-Petor etalon mode structure.
The etalon filters longitudinal modes effectively forcing He-Ne laser to emit only
a single longitudinal mode.
2n0 l/c
tp = . (89)
1 − r 1 r2
The former expression can be interpreted in the following way: the photon
lifetime in a resonator equals to a resonator round-trip time divided by
round-trip losses. On the other hand, photon lifetime can be related with
characteristic resonator mode extinction time:
dE E
=− , (90)
dt tp
where E is the energy of a resonator mode. Resonator quality is usually
defined as a ratio of accumulated energy and dissipated power:
ωE ωE
Q= =− , (91)
P dE/dt
where ω is the frequency of radiation. Comparing both expressions we get
that Q = ωtp . Now let us discuss the causes of optical losses in resonators.
These are produced by several mechanisms:
• Diffraction losses play a double role. On the one hand, when diffrac
tion losses can be decreased or increased by selecting a certain resonator
configuration: mirror radius of curvature and resonator length. Phys
ical apertures of elements inside the resonator are often smaller than
apertures of the resonator mirrors, so these diffraction losses due to
the former must be taken into account. On the other hand, diffrac
tion losses help to suppress higher order TEM modes whose generation
is usually not desirable. Finally, diffraction losses, as a rule, are the
largest in long resonators, but it is the diffraction-related spatial filtra
tion of transverse modes which enables formation of light beams with
very high spatial quality (in simple words, all beam irregularities due to
not ideal optical surfaces, dust, etc are filtered out due to diffraction).
Lasers with short resonators and low diffraction losses (e.g., semicon
ductor lasers and laser diodes) generate low spatial quality beams, their
3.6 Losses in resonators 71
Summary of Chapter 3
• Optical resonator is an indispensable component of any laser,
that is responsible for energy accumulation, feedback and forma
tion of spatial and spectral characteristics of the laser radiation.
Eph = hν = E2 − E1 . (92)
4.1 Einstein’s treatment of spontaneous and induced transitions 73
on the physical state of the matter (e.g., in rarified gases metastable levels
can exist for tens of seconds and even longer). Note that tspont is a material-
related characteristic. In general, lifetime of an energy level depends on the
frequency of radiated photon: the greater is the frequency (also meaning
greater energy difference between the levels), the shorter is the correspond
ing lifetime. This is a logical result since electrons in atoms always try to
occupy the state with a lowest energy.
According to the uncertainty principle, the real energy levels are broadened.
This means that these energy levels emit not infinitely narrow (monochro
matic) light, but a spectrum with finite bandwidth, which is characterized
by a certain linewidth (the so-called natural linewidth) which is described by
a spectral function g(ν):
∞
g(ν)dν = 1. (95)
0
The quantity g(ν)dν indicates the probability that the frequency of emitted
photon is in the frequency interval ν + dν, and the spectral intensity of
radiation is defined as I(ν) = I0 g(ν).
Einstein discovered that besides the spontaneous emission there also ex
ists another form of emission of radiation, which was called the stimulated
emission. This means that excited atoms can return to the ground (un
excited) state not only spontaneously but also can be forced to do that.
Consider an excited atom (or an ensemble of atoms) which is in the field of
external radiation that has spectral energy density ρ(ν). If there are frequen
cies in the external field that match the frequency of transition from excited
to the ground energy level, stimulated (induced) transitions occur, as shown
in Fig.49. An important property of stimulated emission is that emitted
photon has the same direction and energy as the incident photon
of the external radiation field, and the corresponding electromagnetic
wave has the same polarization and phase as the incident one.
′
The rate of stimulated transitions can be expressed as B21 (ν)N2 , where
′
B21 (ν) is a specific function describing the probability of such transition,
′
and is expressed as B21 (ν) = B21 g(ν). Here B21 denotes the Eintein’s
coefficient for stimulated emission. In this case, the evolution of excited
energy level population can be expressed by the following equation
∞
dN2
= −B21 N2 g(ν)ρ(ν)dν = −N2 B21 ρ(ν), (96)
dt 0
for simplicity, assume that ρ(ν) is a narrow line and ρ(ν) = ρ(ν0 ). Quantity
B21 ρ(ν) is called rate of stimulated emission.
Now let us consider absorption. Since absorption is a process which occurs
only in the presence of external field, it is a stimulated (induced) process,
which is depicted schematically in Fig.50.
8πhν 3 1 8πν 2 1
ρ(ν) = hν = × hν × hν . (98)
c3 e kT − 1 c 3
e kT − 1
Here the first multiplier denotes the number of modes within a fre
quency interval, the second denotes the photon energy and the third
denotes the modal population factor, i.e., number of photons per mode;
N2 hν
= exp − ; (99)
N1 kT
• There is a thermodynamic equilibrium between the atom system and
radiation, which means that the population of energy levels does not
change when temperature is stable.
dN2
= N1 B12 ρ(ν) − N2 B21 ρ(ν) − N2 A21 = 0. (101)
dt
Inserting expression of ρ(ν) we get:
� � � �
8πhν 3 8πhν 3
N2 B21 3 hν/kT + A21 = N1 B12 3 hν/kT . (102)
c e −1 c e −1
8πhν 3 A21
= hν/kT
. (103)
3
c ehν/kT −1 B12 e − B21
In this way we get one equation with three unknowns. It is satisfied when:
A21 8πhν 3
B12 = B21 , = . (104)
B21 c3
These are Einstein relations. Let us discuss their meaning. From ex
pressions in Eq.(104) it is clear that Einstein’s coefficients A21 , B21 and B12
are mutually related: if we know at least one of them, we can find the other
two. Another important thing to note is that coefficients for stimulated
emission and stimulated absorption are equal (in the case of non-degenerate
energy levels). This implies that stimulated emission and stimulated absorp
tion are the opposite processes, and not spontaneous emission and stimulated
absorption, as might seem at the first sight. If energy level population is
N1 > N2 , then external radiation field is absorbed and attenuated, while if
N2 > N1 , photons that appear as a result of stimulated emission add to the
external radiation field which is thus amplified. A situation when (N2 > N1 )
is called population inversion. It is clear that in any laser the condition
of population inversion must be satisfied.
Now let us discuss the role of spontaneous emission. The rate of stim
ulated emission depends on the energy density of external radiation field:
W21 = B21 ρ(ν), while the rate of spontaneous emission does not. So if we
ignore transitions related to spontaneous emission when describing the ther
modynamic equilibrium between two energy levels, we get N2 B21 = N1 B12
which, knowing that B21 = B12 , would yield N2 = N1 ; this is obviously
impossible in the case of thermodynamic equilibrium.
An important aspect is the ratio between the rates of spontaneous and
stimulated emissions, which can be defined as
4.2 Gain coefficient 78
A21 hν
R= ≃ exp −1 . (105)
B21 ρ(ν) kT
Let us compare this ratio in the case of microwave and optical range at
equal ambient temperature, T = 300 K. Microwave frequency is roughly
ν ≈ 1010 Hz, then hν/kT ≈ 1.6 × 10−3 and R ≈ 0.0016. This indicates that
in the microwave range the rate of stimulated emission significantly exceeds
the rate of spontaneous emission. The opposite result is achieved in the
optical range, assuming ν ≈ 1015 Hz. Then we get R ≈ 160. This means
that the rate of the stimulated emission in the optical range is very low
compared to that of the spontaneous emission, therefore, it is much harder
to achieve laser amplification in this spectral range. With this in mind, it is
easy to understand why masers were invented much earlier than lasers. It is
also obvious that laser generation and amplification is much harder to realize
in the blue and ultraviolet spectral range. Of course, emission rate ratio for
metastable energy levels may be substantially different in the optical range
but the general tendency is the same.
Although Einstein’s relations were derived in the case of thermodynamic
equilibrium, they are also valid in diverse conditions as transition probabilit
ies are specific characteristics for a given material. In reality, a working laser
is hardly in the thermodynamic equilibrium, but despite this, the meaning
of Einstein relations remains the same.
dIν
= N2 B21 g(ν)ρν − N1 B12 g(ν)ρν hν. (106)
dz
Using Einstein relations Eq.(104) and definition of intensity we get
dIν c2 A21
= N2 − N1 g(ν)Iν . (107)
dz 8πν 2
The solution of this equation is
where Iν (0) is the initial intensity of radiation and γ(ν) is the gain coeffi
cient, which is expressed as:
c2 A21 λ2
γ(ν) = N2 − N1 g(ν) = ∆N g(ν) . (109)
8πν 2 8πtspont
The amplification of light is achieved when γ(ν) > 0, while attenuation
occurs when γ(ν) < 0. The sign of gain coefficient depends only on the differ
ence of energy level population ∆N , since all other coefficients are positive.
If the ensemble of atoms is in a thermodynamic equilibrium, population of
the energy levels obeys Boltzmann‘s distribution, then the following is always
true: N2 < N1 and ∆N < 0 when T > 0. Consequently, in a two-level system
we always have only attenuation (absorption) of the external field and in this
case γ(ν) < 0 is simply equivalent to the absorption coefficient γ(ν) ≡ α(ν).
On the other hand, even if the system is far from thermodynamic equilib
rium and external energy source is present, the population inversion is still
impossible because the rates of stimulated emission and stimulated absorp
tion are equal. In such a two-level system of atoms, amplification is only
possible when T < 0, which has no physical sense according to the classical
(Boltzmann) definition of the absolute temperature. To this end, in the early
studies, the population inversion was related to an effective negative absolute
temperature.
It is clear that in order to achieve amplification of light (laser generation),
a two-level system is insufficient to create the population inversion (∆N > 0).
Leaving the issues how the population inversion is achieved in practice, we
currently only note that population inversion is achieved not only by setting
the system out of thermodynamic equilibrium (with the use of an external
energy source) but also by involving additional energy levels.
4.3 Homogeneous and inhomogeneous spectral linewidth broadening 80
Figure 52: Shape and width of spectral line in the case of homogeneous broadening
of energy levels.
∆ν
gH (ν) = , (110)
2π (ν − ν0 )2 + (∆ν/2)2
where ∆ν is the spectral linewidth at half intensity level as shown in Fig.52
and equals to
A21 1
∆ν = = . (111)
2π 2πtspont
Such spectral line broadening is identical for any atom in the ensemble and
since we cannot distinguish between the radiating atoms, it is called homo
geneous.
Inhomogeneous spectral line broadening is caused by a certain
distribution of features of the individual atoms (and consequently,
their energy levels). In this case we can distinguish between radiating
atoms or their groups, as each atom or a group of atoms emit radiation
with slightly different frequency. There are several physical reasons of in
homogeneous spectral line broadening, which are different for various laser
materials. Here we will discuss several of them which are the most important
and common.
In the case of solid-state lasers, there are two main reasons for inhomo
geneous spectral line broadening. The first is attributed to amorphous struc
ture of laser host material. Various glasses doped with rare-earth metal ions
(Nd, Yb, etc) represent the most illustrative example. Glasses are amorph
ous materials, which from a macroscopic view are isotropic. However, from
4.3 Homogeneous and inhomogeneous spectral linewidth broadening 82
the microscopic view of structure, they are not strictly isotropic but com
posed of many small ordered structures which are ”frozen” in the isotropic
matrix of the glass. Orientation and properties of distinct structures (such
as mechanical stress, concentration of impurities and their energy levels) is
slightly different, which in turn slightly modify the energy levels of ions that
are embedded in the glass matrix. Then the overall radiation spectrum is
a sum of all spectra emitted by these structures which consequently makes
the spectrum relatively broad compared to that of homogeneously broadened
line. For example, the spectral line of Nd ions embedded in a YAG crystal
is homogeneously broadened, whereas the spectral line of the same Nd ions
embedded in a glass is inhomogeneously broadened and its width is roughly
50 times broader.
Inhomogeneous spectral line broadening in crystalline materials occurs
due to various defects and impurities in the crystal lattice, so each radiating
ion is surrounded by a slightly different environment. However, in most cases
such inhomogeneous broadening mechanism is relatively weak, as modern
technologies enable production of laser materials with very high optical and
chemical quality.
Several mechanisms of inhomogeneous spectral line broadening are spe
cifically attributed to gas lasers. The first occurs in gas mixtures due to iso
tope impurities, which emit radiation with slightly shifted frequencies. The
second mechanism is attributed to spectral line broadening due to Doppler
effect. Random velocity distribution of gas atoms or molecules suggests that
a stationary observer due to Doppler effect sees slightly shifted frequency of
radiation, which depends on the direction of the moving particle with re
spect to the observer. Distribution of atom or molecule velocities depends
on the mass and temperature and is described by the Maxwell-Boltzmann
distribution function, which can be converted to the distribution function
of radiation frequency. For an inhomogeneously broadened spectral line the
function is
� �
1/2 2
4 ln 2 1 ν − ν0
gI (ν) = exp −4 ln 2 . (112)
π ∆ν ∆ν
In the ideal case, the shape of the inhomogeneously broadened spectral line
is a Gaussian function whose width is defined as
1/2
8kT ln 2
∆ν = ν0 (113)
M c2
and is a function of atomic mass M and temperature T . In the case of
4.4 Three-level and four-level lasers 83
Figure 53: The shape of the spectral line in the case of inhomogeneous broadening
which consists of multiple homogeneously broadened lines shown by dashed curves.
has to be short, it is also desirable that any radiative transitions from level
(3) do not occur or their rate is very low. In this way, excitations of level
(3) are rapidly transferred by non-radiative transitions to level (2), which
is metastable, i.e. has long lifetime. Levels (2) and (1) are the upper and
lower laser levels, respectively. To achieve population inversion between laser
levels (2) and (1), lifetime of energy level (1) has to be short. i.e. its excit
ations rapidly relax to the ground state (0) by radiative and non-radiative
transitions. Therefore, in a four-level laser, any excited atom in level (3)
quickly transfers its excitation to level (2) and since level (1) is empty (very
quickly relaxing), the population inversion is achieved by very simple means.
Examples of four-level laser systems include Nd-doped solid-state lasers (e.g.
Nd:YAG, Nd:glass), Ti:sapphire, excimer and most of gas lasers.
A different situation occurs in a three-level laser (Fig.55), where the
lower energy level (1) matches the ground level. In this case, to achieve
population inversion it is necessary to excite at least half of the atoms from
the ground level to level (3). Obviously, this can be obtained only using
sufficiently high pumping rate. The most prominent example of three-level
laser is ruby laser.
4.4 Three-level and four-level lasers 85
Figure 57: Energy level structure of Nd ion embedded in YAG crystal, pump and
laser transitions.
Figure 58: Simplified view of excited states of He and Ne atoms and possible
transitions.
There are also some unique laser media, e.g., dye lasers, which possess an
extremely broad emission spectrum. This is related to the energy structure
of dyes. They are usually in liquid phase and their electronic states have
multiple vibrational energy sub-levels. Lasing can occur between any (or all)
of these vibrational energy levels, therefore, lasing is broadband (see Fig.59).
This unique feature is important in femtosecond lasers: the shorter is the
4.6 Pump sources 88
pulse duration, the broader laser gain bandwidth is required to generate such
pulse (see Eq.13). Dye lasers very historically the first femtosecond lasers.
The relevant parameters of several widely used laser media are listed in
Table 2. Gain bandwidth for various laser media differs by more than 6
orders of magnitude. Rhodamine 6G dye and Ti:sapphire crystal are two
unique media with very large gain bandwidth.
Table 2: The main parameters of several laser media. λ0 is the laser wavelength, σ
is the laser transition cross-section, tspont is is the lifetime of the upper laser level,
∆ν is the gain bandwidth, I and H stand for inhomogeneous and homogeneous
line broadening, respectively.
Laser medium λ0 ,µm σ,cm2 tspont ,µs ∆ν Broadening
He-Ne 0.6328 5.8 × 10−13 0.03 1.7 GHz I
Ar+ 0.5145 2.5 × 10−13 0.01 3.5 GHz I
CO2 10.6 3 × 10−18 2.9 × 106 60 MHz I
Nd:YAG 1.064 2.8 × 10−19 230 135 GHz H
Nd:glass 1.054 4.1 × 10−20 350 8 THz I
Ti:sapphire 0.8 3.8 × 10−19 3.9 100 THz H
Rhodamine 6G 0.57 2 × 10−16 0.005 45 THz H/I
(b)
Figure 60: (a): noble gas-filled flashlamps used for optical laser pumping. Image
adapted from Wikimedia under CC BY-SA 2.0 license. (b): solid-state laser res
onator pumping geometries using helical or linear shape flashlamps.
Most laser diodes emit in the near-IR and are of particular importance
for Yb-doped laser pumping, as Yb ion has absorption lines only in
the near-IR. Laser pumping is used only when the absorption lines or
bands of the active laser medium are in the visible spectral range and
excited laser energy levels have short lifetimes. This is the case for the
4.7 Laser oscillation 91
Figure 62: Principal setup of laser oscillator (generator), amplifier and MOPA
system.
1 1
γ0 (ν) ≥ ln =α. (115)
2ll r1 r2
4.7 Laser oscillation 93
Summary of Chapter 4
• The stimulated emission is the fundamental effect that produces
laser radiation.
• Laser generation starts when laser gain exceeds losses per round-
trip in the resonator. This is the laser oscillation condition.
dϕ d∆N
= = 0, (118)
dt dt
we can find the equilibrium values for population and photon flux densities:
1
∆N0 =
Btp
RBtp − 1/tspont
ϕ0 = . (119)
B
From the above expressions it follows that when Rs = (Btp tspont )−1 , ϕ0 =
0. Here Rs marks the threshold pump rate. When it is exceeded, the laser
starts to generate and photon flux density increases. We define a ratio r =
R/Rs indicating how many times the threshold pump rate is exceeded. Then
the second equation in the system of equations (119) can be written as
r−1
ϕ0 = . (120)
Btspont
Let us examine what happens when we have small deviations of popula
tion and photon flux densities from their equilibrium values:
where N ′ (t) ≪ ∆N0 and ϕ′ (t) ≪ ϕ0 are small quantities. Inserting these
expressions into the initial equations and using relations in Eq.(119), we get
a system of equations for small variations of population and photon flux
densities:
dN ′ ϕ′
= −RBtp N ′ −
dt tp
′
dϕ 1
= RBtp − N ′. (122)
dt tspont
We will search for solutions of these equations in the form ϕ′ (t) = eAt sin(ωm t),
i.e., assuming that a perturbation of the photon flux density has the shape of
a damped sine function with a damping parameter (a characteristic time con
stant) A and oscillation frequency ωm . Solving equations (without focusing
on mathematical subtleties) we get:
5.1 Free-running regime and relaxation oscillations 98
r
A= ,
tspont
�
r−1
ωm ≈ . (123)
tp tspont
Figure 64: Relaxation oscillations in a laser and evolution of energy level popula
tion.
quality is very high. Then the relaxation oscillations become chaotic and
emerge as a sequence of non-periodic intensity spikes. The reason for this
behavior is explained as follows. The lowest order mode TEM00 is gener
ated first and strongly depletes population of the laser medium within its
volume. Since the volume of higher order modes is larger, they are easily
excited because the population of laser medium in the peripheral part re
mains large and exceeds the threshold value. Each of the transverse mode
produces a sequence of its own relaxation oscillations, which overlap produ
cing non-periodic intensity spikes. Thereafter the TEM00 mode is generated
again and such mode ”hopping” can continue for quite a long time. The
shape of relaxation oscillations essentially indicates the ability of the laser
resonator to support only the lowest index transverse mode and such simple
observations of relaxation oscillations in free-running mode are readily used
for laser adjustment.
5.2 Q-switching
Q-switching regime is used in order to concentrate laser radiation in a form
of a single pulse and so to increase the peak power of laser output. Q-
switched lasers produce nanosecond pulses with a peak power up to several
hundreds of MW. The first Q-switched laser was demonstrated back in 1962.
At present, Q-switched lasers are used for many scientific and technological
applications such as precision materials processing, laser ranging and detec
tion and pumping ultrashort pulse lasers and amplifiers, to mention a few.
The idea of Q-switching technique is based on the control of losses
in the resonator and consequently the resonator quality. This is
performed as follows: after the pump is turned on, resonator losses are artifi
cially increased. Therefore, the laser generation is prevented, but the active
laser medium at that time efficiently accumulates population acting as an
energy storage. When the population inversion reaches its maximum value,
the resonator quality is rapidly increased (returned to its initial high value
state). Then the gain coefficient rapidly becomes very high (exceeding the
threshold value by many times) and radiation intensity inside the resonator
rapidly (during just several round-trips) increases, ”emptying” the whole ac
cumulated population in the upper laser energy level. Such laser pulse is
often called giant (gigantic) pulse.
Before going into technical details, how Q-switching is performed in prac
tice, we first analyze the basic features of operation of such laser. The prin
cipal scheme of a Q-switched laser is depicted in Fig.65. Initially, let us
estimate the parameters of such laser radiation and factors from which they
depend on. Since the entire radiation (laser generation) process usually oc
5.2 Q-switching 100
curs faster than in 20 ns, we can ignore slow processes, such as change of the
pump rate in time and consequent change of the population.
Figure 65: Principal scheme of a Q-switched laser. r1 and r2 are the mirror
reflection coefficients, l is the length of the resonator, L is the length of laser
medium, R denotes the pump source and Q denotes the resonator Q-switching
unit.
dϕ γcL 1
=ϕ − . (126)
dt n0 l tp
Let us introduce a dimensionless time τ = t/tp , which is normalized to the
photon lifetime in a resonator. Then after multiplying both sides of Eq.(126)
by tp , we get
dϕ γ γ
=ϕ −1 =ϕ −1 , (127)
dτ n0 l/cLtp γt
where γt = n0 l/cLtp is the threshold gain coefficient when number of photons
in the resonator is constant: dϕ/dτ = 0. As the gain coefficient is directly
proportional to population inversion, the former expression can be rewritten
in terms of population inversion:
dϕ n
=ϕ −1 , (128)
dτ nt
where nt = ∆Nt V is the threshold population inversion. The term ϕ(n/nt )
indicates the number of photons generated per unit of normalized time due
to stimulated emission. Because each generated photon is a result of electron
transition between the laser energy levels, after emission of a single photon
the difference of population is reduced by a factor of two: population of
the upper energy level is decreased by one and population of the lower en
ergy level is increased by one. This means that a single optical transition
corresponds to change of population ∆n = −2 and then its evolution is:
dn n
= −2ϕ . (129)
dτ nt
Eq.(128) and Eq.(129) describe the time evolutions of photon number
and population, respectively. Solving these equations yields the radiation
parameters. Dividing Eq.(128) by Eq.(129), we get photon number evolution
which is related to the change of population:
dϕ nt 1
= − . (130)
dn 2n 2
After integration we get:
[
1 n
ϕ − ϕi = nt ln − (n − ni ) , (131)
2 ni
where ϕi and ni are photon number and population, respectively, at the initial
moment of time. Assuming that the initial photon number in the resonator
5.2 Q-switching 102
is infinitely small (ϕi ≈ 0 and this is essentially the case, because at the
very beginning only photons produced by spontaneous emission exist), we
get expression of photon number in the resonator:
[
1 n
ϕ= nt ln − (n − ni ) . (132)
2 ni
The former expression describes the relation between the photon number
and population at any moment of time. Let us take time t ≫ tp , when we
know that there are no photons left in the resonator (ϕ = 0), then we can
evaluate the final population nf :
nf nf − ni
= exp . (133)
ni nt
n −n
The quantity η = i ni f indicates which portion of the accumulated pop
ulation is converted into radiation and η approaches unity when ni /nt → ∞,
i.e., the greater is the initial accumulated population at the moment
of switching to high quality, the more efficient is its conversion to
radiation. The time evolutions of photon number and population, as de
scribed by Eq.(128) and Eq.(129), are plotted in Fig.66.
tp
τd = (134)
ni /nt − 1
and is inversely proportional to the ratio of initial and threshold populations.
The greater is this ratio, the faster is pulse evolution when Q-switching is
turned on. The pulse duration is also a function of this parameter and is
expressed as
η ni /nt
τp = t p . (135)
ni /nt − ln(ni /nt ) − 1
Note that duration of the gigantic pulse is directly proportional to photon
lifetime in the resonator: in order to get shorter pulses, some losses
in the resonator should be introduced, e.g., choosing an output coupler
with greater transmission.
Another important parameter of the laser output is the peak power of
the pulse which is defined as P = ϕhν/tp . According to Eq.(132), the peak
power is expressed as:
[
hν n
P = nt ln − (n − ni ) . (136)
2tp ni
The maximum peak power is achieved at the moment of time when n = nt
(it can be derived from the condition ∂P/∂n = 0). Inserting n = nt to
expression in Eq.(136) we get
[
hν nt
Pmax = nt ln − (nt − ni ) . (137)
2tp ni
Assuming that the initial population significantly exceeds the threshold value,
i.e., ni ≫ nt (this is essentially the case), the expression for peak power sim
plifies to
ni hν
Pmax = . (138)
2 tp
In practical setups Q-switching mode is optimized for either maximum peak
power (minimum pulse duration) or energy, depending on the planned ap
plications. An important feature of the Q-switched regime: the longer is the
photon lifetime in the resonator tp (the greater resonator quality), the lower
is the peak power.
Up to now we have analyzed the case when the quality of resonator is
switched instantaneously. In practice, this means that Q-switching time is
5.2 Q-switching 104
shorter than the development time of the pulse τd . What happens when res
onator quality is switched slowly (the resonator losses α are reduced gradu
ally), is depicted in Fig.67. Then the gain coefficient γ can exceed losses
several times during Q-switching phase (moments t1 and t2 ) and as a result
several pulses are generated, which is often not a desired result.
Figure 67: Dynamics of photon number and gain coefficient in the case of slow
Q-switching. α marks the level of losses in the resonator.
Figure 68: Dynamics of the Q-switching operation in the cases of (a) pulsed pump
and (b) continuous pump.
5.3 Methods of Q-switching 105
When reflection from the prism or mirror aligns with the optical axis of the
resonator, the quality of the resonator suddenly increases and a giant pulse
is generated at this moment. This is a very simple and cheap Q-switching
method, which performs equally well at any laser wavelength. The first Q-
switched ruby laser operated exactly according to this principle. However,
this method has one key drawback: the speed of prism rotation is limited
(in real conditions it is up to several hundred rotations per second), so Q-
switching is slow and, as a consequence, several light pulses are generated.
5.3 Methods of Q-switching 106
Therefore in practice this method is no longer used and has only a historic
value as a simple idea of Q-switching .
Figure 70: The schematic of a Q-switched laser that uses an electro-optic modu
lator. V is the voltage applied to the electro-optical crystal (Pockels cell), PA is
the polarization analyzing element.
Figure 71: Electro-optical modulator based on (a) longitudinal and (b) transverse
Pockels effect.
n3 V
[
ωl
∆ϕ = (no − ne ) − o r63 , (140)
c 2 d
where ne is the refractive index for the extraordinary wave, l is the length of
crystal, d is the width of crystal. Unlike the longitudinal Pockels effect, the
polarization rotation due to transverse Pockels effect depends on the crystal
length. It may seem convenient, since the half-wave voltage can be reduced
significantly, however, excitation of strong high frequency acoustic waves can
considerably modulate the transmission of such electro-optic switch in time.
5.4 Mode-locking
Another laser operation regime is mode-locking. Mode-locking represents
a technique that allows phasing of all frequencies (longitudinal modes)
circulating in the laser resonator. Mode-locking produces the shortest and
the most intense laser pulses with durations varying from several hundred
of ps to < 10 fs (depending on the gain bandwidth of the laser medium
and mode-locking method). The first mode-locked laser was demonstrated
in 196612 . Absolutely all ultrashort pulse lasers employ a certain type of
mode-locking technique. Let us discuss the differences between unphased
and phased (locked) mode laser operation.
Consider a laser medium with inhomogeneously broadened spectral line.
Only longitudinal modes for which the gain is greater than losses are amp
lified. Assuming sufficiently long resonator, a number of such modes can
be very large and the intermodal distance (the difference of adjacent mode
frequencies) which is defined as
πc
∆ω = ωq+1 − ωq = (143)
n0 l
is much smaller than the laser amplification bandwidth ∆ωL , as schematically
illustrated in Fig.73. Only in this case the discussion on mode-locking is
meaningful.
N
N −1
∗
2
I(t) = |E(t) | = EE = E02 ei(ωn t+ϕn ) e−i(ωn t+ϕn ) = N E02 . (146)
n=0
This result shows that the total radiation intensity is a sum of intensities
of the individual modes: I = N Imode . This is the average value because
a certain number of modes may accidentally have the same phase, but the
average result is exactly as specified. Such (unphased mode) radiation is
generated by free-running and Q-switched lasers.
Now let us consider the situation when all longitudinal resonator modes
are phased, i.e., have the same phase ϕn = ϕ0 . Then the the complex field
amplitude is
N
N −1 N
N −1
i(ωn t+ϕ0 ) iϕ0
E(t) = E0 e = E0 e eiωn t . (147)
n=0 n=0
2
2 ∗ 1 − e−iN ∆ωt 2
2 sin (N ∆ωt/2)
I(t) = |E(t) | = EE = E02 = E0 2 . (148)
1 − e−i∆ωt sin (∆ωt/2)
The latter result suggests that phasing of longitudinal modes induces certain
intensity modulation in time with intensity maxima at
∆ωt
= 0, π, 2π, ..., nπ. (149)
2
It is then easy to determine the temporal separation (period) between the
adjacent intensity maxima:
5.4 Mode-locking 113
2 2
2 sin (N ∆ωt/2) N (∆ωt/2)2
I(t)max = lim E0 2 = lim E02 2
= N 2 E02 .
∆ωt/2→0 sin (∆ωt/2) ∆ωt/2→0 (∆ωt/2)
(151)
In the case of mode-locking, intensity of the light pulse is I = N 2 Imode ,
i.e., significantly greater than in the case of modes with random phases.
Therefore, mode-locking produces light pulses with very high intensity. A
comparison of the intensity evolutions I(t) in the cases of unsynchronized
and synchronized modes is shown in Fig.74.
Figure 74: Evolution of light intensity in a resonator in the cases of (a) unsyn
chronized and (b) synchronized mode (mode-locked) cases. The intensities in (a)
and (b) are not to scale.
2π tr 2n0 l 1
τp = = = = . (152)
∆ωN N Nc ∆ωL
This result suggests that the more modes are phased, the shorter is the
generated light pulse. Assuming that all the longitudinal modes are phased,
the minimum pulse duration is inversely proportional to the width of laser
medium amplification band. In a real laser, it is possible to synchronize a
major fraction of longitudinal modes, the relevant examples are presented in
Table 3. Note that with several exceptions (e.g., dye or semicoductor lasers)
only solid-state lasers operate in mode-locking regime as only these lasers
provide sufficiently broad laser amplification bandwidth.
ways and their choice depends on specific properties of the laser medium, i.e.,
the gain coefficient and lifetime of the upper laser energy level. According to
the principle of operation, active mode-locking is performed via amplitude
or phase modulation. The amplitude mode-locking relies on periodic modu
lation of resonator losses, while the phase mode-locking is based on periodic
adjustment of the resonator length. In lasers which have short upper level
lifetime (e.g., dye lasers) active mode-locking is implemented by modulation
of gain coefficient (e.g., by pumping with short pulses from another laser).
Such mode-locking method is called synchronous pumping.
We will now discuss the basic principle of only the most commonly used
active amplitude mode-locking method.
the modulator exactly at the time when resonator losses are reduced again.
Moreover, with each round-trip the pulse shortens due to uneven amplifica
tion: its peak is amplified more efficiently than its fronts. Mode-locking in
the time domain can be understood in a quite straightforward way: when
the pulse gets shorter, its frequency spectrum broadens and this automat
ically implies that more resonator longitudinal modes are made to oscillate
in phase. When such process repeats many times, in an ideal case of in
homogeneous linewidth broadening the minimum pulse duration is defined
as
0.44
tp ≈ . (153)
∆ωL
In the reality, a complete locking of all the longitudinal modes is never
achieved, and the final pulse duration which establishes after many round-
trips is determined by the modulation depth or frequency. In the case of
homogeneous amplification linewidth broadening the minimum pulse dura
tion can be defined as
0.45
tp ≈ √ , (154)
∆ωL fm
where fm is the modulation frequency of resonator losses. The active mode-
locking critically depends on how accurately the modulation frequency matches
the resonator length (or its round-trip time): even a small mismatch between
these two, e.g., length change due to environmental temperature variations,
mechanical vibration, etc, can cause large instabilities in the mode-locking
process, which eventually can be lost.
Using active mode-locking methods, light pulses with durations ranging
between 1 ns and 50 ps are usually generated. This pulse duration interval is
quite unique, as there are no other methods to generate bandwidth-limited
pulses with such duration. Usually, electro-optical modulators are used in
lasers with pulsed pumping, where the gain coefficient is large, while acousto
optical modulators are used in lasers with continuous wave pumping where
the gain coefficient is small.
the remaining noise background. In this way, a single light pulse remains
circulating in the resonator and becomes shorter and shorter with every pass
through the absorber. As a result, its frequency spectrum broadens that is
equivalent to adding more and more longitudinal modes in phase.
We have already mentioned that organic dye saturable absorbers posses
several drawbacks related to their chemical instability (lack of durability).
Furthermore, they are liquid and have to be constantly mixed to avoid un
desired thermal effects. Solid-state saturable absorbers (such as Cr:YAG)
also have drawbacks related to relatively long relaxation time and limited op
erating wavelength range, so they are practically used only for Q-switching.
In 1992 qualitatively new type of saturable absorbers based on semiconductor
media was proposed 13 . These saturable absorbers were called SESAM
(Semiconductor Saturable Absorber Mirror). Unlike organic dye ab
sorbers which are based on saturable transmission, semiconductor absorbers
are based on saturable reflection (Fig.77).
range) and low saturation intensity. Therefore, the use of SESAM is lim
ited to very low pulse energies; such pulses are generated by pulsed lasers
that employ continuous wave pumping (quasi-cw lasers). With high energy
pulses, the coating of the absorber degrades in time as is eventually optically
damaged.
Another method of passive mode-locking is based on the optical Kerr
effect in transparent materials, which results in the intensity-dependent re
fractive index, which is expressed as
n = n0 + n2 I , (156)
where n2 is the nonlinear index of refraction. For all solid-state media that
are transparent in the visible and IR range, this quantity is in the order
of 10−16 cm2 /W. The typical intensities present in laser resonators induce a
refractive index change of the order of 10−5 −10−6 , and this is enough that an
intense light beam starts to self-focus when propagating in such a medium.
Self-focusing occurs because beam in the resonator has intensity distribution
I(r) described by the Gaussian function: beam parts with different intensities
acquire different phase shift, which is the largest at the center of the beam:
2πn2 ln
δϕ(r) = I(r), (157)
λ
where ln is the length of nonlinear medium. Due to that intensity-dependent
phase shift, the initially flat phase front of the beam bends as if the beam
passes through a convex lens and the beam starts to self-focus, as shown in
Fig.78.
Figure 78: Bending of the phase front of intense Gaussian beam due to intensity-
dependent phase shift which causes self-focusing.
The lens induced by nonlinear medium is called Kerr lens and its focal
length
5.5 Methods of mode-locking 120
1 4n2 ln
= Ip (158)
fn n0 w2
is a function of the peak intensity Ip . Here w is the radius of a Gaussian
beam.
Now let us discuss how the self-focusing effect can be exploited for mode-
locking. The method is called Kerr lens mode-locking (KLM) and its prin
ciple of operation is schematically depicted in Fig.79. A nonlinear medium
is placed in the laser resonator together with a diaphragm. The diameter of
the diaphragm is chosen such that only a high intensity beam that carries a
very short light pulse passes the diaphragm without any losses. Only a very
short (mode-locked) light pulse has sufficient intensity to induce a nonlinear
lens, while a noise-like radiation with low intensity experiences large losses.
Hence a Kerr lens operates similarly to a saturable absorber with very fast
relaxation time, which is only few fs (10−15 s), and which is determined by
the response time of electron cloud to the electric field of light. Kerr lens
can be induced in any transparent material, since all media possess electronic
nonlinearity. Moreover, this method is suitable to any wavelength and is non-
resonant: the losses are created without absorption. A small problem is that
often the initial fluctuations appearing in the laser resonator have insufficient
intensity to initialize Kerr lens mode-locking. Then the laser generation is
triggered by an external perturbation e.g., by gently hitting the resonator
mirror mount. Of course, KLM can only be applied in the case of quasi-cw
laser and would be difficult to implement in pulsed pumping regime.
Kerr lens mode-locking was first demonstrated in 1991 by generating
femtosecond pulses from Ti:sapphire oscillator14 . Interestingly, in that con
figuration, sapphire crystal served itself as a nonlinear medium, while the
14
D. E. Spence, P. N. Kean, and W. Sibbett, 60-fsec pulse generation from a self-mode
locked Ti:sapphire laser, Optics Letters 16, 42–44 (1991)
5.5 Methods of mode-locking 121
Summary of Chapter 5
• A unique feature distinguishing lasers from other light sources is
that they can emit light not only in the form of a low divergence
beam but also in the form of very short pulses.
• There are three types of pulsed laser operation regimes: free-
running, Q-switching and mode-locking.
• Free-running is the simplest laser operation regime, which starts
when the laser oscillation threshold is exceeded and produces µs
pulses with kW peak power.
• Q-switching regime is achieved via control of resonator losses and
produces ns pulses with a peak power up to hundreds of MW.
• Mode-locking regime is achieved when longitudinal modes in the
laser resonator are made to oscillate (locked) in phase. Mode-
locking produces ultrashort pulses with durations ranging form
100 picoseconds to a few femtoseconds and peak powers of up to
several hundreds of MW.
6 Generation and amplification of ultrashort light pulses 122
resonator. Let us discuss what these elements are and how they work.
The group velocity dispersion of transparent dielectric media in the vis
ible and near-IR spectral range is normal: the long-wavelength (”red”) spec
tral components propagate faster than the short-wavelength (”blue”) spectral
components and therefore the pulse becomes chirped (phase modulated).
A chirped pulse is ”colored”: its leading front is ”red” and its trailing front
is ”blue”, and the frequency across the pulse changes linearly (the phase
changes quadratically). The pulse duration during propagation changes ac
cording to
�
z2
tp = t0 1+ , (159)
L2d
d2 ϕ λ3 L ∂ 2 n
= , (160)
dω 2 2πc2 ∂λ2
where L is the length of the dispersive medium and n(λ) is the dispersion
law of the refractive index.
Compensation of dispersion requires delaying the spectral com
ponents of the light pulse that propagate faster with respect to
spectral components that propagate slower. There are several ways to
do this.
It would be very simple if one could introduce into the resonator an
optical element, which has anomalous group velocity dispersion
(g0 < 0), where the short-wavelength (”blue”) spectral components propag
ate faster than long-wavelength (”red”) spectral components. Unfortunately,
the range of anomalous group velocity dispersion in transparent dielectric
materials lies close to the mid-IR region, e.g., in glass or sapphire, the dis
persion is anomalous for wavelengths longer than 1.3 µm, while most of the
solid-state lasers emit at shorter wavelengths.
A solution was found by using a pair of prisms, where anomalous group
velocity dispersion occurs due to optical path difference between distinct
6.1 Dispersion control in the resonator 124
Note that after passing the prism pair, the laser beam becomes ”colored”,
as its colors become separated in space due to different refraction angles.
The coloration of the beam is cancelled after back reflection from the mirror
(providing yet another pass), which also doubles the temporal dispersion of
the setup.
If the pulse duration is very short (consequently, its frequency spectrum
is very broad), dispersion cannot be fully compensated using a pair of prisms
because of higher-order dispersion effects (related to higher order derivat
ives: dϕ3 /dω 3 , etc.) that become significant. The effects related to higher
order dispersion manifest themselves as nonlinear distortions of the pulse
chirp preventing optimal compression of the pulse. In that case, optimal
dispersion compensation (pulse compression) is achieved using the so-called
chirped mirrors16 . Chirped mirror is a mirror with a specific dielectric
coating which consists of e.g. alternating SiO2 and TiO2 layers of variable
thickness, deposed in such a way that light penetration depth (where it is
reflected) depends on its wavelength, as illustrated in Fig.82.
Figure 83: Setup for measurement of the autocorrelation function. PD is the beam
splitter, L is the lens, NK is the nonlinear crystal for second harmonic generation,
D is the diaphragm and FD is the photodetector. The inset depicts measured
autocorrelation function of an ultrashort pulse.
∞
IAC (x0 ) = I(x)I(x − x0 )dx. (167)
−∞
The temporal delay in such setup is calibrated very simply: when one
of the pulses is delayed, the image of autocorrelation function on the x axis
shifts accordingly. Since the temporal autocorrelation function is converted
into spatial intensity distribution, the measurement accuracy also depends
on the beam diameter. In this case it is chosen such that it would be much
greater than longitudinal pulse parameters (its duration).
The intensity autocorrelation function has one significant advantage. This
method background free: if the pulses do not overlap in time, the measured
signal is zero. On the other hand, since the autocorrelation function is always
symmetrical even when the intensity profile of the measured pulse is not, or
the pulse shape is complex (i.e. containing several intensity peaks, sub-
pulses, etc.), this kind of information is not retrieved. This shortcoming is
circumvented by the measurement of cross-correlation function. The cross-
correlation function is produced in a similar manner as the autocorrelation
function, just instead of the pulse replica, a pulse with a well-known known
shape (usually much shorter as well), which is called the probe pulse, is used:
∞
ICC (τ ) = I(t)Iz (t − τ )dt, (168)
−∞
where Iz (t) is the intensity of probe pulse. It is obvious that the cross-
correlation function is not symmetrical and the cross-correlation technique
is often applied for the measurement of complex pulses.
The temporal information can also be extracted by measuring higher-
order, e.g., third order correlation function:
∞
I (3) (τ ) = I(t)2 I(t − τ )dt. (169)
−∞
Very often, besides the duration and shape of the ultrashort pulse, the
knowledge of its spectral and temporal phase is required, especially concern
ing optimization of pulse compression setups. Measurement methods that
enable phase retrieval are based on the measurement of instantaneous fre
quency of the autocorrelation function. The measurement setup is basically
identical to that shown in Fig.83, only in the present case instead of a pho
todetector a spectrometer is used. These methods are called FROG; the
abbreviation stands for frequency-resolved optical gating. The obtained res
ults (the so-called frogograms) are processed using mathematical algorithms
which are able to retrieve the actual shape and instantaneous phase of the
pulse.
Figure 85: The layout of a two-stage single-pass laser amplifier. T1 and T2 are
telescopes which increase the size of the beam.
loss of its spectral content by the use of dispersive elements. The stretching
factor may vary from hundreds to thousands of times, in this way significantly
reducing the pulse intensity. The long chirped pulse is then amplified in a
laser amplifier by a factor of 104 − 106 and then recompressed to the original
duration by dispersive elements, which introduce an opposite dispersion as
the stretched does. The CPA technique solved the long-standing problem of
safe and efficient amplification of ultrashort optical pulses without the onset
of optical damage of the amplifier material and other optical components,
enabling a tremendous leap in the peak power and intensity of laser pulses,
boosting an exciting progress in ultrafast laser technology. At present, all
modern ultrashort pulse laser systems are based on this principle19 , while the
inventors of this technique were awarded the Nobel Prize in Physics in 2018.
The number of passes through the active medium is adjusted with the mirrors
and the main advantage of such amplifier is simplicity of the setup which does
not require any complex control elements. However, the mechanical stability
of such amplifier system is low: even a small misalignment of mirrors signi
ficantly alters direction the light beam and small beam distortions occurring
after each pass only increase.
Regenerative amplifier.
The setup of a regenerative amplifier (Fig. 88) is very similar to that
of a Q-switched laser. Here the electro-optic switch (the Pockels cell and
polarization analyzer) performs functions of coupling the seed pulse (the
stretched oscillator pulse) into the resonator and dumping the pulse after
amplification.
The principle of operation of the regenerative amplifier is as follows. Let
us assume that initially the resonator is closed (λ/4 voltage is applied to the
Pockels cell). Then the amplifier medium accumulates population inversion.
When the seed pulse arrives, it is reflected from the polarizer and passes
the Pockels cell which turns its polarization into circular. After reflection
from the mirror and second pass through the Pockels cell, pulse polarization
becomes perpendicular with respect to the initial polarization, therefore, it
passes the polarizer and is coupled into the resonator.
At this moment the voltage applied to the Pockels cell is turned off and
the resonator quality suddenly increases, making the pulse ”trapped” and
amplified in the resonator. The pulse makes several tens to several hundreds
of round-trips in the resonator to reach the gain saturation. When the gain
saturates, the λ/4 voltage is again applied to the Pockels cell which rotates
the polarization by 90◦ (note a double pass) and the amplified pulse is coupled
out (dumped). It is obvious that the regenerative amplifier operates as a Q-
switched laser when there is no seed pulse, but the resonator is opened.
6.4 Pulse stretchers and compressors 135
Figure 88: Layout of the regenerative amplifier. PE is the Pockels cell, PA is the
polarization analyzer.
Figure 89: Principal setup of a diffraction grating pulse compressor. G1, G2 are
the diffraction gratings, V is the mirror.
Figure 90: Layout of a grating pulse stretcher. G1, G2 are the diffraction gratings,
L1, L2 are the telescope lenses, V is the mirror.
If both gratings are placed exactly at the focal planes of the lenses (s = f ),
group velocity dispersion of such system is zero, as the optical paths of all
spectral components are equal, and pulse duration does not change when
passing through such arrangement. Moving the gratings closer to the lenses
(as depicted in Fig.90), the group velocity dispersion of the system becomes
normal: the optical path of ”red” spectral components becomes shorter than
that of the ”blue” spectral components. The closer to the lenses are the
gratings, the larger is the introduced normal group velocity dispersion. The
effective length of such pulse stretcher is expressed as:
hν
Jsat = . (173)
σ
The saturation fluence for dye and most of gas lasers is ∼ 2 mJ/cm2
which means that in order to get, for example, 5 mJ energy ampli
fied pulse, the diameter of the active medium should be around 2 cm.
Therefore, due to low saturation fluence dye and gas lasers are very
inefficient amplifiers. For solid-state lasers the saturation fluencies are
much higher, e.g. 0.9 J/cm2 for Ti:sapphire, 7 J/cm2 for Nd:glass and
32 J/cm2 for Yb:glass. In a 1 cm diameter Ti:sapphire amplifier it is
possible to obtain the amplified pulse with an energy of 1 J.
• Long lifetime of the upper laser energy level, which determines the
ability to accumulate large population inversion.
6.5 Laser amplifier media 139
tspont
M= . (174)
tp Jsat
The evolution of achieved laser powers from the invention of the laser is
illustrated in Fig.91.
Figure 91: The maximum focused laser radiation intensity versus year.
The early rise of the focused intensity was achieved thanks to the invention
of laser pulse generation techniques, Q-switching and mode-locking and sub
sequent amplification of mode-locked pulses. Thereafter, the progress has
slowed down due to physical limits of short pulse amplification, which were
imposed by the optical damage of the laser amplifier material. A major
breakthrough came in 1985, when the chirped pulse amplification technique
was invented, which currently lies on the basis of every modern ultrashort
pulse laser system. Apart from tremendous impact on laser physics, invention
of CPA facilitated rapid development of experimental sciences and opening
new areas of physics, technology and multidisciplinary research. In partic
ular, CPA-based lasers provided experimental access to strong-field physics
6.5 Laser amplifier media 140
and relativistic nonlinear optics, eventually drawing the guidelines for prob
ing vacuum nonlinearity and reaching the fundamental laser intensity limits.
At present, the most powerful laser systems employ Nd:glass and/or Ti:sapp
hire as amplifier media. Fabrication technology of Nd:glass is well developed,
there are many efficient pump sources for this amplifying medium and the
wavelength of Ti:sapphire oscillator can be matched with the Nd:glass emis
sion spectrum. However, Nd:glass is an amorphous medium, its thermal
conductivity is low which significantly reduces repetition rate of the ampli
fied pulses. Moreover, amplification bandwidth of Nd:glass can sustain only
pulses with sub-picosecond duration. The highest peak power of 1.5 PW
was achieved in a hybrid Ti:sapphire and Nd:glass laser system, which pro
duced 0.44 ps pulses with an energy of 660 J, yielding the focused intensity
of > 7 × 1020 W/cm2 20 . More recently, 2 PW femtosecond laser system
based on sole Ti:sapphire laser medium was reported to produce the com
pressed pulses with a duration of 26 fs and an energy of 72.6 J21 . Currently
the highest reported focused beam peak intensity is 1.1 × 1023 W/cm2 22 .
Let us determine what the maximum peak power can be achieved by amp
lifying ultrashort light pulses. Assuming that tp is duration of the amplified
pulse, which is inversely proportional to its spectral width (tp = 0.441/∆ν):
Jsat hν
Pth = = ∆ν. (175)
tp 0.441σ
Assuming that we can focus such a beam into a spot size of the order of a
wavelength, we get the maximum intensity:
Pth hν 3 ∆ν
Ith = = . (176)
λ2 0.441σ c2
It is easy to estimate that for known laser materials Ith ≈ 1024 W/cm2 .
20
M. D. Perry et al, Petawatt laser pulses, Optics Letters 24, 160–162 (1999)
21
Y. Chu et al, High-contrast 2.0 Petawatt Ti:sapphire laser system, Optics Express 21,
29231–29239 (2013)
22
J. Woo Yoon et al, Realization of laser intensity over 1023 W/cm2 , Optica 8, 630–635
(2021)
6.5 Laser amplifier media 141
Summary of Chapter 6
• Dispersive broadening in laser resonator prevents the generation
of ultrashort light pulses carrying very broad spectrum. Disper
sion compensation is performed by introducing a delay of the
fastest (long-wavelength) spectral components with respect to
the slowest (short-wavelength) ones by using prism pairs and/or
chirped mirrors
7 Types of lasers
In this chapter we will briefly discuss the basic features, materials and prin
ciples of operation of solid-state, gas and semiconductor lasers.. The key
difference between these three types of lasers is that the laser radiation is
generated by very different means. In solid-state lasers, the laser radiation is
produced by rare-earth or transition metal ions embedded in a solid-state host
material. In gas lasers the laser radiation is produced by neutral atoms, ions
or molecules, while in semiconductor lasers the laser radiation is generated
at the p-n junction by recombination of electron-hole pairs, i.e. is produced
by the semiconductor crystal itself. These key differences determine the con
struction principles, operation regimes, parameters of output radiation and
the range practical applications of these laser types.
• the laser host material having specific mechanical, optical and thermal
properties;
• the optical pump source with its spectral and temporal properties and
pump geometry.
Solid-state laser materials are crystals or glasses (laser hosts) doped with
small amounts of activator ions which produce laser radiation. In general,
solid-state laser materials have to have sufficient narrow emission lines (ex
cept those used for ultrashort pulse generation), strong absorption bands or
lines and high quantum efficiency of the laser transition. The host material,
where these activator ions are embedded, must be transparent to both pump
and laser radiation, and the absorption of an ion must be in the spectral
7.1 Solid-state lasers 143
Activator ions.
Ions of a certain group of chemical elements termed rare-earth metals,
lanthanide and actinide, serve as activator ions in solid-state lasers. Rare-
earth metal ions (see the two special rows positioned at the very bottom
of the periodic table shown in Fig. 92) possess narrow emission linewidths
(this essentially ensures that the lifetime of the upper laser level is long)
and electronic transitions occur between the inner layers which are not fully
occupied. Since the electronic transitions occur between the inner layers
which are shielded by the outer layers, emission spectrum of such ion is
almost identical to that of a free ion and is almost independent of the laser
host the ion is embedded in.
Xe : 1s2 2s2 2p6 3s2 3p6 3d10 4s2 4p6 4d10 5s2 5p6
In the case of rare-earth metals, electrons further fill the empty orbital of
sublayer 4f :
As a rule, rare-earth metal ions are trivalent; the neutral atoms have lost
three electrons from 6s, 5d (if present) and 4f sublayers:
Laser hosts.
Laser hosts are the materials into which the activator ions are embedded.
Laser hosts have to comply with strict requirements which can be summarized
as follows:
7.1 Solid-state lasers 145
is relatively simple and well elaborated. While doping the sapphire crys
tal, Al3+ ions are replaced with trivalent activator ions. However, there is
a slight problem: the size of Al3+ ion is very small, only ≈ 0.6 Å, while all
rare-earth metal ions are almost twice larger. Therefore, A3+ ions in the
sapphire crystal could be replaced only by transition metal ions: Cr3+ and
Ti3+ , which have similar size (≈ 0.7 Å). Consequently, two doped sapphire
materials are in use: Ti:sapphire and Cr:sapphire (ruby), the latter exists
naturally in nature.
Apart sapphire, other gem class crystals could be doped with trans
ition metal ions: Cr:BeAl2 O4 (alexandrite), Cr:Be3 Al2 (SiO3 )6 (emerald), etc.
Some fluoride crystals, called colquirites, are also doped with transition metal
ions (usually Cr3+ ): LiCaAlF6 (LiCAF), LiSrAlF6 (LiSAF) and LiSrGaF6
(LiSGAF). Compared to sapphire, their thermal properties are significantly
poorer, but they have a very broad amplification bandwidth (they are so
called vibronic materials) that affords either generation of femtosecond pulses
or production of wavelength-tunable laser radiation.
Garnets comprise the most widely used oxide group. Out of these, the
most popular and widespread is yttrium aluminium garnet, Y3 Al5 O12 (YAG).
YAG exhibits outstanding mechanical, thermal, and optical properties, which
makes it a versatile optical material that is widely used in optoelectronics,
laser physics, and nonlinear optics. YAG is one of the most important laser
host materials, which could be doped with various rare earth metal ions, since
radius of yttrium ions (≈ 1.1 Å) is practically identical to rare-earth metal
ions radii. Neodymium doped YAG (Nd:YAG) is one of the most popular
laser materials, which has low generation threshold, large gain and could be
pumped with flashlamps and laser diodes. Currently, several other artificial
garnets were produced: Gd3 Ga5 O12 (GGG), Ga3 Sc2 Al3 O12 (GSGG), etc.
Garnets could also be doped with Yb3+ ions. For example, in the meantime
a lot of research is carried out to develop Yb:YAG based laser systems that are
capable of operating at very high repetition rates and producing ultrashort
light pulses with sub-picosecond duration. Among the other commonly used
laser hosts it is worth to mention yttrium-based crystals: yttrium vanadate
(YVO4 ) and yttrium lithium fluoride YLiF4 (YLF) which all are doped with
Nd3+ ions.
The search for new efficient laser materials continues nowadays, for ex
ample, a recently invented but already commercially used laser material is yt
terbium ion doped potassium gadolinium tungstate KGd(WO4 )2 (Yb:KGW).
Another laser research line is directed toward the development of ultrashort
pulse lasers operating in the mid-IR spectral range. To this end, novel laser
sources based on Ho-doped gain media (Ho:YLF and Ho:YAG), operating
around 2 µm and delivering multimilijoule, multigigawatt pulses with dur
7.2 Gas lasers 147
levels are the same for all laser media, gas media have certain specifics. It
can be briefly described as factors reducing the lifetime of the upper laser
level:
• collisions between excited and not excited atoms. In this way the ex
citations of the upper laser level are lost because they are become re
distributed between the other energy levels;
• Neutral atom gas lasers. The laser media are noble gases and their
mixtures, also metal vapors. All noble gas lasers emit in the infrared,
except the He-Ne laser which emits in the visible spectral range, at
the wavelength of 632.8 nm. Metal vapor lasers emit visible light and
operate only in a pulsed free-running regime. Continuous population
inversion cannot be achieved in these media, because their lower laser
energy level is depleted only through collisions with gas tube walls.
Metal vapor lasers operate only at high (> 1000◦ C) temperature, since
the metal has to be evaporated. The best-known example of metal
vapor laser is copper vapor laser (CVL).
• Ion lasers. The energy levels of an ionized atom are shifted with
respect to those of a neutral atom due to an unbalanced charge. As
a result, the gaps between energy levels increase, and all ion (noble
gas, metal vapor) lasers emit at shorter wavelengths, in the visible and
ultraviolet spectral range. The most important ion gas lasers are argon
ion Ar+ laser and helium-cadmium He-Cd laser.
spectral range. Out of variety of molecular gas lasers, the most import
ant laser is CO2 laser emitting at wavelengths of around 10 µm. CO2
laser is one of the most powerful continuous wave laser (the average
power can exceed kW and more) and one of the most important lasers
for industrial applications, such as laser cutting and welding.
Since the bandgap energy of Ga1−x Alx As depends on the Al mole frac
tion, this double-heterostructure is fabricated so as the potential well for
electrons of height ∆Ec coincides spatially with a well for holes of height
∆Ev . Applying the forward bias of eVa ∼ Eg , large densities of electrons are
7.3 Semiconductor lasers 151
injected from the n side and holes from the p side into the well, providing
the population inversion condition, which is expressed as
EF c − EF v > ℏω,
where EF c and EF v are the quasi-Fermi levels for the conduction and valence
bands, respectively. The GaAs layer where stimulated emission takes place
is called the active region, as shown in Fig. 94(a).
Figure 94: (a) Schematic energy diagram under positive (forward) bias in
GaAlAs/GaAs/GaAlAs double-heterojunction laser diode. (b) The refractive in
dex profile that confines the optical field.
For the maximum gain, the light has to be confined within an active re
gion. This is achieved by a waveguiding effect in a GaAlAs/GaAs/GaAlAs
sandwich, since reduction of the energy gap of a semiconductor by dop
ing, causes increase of its refractive index. The refractive index distribution
and the output profile of a typical double-heterojunction laser is shown in
Fig. 94(b). In this case the gain scales as d−1 , where d is the height of the act
ive region. The concept of double-heterojunction GaAlAs laser is transferred
to other semoconductor lasers, e.g. InAlP.
Conventional heterostructures confine the gain into a small zone (the act
ive region) and prevent electron and hole diffusion. Further miniaturization
of the active region leads to quantized electron motion that is characterized
7.3 Semiconductor lasers 152
by discrete energy levels, giving rise to quantum film, quantum wire and
quantum dot lasers. These lasers are not just simply smaller, but also show
qualitatively novel radiation properties. Structures with reduced dimension
ality offer lower threshold currents, larger gain and lower temperature sens
itivity.
The most unwanted property of a laser diode is the mode hop that results
in the change of laser wavelength and consequently the gain with temperat
ure. To overcomes this, the simplest approach is based on an external cavity
concept, which makes use of a diffraction grating at the Littrow configura
tion, which serves as a cavity mirror. Various modifications of semiconductor
lasers were developed with regard of control the laser wavelength and coher
ence properties by integrating additional components during laser crystal
manufacturing. In such a way, the so-called distributed feedback (DFB),
distributed Bragg reflector (DBR) lasers are designed.
Summary of Chapter 7
• In solid-state lasers the laser radiation is produced by either rare-
earth or transition metal ions embedded in a solid-state host
material (transparent crystals and glasses). Solid-state lasers ef
ficiently operate in all regimes: free-running, Q-switching and
mode-locking.
where ϵ0 is the vacuum dielectric permittivity constant and χ(1) is the linear
optical susceptibility. In general, quantities P (t) and E(t) are vectors and
rapidly oscillate (at the optical frequency) in time. The electric field of a
light wave propagating in a transparent dielectric medium shifts the charges
inducing a dipole moment µ(t), which oscillates at a frequency of the light
wave. The induced polarization is then P (t) = N < µ(t) >, where N is
the number of dipoles and <> denotes the averaging of all dipole moments.
Note that all dipole moments which may naturally exist in the medium are
not accounted for, since they do not emit electromagnetic waves.
8.2 Second harmonic generation 154
The nonlinear optical phenomena occur when the electric field of light is
strong enough to induce the nonlinear polarization, which is however much
smaller that the linear counterpart. This is called perturbative regime of
light-matter interaction or perturbative approximation. In that case, the
polarization can be expressed in power series of the electric field strength:
P (t) = P (1) (t) + P (2) (t) + P (3) (t) + ... = PL + PNL . (181)
Here P (2) (t) = ϵ0 χ(2) E 2 (t) and P (3) (t) = ϵ0 χ(3) E 3 (t) are called second order
(quadratic) and third order (cubic) nonlinear polarizations, respectively. In
what follows, we will further restrict to second order nonlinear optical
phenomena which occur in a certain class of materials: noncentrosym
metric dielectric crystals, i.e. crystals that exhibit birefringence.
P (2) (t) = ϵ0 χ(2) E 2 (t) = 2ϵ0 χ(2) EE ∗ + (ϵ0 χ(2) E 2 e−2iωt + c.c.). (183)
The induced quadratic polarization has two components. The first polariza
tion component (the first term on the right hand side) has zero frequency and
describes the nonlinear phenomenon called optical rectification, i.e. creates
a static electric field in the crystal which exists as long as the optical field is
present. The second polarization component (the second term on the right
hand side) denotes the polarization which oscillates at the doubled frequency
of the incident light wave and is responsible for second harmonic generation.
8.2 Second harmonic generation 155
Figure 95: Energy level diagram depicting second harmonic generation. Eg is the
energy bandgap of dielectric medium. Dashed lines denote virtual energy levels.
Note that virtual energy levels are not real, and they are associated with a
small perturbation of bound electron quantum mechanical state by a light
wave. Therefore, lifetime of virtual energy level can be formally estimated
from the uncertainty principle: ℏ/δE, where δE is the energy difference
between virtual and the closest real energy level. So the lifetime of virtual
energy level is always much shorter than the lifetime of any real energy level,
suggesting that processes involving virtual energy levels are very fast.
Figure 95 also represents the energy conservation law: ℏω1 + ℏω1 = ℏω2 ,
where 2ω1 = ω2 . However, for efficient second harmonic generation, the mo
mentum conservation law also has to be satisfied, i.e., ℏk1 +ℏk1 = ℏk2 , where
k1,2 = n1,2 ω1,2 /c denote fundamental (incident) and second harmonic wave
vector lengths, respectively. Assuming that all interacting waves propagate
at the same direction, we get that the momentum conservation is equivalent
to the condition
with an exciting wave. This in turn guarantees that emitted second harmonic
waves interfere constructively during propagation and the intensity of second
harmonic radiation increases quadratically versus the propagation distance.
The fulfillment of the phase matching condition makes use of crystal bi
refringence and calls for setting appropriate polarizations of the interacting
waves. For instance, in a negative uniaxial crystal (ne < no , where ne is the
principal value of the extraordinary refractive index and no is the ordinary
refractive index), the phase matching condition could be satisfied for the fun
damental wave of an ordinary polarization, and the second harmonic wave
of an extraordinary polarization. This is achieved only for certain propaga
tion direction of the interacting waves with respect to the optical axis of the
crystal, and which is called the “phase matching angle”:
Figure 96: Phase-matching angle for the second harmonic generation in BBO and
KDP crystals.
coupled, and the intensities of the signal and idler waves grow exponentially
with propagation distance, at the cost of the pump wave, whose energy is
continuously depleted.
Figure 98: Phase-matching curve for BBO crystal. The pump wavelength is
400 nm, which corresponds to the second harmonic of Ti:sapphire laser.
1
Is (z) = Is (0) exp(2Γz) (187)
4
ωi
Ii (z) = Is (0) exp(2Γz) (188)
4ωs
here
ωi ωs d2eff Ip2
Γ2 = , (189)
ns ni np ϵ0 c3
and quantity Γ is called the amplification increment, which depends on
the effective nonlinearity of the crystal deff , proportional to certain matrix
elements of the second order nonlinear susceptibility χ(2) , and pump wave
intensity Ip . Here we notice that only radiation intensity of a laser source
is high enough for amplification, while using non-laser sources makes optical
parametric amplification practically impossible.
In this sense the nonlinear crystal performs as a certain catalyzer and
gain coefficient is very high (conversion from pump energy to signal and idler
waves can be more than 80%). The gain coefficient of optical parametric
amplification is defined as:
Is (z) 1
G= = exp(2Γz) , (190)
Is (0) 4
here Is (0) and Is (z) are the signal wave intensities at crystal input and output
respectively. At GW/cm2 level pump intensity, which is easy to achieve with
8.4 Optical parametric amplifiers 160
ultrashort pulse lasers, the gain coefficient can be 103 − 106 for a crystal of
a few mm length. So optical parametric amplification is a unique process
allowing fixed wavelength radiation to be converted into tunable wavelength
radiation with high efficiency and without loss of coherence properties.
The first type of optical parametric amplifiers are made of nonlinear crys
tal which is placed in a resonator as shown in Fig.99. Such optical parametric
amplifiers are called standing wave amplifiers or optical parametric oscillators
(OPO). Here the pump beam is coupled through one of the resonator mir
rors and the mirrors are chosen to reflect either signal or idler waves (singly
resonant OPO, Fig.99(a)) or both (doubly resonant OPO, Fig.99(b)).
For OPO pumping either low power lasers which operate in CW or Q-
switched modes or femtosecond laser oscillators generating ultrashort but low
energy pulses at high repetition rate are usually used. Since the peak power
of such pump laser radiation is low, single-pass amplification is relatively low
as well (from several to a dozen of times). The resonator performs the same
roles of energy accumulation, feedback and formation of high spatial quality
beam (due to diffraction losses) as in a conventional laser. During the first
round-trip parametric superfluorescence is generated in the nonlinear crystal;
a part of which propagating on the resonator axis is amplified during the
following round-trips through the nonlinear crystal. As in conventional laser,
the oscillation condition is defined in the same way: parametric amplification
after one resonator round-trip must exceed losses due to light scattering,
absorption, diffraction and mirror transmission.
In the case of the simplest singly resonant OPO a small fraction (usually
few percent) of resonating signal is coupled out. In addition, the entire idler
wave is coupled out which after each pass is generated again. Notice that
all resonator stability criteria discussed at the beginning of this textbook
also apply to OPO resonator. Wavelength tuning of OPO is performed in
a simple way – by rotating the nonlinear crystal, i.e., by changing pump
beam direction with respect to crystal optical axis. In such case at each
time collinear phase-matching is achieved for a different signal and idler fre
quency pair. Wavelength tuning range of OPO is usually limited by spectral
characteristics of resonator mirrors because fabricating mirrors with uniform
8.4 Optical parametric amplifiers 161
Figure 99: The main configurations of optical parametric oscillators: (a) singly
resonant, (b) doubly resonant.
Figure 100: The principal setup of traveling wave optical parametric amplifier:
SHG - second harmonic generator (optional), BS - beamsplitters, DM - dichroic
mirrors, PSFG - parametric superfluorescence generator, OPA - optical parametric
amplification stages.
Figure 101: Tuning ranges of BBO crystal based optical parametric amplifiers
pumped with commonly used ultrashort pulse lasers and their harmonics.
lifier24 in CPA system has led to novel high power amplification technique,
termed Optical Parametric Chirped Pulse Amplification (OPCPA).
OPCPA favorably combines the advantages of CPA with the advantages of
optical parametric amplification, which offers very high gain, low thermal ef
fects, great wavelength flexibility and intrinsically broad amplification band
width, extending well beyond that afforded by existing solid-state laser amp
lifiers
At present, compact optical parametric chirped pulse amplifiers comprise
a unique class of ultrafast light sources, which currently amplify octave-
spanning spectra and produce few optical cycle pulses with multi-gigawatt
to multi-terawatt peak powers and multi-watt average powers, with carrier
wavelengths spanning a considerable range of the optical spectrum.
24
A. Dubietis, G. Jonušauskas, A. Piskarskas, Powerful femtosecond pulse generation by
chirped and stretched pulse parametric amplification in BBO crystal, Optics Communic
ations 88, 437–440 (1992)
8.4 Optical parametric amplifiers 164
Summary of Chapter 8
• Nonlinear optics offers simple and efficient way to change the
wavelength of laser radiation.
9 EN - LT Glossary
Lasers and laser-related technologies are an integral part of modern science and
are widely used in many areas of everyday life. This textbook provides the basic
principles of laser operation, describes the underlying physical effects a nd gives
an overview of relevant scientific and technological innovations that make possible
generation and amplification o f u ltrashort l ight p ulses. T he m aterial presented
here is intended for Vilnius University Faculty of Physics undergraduate students
of Physics and Light engineering study programmes, but is also suitable to anyone
with a background in natural sciences, interested in this subject.
knygynas.vu.lt
journals.vu.lt