Timedependent nonlinear JaynesCummings dynamics of a trapped ion
Abstract
In quantum interaction problems with explicitly timedependent interaction Hamiltonians, the time ordering plays a crucial role for describing the quantum evolution of the system under consideration. In such complex scenarios, exact solutions of the dynamics are rarely available. Here we study the nonlinear vibronic dynamics of a trapped ion, driven in the resolved sideband regime with some small frequency mismatch. By describing the pump field in a quantized manner, we are able to derive exact solutions for the dynamics of the system. This eventually allows us to provide analytical solutions for various types of timedependent quantities. In particular, we study in some detail the electronic and the motional quantum dynamics of the ion, as well as the timeevolution of the nonclassicality of the motional quantum state.
I Introduction
The verification and quantification of nonclassical effects, that is phenomena which cannot be explained by Maxwell’s equations, is a main concern of theoretical and experimental quantum optics. Many of those effects, like squeezing (1); (2); (3); (4); (5), entanglement (6); (7), and photon antibunching (8), were intensively investigated over many decades. However, there are effects beyond this set, like for example anomalous quantum correlations (9); (10); (11), which arise from the violation of fieldintensity inequalities. In such and related scenarios a subject of interest is the investigation of the interplay of free fields and fields which are attributed to sources, which play an important role in the theory of spectral filtering of light (12); (13). The relationship between field correlation functions of freefield and sourcefield operators were, for example, considered in Refs. (14); (15). Hence, the treatment of a physical system containing contributions from both kinds of fields is an interesting aspect to be studied, especially when the corresponding dynamics is exactly solvable. A suitable model for this purpose is the JaynesCummings model, which contains not only freefield parts but a sourceattributed part as well. In the following we will briefly reconsider its history and possible areas of application.
When the JaynesCummings model was proposed in 1963 (16); (17), its practical relevance was doubted, as it describes an idealized scenario of the resonant interaction of a twolevel system with only a single radiation mode. However, in the 1980s the model’s importance vastly enhanced, since, due to technical progress, it was possible to experimentally prove many of its predictions (18); (19); (20); (21). Remarkably, despite its simplicity, the JaynesCummings model exhibits a plenty of physical effects, e.g. Rabi oscillations (22); (23); (24), collapse and revivals (21); (25); (26), squeezing (27); (28), atomfield entanglement (29); (30); (31), antibunching (32); (33); (34), and nonclassical states such as Schrödinger cat (35); (36) or Fock states (37); (38); (39). Initially intended for describing the interaction of a single atom with a single radiation mode, the JaynesCummings model could be applied to a variety of physical scenarios. Examples are Cooperpair boxes (40); (41), ”flux” qubits (42), and Josephsonjunctions (43); (44); (45). It can also be applied in solid state systems to describe the (strong) coupling of qubits to a cavity mode, for example in quantum dots (46); (47); (48); (49) or superconducting circuits (50); (51); (52). Another recent application of the JaynesCummings model is the description of Rydbergblockaded atomic ensembles (53); (54).
The JaynesCummings model also became relevant for the vibronic dynamics of trapped ions, where the quantized mode of the electromagnetic field is replaced by the quantized centerofmassmotion of the ion (55); (56); (57). Later on, a nonlinear JaynesCummings model was introduced (58), which describes the dynamics of a trapped ion beyond the standard LambDicke regime (59); (60). The motional degrees of freedom are coupled to the electronic states of the ion by a classical pump field in the resolved sideband regime. On this basis it became possible to generate many motional states for trapped ions, such as Fock states, squeezed states (61); (59), even and odd coherent states (62), nonlinear coherent states (63), pair coherent states (64), superpositions of the latter (65), SU(1,1) intelligent states (66), Schrödingercat states (67); (68), entangled coherent states (68), and generalized Kerrtype states (69). As for the standard JaynesCummings Hamiltonian, see for example Ref. (70), the trapped ion dynamics based on the nonlinear JaynesCummings model was also considered beyond the rotating wave approximation (71); (72); (73); (74).
In the present paper we study the vibronic nonlinear JaynesCummings model, when the classical driving laser field is slightly detuned from the th sideband. Such a mismatch yields an explicitly timedependent Hamiltonian in the Schrödinger picture, whose corresponding dynamics is not easily solved due to the relevance of timeordering effects. We will demonstrate that these difficulties can be resolved by extending the Hilbert space of the problem to include the driving field in the quantum description. On this basis, the full dynamics will become exactly solvable. This renders it possible to study sophisticated problems of explicitly timedependent dynamics on the basis of the exactly solvable extended problem. This yields deeper insight in the yet rarely studied quantum dynamics in cases when explicitly timedependent Hamiltonians and the resulting timeordering prescriptions are relevant. Also the quantum effects and the nonclassical correlation properties of the system can be studied by this method in great depth.
The paper is structured as follows. Section II introduces the Hamiltonian used in this paper and briefly discusses its physical meaning as well as timeordering effects. In Sec. III, we solve the dynamics using the eigenstates of the generalized Hamiltonian. Afterwards, in Sec. IV, we use the regularized GlauberSudarshan function to study the nonclassical evolution of the motional quantum properties of the ion. Finally, a summary and some conclusions follow in Sec. V.
Ii Explicitly timedependent nonlinear JaynesCummings model
The nonlinear JaynesCummings model (NJCM) for the vibronic coupling between the electronic and motional degrees of freedom of a trapped ion was introduced for the situation of the exactly resonant interaction of a laser field with the th vibronic sideband of the ion (58); (75). This interaction Hamiltonian can be exactly diagonalized. It was experimentally demonstrated that it properly describes the dynamics of a trapped ion for the case of (59). In the present paper we are interested in more sophisticated timedependent quantum phenomena of such a system. For this reason, in a first step we generalize the NJCM to allow for the explicitly timedependent dynamics. In this case, however, an exact solution of the problem seems to be not feasible and numerical solutions are required.
ii.1 Explicitly timedependent Hamiltonian
Let us start with the following Hamiltonian, describing an ion, trapped in a harmonic trap potential, interacting with a classical laser field, see (58) and Chap. 13 of (75),
(1)  
(2) 
describes the free motion of the vibrational centerofmass and electronic degrees of freedom of the twolevel ion, with the vibrational frequency and the electronic transition frequency . The laser is assumed to be monochromatic and quasiresonant with the electronic transition, , and to have only one nonvanishing wave vector component, such that only one motional degree of freedom appears in the interaction term. The complex amplitude describes the pump laser. The operators () are the creation (annihilation) operators of the vibrational frequency . The electronic flip operators () describe the atomic transitions, is a projection of the electricdipole matrix element on the direction of the electrical field and describes the mode structure of the pump laser at the operatorvalued position of the ion. For a standing wave it reads as
(3) 
where defines the relative position of the trap potential to the laser wave. The LampDicke parameter describes the effects of momentum transfer on the atomic wave packet due to recoil effects.
Applying the BakerCampbellHausdorff formula in together with a power series expansion, we get
(4) 
The interaction Hamiltonian in the interaction picture (indicated by the tilde) reads as
(5) 
If the laser is exactly resonant to the th vibrational sideband, this yields the exactly solvable nonlinear JaynesCummings interaction (58); (75).
For the purpose of the present paper, we are interested in the situation when the laser is slightly detuned from the th sideband:
(6) 
with . We still assume that the ion is in the resolvedsideband limit, i.e., we can resolve the single sidebands very well. This means that the linewidths of the vibronic transitions and the coupling strength are small compared to the vibrational frequency . In this case we can address vibronic transitions of equal transition frequencies. Here we consider the transitions for , cf. Eq. (6). Hence, we perform a vibrational rotatingwave approximation,
(7) 
in Eq. (5). This yields the interaction Hamiltonian
(8) 
where
(9) 
with denoting the generalized Laguerre polynomials. The Hamiltonian (8) describes the nonlinear th sideband coupling of the vibrational mode and the electronic transition, , cf. Fig. 1. It is important that this nonlinear interaction is explicitly timedependent, as long as .
ii.2 Solution of the timedependent interaction and timeordering effects
The most general solution of the time evolution of a quantum system, described by its Hamiltonian , is given by its timeevolution operator,
(10) 
where denotes the timeordering prescription. The latter accounts for the temporal order of the Hamiltonians with different time arguments contained in the exponential function. If the Hamiltonian, however, is not explicitly timedependent or commuting with itself at different times, the ordering symbol becomes superfluous and the standard exponential power series of is recovered.
Alternatively, the representation (10) can also be given in the form of the Magnus expansion (76); (77),
(11) 
which is unitary in each order , with . Herein, the contributions of are referred to as timeordering effects or timeordering corrections (78); (79); (80); (81); (82). However, the contain multiple integrals of nested commutators of the Hamiltonians at different times which are, especially for higher orders, difficult to handle. A possibility to circumvent this problem was presented in Ref. (83), where only one commutator needs to be evaluated. However, in this representation a needed diagonalization of the operatorvalued problem is not trivial, as the different orders of the expansion do not necessarily possess a common eigenbasis. For certain physical models and regimes, the timeordering symbol in (10) can be neglected, for example for parametric downconversion with not too high pump powers (78). Hence, let us begin to study the influence of timeordering effects on the dynamics described by the Hamiltonian in Eq. (8).
For this purpose we use the opensource software package Qutip (84); (85) in Python to obtain numerically the timeordered solutions based on Eq. (10) together with the Hamiltonian (8). To visualize the effects of time ordering, we compare the solutions with those when the time ordering is discarded, ,
(12) 
In this case the integral can be directly evaluated,
(13) 
For convenience we introduce the dimensionless quantities
(14) 
such that
(15) 
Furthermore, we define . Hence, we have
(16) 
and we assume from now on.
The eigenstates of the integrated Hamiltonian (II.2) read as
(17) 
with denoting the electronic () and motional () excitations, see, e.g., Chap. 12 of Ref. (75). Due to normalization we find immediately . These states are often referred to as “dressed states”. Solving
(18) 
yields the parameters
(19) 
where , cf. Eq. (9). The completeness relation of these states reads
(20) 
This yields the time evolution operator (16) in the form
(21) 
since the part cancels, as for .
For further investigations let us consider the population probability of the excited electronic state, which was studied only for in Ref. (58),
(22) 
which is given now in dependence on the scaled time . For the visualization we chose the input state at . The atom is initially prepared in the electronic ground state and the motional state of the ion is a coherent state. Details concerning the coherent state preparation of the motional state of the ion can be found in Refs. (59); (86). This eventually yields
(23) 
Note that the contributions cancel each other. The temporal evolution of is depicted in Fig. 2.
The correct numerical solution significantly differs from the analytical one without time ordering. That is, neglecting the timeordering effects, even for a very small frequency mismatch , strongly falsifies the electronic population dynamics. Hence, the ordering plays an important role and it must not be omitted.
Iii Nonlinear JaynesCummings Model with quantized pump
In this section we will overcome the shortcoming of the nonlinear JaynesCummings model with frequency mismatch by quantization of the pump field. In practice, this can be realized by placing the trapped ion in a high cavity. Now a mode of the quantized cavity field pumps the vibronic transition. We will see that this extension of the Hilbert space allows us to exactly solve the full interaction problem. This opens new possibilities to study problems underlying time ordering analytically, which are not solvable in a semiclassical approach.
iii.1 Quantization of the pump field
As shown above, the explicit time dependence of the Hamiltonian (8) prevents an analytical solution of the dynamics as we cannot discard – even approximately – the timeordering effects. Hence, our aim is to eliminate the timedependence in the Hamiltonian . For this purpose we return to the Hamiltonian (2) and quantize the pump field by replacing
(24) 
where is the annihilation operator of the pump quanta in the Schrödinger picture. The semiclassical timedependence is thus transformed into the free evolution of the operator . In practice this can be realized via “Cavity Quantum Electrodynamics (QED)” with a trapped ion, for theory and corresponding experiments, see (87); (88) and (89); (90), respectively.
The total Hamiltonian in the Schrödinger picture, including the quantized pump field, reads as
(25)  
(26) 
which is now timeindependent. The modified free Hamiltonian, , includes the free evolution of the quantized pump field with the frequency . Here, we again assume that we operate in the resolved sideband regime and we still consider the transitions. As before, only those terms of are relevant which belong to this transition, cf. Eq. (4),
(27) 
with being defined in Eq. (9). Thus, we arrive at
(28) 
The interpretation of the Hamiltonian (28) is the following: A pump photon is absorbed () and the ion is excited (). The vibrational transitions () occur according to our chosen quasiresonance condition, . The Hermitian conjugate (H.c.) term in addition describes the emission of a pump photon (), accompanied by the electronic transition and the vibrational transition . As before, the pump field is not exactly on resonance, as . In the case of interest, , only the wanted transitions significantly contribute to the dynamics. Since the resulting Hamiltonian is not explicitly timedependent anymore, we obtain the timeevolution operator in the form
(29) 
with the definitions of the Hamiltonian according to Eqs. (25) and (28).
iii.2 Solution for the quantized pump field
In this section we solve the timeevolution problem based on the operator (29), i.e. we derive an analytical expression for . Our full Hamiltonian, cf. Eqs. (25) and (28), reads as
(30) 
The eigenstates of this Hamiltonian are
(31) 
where denotes the electronic, pumpphoton (), and motional excitations. The normalization yields . The general procedure resembles that in Subsec. II.2. However, now we have an additional mode and we will solve the problem in the Schrödinger picture.
The parameters are found to be
(32) 
where we used Eq. (6) and defined
(33) 
Here, are the eigenvalues of , associated with the eigenstates , and is the nonlinear quantum Rabi frequency, which was already discussed in Ref. (58). For the present problem there occurs in Eq. (33) the additional factor , which is caused by the quantum treatment of the pump field.
The completeness relation reads
(34) 
Using the latter, we can rewrite the full timeevolution operator, Eq. (29), in the form
(35) 
with . For convenience we use the scaled (dimensionless) parameters:
(36) 
with . In terms of these dimensionless quantities the unitary time evolution operator reads as
(37) 
where
(38) 
iii.3 Semiclassical versus quantized pump
Let us now compare the analytical results obtained from the solution of the quantizedpump dynamics with the numerical solutions, see Sec. II.2, for a classical pump field, when the Hamiltonian is explicitly timedependent, cf. Eq. (8). As an example, we calculate the occupation probability of the excited electronic state,
(39) 
where is the full density matrix of the state. Using the input state at , the analytical treatment yields
(40) 
Note that, due to the dependence on , there is no dependence of on and .
Let us consider the evolution for a relatively weak pump amplitude, cf. Fig. 3. We see that, excluding the shorttime dynamics, the solutions with classical and quantized pump differ significantly from each other. Hence, the used pump amplitude is by far not sufficiently large to be referred to as “quasiclassical”.
In the case of a stronger pump field, cf. Fig. 4, we indeed obtain a dynamics which is almost identical to the numerical solution for a classical pump field. This enables us not only to conclude which pump amplitudes are needed such that the pump can be treated as a classical one on the corresponding time scale. It also reveals, that the solution found via quantization of the pump yields a more general description of the quantum system under study, where the time ordering is contained via the extension of the Hilbert space
Iv Evolution of nonclassicality
In this section we use the solution, Eq. (III.2), to discuss the nonclassical properties of the system. Let us first introduce the notion of nonclassicality to be used in the following.
iv.1 The regularized function
Using the GlauberSudarshan function (91); (92), , any quantum state, given by its density operator at time , can be expressed as a mixture of coherent states ,
(41) 
We call a state nonclassical, if it cannot be expressed as a classical mixture of coherent states. In such cases cannot be interpreted in terms of a classical probability density (93); (94), i.e., it can attain negative values in the sense of distributions. However, for many states is highly singular and, hence, it is not accessible in experiments. To uncover the negativities of it is therefore necessary to use a regularization procedure which yields a regularized version of this function (95). This procedure was successfully applied to experimental data (96); (97); (98) and generalized to different scenarios (99); (100); (101). Here we will only recapitulate the basic idea.
The function is defined by the Fourier transform of the characteristic function with ,
(42) 
The possibly occurring singular behavior of results from the fact that may be unbounded and, hence, not squareintegrable. According to Eq. (42), the function can therefore be highly singular. To get experimental access to the latter, one may introduce a filter function with some filter width , to define the regularized function (95) as
(43) 
The resulting function is a regular and smooth (99) function as long as the following requirements to the filter function are fulfilled:

