Download Modeling of Lithium Niobate Ridge Waveguides - Rays

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts

Harold Hopkins (physicist) wikipedia , lookup

Optical tweezers wikipedia , lookup

Thomas Young (scientist) wikipedia , lookup

Ellipsometry wikipedia , lookup

Fourier optics wikipedia , lookup

Nonimaging optics wikipedia , lookup

Photon scanning microscopy wikipedia , lookup

Magnetic circular dichroism wikipedia , lookup

Dispersion staining wikipedia , lookup

Retroreflector wikipedia , lookup

Birefringence wikipedia , lookup

Anti-reflective coating wikipedia , lookup

Refractive index wikipedia , lookup

Surface plasmon resonance microscopy wikipedia , lookup

Nonlinear optics wikipedia , lookup

Silicon photonics wikipedia , lookup

Transcript
Modeling of Lithium Niobate Ridge Waveguides
Anton Möller
[email protected]
under the direction of
Associate Professor Katia Gallo
Department of Applied Physics,
Quantum Electronics and Quantum Optics
Royal Institute of Technology
Research Academy for Young Scientists
July 10, 2014
Abstract
Waveguides are used in order to guide electromagnetic waves and an example of a waveguide is optical fibers. Problems arise with waveguides when propagation losses are induced, and the aim of this study is to model optical guidance, field profiles and effective
indices in ridge type waveguides and propagation losses in these. This is done by using a model of the waveguides with dimensions experimentally measured, in COMSOL.
Simulations were done for two waveguides with different dimensions at the wavelengths
λ = 1550 nm and its second harmonic λ = 775 nm. The results show that the waveguide
with smaller dimensions have properties which imply a larger confinement with less propagation modes. When studying these waveguides we could also see experimentally that a
second harmonic generation was generated. In the future, more exact modeling should be
done. The asymmetry of the waveguides as well as manufacturing imperfections should
be taken into account when building the simulation which should give even more reliable
data.
Contents
1 Introduction
1
2 Optical Waveguide Theory
1
2.1
General Optics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1
2.2
Waveguides . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
2.2.1
2D Waveguides . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
2.2.2
3D Waveguides . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4
2.3
Refractive Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4
2.4
Second Harmonic Generation . . . . . . . . . . . . . . . . . . . . . . . .
5
2.5
Manufacturing of LiNbO3 Ridge Type Waveguides . . . . . . . . . . . . .
6
2.6
Propagation Modes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
7
2.6.1
Single-Mode Range . . . . . . . . . . . . . . . . . . . . . . . . . .
9
2.6.2
Full Width at Half Maximum . . . . . . . . . . . . . . . . . . . .
9
2.7
Coupling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9
2.8
Wave Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10
2.8.1
Wave Equation for TM Mode in a 2D-Slab Waveguide . . . . . .
11
2.8.2
Wave Equation for TE Mode in 2D-Slab Waveguide . . . . . . . .
11
Propagation Losses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
12
2.9.1
Fresnel losses . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
12
2.9.2
Coupling Efficiency . . . . . . . . . . . . . . . . . . . . . . . . . .
12
2.9.3
Propagation Losses in the Waveguide . . . . . . . . . . . . . . . .
13
2.9.4
Determining the Total Power Loss . . . . . . . . . . . . . . . . . .
13
2.9.5
Determining the Experimental Sampling Rate . . . . . . . . . . .
14
2.9
3 Method
3.1
15
Technical Data of the Waveguides . . . . . . . . . . . . . . . . . . . . . .
4 Result
17
18
2
4.1
Mode Profiles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
18
4.2
Comparison with Experimental data . . . . . . . . . . . . . . . . . . . .
21
5 Discussion
22
6 Acknowledgments
23
A Values of Coefficients for Sellmeier Equations
26
B Code: Full Width at Half Maximum
26
1
Introduction
Waveguides play an important role in the development of optical integrated circuits as it
can be viewed as one of the fundamental building blocks in this area. They are used in
order to guide electromagnetic waves and an example of a system where optical integrated
circuits are active is the fiber communications system. When dealing with waveguides, it is
important that the propagation losses are small to ensure high efficiency and transmission
of light. The ridge type waveguides used in this study are foremost utilized in order to
study the second harmonic generation and non-linear interactions and modeling these
waveguides is an important step in achieving this. The aim of this study was to model
ridge-type lithium niobate waveguides (LiNbO3 ) and propagation losses in these. This has
not been done previously, and this study may provide an approach for modeling other
types of waveguides as well.
2
Optical Waveguide Theory
In this Section, general optics which is necessary in order to fully understand the theory
of waveguides will be presented. After that, the underlying theory of this study will be
covered.
2.1
General Optics
Light travel through different media, and each medium has a specific refractive index,
commonly denoted nmedium and defined as
c
nmedium = ,
v
(2.1)
where c is the speed of light in vacuum and v is the speed of light in that medium. This
definition implies that the refractive index for vacuum is 1 and that the refractive index
describes how much faster light travels in vacuum than in the specific medium [6]. When
1
n2
θr
θi
(a)
θi
θc
θr
(b)
n1
(c)
Figure 1: (a). Light travels from a medium with refractive index n1 to a medium with
refractive index n2 with the incident angle θi and the refracted angle θr . This light is
refracted according to Snell’s law. (b). If the incident angle is greater than the critical
angle, total internal reflection occurs and light is confined in the medium. (c). When the
refracted angle is parallel to the interface (θr = π/2 rad), the incident angle is called the
critical angle, θc .
light travels from one medium to another with different refractive indices, it obeys Snell’s
law,
n1 sin θi = n2 sin θr ,
(2.2)
where n1 and n2 are the refractive indices for medium 1 and 2 respectively, θi is the
incident angle and θr is the refracted angle. The angles are measured from a line normal
to the interface between the two media [5]. When the incident angle causes a refraction
parallel to the interface (θr = π/2 rad), it is denoted the critical angle. If the incident
angle is greater than the critical angle; total internal reflection occurs and light cannot
penetrate the interface between the two media and is therefore confined in one of the
media [5]. These processes are illustrated in Figure 1.
2
Figure 2: A cross-section of a typical optical waveguide. The core has a higher refractive
index than the surrounding cladding, ncore > ncladding which means that the propagating
wave is propagating mostly in the core.
2.2
Waveguides
Waveguides guide electromagnetic waves in a specified direction. They are typically composed of several layers with different refractive indices. Innermost, there is the core in
which the propagation of light takes place. The core has a higher refractive index than
the surrounding layer in order to be able to guide the wave. The surrounding layer, also
known as the cladding or the substrate vary in size but is generally larger than the core [6].
A cross-section of a simple optical waveguide is seen in Figure 2.
2.2.1
2D Waveguides
Waveguides are classified either by how they guide waves or how their refractive indices
change. In two dimensional space; step-index waveguides and graded-index waveguides
are considered. A two dimensional waveguide confines light in one direction only. In a
step-index waveguide, the refractive index changes abruptly in the transition between core
and cladding. The graded-index waveguide on the other hand has a gradual change and
a hard line cannot be drawn between the two different media. This change in refractive
indices result in confinement of the wave as it is not able to escape the core. Because
of confinement in one direction only; the wave expands in the non-confinement direction
due to diffraction as it propagates. This means that the width of the beam increases [6].
3
(a) Buried type
(b) Ridge type
Figure 3: (a) The buried type is characterized by the core being submerged into the
cladding. This type of waveguide is suitable for constructing optical waveguide devices
due to relatively low propagation losses. (b) Ridge type waveguides are on the other hand
characterized by the core being on top of the cladding. Ridge waveguides typically induces
propagation losses due to tight confinement of the wave. [6].
2.2.2
3D Waveguides
When confinement takes place in two directions, x and y, the waveguide is a three dimensional waveguide. These waveguides are further divided into subtypes, depending on
the physical structure of the waveguide. Examples of this are the buried and ridge types.
In a buried type waveguide; the core is submerged into the cladding, and in a ridge type
waveguide the cladding is etched off around the core [6]. This is illustrated in Figure 3.
There are other types of three dimensional waveguides and descriptions of these can be
found elsewhere [6].
2.3
Refractive Index
A wave propagating inside the waveguide along the propagation axis sees the effective
index, neff , and not the refractive index, ncore of the core itself. The effective refractive
index is calculated by solving the wave equation described in Section 2.8 for the spatially
changing refractive index distribution.
n(λ, f, z) = nsubstrate (λ, f ) + ∆n(z, λ),
(2.3)
where ∆n(z, λ) is the refractive index increase induced by annealing process described
in Section 2.5, z is spatial coordinates and λ the wavelength. The refractive index of the
substrate, nsubstrate , is dependent on the wavelength of the light propagating and this
4
index is determined by using the Sellmeier equation for LiNbO3
s
nsubstrate (λ, f ) =
a1 + b 2 f +
λ2
a2 + b 2 f
a4 + b 4 f
− a6 λ 2 ,
+ 2
2
− (a3 + b3 f )
λ − a25
(2.4)
where ai and bi are coefficients determined experimentally [4]. These values are found in
Appendix A. f is a temperature parameter defined as
f = (T − 24.5)(T + 570.82),
(2.5)
where T is the current working temperature expressed in degrees Celsius [4].
In order to account for the calibrated dispersion factor of the annealing process, another Sellmeier equation is used [1],
r
∆n(λ) = c1 +
λ2
c2
,
− c23
(2.6)
where ci are coefficients which are found in Appendix A.
2.4
Second Harmonic Generation
In non-linear optics, generating a second harmonic (SHG) in optical devices is studied.
Suppose an incident wave which approaches an optical device with an angular frequency
of ω. When the wave is propagating through the device, depending on the structure,
two photons with individual energies of ~ω may merge into a single photon with the
combined energy of the two 2~ω. This conversion implies that the newly made photon
has a frequency twice that of the merging photons. This also means that the wavelength
of the second harmony is half than of the original wave. SHG is useful because the
conversion of a wave from a lower frequency to a higher enables more types of optical
integrated circuits to be constructed. One example of this is being able to convert a laser
of a specified wavelength into a laser with half that wavelength [3].
5
2.5
Manufacturing of LiNbO3 Ridge Type Waveguides
The manufacturing process of ridge type waveguides begins with a sample of a material,
in this case lithium niobate, LiNbO3 . Pure benzoic acid is then heated and one side of the
sample is then exposed to the melted substance. This is done in order to replace some of
the lithium ions with protons; hence the name proton exchange (PE) for the process. The
proton exchange process is used in order to locally change the refractive index (∆n(λ)
from Equation 2.6) of a region within the waveguide and this exchange is achieved due
to the melted benezoic acid increasing the optical refractive index of the material.
The depth of which the lithium is replaced by protons is determined by the temperature and duration of the process. When the PE has taken place for a fixed amount of
time; the sample is extracted. The sample is now a slab step-index waveguide with a
guiding layer composed of the volume in which the lithium has been replaced and the
manufacturing process could potentially stop here [2].
In order to refine the waveguide; the next step is annealing. During the annealing the
sample is heated for a fixed amount of time in order to allow the protons which have
replaced lithium to reposition themselves. This process allows the protons to sink further
down into the sample, eliminating the discontinuity of the refractive index creating a
graded-index slab waveguide [2].
In order to create the ridge-type waveguide; the now graded-index slab waveguide is
cut several times along the z-axis (the axis of propagation) and this yield the finished
waveguide. The waveguide itself is now finished, but in order to ensure a high quality
wave propagation its surfaces must be polished carefully. This includes both end faces and
the areas in beteween the individual ridges. When a waveguide has been manufactured
through PE and annealing it is generally called an Annealed Proton Exchange (APE)
waveguide. In Figure 4 a microscopic top-view image of a waveguide is shown.
6
Figure 4: Sample waveguide. If carefully looked upon, imperfections and residue from the
manufacturing process are seen. The waveguide is 7.48 µm wide.
2.6
Propagation Modes
Propagation of a wave in a waveguide only takes place according to different propagation
modes which depend on the wave characteristics and the waveguide. The simplest propagation mode is single-mode propagation. When in single-mode propagation all waves,
even with different wavelengths, travel with equal speed. Single-mode propagation is characterized by the output of the waveguide forming a single sharp dot, compared to the
multi-mode propagation which often forms different patterns due to interference between
waves with different propagation speeds [8].
Guided modes are further classified by the orientation of their electrical and magnetic
fields relative to the plane of propagation (polarization). The most common modes in this
context are called transverse magnetic (TM) and transverse electric (TE). Suppose that
there is an electromagnetic wave propagating in the z-direction. When the electrical field
of the wave is oscillating along the y-axis, the wave is polarized in transverse magnetic
mode. Conversely, when the magnetic field is oscillating in the y-direction, the wave is
polarized in transverse electric mode [8].
When describing a propagating mode the convention is TMmn or TEmn where the
7
(a) TM00
(b) TM10
Figure 5: Two different modes of propagation of transverse magnetic waves. The images
represent the electrical field distribution on a cross section of the waveguide 5 µm wide
and 4 µm high with a refractive index of 2.1401 in the substrate and are to scale. (a)
Single-mode propagation with the mode TM00 . (b) Multi-mode propagation with the
mode TM10 .
TM and TE part states if the mode is transverse magnetic or transverse electric while
m and n are two modal orders which are obtained by solving an eigenvalue equation
from the resulting process of solving the wave equations found in Section 2.8. For the
purpose of this study it is enough to say that m = n = 0 is the fundamental mode which
is equivalent with single-mode propagation and that all other combinations of m and n
imply multi-mode propagation [6].
In Figure 5 two different modes of propagation is seen for transverse magnetic polarized waves at the wavelength λ = 1550 nm.
Waveguides only guide a discrete set of modes, which is limited by a number of factors.
Firstly, the difference between the refractive index of the core and the substrate index,
∆n = ncore − nsubstrate , bounds the range of the solutions to the wave equation. This is
due to the constraints on the effective refractive index, nsubstrate < neff < ncore [6].
When solving the wave equation, the effective mode index is calculated and for a certain value of the index, the electrical field distribution is obtained. Although the number
of solutions already is limited, not all solutions to the wave equation guides a mode which
is seen by looking at their respective electric field distribution [3].
8
2.6.1
Single-Mode Range
The single-mode range of a waveguide is obtained by varying the width of the waveguide
in the COMSOL simulation and measuring the width at which the effective index goes
towards the refractive index of the substrate. This gives the lower boundary of the singlemode range (or cutoff for the TM00 mode). In order to find out the upper boundary value;
the width at which the TM10 mode starts appearing with an effective index meeting the
constraints from Section 2.6 is determined. The range in which single-mode is the only
propagating mode is then defined as wmin < w < wmax where w represents the width of
the waveguide. Observe that fundamental mode propagation is not limited to the singlemode range while the single-mode range specifies in which range the TM00 mode is the
only propagating mode [3].
2.6.2
Full Width at Half Maximum
Full width at half maximum (FWHM) is defined as the absolute value of the distance
between two x-coordinates where the y-coordinates is equal to half of the maximum yvalue. The coordinates were the normalized electrical field distribution versus location in
this study. FWHM is applied to a function representing a pulse or a wave where the data
can be approximated by either a Gaussian fit or a Lorentz fit. Matlab code to calculate
FWHM is found in Appendix B. This value is utilized when comparing and modeling
different modes since it describes in which span the propagation of the transverse side of
the mode takes place.
2.7
Coupling
Coupling is a way to connect a source of waves, typically a laser, to a waveguide. There
are several ways of doing this, and in this study end-fire coupling is described and used. In
order to excite modes in a waveguide; a wave with a similar profile to those modes is fired
directly onto the polished end face of the waveguide. If the wave does not have a similar
9
Figure 6: End-fire coupling. The waves from the light source are focused on the beginning
of the waveguide which lies in the focal point of the lens. The distance between the source,
lens and waveguide can be adjusted in order to obtain the most efficient incident angle.
profile as the waveguide, chances are that there will be no guidance. Often there is a lens
between the source and the waveguide which focuses the light in order to maximize the
efficiency. When utilizing this coupling method; it is important that the incident wave
profile is as similar as possible to that of the guided mode. This is due to the fact that
a waveguide which has an unlike wave profile with respect to the wave may not guide
well [6]. A diagram of an end-fire coupling arrangement is illustrated in Figure 6.
2.8
Wave Equations
From Maxwell’s equations for a lossless dielectric medium, it follows that for an electrical
field variation Ẽ and a magnetic field variation H̃
∇ × Ẽ = −µ0
∂ H̃
,
∂t
and
∇ × H̃ = 0 n2
∂ Ẽ
,
∂t
where µ0 is the magnetic permeability of free space, 0 the electric permittivity of free
space, n the refractive index for the medium, t the time and ∇ the curl operator.
Suppose that a is wave propagating in the orthogonal coordinate system (x,y,z) along
the z-axis with the propagation constant β = 2πneff /λ. The electric and magnetic fields
10
of the wave varies according to
Ẽ = E(x, y) · ei(ωt−βz) ,
(2.7)
H̃ = H(x, y) · ei(ωt−βz) ,
(2.8)
and
with the angular frequency ω = 2πc/λ, and the imaginary unit i, the wavelength λ and
the time t.
These fields are independent of y in a two-dimensional slab waveguide and taking the
partial derivatives yield
∂
∂
∂
= iω,
= iβ,
=0
∂t
∂z
∂y
.
Equations (2.7) and (2.8) result in two different modes with orthogonal polarized
states, the transverse electric and transverse magnetic modes. The TE mode has an electrical field orthogonal to the plane of propagation (the xz-plane), and therefore consists
only of the field components Ey , Hx and Hz . The TM mode would consequently only
have the field components Hy , Ex and Ez [8] [6]. With this information, one could express
the wave equations for both the TM and TE mode.
2.8.1
Wave Equation for TM Mode in a 2D-Slab Waveguide
∂ 2 Hy
+ (k02 n2 − β 2 )Hy = 0
∂x2



