# Enhancement of laser cooling by the use of magnetic gradients

###### Abstract

We present a laser cooling scheme for trapped ions and atoms using a combination of laser couplings and a magnetic gradient field. In a Schrieffer-Wolff transformed picture, this setup cancels the carrier and blue sideband terms completely (up to first order in the Lamb-Dicke parameter), resulting in an improved cooling behaviour compared to standard cooling schemes in the Lamb-Dicke regime (e.g. sideband cooling) and allowing cooling to the vibrational ground state. A condition for optimal cooling rates is presented and the cooling behaviour for different Lamb-Dicke parameters and spontaneous decay rates is discussed. Cooling rates of one order of magnitude less than the trapping frequency are achieved using the new cooling method. Furthermore the scheme exhibits fast rates and low final populations even for significant deviations from the optimal parameters and moreover provides good cooling rates also in the multi-particle case.

###### pacs:

## I Introduction

Cooling of atoms/ions is an important precondition for quantum logic and computation Wineland et al. (2005); Sasura and Buzek (2002), wherein an initial ground state preparation of the vibrational state is required, for high resolution spectroscopy beyond the Doppler limitations and for the analysis of low temperature phenomena (e.g. BEC condensates). Cooling by light (laser cooling) for atoms and ions was originally proposed in Hänsch and Schawlow (1975); Wineland and Dehmelt (1975). For trapped ions, sideband cooling Stenholm (1986); Wineland et al. (1978) is the standard cooling method, achieved by tuning the laser frequency to the red sideband of the bare atomic frequency, which leads to a resonant cooling contribution whereas heating is only induced by off-resonant carrier and blue-sideband excitations. To obtain well-resolved sidebands requires the linewidth of the atomic transition to be much smaller than the trapping frequency (strong confinement). Furthermore, taking into account that off-resonant heating excitations increase with increasing laser intensity the coupling Rabi-frequency has to be smaller than the trapping frequency (weak coupling) Sasura and Buzek (2002). Cooling to lower temperatures was reached using the concept of Raman sideband cooling Lindberg and Javanainen (1986), applicable for atoms Kerman et al. (2000) and ions Lee et al. (1996); Monroe et al. (1995), which allows the realisation of very small cooling linewidths by Raman two-photon transitions. Further improvements on the way to reach the ground state are achieved with dark state concepts based on eliminating off-resonant excitations via quantum interference. In EIT cooling Morigi et al. (2000); Roos et al. (2000) the concept of electromagnetically induced transparency Fleischhauer et al. (2005) is used to eliminate the carrier transition, which gives the main heating contribution in the sideband schemes. It uses a three level system coupled by two lasers, wherein one of the lasers (the stronger one) creates the absorption spectrum, a Fano like profile with a zero at the carrier transition and a peak at the red sideband transition. Therefore stronger couplings can be realised, no strong confinement is required due to the three level nature of the system and a large cooling bandwidth (i.e. frequency range for cooling) allows simultaneous cooling of several modes. An alternative approach, based on the Stark shift gate Jonathan et al. (2000), is the Stark shift cooling Retzker and Plenio (2007), cancelling the carrier transition for a correctly tuned Rabi frequency, which is of the order of the trapping frequency. This is achieved in the dressed state picture, i.e. taking the Stark shift of the laser interaction (originating from the zeroth order in the Lamb-Dicke parameter contributions of the resonant coupling laser) into account. The robust laser cooling scheme Cerrillo et al. (2010) combines the methods of EIT and Stark shift cooling resulting in a dark state for the combination of EIT and Stark shift Hamiltonian contribution and leading to an elimination of both the carrier transition and the blue sideband heating transition by interference.

The scheme presented here also eliminates the carrier and blue sideband transitions in lowest order in the Lamb-Dicke parameter, using an atomic three level system, a combinations of four laser beams and applying a magnetic field gradient. The application of an inhomogeneous magnetic field to trapped ions resulting in a position dependent resonance shift and mediating the coupling between internal and external degrees of freedom (e.g. allowing the coupling by microwaves with a non-vanishing effective Lamb-Dicke parameter) was discussed in Mintert and Wunderlich (2001); Wunderlich and Balzer (2003). For our method only the proper choice of the gradient strength and the phase relations of the laser beams form the conditions for a pure red sideband interaction up to first order in the Lamb-Dicke parameter.

