Next Article in Journal
A New Era for the Early and Non-Invasive Diagnosis of Giant Cell Arteritis: The Use of Fast-Track Ultrasound in Clinical Practice
Next Article in Special Issue
Accurate Characterization of the Properties of the Rare-Earth-Doped Crystal for Laser Cooling
Previous Article in Journal
Detailed and Average Models of a Grid-Connected MMC-Controlled Permanent Magnet Wind Turbine Generator
Previous Article in Special Issue
Saturation, Allowed Transitions and Quantum Interference in Laser Cooling of Solids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Laser Cooling beyond Rate Equations: Approaches from Quantum Thermodynamics

School of Physics, Trinity College Dublin, DN02 PN40 Dublin, Ireland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(3), 1620; https://doi.org/10.3390/app12031620
Submission received: 19 November 2021 / Revised: 27 January 2022 / Accepted: 1 February 2022 / Published: 3 February 2022
(This article belongs to the Special Issue Laser Cooling of Solids: Novel Advances and Applications)

Abstract

:
Solids can be cooled by driving impurity ions with lasers, allowing them to transfer heat from the lattice phonons to the electromagnetic surroundings. This exemplifies a quantum thermal machine, which uses a quantum system as a working medium to transfer heat between reservoirs. We review the derivation of the Bloch-Redfield equation for a quantum system coupled to a reservoir, and its extension, using counting fields, to calculate heat currents. We use the full form of this equation, which makes only the weak-coupling and Markovian approximations, to calculate the cooling power for a simple model of laser cooling. We compare its predictions with two other time-local master equations: the secular approximation to the full Bloch-Redfield equation, and the Lindblad form expected for phonon transitions in the absence of driving. We conclude that the full Bloch-Redfield equation provides accurate results for the heat current in both the weak- and strong- driving regimes, whereas the other forms have more limited applicability. Our results support the use of Bloch-Redfield equations in quantum thermal machines, despite their potential to give unphysical results.

1. Introduction

Laser cooling [1,2,3], in both atomic and solid-state systems, is now a well established technique. In solids, particularly rare-earth-doped glasses, cooling can be achieved by using anti-Stokes fluorescence of the dopants. It provides an example of a quantum thermal machine [4,5,6], in which is a discrete quantum system—in this case, the energy levels of rare earth ion—is the working medium. This working medium couples to two heat baths and a source of work, namely the phonon and photon reservoirs and the driving laser, allowing it to operate as a refrigerator.
Laser cooling is generally modeled using rate equations for the populations of the levels. This approach can also be used for semiconductors, where the rate equations refer to the populations of the electron and hole bands. However, such approaches cannot capture certain effects which, while not expected to be relevant in systems such as rare earths, are increasingly important in quantum thermodynamics more generally. These include the role of coherences in determining heat flows, which have been argued to offer enhanced performances in various quantum thermal machines [7,8,9,10,11]; the effects of strong driving, which can modify the energy levels through the a.c. Stark effect [12,13,14], and so impact on the heat flows [11,15]; and the effects of spectral structure in the heat baths. This last can be considered in two regimes: for strongly structured baths one can expect non-Markovian behavior [16,17,18], whose impact on thermodynamics remains a challenging open topic. However, the spectral structure can be important even where a Markovian description remains appropriate [19]. An important practical target for thermodynamic machines is to maximize their power, and the heat flows to a bath are determined by its spectral density. Thus to achieve maximum power one must consider the spectral structure of baths, if there is any on the energy scales of the working medium. Examples of systems where this occurs include quantum-dot excitons coupled to acoustic phonons [20], colour centres in diamond [21,22], and superconducting circuits [23].
These issues can be treated theoretically by studying models of an open quantum system in which the working medium interacts with its surrounding heat baths. Such models are tractable in the weak-coupling, Markovian regime, where they lead to time-local equations of motion such as the Bloch-Redfield equation [24]. Those approaches can be extended to allow calculations of heat and work in the quantum regime [25]. However, there are several time-local equations which can be obtained, using reasonable approximations, from a given model, and these can make differing predictions for the dynamics [26,27]. This problem has been addressed by several groups, who argue that the Bloch-Redfield equation [10,26,28,29,30,31,32,33] is useful and indeed accurate, despite its potential pathologies [34]. In this paper, we extend such studies to explore the heat flows in a simple laser cooling process, with the aim of identifying an approximate time-local equation that can accurately model them.
In the following, we first review the derivation of the Bloch-Redfield equation for an open quantum system, and outline its extension to calculate heat flows. We also discuss two other time-local equations which can be obtained by making further approximations: a Lindblad form in the energy eigenbasis, obtained by making the secular approximation, and a Lindblad form in the eigenbasis of the undriven system. We use these forms to calculate the cooling spectrum, i.e., the cooling power as a function of driving frequency, in a model of laser cooling. The model allows for strong driving and includes a spectral structure for the environment. We find that a complete description of the cooling spectrum, which covers both the weak-driving and strong-driving regimes, can be achieved using the full Bloch-Redfield equation. We provide further support for the correctness of the Bloch-Redfield master equation—whose use has been controversial because it does not guarantee positivity [34], and can lead to behavior inconsistent with thermodynamic principles [35]—by comparing its predictions to those of an exact numerical method. Our conclusions support the use of Bloch-Redfield equations to model laser cooling and other thermodynamic processes [10,30,32,33,36,37].

2. Materials and Methods

2.1. Laser Cooling Model

