Spin decoherence due to Fluctuating Fields

Spin decoherence due to Fluctuating Fields

Piotr Szańkowski Institute of Theoretical Physics, University of Warsaw, ul. Hoża 69, PL–00–681 Warszawa, Poland Department of Chemistry, Department of Physics and Department of Electro-Optics, and the Ilse Katz Center for Nano-Science, Ben-Gurion University, Beer-Sheva 84105, Israel    M. Trippenbach Institute of Theoretical Physics, University of Warsaw, ul. Hoża 69, PL–00–681 Warszawa, Poland    Y. B. Band Department of Chemistry, Department of Physics and Department of Electro-Optics, and the Ilse Katz Center for Nano-Science, Ben-Gurion University, Beer-Sheva 84105, Israel
Abstract

The dynamics of a spin in the presence of a deterministic and a fluctuating magnetic field is solved for analytically to obtain the averaged value of the spin as a function of time for various kinds of fluctuations (noise). Specifically, analytic results are obtained for the time dependence of the expectation value of the spin, averaged over fluctuations, for Gaussian white noise, Guassian colored noise, as well as non-Gaussian telegraph noise. Fluctuations cause the decay of the average spin vector (decoherence). For noise with finite temporal correlation time, a deterministic component of the field can suppress decoherence of the spin component along the field. Hence, decoherence can be manipulated by controlling the deterministic magnetic field. A simple universal physical picture emerges which explains the mechanism for the suppression of decay.

pacs:
05.10.Gg, 03.65.Yz, 05.30.-d

I Introduction

Quantum systems are almost invariably coupled to an environment that degrades their coherence. Decoherence of quantum systems is a fundamental problem, with implications across all branches of physics. The temporal evolution of a spin system in the presence of coupling to an environment is perhaps the simplest phenomenon involving decoherence. For example, a particle of spin 1/2, which can be mapped onto any two-level system (TLS) Feynman_63 (), and vice versa, has served as an important paradigm Itakura_03 (); Altshuler_09 (); Cheng_08 (); Aharony_10 (). Specific examples include the decoherence of a qubit in a quantum-computer or quantum-information processor QI (), spin dynamics in a magnetometer Budker_07 (), accuracy and stability limitations of atomic clocks due to interactions with an environment clocks (), decoherence of two-level quantum dot systems QDots (); TaylorQDots (), nitrogen vacancy centers in diamond NV (); NV_fluc ().

Under certain assumptions, which we discuss in the next section, the evolution of a spin in the presence of an environment can be represented by evolving the spin in an effective magnetic field, . Here is a deterministic magnetic field (which can be time-dependent) that can be used to exert control over the spin. The magnetic field models the influence of an environment on the system, and is represented by a vector stochastic process . Averaging over fluctuations corresponds to tracing out the environmental degrees of freedom. This yields a reduced, non-unitary dynamics wherein the averaged spin decoheres in time.

The physical properties of the environment determine the statistical properties of , which in turn determine the type of stochastic process . Here we consider several types of stochastic fluctuations. A prototype model for fluctuations is Gaussian white noise vanKampenBook (); Kloeden (), a random process having vanishing correlation time. As a consequence, the system driven by it evolves without memory, since it has no way to probe the past. A natural generalization is a Gaussian colored noise Valsakuma_85 (); vanKampen_89 (), wherein the random process has a finite correlation time, hence, the system evolves with memory. The ubiquitousness of this type of noise can be attributed to central limit theorem which states that the superposition of a very large collection of independent sources of fluctuation tends to Gaussian noise. Moreover, in some spin systems there are more than one source of noise which differ in their intensity and correlation time scale. For example, in the nitrogen-vacancy diamond system Bar-Gill_12 (), NV centers experience slowly fluctuating noise from C impurity nuclear magnetic moments, and extremely fast fluctuations from nearby electronic nitrogen spins. The former can be modeled by colored noise, which is characterized by a very long correlation time, while the latter can be treated as a white noise, since it sets the fastest time scale in the system. A particular case of colored noise that is not Gaussian, telegraph noise, has been studied in connection with quantum dot qubit systems Aharony_10 (); Altshuler_09 (). Reference Cheng_08 () claims that telegraph noise can be used to model the environment of systems where a spin can interact with relatively few random fluctuators in its neighborhood, and these fluctuators go back and forth between only two states.

