Strong excitation of emitters in an impedance matched cavity: the area theorem, \pi-pulse and self-induced transparency

I theoretically study the behavior of strong pulses exciting emitters inside a cavity. The ensemble is supposed to be inhomogeneously broadened and the cavity matched finding application in quantum storage of optical or RF photons. My analysis is based on energy and pulse area conservation rules predicting important distortions for specific areas. It is well supported by numerical simulations. I propose a qualitative interpretation in terms of slow-light. The analogy with the free space situation is remarkable.


Introduction
Electromagnetic resonators offer the possibility to enhance the interaction of waves with matter.In the RF or optical domain, they are the basis of many commercial devices, going from the microwave-oven to lasers.At the most fundamental level, they permit to explore the boundary between the classical and the quantum world allowing experimental tests of quantum electrodynamics (QED) [1].
The practical interest for active resonators, i.e. filled with emitters, has been reactivated in the context of quantum storage for which a complete mapping of the field carrying the information into a long-lived atomic system is necessary.Between the two extreme situations, a single emitter in an ultra-high-finesse cavity (Cavity QED) on one side and a strongly absorbing medium in free-space on the other side, an intermediate regime exists.A weakly absorbing sample can be placed in a medium-finesse cavity to obtain a significant interaction.
I will consider the specific situation of an impedance matched ring cavity.The matching condition means that the input mirror transmission equals the losses.I assume the loss mechanism being dominated by the emitters absorption of the active medium inside the resonator.Conversion into the atomic excitation either way represents a loss or a mapping of the incoming field depending on the point of view.This approach has been particularly fruitful in non-linear optics for second harmonic generation or more generally for frequency mixing [2].The incoming power is then fully converted into the target frequency.Following the same idea, the use of a matched cavity has been proposed for quantum light storage [3,4] and successfully implemented in a luminescent crystal [5].The incoming signal is then fully mapped into the atomic coherences even if the single path absorption is moderate as soon as it is compensated by the cavity quality factor.This situation is not restricted to optics but is also considered very actively in the RF domain with spins [6,7,8,9,10].Similar proposals have been made to store and retrieve microwave photons [11,12].A resonator, in the optical or in the RF domain, is indeed advantageous as compared to free-space propagation.That's the reason why this technique is actively investigated in a topical context.It allows the use of low absorption samples reducing the constraint on the physical size of the medium or the concentration of emitters.Highly doped samples can indeed exhibit a important coupling between individual dipoles thus reducing the coherence time.It is the case for example for NV centers in diamond [13] and potentially for different kinds of impurities in insulators as rare-earth or transition ions.This statement can be generalized to atomic vapors.Extreme optical depth can be obtained in Bose-Einstein condensate at the price of interaction or collisions between dipoles.Low doping samples in a resonator are clearly an alternative.Using optically thin samples also facilitates the preparation stage as described by Afzelius et al. [3] when optical pumping can be operated from the side of the cavity.In the optical or RF domain, involving optical dipoles or spins, matched resonators offer interesting prospects.
The different protocols [3,4,11,12,14,15] to store and retrieve quantum signals are based on the coherent manipulation of spins or optical dipoles.A π-pulse is a tool of choice in that case because it permits the delayed emission of the stored signal in the lineage of the spin or photon echo experiments.The propagation of such strong pulses in free-space absorbing samples is a complex problem known for decades [16,17].More recently, the use of the area theorem has been proposed in the context of quantum memories [18].This general approach allows the derivation of analytic solutions for weak and strong area pulses.Later on , the distortion though the sample of the rephasing π-pulse has been identified experimentally as critical to explain the observed storage efficiency [19].My main motivation is to transpose this propagation analysis in a resonator filled with absorbers.
For the rest of the paper, I study the ring cavity design (see fig. 1).The ring configuration has the advantage of the theoretical simplicity.When the round-trip absorption is small, the traveling wave inside interacts only with one atomic mode.As a consequence, inside the resonator, the field and the atomic variables (population and polarization) are fully described by spatially independent parameters (the amplitudes of the field and atomic modes).This model has been extensively and accurately discussed in the context of optical bistability [20,21,22] based on previous work for the semi-classical description of lasers [23].
Even if my paper is restricted to the ring resonator design, it should be noted that the standingwave configuration has the advantage of the experimental simplicity.It only involves two mirrors thus reducing the unwanted passive losses, allowing more compact and stable setup.Semimonolithic linear resonators are indeed widely used in optics for pulse compression [24] or non-linear frequency conversion [2].The theoretical treatment is more complex because the standing wave induces a population grating in the non-perturbative regime.It couples the forward and backward propagating modes fundamentally differing from the traveling wave description [20,21].It is beyond the scope of the present paper but it should deserve a specific study because of its experimental simplicity and its current extensive usage for electronic spins coupled to planar waveguides [7,8,9,10].
I finally assume the atomic transition to be inhomogeneously broadened (local shift of the transition).It allows the use of photon echo techniques to control the sample and possibly offers a larger intrinsic interaction bandwidth.I assume the coupling constant being identical for all atoms.The assumption is well suited in optics because the orientation of the transition dipoles is well defined with respect to the incoming polarization.It is generally not appropriate for spins in planar waveguides because the magnetic dipoles can have different orientations and the magnetic field orientation around the waveguide has a radial dependency [25].Even if I employ the optics terminology and use appropriate assumptions my approach can hopefully open some perspectives for RF or microwave resonators offering a whole bestiary of integrated designs with waveguides on different substrates [25].
The aim of this paper is manifold.I primarily study the input/output relations of strong pulses in the matched cavity and point out the possible distortions.I also draw an analogy with the free space propagation through thick sample giving a solid basis for the intuition.I use the intracavity version of the McCall & Hahn area theorem that has been already derived by Zakharov in the context of optical bistability and superraddiance [21,22,26,27].The area theorem is extremely helpful to predict the qualitative behavior of strong pulses, π/2 or π for example.π-pulses for example are fundamental tools for the manipulation of quantum systems.It is still widely unknown in the quantum information community.My goal is to replace the work of Zakharov in this context and show how helpful it is to estimate the possible distortions when the pulse travels through the medium.The analytic part is then essentially a recontextualization of previous results.My numerical calculations of the pulse shapes are strongly supported by the area and energy conservation rules whose simplicity reinforce this kind of analysis.I propose a qualitative interpretation in terms of slow-light.This explanation has the advantage to reintroduce the dispersion properties of the medium particularly useful in the weak signal limit.My approach is clearly inspired by the problematic of optically active emitters in a resonant cavity.Input/output relations are well-known in optics and are formally equivalent to a propagation problem in free space.It is less natural in the microwave domain (mostly in EPR spectroscopy for which relatively large coupling strength can be obtained), because the wavelength is usually larger than the sample.Considering the complete mapping between a microwave or optical photon into spins or optical dipoles respectively brings two fields of study closer.Bridging this gap between optics and RF is also an objective of this paper and is then intended to a larger audience.

