Autocorrelation of An Ultrashort Laser Pulse To Determine Its Duration
Autocorrelation of An Ultrashort Laser Pulse To Determine Its Duration
Autocorrelation of An Ultrashort Laser Pulse To Determine Its Duration
net/publication/263351073
CITATIONS READS
0 1,054
2 authors, including:
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Autocorrelation of an ultra short laser pulse to determine its duration. View project
All content following this page was uploaded by Tamba William Saa on 25 June 2014.
22 May 2014
Submitted in partial fulfillment of a structured masters degree at AIMS South Africa
Abstract
Various techniques to find the duration of an ultrashort laser pulse have been developed. This work is a
report of an experiment to determine the duration of an ultrashort laser pulse with 770 nm wavelength
in both the temporal and spectral domain using autocorrelation technique. Autocorrelation technique
is a technique based on the interaction of a pulse with a time delayed copy of itself within a medium. In
this experiment I will exploit the phenomenon of second harmonic generation (SHG) and amplification
of ultrashort laser pulse. Using measured values from the plots I will produce from the autocorrelation
trace and spectrum respectively I will subsequently determine the pulse duration and spectra width in
both frequency and wavelength. Lastly I will compute the spectral bandwidth product and determine if
the pulse is transform limited or not.
Declaration
I, the undersigned, hereby declare that the work contained in this research project is my original work, and
that any work done by others or by myself previously has been acknowledged and referenced accordingly.
i
Contents
Abstract i
1 Introduction 1
3 Experiment 10
3.1 Experiment Set Up . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.2 White Light generation Stage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.3 First amplification stage of white light . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.4 Second Amplification Stage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.5 The Output and Compressor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.6 Autocorrelator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
5 Conclusion 21
References 22
ii
1. Introduction
Atomic nuclei vibrates around their equilibrium positions, binding electrons move into outer orbitals,
reactants form new compound due to the arrangement of atomic nuclei and electrons. All of these pro-
cesses takes place within only a few femtoseconds (10− 15), fractions of a second so small like comparing
one second to the age of the universe. This temporal window is experimentally accessible by ultrafast
laser spectroscopy, which utilizes laser pulses of durations down to 10 fs to start photoinduced dynamics
in matter and to probe the evolution of the excited system at later times. However, temporal resolu-
tion of these experiments cannot be achieved by any electronic detector, but must be accomplished by
purely optical techniques which in essence take laser induced snapshots of the dynamics at mechanically
determined delayed times.
My goal in this experiment is to introduce the formation of ultrashort pulse, describe the autocorrelation
operation, measure an ultrashort laser pulse widths in frequency and time using autocorrelation tech-
nique. In the first part of my work I will explain the theory of the time and frequency domains than use
the Fourier transform method to measure the pulse duration.The second part of this work describes the
experimental set up with detail diagrams at the Laser Research Center at Stellenbosch University where
the experiment was conducted. Data of the autocorrelation out put signal was collected and plotted
to determine the pulse duration and spectral bandwidth. These results were use to calculate the time
bandwith product and compare to the theoretical value of a transform limited pulse.
1
2. Theory: The Ultrashort Laser Pulse
A laser is an optical source that emits coherent radiation at a specific wavelength. Ultrashort laser pulse
is an optical pulse with duration of picoseconds or less. It is one billionth of one millionth of a second.
10−15 10−12 10−9 10−6 10−3 100 103 106 109 1012 1015 1018
1picosecond
Femtosecond (fs) light pulses are electromagnetic wave packets and as such are fully described by the
time and space dependent electric field.
A pulse is characterized by measurable quantities that can be directly related to the electric field.
Consider the temporal dependence of the electric field neglecting its spatial and polarization dependence,
i.e, E(x, y, z, t) = E(t). We can completely describe a pulse in either the time or the frequency domain
as we will do in the next section using Fourier transform to achieved this.
The frequency domain is used to measure the spectrum and the spectral phase (phase of the electric
field in the frequency domain) while in the Time-Domain we measure the intensity and the phase.
2.1.1 Time-Domain. In order to measure the duration of an event, you must use a shorter event. But
than to measure the duration of a shorter event, you must use an even shorter one and so on. So, how
do we measure the duration of an ultrashort laser pulse? We will use a pulse to measure itself. Now
our problem becomes measuring a pulse by itself.
Consider the energy time uncertainty principle given by
~
∆E∆t ≥ , (2.1)
2
∆E represents the standard deviation in energy and ∆t is the amount of time it takes the expectation
value of some operator to change by a standard deviation.
2
Section 2.1. Time- and Frequency-Domain Measurements Page 3
2
t
E(t) = exp − 2 cos(ω0 t), (2.3)
τ
c
which is plotted below using the values τ = 30 fs , ω = 2πf where f = ; λ = 530 nm is the pulse
λ
2πc
wavelength and c is the speed of light in vacuum. Hence, ω = = 3.554 · 1015 Hz.
λ
The frequency domain sometimes refers to as the domain of the spectrometer is generally used to
measure the spectrum, S(ω), of a laser pulse.
The measurement of the duration τ of this laser pulse is achieved by finding the Fourier transform of
the electric field S(ω) of the pulse.
2.2.1 Calculating the Spectrum. We will analytically calculate the Fourier transform of the spectrum
of the Gaussian shape light pulse using Fourier transform technique.
The Fourier transform, is a mathematical transformation used to transform signals from time (or spatial)
domain to frequency domain. It is reversible, being able to transform from either domain to the other.
The term itself refers to both the transform operation and to the function it produces.
We therefore calculate the spectrum S(ω), from the electric field, E(t), using the Fourier transform
method. Z ∞ 2
t
E(ω) = exp − 2 cos(ω0 t) · cos(ωt)dt. (2.4)
−∞ τ
Z ∞ 2 Z ∞ 2
t t
E(ω) = exp − 2 cos[(ω0 − ω)t]dt + exp − 2 cos[(ω0 + ω)t]dt.
0 τ 0 τ
Now we expand and simplify the above expression using the Maclaurin series, a special case of Taylor
expansion. By definition
∞
X (−1)n 2n
cos(x) = ·x . (2.6)
(2n)!
n=0
Therefore, for x = (ω0 − ω)t for the first integral and x = (ω0 + ω)t for the second integral we have
∞ ∞ ∞ ∞
(−1)n −t2 (−1)n −t2
X Z X Z
2n 2n 2n
E(ω) = (ω0 − ω) exp t dt + (ω0 + ω) exp t2n dt.
(2n)! 0 τ2 (2n)! 0 τ2
n=0 n=0
Let
t √
= x
τ
Section 2.2. The Frequency Domain Page 5
√
t = τ x,
τ
dt = √ dx.
2 x
Therefore
∞ ∞
" #Z
X (−1)n 2n
X (−1)n 2n
∞ √ τ
E(ω) = (ω0 − ω) + (ω0 + ω) exp(−x) · (τ x)2n · √ · dx,
(2n)! (2n)! 0 2 x
n=0 n=0
∞ ∞
" #
(−1)n (−1)n τ 2n+1 ∞
Z
X X 1
2n 2n
E(ω) = (ω0 − ω) + (ω0 + ω) xn− 2 · exp(−x) · dx.
(2n)! (2n)! 2 0
n=0 n=0
similarly Z ∞
1 1
Γ n+ = xn− 2 exp(−x)dx,
2 0
hence,
∞ ∞
" #
(−1)n (−1)n τ 2n+1
X
2n
X
2n 1
E(ω) = (ω0 − ω) + (ω0 + ω) ·Γ n+ .
(2n)! (2n)! 2 2
n=0 n=0
By definition √
1 π · (2n)!
Γ n+ = n . (2.8)
2 4 · (n!)
Therefore,
" ∞ ∞
# √
X (−1)n 2n
X (−1)n 2n τ 2n · τ π · (2n)!
E(ω) = (ω0 − ω) + (ω0 + ω) ,
(2n)! (2n)! 2 4n · (n!)
n=0 n=0
√ "∞ ( 2 )n ∞
( 2 )n #
τ π X ω0 − ω 1 X ω 0 + ω 1
E(ω) = (−1)n τ + (−1)n τ ,
2 2 n! 2 n!
n=0 n=0
√ " #
τ (ω0 − ω) 2 τ (ω0 + ω) 2
τ π
E(ω) = exp − + exp − .
2 2 2
For very large ω and ω0 the second term vanishes. Hence our spectrum simplifies to
√
τ (ω0 − ω) 2
τ π
E(ω) = · exp − . (2.9)
2 2
P∞ xn τ (ω0 ± ω) 2
x
where e = n=0 for x = . ω0 and ω are optical frequencies in the order of 1014th ,
n! 2
τ is the pulse duration assume to be longer than a few periods of the wave, being a few tens of
femtoseconds.
Section 2.2. The Frequency Domain Page 6
1 1
· g(x) = · C.
2 2
If we let h = x at half maximum we have
h2
1
· C = C · exp − 2 ,
2 2σ
h2
1
= exp − 2 .
2 2σ
Taking natural log on both sides gives us
h2
1
ln = ln exp(− 2 ) ,
2 2σ
1 h2
ln = − 2,
2 2σ
1
−2σ 2 ln = h2 ,
2
−2σ 2 (ln 1 − ln 2) = h2 ,
−2σ 2 (− ln 2) = h2
√
h = σ 2 ln 2,
hence √
FWHM = 2 · σ · 2 ln 2 (2.11)
Section 2.2. The Frequency Domain Page 7
y
C
1
C
2
x
3σ 2σ σ σ 2σ 3σ
√
2·σ· 2ln2
2.2.3 Product of Pulse Durationand Spectral Width. If we compare equation 2.9 with the definition
1 (x−µ)2
of a Gaussian P(x) = √ · exp − 2σ2 ,
σ 2π
we have
1 τ2
= ,
2σ 2 4
2
σ2 = 2 ,
τ
√
± 2
σ= . (2.12)
τ
√
2 2
Only the positive root is physically relevant here. Hence our spectral width is .
τ
Now our time-bandwidth product is
∆ω · ∆τ ≥ k, (2.13)
∆ω · ∆τ ≥ 0.441, (2.14)
where ω is the frequency bandwidth measured at full width at half maximum (FWM) and τ is FWHM
in time and, k is a number which depends only on the pulse shape. Therefore,
Section 2.2. The Frequency Domain Page 8
Figure 2.4: Simulated intensity of a Gaussian laser pulse plotted against frequency
The central frequency is 3.55 Hz and the full width at half maximum of the curve is 0.1 fs.
Section 2.3. Intensity Autocorrelation Page 9
Figure 2.5: Simulated intensity of a Gaussian laser pulse plotted against its wavelength
The intensisty of this field is proportional to the product of the intensities of the two input pulses:
Since common detectors are too slow to resolve this beam in time, they will measure
Z ∞
(2)
A(τ ) = I(t)I(t − τ )dt. (2.17)
−∞
(2)
A(τ ) is second order autocorrelation, the subscript (2) indicate second order autocorrelation, third order
generation is also possible.
3. Experiment
3.1 Experiment Set Up
The following set up is a schematic diagram of a Nonlinear Optical Parametric Amplifier (NOPA). A
NOPA is an optical device or set up which receives an input signal and generate an output signal with
higher optical power based on nonlinear interactions.
BS
CPA Output
λ = 775nm, τ = 150f s
m
SHG
λ = 387nm BS BS
BS
1:9 BBO 30:70
R:T R:T
m
neutral focusing Collimating
density m
filter lens lens
OPA m
m
iris Short pass
diaphragm saphire BBO
crystal filter
Prism Compressor
D-shape
m
to autocorrelator
m
m m
m
m
m m
OPA
BBO
Input stage
In this stage a beam of a 775 nm wavelength laser, is emitted by Ti:sapphire laser (also known as
T i : Al2 O3 laser, titanium sapphire laser, or simply Ti:sapphires) and is splitted by a beam splitter. The
pulse is further splitted with approximately 4% of it been reflected. The frequency of our transmitted
10
Section 3.2. White Light generation Stage Page 11
pulse is doubled in a Second Harmonic Generation (SHG), a 0.77mm Beta Barium Borate (β −BaB2 O4 )
crystal and is used as our pump pulse.
Figure 3.2: Frequency is doubled inside a nonlinear crystal Troiano (May 14, 2013)
The reflected pulse at this point, low in intensity, compare to our pump pulse at a ratio of 1 : 9 is
reflected by a broadband mirror into the white light generation stage.
The term white light refers to a continuous spectrum, including all wavelengths from green to red and
way into the infrared spectral range.
To generate white light we start by intensity control of the beam, that is the modulation of a 775 nm
pulse by passing it through a neutral density filter. Upon exiting the filter the beam pass through an
iris diaphram, which controls the numerical aperture of the incident beam. At this point it is important
to note that the beam intensity and aperture are controlled in order to create white light filament with
a discretized energy of approximately 2µJ. After generating white light, a 30 mm plano-convex fused
silica lens is used to align and collimate the beam.Here again a filter is used to improve relative intensity
distribution of the white generated light. The filter blocks all light above 700 nm of wavelength.
Figure 3.3: White generation inside a sapphire crystal (Troiano, May 14, 2013)
Section 3.3. First amplification stage of white light Page 12
Figure 3.4: A focusing lens, sapphire crystal and, collimating lens in the white light generation stage(?)
22% of the frequency doubled blue pump pulse from the second harmonic generator (SHG) is reflected
by a beam splitter on to a focusing mirror with focusing length of 150 mm. This highly reflective
focusing mirror than reflects the blue laser light into a 1 mm BBO crystal where it is overlapped with
the seed light from the white light stage. Diagram and explaination of white light amplification inside
the optical parametric amplifier (OPA) is given in later section.
Again, it is important to note that in order to avoid damaging of the crystal, the pump light is focus in
front of the crystal such that the diameter of the pump and the seed pulses match within the crystal.
In this stage the amplified seed pulse is directed past two silver mirrors and focus by a silica lens onto a
BBO crystal. Coming from another path is the transmitted pump light which is reflected via two steering
mirrors and a focusing mirror at an adjustable angle unto the BBO crystal. Here the pre-amplified white
light from the first amplification stage and the pump light are carefully phase-matched to avoid getting
an output power that is the sum of two individual powers from the amplification of different spectral
regions. The output is therefore, a multiplicative amplification.
Section 3.5. The Output and Compressor Page 13
χ2
ω2
ω2 seed
seed ω2 ω2
ω2 signal
ω1 ω3 ω3
pump ω1 idler
Process of white light amplification. In the BBO crystal a pump pulse of frequency ω1 is excited by a
weaker seed pulse of frequency ω2 . This stimulation results in the creation of a signal pulse of frequency
ω2 and an idler pulse of frequency ω3 due to the anihilation of the pump pulse. That is, ω1 = ω2 + ω3 .
Different colors propagate at different velocities in a medium experiencing an effect called group-delay
dispersion (GDD). Since ultrashort pulses have broad spectra, GDD results in significant pulse length-
ening or intensity increase.
To compensate for the GDD of the NOPA output, a pulse compressor, a silica Brewster prism pair is used.
A prism pulse compressor is an optical device used to shortened the duration of a positively chirped
ultrashort laser pulse by giving different wavelength components different time delays. A positively
chirped ultrashort laser pulse is a pulse in which light with longer wavelength travels faster than shorter
wavelength light. Overall the prism compressor makes wavelength components of a laser pulse overlap
with each other, thus causing a shorter pulse.
Prism 1 in Figure (3.1) refracts different wavelength components of light to slightly different angles
causing these different wavelength components to travel in different paths but in parallel directions.
Prism 2 also refracts all components as was done by prism 1 onto a reflecting mirror.These refracted
wavelength components are retro reflected by the mirror, again to let them propagate in parallel directions
retracing their paths.
3.6 Autocorrelator
The amplified pulse from the NOPA is directed to a BBO crystal via reflecting mirrors and a pair of
irises. The irises align the beam to hit the two D-shape mirrors precisely at the middle of the mirrors.
The two D-shape mirrors split the incoming NOPA output pulse into two identical halves and reflect
the split pulses onto a focusing mirror which in turn focus the pulses into a BBO crystal. In the BBO
crystal the pulses are convoluted and frequency is doubled. Emerging from the BBO crystal we observe
three outcomes; pulse a interact with itself, pulse a interact with pulse b and, pulse b interacts with
itself. The use of a and b here is only meant to indicate two pulses and don’t refer to chronological
Section 3.6. Autocorrelator Page 14
sequence. The frequency doubled pulses are then filtered to block visible light and transmit UV light
pulses and remaining light is then focus on to a detector by a plano convex focusing lens.
m from
NOPA
iris iris
m
m
filter m
m BBO
crystal
iris
•
m
focusing detector
lens
Figure 3.6: A Sketch of the autocorrelator. Depicted is the path of the NOPA output pulse to the
autocorrelator. A photodiode detector register the signal from the autocorrelation.
4. Results and Measurements
We used an autocorrelator to measure the intensity of the convoluted pulse. We achieved this by
splitting the frequency doubled pulse using an autocorrelator, time delaying one of the splitted pulses
than recombining them in a SHG crystal. The out put signal current from the autocorrelator was than
converted into voltage and read on a digital oscilloscope. This signal out put is plotted in the figure
below against time.
Figure 4.1: Second order Intensity (a.u.) autocorrelation signal of a 530 nm pulse plotted against time
(fs).
This autocorrelation signal can be fitted with a Gaussian as shown in the following plot. The accuracy
of this measurement is limited to about 25% due to several experimental parameters, which are not
trivial to quantify.
15
Section 4.2. Calculating the Pulse Duration. Page 16
Figure 4.2: Intensity (a.u.) of a 530nm pulse of a Second order autocorrelation signal from the au-
tocorrelator plotted against time (fs) and fitted against a Gaussain function. Horizontal dashed lines
intersect vertical dashed lines at points where the intensity value drop to 1/e.
From the plot of the pulse intensity against time in Figure (4.2) we calculate the full width at half
maximum of the pulse (τF W HM ),
√
Using equation (2.12) FWHM = τ = 2 · σ · 2ln2 we have
√
τAC = 2 · σ · 2ln2,
√
τAC = 2 · 18.1 · 2ln2,
Section 4.2. Calculating the Pulse Duration. Page 17
τAC = 42.6 fs
1
DAC = √ , (4.2)
2
Therefore,
1
τ = √ · τAC ,
2
1
τ = √ · 42.6 fs.
2
4.3 Spectrum
As an extra we calculate the bandwidth spectrum and also fit it with a Gaussian function.
Figure 4.3: Spectrum of a < 40 fs, 530 nm pulse plotted as intensity (a.u) versus wavelength (nm).
The uncertainty of the experimentally determined spectral width is about 10% again due to several
experimental parameters, of which probably the most influencial one is the flutuation of the pulse
spectrum from pulse to pulse.
Section 4.3. Spectrum Page 19
Figure 4.4: The Spectrum above in Figure 4.3 fitted against a Gaussian function. Horizontal dashed
lines intersect vertical dashed lines at points where the intensity value drop to 1/e.
From spectrum plot in Figure 4.4 we extract values of the amplitude of the signal, central wavelength
of the pulse, and the standard deviation.
amplitude = 4.02e4, central wavelength = 530e2, standard deviation = 7
4.3.1 Calculating the Bandwidth.. Bandwidth is the width of some frequency or wavelength range.
A common definition of bandwidth is the full width at half maximum (FWHM).
Spectral width in terms of wavelength of the ultrashort laser pulse
√ plotted in Figure (4.4) is calculated
below. Again we use equation 2.12 (FWHM = τ = 2 · σ · 2ln2) to calculate the FWHM of the
spectrum √
∆λ = 2 · σ · 2ln2,
√
∆λ = 2 · (7) · 2ln2,
∆λ = 17.4 nm.
Section 4.3. Spectrum Page 20
c
ω= , (4.3)
λ
but for small wave length interval we have
c
∆ω = · ∆λ
λ2
2.99 · 108 m/s
= · 17.4 nm (4.4)
(5.29 · 102 nm)2
= 1.87 · 1013 Hz
= 18.7 THz.
5. Conclusion
Using data from the experiment, I plotted intensity against time and extracted the pulse duration of the
ultrashort laser pulse as full width at half maximum of the autocorrelation trace. We determined the
pulse duration to be 30.1 fs. We also plotted intensity against wavelength and calculated the wavelength
spectral width of the pulse to be 17.4 nm and the frequency spectral width as 18.t THz.
Now, we check if our pulse is transform limited. Transform limit or Fourier limit is the lowest possible
limit for a pulse duration for a given optical spectrum of a pulse. The minimum time bandwidth product
as shown in equation 2.14 is 0.441 for a Gaussian pulse. This time bandwidth product indicates if a
pulse is transform limited.
For my pulse the time bandwidth product is
This result means that our ultrashort laser pulse is nearly transform limited.
Our measured pulse is longer than the transform limited pulse due to chromatic dispersion. Frequency
is dependent on phase velocity in transparent medium. A related quantitative measure of chromatic
dispersion is the group velocity dispersion. Optical filters can also affect the spectral width or shape
of an ultralshort pulse. Reduction of the spectral width may result into temporal broadening. Light
reflection from mirrors also causes dispersion. This temporal broadening occurs because the reflected
pulse has specific reflectivity at different wavelengths of the spectrum.
21
References
C. Antoncini. Ultrashort laser pulses. Department of Physics University of Reading, pages 1–35, May
2014.
J.-C. Diels and W. Rudolph. Ultrashort laser pulse phenomena. Academic press, 2006.
H. Giessen, L. F. Hoos, and J. Teipel. White light lasers and their applications. pages 1–2, May 2014.
O. Gobert, D. Rovera, G. Mennerat, and M. Comte. Linear electro optic effect for high repetition rate
carrier envelope phase control of ultra short laser pulses. Applied Sciences, 3(1):168–188, 2013.
E. J. Rohwer. Development of non-collinearly phase matched optical parametric amplifier and application
in pump-probe spectroscopy. Department of Physics Stellenbosch University, March 2011. Master
Thesis.
J. Troiano. Laser science. light can do way more than just bend. Sustainable Nano, page 1, May 14,
2013.
22