The idea of using stochastic methods for modeling an environment or other complex interactions within the spin system is not new. For example, the problem of the line shape in paramagnetic resonance is discussed in Ref. Anderson_53 (). It was assumed that the electrons in an atom are affected by a perturbation which varies randomly in time at a rate determined by exchange interactions. As a result, the fluctuations around the stationary state are induced and their effect on the absorption spectra can be calculated. Another example is presented in Ref. Anderson_62 (), where the authors consider an ensemble of spins interacting with an environment as well as each other. The problem of calculating the effects of these complex interactions on the observed line shape is circumvented by introducing stochastic fluctuations of the precession frequency of the observed spins.

The goal of this work is to develop methods to describe and analyze the process of relaxation of a single spin induced by contact with an environment. An important result of our analysis is that, when the spin evolves with memory, it can be manipulated by applying a sufficiently strong deterministic magnetic field to significantly suppress decoherence. We hope that the level of generality and overall simplicity of our approach allows for a straightforward application to appropriate experimental systems.

The outline of the this paper is as follows. Section II develops the equations of motion for a spin in the presence of a deterministic magnetic field and a stochastically fluctuating magnetic field. Section III describes evolution of the system with Gaussian colored noise fluctuations. In that section we present a simple, intuitive view of the nature of the decay of the averaged spin due to field fluctuations. We also show that this decay can be suppressed when the deterministic magnetic field is strong enough, and the correlation time is long enough (i.e., when the system has memory). Section IV develops the white noise limit of colored noise. Section V explains the effects of memory on decoherence, and Sec. VI considers telegraph noise, which is a specific from of non-Gaussian noise. Finally, a summary and a conclusion are presented in Sec. VII.

Ii Spin dynamics in a fluctuating magnetic field

We consider fluctuations in an otherwise deterministic system that results from the application of a effective random field generated by an environment with which it interacts. The back-coupling of the spin to the source of the noise is not taken into account (Van Kampen vanKampenBook () calls this external noise). Hence, the stochastic properties of the noise result only from the environment which is unaffected by the system. We assume that these properties can be measured or otherwise deduced. For each particular system, one should carefully check whether the back-action can indeed be neglected. Examples of systems for which back-action was neglected include: nitrogen-vacancy centers in diamond affected by magnetic moments of impurity spins Lange (); NV_fluc (), electron spin in a quantum dot affected by nuclear magnetic moments that fluctuate due to crystal lattice vibrations TaylorQDots (); Dobrovitski (), and magnetic noise in atom chips caused by fluctuations of electron currents in wires that make up the atom chip Folman_05 ().

Consider a particle of spin with magnetic moment in the presence of a time-dependent field . The dynamics is determined using a Zeeman-type Hamiltonian,

(1)

where is the vector of spin operators satisfying commutation relations , is the Rabi frequency vector, is a total magnetic field felt by the particle, is the magnetic moment of the particle, is the -factor and is the Bohr (or nuclear) magneton. The magnetic field is the sum of a deterministic field, , whose direction defines the -axis, and a fluctuating field, . Hence, the Rabi frequency is given by

(2)

where and . For simplicity, we have taken the deterministic field to have a component only along the -axis. If the deterministic field varies slowly in time (adiabatically), our conclusions below remain valid, but for deterministic fields that vary more quickly, the simple version of the deterministic part of the evolution matrix that satisfies Eq. (7) requires revision, as detailed below.

In the Heisenberg representation, the evolution of the spin operators is given by the equations,

(3)

where the superscript denotes the Heisenberg picture. In our case, the Hamiltonian contains a stochastic term, :

(4)

The commutator on the right hand side of the Heisenberg equation of motion, Eq. (3), is easy to evaluate, and leads to the vector form of the equation,

(5)

and the initial condition is simply

(6)

The solution to Eq. (5) with the initial condition (6) can be written as a 33 evolution matrix, , acting on the initial vector, . It is convenient to write the equations, , in vector notation as , where the hat indicates a 33 matrix. For convenience, let us take the deterministic field to be constant in time. Substituting into the equations of motion yields,

(7)

Here the matrix elements of are defined as , where is the Levi-Civita symbol. The matrices satisfy commutation relations similar to those of angular momentum operators,

(8)

Consequently, they are generators of rotations in a real, three-dimensional vector space. For example, generates rotations around axis,

(9)

The evolution operator is a rotation operator that causes precession of the spin vector around the instantaneous rotation axis defined by the total field . The equation of motion for the spin does not depend upon the details of the representation of the spin, hence the solution to the dynamics for an arbitrary spin is also given by (7). The only indication of the spin representation, i.e., the dimension of the representation, is in the initial condition.

Experimental measurements of the spin will of necessity correspond to quantum expectation values and averages over the stochastic fluctuations of the magnetic field. Hence the quantity corresponds to experimental measurements of the spin. Here, the symbol means quantum average, and the symbol means average over stochastic variables. Note that the average over fluctuations and the average over the initial quantum state factorizes in the following way: .