We consider a simple laser-cooling scheme, depicted in Figure 1, involving an impurity with two states forming a ground-state manifold, and a single state in an excited-state manifold. A driving laser of frequency ω l excites the transition from the upper state of the ground-state manifold to an excited state | e and energy E 0 above the top of the ground-state manifold. We assume this state decays by radiative emission to the ground state. Crucially, the two states within the ground-state manifold, | g l , | g u , which are split by an energy E man , are coupled by the emission and absorption of lattice phonons.
We will analyze this problem within the standard framework of open quantum systems theory, separating it into parts corresponding to a ‘system’, one or more baths, and the interactions between the system and the baths. In our case, the baths are the phonons, and the continuum of photon modes that give the radiative decay. We describe those aspects of the model in Section 2.4 below. The system, meanwhile, comprises the states of the impurity and their interaction with the laser, which can be treated as a classical driving field [38]. The system Hamiltonian is then ( = 1 )
H S = E 0 | e e | E man | g l g l | + Ω cos ω l t | g u e | + | e g u | ,
where Ω is the Rabi frequency given by the product of the electric field amplitude of the driving laser and the dipole moment of the transition, d . E 0 . The time-dependence of the driving field can be removed, in the rotating wave approximation, by using a unitary transformation
U = exp ( i ω l t | e e | ) .
In this frame the driving field is time-independent, and the Hamiltonian for the system is
H S = δ Ω / 2 0 Ω / 2 0 0 0 0 E man .
Here δ = ω l E 0 denotes the driving laser frequency relative to the transition. The diagonal terms in Equation (2) are the energies of the electronic states, in the rotating frame. The off-diagonal terms are the coupling between those states produced by the electric-dipole interaction with the driving field [38]. The coupling between the impurity states and the lattice phonons is considered in the following, where we discuss master equations for open quantum systems and specify the system-bath interaction Hamiltonian.

2.2. Master Equations for Open Quantum Systems

Figure 1 depicts an open quantum system: one which interacts, explicitly or implicitly, with a wider environment. These interactions lead to an exchange of energy between system and environment, and dephasing and decoherence effects. Here, we have an environment comprising the phonons in the host crystal of the impurity and the photons associated with the radiative decay of the upper level.
The dynamics of an open quantum system can, in certain circumstances, be described by a time-local master equation for its reduced density matrix [24]. Such equations can be obtained from microscopic models which consider the environment explicitly in making the weak-coupling and Markovian approximations. They are also often postulated phenomenologically, based on the observation that the most general equation of motion is one of Lindblad form. However, there are several different forms of equations that can result from a microscopic model, depending on the details of the approximations made. The predictions of these forms can, furthermore, differ from those based on phenomenological Lindblad forms.
These issues have been discussed in previous works [26,28] which suggest that the full Bloch-Redfield equation—obtained by using the weak-coupling and Markovian approximations, but without making the secular approximation—gives a good description of the dynamics. This is in spite of the fact that the Bloch-Redfield equation does not guarantee that the eigenvalues of the reduced density matrix remain positive [24,34]. For a system where there are no degeneracies, or near-degeneracies, that issue can be cured by secularization [34,39], which corresponds to eliminating oscillating terms in the dissipator that average to zero over time. This leads to a Lindblad form [40,41] with positive rates. It is, however, a priori invalid for the laser-cooling protocol considered here, where weak driving near resonance means we have Ω 0 and δ 0 , so that two of the eigenstates of Equation (2), | g u and | e , are almost degenerate in the frame where the laser field is time independent.
A fairly generic form for the Hamiltonian of an open quantum system is
H = H S + H B + H S B
H S B = k g k O ( b k + b k ) .
Here H S is the Hamiltonian for the system, H B for its environment, or bath, and H S B is the system-bath coupling. We consider the common situation in which the bath comprises a set of harmonic oscillators [24], which we index using a quantity or quantities labeled r. Note that r denotes the full set of quantum numbers required to label the modes. The oscillators have frequencies ω r , and ladder operators b r and b r . The displacement of the rth bath mode is coupled to the system operator O, with coupling strength g r . The dissipative effects of the bath depend on its spectral density, J ( ω ) = r g r 2 δ ( ω ω r ) .
To fix notation we recall the standard procedure for deriving a Bloch-Redfield master equation [24,42,43]. We work in the interaction picture with respect to H S + H B , so that O ( t ) = e i H S t O e i H S t . Note that where necessary we will distinguish operators in the interaction and Schrödinger pictures as, for example, O ( t ) and O. From the von Neumann equation we obtain the form
d ρ ( t ) d t = i [ H S B ( t ) , ρ ( 0 ) ] 0 t d t [ H S B ( t ) , [ H S B ( t ) , ρ ( t ) ] ]
where ρ ( t ) is the full density operator of the system and environment. For weak coupling to a bath one can replace ρ ( t ) ρ S ( t ) ρ B ( t ) on the right-hand side, where ρ S is the reduced density matrix of the system, and ρ B that of the bath. Since the bath is macroscopic it can be assumed to be unperturbed by the system, and ρ B taken to be a thermal state at inverse temperature β . For a Markovian system one may, furthermore, approximate ρ S ( t ) ρ S ( t ) . We can write the coupling operator in the eigenbasis of H S as
O ( t ) = i j e i ( E i E j ) t i | O | j | i j | i j O ^ i j ( t ) .
Taking the trace of Equation (5) over the environment’s degrees-of-freedom we find
d ρ S ( t ) d t = i j { A i j [ O ^ j i ( t ) ρ S ( t ) O ( t ) + O ( t ) ρ S ( t ) O ^ i j ( t ) ρ S ( t ) O ^ i j ( t ) O ( t ) O ( t ) O ^ j i ( t ) ρ S ( t ) ] i B i j [ O ^ j i ( t ) ρ S ( t ) O ( t ) O ( t ) ρ S ( t ) O ^ i j ( t ) + ρ S ( t ) O ^ i j ( t ) O ( t ) O ( t ) O ^ j i ( t ) ρ S ( t ) ] } .
The quantities A i j and B i j are related to the the real-time Green’s functions of the environment at the transition frequency ν i j = E i E j connecting levels i and j. The quantities A i j are associated with dissipation, and are
A i j = π { [ n ( ν i j ) + 1 ] J ( ν i j ) + n ( ν j i ) J ( ν j i ) } .
Here n ( ν > 0 ) = 1 / ( exp ( β ν ) 1 ) is the Bose function describing the bath occupation, and J ( ν ) = 0 for ν < 0 . The first term in A i j corresponds to the creation of a bath quantum as the system transitions from a state i to j with E i E j > 0 , whereas the second corresponds to the absorption of a bath quantum in the opposite case, E i E j < 0 . The quantities B i j are associated with energy shifts, and are given by the principal value integral
B i j = P J ( ω ) ω + ( 2 n ( ω ) + 1 ) ( E i E j ) ω 2 ( E i E j ) 2 d ω .
Equation (7) can be used directly, but is often further approximated, leading to other forms of equation-of-motion for an open quantum system. One very common approximation is to drop the principal value terms proportional to B i j . Another common approximation is to secularize the equation-of-motion. This is done by decomposing the remaining coupling operators, O ( t ) , into the energy eigenbasis: O ( t ) = k l O ^ k l ( t ) . Every term in Equation (7) then involves a product of operators corresponding to two transitions, one involving the pair of levels i and j, and one involving the pair k, l. If the levels are non-degenerate these products of operators are, in general, time-dependent in the interaction picture, and average to zero. The exception is where a transition in one direction is paired with the same transition in the opposite direction, so that the time-dependence cancels out. Retaining only those terms the dissipative part of Equation (7) becomes
d ρ S ( t ) d t = i j 2 A i j O ^ j i ( t ) ρ S ( t ) O ^ i j ( t ) 1 2 [ ρ S ( t ) , O ^ i j ( t ) O ^ j i ( t ) ] + ,
where [ A , B ] + = A B + B A is an anticommutator. This is of Lindblad form, and therefore guarantees the positivity of the density operator. It has a straightforward physical interpretation: the environment causes transitions from the system state i to the system state j at rate 2 A i j .