E = β H
x


Ez =
2.8.2
ω0 n2
(2.9)
y
(2.10)
∂Hy
1
iω0 n2 ∂x
Wave Equation for TE Mode in 2D-Slab Waveguide
∂ 2 Ey
+ (k02 n2 − β 2 )Ey = 0
∂x2
11
(2.11)



Hx = − β Ey
ωµ0
(2.12)


Hz = − 1 ∂Ey
iωµ0 ∂x
where k0 = 2π/λ.
When solving either of these wave equations; the result is an electrical or magnetic
field as a function of position in the x,y plane. From this the field distribution is utilized
in order to check whether a waveguide actually guides a mode or not. This is achieved by
determining if the field distribution stabilizes and forms a guiding shape such as a singleor multi-mode [6].
2.9
2.9.1
Propagation Losses
Fresnel losses
Between any two interfaces of two different refractive indices there is a Fresnel loss, R.
Fresnel losses occur due to light being reflected back towards the source by the interface.
This loss is estimated for weakly guiding structures by
R=
(neff − 1)2
,
(neff + 1)2
(2.13)
where neff = βλ/(2π) is the effective refractive index and β is the propagation constant.
This is a good approximation for waveguides with perfectly polished end-faces, but in
reality there are imperfections caused by the manufacturing process and polishing which
increases the losses [7][3].
2.9.2
Coupling Efficiency
During the coupling losses occur due to the waveguide and the source being misaligned.
Losses are also induced by misalignment of the focal point causing the light beam to
not being focused on the coupling point and by the shape of the incoming beam not
matching the mode in the waveguide. This loss is minimized by optimizing the coupling
12
experimentally [6].
The coupling loss is approximated by comparing a profile image of the waveguide to
a simulated one and compare the difference in area of the two. This yields an estimate of
how well the coupling is performed.
2.9.3
Propagation Losses in the Waveguide
When the light has entered the waveguide itself, a propagation loss may occur due to imperfections on the surface. This loss is interesting when trying to optimize the waveguides
because of the length dependent exponential power drop. The formula for calculating the
power at the end of the waveguide, Pend , with a starting power Pstart is
Pend = Pstart · e−αL ,
(2.14)
where α is the attenuation constant and L being the length of the waveguide [6].
2.9.4
Determining the Total Power Loss
Define transmission as T := 1 − R, i.e. how much power that is transmitted through an
interface due to some loss factor. Consider a waveguide for which the losses in power are
to be determined. If the incoming power is denoted Pin , the power which continues after
the interface, P1 , is calculated as follows where T is determined from R in Equation 2.13.
P1 = T Pin ,
(2.15)
During propagation through the waveguide the power is decaying according to Equation (2.14). This equation is rewritten as
P2 = P1 · e−αL ,
13
(2.16)
Figure 7: With an incoming power Pin at the first interface, the power P1 is the output
power from that interface due to a Fresnel loss. During propagation through the waveguide, P1 decreases to P2 according to Equation (2.16) and P2 hitting the end interface
which results in an output power of Pout which is calculated by Equation 2.18.
where P2 is the power at the end interface of the waveguide. At the end interface there
is once again a Fresnel loss which is calculated in a similar way of Equation (2.15).
Pout = T P2 ,
(2.17)
Combine equations (2.15), (2.16) and (2.17) as follows and add a term, ηcoupling , to
account for the coupling losses
Pout = (1 − R)2 Pin e(−αL) ηcoupling = T 2 Pin e(−αL) ηcoupling ,
(2.18)
in order to get an expression for the final output power. The different spots where the
power is considered are illustrated in Figure 7.
2.9.5
Determining the Experimental Sampling Rate
When sampling power as a function of the wavelength λ, it is important to use the correct
sampling step ∆λ in order to get useful data. The power as a function of wavelength is
sampled in order to be able to calculate the attenuation coefficient which describes how
the power attenuates as a wave propagates through a waveguide. Since the graph of this
sampling resembles the graph of a sinusoidal function, using a too big sampling step
would result in a data set where the crests and troughs are not visible whereas a too
14
small step would result in unnecessary data.
There are two ways of calculating this; and the differential method is used here.
Assume that
∆φ
∂φ
=
,
∂λ
∆λ
where φ = 2βL = 4πneff L/λ,. From the Taylor series of φ it follows that
∆φ =
∂φ
∆λ,
∂λ
where ∆φ = 2β(λ0 + ∆λ)L − 2β(λ0 )L = 2π and λ0 is an arbitrary point in space if the
higher order terms are neglected due to their minimal impact for this rough measurement.
Consequently,
∂φ
∆λ = 2π,
∂λ
∂
[4πneff Lλ−1 ]∆λ = 2π,
∂λ
−4πneff L
∆λ = 2π,
λ2
∆λ =
−1 λ2
,
2 neff L
(2.19)
The results given by this calculation is only an approximation of the distance between
two fringes of a curve. Therefore, the value given here cannot be used directly as the
sampling step while this would result in only getting measurements for every fringe and
nothing in between. The calculated value must therefore be reduced by a factor in order
to get good experimental results.
3
Method
In this study, two selected waveguides were studied at λ = 1550 nm and at the second
harmonic generation (λSHG = 775 nm). The two waveguides used are seen in Figure 8.
By utilizing Matlab and COMSOL, numerical solutions to Equation (2.9), the wave
15
equation, are obtained for a modeled waveguide. The COMSOL model consists of a ridge
waveguide defined by areas of different refractive indices. Matlab is used simultaneously
and provides COMSOL with nessecary data about how the refractive index of the core
changes depending on position, as the waveguide is a graded one which is seen in Section 2.5. The effective index is calculated in Matlab by using the Sellmeier equations
found in Section 2.3. COMSOL returns the mode electric field distribution and the mode
effective index. Therefore; multiple solutions are obtained for every simulation, one field
distribution per effective index.
By looking at the distribution by hand, the TM00 mode and its associated effective
index are selected. The electric field distribution is then compared with experimentally
found mode profiles of the waveguide to determine the coupling efficiency. The profile
image is taken by utilizing an infrared camera at the end of a coupled waveguide setup.
To get accurate results from experimental trials; the sampling step ∆λ is obtained
by using Equation (2.19) with approximate values for the length, effective index (which
is approximated as the refractive index of the core in this context) and the wavelength
around which the sampling takes place. This is done in order to model the power losses.
In order to properly model the waveguides, the dimensions of the waveguides must
be known. These values are measured using a surface imaging camera. Furthermore, the
refractive index of the substrate and the wavelength of the electromagnetic wave used
must be known or determined to simulate propagation.
For this study, the single-mode range, mode profiles and effective indices were calculated for the mode TM00 at the wavelength λ = 1550 nm and at the second harmonic
generation (λSHG = 775 nm) in order to model the waveguide. The calculations were based
upon measured dimensions of the waveguide. For the single-mode, FWHM is calculated
by using a Matlab function found in Appendix B.
16
(a) Waveguide 4
(b) Waveguide 6
Figure 8: Both images are taken with a microscope. (a) Microscopic image on the fourth
waveguide with a width of 8.47 µm. (b) Microscopic image of the sixth waveguide with
a width of 7.48 µm
3.1
Technical Data of the Waveguides
The technical specifications of the waveguides used in these experiments are found in
Table 1.
Table 1: Technical data of the waveguides used. Entries beginning with PE are specifications for the proton exchange process and entries beginning with APE are details about
the annealing process.
Property
Waveguide 4 Waveguide 6
Type
Ridge
Ridge
Length
10 mm
10 mm
Width
8.47 µm
7.48 µm
PE Environment
Benzoic acid
Benzoic acid
PE Duration
14 hours
14 hours
PE Temperature
160◦ C
160◦ C
PE Depth of exchange 700 nm
700 nm
APE Environment
Air
Air
APE Duration
26 hours
26 hours
◦
APE Temperature
330 C
330◦ C
17
Table 2: The effective refractive index for different propagation modes in the two
waveguides at the fundamental (λ = 1550 nm) and the second harmonic generation
(λSHG = 775 nm). These values were acquired from COMSOL simulations where the
waveguides had been modeled. A dash indicates that the mode is not guided and therefore no effective index exists for that specific waveguide at a certain wavelength.
Waveguide 4
Wavelength λ = 1550 nm λ = 775 nm
TM00
2.1554047
2.213887
TM10
2.148547
2.212509
TM20
−
2.210211
TM30
−
2.20699
TM40
−
2.202841
TM50
−
2.19776
TM60
−
2.19187
4
Waveguide 6
λ = 1550 nm λ = 775 nm
2.153548
2.21376
2.146548
2.212
−
2.209063
−
2.204945
−
2.19964
−
2.193139
−
−
Result
The single-mode range of the waveguides are 2.70 µm < w1550 < 5.70 µm and 0.895 µm <
w775 < 1.92 µm for the fundamental and second harmonic generation respectively. At the
fundamental wavelength, the refractive index of the substrate was determined according
to Equation (2.3) to be 2.1401 and at the second harmonic 2.1814. Table 2 shows the
effective refractive indices from the COMSOL simulations for the different propagation
modes and wavelength.
4.1
Mode Profiles
Waveguide 4 and 6 both theoretically support the TM00 and TM10 modes at a wavelength
of λ = 1550 nm. The mode profile images are seen in Figure 9. For the second harmonic
generation (λSHG = 775 nm), waveguide 4 supported modes ranging from TM00 to TM60 .
Waveguide 6 supported all modes that waveguide 4 did, except TM60 .
The cross section profiles of some of the modes are shown in Figure 10. For the TM00
modes, the full width at half maximum values are calculated and the results are shown
in Table 3.
18
(a) Waveguide 4 - TM00
(b) Waveguide 4 - TM00 (SHG)
(c) Waveguide 6 - TM00
(d) Waveguide 6 - TM00 (SHG)
Figure 9: Electric field distribution cross section of waveguide 4 and 6 at the fundamental
and second harmonic generation in the x,y-plane. (a) and (b) show the field distribution at
single-mode (TM00 ) for the fundamental and second harmonic generation respectively of
waveguide 4. (c) and (d) show the field distribution at the TM00 mode for the fundamental
and second harmonic generation of waveguide 6. All images are to scale.
Table 3: Calculated values for full width half maximum of the two different waveguides
at the fundamental and second harmonic generation. The data was retrieved from COMSOL simulations and a cross-section of the electrical field distribution used. FWHM is
expressed terms of micrometers.
Waveguide 6
Waveguide 4
Wavelength [nm]
1550
775
1550
775
Full Width Half Maximum [µm]
5.1826
5.1293
5.8226
5.7780
19
60
Electrical field norm [V/m]
Electrical field norm [V/m]
60
50
40
30
20
10
50
40
30
20
10
0
0
0
2
4
6
8
10
0
2
Position [µm]
(a) Waveguide 4 - TM00
6
8
10
(b) Waveguide 4 - TM00 (SHG)
60
Electrical field norm [V/m]
60
Electrical field norm [V/m]
4
Position [µm]
50
40
30
20
10
0
50
40
30
20
10
0
0
2
4
6
8
10
0
Position [µm]
2
4
6
8
10
Position [µm]
(c) Waveguide 6 - TM00
(d) Waveguide 6 - TM00 (SHG)
Figure 10: Cross section profiles of the electric field distribution in the x,y-plane of the two
different waveguides at the fundamental and second harmonic generation. These graphs
were then used in order to calculate the FWHM value.
20
Figure 11: Mode profile of waveguide 6 at the fundamental mode taken with an IR imaging
camera. The camera was an Alpha NIR InGaAs with a field of view of 30 × 30 µm with
an objective of 10x and with a focal length of 16.9 mm.
4.2
Comparison with Experimental data
The absolute value of the experimental sampling rate, |∆λ| was determined to be 0.056 nm
according to Equation (2.19) with the following values:




