Perturbative dissipation dynamics of a weakly driven Jaynes-Cummings system

Perturbative dissipation dynamics of a weakly driven Jaynes-Cummings system

S. M. Yu School of Applied Sciences, Beijing Univertiy of Technology, Beijing, 100124, China    Y. B. Gao School of Applied Sciences, Beijing Univertiy of Technology, Beijing, 100124, China    H. Ian Institute of Applied Physics and Materials Engineering, University of Macau, Macau UMacau Research Institute, Zhuhai, Guangdong, China

We generalize a microscopic master equation method to study the dissipation dynamics of Jaynes-Cummings two-level system with a weak external driving. Using perturbative analysis to extend the damping bases theory, we derive the corrected Rabi oscillation and vaccum Rabi splitting analytically. The evolution of the decoherence factor of the weakly driven system reveals that the off-diagonal density matrix elements are oscillating at a frequency dependent on the driving strength and the initial population inversion. For highly-inverted systems at the weak-driving limit, this frequency reduces to twice the value for the non-driven system, showing the dissipation dynamics unable to be discovered using more conventional approaches.

I introduction

The celebrated Jaynes-Cummings (JC) model jaynes63 () predicts the splitting of the energy levels of two-level atoms under their coupling to a quantum field and provide an analytical basis for the full-quantum oscillation dynamics of two-level systems. The simple yet elegant model has paved the way to the much ramified research subject of cavity quantum electrodynamics (cQED) Raimond01 () for the atom-field interaction in optical cavities. Recently, it is incarnated in the cQED studies of superconducting circuit cavities wallraff04 (); majer07 () where the two-level system is a superconducting qubit comprising two macroscopic persistent-current states, which has contingent applications in quantum logic gate implementation devoret13 () for quantum computation divincenzo00 ().

It is well-known that one of the major challenges in realizing a functional quantum computer is the decoherence inherent to a superconducting qubit, whose sources are the multiple environmental couplings such as the electron-phonon coupling and the electron-electron coupling that results in 1/f noises and phase noises innate to a solid-state system. Therefore, it is imperative to describe the qubit dynamics that accounts for the decoherence effects more accurately under the JC-model.

Earlier studies that accounts the decoherence effects employ the so-called phenomenological master equation that extends the quantum Liouvill equation for the density matrix to include terms consisting of products of the density matrix and operators from a semigroup. It is the simplest implementation of Lindblad’s modellindblad76 () for density matrix evolution in the sense that the coefficients of these terms responsible for decoherences are assumed a uniform constant extracted phenomenologically from experiments. Later, progresses in the direction of nonequilibrium thermodynamics davies74 (); spohn78 () has accumulated to the Kubo-Martin-Schwinger condition that permits the depiction of a state-dependent decay coefficient for different terms in the master equation, where is the eigenfrequency associated with the operators in the semigroup for decoherence dynamics. Consequently, the set of decoherence operators can be expanded under an exact basis of the environmental Hilbert space pertinent to a particular form of heat-bath interaction Hamiltonian, giving rise to a precise microscopic master equation derived exactly from a given environmental coupling BREUER ().

Here, we employ the microscopic master equation to study the quantum dynamics of a cavity QED system with the cavity mode coupled to a two-level system under the JC-model. The cavity itself is weakly driven by a resonant laser field. The decay dynamics without an external driving previously studied by Scala et al. shows that the populations in the quantized levels are oscillating around an averaged exponential curve derived from the phenomenological theory. When the cavity is weakly driven, we amend the exactly diagonalizing eigenvectors with perturbative terms to derive a set of nine damping bases for the Liouville superoperator, under which decay sources previously underdiscovered gao13 (); gao16 (); gao09 () are found. Using these damping bases, we find minute-scale oscillations of higher frequency on top of the original Rabi oscillations eberly80 (); eberly81 (); rempe87 () that is envelopped by a decay due to the relaxation of the two-level system. In addition, the splitting of the dressed levels at vacuum is contracted compared to the non-driven case. More importantly, we observe the decoherence factor would follow different oscillation frequencies under different initial population inversions, which can be explained by the distinct damping bases dominating under particular inversion cases.