Ring cavity model
I consider a ring cavity uniformly filled with emitters (fig.1).In this geometry, the outgoing field Ω out and the incoming one Ω in can be separated easily.The fields interact with only one mode of the cavity whose amplitude Ω uniquely characterizes the traveling-wave inside the resonator.All the fields are then completely defined by their time-varying Rabi frequencies.Fig. 1.Active ring cavity uniformly filled with emitters.The entrance mirror is partially reflecting.The different fields are defined by their time-varying Rabi frequencies: Ω in , Ω out represent the incoming and the outgoing amplitudes on the entrance mirror and Ω the intracavity mode amplitude.The spatial dependence can be neglected when the round-trip absorption is small.The local interference because of a partial beam overlap is supposed to be negligibly within the cavity volume.

Atoms dynamics and cavity master equations
The coupled system between the cavity field and the emitters is described on one side by the input/output relations of the resonator including the atomic polarization (source term) and on the other side by the Bloch equation characterizing the dipole dynamics.This model is wellestablished in quantum optics [28,29,30,11,12] including sophisticated description of multilevel atomic level scheme [31] and complete quantum description of QED cavity with inhomogeneous broadening [4].Since I specifically study the strong pulse distortions, I use a semi-classical formalism for the cavity master equation and the atomic variables dynamics.It can be naturally extended to describe the interaction at the quantum level [4,11,12,31,32].
Concerning the dynamics of the emitters, the transition is assumed to be only inhomogeneously broadened.The dipole dephasing and the population lifetime are neglected, in other words the atomic response is examined in the coherent transient regime as in echo type experiments [16].In the Bloch vector formalism, the transverse and longitudinal decays are then neglected.In the rotating frame the equations read as: ∆ is the detuning from the excitation frequency and is used as an index accounting for the inhomogeneous broadening.U ∆ and V ∆ are the in-phase and out-of-phase components of the Bloch vector.The real Rabi frequency Ω is the slowly time-varying envelope of the monochromatic field.The last equation describes the population dynamics W ∆ .
Using the complex amplitude R ∆ = U ∆ + iV ∆ , the dipole moment's equation may be written in a more compact manner: The cavity master equations also known as input/output relations are slightly modified as compared to the empty ring cavity case [30].They include a source term corresponding to the atomic emission: The source term is the sum of the polarizations R ∆ over the inhomogeneous broadening g(∆).I take g(∆) = 1 over the spectrum of interest as the normalization assuming the inhomogeneous broadening is much larger than the interaction bandwidth.It is an important assumption.Even if numerical simulations can be implemented without this condition, it is required by the area theorem as we will see in 3.2.Since my paper precisely compares the intuition from the area conservation law and numerical simulations of the temporal shape, I make this assumption from now.To summarize, the dynamics is assumed to verify: homogeneous linewidth ≪ interaction bandwidth ≪ inhomogeneous broadening [33, p.3].
It is clearly a restriction because to optimize the resources one would naturally match the bandwidth and the inhomogeneous linewidth as in [7,10] for example.In that case, the analytic area theorem is not valid anymore, one can only rely on numerical simulations.
Coming back to the model description (eq.3), the atoms-light coupling constant is directly included into α, the measured absorption coefficient.The convergence of the integral is ensured by R ∆ whose support is the excitation bandwidth.
The term αL equals the round-trip absorption of the uniformly filled resonator.It can be rescaled by the filling ratio if the sample is smaller than the cavity.
The cavity is described by the free-spectral range defined as D = c L where L is its round-trip length and the intensity transmission coefficient κ of the entrance mirror.
I now focus on the specific condition where weak incoming pulses are completely mapped into the ensemble.