2.3. Heat Flows from Master Equations

The method of full counting statistics [25] allows one to extend the approaches above so as to compute the heat transferred to the bath. It has been used, often with the secular approximation [44], to obtain master equations and study heat statistics in various systems, including driven quantum-dot excitons [11,15], a driven two-level system [45], a steady-state (absorption) refrigerator [10,30,36], and a two-bath spin-boson model [32,33]. The absorption refrigerator and spin-boson model have been studied using the full Bloch-Redfield approach, without the secular approximation, which highlights the role of coherences [10,30,36]. Here we give an outline of the method and present a complete form for the full counting-field Bloch-Redfield equation, which we shall use to calculate laser cooling spectra.
The heat transferred to a bath is, by definition, the change in its energy between two times. Thus we consider a process involving projective measurements of the bath energy at two times. We take the initial time to be t i = 0 , and suppose that at this time the system and bath are in a product state, ρ S ( 0 ) ρ B . We can then consider the probability distribution of the heat, P ( Q , t ) , which is the probability that the energy measurements of the bath at times t i and t give results differing by Q. It is convenient also to introduce the characteristic function of the heat distribution, χ ( u , t ) = d Q P ( Q , t ) e i u Q . The variable u is known as the counting field. (This term should not be taken to imply that heat is necessarily a discrete, countable quantity. It arises from other uses of the method, such as calculations of the number of electrons transferred across a tunnel junction [25].).
One can evaluate χ ( u , t ) by introducing an annotated density operator, ρ u ( t ) , such that χ ( u , t ) = Tr ρ u ( t ) . ρ u ( t ) has a non-unitary time evolution given by
ρ u ( t ) = U u / 2 ρ u ( 0 ) U u / 2 ,
where U u is related to the normal time-evolution operator, U = e i H t , by
U u = e i u H B U e i u H B .
Note the similarity between these phase factors and the factor e i u Q in the definition of the characteristic function; it is these factors that incorporate the results of the measurements of the bath energy, H B , into ρ u ( t ) . At the initial time the annotated density matrix is given by ρ u ( 0 ) = ρ ( 0 ) .
For a general operator P we define the annotated version P u = e i u H B P e i u H B , which obeys the Heisenberg-like equation
i d P u d u = [ P u , H B ] .
For the lowering operator appearing in Equation (4) we have b u , k = e i ω k u b 0 , k . Thus, the time-evolution operators, U ± u / 2 , can be obtained from the standard form, e i H t , by replacing the coupling Hamiltonian, Equation (4), with H S B ± = g k O ( b k e i ω u / 2 + b k e ± i ω u / 2 ) .
A master equation for the reduced annotated density matrix, ρ u , S ( t ) can now be obtained, following the steps above. The essential difference is that the von Neumann equation for ρ ( t ) , in the interaction picture, must be replaced by
d ρ u ( t ) d t = i ( H S B + ρ u ( t ) ρ u ( t ) H S B ) .
The result is
d ρ u , S ( t ) d t = i j { A i j [ e i u ( E i E j ) ( O ^ j i ( t ) ρ u , S ( t ) O ( t ) + O ( t ) ρ u , S ( t ) O ^ i j ( t ) ) ρ u , S ( t ) O ^ i j ( t ) O ( t ) O ( t ) O ^ j i ( t ) ρ ( t ) ] i B i j [ e i u ( E i E j ) ( O ^ j i ( t ) ρ u , S ( t ) O ( t ) O ( t ) ρ u , S ( t ) O ^ i j ( t ) ) + ρ u , S ( t ) O ^ i j ( t ) O ( t ) O ( t ) O ^ j i ( t ) ρ ( t ) ] } .
This form differs from Equation (7) by the addition of phase factors in the four terms that cause transitions between the system eigenstates. It can be approximated as discussed above, by dropping the principal value terms, or by making the secular approximation.
The mean heat is
Q = Q P ( Q ) d Q = i d χ d u u = 0 = i Tr d ρ u , S ( t ) d u u = 0 .
From Equation (14) we find that the heat current is
d Q d t = i j { A i j [ ( E i E j ) Tr ( O ^ j i ( t ) ρ S ( t ) O ( t ) + O ( t ) ρ S ( t ) O ^ i j ( t ) ) ] i B i j [ ( E i E j ) Tr ( O ^ j i ( t ) ρ S ( t ) O ( t ) O ( t ) ρ S ( t ) O ^ i j ( t ) ) ] } .
This can be used to calculate the heat current from the density matrix, ρ u = 0 , S ( t ) = ρ S ( t ) , obtained by solving the standard Bloch-Redfield Equation (7).