For spin 1/2 (i.e., for a TLS), the average expectation value of the spin fully determines the state of the system. This is because the density matrix can be expressed using the Pauli matrices which, together with identity operator , form a basis in the space of 22 hermitian matrices. In the Schrödinger representation,

(10)

and the expectation values of the spin at time is given by vector parameter . Thus, at any time , there are three parameters, , , which completely specify the density matrix. This does not mean that the dynamics of the density matrix for any -level system can be represented in terms of the dynamics of a particle with spin coupled to a magnetic field, since the number of parameters necessary to specify the density matrix larger than the three components of the spin No_indep_param ().

Iii Gaussian Colored Noise

First, let us consider a process that is Gaussian, isotropic and Markovian vanKampenBook (). The field fluctuates in all three dimensions with independent components (isotropy) and statistical properties of each component is completely determined by the first two moments (Gaussian), the average and the correlation function. The stationary condition means that the correlation function depends on the time difference and the average is time independent. According to Doob’s theorem Doob (), if the process is also Markovian, the correlation function is exponential and the average is zero, i.e.,

(11)
(12)

Here the bar indicates average over the fluctuations, is the correlation function with correlation time , the average fluctuating Rabi frequency vanishes, , and the variance of the Rabi frequency, , specifies the strength of the fluctuations.

Several comments regarding Guassian colored noise are in order. First, note that the Gaussian process can model a wide variety of environments due to the central limit theorem; if the environment consists of many individual and independent fluctuating elements, the central limit theorem substantiates the approximation of using a Gaussian process as the effective field affecting the spin. Secondly, if there is no back-action of the spin on the environment, and the environment is in a stationary state (e.g., is in equilibrium) then the process is stationary as well. Finally, the isotropy assumption regarding the random magnetic field produced by the environment is appropriate if no particular special direction can be associated with the environment.

To solve the equation of motion (7), first we eliminate the so called drift term due to the deterministic part of the magnetic field, , from the right hand side of the equation by transforming to the rotating frame. Then Eq. (7) takes the form

(13)

In this equation we introduced transformed matrices, denoted by a prime, and , which are given by

(14)
(15)

where is a deterministic Euclidean rotation matrix that results in counterclockwise rotation around the -axis by the angle [see Eq. (9)]. Note that the scalar product in Eq. (13) can be rewritten as

(16)

The time-dependent rotation matrix “mixes” field fluctuations in the and component. In the Lab frame, and are independent, however, because of the mixing, they are correlated in the rotating frame and the correlation functions are modified in the following way:

(17)
(18)
(19)

The stochastic field in the rotating frame, , remains a Gaussian process, since it is a linear combination of Gaussian processes. Note that any time-independent rotation does not change the correlation functions of the isotropic Gaussian vector process .

The formal solution to Eq. (13) is given in the terms of the time-ordered exponential function,

(20)

defined by the power series, , where is the time-ordering operator vanKampen_Cumulant_Exp (). The average value of , , can be computed by noting that it has the form of a moment generating functional {} for matrix processes, , , where is the vector of trial functions 111The moment generating functional is an alternative means of specifying the probability distributions of stochastic processes. The coefficients in the expansion of in powers of the trial function define the moments of the process.. Then, the average can be “absorbed” inside the exponential function by expressing through a cumulant generating functional {} in terms of cumulant averages denoted here by a double bar (see Ref. Kubo:Cumulants () for a discussion of the cumulant expansion methods used here):

(21)

In the last step, we took advantage of the fact that all cumulants beyond the second vanish for Gaussian processes, and if the average is zero.

Substituting the correlation functions given by Eqs. (17), (18) and (19), we find the following expession for the cumulant generating functional evaluated at ,

(22)

where the rates are given by integrals of the correlation functions:

(23)
(24)
(25)

can be further simplified if we take into account properties of the rotation generators . From the commutation relations (8), the commutator in the third term of Eq. (22) is equal to . It is easy to verify that the sum of squares of and in the second term of can be written as , hence, the parts of Eq. (22) commute with each other. Thus, the average evolution matrix, which is equal to the exponential of , factorizes into three parts,

(26)

where

(27)
(28)
(29)

The first term on the right hand side of Eq. (26) describes the decay caused by the fluctuations in the direction parallel to the applied field . These fluctuations affect only spin components perpendicular to the axis, and the decay rate does not depend on the intensity of the constant field. The second term originates from the fluctuations in the plane, which is perpendicular to the applied field. Fluctuations in the plane couple to all the components of the spin and cause their decay. The key feature is that the rate of decay due to these fluctuations is dependent on . Finally, the third term, also related to the fluctuations in the plane, describes the modification of the precession frequency of the system caused by the noise.