Matched cavity condition
The matching condition is defined in the perturbative limit: the populations are not modified (weak pulse) and are supposed to keep its initial value W ∆ = −1.I define the reflection coefficient in the spectral domain Ω out = r (ω) Ω in where Ω out and Ω in are the Fourier transform of the incoming and outgoing pulses respectively.r (ω) completely characterizes the response of the cavity in the weak signal limit.It reads as The cavity is matched when the reflection is zero on resonance (ω = 0): i.e when the round-trip absorption equals the inverse of the finesse 2π κ .
For a matched cavity, the intensity reflection |r (ω) | 2 exhibits a dip whose FWHM defines the cavity linewidth ∆ω cav .It scales the possible interaction bandwidth as previously demonstrated by Moiseev et al. [4,15].It is then an important parameter: ∆ω cav = 2κD.
Because I neglect the atomic decay, the field is completely mapped into the dipoles without damping.It allows the implementation of the field-matter interface.The system is clearly conservative so it is particularly helpfull to use the energy conservation rule to evaluate the conversion from field to atoms.

Energy conservation
In the coherence propagation regime, there is a conservation of quanta between the field (photons) and the excitation of the ensemble (population).Writing the equations by using the Rabi frequencies to describe the field envelope and the absorption coefficient (eqs. 1 and 3) hides the microscopic parameters (energy quanta) but has the advantage to exhibit experimentally measurable macroscopic parameters.The conservation of quanta can be straightforwardly retrieved by integrating the equation of motion [16].
On the onde side, the energy contained in the incoming field (electromagnetic) is given by the integration of the intensity i.e the incoming pulse energy On the other side, the sum of the population gives the energy stored into the atomic excitation.
It reads with the same units U w = αL Comparing the two forms of energy is a simple way to predict the qualitative shape of the output pulse as we'll see later.This analysis tells us how much energy is potentially left behind by the incoming pulse.The area conservation rule offers a second powerful tool.