We first derive the microscopic master equation and the associated eigenoperators in Sec. II. Using these operators as bases, we compute the correlation functions of the cavity mode quadrature to show the effective decay rates and the noise density spectrum in Sec. III. The evolution of the level populations is shown in Sec. IV and conclusions are given in Sec. V.

Ii microscopic master equation

We start the discussion with the model Hamiltonian where ()


stands for the JC-model Hamiltonian for one two-level system of eigenfrequency and one cavity quantum field of frequency under the rotating wave approximation for their interaction. The weak external driving is represented by the Hamiltonian


Under weak-excitation, we consider only the lowest three eigen-levels for at atom-cavity resonance: and . Further considering as a weak perturbation to , we can find the perturbed eigen-basis vectors


in terms of the bare basis vectors under the first-order perturbation expansion. The perturbation contributes a quadratic term correction to the eigenvalues from the driveless JC-model, which gives


Pairing the vectors of Eqs. (3)-(4) with their conjugates form a set of nine damping bases briegel93 (); scala07 () for the Liouville operator of the master equation BREUER ()


where we have split the environmental contribution into two parts with each having the effect


on the density matrix. The damping rate obeys the Kubo-Martin-Schwinger (KMS) condition for the quasi-equilibrium system at temperature . Using and for all other to denote the nine damping bases, the corresponding eigenvalues associated with the eigen-equation read


where the abbreviation has been adopted. We note that since the damping basis has a zero eigenvalue, it is obviously identical to the steady-state density matrix for . Thus, for any initial state expanded in the damping bases, the master equation (7) admits the formal solution


Iii Corrected Rabi oscillation and vacuum splitting

The damping bases with their associated eigenvalues (9) to the master equation (7) leads to a decaying envelop to the Rabi oscillation of the two-level system. The effect of the weak driving can be examined in the temporal domain where we assume the two-level system is fully inverted initially, i.e. . In the dressed state bases, one can expand the state to first-order in the perturbative series as


making the corresponding density matrix be expandable in the damping bases as


The excited-state population can then be computed from Eq. (10) as a function of the driving strength ,


where the first term indicates the spontaneous decay from the two dressed levels and while the second term indicates the Rabi oscillation of population. The Rabi frequency is corrected from its original for a bare system by a correction term due to the driving, which makes the oscillation slower. In addition, the driving leads to a third term with , which breaks the oscillation symmetry of the dressed levels.

Figure (1) shows the comparison of the evolution over a dimensionless time scale between the non-driven systems and the driven systems, using the experimental parameters typical of a superconducting qubit where GHz and , and . As one expects from Eq. (13), the non-driven case has follow an enveloped Rabi oscillation of frequency contributed by the second term. The third term proportional to squared driving strength has the effect of minute oscillations of higher frequencies and imposed on top of the Rabi oscillation, which are themselves enveloped by the decay rates and respectively. Note that the summation of the two sinusoidals has the apparent effect that the minute oscillation has a carrier frequency of and an amplitude oscillation frequency of . The first term of Eq. (13) contributes to a non-oscillating offset, which is only visible in the plot over long-term.

Figure 1: (a) Plot of the excited state population of a bare two-level system with no external driving, i.e. , over dimensionless time. (b) The change over when the system is weakly driven with driving strength (blue curve) and (red curve). 

To find the spectrum of the driven dressed system, we employ the so-called input-output theory, whereby the photon statistics of a probed output field is dependent on the noise input through a correlation Green function of two operators. The time-dependent operator is meant to represent and the noise density spectrum is the Fourier transform


For the cavity mode considered in the JC-model, we have so that the resulting correlation functions are essentially the perturbation-corrected fluctuation-dissipation (FD) relations of the single-mode field under the framework of microscopic master equation. At the first-order perturbative expansion given by Eqs. (3)-(4), is also an eigenvector of , the corrections contributed by carrying a negligible coefficient on the order of . Consequently, when we consider the system evolution begin with the equilibrium state , it becomes easy to find that