can be Fourier transformed for all filter widths , with .

The Fourier transform of is a probability density, so that it is nonnegative.

For a filter which is infinitely broad, , we obtain the original function, .
For an overview and the discussion of different filter functions we refer to Ref. (102).
iv.2 Calculation of in Fock basis
In practical calculations, to obtain the full information on the quantum state, the density matrix of the state is calculated. In the following, we implement a suitable procedure to calculate directly out of . Let us first rewrite the definition of , cf. Eq. (43), in Fock state basis,
(44) 
with and . The functions are the regularized elements of the function in the Fock basis. In Eq. (IV.2), the complete time evolution is contained in the density matrix elements . The functions however, only depend on the fixed parameter and the phasespace coordinate . Hence, it is possible to calculate these elements only once and after that we apply them to , for arbitrary . Let us therefore find a suitable expression of .
We make use of (75)
(45) 
with
(46) 
where are the generalized Laguerre polynomials. Using a radial symmetric filter function (102), with for and , then may be rewritten as
(47) 
The phaseintegral can be evaluated via substitution of the limits of integration
(48) 
where are the Bessel functions of the first kind.
Finally, we arrive at the expression
(49) 
This relation holds true for all radial symmetric filters . We will use the filter (102)
(50) 
with if and elsewhere. Inserting Eq. (50) in Eq. (IV.2) and using the substitution yields
(51) 
The integral needs to be evaluated numerically in general. Note that here holds. We stress that this procedure applies to any time evolution.
iv.3 Nonclassicality in the nonlinear JaynesCummings model
We are interested in the nonclassical properties of the vibrational states. Hence we calculate the reduced density matrix
(52) 
where the trace over the electronic states and the pump states is evaluated. Here we use with the timeevolution operator given in Eq. (III.2) and at . This yields
(53) 
which does not depend on .
The surface plot of the regularized function is given in Fig. 5. Due to the vibronic coupling the initial motional coherent state at evolves into a nonclassical state. For a rather small time, cf. Fig. 5 (a), the state is still close to a coherent one. For larger times, one obtains more distorted states, cf. Figs. 5 (b) and (c). The quantum character is displayed by the clearly visible negativities of the regularized functions at the corresponding times. Note that the choice of does not affect the nonclassical properties of the state but leads to a rotation in phase space.
Finally we would like to note that the nonclassicality quasiprobabilities shown in Fig. 5 can be determined straightforwardly in experiments. This can be done by the method introduced in (103) and realized in (104); (105), which allows the direct measurement of the characteristic function of the function, cf. Eq. (42). This technique can be readily combined with the direct sampling approaches as developed for the nonclassicality quasiprobabilities of radiation fields (97); (98).
V Summary and Conclusions
In this work we considered the timedependent Hamiltonian of a nonlinear JaynesCummings system that is driven in quasiresonance. We showed, that timeordering effects have a crucial impact on the system and can therefore not be omitted. As the general solution of a timedependent Hamiltonian can become a cumbersome task, we introduce a method to circumvent this issue via quantizing the pump field. By extending the Hilbert space of the system, the dynamics becomes exactly solvable. Using the resulting timeindependent Hamiltonian we derived an analytical expression of the timeevolution operator.
For a pump field prepared in a coherent state the solutions were shown to converge to the classicalpump scenario where the discrepancies shrink with increasing coherent amplitude. Furthermore, we visualized the temporal evolution of the nonclassicality quasiprobability of the motional states of the ion. This regularized version of the often strongly singular GlauberSudarshan function has the advantage that it can be determined in experiments. Their negativities certify the quantum nature of the system under study. The introduced method to calculate this quasiprobability out of the density matrix applies to any time evolution.
In general, the derived algebra of the quasiresonantly driven trapped ion renders it possible to investigate complex scenarios where the interaction of the vibrational and the atomic (source) degrees of freedom is of interest. This may include the study of timedependent motional quantum correlation effects. Furthermore, our analytical approach may yield a deeper insight into the properties of nonequaltime commutation rules, in cases with explicitly timedependent interactions.
Acknowledgements.
We thank Ruynet Lima de Matos Filho for helpful comments. W. V. acknowledges funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 665148.References
 D. F. Walls, Squeezed states of light, Nature (London) 306, 141 (1983).
 R. E. Slusher, L. W. Hollberg, B. Yurke, J. C. Mertz, and J. F. Valley, Observation of Squeezed States Generated by FourWave Mixing in an Optical Cavity, Phys. Rev. Lett. 55, 2409 (1985).
 L.A. Wu, H. J. Kimble, J. L. Hall, and H. Wu, Generation of Squeezed States by Parametric Down Conversion, Phys. Rev. Lett. 57, 2520 (1986).
 H. Vahlbruch, M. Mehmet, S. Chelkowski, B. Hage, A. Franzen, N. Lastzka, S. Goßler, K. Danzmann, and R. Schnabel, Observation of Squeezed Light with 10db QuantumNoise Reduction, Phys. Rev. Lett. 100, 033602 (2008).
 H. Vahlbruch, M. Mehmet, K. Danzmann, and R. Schnabel, Detection of 15 dB Squeezed States of Light and their Application for the Absolute Calibration of Photoelectric Quantum Efficiency, Phys. Rev. Lett. 117, 110801 (2016).
 A. Einstein, N. Rosen, and B. Podolsky, Can QuantumMechanical Description of Physical reality Be Considered Complete?, Phys. Rev. 47, 777 (1935).
 E. Schrödinger, Die gegenwärtige Situation in der Quantenmechanik, Naturwiss. 23, 807 (1935).
 H. J. Kimble, M. Dagenais, and L. Mandel, Photon Antibunching in Resonance Fluorescence, Phys. Rev. Lett. 39, 691 (1977).
 W. Vogel, Squeezing and anomalous moments in resonance fluorescence, Phys. Rev. Lett. 67, 2450 (1991).
 B. Kühn, W. Vogel, M. Mraz, S. Köhnke, and B. Hage, Anomalous Quantum Correlations of Squeezed Light, Phys. Rev. Lett. 118, 153601 (2017).
 S. Gerber, D. Rotter, L. Slodicka, J. Eschner, H. J. Carmichael, and R. Blatt, IntensityField Correlation of SingleAtom Resonance Fluorescence, Phys. Rev. Lett. 102, 183601 (2009).
 L. Knöll, W. Vogel, and D.G. Welsch, Quantum noise in spectral filtering of light, J. Opt. Soc. Am. B 3, 1315 (1986).
 J. D. Cresser, Intensity correlations of frequencyfiltered light fields, J. Phys. B: At. Mol. Phys. 20, 4915 (1987).
 L. Knöll, W. Vogel, and D.G. Welsch, Action of passive, lossless optical systems in quantum optics, Phys. Rev. A 36, 3803 (1987).
 A. Stokes, Vacuum sourcefield correlations and advanced waves in quantum optics, Quantum 2, 46 (2018).
 E. T. Jaynes and F. W. Cummings, Comparison of Quantum and Semiclassical Radiation Theories with Application to the Beam Maser, Proc. IEEE 51, 89 (1963).
 H. Paul, Induzierte Emission bei starker Einstrahlung, Ann, Phys. (Leipzig) 11, 411 (1963).
 S. Haroche, Les Houches Lecture Notes, Session XXXVIII ed, New Trends in Atomic Physics, eds G. Grynberg and R. Stora (Elsevier, New York), 195 (1984).
 S. Haroche and J.M. Raimond, Radiative Properties of Rydberg States in Resonant Cavities, Adv. At. Mol. Phys. 20, 347 (1985).
 D. Meschede, H. Walther, and G. Müller, OneAtom Maser, Phys. Rev. Lett. 54, 551 (1985).
 G. Rempe, H. Walther and N. Klein, Observation of quantum collapse and revival in a oneatom maser, Phys. Rev. Lett. 58, 353 (1987).
 I. I. Rabi, On the Process of Space Quantization, Phys. Rev. 49, 324 (1936).
 I. I. Rabi, Space Quantization in a Gyrating Magnetic Field, Phys. Rev. 51, 652 (1937).
 D. Esteve, J.M. Raimond, and J. Dalibard, Entanglement and Information Processing, 1st ed., (Elsevier, Amsterdam, 2003).
 J. H. Eberly, N. B. Narozhny, and J. J. SanchezMondragon, Periodic Spontaneous Collapse and Revival in a Simple Quantum Model, Phys. Rev. Lett. 44, 1323 (1980).
 N. B. Narozhny, J. J. SanchezMondragon, and J. H. Eberly, Coherence versus incoherence: Collapse and revival in a simple quantum model, Phys. Rev. A 23, 236 (1981).
 J. R. KukliÅski and J. L. Madajczyk, Strong squeezing in the JaynesCummings model, Phys. Rev. A 37, 3175(R) (1988).
 M. Hillery, Squeezing and photon number in the JaynesCummings model, Phys. Rev. A 39, 1556 (1989).
 S. J. D. Phoenix and P. L. Knight, Establishment of an entangled atomfield state in the JaynesCummings model, Phys. Rev. A 44, 6023 (1991).
 S. Furuichi and M. Ohya, Entanglement Degree in the Time Development of the JaynesCummings Model, Lett. Math. Phys. 49, 279 (1999).
 S. Bose, I. FuentesGuridi, P. L. Knight, and V. Vedral, Subsystem Purity as an Enforcer of Entanglement, Phys. Rev. Lett. 87, 050401 (2001)
 H. J. Carmichael, Photon Antibunching and Squeezing for a Single Atom in a Resonant Cavity, Phys. Rev. Lett. 56, 539 (1986).
 M. Hennrich, A. Kuhn, and G. Rempe, Transition from Antibunching to Bunching in Cavity QED, Phys. Rev. Lett. 94, 053604 (2005).
 F. Dubin, D. Rotter, M. Mukherjee, C. Russo, J. Eschner, and R. Blatt, Photon Correlation versus Interference of SingleAtom Fluorescence in a HalfCavity, Phys. Rev. Lett. 98, 183003 (2007).
 M. Brune, S. Haroche, J. M. Raimond, L. Davidovich, and N. Zagury, Manipulation of photons in a cavity by dispersive atomfield coupling: Quantumnondemolition measurements and generation of ”Schrödinger cat” states, Phys. Rev. A 45, 5193 (1992).
 G.C. Guo and S.B. Zheng, Generation of Schrödinger cat states via the JaynesCummings model with large detuning, Phys. Lett. A 223, 332 (1996).
 J. J. Slosser, P. Meystre, and S. L. Braunstein, Harmonic oscillator driven by a quantum current, Phys. Rev. Lett. 63, 934 (1989).
 M. Weidinger, B. T. H. Varcoe, R. Heerlein, and H. Walther, Trapping States in the Micromaser, Phys. Rev. Lett. 82, 3795 (1999).
 S. Brattke, B. T. H. Varcoe, and H