2.4. Master Equations for Laser Cooling

We consider a model in which the system Hamiltonian is given by Equation (2). We suppose that there is a continuum of phonons responsible for transitions between the states of the ground-state manifold. This phonon bath will be described by Equations (3) and (4), with coupling operator O = | g l g u | + | g u g l | . For the spectral density of this bath, we take the super-Ohmic form with an exponential high-frequency cut-off, J ( ω ) = 2 α ( ω 3 / ω c 2 ) exp ( ω / ω c ) . We are not targetting a detailed model of a real system, and this form is chosen largely for illustrative purposes. It may, however, be noted that it corresponds to that for acoustic phonons coupling to localized impurities such as the silicon-vacancy center in diamond [22] or a quantum-dot exciton [20]. α is a dimensionless measure of the coupling strength, and ω c is a high-frequency cut-off. Such cut-offs arise from the size of the electronic states and correspond roughly to the phonon frequency at a wavelength given by that size.
We also consider, in the following, an alternative form of dissipator, of standard Lindblad form. For a transition caused by a jump operator A, with rate γ A , the standard Lindblad form is
d ρ S ( t ) d t = γ A L A ρ S ( t ) = γ A A ρ S ( t ) A 1 2 [ ρ S ( t ) , A A ] + .
Thus the natural phenomenological form, capturing the processes shown in Figure 1, is to combine two of these dissipative terms, one for phonon absorption, with rate γ + and jump operator σ + = | g u g l | , and one for phonon emission, with rate γ and jump operator σ = σ + = | g l g u | . Such a form corresponds to Equation (10) when the eigenstates of H S are simply | g l and | g u , which is resonable for weak driving. This comparison allows us to identify the appropriate rates, from Equation (8), as γ = 2 π ( n ( E man ) + 1 ) J ( E man ) and γ + = 2 π n ( E man ) J ( E man ) .
In addition to the phonon dissipation, our model involves the radiative decay of the excited state, | e to the ground state | g l . We model this as a Lindblad form with jump operator | g l e | , and rate γ .

2.5. Exact Methods

As well as results of master equations, we shall present, in the following, calculations of the heat flows obtained by numerically-exact simulations [46,47,48] of the model open quantum system described above. The technique, known as TEMPO, calculates the path integral for the evolution of an open quantum system, discretizing time into a series of steps [49]. It uses a matrix-product state representation to efficiently store the augmented density tensor, which allows it to consider large memory times for the bath [47,48]. Combining path-integral methods with the counting-field technique [32,46], allows calculations of the total heat transferred to the phonon bath up to a particular time. Details of the method, and the associated code, are given in Ref. [46]. We use it to calculate the heat currents to the phonons by taking the difference of the total heat transferred to the bath between two times, separated by a single timestep. In these calculations, the dynamics of the system and the effects of the phonon bath are treated exactly. We do not treat the radiative decay in this first-principles fashion, but rather include it using the same Lindblad form we use for the master equation approach. We believe this is appropriate, in as much as the bath associated with radiative decay has no spectral structure, in contrast with that associated with the phonons. The TEMPO approach has recently been extended to simulations with multiple baths [50], which would allow it to treat laser cooling with structured photon environments, e.g., in optical resonators.

3. Results