The paper is organized as follows: In section II the cooling scheme is presented and the influence of the magnetic field gradient is analysed. It turns out, that in a Schrieffer-Wolff transformed picture, both the carrier as well as the blue sideband transition are eliminated by interference and thus cooling to the vibrational ground state in zeroth order is achieved. In section III the resonance condition for the cooling process is derived. Moreover it is shown that the maximal cooling amplitude is limited by the trapping frequency, a fact that arises from the amplitudes of the dressed states, which are tuned to resonance. Section IV presents an analytic formula for the cooling rate obtained by the method of adiabatic elimination and valid in the weak coupling case. Indeed a comparison to numerical results shows a good agreement in the range of small Rabi frequencies. In section V the dependence of the cooling rate on the spontaneous decay rate as well as on the magnitude of the Lamb-Dicke parameter is analysed. Furthermore the dependence of the cooling rate and final population on phase fluctuations of the coupling lasers is presented resulting in the conclusion that the scheme, albeit becoming slower, still provides fast cooling and low final populations even for relatively large deviations from the optimal phase relation. Section VI deals with the cooling of multiple vibrational modes, i.e. with multiple trapped ions. Only for the center of mass mode a complete elimination of heating terms (carrier and blue sideband) is possible; however the strength of the heating terms for other modes is smaller than in standard cooling schemes.

## Ii Cooling scheme and Hamiltonian

The scheme presented here is based on a complete elimination of the blue sideband heating terms in a Schrieffer-Wolff transformed picture (see eq. (6) and (7) below) in lowest order in the Lamb-Dicke parameter, which is achieved by a magnetic gradient field together with a four laser coupling scheme. It requires an atomic three level system consisting of a (fast decaying) excited state and a Zeeman split ground state and induced by the external position dependent magnetic field. Each of the Zeeman split levels is coupled by two lasers with Rabi frequency and opposite Lamb-Dicke parameters (, with the angle of the laser beam to the trap axis and the absolute value of the wavevector), obtained by irradiating the beams from opposite directions, to the excited level having a phase difference (see fig.1).

The ion is assumed to be trapped in an harmonic trapping potential with frequency and for simplicity only the motion along the trap axis is taken into account. Moreover, the expansion up to first order in terms of the Lamb-Dicke parameter requires the atom/ion to be localized on the scale provided by the laser wavelength (Lamb-Dicke limit, ). The position dependent magnetic field has the form

(1) |

and higher terms are not present or can be neglected in the following expansion in terms of the effective Lamb-Dicke parameter respectively. The corresponding energy shift due to the magnetic Zeeman splitting is (with the position operator , wherein and are the phonon creator and annihilator operators respectively) given by

(2) |

with the magnetic moment, the magnetic quantum number and , providing a coupling between the spatial position and the internal energy, i.e. phonons and internal atomic states. The effective Lamb Dicke parameter has a physical meaning as the energy difference between the left and right bound of the ground state wavefunction extent in units of . The energy shift of the homogeneous part is given by , wherein denotes the g-factor and the Bohr magneton. Thus the Hamiltonian of the system (cp. to fig. 1) is given by (setting )

(3) |

with being the energy of the upper level . In order to end up with a time independent Hamiltonian we assume for the detuning: , wherein is the frequency of the first laser pair (coupling to ) and the frequency of the second pair respectively (coupling to ). In an interaction picture, expanding the Hamiltonian up to the first order in the Lamb-Dicke parameter and applying the rotating wave approximation, the Hamiltonian takes the form

(4) |

In order to cancel the blue sideband heating contributions completely, it is important to adjust the magnetic gradient field in a way so that

(5) |

as will be explained later. Values for the necessary magnetic gradient to fulfill condition (5) for Ytterbium-172 and Calcium-40 ions are given in table 1.

Trapping frequency | ||||

Magnetic gradient /(T/m) | 607.6 | 1215 | 564.7 | 1130 |

Lamb-Dicke parameter | 0.13 | 0.09 | 0.25 | 0.18 |

The basic idea of the cancellation mechanism is that generates the red sideband and cancels the blue one. Since the phonon coupling introduces a term that is proportional to and the magnetic gradient introduces a term that is proportional to we can expect that if the phases are adjusted appropriately we can get a complete interference of the blue sideband.