Intracavity area theorem
The McCall-Hahn area theorem is extremely helpful in free space to guide the intuition when propagation in absorbing media is considered [17].It gives a simple analytic law for the pulse area defined as Θ = t Ωdt.
Its intracavity version has been derived by Zakharov [21,22].It has been studied experimentally with doped solids [26,27].I briefly rederive it here to make the present article more self-contained.I recommend the reading of [33] for a pedagogical introduction and accurate derivation to the area theorem in free space (see also [16, p.19] and [34, p.816]) .The intracavity version follows the same formal demonstration [21].
It is given by a time integration of the equation of motion (eq.3): The source term is resolved by rewriting the Bloch equation ( 2) in its integral form [33].
Following [33] because the calculation of the source term contribution is the same as in free space, the distribution g(∆) is assumed much broader than the interaction bandwidth so the term exp(−i∆(t ′ − t)) is rapidly oscillating.This simplification is possible because of the three previously mentioned well-separated temporal dynamics: homogeneous linewidth ≪ interaction bandwidth ≪ inhomogeneous broadening.The integration of the field Θ = t Ωdt introduces a time when the field disappears and the macroscopic polarization vanishes because of the inhomogeneous dephasing but the coherences still freely oscillate.This oscillating term is a formal representation of the Dirac distribution centered on ∆ = 0 [16,34].After integration over the inhomogeneous profile, only remains W 0 (t), the on-resonance population [33, eq. ( 7)].It is the solution of the Rabi flopping problem: I finally obtain the area theorem It gives input/output relations for the pulse area whatever is the exact incoming temporal shape.
The sin (Θ) term introduces singular propagation effects for specific area integer of π as in the free-space situation [35].Assuming a matched cavity (eq.5) further simply the expression: In the small area limit (perturbative regime), one verifies that Θ out = 0.No light escapes the cavity (matching condition).For incoming area ranging from 0 to √ κπ, I represent the outgoing and intracavity area (fig.2).The latter ranges from 0 to 2π accordingly.
Θ/π Fig. 2. Outgoing Θ out (top) and intracavity area Θ (bottom) as a function of the incoming area Θ in calculated from the matched cavity area theorem (eqs.10).Squares correspond to the numerical simulation of the pulse temporal shapes that will be detailed later on (see 4.3).It serves as a validation of the numerically calculated area that can be compared to the analytic result of the area theorem.
I now analyze the two singular situations for the area theorem: Θ = π and Θ = 2π.For both and more generally as any integer of π, the area is conserved.They are essential for the manipulation of coherences in the context of quantum information storage.

Pulse distortion
I here show that energy and area conservation rules can predict the general behavior of strong pulses through the cavity.The intuition is well supported by numerical results.
In other words, the area is conserved.It may be surprising at very first sight because the intracavity π-pulse is inverting the medium and is then leaving some energy behind.This apparent contradiction actually tells us that the pulse stretches to conserve its area but to reduce its energy.My point is actually to look at the energy conservation to evaluate qualitatively these distortions.If the pulse looses are a neglegible fraction of its energy, it will keep its shape.κ∆ω for the incoming energy.
On the other side, if the medium is inverted (W ∆ = 1) over a bandwidth ∆ω, the population αL∆ω.For a matched cavity (eq.5), U Ω and U w have precisely the same scaling with κ∆ω.The ratio doesn't depend on the atomic or cavity parameters.
It should be noted that the two forms of energies have the same order of magnitude precisely because the matching condition is assumed.On the contrary, the area theorem in general (eq.9) doesn't depend on the absorption for Θ = π.In other words, a π-pulse conserves its area if the cavity is matched or not.It offers a certain degree of freedom if the pulse distortions have to be minimized.If the transmission of the input mirror can be changed for example.Even if the matching is required for a complete storage of the incoming signal, a dynamical control of the resonator parameters as it is the case for RF waveguides would certainly open some perspectives.

Discussion
An equivalent energy comparison can be done for an arbitrary area.It tells that the pulses have to leave a constant fraction of energy behind independently from their duration.Let's take weak pulses as a reference.They reduce both area and energy by simply reducing the pulse amplitude keeping the duration constant.This option is not possible for a π-pulse because the area should be conserved.The only possibility here is to associate amplitude reduction and pulse elongation.A π-pulse reduces the energy but roughly keeps its area.Major distortions of the incoming pulse are then expected.The analogy with strongly absorbing media in free-space situation is remarkable [35].
Before considering numerical simulations which are necessary to obtain the exact shape of the outgoing pulse, I treat the case of Self-Induced Transparency (SIT).