The parameters in our model are the energy splitting of the ground-state manifold, the detuning and Rabi frequency of the driving, the radiative decay rate, γ , the cut-off frequency, ω c , the dimensionless coupling, α , and the temperature T. We choose the energy and time units such that E man = 2 . For the remaining parameters we take γ = 0.5 , ω c = 1 , α = 0.01 , and T = 3 . These parameters are not intended to be realistic but are chosen to allow us to compute the exact solutions with a reasonable effort, and compare the results of the different master equations. In particular, we choose a large value for the radiative decay rate, γ , to increase the magnitude of the heat current. It may be noted that for these parameters the phonon absorption rate, γ + 0.14 , is comparable to, but smaller than, the radiative decay rate. This differs from the situation for conventional laser cooling, appropriate in systems such as rare-earth ions, where the phonon rates are much larger than those for radiative decay [1], and the electronic populations are very close to equilibrium. It implies that, in our case, the heat current will be limited by the driving strength (for weak driving) or the phonon rate (for strong driving), and not the radiative lifetime.
Figure 2 shows the calculating cooling power as a function of the detuning, δ , for four different strengths of the driving field. The different curves are computed using the full Bloch-Redfield Equation, (7), the phenomenological Lindblad form, Equation (17), and the secular Bloch-Redfield Equation, (10). Considering first weak driving, in Figure 2a, we see that the Bloch-Redfield and phenomenological theories agree well, and give a cooling profile that appears to be Lorentzian, as one would expect. While the secular Bloch-Redfield equation agrees away from the resonance, we see that it fails close to it, massively overestimating the cooling power. The secular approximation is, of course, not justified here, because there are near degeneracies in the Hamiltonian. Nonetheless, the level of disagreement seems surprising, given the agreement away from resonance.
In the converse, strong-driving region, Figure 2d, all three methods give similar results. However, there is a noticeable difference on the high-energy side of the transition, with the phenomenological theory giving, as before, a Lorentzian profile, while the other theories predict the heat current drops off more rapidly, and indeed switches direction, from cooling to heating, in the range of detunings shown.
Figure 3 shows the temperature dependence of the cooling power for weak resonant driving. As noted above, the secular approximation is inappropriate in this regime and massively overestimates the cooling power. The other theories agree closely and appear physically reasonable, predicting that the cooling power drops rapidly once the temperature is lowered below the splitting E man . This is the expected physical behavior for laser cooling in a discrete level structure, caused by the vanishing of the phonon occupation at temperatures much less than E man . The precise behavior at very low temperatures is not relevant to laser cooling since it is not expected to operate there. Nonetheless, it may be noted that, for these parameters, the Bloch-Redfield theory predicts a very small but negative cooling power, i.e., net heating, at very low temperatures ( T < 0.38 ). However, this is an approximate theory whose accuracy is not sufficient to discern the true behavior of the heat current in this regime. Indeed, we find that at these very low temperatures the theory does not predict a physical density matrix, giving one which has a negative, albeit very small, eigenvalue.
In Figure 4 we compare the cooling spectra predicted by the full Bloch-Redfield equation with those obtained from the exact numerical method [46]. The numerical method simulates the time-evolution of the open quantum system, using discrete timesteps. For these simulations, we have taken a timestep d t = 0.05 , and computed the heat current, at a time t = 30.0 , from the difference in the heat transfer at two times. The mean heat transfer is computed by evaluating the annotated reduced density matrix, ρ u , S ( t ) , and computing the finite difference approximation to the derivative in Equation (15) from the values of Tr ρ u , S ( t ) at u = 0.05 and u = 0.0 . The numerical accuracy and convergence of these simulations involve two further parameters: a maximum number of timesteps retained in the influence functional, K, and a cut-off parameter controlling the truncation of the singular-value decompositions. We take K = 100 , and use a cut-off of 10 7 .
We see from Figure 4 that the Bloch-Redfield equation is in excellent agreement with the numerical results. The non-Lorentzian behavior of the cooling profile on the high-energy side, predicted by the full Bloch-Redfield and secular equations, is present. There is a slight overestimate of the cooling power in the tails of the profiles, which we believe is because of the heat current has not yet reached its steady-state value at those small values of the cooling power.

4. Discussion

Figure 2 and Figure 4 suggest that the full Bloch-Redfield equation gives an accurate account of the cooling profile, in both the weak- and strong-driving cases. In the weak-driving case, it agrees with the phenomenological theory, which is well-justified for weak-driving, while in the strong-driving case it agrees with the secular theory, which is well-justified there. Furthermore, it agrees with exact numerical results, in the strong-driving regime where such simulations are possible. Thus the Bloch-Redfield equation allows a complete treatment of both regimes, using a single equation. This conclusion is similar to previous conclusions on the dynamics and thermodynamics of other open quantum systems, where the Bloch-Redfield equation has similarly been argued to provide the most accurate description [10,26,28,29,30]. This is in spite of the possibility that it produces unphysical density matrices with non-positive eigenvalues. That possibility does not occur in our results, except at very low temperatures where laser cooling would not, in any case, be expected to operate.
In previous works, it has been noted that the secular approximation does not allow for the presence of bath-induced or noise-induced coherence [51] in multilevel systems where near-degenerate levels have different couplings to the bath [26]. This phenomenon can play an important role for the heat currents, as has been pointed out previously for quantum absorption refrigerators [10,30]. Its significance in our case can be seen by comparing the secular result (where there is no bath-induced coherence) and the Bloch-Redfield result (where there is) in Figure 2. When δ = 0 , Ω = 0 our Hamiltonian has two degenerate eigenstates, but the form of those eigenstates depends on how the limit is taken: for δ = 0 , Ω 0 they are | ± = | g u ± | e , but for Ω = 0 , δ 0 they are | g u and | e . The naive form of secular approximation in Equation (10) produces a dissipator which populates the states | + and | , and destroys coherences between them. However, we observe that the phonon bath couples only to | g u | + + | , and not to | e | + | , so in the weak-driving case the correct dissipator should affect the population of the first combination, while leaving that of the second undamped. This means that there are undamped coherences in the | ± basis, which survive in the steady-state [26], and produce corrections to the heat currents relative to the results of the secular approximation.

Author Contributions

Conceptualization, methodology, software, analysis, investigation, C.N.M., L.T.T. and P.R.E. Writing—original draft preparation, P.R.E.; writing—review and editing, C.N.M. and L.T.T.; visualization, supervision, P.R.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Irish Research Council grant number GOIPG/2017/1091. Some calculations were performed on the Boyle cluster maintained by the Trinity Centre for High-Performance Computing. This cluster was funded through grants from the European Research Council and Science Foundation Ireland.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data generated in this study are available in Zenodo at https://dx.doi.org/10.5281/zenodo.5715836 (accessed on 18 November 2021).

Acknowledgments