The effect of the magnetic field is best seen by applying a unitary Schrieffer-Wolff transformation Schrieffer and Wolff (1966) to the Hamiltonian (4)

(6) |

resulting in the following transformations

(7) |

Application of the Schrieffer-Wolff transformation to the Hamiltonian (4) and keeping only terms up to first order in we get (fig. 1b, note that is equal to ):

(8) |

Herein we have defined:

(9) |

Note that the first order expansion terms of the Schrieffer-Wolff transformation () of the zeroth order term in (Lamb-Dicke-Parameter) Hamiltonian contributions (see 4) lead to a cancellation of the first order in blue sideband laser interaction contributions under the assumption (5).

The complete time evolution in the Schrieffer-Wolff transformed picture is given by

(10) |

The Liouvillian Superoperator describes the spontaneous decay from the excited level to the levels and . It exhibits the following form:

(11) |

describes the general spontaneous decay terms in zeroth order in the Lamb-Dicke Parameter:

(12) |

is the Schrieffer-Wolff transformation contribution of the spontaneous decay term (in first order in ):

(13) |

Note that recoil only contributes in order , leading to a very small correction to the cooling rate (decreasing it) due to the dark state nature of the cooling scheme and therefore can be neglected for small (which was tested by numerical simulations). This is due to the fact, that the probability for being in the excited state is given by and the spontaneous decay including recoil from this state is a process in second order in and therefore recoil only contributes with a strength .

If the final phonon population is reached, the system will be decoupled from the laser interaction and will end up in the state (“ dark state”). Note that due to the presence of only red sideband terms the phonon state will always be reached (unless ). In the original (i.e. non Schrieffer-Wolff transformed) picture, this final state has the form:

(14) |

Therefore the presented scheme cools to the vibrational ground state (n=0) in zeroth order in . Furthermore the Schrieffer-Wolff transformation preserves pure states and thus does not increase the entropy (which remains zero). A simple (projective) measurement of the internal state as well as a unitary rotation could further improve the achieved final population, that is they would result in a state with up to second order in .

## Iii Conditions for optimal cooling rates and limitations

At first sight it might seem possible to increase the cooling rate by an arbitrary amount just by increasing the laser Rabi frequency . However a detailed analysis of the -coupling terms (cp. Hamiltonian (8)) shows that this is not the case. These terms are changing the -state (shifting it and mixing it with the -state) and therefore crucially influence the cooling term () . Optimal rates are achieved for tuning the cooling interaction in resonance with one of the dressed states of the interaction. The Hamiltonian (8) in the dressed state basis takes the form

(15) |

wherein and is the dressed state basis (each having contributions of and states), and are the energies of those states (i.e. eigenvalues of the diagonalization procedure)

(16) |

and and the amplitudes of and in respectively (see also eq. (17)), i.e. the excited level expressed in the dressed state basis is given by

(17) |

In the interaction picture with respect to the first three Hamiltonian contributions, the resonance condition (for tuning the state to resonance) is easily seen to be fulfilled by setting . This procedure leads to the following resonance condition for optimal cooling rates:

(18) |

Consequently the detuning is - for a given trapping and laser Rabi frequency - uniquely determined in the optimal cooling case.

To see how the magnitude of the cooling term increases with increasing it is useful to define an “effective Rabi frequency“ .
If corresponds to the resonant state (the one that has been tuned to resonance), then is defined as

(19) |

and is a direct measure of the strength of the resonant (and therefore relevant) cooling contribution. In case of fulfilling condition (18) it follows that

(20) |

(see figure 2a).

From equation (20) it is clearly seen, that it is impossible to increase the cooling term contribution above a limiting upper value, or more precisely, is limited to a maximal value, which corresponds to the trapping frequency (). Therefore the trapping frequency is responsible for an upper limit to the maximal achievable cooling rate.