Self-induced transparency of 2π-pulses
The situation is clearly different for Θ in = √ κπ.The area should be still conserved but the energy conservation is relaxed.The on-resonance atoms are returned back to the ground state.Even if off-resonant atoms can still be excited because they don't undergo a 2π area, less energy should be left behind as compared to a π-pulse.Because the energy conservation is less drastic, less distortions and a minimized absorption are likely.Soliton like behavior are expect in analogy with free-space defining the SIT phenomenon [16,17].Since the present paper is primarily focused on quantum storage, 2π rotations of the atomic state are not really interesting because they essentially leave the system unchanged.π-pulses are more relevant in that sense.They can indeed rephase the inhomogeneous broadening (photon or spin echo) and be used in series for dynamical decoupling sequences [36,37] 4.3.Numerical simulation of the pulse temporal shapes The numerical simulation first requires a discretization of the detuning ∆.I write the Bloch system of eq. ( 1) for n evenly-spaced detunings ∆ n (spacing d∆) leading to 3n equations for The cavity master (eq. 3) becomes It gives (3n + 1) linear first-order differential equations for the variables The boundary conditions are given on the one hand by assuming the atoms in the ground state (U ∆ n ,V ∆ n ,W ∆ n ) = (0, 0, −1) and on the other hand by giving an incoming temporal pulse shape Ω in (t).I solve the (3n + 1) linear system numerically (n = 256) by applying a forth-order Runge-Kutta method.
For numerical application, I choose a finesse of 2π κ = 500 and a free-spectral-range D = 3GHz.These common parameters covers different physical realities.In the optical domain, they correspond to a 10 cm long cavity for which finesse of 1000 have been observed [38,39] using monolithic samples.In the RF domain, a GHz range can be associated with the hyperfine splitting of NV centers in diamond [7,8,10] or with the electron spin of paramagnetic impurities in luminescent crystals [6,9].The finesse in that case is not limited by the ultimate performance of the superconducting resonator but by the atomic linewidth (typical in the MHz range) then limiting the quality factor (few hundreds).A finesse of 2π κ = 500 is then realistic and covers a wide range of physical systems.The cavity is supposed to be matched (eq 5).The different parameters are now fixed.

Outgoing pulse shapes
I choose gaussian shaped incoming pulses for the simulation with a 2µs duration.The corresponding bandwidth is 2π × 80kHz much smaller than the cavity linewidth ∆ω cav = 2π × 12MHz.I keep the duration constant and vary the pulse amplitude ranging the incoming area from zero to Θ in = √ κπ (2π-pulse).I then obtain the temporal shape of the intracavity and outgoing pulses in fig. 3. I also perform a rudimentary test of my simulation by computing the intracavity Θ and outgoing Θ out areas.They can be easily compared (squares on fig.2) to the analytic result of the area theorem (eq.10).
The pulse for Θ in = 0.4 √ κ 2 π behaves as a weak area pulse meaning it is mostly absorbed without distortion.No reflection at all is expected for areas much smaller than √ κ 2 π.The 2π-pulse with Θ in = √ κπ shows remarkable similarities with SIT in strongly absorbing media.It conserves its amplitude and shape.It is essentially delayed.The delay scales as the incoming pulse duration, 2µs in that case [16].The situation is particular and has been observed experimentally [26,27].It is less interesting in the context of quantum manipulation of qubits but it would certainly deserve further investigation.I now consider more specifically the π-pulse situation.
The outgoing shape Ω out for Θ in = √ κ 2 π is characterized by a long tail, much longer than the pulse duration.This strong distortion is precisely expected for π-pulse which have to leave most of its energy behind by conserving its area (see section 4).To make this comparison more clear, I plot the normalized outgoing pulse (fig.4, left).
As previously mentioned, the 2π outgoing pulse (SIT) is remarkably similar to the incoming pulse (fig.4, left, blue curve).It is essentially delayed.The almost weak area pulse Θ in = 0.4 4, left, black curve) is not distorted neither.For a matched cavity, it is strongly absorbed, the observed shape and delay should be taken with precaution because they strongly depend on the numerical parameters.Temporal discretization may typically explain the pedestal on the rising edge of outgoing pulse (fig.4, black curve).Anyway, it serves as a reference as compared to π-pulse whose behavior is completely different.A long tail clearly appears.red and blue respectively (as in fig.3).The normalised incoming pulse is plotted as a reference (dashed red).The long tail of the outgoing π-pulse (red) is prominent.Right: root mean square (rms) temporal widths of the calculated pulse (squares, a dashed line is used to guide the eye).

Outgoing pulse width
To quantify the distortion, I plot the rms width σ of the pulse (fig.The pulse distortion is usually associated with a steep dispersion profile.The slow-light phenomenon recently reappeared in the context of impedance matched cavity [40] with particularly strong and observable effects.Considering the reflection coefficient (eq.4) as an effective susceptibility whose imaginary part is the equivalent refractive index is tempting.It would be only justify in the weak signal perturbative limit and cannot be used in general to explain the outgoing shape of strong area pulses.I nevertheless propose to consider the effective group delay for a qualitative analysis essentially demonstrating that extreme dispersion is expected for an inverted medium.One can push the perturbative approach by looking at the reflection coefficient (eq.4).I assume for example that a given time the atomic population is excited W ∆ (t) ≥ −1, and from that point calculate an effective reflection coefficient.I now assume the population constant W over the bandwidth of interest for the sake of simplicity since my goal is to derive an order of magnitude.The population simply rescales the absorption coefficient so to say taking the form −W × α.It leads to gain or absorption if the medium is inverted or not.Starting from the generalized expression of the reflection coefficient One can derive the group delay T g for a matched cavity (eq.5) by a first order expansion r W (ω) = r W (0) + iωT g + ... One simply finds I now compare the initial situation in the ground state (the population is weakly modified by small area pulses) and a inverted medium as expected for π-pulses.For atoms in the ground state, the reflection tends toward zero and the group delay toward 1 ∆ω cav , that can be compared to the empty ring cavity case with a 4 ∆ω cav delay.When the medium is inverted, both the reflection r W (0) and the group delay diverge.It is extremely surprising at first sight because an infinite reflection seems to break the energy conservation rule.It is not the case because it is precisely associated with an infinite group delay.It tells that the weak pulse stays for an infinitely long time in the cavity.Because the ensemble is supposed inverted (and not depleted), it gives constantly its energy to the pulse leading to an endless amplification r W (0) → ∞.This analysis should not be confused with the distortion of a π-pulse.A strong pulse enters an notinverted medium but creates the inversion through which it propagates.It cannot be accurately described by my perturbative expansion and precisely requires a numerical simulation to solve the coupled equations.We nevertheless understand qualitatively that the more the pulse inverts the medium, the slowest is goes through the sample.It is not so surprising at the end that a long tail appears following the outgoing π-pulse.

Conclusion
The intracavity and output shapes of a strong exciting pulse are extremely sensitive to its incoming area close to Θ in = √ κ 2 π.It corresponds to a singularity for the intracavity area theorem.This effect can be interpreted as a competition between the area and energy conservation rules.The pulse stretches because it reduces its energy by keeping its area.The parallel with the free space propagation is remarkable.This effect appears to be critical to explain the observed efficiency of the two-pulse echo in an optically-thick sample [19].A similar result is expected for echo type experiments in a matched cavity.
A long pulse tail can be problematic when a weak signal has to be isolated from a strong control pulse.To get around this potential limitation, monochromatic π-pulses can be replaced advantageously by adiabatic inversion as proposed in free-space [41].Frequency swept pulses are not only more robust in terms of power fluctuation for example but they are known to undergo less distortion when propagating though absorbing samples [42].A similar behavior is expected in a matched-cavity.Frequency swept pulses would deserves a specific study.The Bloch-Maxwell model can be extended to include the time-dependent phase of the field in a straightforward manner [16].It should be noted that the area theorem is not valid anymore, in a certain sense relaxing the constraints on the pulse dynamics [33].

4. 3 . 4 .
Slow-light interpretation It is not the case here.Let's compare the incoming energy and what's left in the atomic excitation by a π-pulse of bandwidth ∆ω (see 3.1).Its amplitude Ω in √ κ 2 π by definition of a π-pulse.So I obtain U Ω ∼ π 2 4