λ = 1550 nm




neff = 2.1401






L = 10 mm
and a sampling step of 0.005 nm was therefore used in the experiments.
With the effective refractive index the transmission factor for the Fresnel losses were
calculated according to (2.13) yielding T = 0.8682 at the fundamental mode. The coupling
efficiency, ηcoupling , was determined to be 0.41 by using Equation (2.18) with the following
21
values




Pout = 3.09585 mW








Pin = 10.0 mW




T = 0.8681








L = 10 mm





α = 0.8128
obtained by calculations based on experimentally found values.
5
Discussion
Both waveguides have a single-mode range of 2.70 µm < w < 5.70 µm at the fundamental
mode. This means that they should guide a combination of the TM00 and TM10 modes due
to the actual widths of the waveguides being above 7 µm. In Figure 11 we see the guided
mode profile and the single-mode propagation seems to be dominating which prohibits
us from seeing the multi-mode profile. Below the main concentration of the beam there
is a small dot, which probably is a diffracted part of the propagating waves or scattered
light.
Waveguide 4 is wider and therefore the refractive indices are greater than in waveguide 6 at the fundamental mode. Due to its width, waveguide 4 is theoretically able to
guide up to TM60 at the second harmonic generation but this is difficult to observe experimentally. Waveguide 6 has a smaller width which in this case only affects the second
harmonic generation limiting the guided modes to up until TM50 . Both waveguides supports the same modes at the fundamental wavelength. In Figure 9 the guided modes for
the fundamental and second harmonic generation are shown. Once again, observe that
the modes in waveguide 4 are wider and not as confined as in waveguide 6. If we look at
the FWHM, a narrower beam is observed at waveguide 6 due to the larger confinement
which is observed both at the fundamental and second harmonic generation.
22
In order to model a waveguide and its losses it is useful to simulate propagation
and determine characteristics such as the single-mode range and which modes that are
guided. These characteristics serves as a reference of what one should expect to see
in experiments but the modeling and experiments are dependent on one another. For
example, the simulations generate a single-mode range and an effective index while the
experiments determine dimensions and power measurements. One could determine the
effective index and single-mode range experimentally but the process is time-consuming
and inefficient. Improvements of this study would include an automated way of calculating
the single-mode range and some type of automation for the extraction of information from
COMSOL and Matlab when simulating the waveguides.
For future studies it would be interesting to model the second harmonic generation
more thoroughly and study the non-linear interactions in these waveguides. Another
study would be to model the imperfections and the end-face reflections in COMSOL and
to model the exact structure of the ridge and not assume that it is homogeneous. This
would yield a more accurate model and hence more accurate results.
6
Acknowledgments
First, I would like to thank my mentor Associate Professor Katia Gallo for all support
given and all time spent guiding me and for allowing me to work with her. I would then
like to thank all people at the department of Quantum Electronics and Optics for the
encouragement and helpfulness. Hampus Gummesson Svensson, another Rays student,
worked with a project similar to mine but with a different point of view. I would therefore
like to thank Hampus for this amazing journey we’ve been through together and for always
being there for me to bounce ideas off.
All the people making Rays possible have my sincere gratitude, especially my counselor
Dennis Alp, project director Philip Frick and of course counselors Mariam Andersson and
Agnes Nordquist.
23
I would also like to thank ABB, Teknikföretagen and The Royal Patriotic Society for
making Rays possible.
24
References
[1] Bortz, M., Arbore, M., and Fejer, M. Quasi-phase-matched optical parametric
amplification and oscillation in periodically poled LiNbO3 waveguides. Optics letters
20, 1 (1995), P. 59–60.
[2] Bremer, M. Periodically poled ridge waveguides in LiNbO3 . KTH physics, TritaFYS 2012:30 (2012), 17ff, 44ff.
[3] Hecht, E. Optics 4th edition. Optics, 4th Edition, Addison Wesley Longman Inc 1
(1998), P. 13f, 27f, 396ff, 641ff.
[4] Jundt, D. H. Temperature-dependent sellmeier equation for the index of refraction
in congruent lithium niobate. Optics Letters 22, 20 (1997), P. 1553–1555.
[5] Nilsson, K., Pålsgård, J., and Kvist, G. Ergo Fysik 2. Liber Förlag, 2010.
[6] Nishihara, H., Haruna, M., and Suhara, T. Optical integrated circuits. New
York (1989), P. 7–8, 10–15, 27ff, 225ff.
[7] Regener, R., and Sohler, W. Loss in low-finesse ti: LiNbO3 optical waveguide
resonators. Applied Physics B 36, 3 (1985), P. 143–147.
[8] Yariv, A., and Yeh, P. Photonics: Optical electronics in modern communications
(the oxford series in electrical and computer engineering). Oxford University Press,
Inc. (2006), P. 110–136.
25
A
Values of Coefficients for Sellmeier Equations
Table 4: Coefficients for the Sellmeier equation (equation 2.4) [4].
Coefficients
Value
a1
5.35583
a2
0.100473
a3
0.20692
a4
100
a5
11.34927
a6
1.5334 · 10−2
b1
4.629 · 10−7
b2
3.862 · 10−8
b3
−0.89 · 10−8
b4
2.657 · 10−5
Table 5: Coefficients for the Sellmeier equation (equation 2.6) [1].
Coefficients
B
Value
c1
3.43 · 10−5
c2
1.01 · 10−5
c3
0.326
Code: Full Width at Half Maximum
function [ fwhm ] = fw_hm( data )
%fw_hm Calculates the full width half maximum value for
% an input ’data’ with two column vectors x,y by fitting a gaussian
26
% curve to the vectors. By Anton Moeller and Philip Frick 2014-07-08.
load(data,’x’,’y’);
if isequal(size(x), size(y))
%Make them into orderd pairs
values = [x,y];
%Half max of y
half_max = max(y)/2;
%Curve fitting
[curve,gof]
= fit(x,y,’gauss8’);
%Output fitting data
rsquared = gof.rsquare
adjrsquared = gof.adjrsquare
%Prepare equation environment
syms f(x)
f(x) = formula(curve);
curve_coeffs = coeffvalues(curve);
%Load coeffficient values into the function
f(x) = subs(f(x),’a1’,num2str(curve_coeffs(1)));
f(x) = subs(f(x),’a2’,num2str(curve_coeffs(4)));
f(x) = subs(f(x),’a3’,num2str(curve_coeffs(7)));
f(x) = subs(f(x),’a4’,num2str(curve_coeffs(10)));
f(x) = subs(f(x),’a5’,num2str(curve_coeffs(13)));
27
f(x) = subs(f(x),’a6’,num2str(curve_coeffs(16)));
f(x) = subs(f(x),’a7’,num2str(curve_coeffs(19)));
f(x) = subs(f(x),’a8’,num2str(curve_coeffs(22)));
f(x) = subs(f(x),’b1’,num2str(curve_coeffs(2)));
f(x) = subs(f(x),’b2’,num2str(curve_coeffs(5)));
f(x) = subs(f(x),’b3’,num2str(curve_coeffs(8)));
f(x) = subs(f(x),’b4’,num2str(curve_coeffs(11)));
f(x) = subs(f(x),’b5’,num2str(curve_coeffs(14)));
f(x) = subs(f(x),’b6’,num2str(curve_coeffs(17)));
f(x) = subs(f(x),’b7’,num2str(curve_coeffs(20)));
f(x) = subs(f(x),’b8’,num2str(curve_coeffs(23)));
f(x) = subs(f(x),’c1’,num2str(curve_coeffs(3)));
f(x) = subs(f(x),’c2’,num2str(curve_coeffs(6)));
f(x) = subs(f(x),’c3’,num2str(curve_coeffs(9)));
f(x) = subs(f(x),’c4’,num2str(curve_coeffs(12)));
f(x) = subs(f(x),’c5’,num2str(curve_coeffs(15)));
f(x) = subs(f(x),’c6’,num2str(curve_coeffs(18)));
f(x) = subs(f(x),’c7’,num2str(curve_coeffs(21)));
f(x) = subs(f(x),’c8’,num2str(curve_coeffs(24)));
%Divide graph into two different sections and solve for the bounds
[trash index] = max(y);
center_x = values(index,1);
max_x = max(values(:,1));
%Solve for the x-coorinates of half max in both areas
solution_left = vpasolve(f(x) == half_max, [0 center_x]);
solution_right = vpasolve(f(x) == half_max, [center_x max_x]);
28
else
warning(’Input data does not have compatible dimensions.’);
end
fwhm = abs(solution_left-solution_right);
end
29