Abstract
Photon echoes in rare-earth-doped crystals are studied to understand the challenges of making broadband quantum memories using the atomic frequency comb (AFC) protocol in systems with hyperfine structure. The hyperfine structure of Pr3+ poses an obstacle to this goal because frequencies associated with the hyperfine transitions change the simple picture of modulation at an externally imposed frequency. The current work focuses on the intermediate case where the hyperfine spacing is comparable to the comb spacing, a challenging regime that has recently been considered. Operating in this regime may facilitate storing quantum information over a larger spectral range in such systems. In this work, we prepare broadband AFCs using optical combs with tooth spacings ranging from 1 MHz to 16 MHz in fine steps, and measure transmission spectra and photon echoes for each. We predict the spectra and echoes theoretically using the optical combs as input to either a rate equation code or a density matrix code, which calculates the redistribution of populations. We then use the redistributed populations as input to a semiclassical theory using the frequency-dependent dielectric function. The two sets of predictions each give a good, but different account of the photon echoes.
© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Quantum memory is a crucially important technology. Perhaps most prominently, it will play an important role in long-range quantum networks [1,2], which will be required to fully realize the power of quantum computers and other quantum technologies: quantum repeaters will be necessary to overcome transmission losses in these networks [3], and most current approaches to quantum repeaters depend on optical quantum memory [4,5], a device for storing quantum states and recovering them at a later time. A practical quantum memory would enable long-lived storage and high-fidelity, high-efficiency retrieval of broadband single photons. Rare-earth-ion-doped crystals are promising candidates for implementing quantum memory due to their long ground-state lifetimes, large inhomogeneous bandwidths, and narrow homogeneous bandwidths [6]. Various approaches to quantum memory have been developed for rare-earth-crystal platforms. The atomic frequency comb (AFC) protocol is one such scheme [7]. In the simplest case, the central idea is to store single photons in a material with a frequency-comb absorption spectrum, with comb-tooth spacing $f_{\rm rep}$, and then to recover photon echoes after a time $T_{\rm rep}=1/f_{\rm rep}$. On-demand storage, i.e., with selectable storage time, is also possible using control signals.
Pr$^{3+}$:Y$_2$SiO$_5$ (Pr:YSO) is suitable for various optical storage techniques, and it was the first material used to demonstrate on-demand AFC quantum memory [8]. It has been studied extensively, including work on photon echoes [9,10], spectral-hole burning [11], optical filtering [12], electromagnetically induced transparency [13], and stimulated Raman adiabatic passage [14,15]. The latter technique was used for selective retrieval of pulses [16], and further studies investigated the rephasing of these coherent populations [17] and on-demand retrieval based on spontaneous Raman scattering [18]. Storage of optical pulses has been demonstrated with stopped light using electromagnetically induced transparency [19,20] and a memory has been realized at the single-photon level using stopped light in a spectral hole [21]. Gradient-echo quantum memory was demonstrated [22], and AFC quantum memory has also been demonstrated at the single-photon level [23,24]. Recently, Pr:YSO has been used for high-rate entanglement distribution as a step towards the implementation of quantum repeaters [25] and on-demand quantum memory [26].
The AFC protocol has also been studied in other rare-earth-ion-doped crystals, which have similarly advantageous optical properties. Reference [27] studied the effect of comb-tooth shape on echo efficiency. Storage and retrieval of single photons has been demonstrated: time-bin qubits [28,29] and entanglement with another photon [30] storage. Temporal multimode storage has been studied [31]. Pulse storage and retrieval has also been demonstrated using Stark shifts [32], superhyperfine levels [33] and hybridized electron-nuclear hyperfine levels [34]. Gigahertz bandwidths were demonstrated using Tm$^{3+}$ [30] and Er$^{3+}$ [35] ions. Er$^{3+}$ has also been studied for use in quantum repeaters [36].
One attractive feature of the AFC protocol highlighted in Ref. [7] is that in principle it enables retrieval efficiency arbitrarily close to unity. In practice, however, experimental demonstrations have achieved limited efficiency: Refs. [37,38] used cavity enhancement to reach efficiencies just over 50%.
Many rare-earth species and isotopes, including Pr:YSO, have hyperfine structure, which presents an obstacle to preparing high-bandwidth AFCs. To circumvent this challenge, past works on AFCs in Pr:YSO have either (1) used an AFC bandwidth smaller than the hyperfine spacing [18,21,23], (2) used an AFC tooth spacing larger than the entire 37 MHz spectral width of the hyperfine levels [39], or (3) both [24]. Here, we explore the intermediate regime in which the AFC bandwidth is larger than 37 MHz and the AFC tooth spacings are comparable to the hyperfine spacing. While some past studies on AFCs worked in this regime [33,40], we systematically investigate AFC formation and echo retrieval for a range of comb-tooth spacings. We experimentally demonstrate AFC formation and echo retrieval across a range of comb-tooth spacings. Through this systematic study, we hope to elucidate the relationship between AFC tooth spacing and hyperfine structure, which could aid in the design of future high efficiency AFC-based quantum memories.
We also develop a semiclassical theory of photon echoes which calculates the effect of the experimental optical fields applied to the crystal. Historically, photon echo phenomena have been understood as collective emission by atoms prepared in a Dicke state [41], whereas here we achieve good results with a semiclassical model in which individual Pr$^{3+}$ ions, treated quantum mechanically, interact with a classical electromagnetic field. Although several authors have noted the presence of a semiclassical limit [7,42,43], and the formalism of this limit was presented in a recent review article [44], there does not appear to be a realistic calculation of photon echoes for rare-earth-doped crystals in the literature, particularly not one which calculates the AFC based on an observed optical signal as well as the photon echoes. We calculate the atomic population shifts that occur during AFC preparation; the change to the dielectric response due to these shifts; and the echoes that arise from a pulse interacting with the highly dispersive medium. We present two sets of results, one corresponding to AFCs as predicted by the rate equations and a second set as predicted by the density matrix formalism. Given the redistribution of population which is a quantitative description of the AFC, the photon echo is described within the same formalism of highly dispersive linear optical response.
Regarding other calculations in this area, Thiel et al. [45] gave a theoretical account matching experimental results of pulse echo linewidths in three rare-earth crystals Tm$^{3+}$:YAG, Tm$^{3+}$:LiNbO3 and Tm3$^{3+}$:YGG. More frequently, the calculations are models [46] or done to propose experiments [47,48] rather than to do a careful comparison.
The YSO crystal structure has two distinct Y sites where Pr$^{3+}$ impurities can occur. Here we study Pr$^{3+}$ ions at site 1 [49,50], which have $^3H_4\,\rightarrow \,^1D_2~(\mathbf {g}\rightarrow \mathbf {e})$ optical transitions at 605.977 nm. The homogeneous linewidth of these optical transitions is less than 1 kHz and can be observed when the crystal is cooled to liquid helium temperatures to eliminate phonon broadening; the inhomogeneous broadening, arising from local variations in the crystal structure surrounding each Pr$^{3+}$ ion, is typically several gigahertz [51]. In the absence of an external magnetic field, the hyperfine states are all two-fold degenerate, and the ground- and excited-state manifolds each have three distinct energy levels (Fig. 1). Therefore, each ion has nine transition energies labeled $a-i$ in ascending order; conversely, for any laser frequency within the inhomogeneous bandwidth, there are nine different classes of ions $A-I$ for which one of these transitions is resonant with the laser, with the laser addressing transition $x$ for class $X$.
As a result of the hyperfine structure, spectral-hole burning at even a single frequency results in multiple features (Fig. 2). The large transmittance peak at zero detuning occurs because all 9 classes become transparent at this frequency. The smaller peaks occur because some subset of the classes become transparent: for example, only 3 of the 9 classes that would normally absorb at 9.42 MHz detuning are depopulated. Due to power broadening, our measurement cannot distinguish features separated by less than 0.2 MHz, but such closely-spaced peaks may appear as broader peaks or shoulders. The transmission peaks in Fig. 2 correspond to population depletion at detuning frequencies shown in bold in Fig. 1(c). The other frequencies correspond to population enhancement, hence transmission dips in most cases.
The excited-state spacings in Pr:YSO are nearly equal, leading to a comb-like response to spectral-hole burning at a single frequency. This occurs because (1) the transparency for classes $A,D,G$ at 4.58 MHz detuning overlaps with the transparency for $B,E,H$ at 4.84 MHz, and (2) an additional transparency occurs for $A,D,G$ at 9.42 MHz, further extending the comb-like response (similar reasoning for classes $C,F,I$ and $B,E,H$ explains the peaks at $-9.42$ MHz, $-4.84$ MHz, and $-4.58$ MHz). This observation suggests that the material has a naturally periodic spectrum. Although it has been noted in previous works, for example Ref. [33], that this periodicity dictates which tooth spacings are compatible with AFC formation, here we carefully examine tooth spacings close to the periodicity. In fact, we find that AFC formation is optimal not for tooth spacings exactly matching the excited-state hyperfine spacings, but rather a few hundred kilohertz detuned from those values.
2. Experiment
We use our experimental setup for two different measurement protocols, a comb-measurement protocol and a pulse-echo protocol, each of which proceeds in two stages as shown in Fig. 4. In the first stage of both protocols, we prepare an AFC. In the second stage, we either probe the transmission spectrum (the comb-measurement protocol), or observe the transmission and echoes of a pulse sent to the AFC (the pulse-echo protocol).
We use a tunable diode laser near 1212 nm, which is amplified in a fiber Raman system and then doubled, producing a beam at the $^3H_4\rightarrow {^1}\!D_2$ optical transition. To stabilize the laser frequency, we use Pound-Drever-Hall feedback [52] to lock the infrared light to a reference cavity with a Zerodur (Schott Glass), spacer, held in a vacuum chamber. This beam is modulated using double-passed acousto-optic modulators (AOMs) to prepare the three different beams used in our measurements. A burn beam is produced by applying a frequency-modulated (FM) radio-frequency (RF) sine wave to the AOM. This signal has a frequency-comb structure, centered at 100 MHz, with 60 MHz bandwidth, and a tooth spacing equal to the FM modulation frequency. After the double pass, this results in an optical frequency comb with 120 MHz bandwidth and tooth spacing equal to the FM modulation frequency (Fig. 3). A probe beam is produced by (slowly) sweeping the AOM frequency, resulting in a frequency-swept beam covering the same bandwidth as the burn beam. Alternatively, 10 ns pulses are generated by applying RF pulses to the AOM.
In our setup, we use a Pr:YSO crystal with 0.05% substitution of Pr for Y, housed in a cryostat maintained at 3.5 K. Two beams are focused to overlapping waists inside the crystal. Beam 1 is the burn beam described above: we vary the tooth spacing of this optical frequency comb in 0.1 MHz steps, from 1.0 MHz to 16 MHz. The first stage of both measurement protocols consists of 30 seconds of nearly continuous burning with Beam 1, aiming to imprint the optical frequency comb on the atomic transmission spectrum and establish an AFC. In our comb-measurement protocol, Beam 2 is the probe beam described above. During the 4 ms frequency sweep of the probe beam, we acquire the transmission spectrum of the crystal. In the pulse-echo protocol, Beam 2 is the 10 ns pulse described above. We monitor emission from the crystal to detect the transmitted pulse and its echoes. Due to the different natures of the output signals from the two measurements, detectors with different response bandwidths are required. In the second stage of each measurement protocol, the 1 s measurement cycle is synchronized to an electrical pulse indicating the cryostat pump cycle, in order to minimize effects of vibrations. We repeat the measurements for approximately 30 cycles [53] and average the results. In addition to the 30 s burn in stage 1, in stage 2 of each cycle, 930 ms is spent burning the crystal to reinforce the AFC.
3. Theory
Most discussions of the AFC protocol rely on preparation of a Dicke state of the atomic ensemble [7,8,31,33,35,39,40,42,43,54–58]; in this paper, however, we show that the AFC transmission spectra and echo pulses can be obtained using only semiclassical concepts, provided the beams are in coherent states. Not all authors invoke the Dicke state. Chanelière et al. [59] also outlines a semiclassical treatment of photon echoes in AFCs, but the calculation of the susceptibility is not explicit. Burman and Le Gouet [60] have made an analysis of both uniform and random AFCs in highly dispersive media using the Maxwell-Bloch equations for a two-level system, with an emphasis on analytic modelling.
In Section 3.1, we describe the theory required for calculating the dielectric response function. In Section 3.2 we show that this function can be used to calculate the propagation of light pulses through the crystal. Finally, in Section 3.3 we describe our method for determining how an incident field changes the ground-state population distribution.
3.1 Calculating the dielectric response function
The Clausius-Mossotti relation gives the macroscopic dielectric function $\varepsilon (\omega )$ of a solid in terms of the polarizabilities $\alpha _k(\omega )$ of its constituent atoms. The relation remains valid for particles embedded in a homogeneous medium if $\varepsilon (\omega )$ is interpreted as the dielectric constant relative to the undoped crystal dielectric constant $\epsilon _0\varepsilon _{\rm bkg}$ [61]:
The validity of the Clausius-Mossotti equation depends on our ability to define regions within the solid where charge does not pass. If such units do not exist [65], then intrinsically solid-state approaches such as band theory [66] or the modern theory of polarizability [67] are required to calculate the dielectric function. In the present case, the Pr$^{3+}$ impurity ions (a) have open 4f shells which are located in the interior of the ions, (b) are dilute, hence located far from each other on the scale of an atomic bond length, and (c) individually have very narrow resonances, whose energies differ from ion to ion by many times the resonance width. Each of these circumstances inhibits any conventional quantum mechanical interactions such as overlap of wave functions leading to bonding-antibonding splittings or band structure. Hence, Martin’s criterion [65] of a well-defined surface to separate the various oscillating charges is satisfied, building our confidence in the applicability of the Clausius-Mossotti equation. The surface in question is simply a sphere drawn around each individual Pr$^{3+}$ ion well outside of the 4f electrons.
The sum over $k$ in Eq. (1) may be written as a sum over the ground-state hyperfine levels $i$, integrated over the detuning $\Delta _0$ of their lowest-energy hyperfine transition $a$ (Fig. 1). We measure $\Delta _0$ relative to the center of the inhomogeneous band $\omega _0$. We define
In this work, $i$ runs over the three ground states $|\,\mathbf {g},{1}/{2}\rangle$, $|\,\mathbf {g},{3}/{2}\rangle$, and $|\,\mathbf {g},{5}/{2}\rangle$, and we consider transitions to the three excited states $|\,\mathbf {e}, {1}/{2}\rangle$, $|\,\mathbf {e}, {3}/{2}\rangle$, and $|\,\mathbf {e}, {5}/{2}\rangle$. However, our formalism applies to arbitrary numbers of ground and excited states: this makes it readily adaptable to cases where the degenerate hyperfine levels are split due to external fields, and even to calculations for other dopant ions.
In the rotating-wave approximation (RWA), the atomic polarizability has the Drude-Lorentz form
Reference [11] gives oscillator strengths for the 9 hyperfine transitions in Pr:YSO. Within a manifold, the transition rates are proportional to these oscillator strengths. Oscillator strengths differ from squares of matrix elements by a factor proportional to the transition frequencies. This factor is $1+O(10^{-7})$ because the hyperfine splittings are less than 50 MHz and the reference transition frequency is 495 THz.
Equation (4) allows us to write
Using the convention $\omega =\omega _0+\Delta$, the polarizability can be written as
Similarly, we define $\tilde \varepsilon (\Delta )=\varepsilon (\omega )$ and $\tilde B_i(\Delta ) = B_i(\omega )$. We may rewrite Eq. (2) as
To find the dielectric response of the unperturbed impurity band in our study, the three ground-state hyperfine levels were each assumed to have an identical Gaussian density of states with a full width at half maximum (FWHM) of 4400 MHz (standard deviation $\sigma _{\rm inh}= 1869$ MHz). When the crystal is unperturbed, all three states have equal occupancy. The frequency step is 50 kHz with $2^{20}$, i.e., over 1 million, points used in the FFT, leading to a bandwidth of 52 GHz, which is more than $28\,\sigma _{\rm inh}$. The polarization sums were formed on this grid leading to a function proportional to the polarizability, or equivalently $\varepsilon (\omega )-1$. The proportionality constant was set at the peak absorption as described in Section 3.2.
3.2 Pulse propagation
Taking the beam direction in the crystal to be the $+z$ direction, a scalar wave propagating in such a medium obeys
where $U_0$ is the incident optical amplitude, $k=2\pi n_{\rm bkg} / \lambda$ is the wavevector in an undoped YSO crystal, $n=\sqrt {\varepsilon }=n_1 + \text {i} n_2$ is the adjustment to the refractive index $n_{\rm bkg}$ due to the dopant ions, and $\lambda$ is the free-space wavelength. The transmission is given byAssuming $T=9.9$ dB, as observed in our single-frequency-burn measurements in the 10 mm long crystal, we evaluate Eq. (10) at the center of the inhomogeneous band, finding $n_2(\omega _0)=1.1 \times 10^{-5}$. Here $n_1(\omega _0)=1$, leading to $\varepsilon (\omega _0)-1 = 2 \, \text {i} \, n_2(\omega _0)$ where we neglect the term $[n_2(\omega _0)]^2$.
In practice, we find the sum of convolution integrals $\sum _i\tilde B_i(\Delta )$, divide by $\sum _i\tilde B_i(0)$ for the unperturbed system and multiply by $\tilde \varepsilon (0)-1$, found experimentally, to obtain $\tilde \varepsilon (\Delta )-1$. Thus, we avoid an explicit value for the density of impurity ions.
Our model of photon echoes consists of applying Eq. (9) to each Fourier component of a pulse incident on the crystal with a waveform given by $U_0(t)$. Because the pulse has a narrow bandwidth about a central frequency, the transform variable is the detuning $\Delta$. Because the crystal has an antireflective (AR) coating, reflections are neglected. The dielectric response is given by $\tilde \varepsilon (\Delta )$. At the exit face, the frequency components are back transformed to the time domain yielding $U_1(t)$ at the exit face. Symbolically,
3.3 Calculations of ground-state population redistribution
We calculate the AFC spectra, i.e., the redistribution of the ground-state population using two different methods, namely by solving the rate equations and by using a density matrix formalism. Each approach has strengths and weaknesses. The density matrix formalism is more generally applicable, while the rate equations may be derived from the density matrix under certain physical conditions. When the rate equations are valid, only the populations, and not the coherences, need be calculated. In both cases, we treat each ion class independently: a more thorough analysis could consider the interplay between different classes, to include spectral diffusion effects, for example, but this is beyond the scope of our present work. In addition, we neglect relaxation of the ground-state which occurs on a time scale of 100 s vs. 164 $\mu$s for relaxation of the excited state [11].
The rate equations, introduced by Einstein in 1905, have played a huge role in understanding optical systems such as lasers that depend on redistribution of population of quantum levels. The density matrix is a more sophisticated level of theory which also includes quantum coherence. While the rate equation captures the effects of absorption, stimulated emission, and spontaneous emission, it cannot capture coherent phenomena such as Rabi oscillations and coherent population redistribution mechanisms such as Stimulated Raman Adiabatic Passage (STIRAP) [69]. Indeed, the name “counterintuitive pulse sequence” associated with STIRAP is “counterintuitive” under the unstated assumption that the intuition was built on the rate equations.
The text of Grynberg et al. [70] discusses the validity conditions for the rate equations. Neither of the two conditions they discuss apply: neither do the coherences relax substantially faster than the populations are transferred, nor is the correlation time of the laser small compared to the decay time associated with the (power-broadened) absorption line. Despite a question of the region of validity, we present results of atomic frequency combs based on the rate equations. We also present a calculation of the redistribution of population using the density matrix.
We considered using the stochastic Schrödinger equation (SSE) [71] as well. An average of the SSE is formally equivalent to the density matrix [72], although the SSE has advantages for single-particle dynamics [71]. Reference [72] suggests that the density matrix is favored algorithmically for small matrices, and the SSE is favored for large ones. Reference [72] describes a cross-over when the dimension of the Hamiltonian exceeds 200. Here, the Hamiltonian is only of dimension 6, suggesting the density matrix formalism will be more efficient. This is true, but it understates the case. Because our Hamiltonian models a laser field which is periodic in time, it is possible to obtain a solution in the logarithm of the number of periods. Attal et al. [73] make a similar observation to accelerate a similar calculation using the rate equation. We were not able to find a comparable algorithm for the SSE, making the density matrix the better method for this application.
3.3.1 Rate equation formalism
In our rate equation calculations, we assume the entire population is initially distributed equally among the 3 ground-state levels, and describe their occupancies with a vector
Note that the off-diagonal elements are necessarily non-negative, since ions excited from $|\,\mathbf {g},j\rangle$ can only increase the population of $|\,\mathbf {g},\ell \neq j\rangle$; since these ions do not change the total population of all the ground states, the diagonal elements $R_{jj}\leq 0$. The time evolution of the populations is calculated by
3.3.2 Density matrix formulation
The Lindblad master equation for the density matrix $\rho$ is
The further development of this equation requires an explicit form for the Hamiltonian. Unfortunately, published Hamiltonian parameters are in conflict and are internally inconsistent. Parameters for the Hamiltonian determined experimentally by Longdell et al. [74] were later found to be in conflict with the empirical optical oscillator strengths of Nilsson et al. [11]. In an attempt to resolve this conflict, Lovrić et al. [75] remeasured the Hamiltonian parameters. Unfortunately, these new parameters result in an oscillator-strength matrix with two columns swapped when compared with Nilsson et al. and even with their own work, as has been pointed out previously [76].
Since we were unable to find the Hamiltonian elsewhere in the literature, we derived a Hamiltonian that is consistent with the hyperfine level structure and oscillator strengths [75] and that takes the conventional form $\frac {1}{2} I\cdot Q \cdot I$, where $I$ is a vector of length 3 representing the spin $\frac {5}{2}$ nuclear spin states and $Q$ is a rank-2 traceless symmetric tensor giving the orientation of these spins in a given manifold relative to the crystal axes. There are 10 degrees of freedom to be determined: the ground- and excited-state manifolds each have 2 energy spacings (Fig. 1) and 3 Euler angles indicating their orientation with respect to the crystal axes. However, it is actually sufficient to know only the relative orientation of the manifolds, so we only need to determine 3 Euler angles. The 3$\times$3 oscillator-strength matrix given in Ref. [51] provides four constraints on these Euler angles. Using a least-squares fit we find Euler angles $\alpha _E=10.3(15)^\circ$, $\beta _E = -164.4(15)^\circ$, and $\gamma _E=-130.7(15)^\circ$ (digits in parentheses, here and throughout, are total uncertainties at one standard error, and represent the uncertainty of the least significant digits). The RMS deviation of the squares of the matrix elements in Table 1 from the oscillator strength values given in Ref. [51] is 0.013. This value may be compared with the quadrature sum of the experimental uncertainty of 0.01 in Ref. [51] and the RMS uncertainty from the above fitting procedure, leading to a combined uncertainty of 0.023.
Returning to the Lindblad Master equation, we proceed to a rotating frame
In practice, the matrix exponential is found through its Taylor series in the form used in Horner’s method, namely
for any matrix $M$. The terms are evaluated from the innermost parenthesis, starting from $\frac {M}{30}$. This method proves to be more accurate than evaluating the eigenvalues and eigenvectors with the linear algebra package LAPACK [77] and taking the exponential, stabilizing the solutions.To evaluate the solution for a large number of frequency-modulated periods, we make use of the identity
Finally, the Rabi overall coupling strength $\Omega _0$ is a fitting parameter. We compare theory to experiment for the single-frequency burn given in Fig. 2. The linewidth is a function of $\Omega _0$, and the value that matches the FWHM of the central peak is shown in the plot. The five prominent transmission features are in good agreement with the experiment as well as the antiholes at positive detuning. The more complex pattern at negative detuning is in somewhat less satisfactory agreement with the experiment. For the AFCs, we present values for $\Omega _0=0.0125$. We also calculated for 8 other values between $\Omega _0=0.0032$ and $0.02$. The results are quite similar in the range of 0.005 to 0.02, although the peak coefficient of correlation with the experimental results below is maximized with the value presented.
The techniques developed in this section can be applied to periodic optical excitations of the crystal. In practice, the burn beam produced in the experiment is not exactly periodic, and the procedure for obtaining a periodic approximation to it is fairly subtle (Fig. 5). Details about this signal averaging step are in the Appendix.
In principle, spatial effects play a role in spectral-hole burning. For example, along the burn beam propagation direction, transmission peaks will become sharper as the beam is progressively filtered by the crystal. We leave these refinements to future work.
To simulate our experimental AFC preparation, we run the density matrix code at 7001 frequencies over a bandwidth of 350 MHz (in 50 kHz steps which are smaller than the narrowest features experimentally observed when applying a single-frequency burn to the crystal, as shown in Fig. 2), yielding predicted ground-state population fractions for each of three hyperfine levels. Outside of the 350 MHz band, the populations are assumed to be undisturbed. We multiply all population fractions by the Gaussian density of states given in Section 3.1 to obtain the absorbing populations throughout the inhomogeneous band. The highly dispersive dielectric function of the AFC was calculated from these populations using Eq. (8) and the equations leading to it in Section 3.1.
4. Results
4.1 AFC transmission
A typical dielectric function calculated in the presence of an AFC is shown in Fig. 6. The AFC modifies only a small central region of the dielectric function: at large detunings, the AFC dielectric function approaches the unperturbed dielectric function. The key feature of the AFC dielectric function is its rapid oscillations with frequency.
The transmission spectra can be calculated from the dielectric function using Eq. (10). Figure 7 shows experimental results along with the corresponding theoretical spectra for both the rate equation and density matrix. The theories capture many features of the experimental transmission spectra, such as the widths and shapes of transmission windows between comb teeth, as well as transitions to non-periodic forms. The rate equations tend to predict more regular combs than the density matrix theory. We see that combs with large fundamental Fourier components occur a few hundred kHz above and below 4.7 MHz, as opposed to at 4.7 MHz where the combs are not dominated by the fundamental [78,79]. The effect of these features is discussed further in Section 4.2.
4.2 Photon echoes
The input pulse was modeled as a Gaussian in the time domain with 10 ns FWHM, corresponding to 88 MHz FWHM in the frequency domain. The pulse bandwidth is compared with the AFC bandwidth in Fig. 6, showing that it is a good match for the oscillatory region of the AFC dielectric function. After discarding spurious points arising from the artificial periodicity of the FFT, we obtain a time span of 10 $\mu$s, several times longer than the 1 $\mu$s time period in which the pulses appear. Figure 7 shows the experimental AFCs (top center) alongside density matrix (top left) and rate equation (top right) calculations, and the corresponding pulses (bottom). Figure 8 shows the echo pulses over a wider range of $f_{\rm rep}$ values. Although Fig. 6 is given for one of the density matrix calculations at a single repetition rate, we find the dielectric function separately for each theory at each repetition rate.
The initial pulse near 0 $\mu$s is present in both experiment and theory. The echo pulses appear on hyperbolas, with the $N^{\rm th}$-order echo emission occurring at $T=N/f_{\rm rep}$. Figure 7 highlights a number of features common to the experiment and the theory. Most notably, $f_{\rm rep}$ the largest first-order echoes do not align well with the excited-state hyperfine energy spacing values (dashed lines in panel 5): the largest echoes occur for $f_{\rm rep}$ detuned slightly above and below these energy spacings. Indeed, the echo efficiency at one of those spacings is near its lowest value. The patterns of drop-out of the first and second photon echoes is predicted well: for example, the absence of a photon echo for $f_{\rm rep}=2.5$ MHz is predicted by both theories and it is experimentally observed, despite photon echoes for both $f_{\rm rep}=2.4$ MHz and 2.6 MHz. Also common to both theories and experiment are valleys near 3.5 MHz and 4.5 MHz, although the drop-out is over-predicted by theory for the 4.5 MHz case. There are also small echoes in both experiment and density matrix theory for $f_{\rm rep}$ between 2.0 MHz and 3.5 MHz, at nearly fixed delay times between 0.20 $\mu$s and 0.25 $\mu$s. These are much less prominent in the rate equation calculation.
5. Conclusions
We present a semiclassical theory of photon echoes that gives a good account of the pulses arising in an AFC experiment in Pr:YSO with fine-grained systematic spectral coverage. The theory has two variants, one in which the AFC is predicted by the Rate Equations and a second in which the AFC is predicted by a Density Matrix. In contrast to most descriptions of AFC’s in the literature, our approach does not invoke the Dicke state [41]. The theory involves obtaining the dielectric function of the contribution of the Pr$^{3+}$ ions to the crystal from (a) a rate-equation or density-matrix description of the relevant hyperfine levels in the Pr$^{3+}$ impurity ions leading to the atomic polarizability (up to a constant) along with the input optical signals that prepared the AFC, (b) the Clausius-Mossotti relation, and (c) a fit to the maximum absorption in the case of a crystal without an AFC (thus fixing the constant). The photon echoes are modeled by letting the probe pulse pass through the AFC, here modeled as a highly dispersive medium. There does not appear to have been any similar attempt to simulate realistic atomic frequency combs based on the actual optical fields used to prepare them. Further, by generating over 100 AFCs for our systematic study, we demonstrate the predictive power of our theoretical model.
Looking ahead, using our predictive model and its natural extensions will enable tailoring the AFC preparation protocol to optimize the efficiency of the system as a quantum memory. The fact that our model is relatively simple in that it relies on long-established principles of the response of an oscillator to classical electromagnetic fields will hopefully enable applications.
Appendix
Signal averaging
The theory requires knowledge of the periodic signal entering the crystal. Here, we describe how that periodic signal is obtained from the RF recorded after heterodyne mixing, shown in Fig. 4(a). Signal averaging is used to extract the periodic signal from 20 $\mu$s samples of data discretized in 0.2 ns steps. The algorithm for averaging the signals is given below.
Results for the case of a modulation frequency $f_{\rm rep} = 4.8$ MHz are shown in Fig. 5, including the data before and after the averaging procedure. The instantaneous frequency shown in Fig. 5(c) may be compared to the intended waveform, an ideal triangle-wave modulation: the discrepancy is due to limitations of the wave form generator.
The quality of the averaged signal may be assessed by comparing its Fourier transform to the Fourier transform of the input signal, as shown for the example of 4.8 MHz in Fig. 3. The averaged spectrum is necessarily discrete. The experimental spectrum shows strong peaks that are broadened in analogy with x-ray-diffraction Bragg peaks broadened through thermal motions of the atoms, as described by the Debye-Waller factor. The alignment of the discrete points to the peaks of the black curve validates the averaging process. Matches of similar quality were obtained for each of the 161 modulation frequencies from 1 MHz to 16 MHz in steps of 0.1 MHz. The signal averaging algorithm has the following steps:
- 1. The mean is subtracted from the sample and the sample is smoothed with a 7-point symmetric binomial filter.
- 2. Locations of zero crossings (i.e., sign switches in adjacent smoothed data points) are found. If there are two zero crossings in consecutive time steps, the pair is not considered in the next step.
- 3. The instantaneous frequency $f\approx \frac {1}{2\pi } \frac {\Delta \phi }{\Delta t}$ is calculated by assigning a phase shift of $\Delta \phi =\pi$ to adjacent zero crossings which also give the time interval $\Delta t$.
- 4. The instantaneous frequency is averaged stroboscopically given the known driving period.
- 5. The envelope is obtained by extracting the analytic signal from the smoothed data of step 2 above. The analytic signal is obtained through discrete FFT followed by zeroing out the negative frequencies, followed by the inverse FFT. The envelope is the absolute value of the analytic signal.
- 6. The envelope is averaged stroboscopically to create a periodic function.
- 7. The periodic phase is obtained by integrating the averaged instantaneous frequencies obtained in step 3 over a period. A small linear adjustment in the phase is made to enforce strict periodicity of both amplitude and phase.
- 8. The signal, a complex quantity, is formed as the product of the envelope from step 6 and the phasor of the periodic phase from step 7. (The phasor of a phase $\phi$ is $\text {e}^{\text {i} \phi }$.)
- 9. The signal is shifted in the Fourier domain to remove the central modulating frequency. This corresponds physically to incorporating the central frequency of the AOM into the optical reference signal. The effect is to create a function that may be sampled with more widely spaced points. In practice, there is a 5:1 downsampling at this step from 0.2 ns to near 1 ns. The step sizes are not exactly equal for different repetition rates because each time step must be an integer divisor of the corresponding period.
Acknowledgments
We are grateful to Alexey Gorshkov, Robinjeet Singh, and Eite Tiesinga for helpful discussions and to Ivan Burenkov, and Sergey Polyakov for helping to build the experimental setup. Mention of commercial products does not imply endorsement by the authors’ institutions.
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. A. Acín, I. Bloch, H. Buhrman, T. Calarco, C. Eichler, J. Eisert, D. Esteve, N. Gisin, S. J. Glaser, F. Jelezko, S. Kuhr, M. Lewenstein, M. F. Riedel, P. O. Schmidt, R. Thew, A. Wallraff, I. Walmsley, and F. K. Wilhelm, “The quantum technologies roadmap: a European community view,” New J. Phys. 20(8), 080201 (2018). [CrossRef]
2. M. Pompili, S. L. N. Hermans, S. Baier, H. K. C. Beukers, P. C. Humphreys, R. N. Schouten, R. F. L. Vermeulen, M. J. Tiggelman, L. dos Santos Martins, B. Dirkse, S. Wehner, and R. Hanson, “Realization of a multinode quantum network of remote solid-state qubits,” Science 372(6539), 259–264 (2021). [CrossRef]
3. H.-J. Briegel, W. Dür, J. I. Cirac, and P. Zoller, “Quantum repeaters: The role of imperfect local operations in quantum communication,” Phys. Rev. Lett. 81(26), 5932–5935 (1998). [CrossRef]
4. L.-M. Duan, M. D. Lukin, J. I. Cirac, and P. Zoller, “Long-distance quantum communication with atomic ensembles and linear optics,” Nature 414(6862), 413–418 (2001). [CrossRef]
5. D. D. Awschalom, R. Hanson, J. Wrachtrup, and B. B. Zhou, “Quantum technologies with optically interfaced solid-state spins,” Nat. Photonics 12(9), 516–527 (2018). [CrossRef]
6. C. Simon, H. De Riedmatten, M. Afzelius, N. Sangouard, H. Zbinden, and N. Gisin, “Quantum repeaters with photon pair sources and multimode memories,” Phys. Rev. Lett. 98(19), 190503 (2007). [CrossRef]
7. M. Afzelius, C. Simon, H. de Riedmatten, and N. Gisin, “Multimode quantum memory based on atomic frequency combs,” Phys. Rev. A 79(5), 052329 (2009). [CrossRef]
8. M. Afzelius, I. Usmani, A. Amari, B. Lauritzen, A. Walther, C. Simon, N. Sangouard, J. Minář, H. de Riedmatten, N. Gisin, and S. Kröll, “Demonstration of atomic frequency comb memory for light with spin-wave storage,” Phys. Rev. Lett. 104(4), 040503 (2010). [CrossRef]
9. F. R. Graf, A. Renn, and U. P. Wild, “Site interference in stark-modulated photon echoes,” Phys. Rev. B 55(17), 11225–11229 (1997). [CrossRef]
10. F. R. Graf, A. Renn, G. Zumofen, and U. P. Wild, “Photon-echo attenuation by dynamical processes in rare-earth-ion doped crystals,” Phys. Rev. B 58(9), 5462–5478 (1998). [CrossRef]
11. M. Nilsson, L. Rippe, S. Kröll, R. Klieber, and D. Suter, “Hole-burning techniques for isolation and study of individual hyperfine transitions in inhomogeneously broadened solids demonstrated in Pr3+:Y2SiO5,” Phys. Rev. B 70(21), 214116 (2004). [CrossRef]
12. S. E. Beavan, E. A. Goldschmidt, and M. J. Sellars, “Demonstration of a dynamic bandpass frequency filter in a rare-earth ion-doped crystal,” J. Opt. Soc. Am. B 30(5), 1173–1177 (2013). [CrossRef]
13. H. Q. Fan, K. H. Kagalwala, S. V. Polyakov, A. L. Migdall, and E. A. Goldschmidt, “Electromagnetically induced transparency in inhomogeneously broadened solid media,” Phys. Rev. A 99(5), 053821 (2019). [CrossRef]
14. J. Klein, F. Beil, and T. Halfmann, “Robust population transfer by stimulated Raman adiabatic passage in a Pr3+:Y2SiO5 crystal,” Phys. Rev. Lett. 99(11), 113003 (2007). [CrossRef]
15. H. Goto and K. Ichimura, “Observation of coherent population transfer in a four-level tripod system with a rare-earth-metal-doped crystal,” Phys. Rev. A 75(3), 033404 (2007). [CrossRef]
16. H.-H. Wang, L. Wang, X.-G. Wei, Y.-J. Li, D.-M. Du, Z.-H. Kang, Y. Jiang, and J.-Y. Gao, “Storage and selective release of optical information based on fractional stimulated Raman adiabatic passage in a solid,” Appl. Phys. Lett. 92(4), 041107 (2008). [CrossRef]
17. S. Mieth, D. Schraft, T. Halfmann, and L. P. Yatsenko, “Rephasing of optically driven coherences by rapid adiabatic passage in Pr3+:Y2SiO5,” Phys. Rev. A 86(6), 063404 (2012). [CrossRef]
18. E. A. Goldschmidt, S. E. Beavan, S. V. Polyakov, A. L. Migdall, and M. J. Sellars, “Storage and retrieval of collective excitations on a long-lived spin transition in a rare-earth ion-doped crystal,” Opt. Express 21(8), 10087–10094 (2013). [CrossRef]
19. J. J. Longdell, E. Fraval, M. J. Sellars, and N. B. Manson, “Stopped light with storage times greater than one second using electromagnetically induced transparency in a solid,” Phys. Rev. Lett. 95(6), 063601 (2005). [CrossRef]
20. G. Heinze, C. Hubrich, and T. Halfmann, “Stopped light and image storage by electromagnetically induced transparency up to the regime of one minute,” Phys. Rev. Lett. 111(3), 033601 (2013). [CrossRef]
21. K. Kutluer, M. F. Pascual-Winter, J. Dajczgewand, P. M. Ledingham, M. Mazzera, T. Chanelière, and H. de Riedmatten, “Spectral-hole memory for light at the single-photon level,” Phys. Rev. A 93(4), 040302 (2016). [CrossRef]
22. M. P. Hedges, J. J. Longdell, Y. Li, and M. J. Sellars, “Efficient quantum memory for light,” Nature 465(7301), 1052–1056 (2010). [CrossRef]
23. D. Rieländer, K. Kutluer, P. M. Ledingham, M. Gündoǧan, J. Fekete, M. Mazzera, and H. de Riedmatten, “Quantum storage of heralded single photons in a praseodymium-doped crystal,” Phys. Rev. Lett. 112(4), 040504 (2014). [CrossRef]
24. A. Seri, D. Lago-Rivera, A. Lenhard, G. Corrielli, R. Osellame, M. Mazzera, and H. de Riedmatten, “Quantum storage of frequency-multiplexed heralded single photons,” Phys. Rev. Lett. 123(8), 080502 (2019). [CrossRef]
25. K. Mannami, T. Kondo, T. Tsuno, T. Miyashita, D. Yoshida, K. Ito, K. Niizeki, I. Nakamura, F.-L. Hong, and T. Horikiri, “Coupling of a quantum memory and telecommunication wavelength photons for high-rate entanglement distribution in quantum repeaters,” Opt. Express 29(25), 41522–41533 (2021). [CrossRef]
26. S. P. Horvath, M. K. Alqedra, A. Kinos, A. Walther, J. M. Dahlström, S. Kröll, and L. Rippe, “Noise-free on-demand atomic frequency comb quantum memory,” Phys. Rev. Res. 3(2), 023099 (2021). [CrossRef]
27. M. Bonarota, J. Ruggiero, J. L. L. Gouët, and T. Chanelière, “Efficiency optimization for atomic frequency comb storage,” Phys. Rev. A 81(3), 033803 (2010). [CrossRef]
28. I. Usmani, M. Afzelius, H. De Riedmatten, and N. Gisin, “Mapping multiple photonic qubits into and out of one solid-state atomic ensemble,” Nat. Commun. 1(1), 12–17 (2010). [CrossRef]
29. J. H. Davidson, P. Lefebvre, J. Zhang, D. Oblak, and W. Tittel, “Improved light-matter interaction for storage of quantum states of light in a thulium-doped crystal cavity,” Phys. Rev. A 101(4), 042333 (2020). [CrossRef]
30. E. Saglamyurek, N. Sinclair, J. Jin, J. A. Slater, D. Oblak, F. Bussieres, M. George, R. Ricken, W. Sohler, and W. Tittel, “Broadband waveguide quantum memory for entangled photons,” Nature 469(7331), 512–515 (2011). [CrossRef]
31. P. Jobez, N. Timoney, C. Laplane, J. Etesse, A. Ferrier, P. Goldner, N. Gisin, and M. Afzelius, “Towards highly multimode optical quantum memory for quantum repeaters,” Phys. Rev. Appl. 93(3), 032327 (2016). [CrossRef]
32. T. Zhong, M. Kindem, J. G. Bartholomew, J. Rochman, I. Craiciu, E. Miyazono, M. Bettinelli, E. Cavalli, V. Verma, S. W. Nam, F. Marsili, M. D. Shaw, A. D. Beyer, and A. Faraon, “Nanophotonic rare-earth quantum memory with optically controlled retrieval,” Science 357(6358), 1392–1395 (2017). [CrossRef]
33. M. F. Askarani, M. G. Puigibert, T. Lutz, V. B. Verma, M. D. Shaw, S. W. Nam, N. Sinclair, D. Oblak, and W. Tittel, “Storage and reemission of heralded telecommunication-wavelength photons using a crystal waveguide,” Phys. Rev. Applied 11(5), 054056 (2019). [CrossRef]
34. M. Businger, A. Tiranov, K. T. Kaczmarek, S. Welinski, Z. Zhang, A. Ferrier, P. Goldner, and M. Afzelius, “Optical spin-wave storage in a solid-state hybridized electron-nuclear spin ensemble,” Phys. Rev. Lett. 124(5), 053606 (2020). [CrossRef]
35. E. Saglamyurek, M. G. Pugibert, Q. Zhao, L. Giner, F. Marsili, V. B. Verma, S. W. Nam, L. Oesterling, D. Nippa, D. Oblak, and W. Tittel, “A multiplexed light-matter interface for fibre-based quantum networks,” Nat. Commun. 7(1), 11202 (2016). [CrossRef]
36. I. Craiciu, M. Lei, J. Rochman, J. G. Bartholomew, and A. Faraon, “Multifunctional on-chip storage at telecommunication wavelength for quantum networks,” Optica 8(1), 114–121 (2021). [CrossRef]
37. M. Sabooni, Q. Li, S. Kröll, and L. Rippe, “Efficient quantum memory using a weakly absorbing sample,” Phys. Rev. Lett. 110(13), 133604 (2013). [CrossRef]
38. P. Jobez, I. Usmani, N. Timoney, C. Laplane, N. Gisin, and M. Afzelius, “Cavity-enhanced storage in an optical spin-wave memory,” New J. Phys. 16(8), 083005 (2014). [CrossRef]
39. M. Nicolle, J. N. Becker, C. Weinzetl, I. A. Walmsley, and P. M. Ledingham, “Gigahertz-bandwidth optical memory in Pr3+:Y2SiO5,” Opt. Lett 46(12), 2948–2951 (2021). [CrossRef]
40. G. P. Teja, C. Simon, and S. K. Goyal, “Photonic quantum memory using an intra-atomic frequency comb,” Phys. Rev. A 99(5), 052314 (2019). [CrossRef]
41. R. H. Dicke, “Coherence in spontaneous radiation processes,” Phys. Rev. 93(1), 99–110 (1954). [CrossRef]
42. N. Sangouard, C. Simon, M. Afzelius, and N. Gisin, “Analysis of a quantum memory for photons based on controlled reversible inhomogeneous broadening,” Phys. Rev. A 75(3), 032327 (2007). [CrossRef]
43. A. V. Gorshkov, A. André, M. D. Lukin, and A. S. Sørensen, “Photon storage in λ-type optically dense atomic media. ii. free-space model,” Phys. Rev. A 76(3), 033805 (2007). [CrossRef]
44. T. Chanelière, H. Hétet, and N. Sangouard, “Quantum optical memory protocols in atomic ensembles,” Adv. Atom. Molec. and Opt. Phys. 67, 77–150 (2018). [CrossRef]
45. C. Thiel, R. Macfarlane, Y. Sun, T. Böttger, N. Sinclair, W. Tittel, and R. Cone, “Measuring and analyzing excitation-induced decoherence in rare-earth-doped optical materials,” Laser Phys. 24(10), 106002 (2014). [CrossRef]
46. J. Xiong, M. Colice, F. Schlottau, K. Wagner, and B. Fornberg, “Numerical solutions to 2D Maxwell–Bloch equations,” Opt. Quantum Electron. 40(5-6), 447–453 (2008). [CrossRef]
47. N. M. Arslanov and S. A. Moiseev, “Optimal periodic frequency combs for high-efficiency optical quantum memory based on rare-earth ion crystals,” Quantum Electron. 47(9), 783–789 (2017). [CrossRef]
48. W. Tittel, M. Afzelius, T. Chaneliere, R. L. Cone, S. Kröll, S. A. Moiseev, and M. Sellars, “Photon-echo quantum memory in solid state systems,” Laser Photonics Rev. 4(2), 244–267 (2009). [CrossRef]
49. B. Maksimov, Y. A. Kharitonov, V. Ilyukhin, and N. Belov, “Crystal structure of the Y-Oxysilicate Y [SiO 4] O,” Sov. Phys. Doklady 13, 1188 (1969).
50. K. Holliday, M. Croci, E. Vauthey, and U. P. Wild, “Spectral hole burning and holography in an Y2SiO5:Pr3+ crystal,” Phys. Rev. B 47(22), 14741–14752 (1993). [CrossRef]
51. R. W. Equall, R. L. Cone, and R. M. Macfarlane, “Homogeneous broadening and hyperfine structure of optical transitions in Pr3+:Y2SiO5,” Phys. Rev. B 52(6), 3963–3969 (1995). In order to enforce the f-sum rule exactly, we make the following small adjustments 0.55 → 0.552, 0.38 → 0.382, 0.07 → 0.066, 0.40 → 0.396, 0.60 → 0.596, 0.01 → 0.008, 0.05 → 0.052, 0.02 → 0.022, and 0.93 → 0.926. These corrections are smaller than the measurement uncertainty which is 0.01. [CrossRef]
52. E. D. Black, “An introduction to Pound-Drever-Hall laser frequency stabilization,” Am. J. Phys. 69(1), 79–87 (2001). [CrossRef]
53. We have automated the experiment: due to variations in the computer speed over the course of the hours-long measurements, there is a small variation in the number of measurements for each modulation frequency.
54. J. Minář, N. Sangouard, M. Afzelius, H. de Riedmatten, and N. Gisin, “Spin-wave storage using chirped control fields in atomic frequency comb-based quantum memory,” Phys. Rev. A 82(4), 042309 (2010). [CrossRef]
55. J. Jin, E. Saglamyurek, M. l. G. Puigibert, V. Verma, F. Marsili, S. W. Nam, D. Oblak, and W. Tittel, “Telecom-wavelength atomic quantum memory in optical fiber for heralded polarization qubits,” Phys. Rev. Lett. 115(14), 140501 (2015). [CrossRef]
56. K. Furuya, A. Nandi, and M. Hosseini, “Study of atomic geometry and its effect on photon generation and storage [invited],” Opt. Mat. Express 10(2), 577–587 (2020). [CrossRef]
57. A. Holzapfel, J. Etesse, K. T. Kaczmarek, A. Tiranov, N. Gisin, and M. Afzelius, “Optical storage for 0.53 s in a solid-state atomic frequency comb memory using dynamical decoupling,” New J. Phys. 22(6), 063009 (2020). [CrossRef]
58. J. Etesse, A. Holzäpfel, A. Ortu, and M. Afzelius, “Optical and spin manipulation of non-kramers rare-earth ions in a weak magnetic field for quantum memory applications,” Phys. Rev. A 103(2), 022618 (2021). [CrossRef]
59. T. Chanelière, J. Ruggiero, M. Bonarota, M. Afzelius, and J. LeGouet, “Efficient light storage in a crystal using an atomic frequency comb,” New J. Phys. 12(2), 023025 (2010). [CrossRef]
60. P. R. Berman and J.-L. Le Gouët, “Pulsed field transmission by atomic frequency combs and random spike media: The prominent role of dispersion,” Phys. Rev. A 103(4), 043723 (2021). [CrossRef]
61. O. Levy and D. J. Bergman, “Clausius-Mossotti approximation for a family of nonlinear composites,” Phys. Rev. B 46(11), 7189–7192 (1992). [CrossRef]
62. P. J. Mohr and W. D. Phillips, “Dimensionless units in the SI,” Metrologia 52(1), 40–47 (2015). [CrossRef]
63. A. Zangwill, Modern Electrodynamics (Cambridge University Press, 2013) p. 789.
64. R. Beach, M. D. Shinn, L. Davis, R. W. Solarz, and W. F. Krupke, “Optical absorption and stimulated emission of neodymium in yttrium orthosilicate,” IEEE J. Quantum Electron. 26(8), 1405–1412 (1990). [CrossRef]
65. R. M. Martin, “Comment on the calculation of electric polarization in crystals,” Phys. Rev. B 9(4), 1998–1999 (1974). [CrossRef]
66. Z. H. Levine and D. C. Allan, “Quasiparticle calculation of the dielectric response of silicon and germanium,” Phys. Rev. B 43(5), 4187–4207 (1991). [CrossRef]
67. R. Resta, “Macroscopic polarizability in crystalline dielectrics,” Rev. Mod. Phys. 66(3), 899–915 (1994). [CrossRef]
68. H. Sõnajalg and P. Saari, “Diffraction efficiency in space-and-time-domain holography,” J. Opt. Soc. Am. B 11(2), 372–379 (1994). [CrossRef]
69. R. Unanyan, M. Fleischhauer, B. W. Shore, and K. Bergmann, “Robust creation and phase-sensitive probing of superposition states via stimulated raman adiabatic passage (stirap) with degenerate dark states,” Opt. Commun. 155(1-3), 144–154 (1998). [CrossRef]
70. G. Grynberg, A. Aspect, and C. Fabre, Introduction to Quantum Optics: From the Semi-Classical Approach to Quantized Light (Cambridge University Press, 2010).
71. M. B. Plenio and P. L. Knight, “The quantum-jump approach to dissipative dynamics in quantum optics,” Rev. Mod. Phys. 70(1), 101–144 (1998). [CrossRef]
72. J. R. Johansson, P. D. Nation, and F. Nori, “QuTiP: An open-source Python framework for the dynamics of open systems,” Comp. Phys. Comm. 183(8), 1760–1772 (2012). [CrossRef]
73. Y. Attal, P. Berger, L. Morvan, P. Nouchi, D. Dolfi, T. Chanelière, and A. Louchet-Chauvet, “Rate equation reformulation including coherent excitation: application to periodic protocols based on spectral hole-burning,” J. Opt. Soc. Am. B 35(6), 1260–1270 (2018). [CrossRef]
74. J. J. Longdell, M. J. Sellars, and N. B. Manson, “Hyperfine interactions in ground and excited states of praseodymium-doped yttrium orthosilicate,” Phys. Rev. B 66(3), 035101 (2002). [CrossRef]
75. M. Lovrić, P. Glasenapp, and D. Suter, “Spin Hamiltonian characterization and refinement for Pr3+:YAlO3 and Pr3+:Y2SiO5,” Phys. Rev. B 85(1), 014429 (2012). [CrossRef]
76. J. G. Bartholomew, R. L. Ahlefeldt, and M. J. Sellars, “Engineering closed optical transitions in rare-earth ion crystals,” Phys. Rev. B 93(1), 014401 (2016). [CrossRef]
77. E. Anderson, Z. Bai, C. Bischof, S. Blackford, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, and D. Sorensen, LAPACK Users’ Guide (Society for Industrial and Applied Mathematics, 1999), 3rd ed.
78. A. M. Weiner and D. E. Leaird, “Generation of terahertz-rate trains of femtosecond pulses by phase-only filtering,” Opt. Lett. 15(1), 51–53 (1990). [CrossRef]
79. K. Slattery and Z. H. Levine, “Simulating photon echoes for quantum memory,” https://discretionary/demonstrations.wolfram.com/SemiclassicalSimulationOfPhotonEchoesForAtomicFrequency/discretionary-Comb; https://data.nist.gov/od/id/mds2-2454 (2021).