We thank D. Segal for helpful comments on the manuscript, and G. Fux, B. Lovett, and J. Keeling for discussions and assistance with the TEMPO and PT-TEMPO codes.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Seletskiy, D.V.; Epstein, R.; Sheik-Bahae, M. Laser Cooling in Solids: Advances and Prospects. Rep. Prog. Phys. 2016, 79, 096401. [Google Scholar] [CrossRef] [PubMed]
  2. Nemova, G.; Kashyap, R. Laser Cooling of Solids. Rep. Prog. Phys. 2010, 73, 086501. [Google Scholar] [CrossRef] [Green Version]
  3. Epstein, R.; Sheik-Bahae, M. (Eds.) Optical Refrigeration: Science and Applications of Laser Cooling of Solids, 1st ed.; Wiley: Hoboken, NJ, USA, 2009. [Google Scholar] [CrossRef]
  4. Scovil, H.E.D.; Schulz-DuBois, E.O. Three-Level Masers as Heat Engines. Phys. Rev. Lett. 1959, 2, 262–263. [Google Scholar] [CrossRef]
  5. Geusic, J.E.; Schulz-DuBios, E.O.; Scovil, H.E.D. Quantum Equivalent of the Carnot Cycle. Phys. Rev. 1967, 156, 343–351. [Google Scholar] [CrossRef]
  6. Linden, N.; Popescu, S.; Skrzypczyk, P. How Small Can Thermal Machines Be? The Smallest Possible Refrigerator. Phys. Rev. Lett. 2010, 105, 130401. [Google Scholar] [CrossRef]
  7. Creatore, C.; Parker, M.A.; Emmott, S.; Chin, A.W. Efficient Biologically Inspired Photocell Enhanced by Delocalized Quantum States. Phys. Rev. Lett. 2013, 111, 253601. [Google Scholar] [CrossRef] [Green Version]
  8. Fruchtman, A.; Gómez-Bombarelli, R.; Lovett, B.W.; Gauger, E.M. Photocell Optimization Using Dark State Protection. Phys. Rev. Lett. 2016, 117, 203603. [Google Scholar] [CrossRef] [Green Version]
  9. Dorfman, K.E.; Xu, D.; Cao, J. Efficiency at Maximum Power of Laser Quantum Heat Engine Enhanced by Noise-Induced Coherence. Phys. Rev. E 2018, 97, 042120. [Google Scholar] [CrossRef] [Green Version]
  10. Kilgour, M.; Segal, D. Coherence and Decoherence in Quantum Absorption Refrigerators. Phys. Rev. E 2018, 98, 012117. [Google Scholar] [CrossRef] [Green Version]
  11. Murphy, C.N.; Eastham, P.R. Quantum Control of Excitons for Reversible Heat Transfer. Commun. Phys. 2019, 2, 120. [Google Scholar] [CrossRef] [Green Version]
  12. Brash, A.J.; Martins, L.M.P.P.; Barth, A.M.; Liu, F.; Quilter, J.H.; Glässl, M.; Axt, V.M.; Ramsay, A.J.; Skolnick, M.S.; Fox, A.M. Dynamic Vibronic Coupling in InGaAs Quantum Dots. J. Opt. Soc. Am. B 2016, 33, C115–C122. [Google Scholar] [CrossRef]
  13. Eastham, P.R.; Spracklen, A.O.; Keeling, J. Lindblad Theory of Dynamical Decoherence of Quantum-Dot Excitons. Phys. Rev. B 2013, 87, 195306. [Google Scholar] [CrossRef] [Green Version]
  14. Ramsay, A.J.; Gopal, A.V.; Gauger, E.M.; Nazir, A.; Lovett, B.W.; Fox, A.M.; Skolnick, M.S. Damping of Exciton Rabi Rotations by Acoustic Phonons in Optically Excited InGaAs/GaAs Quantum Dots. Phys. Rev. Lett. 2010, 104, 017402. [Google Scholar] [CrossRef] [Green Version]
  15. Gauger, E.M.; Wabnig, J. Heat Pumping with Optically Driven Excitons. Phys. Rev. B 2010, 82, 073301. [Google Scholar] [CrossRef] [Green Version]
  16. Abiuso, P.; Giovannetti, V. Non-Markov Enhancement of Maximum Power for Quantum Thermal Machines. Phys. Rev. A 2019, 99, 052106. [Google Scholar] [CrossRef] [Green Version]
  17. Thomas, G.; Siddharth, N.; Banerjee, S.; Ghosh, S. Thermodynamics of Non-Markovian Reservoirs and Heat Engines. Phys. Rev. E 2018, 97, 062108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Bylicka, B.; Tukiainen, M.; Chruściński, D.; Piilo, J.; Maniscalco, S. Thermodynamic Power of Non-Markovianity. Sci. Rep. 2016, 6, 27989. [Google Scholar] [CrossRef] [Green Version]
  19. Correa, L.A.; Palao, J.P.; Alonso, D.; Adesso, G. Quantum-Enhanced Absorption Refrigerators. Sci. Rep. 2015, 4, 3949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Nazir, A.; McCutcheon, D.P.S. Modelling Exciton–Phonon Interactions in Optically Driven Quantum Dots. J. Phys. Condens. Matter 2016, 28, 103002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Wrachtrup, J.; Jelezko, F. Processing Quantum Information in Diamond. J. Phys. Condens. Matter 2006, 18, S807–S824. [Google Scholar] [CrossRef]
  22. Norambuena, A.; Reyes, S.A.; Mejía-Lopéz, J.; Gali, A.; Maze, J.R. Microscopic Modeling of the Effect of Phonons on the Optical Properties of Solid-State Emitters. Phys. Rev. B 2016, 94, 134305. [Google Scholar] [CrossRef] [Green Version]
  23. Basilewitsch, D.; Cosco, F.; Lo Gullo, N.; Möttönen, M.; Ala-Nissilä, T.; Koch, C.P.; Maniscalco, S. Reservoir Engineering Using Quantum Optimal Control for Qubit Reset. New J. Phys. 2019, 21, 093054. [Google Scholar] [CrossRef]
  24. Breuer, H.P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002. [Google Scholar]
  25. Esposito, M.; Harbola, U.; Mukamel, S. Nonequilibrium Fluctuations, Fluctuation Theorems, and Counting Statistics in Quantum Systems. Rev. Mod. Phys. 2009, 81, 1665–1702. [Google Scholar] [CrossRef] [Green Version]
  26. Eastham, P.R.; Kirton, P.; Cammack, H.M.; Lovett, B.W.; Keeling, J. Bath-Induced Coherence and the Secular Approximation. Phys. Rev. A 2016, 94, 012110. [Google Scholar] [CrossRef] [Green Version]
  27. Hofer, P.P.; Perarnau-Llobet, M.; Miranda, L.D.M.; Haack, G.; Silva, R.; Brask, J.B.; Brunner, N. Markovian Master Equations for Quantum Thermal Machines: Local versus Global Approach. New J. Phys. 2017, 19, 123037. [Google Scholar] [CrossRef]
  28. Hartmann, R.; Strunz, W.T. Accuracy Assessment of Perturbative Master Equations: Embracing Nonpositivity. Phys. Rev. A 2020, 101, 012103. [Google Scholar] [CrossRef] [Green Version]
  29. Purkayastha, A.; Dhar, A.; Kulkarni, M. Out-of-Equilibrium Open Quantum Systems: A Comparison of Approximate Quantum Master Equation Approaches with Exact Results. Phys. Rev. A 2016, 93, 062114. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, J.; Segal, D. Coherences and the Thermodynamic Uncertainty Relation: Insights from Quantum Absorption Refrigerators. Phys. Rev. E 2021, 103, 032138. [Google Scholar] [CrossRef]
  31. Jeske, J.; Ing, D.J.; Plenio, M.B.; Huelga, S.F.; Cole, J.H. Bloch-Redfield Equations for Modeling Light-Harvesting Complexes. J. Chem. Phys. 2015, 142, 064104. [Google Scholar] [CrossRef] [Green Version]
  32. Kilgour, M.; Agarwalla, B.K.; Segal, D. Path-Integral Methodology and Simulations of Quantum Thermal Transport: Full Counting Statistics Approach. J. Chem. Phys. 2019, 150, 084111. [Google Scholar] [CrossRef] [Green Version]
  33. Boudjada, N.; Segal, D. From Dissipative Dynamics to Studies of Heat Transfer at the Nanoscale: Analysis of the Spin-Boson Model. J. Phys. Chem. A 2014, 118, 11323–11336. [Google Scholar] [CrossRef] [PubMed]
  34. Dümcke, R.; Spohn, H. The Proper Form of the Generator in the Weak Coupling Limit. Z. Phys. B 1979, 34, 419–422. [Google Scholar] [CrossRef]
  35. González, J.O.; Correa, L.A.; Nocerino, G.; Palao, J.P.; Alonso, D.; Adesso, G. Testing the Validity of the ‘Local’ and ‘Global’ GKLS Master Equations on an Exactly Solvable Model. Open Syst. Inf. Dyn. 2017, 24, 1740010. [Google Scholar] [CrossRef]
  36. Friedman, H.M.; Agarwalla, B.K.; Segal, D. Quantum Energy Exchange and Refrigeration: A Full-Counting Statistics Approach. New J. Phys. 2018, 20, 083026. [Google Scholar] [CrossRef]
  37. Trushechkin, A.S.; Merkli, M.; Cresser, J.D.; Anders, J. Open Quantum System Dynamics and the Mean Force Gibbs State. arXiv 2021, arXiv:2110.01671. [Google Scholar]
  38. Allen, L.; Eberly, J.H. Optical Resonance and Two-Level Atoms; Dover Publications: New York, NY, USA, 1987. [Google Scholar]
  39. Wangsness, R.K.; Bloch, F. The Dynamical Theory of Nuclear Induction. Phys. Rev. 1953, 89, 728–739. [Google Scholar] [CrossRef]
  40. Lindblad, G. On the Generators of Quantum Dynamical Semigroups. Commun. Math. Phys. 1976, 48, 119–130. [Google Scholar] [CrossRef]
  41. Gorini, V.; Kossakowski, A.; Sudarshan, E.C.G. Completely Positive Dynamical Semigroups of N-level Systems. J. Math. Phys. 1976, 17, 821–825. [Google Scholar] [CrossRef]
  42. Redfield, A.G. On the Theory of Relaxation Processes. IBM J. Res. Dev. 1957, 1, 19–31. [Google Scholar] [CrossRef]
  43. Bloch, F. Generalized Theory of Relaxation. Phys. Rev. 1957, 105, 1206–1222. [Google Scholar] [CrossRef]
  44. Silaev, M.; Heikkilä, T.T.; Virtanen, P. Lindblad-Equation Approach for the Full Counting Statistics of Work and Heat in Driven Quantum Systems. Phys. Rev. E 2014, 90, 022103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Gasparinetti, S.; Solinas, P.; Braggio, A.; Sassetti, M. Heat-Exchange Statistics in Driven Open Quantum Systems. New J. Phys. 2014, 16, 115001. [Google Scholar] [CrossRef]
  46. Popovic, M.; Mitchison, M.T.; Strathearn, A.; Lovett, B.W.; Goold, J.; Eastham, P.R. Quantum Heat Statistics with Time-Evolving Matrix Product Operators. PRX Quantum 2021, 2, 020338. [Google Scholar] [CrossRef]
  47. Strathearn, A.; Kirton, P.; Kilda, D.; Keeling, J.; Lovett, B.W. Efficient Non-Markovian Quantum Dynamics Using Time-Evolving Matrix Product Operators. Nat. Commun. 2018, 9, 3322. [Google Scholar] [CrossRef]
  48. Fux, G.E.; Butler, E.P.; Eastham, P.R.; Lovett, B.W.; Keeling, J. Efficient Exploration of Hamiltonian Parameter Space for Optimal Control of Non-Markovian Open Quantum Systems. Phys. Rev. Lett. 2021, 126, 200401. [Google Scholar] [CrossRef] [PubMed]
  49. Makri, N.; Makarov, D.E. Tensor Propagator for Iterative Quantum Time Evolution of Reduced Density Matrices. II. Numerical Methodology. J. Chem. Phys. 1995, 102, 4611–4618. [Google Scholar] [CrossRef]
  50. Gribben, D.; Rouse, D.M.; Iles-Smith, J.; Strathearn, A.; Maguire, H.; Kirton, P.; Nazir, A.; Gauger, E.M.; Lovett, B.W. Exact Dynamics of Non-Additive Environments in Non-Markovian Open Quantum Systems. arXiv 2021, arXiv:2109.08442. [Google Scholar]
  51. Tscherbul, T.V.; Brumer, P. Long-Lived Quasistationary Coherences in a V-type System Driven by Incoherent Light. Phys. Rev. Lett. 2014, 113, 113601. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Energy levels of an impurity in a model laser-cooling process. The two states of a ground-state manifold, | g u , and | g l , are coupled by the emission and absorption of lattice phonons (vertical solid lines). Laser driving occurs on the transition from the upper level of the ground-state manifold to an excited state | e (block arrow). This state decays radiatively to the ground state (wavy arrow).