The averaged expectation values of the spin are given by:

(30)
(31)
(32)

In the limit , we obtain and , hence Eqs. (30)-(32) reduce to

(33)

Substituting Eq. (29) into (33) yields the analytic form (33) of the isotropic decay of the average spin, regardless of the initial state of the spin, as can be seen in Fig. 1. For times , the decay is simply exponential with the rate of decay given by . For times , the rate of decay is . In the absence of an external field, the fluctuations in all three directions remain independent and the average evolution matrix decomposes in the following way:

(34)

Here the first part comes from fluctuations in and causes the decay of the and components of the spin, the second part comes from fluctuations in and causes the decay of the perpendicular components, and the third part originates from fluctuations in direction causes the decay of the and components of the spin.

Figure 1: The expectation value of the spin, , subject to Gaussian colored noise, versus time for and . The analytic formula for is given in Eq. (33).

Now let us consider finite . The dependence of , the component parallel to the external field , versus time for three values of dimensional parameter is plotted in Fig. 2. As we discussed previously, component of the average spin decays due to fluctuations in the plane, which are affected by the presence of the constant field. For , the quantity is clearly very similar to the results plotted in Fig. 1, since . As increases, the decay of is significantly slowed, and in the limit of , i.e., the decay is fully suppressed.

Figure 2: subject to Gaussian colored noise versus time for nonvanishing constant (see legend) and . Equation (32) gives the analytic formula for .

Figure 4 plots versus time for the case when and Fig. 4 presents the time dependence of the length of the vector . The decay of the components of the average spin perpendicular to the constant field is not only induced by the fluctuations in the plane, but also by those in the direction [see Eq. (26)]. We have previously seen that the effects of the fluctuations in and are inhibited if is large enough. However, the influence of the fluctuations in the direction parallel to applied field is unperturbed by the presence of . As long as the component of the field is fluctuating the perpendicular component of the average spin will decay to zero for any value of as . For small the decay is, again, similar to the isotropic case shown in Fig. (2), but when increases, the decay is modulated by oscillations which result from the precession of the spin about the external field with frequency modified by the fluctuations. For , the oscillation frequency is given by [see Eq. (28)]. When the frequency of precession tends to as the frequency modification due to fluctuations, vanishes, in a similar manner as that of .

For spin 1/2 system (a TLS), the decay of and correspond to the decay of the coherence of the system, whereas the decay of corresponds to the vanishing of population imbalance. If the constant field is taken along (as opposed to along ), decay of is suppressed, and the decoherence time , which can be made arbitrarily large (suppressed decay) by increasing .

The averaged spin vector solution obtained in this section can be extended to the case of a time dependent deterministic field. If is a function of time then the frequency of oscillations in the and functions appearing in the definitions of the decay rates (24) and (25) are to be replaced by an integral, .

Figure 3: The average value , when the spin is subject to Gaussian colored noise for nonvanishing constant (see legend) and . Equation (30) together with initial condition yield the analytic formula for the function plotted. The average value for initial condition is identical to the results shown here.
Figure 4: Average value of subject to Gaussian colored noise for nonvanishing constant (see legend) and . Equations (30) and (31) yield the analytic formula for the function plotted.

Iv The White Noise Limit

The white noise limit of colored noise involves going to the regime of extremely rapid fluctuations of the field, i.e., vanishing correlation time . Unfortunately, simply taking does not yield the correct limit. If we put in Eq. (29), we obtain that the decay rate of the average spin equals zero, and fluctuations do not affect the evolution of the spin. The proper limit can be obtained if the strength of fluctuations, as given by the variance of the field, , is increased as the correlation time goes to zero, so that remains constant. Thus, we arrive at the limit when the correlation function of the colored noise becomes proportional to a Dirac delta function,

(35)

In the limit of the decay rates and . The decay of the expectation value of the spins are then isotropic and are not affected by the presence of a constant field , i.e.,

(36)

Thus, the spin decay is isotropic, purely exponential with decay rate , and is independent of . Moreover, the precession due to is unaffected by the fluctuations. In the white noise limit, the precession due to deterministic part of the field and the decoherence due to field fluctuations are completely decoupled.

V Memory effects