where is a dimensionless driving-dependent quantity, which vanishes when the perturbation is taken to its weak-driving limit. Following the same approach, we can find


recovering the familiar FD-relations for the cavity vacuum under the normal approach to dissipation dynamics using phenomenological theories. The difference stemming from the microscopic approach lies in the correlation because is not an eigenvector of , rather




where we have omitted the contribution from in the second equality since its coefficient also vanishes in the weak-driving limit, similar to that of Eq. (15).

Collecting the four FD-relations above, we have the correlation function for the real quadrature identical to Eq. (18), i.e. composing of two exponentials. Then it is straightforward to apply Eq. (14) to find each exponential contribute a Lorentzian spectral line in the frequency space, i.e.


Figure 2 shows a plot of the spectrum where we have adopted an ohmic dissipation spectrum leggett87 ()


with cutoff frequency and dissipation coefficient . Then the damping channels along the two levels and become evidently unequal. With and set to 1 GHz, we have , making the right Lorentzian at in the plot exhibit a higher peak and a wider half-width than the left one at .

Figure 2: (a) Plot of the spectrum for the weak-driving scenario (red dashed) and the non-driven scenario (blue solid), showing two peaks of uneven heights and widths at and . Magnified views of the left and the right peaks are shown in (b) and (c), respectively. 

The driving results in red-shifts of the two Lorentzians towards to the left-end of the plot with unequal shift amplitude. As a whole, the driving as a perturbation has the effect of reducing the vacuum Rabi split, following the formula


Iv Decoherence in Driven Jaynes-Cummings Model

The microscopic master equation approach to the driven JC model also gives a more accurate prediction of the decoherence dynamics of a qubit or two-level system in general. Consider an arbitrary initial state of the qubit and a vacuum state for the resonator, i.e. where we let and be real and the complex phase be accounted by the phase factor . Under the truncated dressed bases, the initial state becomes


The corresponding pure-state density matrix can be decomposed into the damping bases like we did in last section and here the coefficients are


, and .

Consequently, through the formal solution Eq. (10), we can find the off-diagonal element of the density matrix to be


for calculating the decoherence factor  zurek05 (). Note that when the driving vanishes, the terms stemmed from the perturbation in the equation above vanishes and the decoherence factor reduces to gao16 ()


showing that the decoherence in the absence of driving is insensitive to the initial population distribution in the qubit.

Figure (3) shows the comparison of decoherence factors of the driven cases to the non-driven case, which demonstrates the inadequacy of the non-driven model to describe the relaxation dynamics of a driven model where the initial population distribution plays a crucial role, particularly when the system is strongly coupled with . The parameters are typical of superconducting qubits with GHz and . The decay rates and are chosen unequal to reflect the uneven damping channels they associated with for a given environmental reservoir, such as that of Eq. (20). The influence of population distribution is shown using three different ratios of at (equal distribution, blue curve), (ground populated, red curve), and (excited-state populated, yellow curve).

Figure 3: (a) Plot of the decoherence factor as a function of time in a non-driven two-level system. (b) Plots of the decoherence factor under strong coupling and weaking driving when the system is populated at different proportions: (red curve), (blue curve), and (yellow curve). 

We observe that only when the two-level qubit is initially equally populated would the decoherence factor resemble that of a non-driven case. This shows that the driving has minimal effect on a thermally equilibrated system. On the other hand, when the system is initially either inverted or condensed to ground, the driving would set the decoherence into a more complicated oscillation cycle. From the case, one can see that the grounded population first experiences a much larger decoherence, about twice, and undergoes a much quicker decay before it converges with the curve of the non-driven case. It is also noticeable that the offset, which is mainly due to the later terms in , becomes nonzero and the phase of the decoherence oscillation is inverted from the non-driven case. Since is small, the oscillation is still dominated by the terms led by the coefficients and , making the frequency of similar to that of given in Eq. (28), i.e. at .