A second limitation arises from the decay rate and the contributions of the excited state in the new dressed state basis. Note that spontaneous decay to both dressed states and the state takes place. The cooling process relies on the decay to the state, and thus - to keep the process going - it is necessary that the other states decay fast enough, i.e. have a sufficient large contribution. Otherwise the system will be trapped (at least for a certain amount of time) in one of the dressed states and thus the cooling rate is reduces significantly. For equal decay rates in (12),(13), equal contributions in the dressed states would lead to the optimal case when taking only the spontaneous decay influence into account. The -amplitudes in the dressed states for the resonance case (18) are plotted in figure 2b. Note that the minimal amplitude determines the dynamics. After a point of equal -amplitudes in the dressed states, the minimal -amplitude decreases again for further increasing , which sets the optimal rate to a value , depending on the relation between and .

## Iv Rate equation formula and comparison to the numerical results

The procedure described in Cirac et al. (1992); Lindberg and Javanainen (1986) allows to calculate an analytic formula for the cooling rate by adiabatically eliminating the internal degrees of freedom under the assumption that the cooling dynamics is much slower than the internal dynamics. Moreover it is assumed that the trapping frequency is large compared to the cooling dynamics (i.e. ), which allows a projection on the phonon subspace; i.e., a derivation of a rate equation. This second condition is only fulfilled for small values of in our case (otherwise the cooling term is of the order of the trapping frequency or for further increasing even bigger).
The rate can be written in the form and the final steady state population is given by , wherein /- are the transition rates for increasing or lowering the phonon number respectively. Those quantities consist of two time correlation functions of internal operators, which can be calculated using the quantum regression theorem and the optical Bloch equations. It follows that, due to the fact that there exist no heating terms in first order in the Lamb-Dicke Parameter, is equal to zero and thus and the final phonon population is given by .

Following the procedure given in Cirac et al. (1992), we get for the rate (see appendix for the calculation)

(21) |

where the rate is defined via:

(22) |

This rate (21) gets optimal if the term in brackets is equal to zero and the resulting condition is equivalent to the resonance condition given in eq. (18). In that case, the rate takes the following form:

(23) |

The fact that the rate is proportional to is a direct consequence of the adiabatic elimination in high . More precisely, the rate is proportional to the spontaneous decay rate of the excited level times the probability for achieving this level. The latter quantity is proportional to , i.e. the probability for reaching the excited dressed state (here: ) decreases with increasing dissipation. Therefore the rate is given by: . However it does not mean that the optimal achievable rate is getting slower for increasing , because it is possible to reach higher final values for due to the fact, that the (internal) dynamics is less limited by the decreasing -amplitude of the resonant state, which occurs for increasing beyond a point of equal dressed state amplitudes. Figure 3 shows the rate obtained by numerical simulations and that one calculated using eq. (23) for different Rabi frequencies. As expected, the rate formula describes the rate quite well in the limit of small values whereas it fails completely in the opposite limit (when the cooling contribution is getting stronger and the conditions used for the formula derivation are no longer fulfilled).

## V Simulation and results

Figure 4 shows optimal cooling cycles for different spontaneous decay rates ( and ). In the first case (, the rate is of the magnitude of ) the optimal cooling rate is achieved for (cp. to eq. (18)) , i.e. the case of equal -amplitudes in the dressed states. That is what one would expect, because with further increasing to its maximal value there is one dressed state with a small amplitude and therefore the slower spontaneous decay from this state leads to a decreasing rate despite the increasing strength of the cooling term contribution. In the , (cooling rate) case, the spontaneous decay rate is much faster than the cooling dynamics. Therefore, as expected, it is possible to increase beyond the ”equal amplitude“ situation and thus a slightly higher effective Rabi frequency can be realized. The spontaneous decay time of the dressed states is much less limiting in this case. Note that this higher is needed to compensate the influence of the quantum Zeno effect, which reduces the rate compared to the case for a given . The optimal achieved rates (see also eq. (22)) are () for the case and () for the case respectively. Rates for different values of are presented in figure 3.

The dependence of the cooling rate on the spontaneous decay rate is plotted in fig. 5. As discussed before, for higher spontaneous decay rates higher values of are possible because of less restrictions due to the spontaneous decay and therefore the rate increases with increasing . For small values of , when the cooling dynamics is limited by the spontaneous decay rate and not by the coupling, the rate increases significantly with increasing . For high values of the rate converges to a limiting value for the rate. In this case the cooling dynamics is clearly limited by the fact, that reaches its maximal value (the trapping frequency (20), see right inset of fig. 5 for the corresponding values). Thus further increasing in the region for increasing (which is indeed necessary to compensate for the higher quantum Zeno effect) leads to a very small increasing of and therefore leads to a saturation of the rate with converging to its maximal value.

The dependence on (for otherwise fixed parameters) is shown in fig. 6. For , i.e. in the regime where heating terms of higher orders in can be neglected, the rate scales approximately linearly with increasing Lamb-Dicke parameter . Note that the rate equation (23) includes a quadratic contribution of , but also the optimal Rabi frequency depends on , scaling approximately and thus leading to an approximately linear scaling^{1}^{1}1More exactly, numerical simulations show that the approximate scaling is given by and . of the rate dependence on . Thus decreasing leads to increasing for the optimal condition, which can be explained as compensating for the quantum Zeno effect although it deteriorates the spontaneous decay amplitude behaviour. By increasing to higher values it reaches some optimal point and decreases again afterwards due to the fact that the heating terms in second order (and by further increasing also in higher orders) in become more and more important.

The cooling scheme concept is designed for the Lamb-Dicke regime (), because only in this regime higher order contributions can be neglected. Numerical simulations for higher Lamb-Dicke parameters show that the (minimal) final population always increases with increasing Lamb-Dicke parameter. While the final population for is given by corresponding to as expected according to (14), it is given by for , for and for (, ). Therefore only in the Lamb-Dicke regime the goal to achieve the vibrational ground state can be achieved. With respect to the new final population the cooling exhibits for higher Lamb-Dicke parameters in average again an exponential decay, leading e.g. for the above described case to a cooling rate of . Note however that the cooling rate and the final temperature are not optimized simultaneously for the same parameters, which especially for high Lamb-Dicke parameters leads to significant differences in the final populations and rates. The fact that both the final population and rate are worse for increasing Lamb-Dicke parameters underlines the restriction to the Lamb-Dicke regime.

Robustness to phase fluctuations — The presented scheme relies on destructive interference and thus may be fragile to dephasing, i.e., fluctuations in the phase difference between the laser pairs around the optimal value. The numerical analysis of the phase fluctuation is shown in fig. 7. It can be seen that even at fluctuations which are larger than the population is still less than and the rates decrease only by a factor of

The proposed cooling scheme is faster than previously proposed continuous cooling schemes. The fastest continuous scheme that was implemented is the EIT cooling scheme Roos et al. (2000) which cooled at a rate of . The robust cooling scheme Cerrillo et al. (2010) can in principle reach a rate of In contrary the proposed cooling scheme can reach rates which are one order of magnitude less than the trap frequency.

## Vi Multimode cooling

The presented cooling scheme was tested for the case of several trapped ions. Furthermore the cooling scheme, i.e. the laser coupling of internal and external degrees of freedom as presented in section II, is applied to one of the trapped ions (named ion 1 in the following). It is advantageous to express the Hamiltonian in terms of normal modes, wherein the creator and annihilator of the n-th normal mode is denoted by and respectively. In contrast, the ones of the i-th local mode are denoted by and . Expanding equation (4) to the multi particle case leads to the Hamiltonian

(24) |

Herein describes the frequencies of the corresponding normal modes, the summation over equals a summation over all normal modes and all internal states correspond to those of ion 1. With the transformation between the local coordinates and the normal coordinates the local operators in (24) can be expressed in terms of normal ones using the relation ()

(25) |

and thus the Hamiltonian (24) reads

(26) |

wherein .

Again, the effect of the magnetic gradient terms is best seen by moving to a Schrieffer-Wolff transformed picture in which the magnetic gradient displacements of the individual normal mode oscillators are removed. This transformation is given by

(27) |

with . Note that this transformation differs from the transformation used in sect. II what can easily be checked by inserting (25) into (6). Applying this transformation (cp. to eq. (7)) to the Hamiltonian (24) leads to

(28) |

Consequently, the amplitudes of the blue and red sideband terms for a specific mode depend on the relation of the trapping frequency to the normal mode frequency . A complete cancellation of blue sideband heating terms is only possible for the center of mass mode (COM), because only in this case . For higher modes heating terms are always present, although they are off-resonant by tuning the laser frequency to the red sideband. In general, for increasing the heating amplitude increases and the cooling amplitude decreases, which leads to equal amplitudes for . The latter case corresponds to the situation in cooling schemes that do not cancel blue sideband terms, e.g. EIT cooling and sideband cooling. Therefore the presented cooling scheme results in an improved cooling behaviour with respect to those schemes also in the multiparticle case albeit not being able to cool other modes than the COM-mode without heating terms.

