Abstract
We theoretically investigate the quantum statistical properties of light transmitted through an atomic medium with strong optical nonlinearity induced by Rydberg–Rydberg van der Waals interactions. In our setup, atoms are located in a cavity and nonresonantly driven on a two-photon transition from their ground state to a Rydberg level via an intermediate state by the combination of the weak signal field and a strong control beam. To characterize the transmitted light, we compute the second-order correlation function . The simulations we obtained on the specific case of rubidium atoms suggest that the bunched or antibunched nature of the outgoing beam can be chosen at will by tuning the physical parameters appropriately.
Export citation and abstract BibTeX RIS
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
In an optically nonlinear atomic medium, the dispersion and absorption of a classical light beam depends on the powers of its amplitude [1]. At the quantum level, dispersive optical nonlinearities translate into effective interactions between photons. The ability to achieve such strong quantum optical nonlinearities is of prominent importance in quantum communication and computation, as it allows us to implement photonic conditional two-qubit gates. The standard Kerr dispersive non linearities obtained in noninteracting atomic ensembles, either in off-resonant two-level or resonant three-level configurations involving electromagnetically induced transparency (EIT), are too small to allow for quantum nonlinear optical manipulations. To further enhance such nonlinearities, EIT protocols were put forward in which the upper level of the ladder is a Rydberg level. In such schemes, the strong van der Waals interactions between Rydberg atoms result in a cooperative Rybderg blockade phenomenon [2–4], where each Rydberg atom prevents the excitation of its neighbours inside a 'blockade sphere'. This Rydberg blockade deeply changes the EIT profile and leads to magnified nonlinear susceptibilities [5–8]. In particular, giant dispersive nonlinear effects were experimentally obtained in an off-resonant Rydberg–EIT scheme using cold rubidium atoms placed in an optical cavity [9, 10]. In this paper, we theoretically investigate the quantum statistical properties of the light generated in the latter protocol. Note that, contrary to other theoretical works, e.g. [11, 12], here, we are interested in the dispersive region. Moreover, since we place the atoms in a cavity rather than in free-space, the theoretical framework and calculations we perform also differ from [11, 12]. In particular, a technical benefit of our approach is that we are not restricted to considering only photon pairs but could, in principle, investigate higher-order correlations.
We first write the dynamical equations for the system of interacting three-level atoms coupled to the strong control field and the nonresonant cavity mode fed by the probe beam. We show that, under some assumptions, the system effectively behaves as a large spin coupled to the cavity mode [13]. We then compute the steady-state second-order correlation function to characterize the emission of photons out of the cavity. Our numerical simulations suggest that the bunched or antibunched nature of the outgoing light, as well as its coherence time, may be controlled through adjusting the detuning between the cavity mode and probe field frequencies.
The paper is structured as follows. In section 2, we present our setup and the assumptions we make to compute its dynamics. We also explain the analytical and numerical methods we employ to calculate the second-order correlation function of the outgoing light beam. In section 3, we present and interpret the results of the simulations we obtained for and on the specific experimental case considered in [9]. Finally, we conclude in section 4 by evoking open questions and perspectives of our work. The appendices address the supplementary technical details that are omitted in the text for readability.
2. Model and methods
The system we consider comprises N atoms that present a three-level ladder structure with a ground , intermediate , and Rydberg states (see figure 1). The energy of the atomic level is denoted by ℏωk (by convention ), and the dipole decay rates from the intermediate and Rydberg states are denoted by γe and γr, respectively. The transitions and are, respectively, driven by a weak probe field of frequency ωp and a strong control field of frequency ωcf. To limit absorption, both fields are off-resonant; the respective detunings are given by and . Moreover, to enhance dispersive effects while keeping a high input–output coupling efficiency, the atoms are placed in an optical low-finesse cavity. The transition is supposed to be in the neighbourhood of a cavity resonance. The frequency and annihilation operator of the corresponding mode are denoted by ωc and a, respectively; the detuning of this mode with the probe laser is defined by , and α denotes the feeding rate of the cavity mode with the probe field, which is assumed real for simplicity. Finally, we introduce g and Ωcf, which are the single-atom coupling constant of the transition with the cavity mode and the Rabi frequency of the control field on the transition , respectively. In the following paragraphs, we study the dynamics of the system, which, under some assumptions, is equivalent to a damped harmonic oscillator, i.e. the cavity mode, coupled to an assembly of spins , and the Rydberg bubbles corresponding to the 'super-atoms' delimited by the Rydberg blockade spheres.
Starting from the full Hamiltonian, we perform the rotating wave approximation and adiabatically eliminate the intermediate state as described in appendix
coupled to the cavity mode of the effective decay rate
increased by the coupling to the atomic ensemble. The Hamiltonian reads
In this expression, we have introduced the atomic operators for as well as the effective detunings
and
respectively shifted from Δr and Δc by the light-induced Stark shift of the control beam and by the linear atomic susceptibility. The quantity is the van der Waals interaction between atoms in their Rydberg level—when atoms are in the ground or intermediate states, their interactions are neglected—while
is the effective coupling strength of the two-photon transition driven by the cavity mode and the control laser.
At this point, following [13], we introduce the Rydberg bubble approximation. In this approach, the strong Rydberg interactions are assumed to effectively split the sample into bubbles , each of which contains atoms but can only accomodate a single Rydberg excitation, delocalized over the bubble. Note that, within this approximation, all bubbles have the same radius, which is fixed by the coefficient C6 and the detunings Δr and Δe, as well as the control field Rabi frequency Ωcf; assuming that the atomic ensemble is homogeneous, the number of atoms per bubble nb is approximately given by [9]
where is the atomic density. Each bubble can, therefore, be viewed as an effective spin whose Hilbert space is spanned by
the ground state of the bubble and its symmetric singly Rydberg excited state, respectively. Introducing the bubble spin- operators , where the operator corresponds to the lowering operator of the spin and the annihilation of a Rydberg excitation, one can write the Hamiltonian under the approximate form (see appendix
where we introduced the collective angular momentum . The system is, therefore, equivalent to a large spin, i.e. the assembly of spin- Rydberg bubbles, coupled to a harmonic oscillator. Its density matrix satisfies the master equation
One can also write the Heisenberg–Langevin equations for the time-dependent operators
where are the Langevin forces associated with a and , respectively. Note that we neglected the effect of extra dephasing due to, for example, collisions or laser fluctuations.
To study the quantum properties of the light transmitted through the cavity, we shall compute the function , which characterizes the two-photon correlations. In the input–output formalism [14], one shows that this function simply equals the function for the intracavity field (see appendix
where ρss denotes the steady state of the system defined by ; see equation (1).
In the region of the small feeding parameter α, one can compute ρss numerically by propagating in time the initial state , where represents the symmetric state in which bubbles are excited, and are the Fock states of the cavity mode. To this end, one applies the Liouvillian evolution operator in a truncated basis, restricted to states of low numbers of excitations (typically with ). The steady state is reached in the limit of large times,ideally when . The denominator of the ratio in equation (2) is directly obtained from ρss. To compute its numerator, one first propagates in time from t = 0 to τ, using the same procedure as earlier, then applies the operator and takes the trace.
In the region of weak feeding, it is also possible to get a perturbative expression for by computing the expansion of and in powers of α. To this end, one uses the Heisenberg equations of the system in equations (2) and (3) to derive the hierarchy of equations relating the different mean values and correlations up to the fourth order in α. After straightforward, though lengthy, algebra, one gets an expression for that is too cumbersome to be reproduced here, but allows for faster calculations than the numerical approach. Such a fully analytical treatment, however, cannot, to our knowledge, be extended to ; for , therefore, we rely entirely on numerical simulations.
To conclude this section, we consider the region of a large number of bubbles and low number of excitations, i.e. and . As shown in appendix
Finally, we get the following approximate expression for the effective Hamiltonian
where . In this region, the system thus behaves as two coupled oscillators: one is harmonic, the cavity field; the other is anharmonic, the Rydberg bubble field. The cavity resonance is therefore shifted depending on the number of excitations in the system.
In the following section, we present and discuss the results we obtained with the specific system used in [9]. It appears that one can choose the bunched or antibunched behaviour of the light transmitted through the cavity by adjusting the detuning Δc. We also show that the time behaviour of the function depends on the region considered, and can be roughly understood as resulting from the damped exchange of a single excitation between the atoms and field.
3. Numerical results and discussion
We now consider the physical setup presented in [9], i.e. an ensemble of atoms whose state space is restricted to the levels , , and with the decay rates MHz and MHz. The other physical parameters must be designed so that strong nonlinearities may be observed at the single-photon level. In the specific system considered here, we find that this is achieved for a cavity decay rate MHz, a volume of the sample , a sample density , a control laser Rabi frequency , a cooperativity C = 1000, a detuning of the intermediate level , a detuning of the Rydberg level , and a cavity feeding rate . For these parameters, the cavity detuning corresponds to the maximal average number of photons in the cavity. Note that these physical parameters are experimentally realistic and feasible.
Let us first focus on the second-order correlation function at zero time , represented in figure 2(a) as a function of the reduced detuning . The numerical and analytical results are in such good agreement for the region considered that the corresponding curves cannot be distinguished. One notes a strong bunching peak (B) and a deep antibunching area centered on (A) . This suggests that around (A), photons are preferably emitted one by one, while around (B) they are preferably emitted by pairs. Note, however, that as a ratio, only gives information on the relative importance of pair and single-photon emissions. Its peaks, therefore, do not correspond to the maxima of photon pair emission, but to the best possible compromises between and , which can be checked by comparing figures 2(a) and (b). Hence, pair emission might dominate in a region where the number of photons coming out from the cavity is actually very small.
Download figure:
Standard image High-resolution imageWe now investigate the behaviour of for two different values of the detuning, i.e. and , which respectively correspond to the peak (B) and minimum (A) of . The numerical simulations we obtained are given in figure 3. The plot relative to (B) exhibits damped oscillations, alternatively showing a bunched or antibunched behaviour. The plot corresponding to (A) always remains on the antibunched side, though asymptotically tending to 1.
Download figure:
Standard image High-resolution imageThe dynamical features observed can be understood and satisfactorily accounted for by a simple three-level model. Indeed, due to the weakness of α, the system in its steady state is expected to contain, at most, two excitations (either photonic or atomic). After a photon detection at t = 0, it contains, at most, one excitation, which can be exchanged between the cavity field and atoms, as it has been known for a long time [15, 16]. In other words, the operator can be expanded in the space restricted to the three states and the effective non-Hermitian Hamiltonian for the system, in this subspace, takes the following form:
The order of magnitude of the frequencies and decay rates of the oscillations observed for in the specific cases (A,B) are satisfactorily recovered by this Hamiltonian, which validates the schematic model we used and suggests it comprises the main physical processes at work.
To conclude this section, it is worth mentioning that the two-boson approximation, though strictly speaking not applicable here—the parameters considered in this section indeed correspond to a number of bubbles —yields, however, the qualitative behaviour for . The minimum is correctly located, though slightly higher than in the spin model; the antibunching peak is slightly shifted towards positive detunings and is weaker than in the previous treatment. These discrepancies result from too low a value of the nonlinearity parameter ; they can be corrected by replacing with in the two-boson Hamiltonian. We first note that and coincide in the region of a large number of bubbles. Moreover, makes sense in the region of a low number of bubbles: in particular, when , i.e. when only one bubble is available, the nonlinearity proportional to diverges accordingly, thus forbidding the boson field to contain more than one excitation. Finally, let us mention that can also be recovered via a perturbative treatment of the full model, which will be presented in a future paper.
4. Conclusion
In this work, we studied how the strong Rydberg–Rydberg van der Waals interactions in an atomic medium may affect the quantum statistical properties of an incoming light beam. In our model, atoms are located in a low-finesse cavity and subject to a weak signal beam and a strong control field. These two fields nonresonantly drive the transition from the ground to a Rydberg level. The system was shown to effectively behave as a large spin coupled to a damped harmonic oscillator, i.e. the assembly of Rydberg bubbles and the cavity mode, respectively. The strong anharmonicity of the atomic spin affects the quantum statistics of the outgoing light beam. To demonstrate this effect, we performed analytical and numerical calculations of the second-order correlation function . The results we obtained on a specific physical example with rubidium atoms indeed show that the transmitted light presents either bunched or antibunched characters, depending on the detuning between the cavity mode and the probe field. This suggests that in such a setup, one could design light of arbitrary quantum statistics through appropriately adjusting the physical parameters.
In this work, we performed the Rydberg bubble approximation, which allowed us to derive a tractable effective Hamiltonian. This scheme is, however, questionable: interactions between bubbles are indeed neglected, and the different spatial arrangements of the bubbles in the sample are not considered. Though challenging, it would be interesting to run full simulations of the system, rejecting those states that are too far off-resonant due to Rydberg–Rydberg interactions. Besides validating the assumption of the present work, this would indeed enable us to consider other regions, such as, for instance, the case of resonant transition towards the Rydberg level. We also implicitly made the assumption that the cavity mode and control beam were homogeneous. Spatial variations should be included in the model, and their potential influence studied in future work. Finally, due to the very weak probe field region considered in this paper, we only presented results on the function : the production of correlated photons is indeed very unlikely. In principle, we can, however, numerically compute for any , which might be relevant in a future work, if addressing stronger probe fields.
Acknowledgments
This work was supported by the EU through the ERC Advanced Grant 246669 'DELPHI' and the Collaborative Project 600645 'SIQS'.
Appendix A.: Derivation of the effective Hamiltonian
A.1. Rotating wave approximation
The full Hamiltonian of the system can be written under the form
where , is the energy of the atomic level for (with the convention ), and denotes the van der Waals interaction between atoms in the Rydberg level. When atoms are in the ground or intermediate states, their interactions are neglected.
We switch to the rotating frame defined by where
and perform the rotating wave approximation to get the new Hamiltonian , where
with the detunings , , and .
The corresponding Heisenberg–Langevin equations are:
where ain and denote Langevin forces.
A.2. Elimination of the intermediate state
Let us now simplify the system. First, one deduces from equation (A.4) that σer is of second order in the small feeding constant α. The term can, therefore, be neglected in equation (A.3). Moreover, since the ground-state population remains dominant during the evolution of the system, we can write ; from equation (A.2), the steady-state solution for in the far detuned region is thus
Finally, substituting this relation into equations (A.1) and (A.3) one gets
where
are the parameters for the effective two-level model and are the modified Langevin noise operators
Note that, in the absence of collisional terms, one simply recovers the standard three-level EIT susceptibility in the far-detuned region
Finally, we get the effective Hamiltonian
A.3. Rybderg bubble approximation
As described in the main text, we introduce the Rydberg bubble approximation. In this approach, the strong Rydberg interactions are assumed to effectively split the sample into bubbles , each of which contains atoms but can only accomodate a single Rydberg excitation, delocalized over the bubble. Note that the number of atoms per bubble nb is approximately given by [9]
where is the atomic density. Each bubble can, therefore, be viewed as an effective spin whose Hilbert space is spanned by
the ground state of the bubble , and its symmetric singly Rydberg excited state, respectively. Introducing the bubble Pauli operators , the operator corresponds to the lowering operator of the spin and the annihilation of a Rydberg excitation, one can write
where we introduced the collective angular momentum . In the same way,
where we used . Finally, the Hamiltonian of the system takes the approximate form
which represents the interaction of the large spin with the cavity mode a.
A.4. Region of a large number of bubbles and a low number of excitations
From the well-known relation , we deduce the second-order operator equation
In the region of a large number of bubbles and for low excitation numbers, i.e. eigenstates of the total angular momentum with , the solution of this equation is approximately given by
whence, at the lowest order in the excitation number,
Injecting equation (A.7) into the previous form of the Hamiltonian, we get
Moreover, from equation (A.8), we deduce that the operator is approximately bosonic and, therefore, the Hamiltonian can finally be put under the form
where .
Appendix B.: Calculation of
By definition, the second-order correlation function for the outgoing field is
Using the relations [14]
and keeping only nonzero terms (all terms like and equal zero), one obtains in the numerator four nonzero terms
Let us consider the first term. Using the standard commutation relations between a and ain operators, we have:
Here we used the relation
where X is any system operator [14] and where is the Heaviside step-function (with ). Evaluating the other terms in the same way, one finally obtains
where and .