Figure 1. Energy levels of an impurity in a model laser-cooling process. The two states of a ground-state manifold, | g u , and | g l , are coupled by the emission and absorption of lattice phonons (vertical solid lines). Laser driving occurs on the transition from the upper level of the ground-state manifold to an excited state | e (block arrow). This state decays radiatively to the ground state (wavy arrow).
Applsci 12 01620 g001
Figure 2. Rates of heat absorption from the phonon bath, as a function of the detuning δ = ω l E 0 of the driving laser from resonance, for four different Rabi frequencies. For each Rabi frequency, we show results computed using the full Bloch-Redfield equation (solid black curve), a phenomenological Lindblad equation (dashed orange curve), and the Bloch-Redfield equation in the secular approximation without the principal value terms (dashed black curve). The Rabi frequencies Ω are: (a) 0.01, (b) 0.1, (c) 0.5, and (d) 1.0.
Figure 2. Rates of heat absorption from the phonon bath, as a function of the detuning δ = ω l E 0 of the driving laser from resonance, for four different Rabi frequencies. For each Rabi frequency, we show results computed using the full Bloch-Redfield equation (solid black curve), a phenomenological Lindblad equation (dashed orange curve), and the Bloch-Redfield equation in the secular approximation without the principal value terms (dashed black curve). The Rabi frequencies Ω are: (a) 0.01, (b) 0.1, (c) 0.5, and (d) 1.0.
Applsci 12 01620 g002
Figure 3. Rate of heat absorption from the phonon bath, as a function of temperature, for the different approaches. Results are shown for weak driving, Ω = 0.1 , at resonance, δ = 0 . The heat current is computed using the full Bloch-Redfield equation (solid black curve), the phenomenological Lindblad equation (dashed orange curve), and the Bloch-Redfield equation in the secular approximation without the principal value terms (dashed black curve).
Figure 3. Rate of heat absorption from the phonon bath, as a function of temperature, for the different approaches. Results are shown for weak driving, Ω = 0.1 , at resonance, δ = 0 . The heat current is computed using the full Bloch-Redfield equation (solid black curve), the phenomenological Lindblad equation (dashed orange curve), and the Bloch-Redfield equation in the secular approximation without the principal value terms (dashed black curve).
Applsci 12 01620 g003
Figure 4. Heat absorption rates, as a function of the detuning, computed using the full Bloch-Redfield equation (solid black curves) and a numerically exact method (red squares), for two different Rabi frequencies. The Rabi frequencies Ω are (a) 0.5 and (b) 1.0.
Figure 4. Heat absorption rates, as a function of the detuning, computed using the full Bloch-Redfield equation (solid black curves) and a numerically exact method (red squares), for two different Rabi frequencies. The Rabi frequencies Ω are (a) 0.5 and (b) 1.0.
Applsci 12 01620 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Murphy, C.N.; Toledo Tude, L.; Eastham, P.R. Laser Cooling beyond Rate Equations: Approaches from Quantum Thermodynamics. Appl. Sci. 2022, 12, 1620. https://doi.org/10.3390/app12031620

AMA Style

Murphy CN, Toledo Tude L, Eastham PR. Laser Cooling beyond Rate Equations: Approaches from Quantum Thermodynamics. Applied Sciences. 2022; 12(3):1620. https://doi.org/10.3390/app12031620

Chicago/Turabian Style

Murphy, Conor N., Luísa Toledo Tude, and Paul R. Eastham. 2022. "Laser Cooling beyond Rate Equations: Approaches from Quantum Thermodynamics" Applied Sciences 12, no. 3: 1620. https://doi.org/10.3390/app12031620

APA Style

Murphy, C. N., Toledo Tude, L., & Eastham, P. R. (2022). Laser Cooling beyond Rate Equations: Approaches from Quantum Thermodynamics. Applied Sciences, 12(3), 1620. https://doi.org/10.3390/app12031620

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop