1. Introduction
In recent times, a lot of attention has been directed at the subject of excitons in bulk crystals due to the experimental observation of the so-called yellow exciton series in Cu
2O, up to a large principal quantum number of
[
1,
2,
3,
4]. In analogy to atomic physics, such excitons have been named Rydberg excitons (RE). By virtue of their specific properties, Rydberg excitons are of great interest in solid-state and optical physics. These objects, whose size scales as
and where the energy spacing of neighbouring states decreases as
, are well-suited for performing experiments in parameter ranges that are quite different from other quantum systems, such as atomic vapours. The small binding energy of RE makes them more sensitive to external magnetic and electric fields than other systems. The observation and detailed description of Rydberg excitons has opened a new field in condensed matter spectroscopy. Their specific properties have motivated both theoretical and experimental interest in this field, with studies ranging from their spectroscopy—that is, optical [
2], electro-optical [
3], and magneto-optical spectra [
4,
5] through non-atomic scaling laws [
6] to quantum chaos [
7].
In the last two years, particular interest has been directed to the study of the dynamical properties of systems with RE. Due to a large orbital radius, their dipole moments are exceptionally large for higher values of
n, so that the interaction between higher Rydberg excitons is exceptionally strong and leads to the so-called Rydberg blockade. The appearance of a Rydberg blockade separates linear and nonlinear regimes of dynamical phenomena with RE. In the nonlinear regime, the possibility of observing giant optical nonlinearities of RE in microcavities has been considered [
8], while the paper [
9] deals with strong interaction phenomena at very low densities of RE, which enables one to determine the contribution of the nonlinear optical response of the medium. In the case of a smaller exciton density, it is possible to remain in the linear range, for which a single-photon source has been proposed [
10]. RE mediums could be used as a gain system for solid-state, highly tunable masers, allowing to achieve output powers of order ranging from
to
W in continuous and pulse mode, respectively [
11]. The realization of light slowing or storing and retrieving experiments in crystal with RE could unlock a plethora of dynamical effects, which might be observed in such media. These phenomena are based on the electromagnetically induced transparency (EIT), which has various potential applications in both classical and quantum optics. Słowik et al. [
12] theoretically proposed to implement single-qubit gates and a two-qubit controlled not (CNOT) gate operating on polarized photons based on light storage. Recently, a highly efficient quantum memory protocol based on EIT has been experimentally demonstrated for cold rubidium [
13] and cesium atoms [
14]. Experimental evidence of spinor slow light as a quantum memory for two-color qubits in cold rubidium has also been presented [
15].
EIT [
16] is one of the most important effects utilized in quantum optics as it allows for the coherent control of materials’ optical properties. The generic EIT utilizes extraordinary dispersive properties of an atomic medium with three active states in the
configuration. This phenomenon leads to a significant reduction of absorption of a weak, resonant laser field (probe field) by irradiating the medium with a strong control field, and thus making an otherwise opaque medium transparent. This leads to dramatic changes of dispersion properties of the system. Absorption forms a dip called the transparency window and approaches zero, while the dispersion in the vicinity of this region becomes normal with a steep slope, which increases for a decreasing control field—the resonant probe beam is transmitted almost without losses. It should be stressed that there is a significant distinction between Electromagnetically Induced Transparency and Autler–Townes Splitting (ATS) [
17]. The Autler–Townes effect results in a “split” transition within a coupled three-level system due to the AC Stark effect and emerges in the strong coupling limit of the more recently discovered electromagnetically induced transparency [
17]. The EIT effect is described by the formation of a dark state due to a destructive quantum interference between the transition pathways. Both ATS and EIT result in a transparency feature, which is qualitatively identified as a wide spectral region between the split-absorption peaks for ATS, but a narrow transmission window within a single-absorption peak for EIT. This common feature has been at the centre of long-standing confusion as to whether an observed transparency is due to EIT or ATS, and as such, the distinction between the two is an active topic of research in the quantum optics community. It has recently been shown theoretically [
18], as well as experimentally [
19] that it is possible to objectively distinguish between the regime where EIT dominates (EIT regime) and the one where ATS dominates (ATS regime). Usually, EIT is explained by the destructive quantum interference between different excitation pathways of the excited state, or alternatively, in terms of a dark superposition of states. For at least 20 years, there has been a considerable level of activity devoted to the EIT [
20], which was motivated by the recognition of a number of its potential applications, such as slowing and storing the light (see e.g., [
21,
22,
23] and the References therein). Some of the most exciting developments on EIT refer to their manifestations in artificial atoms [
24,
25,
26,
27]. A remarkable quenching of absorption due to EIT in an undoped bulk crystal of Cu
2O in a
configuration involving lower levels was examined in [
28], and evidence of quantum coherences has been pointed out to exist in Cu
2O [
29].
In Rydberg systems, a ladder configuration enables one to couple long-living metastable, initially empty upper levels with the lower levels coupled by the probe field. Demonstration of EIT in Rydberg atoms involving the ladder levels scheme by Mohapatra et al. [
30] has taken advantage of their unique properties and entailed investigations concerning transient and steady-state EIT spectra in Rydberg systems [
31]. The Rydberg excitons offer an unprecedented potential to study the above-mentioned phenomena in solids—excitons in Cu
2O offer a great variety of accessible states for creating the ladder configuration, which enables such a choice of coupling that could be realized by accessible lasers, or which could eventually to be suitable for desirable coherent interaction implementations. The first step toward the study of the photon blockade, many-body physics with Rydberg excitons and quantum non-linear optics is the realization of Rydberg EIT. The optical nonlinearity can be significantly enhanced in the presence of quantum interference in EIT systems, so examination of this phenomena in RE may be the first step to proposing an eventual application of this medium in quantum information processing.
Here, we focus our interest on the dynamical aspects and properties of Rydberg excitons in Cu
2O, and we show that the Rydberg excitonic states can be used to realize the EIT. We indicate excitonic states which guarantee the most efficient realization of the experiment, taking into account our previous results concerning excitonic resonances [
2,
3,
4,
5], as well as damping parameters and the matrix elements of the interband dipole operator. The influence of nonlinear interactions will also be taken into account. Moreover, a study of Rydberg EIT spectra will be performed in the presence of dipolar interactions, and their dependence on excitonic density, principal quantum numbers of Rydberg states, and control beam intensities will be discussed.
Our paper is organized as follows. In
Section 2 we present the assumptions of the considered model and solve the time evolution equations, obtaining an analytical expression for the susceptibility and the group index. We use the obtained expression to compute the real and imaginary part of the susceptibility and the group index (
Section 3) for a Cu
2O crystal slab. We examine in detail the influence of the control field on the dispersion, absorption, and the group index of the medium. In
Section 4, we draw conclusions based on the model studied in this paper and indicate the optimal choice of Rydberg states to realize the EIT and light slowing in the linear regime.
2. Theory
The phenomenon of EIT, described qualitatively above, has been studied theoretically in various configurations of the transitions, probe, and control beams. Below, we propose a theoretical description for the case when the atomic transitions are replaced by intra-excitonic transitions in Rydberg excitons. The condensed matter exhibits quite a variety of three-level systems where an induced transparency could be achieved in much the same way as done in atomic mediums. Yet, dephasings, which can easily break the coherence of the population trapping state, are typically much faster in solids than in atomic vapours; it has caused great difficulty in observing a large, electromagnetically induced transparency effect in solids. We believe that this difficulty could be overcome by using the Rydberg exciton states. In particular, the higher states have exceptionally long lifetimes (which are proportional to
[
1]) and the dephasing can reach values which enables observation of the EIT effect.
Below, we consider a Cu
2O crystal as a medium where the EIT phenomenon can be realized. We use a ladder configuration (
Figure 1) consisting of three levels,
, and
c. As in previous works on Rydberg excitons, we focus our attention on the so-called yellow series associated with the lowest inter-band transition between the
valence band and the
conduction band. Because both band-edge states are of even parity, the lowest
exciton state is dipole-forbidden, whereas all the
P states are dipole-allowed; the
to
transition is also allowed. We have chosen the valence band as the
b state. As a practical example, the
and
excitonic states are chosen. The proposed ladder system contains only dipole-allowed
transitions.
Let the probe/signal field of frequency
and amplitude
couple the ground state
b of energy
with an excited state
a of energy
. The control field of frequency
and amplitude
couples the state
c of energy
with the state
a, as it is illustrated in
Figure 1. The Hamiltonian of such a three-level system interacting with an electromagnetic wave in the rotating wave approximation leads to
where h.c. stands for hermitian conjugate. The
,
,
,
are the wave vectors and real Rabi frequencies corresponding to the particular couplings, respectively. The
are the dipole transition moments related to the specific transitions, and
describes the nonlinear interactions between Rydberg excitons. One of the characteristic features of Rydberg media is the existence of the Rydberg blockade. Strong dipolar interaction between Rydberg excitons, which strongly depends on their distance, shifts the Rydberg levels, preventing the optical excitation of nearby excitons by shifting their corresponding levels out of the resonance with the exciting electromagnetic field. The term Rydberg blockade refers to the case where the interaction-induced shift is much larger than the EIT linewidth. In this case, the resulting nonlinearity can be interpreted as a switch from three-level EIT susceptibility to two-level susceptibility [
32]. Resonant absorption and exciton creation are no longer possible inside the blockade volume
, in which dipole interaction energy is larger than the EIT linewidth. This effect influences the generic, linear EIT, so that the nonlinear modification should be taken into account [
32]. In the presented manuscript, we consider the situation when one remains in the linear regime, that is, for exciton density small enough to avoid the Rydberg blockade. The presented calculations allow one to make such a choice of excitonic states where linear approximation still holds, such as the blockade-induced shift
; this is ensured by the exciton density which is two to three orders of magnitude below the limit where the excitons are packed closely together and
.
The time evolution of the system is governed by the von Neumann equation with a phenomenological relaxation contribution
where
denotes the density matrix for an exciton at position
z and time
t, and
R is the relaxation operator accounting for all relaxation processes in the medium.
By neglecting propagation effects for the control field and denoting the probe and control beam detunings by
and
, respectively, the evolution of the system can be described by the following equations:
At large separations between RE, the dipole–dipole interaction is dominant, and the potential
strongly depends on the separation distance between Rydberg excitons [
1]. The parameters
,
describe the damping of exciton states and are determined by temperature-dependent homogeneous broadening due to phonons, and broadening due to structural imperfections and eventual impurities. The relaxation damping rates for the coherence are denoted by
,
[
28]. It should be noted that in the above equations, only the relaxations inside the three level system are considered, so that the total probability for the populations of the three levels is conserved:
.
Due to the fact that at
,
for a weak probe field (i.e.,
, in the first-order perturbation with respect to the probe field, the evolution of our system, given by Bloch equations, reduces to a set of the following equations for the density matrix:
Taking into account the fact that for a weak probe pulse, for a given frequency, the polarization of the medium
is proportional to the signal field
and susceptibility
, it takes the form
where
N is the density of excitons,
is the transition dipole matrix element [
11], and
is the vacuum dielectric permittivity. In the limit of the low probe intensity, the susceptibility, which follows from the steady-state solutions of Equations (4) to (5), can be expanded into the Taylor series
where
is the linear, and
is the nonlinear part of electric susceptibility. The formula for the linear part
follows from the steady-state solution of Equations (
4) to (5) and
is treated as a correction, which depends on interaction potential
and causes the nonlinear modification of EIT due to the strong interactions between Rydberg excitons. For a large exciton separation
R, this potential has the form
, with
being a coefficient proportional to
, which drastically increases for higher Rydberg excitons. On the other hand, for closer separation of excitons, when the dipole–dipole interaction becomes comparable or even exceeds the energy spacing between excitonic levels, this interaction has the Foerster type and is proportional to
[
1]. In our system, we remain in the low-density regime where the
interaction is dominant.
The susceptibility is a complex, rapidly varying function of
, with its real part responsible for the dispersion and the imaginary part describing the absorption. If the excitons are driven on a single-photon resonance
, the main nonlinear effect will be the nonlinear absorption. Neglecting the transverse probe beam dynamics for lower excitonc states, that is, when van der Waals is dominant, it is possible to analytically calculate this nonlinear part of susceptibility [
33],
It should be stressed that the interaction between Rydberg excitons modifies the system in a nonlinear way by shifting the Rydberg levels, effectively changing the control field detuning and influencing the absorption and the width of the transparency window. Note that the equality of real and imaginary parts of
holds in a limited frequency range around the
resonance [
33].
Due to the dependence of the refractive index
on the medium dispersion, we define the group index
which accounts for the time delay of a pulse propagating in a medium with the group velocity
The slope of dispersion relation inside the transparency window determines the velocity of the pulse propagating inside the medium. Note that the derivative of the real part of the susceptibility may be positive (normal dispersion) or negative (anomalous dispersion); in the latter case, the group velocity may even become negative.
3. Numerical Results
We considered a Cu
2O crystal slab with a thickness of 30 μm as a medium where the EIT phenomenon can be realized. We used the ladder configuration (
Figure 1), consisting of three levels. Let the probe/signal field of frequency
couple the ground state (valence band) with the excited state
. The control field of a much lower frequency
couples states
and
. Using the Formulas (
7) and (
9), we calculate the real and imaginary part of the susceptibility and the group index. For the state energies and damping parameters, we use the values obtained from fitting to experimental data, as described in [
3]. The transition dipole moments were calculated from hydrogen-like functions, as described in detail in [
11]. In particular, we used Hydrogen P-orbital with the radius proportional to
[
1], which results in dipole moments in the range of
ea
(1S state) up to over 150 ea
(10P state) [
11]. Note that instead of
, we used the function
, which is an analogue of the dipole moment element, but takes into account the smeared-out transition dipole density characteristic for coupling between the valence band and the 2
P excitonic state [
4]. Application of this smeared-out dipole density is justified by a good agreement with an experiment of our previous results for electro-optical properties of Rydberg excitons [
3]. We chose the exciton density
, which is two orders of magnitude below the maximum value of
cm
−3 for
[
1].
The system parameters are summarized in
Table 1. We assumed that
[
28]; the damping rates are calculated for
K, where the contribution of exciton scattering with acoustic phonons is relatively small, and the main factor is spontaneous, with an emission rate proportional to
[
11].
The choice of states
and
has a fundamental effect on the properties of the system. The interaction coefficient
shows the strongest dependence on the state number and puts an upper limit on the lower state
. Moreover, the maximum exciton density is inversely proportional to the Rydberg blockade volume
(note that the upper state
is empty, so that the blockade does not apply to it). The state
affects the susceptibility by relaxation rates. All these factors indicate that a successful light slowdown can be achieved only for some parameters and the state configurations, where the linear part of susceptibility is dominating, and thus the group velocity is well-defined. Therefore, as a first step, we have investigated the range of states and control field intensities where the linear part of susceptibility is dominating, as shown in the
Figure 2. As pointed out above, the nonlinearity is very sensitive to the state number and varies by many orders of magnitude. One can also see that, in general, linear response is obtained for weaker control fields, which is advantageous for a large slowdown factor, but results in a relatively narrow transparency window. In the range where
is negligible, there is very little absorption inside the transparency window, and the group velocity
is a well-defined quantity that can be calculated in a linear regime from Equation (11).
The linear region is outlined on the insets of the
Figure 2. One can see that, practically, only the lower states
= 1, … 4 are usable, except for extremely small values of control field intensity which are impractical due to the fact that in EIT, the control field needs to be much stronger than the probe field, that is,
. On the other hand, the
should be large; due to the fact that this level remains empty, it is not constrained by the Rydberg blockade.
The
Figure 3 shows the medium susceptibility and the group index
. For a non-zero Rabi frequency
of the control field, the imaginary part of the system’s susceptibility reveals a dip in the Lorentzian absorption profile, called a transparency window. This means that the resonant probe beam which otherwise would be strongly absorbed is now transmitted almost without losses. The width of the transparency window is proportional to the square of the control field amplitude, and therefore, by increasing the control field strength, it is possible to open it out. This is shown in the
Figure 4, where one can observe a dip in absorption (Im
) around the zero-detuning part of the spectrum, which is accompanied by a very steep slope of the Re
. That strong, normal dispersion inside the window is responsible for the reduction of the group velocity. The absorption at the resonance does not reach zero due to the finite value of the relaxation rate
, but it is indeed very small (absorption coefficient
3 cm
−1). This is easily visible in the
Figure 3, where the control field is set to
GHz. This means that the medium has become almost transparent for a probe pulse which now travels with a reduced velocity. Indeed, the group index shows that the impulse slowdown factor of about 3700 is expected. It should be stressed that the value of
is directly proportional to the magnitude of susceptibility changes, and thus, to the exciton density. Clearly, high values are favourable; for the lowest state, the densities on the order of
cm
−3 are possible [
34,
35]. In our calculations, we have chosen
cm
−3, which corresponds to susceptibility values on the order of
(
Figure 3), comparable with values presented in [
28] for the low excitonic states and well below the limit imposed by the Rydberg blockade for the chosen state. It should be stressed that while these densities are orders of magnitude greater than in atomic EIT systems, the possible impulse slowdown is limited by much larger dissipation rates, particularly
.
As mentioned previously, there are many Rydberg state combinations that provide a working EIT system. The maximum possible slowdown in all sets of states
,
, a strong control field
GHz, and probe field
GHz is shown in the
Figure 5a. The blue dots mark the systems where the linear part of susceptibility is much greater than the nonlinear, so that the group velocity is a well-defined quantity. The red dots mark the range where
, so that the calculated slowdown is only a rough estimation. Remarkably, a very large slowdown can be achieved for
due to the small size of excitons and, correspondingly, large exciton density. In the strong control field regime, the linear systems can be based only on
; for higher states, the effect of a Rydberg blockade becomes significant. Moreover, it is beneficial to use the high upper state
so that the probe and control field frequencies are similar, for example,
. It should be pointed out that a high value of
is only possible at highly cryogenic temperatures. Additionally, the probe beam may excite a group of closely-lying upper states, further degrading the performance [
11].
By reducing the probe and control field strengths to
GHz and
GHz, correspondingly, one can extend the range of linear systems to
, as shown in the
Figure 5b. In such a case, the transparency window is much narrower. Due to this, the slowdown in the systems with
is limited by the large dissipation constant
, which prevents the transparency window from fully opening. On the other hand, the weaker probe and control fields result in a smaller density, so that the Rydberg blockade effect is delayed to
.