Whereas for the inverted case with large, we observe from Eqs. (24)-(26) that is dominated by and terms. Then from their associated eigenvalues and of the damping bases given in Eq. (9), the oscillation frequency becomes close to . This explains the behavior of the yellow curve for in Fig. 3.

Moreover, one sees from Eqs. (24)-(26) that the complex phase factor of the initial state also affects the decoherence. Fig. 4 shows the change in the decoherence factor from the non-driven case when the phase is set to five values from to . We can discern clearly a period at in all subplots and the inversion of decoherence between a half-period difference, that is .

Figure 4: Plots of the change in the decoherence factor due to different initial phases: (a) , (b), (c) , (d) , and (e)

V Conclusions and discussions

We have developed a weakly driven Jaynes-Cummings model under the microscopic master equation. The effects of the weak drivings are given as second-order perturbative terms in the eigenvalues associated with the damping bases of the Liouville superoperator. The exemplified effects due to the weak driving include corrected Rabi oscillations and contracted vacuum Rabi splitting. More noticeable is the sensitivity of the decoherence factor to initial population of the two-level system given the external driving, which were unaccounted for in previous non-driven JC-models. Especially when the system is initially inverted, the decoherence factor would oscillating at twice the Rabi frequency, showing the decoherence dynamics is dominated by a set of damping bases different from the usual set for its non-driven counterpart. This finding demonstrates that the microscopic master equation approach is not mere corrections to the simpler Langevin equation approach and Maxwell-Bloch equation approach for analyzing the relaxation process of a quantum system. On the contrary, it provides a better ground for analysis that leads to previous undiscovered phenomena.

H. I. acknowledges the support of FDCT Macau under grant 013/2013/A1, University of Macau under grant MRG022/IH/2013/FST and MYRG2014-00052-FST, and the National Natural Science Foundation of China under Grant No. 11404415.


  • (1) E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963).
  • (2) J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. Phys. 73, 565 (2001).
  • (3) A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 431, 162 (2004).
  • (4) J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, A. Wallraff, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, Nature 449, 443 (2007).
  • (5) M. H. Devoret and R. J. Schoelkopf, Science 339, 1169 (2013).
  • (6) D. DiVincenzo, Fortschr. Phys. 48, 771 (2000) .
  • (7) G. Lindblad, Commun.Math. Phys. 48, 119 (1976).
  • (8) E. B. Davies, Commun.Math. Phys. 39, 91 (1974).
  • (9) H. Spohn and J. L. Lebowitz, Adv. Chem. Phys. 38, 109 (1978).
  • (10) Heinz-Peter Breuer and Francesco Petruccione, The Theory of Open Quantum Systems, (Oxford Press, Oxford 2002).
  • (11) C. Chen and Y. B. Gao, Commun. Theor. Phys. 60, 531 (2013).
  • (12) X. Xiao, M. Y. Zhao, S. M. Yu and Y. B. Gao, Commun. Theor. Phys. 65, 273 (2016).
  • (13) Y. B. Gao, S. Yang, Yu-xi Liu, C. P. Sun, and Franco Nori, arxiv: 0902.2512.
  • (14) J. H. Eberly, N. B. Narozny, and J. J. Sanchez- Mondragon, Phys. Rev. Lett. 44, 1323 (1980).
  • (15) N. B. Narozny, J. J. Sanchez-Mondragon, and J. H. Eberly, Phys. Rev. A. 23, 236 (1981).
  • (16) G. Rempe, H. Walther, and N. Klein, Phys. Rev. Lett. 58, 353 (1987).
  • (17) H.-J. Briegel and B.-G. Englert, Phys. Rev. A 47, 3311 (1993).
  • (18) M. Scala, B. Militello, A. Messina, J. Piilo, S. Maniscalco, Phys. Rev. A 75, 013811 (2007).
  • (19) A. J. Leggett, S. Chakravarty, A. T. Dorsey, Matthew P. A. Fisher, Anupam Garg, and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987).
  • (20) F. M. Cucchietti, J. P. Paz, and W. H. Zurek, Phys. Rev. A 72, 052113 (2005).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description