In the expression for for the case of Gaussian colored noise, the decay rates of the spin depend on the history of the fluctuating field through the integrals of correlation functions. But the decay rate for the white noise case does not depend on history. Based on this observation, we say that spin system has memory if the fluctuating field driving it has finite correlation time. This definition is consistent with a rigorous notion of memory, since , understood as a stochastic process, is Markovian vanKampenBook () (i.e. memoryless) if and only if the noise driving it is white noise. Memory in the system has important consequences, because it allows control of spin decay, through the application of an external field.

The role of memory can be understood as follows. Consider first the case of . Each realization of the fluctuating field determines corresponding realization of via the equation of motion, where defines an instantaneous rotation axis for the spin vector. Hence, each realization draws a trajectory of the spin on a sphere of radius . The evolution can be viewed in a simplified “stroboscopic” picture in which, within a correlation time , the axis does not change much, and subsequently, it jumps to a new random value within a range on the order of . Eventually such a sequence of rotations around randomly selected axes “smears” the trajectories of the spin over the whole sphere. This results in the decay of the spin averaged over all realizations.

The situation is different when a constant field along, say, the -axis is present. This becomes clear when we view the evolution in the reference frame rotating with frequency around the -axis. Consider a particular realization of the fluctuating field , and denote it in the rotating frame as . Since , the fluctuations in the direction parallel to the constant field are unchanged when shifted to the rotating frame. On the other hand, the projection of the fluctuating field on the plane perpendicular to -axis, , is modified when viewed in the rotating frame. Within a time of order the correlation time , the projection revolves around the -axis with frequency , and if , is able to perform many revolutions between jumps. The effect of the jumps becomes negligible in comparison with the systematic rotation. Thus the stochastic character of is lost, and the only contribution to the decay caused by fluctuations comes from . A similar effect provides stable motion of a spinning top, and for that reason we will refer to it as a gyroscopic effect. If , then revolution in between jumps is negligible, so , and by necessity the decay induced by the fluctuations must be almost the same as if there was no constant field at all.

The suppression mechanism of the decay described above is inoperative in the white noise limit where the correlation time vanishes, so the fluctuations of the rotation axis are instantaneous. Hence, the decay caused by white noise cannot be affected by a constant field of finite intensity. Moreover, even if the strength of the fluctuations is small in comparison to , the effect of fluctuations cannot be neglected if the correlation time is too short.

Note that the important dimensionless parameter here is , which can be interpreted as an average angle of precession due to the deterministic field within the time period of . Bearing this in mind, it is clear that the conclusion of this paragraph remains valid even if the field varies in time as long as . Moreover, the reasoning used above shows that the same decay suppression effects apply even for non-Gaussian noise as long as it has non-vanishing correlation time. We verify this assessment in the following section.

Vi Non-Gaussian process: telegraph noise

In this section we present an example of a non-Gaussian process, telegraph noise, i.e., dichotomic noise corresponding to a memoryless continuous-time stochastic process that jumps between two distinct values, sometimes called burst noise vanKampenBook (); GardinerBook (). We consider telegraph noise in one component of the field, . The Rabi frequency of the fluctuating field, , jumps randomly between and , with a jump rate . The average value and the correlation function of such a process is given by

(37)

Because telegraph noise is non-Gaussian, the correlation function does not contain all the information regarding the process.

Telegraph noise can model the effect of a spin 1/2 impurity atom near the spin Itakura_03 (); Cheng_08 (); Aharony_10 (). The impurity spin has probabilities and for being in the spin-up and spin-down state, and a hopping rate for transferring from one spin state to the other. Neglecting the back-action of the magnetic moment on the impurity, this model reduces to telegraph noise.

Our goal is to find the average evolution matrix for the telegraph noise case. The technique we have used previously to find the average does not work in the case of telegraph noise, because unlike with Gaussian processes, the closed form of the cumulant generating functional is not known. Fortunately, it is still possible to find an exact, analytical formula for if the fluctuations are parallel to the constant field. One way of solving for is to follow Ref. VanKampen_Tele (). We can write the average evolution operator as , where are the averages conditional on being equal to . These averages vary in time because of equation of motion (7), and because jumps. Thus,

(38)
(39)

Adding and subtracting these equations, we obtain

(40)
(41)

where . The next step is to eliminate from Eq. (40) by substituting the solution to Eq. (41), with the initial condition . This yields a closed equation for average evolution matrix:

(42)

Note the particular form of the “memory integral” in the equation (42), where the dynamical variable is evaluated at times previous to . In this case the system clearly has memory and it remembers not only the history of the noise but also its own history. It is easy to check that the solution to this equation that satisfies the initial condition is given by

(43)

where

(44)

The average expectation values of the spin are obtained by applying the average evolution matrix (43) to the initial vector of spin operators, , and this yields,