ABSTRACT
SLOW AND STOPPED LIGHT WITH MANY ATOMS, THE ANISOTROPIC RABI MODEL AND A PHOTON COUNTING EXPERIMENT ON A DISSIPATIVE OPTICAL LATTICE
by Tyler Thurtell
First, we study electromagnetically induced transparency (EIT) in the case that interactions between the atoms are allowed. Under the right circum- stances EIT can lead to a dramatic slowing or even stopping of a pulse. We consider EIT in the case of incoherent collective emission. We discuss how collective emission can be used to enhance the similarity between the input and output pulse in stopped light. Next, we introduce coherent, ’spin flip’ interactions that allow the atoms to ’trade’ excitations. For coupling on the probe transition we find that the minimum speed to which the light can slowed depends on how homogeneous the behavior of the atoms is. In the case of probe transition coupling, we find that EIT may not even be possible and if it is a frequency shift of the atomic resonance will occur. We then move on to discuss some generalizations of the Jaynes-Cummings model. The Jaynes-Cummings model is actually a simplification of the Rabi model. We review these two models and then discuss the anisotropic Rabi model in which co-rotating and counter-rotating terms may have distinct coupling strengths. We find eigenstates of the interaction Hamiltonian with zero eigenvalue for the two couplings not equal and investigate the behavior as the couplings become more similar. Finally, we discuss the beginnings of a photon counting experiment on a dissipative optical lattice. The photon counting apparatus will be used to study the well depths of an optical lat- tice. Unfortunately, due to time constraints, we were unable to complete the experiment and so we will discuss only the experimental equipment and data analysis techniques. SLOW AND STOPPED LIGHT WITH MANY ATOMS, THE ANISOTROPIC RABI MODEL AND A PHOTON COUNTING EXPERIMENT ON A DISSIPATIVE OPTICAL LATTICE
A Thesis
Submitted to the Faculty of Miami University in partial fulfillment of the requirements for the degree of Master of Science Department of Physics by Tyler Thurtell Miami University Oxford, Ohio 2018
Advisor: (Perry Rice)
Reader: (Samir Bali)
Reader: (Carlo Samson) Contents
List of Figures v
1 Introduction 1
2 Open Quantum Systems 6
3 Numerical Techniques 19
4 Electromagnetically Induced Transparency 24 4.1 Fundamental Processes Approach ...... 24 4.2 Adiabatic Elimination ...... 27 4.3 Dark States ...... 28 4.4 Electromagnetically Induced Transparency ...... 31 4.5 Slow Light ...... 34 4.6 Stopped Light ...... 35
5 Superradiance, Subradiance, and Selective Radiance 37 5.1 Standard Formalism ...... 39 5.2 Green’s Function Approach ...... 44
6 EIT and Slow Light with Spin Flip Interactions 50 6.1 Control Transition Interactions ...... 50 6.2 Probe Transition Interactions ...... 53
ii 7 The Jaynes-Cummings Model and its Generalizations 56 7.1 Cavity QED Basics ...... 56 7.2 The Jaynes-Cummings Model ...... 58 7.3 The Rabi Model ...... 60 7.4 The Anisotropic Rabi Model ...... 60
8 Photon Counting in a Dissipative Optical Lattice 66 8.1 Theory of Optical Lattices ...... 67 8.2 Numerical Simulations ...... 70 8.3 Experimental Techniques ...... 76 8.4 Data Analysis ...... 79
9 Summary 85
A Properties of the Density Operator 88
B Alternative Derivation of the Master Equation 90
C Guide to the Use and Abuse of QuTiP 95
D Green’s Functions 109
E Pulse Propagation Programs 112
F Jaynes-Cummings Model Generalization Program 122
G Semi-Classical Diffusion Programs 125
H Photon Counting Experimental Apparatus Details 134 H.1 Single Photon Counting Module Details ...... 134 H.2 FPGA Multichannel Acquisition Board ...... 135
I Proof of the Cross Correlation Theorem 138
iii J Photon Counting Data Analysis Programs 140 J.1 Cross Correlation Theorem Program ...... 140
K Bibliography 153
iv List of Figures
3.1 A flow chart indicating the algorithm at work in quantum trajectory theory...... 22
4.1 A lambda system with three decay channels interacting with two near resonant laser beams...... 25 4.2 A diagrammatic representation of the sum over processes in- volved in EIT...... 26 4.3 (a) A plot of Re[α(ω)] versus δ/Γ for an EIT atom. (b) A plot of Im[α(ω)] versus ∆/Γ for an EIT atom. Both (a) and (b) are for the case of Raman resonance. (c) A plot of Re[α(ω)] versus δ/Γ as predicted by Lorentz oscillator theory for comparison with (a). (d) A plot of Im[α(ω)] versus δ/Γ as predicted by Lorentz oscillator theory for comparison with (b)...... 33
5.1 Heat plots of the probe Rabi frequency for a pulse propagat- ing through a gas of EIT atoms with time on the x-axis and 0 distance on the z-axis for the case Γ1 = Γ = 0. From left to right the storage times are 0.6T , 1.2T , 1.6T , 2.0T where
T = 1/Γr for the top row. (a) The case where Γr = 10Γ3. (b)
The case enhanced by selective radiance where Γr = 20Γ3. . . 38
v 5.2 The relationship of the spin and the magnetic field at (a) the beginning of the radiative process, (b) halfway through the radiative process when the intensity is the greatest, and (c) at the end of the process...... 43 5.3 Illustration of atoms coupled to a tapered optical nanofiber. . 48
7.1 An illustration of the situation studied in Cavity QED. . . . . 57 7.2 The atomic inversion vs time in the Jaynes-Cummings model. 59 7.3 The atomic inversion (blue curves) and photon number (green curves) vs time in the Rabi mode. In (a) the atom is initially in the excited state and the initial photon number is 3. In (b) the atom is initially in the ground state and the initial photon number is 0...... 61 7.4 The atomic inversion (blue curves) and photon(green curves)
number vs time in the anisotropic Rabi model for g1 = ω = ωeg
with (a) g2 = 0.01g1, (b) g2 = 0.1g1, (c) g2 = 0.5g1, (d)
g2 = 2g1, (e) g2 = 5g1 and (f) g2 = 10g1. Notice that for g2 ≈ 0
Jaynes-Cummings dynamics are recovered, for g2 = g1 the
Rabi model dynamics recovered, and for g2 >> g1 a distinct beating appears...... 63
8.1 Illustration of the lattice created by two linearly polarized beams perpendicular to each other and propagating in oppo- site directions. The result is periodic net right and left circu- lar polarization. In between the lattice sites the polarization varies continuously between these two extremes. The atoms are trapped in the locations of pure circular polarization and oscillate about these points as shown, with a characteristic
vibrational frequency denoted by ωv...... 67
8.2 Energy level diagram of a six level,Fg = 1/2 → Fe = 3/2, atom. 67
vi 8.3 Illustration of the five possible processes for an atom in the +1/2 ground state. A blue arrow signifies that the process
involves a right circularly polarized photon, σ+, red arrow indicates a linearly polarized photon, π, and a green arrow
indicates a left circularly polarized photon, σ−. (a) A σ+ is
absorbed, another σ+ is emitted and the atom remains in the same state. (b) A π is absorbed, another π is emitted and the
atom stays in the same state. (c) A π is absorbed, a σ− is
emitted and the atom changes state. (d) A σ− is absorbed, a
π is emitted and the atom changes state. (e) A σ− is absorbed,
a σ− is emitted and the atom changes state...... 68 8.4 Examples of position vs. time and momentum vs. time be- havior in an optical lattice as predicted by the semi-classical algorithm for an atom beginning at rest at origin in the -1/2 ¯h2k2 state with (a) U0 = 20ER and (b) U0 = ER where ER = m2 is the photon recoil energy...... 73 8.5 Examples of position vs. time and momentum vs. time be- havior in an optical lattice as predicted by the modified semi- classical algorithm for an atom beginning at rest at origin in
the -1/2 state with (a) U0 = 20ER and (b) U0 = ER where ¯h2k2 ER = m2 is the photon recoil energy...... 75 (2) 8.6 An example of the σ+ − σ+ g (τ) predicted by the modified semi-classical algorithm...... 76 8.7 Example of position vs. time and momentum vs. time be- havior in an optical lattice as predicted by the modified semi- classical algorithm for an atom beginning at rest at origin in
the -1/2 state with U0 = ER allowed to run for long enough for the atom the change lattice sites many times...... 77
vii 8.8 An example of position vs time and momentum vs time be- havior with a constant applied force as simulated based on the modified semi-classical algorithm...... 78 8.9 A picture of the single photon counting module. Connected to the left side of the detector is multimode optical fiber. Con- nected to the right side is the 5V DC power supply and two BNCs. One of the BNC carries the TTL output pulse when a photon is detected. The other BNC is a gating input. When the gate has a TTL high pulse it should turn the detector off. 79 8.10 A picture of the (a) outside and (b) inside of the FPGA mod- ule. The BNC ports on the front of the FPGA module are well labeled. The output from the photon counting module should be connected to one of the ports labeled detector. The other important port is the one labeled external clock. It may be used to have the FPGA module keep track of shorter timescales. 80 8.11 The averaged white light data. It should show a flat line at g(2)(τ) = 1 but it shows a distinct decay...... 81 8.12 The unaveraged white light data for (a) 6446 counts and (b) 64460 counts. In case (a) th value is clearly elevated above one. In case (b) it is not so obviously above one but it most likely is...... 82 8.13 The unaveraged white light data for (a) 6446 counts and (b) 64460 counts as calculated by the shift register program. Again, in case (a) th value is clearly elevated above one. In case (b) it is not so obviously above one but it most likely is...... 84
C.1 The atomic inversion vs time in the Jaynes-Cummings model. 99 C.2 The atomic inversion vs time in the Jaynes-Cummings model for an initial state that is mixed...... 101 C.3 The atomic inversion vs time in the damped Jaynes-Cummings model for an initial state that is mixed...... 102
viii C.4 The atomic inversion vs time in the damped Jaynes-Cummings model for a coupling with a sinusoidal time dependence for an initial state that is mixed...... 103 C.5 The atomic inversion vs time in the damped Jaynes-Cummings model for a coupling with a sinusoidal time dependence for an initial state that is mixed produced using quantum trajec- tories. In (a) only one trajectory was used and in (b) 100 trajectories were averaged over...... 104 C.6 The steady state atomic inversion vs coupling strength in the driven Jaynes-Cummings model...... 106 C.7 The steady state population of level |3i vs detuning for a lambda system showing the famous EIT dip...... 108
H.1 The light tighting apparatus from the (a) front and (b) back. The opening in the front is covered with lens tissue and marked with pencil. Even with the light tighting the flashlight was not shown directly at the apparatus but off at an angle. The laser light was shown at the very edge of the opening with a current of 90mA...... 136
ix ACKNOWLEDGEMENTS
Foremost, I would like to thank my thesis advisor, Dr. Perry Rice for all the important lessons he has taught me and his constant encouragement. I must also acknowledge the many great teachers I have add, including my other research advisor Dr. James Clemens and committee members Dr. Samir Bali and Dr. Carlo Samson for teaching me the methods and joys of quantum mechanics. I also must thank Dr. Xiaodong Qi and Dr. Ivan Deutsch from the University of New Mexico who helpfully supplied the code used to solve for the decay rates and frequency shifts of atoms near an optical nanofiber. I am also indebted to my friends and fellow students, including but not limited to: Nazar Al-Aayedi, Zeeshan Ali, Dharma Raj Basaula, Subhash Bhatt, Ben Blankartz, Patrick Carroll, Lyndon Cayayan, Jijun Chen, Ken DeRose, Ajithamithra Dharmasiri, Billy Drake, Arkan Hassan, Martin Heidelman, Patrick Janovick, Daniel King, Eitan Lees, Mitch Mazzei, Anthony Rapp, Jayson Rook, Micheal Saaranen, Joshua Schussler, Dinesh Wagle, Anthony Young, and Sara Zanfardino for many helpful discussions and for making the past two years an enjoyable experience. Above all, I must thank my family and loved ones, especially, my parents and Elizabeth, without whom I would never have achieved this.
x Chapter 1
Introduction
First, we will examine electromagnetically induced transparency (EIT)[1][2][3][4][5][6][7] in the case that additional interactions are allowed between the atoms. We will examine EIT in the case that ’spin flip’ interactions are allowed and in the case that the atoms radiate collectively. By ’spin flip’ interactions we essentially mean interactions in which an excitation in one atom can become an excitation in another atom with negligible loss. When atoms radiate col- lectively they may effectively radiate more quickly or more slowly than an isolated atom would. These questions are of particular interest because of their applications to quantum information processing. Many quantum information processing devices will require some physical hardware to function as a memory. However, this is difficult since measuring a single quantum state does not reveal the exact state due to the proba- bilistic nature of measurement in quantum mechanics. In other words, this is a challenge because quantum tomography, the process of determining the quantum state based on repeated measurements on an ensemble of identi- cally prepared states, is an open and difficult problem[8] when only a finite number of samples is available. Even if good procedures for quantum state determination were available it could take many classical bits to record a single qubit of information. The alternative solution is to develop hardware
1 that can store a qubit directly. One physical realization of a qubit is the polarization state of a photon. A memory that can store a photon state is call a photonic memory. One proposed photonic memory is based on the phenomenon of stopped light or light storage[4]. In this phenomenon the information about the state of the photon is transfered to the matter where it is held until some later retrieval time. Specifically, it is transfered to atoms exhibiting EIT. EIT is a non- perturbative effect in which the presence of a strong control laser beam can make a medium transparent to a weak probe beam it was previously opaque to. This effect occurs in three level atoms and is a result of coherent interfer- ence between the levels. Associated with this effect is a rapid change of the index of refraction of the medium with frequency and therefore a low group velocity for the light propagating through the medium. This is referred to as slow light[4][6][7][9]. If the control beam is manipulated properly, the group velocity of the light can be reduced to zero resulting in stopped light[4][6][10]. One problem with this photonic memory is that the fidelity of the stor- age is limited by the spontaneous emission rate of the atoms[11]. However, while the spontaneous emission rate of isolated atoms is determined solely by properties of the atom the spontaneous emission rate of a collection can be modified by coherent interference between the atoms. When this inter- ference increases the spontaneous emission rate, it is referred to as suppradi- ance and when it decreases the spontaneous emission rate it is referred to as subradiance[11]. The phenomena of superradiance was first investigated by Dicke in 1954[12]. Since then it has been extensively studied[6][13][14]. These effects usually only play a role when the atomic separation is less than the wavelength of the emitted light[13]. However, recent advances in nanotech- nology have made nanophotonic waveguides available. Since the field inside the waveguide does not decay, or decays much more slowly, the atoms can un- dergo the same interference effects when separated by distances significantly greater than a wavelength. For example, this effect has been used for appli-
2 cations in quantum metrology[15][16]. This suggests that when the modes of a nanophotonic waveguide are properly engineered, three level atoms around the waveguide can serve as a higher fidelity photonic memory[11]. However, consideration of this problem also leads one to consider the problem of EIT with spin flip interactions since any situation in which the atoms may radiate collectively is likely to be a situation in which they may trade excitations because the atoms will likely either be close together or be linked by some dielectric. Rather than examine particular atomic geometries we will introduce an effective coupling between the atoms, the form of which would have to be derived from first principles for any particular circumstance. Next, we study the anisotropic Rabi model which may be viewed as a gen- eralization of the Jaynes-Cummings model. The Jaynes-Cummings model is the central model of cavity quantum electrodynamics (QED)[7][17][18][19]. It models the interaction of a two level atom with a single field mode. This model is quite general and emerges outside of atomic physics in a number of situations. In atomic physics, it actually emerges as a simplification of the Rabi model which includes so called non-energy conserving terms[20][21]. The Rabi model contains some new phenomena. For example, as the cou- pling strength is tuned, this system exhibits a phase transition[22]. We re- view these models and then discuss the anisotropic Rabi model in which the relative coupling of the energy conserving and energy non-conserving terms becomes arbitrary[23]. This model has been physically realized, for example in strongly coupled superconducting circuits[24]. In particular, in the in- teraction picture, we find Hamiltonian eigenstates with zero eigenvalues for non-equal couplings. Chapter 8 is conceptually separate from the rest of this manuscript. It discusses an experiment preformed under the supervision of Dr. Samir Bali on photon counting in a dissipative optical lattice. Optical lattices, which allow neutral atoms, usually alkalis, to be trapped at temperatures below the Doppler limit in specific locations known as lattice sites, were first discussed
3 by C. Cohen-Tannoudji in 1989[25]. Optical lattices are interesting in their own right but have also found wide applications, for example in quantum information[26]. Dr. Bali’s group is interested in using the optical lattice to produce thermally directed motion of the atoms. That is, to create a ratchet. Such cold atom ratchets have been studied previously[27]. Dr. Bali’s group hopes to use the ratchet to simulate biomolecular motors in an effort to understand the origin of their surprisingly high efficiency[28]. Correlations between the types of photons emitted from the lattice carries important information about the lattice[29]. To understand how this works think of each lattice site as a potential well the atom may become trapped in. When trapped in some of the wells the atom will emit only right circularly polarized, σ+, light. When the atom is trapped in an adjacent site, it will radiate only left circularly polarized, σ−, light. So that if the observation of a
σ+ photon is only correlated with the observation of a σ+ photon some ’long’ time later it indicates that the atoms are well trapped. On the other hand, if the observation of a σ+ is correlated with the observation of a σ− photon a short time later this indicates that the atoms are not well trapped. These types of experiments have been preformed before but many details remain to be worked out[29][30][31]. For example, the behavior in the limit of shallow wells has not been extensively studied. These ideas will be explored further and the experimental techniques will be discussed in Chapter 8. The layout of the manuscript is as follows: In Chapter 2 the density opera- tor is introduced, the master equation is derived in a manner that emphasizes tracing over the reservoir degrees of freedom, and quantum trajectory theory is discussed. In Chapter 3 the numerical techniques used throughout are briefly discussed with a special emphasis on the use of the Quantum Tool- box in Python (QuTiP)[32][33], although readers interested in the details of any particular simulation will have to refer to the appropriate appendix. Chapters 4 though 6 discuss EIT in the situations mentioned above. Chap- ter 7 discusses the generalizations of the Jaynes-Cummings model. Finally,
4 chapter 8 discusses photon counting in a dissipative optical lattice.
5 Chapter 2
Open Quantum Systems
In the normal formulation of quantum mechanics[34], the state of a system is described by an element of a Hilbert space and is denoted by
|ψi. (2.1)
These elements are referred to as kets. Properties of the Hilbert space imme- diately provide an alternative representation as elements of the dual vector space referred to as bras. These are denoted by hψ|
hψ| = |ψi†, (2.2) where the dagger indicates Hermitian conjugation. Except for in the case of measurement, the time evolution of the kets is unitary and given by the Schrodinger equation i |ψ˙ i = H|ψi, (2.3) h¯ whereh ¯ is Planck’s constant divided by 2π and H is the Hamiltonian of the system. The time evolution of the bra representation is given by the Hermitian conjugate of this equation. Unfortunately, these representations have some issues. Notably, they cannot handle statistical uncertainty in the
6 state. As we shall see this also means they cannot handle open quantum systems[19][35]. If we have perfect knowledge of a quantum system we assign it a state exactly as above. The state is then referred to as a pure state. If we have imperfect knowledge of a quantum system we require a new way to assign it a state that reflects this lack of information. Let’s review how this is done in classical mechanics[36]. If we have perfect knowledge of a classical system, its canonical coordinates and momenta are assigned exact values and their time evolution is given by Hamilton’s equations
q˙i = (qi,H) (2.4)
p˙i = (pi,H), (2.5) where the qi are the canonical coordinates, the pi are the canonical momenta, and the () are the Poisson brackets. If we have imperfect knowledge of the system it is described by Liouville’s equation of motion
∂ D = (H,D), (2.6) ∂t where D is a probability distribution on phase space. There is a representation is called the density operator, which can be appropriately generalized[34]. It is given by
X ρ = Pi|ψiihψi|, (2.7) i where Pi is the probability that the system is in the state represented by
|ψii. To understand why this is a good generalization of the probability den- sity consider the function of a probability density in calculating expectation values Z hOi = DOdV, (2.8) V
7 where O is an observable and V is the volume of phase space. For a pure quantum state, expectation values are calculated according to
hOi = hψ|O|ψi. (2.9)
However, since the trace of a scalar is simply the original scalar this may be written as
T r (hOi) = T r (hψ|O|ψi) = T r (|ψihψ|O) = T r(ρO), (2.10) where the cyclic property of the trace has been used. This equation is of the same form as the classical equation if the integral generalizes to the trace and D generalizes to ρ. This calculation has also shown how to generalize the calculation of expectation values from pure states to mixed states. The time evolution of the ρ follows from the Schrodinger equation
X h i i X ρ˙ = P |ψ˙ i hψ| + |ψi hψ˙ | = − P (H|ψihψ| − |ψihψ|H) . i h¯ i (2.11) This may be written compactly as
i ρ˙ = − [H, ρ]. (2.12) h¯
This equation is usually referred to simply as the Schrodinger equation but it is sometimes called the von Neumann equation or the Liouville equation. Notice, there is an important physical distinction between coherent su- perpositions and statistical mixtures[7][34]. The former is a fundamental property of quantum mechanics. It is responsible for the famous double slit experiment and electromagnetically induced transparency which will be dis- cussed below. The latter is not fundamental but rather arises out of our ignorance about the state of the system. The two kinds of uncertainty can also lead to different results when calculating expectation values. Consider
8 the following two states in the density operator representation
1 ρ = (|1ih1| + |1ih2| + |2ih1| + |2ih2|) (2.13) sup 2 1 ρ = (|1ih1| + |2ih2|) . (2.14) mix 2 The first is a coherent superposition of the states |1i and |2i while the latter is a statistical mixture. They both give the same probabilities for measuring the system to be in states |1i or |2i but they do not give the same expectation value for all observables. As an example, consider the observable
O = |1ih2| + |2ih1|. (2.15)
The expectation values of this observable for the two states above are
hOisup = T r(ρsupO) = 1 (2.16)
hOimix = T r(ρmixO) = 0. (2.17)
Some important properties of the density operator are discussed in more detail in Appendix A. This discussion is important for our purposes because open quantum sys- tems will naturally become mixed even if they are initially pure. To un- derstand why this happens consider a classical system of interest coupled to an environment that has many more degrees of freedom than the classical system. Clearly the environment can affect the system but any energy that flows from the system to the environment will be spread out among all of the many degrees of freedom. This means two things. First, the environment will not be appreciably affected by the system. Second, energy and information that flow into the environment will not flow back to the system. A system that has this property is described as Markovian[35]. As concrete example consider an ice cube placed in the ocean. The ocean will cause the ice cube
9 to melt but the ice cube will not appreciably decrease the temperature of the ocean. The ice cube will also never spontaneously reassemble itself in part or in whole. Clearly, if we only keep track of the system degrees of freedom in- formation about the state of the system will be lost in time. In the quantum generalization, the system becomes entangled with the environment. Since the environment degrees of freedom are not kept track of, this leads to a loss of information and system initially in a pure state evolving to a mixed state. The appropriate generalization of the Schrodinger equation was worked out in the early 1960s[37][38]. To see how it arises we follow a well trodden path[39][40], consider the combined state of a system of interest and a reser- voir that the system interacts with described by a density operator ρ and assume that initially the system and the reservoir are not entangled so we may write
ρ = ρs ⊗ R, (2.18) where ρs depends only on the state of the system and
X R = Pn|nihn|, (2.19) n describes the state of the reservoir. In order to find a description that only references the state of the system we trace over the reservoir degrees of free- dom. The system density operator at a time t is then related to the initial system density operator by
† X X p † p ρs(t) = T rR(Uρs ⊗ Pn|nihn|U ) = Pnhµ|U|niρs(0)hn|U |µi Pn, µ n (2.20) where the |µi span the reservoir section of the Hilbert space. Defining the Kraus[41] operators to be
Mµ,n(t) = Mν(t) = hµ|U|ni, (2.21)
10 where the first equal sign is a relabeling. The system density operator at time t is then given by
X † ρs(t) = Mν(t)ρsMν (t) ≡ a(t)[ρs(0)] (2.22) ν
This map is referred to as completely positive since it will always map a positive operator to a positive operator even if there are other subsystems present on which it must act as the identity. An important property of this map is that it is trace preserving.
! ! X † X † X T rs (ρs(t)) = T rs Mµ,nρs(0)Mµ,n = T rs Pnhn|U |µihµ|U|niρs(0) = T rs (ρs(0)) . µ,n n µ (2.23) It is important to note that this property held true because
X † Mν Mν = 1s. (2.24) ν
This is sometimes referred to as the completeness property of the Kraus operators. One subtlety here is that the map is only trace preserving if the |ni and the |µi are capable of forming resolutions of the identity in some section of the Hilbert space. To first order in dt we may write one of the Kraus operators as
M0 = 1s + G(t)dt, (2.25)
and define a new set of operators Lν(t) by √ Mν(t) = Lν(t) dt. (2.26)
11 The completeness property of the Kraus operators implies
† X † M0 M0 = 1s − Mν Mν. (2.27) ν=1
To first order in dt this implies