
Quadratic Zeeman Spectral Diffusion of Thulium Ion Population in a Yttrium
Gallium Garnet Crystal
Jacob H. Davidson,1, ∗Antariksha Das,1, †Nir Alfasi,1Rufus L. Cone,2Charles W. Thiel,2and Wolfgang Tittel1, 3, 4
1QuTech and Kavli Institute of Nanoscience, Delft University of Technology, 2628 CJ Delft, The Netherlands
2Department of Physics, Montana State University, Bozeman, Montana 59717, USA
3Department of Applied Physics, University of Geneva, 1211 Geneva 4, Switzerland
4Schaffhausen Institute of Technology - SIT, 1211 Geneva 4, Switzerland
(Dated: October 12, 2022)
The creation of well understood structures using spectral hole burning is an important task in the
use of technologies based on rare earth ion doped crystals. We apply a series of different techniques
to model and improve the frequency dependent population change in the atomic level structure of
Thulium Yttrium Gallium Garnet (Tm:YGG). In particular we demonstrate that at zero applied
magnetic field, numerical solutions to frequency dependent three-level rate equations show good
agreement with spectral hole burning results. This allows predicting spectral structures given a
specific hole burning sequence, the underpinning spectroscopic material properties, and the relevant
laser parameters. This enables us to largely eliminate power dependent hole broadening through the
use of adiabatic hole-burning pulses. Though this system of rate equations shows good agreement
at zero field, the addition of a magnetic field results in unexpected spectral diffusion proportional to
the induced Tm ion magnetic dipole moment and average magnetic field strength, which, through
the quadratic Zeeman effect, dominates the optical spectrum over long time scales. Our results
allow optimization of the preparation process for spectral structures in a large variety of rare earth
ion doped materials for quantum memories and other applications.
I. INTRODUCTION
Rare-earth ion doped crystals (REICs) are interesting
materials due to their long-lived excited states and their
exceptionally long optical coherence times at cryogenic
temperature [1, 2]. In particular, along with the pos-
sibility for spectral tailoring of their inhomogeneously
broadened 4fN-4fNtransitions, this makes them prime
candidates for a number of applications in classical and
quantum optics. Examples include laser stabilization,
RF spectrum analysis, narrow band spectral filtering,
and quantum information storage and processing [2–7].
Thulium-doped Yttrium Gallium Garnet (Y3Ga5O12,
Tm:YGG) is one such material. Its 3H6↔3H4transition
at 795 nm wavelength features an optical coherence time
of more than 1 ms [8–10], which is one of the longest
among all studied REICs. In combination with the ac-
cessibility of this transition—within the range of com-
mercial diode lasers—this makes it a natural candidate
for applications.
The quality of created features and the resulting con-
sequences for associated applications, are dependent on
the spectroscopic properties of the dopant ions and their
numerous interactions with other atomic components in
their local crystalline environment [11, 12], the details of
the optical pumping process, and the spectral and tempo-
ral profile of the applied laser pulses[13, 14]. Deep under-
standing of the relation between spectroscopic properties,
∗These authors contributed equally to this work. Present Address:
National Institute of Standards and Technology (NIST), Boulder,
Colorado 80305, USA
†These authors contributed equally to this work.
optical control fields, and spectral diffusion dynamics has
resulted in improvements of this process in a number
of other rare-earth-doped materials including Tm:YAG,
Eu:YSO, and Pr:YSO [15–17]. However, this important
connection has thus far not been made for Tm:YGG.
In this paper we track the evolution of population
within the electronic levels of Tm3+ ions in YGG (see
Fig. 2 for simplified level scheme) by semi-continuous
monitoring of spectral holes for many sequences of ap-
plied spectral hole burning pulses. The characteristic
shapes and sizes of these spectral features are matched
to a rate equation model that encompasses the ground
(3H6), excited ( 3H4), and bottleneck (3F4) levels in this
material with associated lifetimes and branching ratios.
At zero magnetic field we see good agreement between
our numerical model and measured results across many
different pump sequences of varying duration, power, and
spectral shape. With the addition of an external mag-
netic field the agreement with our numerical model dis-
appears as spectral diffusion from local host spins begins
to dominate the shape of all spectral features over long
timescales. We characterize the nature of this unexpected
behavior and expand our model accordingly by adding a
spectral diffusion term to account for a quadratic Zeeman
interaction with present noisy magnetic fields [18, 19].
The letter is structured as follows: In section II we
describe the experimental setup used to collect our mea-
surements. In section III we detail the atomic level struc-
ture in Tm:YGG and introduce spectral hole burning, the
workhorse of our investigations, to select a known set of
atomic population. In section IV we introduce and apply
a rate equation model which shows good agreement to the
measured spectral hole features. In section IV A we detail
the use of adiabatic pulses to shape spectral holes at zero
arXiv:2210.05005v1 [quant-ph] 10 Oct 2022