The cooling of the individual modes (COM mode (blue), first excited mode (breathing mode, red) and second excited mode (egyptian mode, green)) for the case of three ions is presented in fig. 8. For three ions Sasura and Buzek (2002) , , , . The left plot shows the cooling starting with an initial population for each mode. Due to the fact that in the higher excited modes heating terms are always present while in the COM mode those contributions are eliminated completely in first order in , the COM mode is always cooled with a faster rate also in the case when the cooling terms of any of the excited modes is set to resonance, i.e. the COM cooling transitions are off-resonant. The right plot shows the heating of the non-resonant excited modes originally in during the resonant cooling process of the center of mass mode starting from . The modes are cooled/heated to final populations of , and . That is the heating turns out to be very small and increases with as expected from eq. (28).

## Vii Conclusion and comments

We presented a cooling scheme in the Lamb-Dicke regime, which completely eliminates the influence of the blue sideband heating transitions as well as the carrier transitions in lowest order in the Lamb-Dicke parameter, resulting in a good cooling behaviour and providing cooling rates of the order of one magnitude below the trapping frequency. Albeit the Lamb-Dicke regime required for the Hamiltonian expansion and the resonance condition limit the cooling rate well below the trapping frequency, the achieved rates are faster than in comparable existing schemes working in the same regime (e.g. EIT cooling Morigi et al. (2000), Robust cooling scheme Cerrillo et al. (2010)). This is achieved by a proper choice of the magnetic gradient strength whose influence was analysed in a Schrieffer-Wolff transformed picture. Optimal cooling conditions were derived together with a formula valid in the weak laser coupling case. The dependence on the Lamb-Dicke parameter as well as on the spontaneous decay rate was presented. Moreover the cooling behaviour under deviations around the optimal phase differences between the coupling laser pairs was analysed showing that deviations of about 20 percent only reduce the rates by a factor of two and the final population is still of the order of , both very satisfying values compared to other cooling schemes working in the Lamb-Dicke regime. Furthermore an analysis of the multiparticle case revealed good cooling behaviour also for several modes albeit not being able to eliminate the carrier and blue sideband terms completely up to first order in in other modes than the center of mass mode.

###### Acknowledgements.

Support by the European Commission (STREPs HIP and PICC), the Alexander von Humboldt foundation, the Deutsche Forschungsgemeinschaft, and secunet AG is acknowledged.## Appendix A Rate equation derivation

In this appendix, the calculation leading to the rate equation formula (21) is presented. The master equation (10) up to first order in can be written in the following form

(29) |

wherein the lower index indicates the order in the Lamb-Dicke parameter. The Liouvillians are given by

(30) |

(31) |

with the definitions of the spontaneous decay contributions and given in eq. (12) and eq. (13) respectively. The complete first order Liouvillian is defined as the sum of the two first order contributions, i.e. . The zero and first order Hamiltonian contributions (cp. to eq. (8)) are

(32) |

(33) |

where are pure external operators and pure internal ones: , , and .

The projector on the relevant subspace is given by

(34) |

with the external/vibrational density matrix contribution , the internal steady state and and the projectors on the internal steady state (assuming that the internal dynamics is much faster than the cooling dynamics) and the vibrational subspace respectively. denotes the trace over the internal part. The master equation for the relevant external part up to first order in the Lamb-Dicke parameter is given by (see e.g. Cirac and Zoller (1995))

(35) |

Tracing equation (35) over the internal part, the master equation for the external part is obtained

(36) |

wherein it was used that , which follows directly from eq. (34). The last step takes into account that the terms (originating from the Schrieffer-Wolff transformation of the spontaneous decay terms) do not contribute in the master equation, because

(37) |

with and furthermore assuming equal -rates in .

Inserting (33) into eq. (31) the first order Liouvillian can be written as

(38) |

using the definitions and . Thus, inserting (38) into (36) and noting that

(39) |

one obtains

(40) |

with the internal operator two time-correlation functions and .

Furthermore, using that and defining we end up with

(41) |

Inserting the concrete expressions for the external and internal operators (see definitions below eq. (33)) and applying the projection (eq. (34)) the master equation for the external part turns out to be

(42) |

wherein

(43) |

The time evolution of the vibrational state population () following from eq. (42) takes the form

(44) |

with the phonon number increasing / decreasing rates

(45) |

Therefore the equation of motion for the phonon population expectation value is found to be

(46) |

with the steady state solution

(47) |

Thus it remains calculating the rates and , i.e. calculating the two-time correlation functions (43). Using the Bloch equations (the internal equation of motion in the zeroth order in the Lamb-Dicke parameter) and the quantum regression theorem, one obtains the results in section IV.

## References

- Wineland et al. (2005) D. J. Wineland, D. Leibfried, M. D. Barrett, A. Ben-Kish, J. C. Bergquist, R. B. Blakestad, J. J. Bollinger, J. Britton, J. Chiaverini, B. Demarco, D. Hume, W. M. Itano, M. Jensen, J. D. Jost, E. Knill, J. Koelemeij, C. Langer, W. Oskay, R. Ozeri, R. Reichle, T. Rosenband, T. Schaetz, P. O. Schmidt, and S. Seidelin, Proceedings of the International Conference on Laser Spectroscopy (ICOLS) (2005).
- Sasura and Buzek (2002) M. Sasura and V. Buzek, J. Mod. Opt., 49, 1593 (2002).
- Hänsch and Schawlow (1975) T. W. Hänsch and A. L. Schawlow, Opt. Commun., 13, 68 (1975).
- Wineland and Dehmelt (1975) D. J. Wineland and H. Dehmelt, Bull. Am. Phys. Soc., 20, 637 (1975).
- Stenholm (1986) S. Stenholm, Rev. Mod. Phys., 58, 699 (1986).
- Wineland et al. (1978) D. J. Wineland, R. E. Drullinger, and F. L. Walls, Phys. Rev. Lett., 40, 1639 (1978).
- Lindberg and Javanainen (1986) M. Lindberg and J. Javanainen, J.Opt.Soc.Am.B, 3, 1008 (1986).
- Kerman et al. (2000) A. J. Kerman, V. Vuletic, C. Chin, and S. Chu, Phys. Rev. Lett., 84, 439 (2000).
- Lee et al. (1996) H. J. Lee, C. S. Adams, M. Kasevich, and S. Chu, Phys. Rev. Lett., 76, 2658 (1996).
- Monroe et al. (1995) C. Monroe, D. M. Meekhof, B. E. King, S. R. Jefferts, W. M. Itano, and D. J. Wineland, Phys. Rev. Lett., 75, 4011 (1995).
- Morigi et al. (2000) G. Morigi, J. Eschner, and C. H. Keitel, Phys. Rev. Lett., 85, 4458 (2000).
- Roos et al. (2000) C. F. Roos, D. Leibfried, A. Mundt, F. Schmidt-Kaler, J. Eschner, and R. Blatt, Phys. Rev. Lett., 85, 5547 (2000).
- Fleischhauer et al. (2005) M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. Mod. Phys., 77, 633 (2005).
- Jonathan et al. (2000) D. Jonathan, M. B. Plenio, and P. L. Knight, Phys. Rev. A., 62, 042307 (2000).
- Retzker and Plenio (2007) A. Retzker and M. B. Plenio, New J. Phys., 9, 279 (2007).
- Cerrillo et al. (2010) J. Cerrillo, A. Retzker, and M. B. Plenio, Phys. Rev. Lett., 104, 043003 (2010).
- Mintert and Wunderlich (2001) F. Mintert and C. Wunderlich, Phys. Rev. Lett., 87, 257904 (2001).
- Wunderlich and Balzer (2003) C. Wunderlich and C. Balzer, Adv. Atom. Mol. Opt. Phys., 49, 293 (2003).
- Schrieffer and Wolff (1966) J. R. Schrieffer and P. A. Wolff, Phys. Rev., 149, 491 (1966).
- Cirac et al. (1992) J. I. Cirac, R. Blatt, P. Zoller, and W. D. Phillips, Phys. Rev. A, 46, 2668 (1992).
- Cirac and Zoller (1995) J. I. Cirac and P. Zoller, Phys. Rev. Lett., 74, 4091 (1995).