Abstract
We present an active feedback scheme acting continuously on the state of a quantum gas dispersively coupled to a high-finesse optical cavity. The quantum gas is subject to a transverse pump laser field inducing a self-organization phase transition, where the gas acquires a density modulation and photons are scattered into the resonator. Photons leaking from the cavity allow for a real-time and non-destructive readout of the system. We stabilize the mean intra-cavity photon number through a micro-processor controlled feedback architecture acting on the intensity of the transverse pump field. The feedback scheme can keep the mean intra-cavity photon number nph constant, in a range between nph = 0.17(4) and nph = 27.6(5), and for up to 4 s. Thus we can engage the stabilization in a regime where the system is very close to criticality as well as deep in the self-organized phase. The presented scheme allows us to approach the self-organization phase transition in a highly controlled manner and is a first step on the path towards the realization of many-body phases driven by tailored feedback mechanisms.
Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Feedback on quantum gases
Ultracold atomic quantum gases are a well-suited platform to study transitions and crossovers between different phases of matter. Prominent examples are the phase transition from a thermal gas to a Bose–Einstein condensate [1, 2] and between a superfluid and a Mott insulator [3], or the crossover between a Bose–Einstein condensate of molecules and a Bardeen–Cooper–Schrieffer superfluid of loosely bound pairs in quantum degenerate Fermi gases [4–6]. Another well studied example is the transition to a superradiant or self-organized phase in the driven-dissipative Dicke model [7].
A characteristic of cold atom experiments is that most probing techniques are inherently destructive [8]. This requires experiments to be repeated many times with the same initial conditions, such that statistically significant findings can be derived. The successive preparation of the system with identical parameters proves challenging and therefore often requires postselection or additional procedures during state preparation. In a recent experiment, non-destructive Faraday imaging was employed to measure the atom number and adjust the subsequent optical evaporation to prepare a number stabilized ultracold atomic cloud [9]. In other experiments, the transmission spectrum of a cavity was monitored to trigger the start of the experiment once a set atom number, detected via the dispersive shift of the cavity resonance, was reached [10, 11].
Dispersively coupling an atomic gas to a cavity comes with a non-destructive measurement channel. The continuous photon leakage through the cavity mirrors can be used to monitor the evolution of the system in real time. Starting point for our feedback scheme is an experimental setup in which we prepare a degenerate Bose gas inside a high-finesse optical cavity. By tuning the strength of an external drive field, the combined atom-cavity system is undergoing a self-organization phase transition which can be mapped onto the Dicke model, which is of fundamental importance in quantum optics [12–15] and has been recently experimentally realized in cold atoms [7]. The strong collective coupling between a quantum gas and the cavity field effectively mediates interactions between the atoms. This mechanism has been used to study e.g. competing short- and long range interactions [16, 17], dissipation-induced instabilities [18] and to engineer dynamical spin–orbit coupling [19]. The intra-cavity light field is related to the order parameter of the self-organization phase transition, providing in principle real-time access to critical properties such as the fluctuations of the order parameter. However, since the system close to the critical point is highly susceptible to minute drifts of any parameter, such a measurement is challenging and requires many experimental runs and a sophisticated data analysis [20, 21]. Here, we stabilize the state of the system via detecting the photons leaking through the cavity mirrors and applying an according feedback signal to the control parameter. This way, the lifetime of the self-organized phase close to the phase transition can be dramatically extended.
Going a step further, feedback could be used to control phase transitions and engineer novel non-equilibrium phases in driven many-body systems. The combination of feedback and weak measurements promises to realize new feedback-induced phase transitions and control of their critical exponent [22]. The simulations of interesting many-body problems like spin-bath models, Ising type interactions, the Lipkin–Meshkov–Glick model and Floquet time crystals through specifically engineered feedback schemes have also been proposed [22]. Applying time-delayed feedback to the driven-dissipative Dicke model opens the prospect to study non-equilibrium dynamics with fixed points and limit cycles in the superradiant regime [23].
In this work, we report on the basic building block to implement such feedback schemes and present a feedback architecture designed to stabilize the mean intra-cavity photon number. We first review the self-organization phase transition in our system. Following, we present the technical implementation of our feedback scheme. We then compare experimental realizations with and without the feedback and demonstrate stabilization of the mean intra-cavity photon number in a wide range.
2. Self-organization phase transition of an atomic gas in a cavity
We prepare a 87Rb Bose–Einstein condensate (BEC) inside a high-finesse optical cavity and illuminate it with a retroreflected standing wave transverse laser pump (TP). The frequency of the TP is red detuned from atomic resonance. Both the exact detuning between the TP frequency and the cavity resonance frequency as well as the intensity of the TP are adjustable. The intensity of the TP controls the depth of the generated standing wave lattice. The setup is sketched in figure 1.
Figure 1. Sketch of the experimental setup. An atomic gas (blue) inside an optical cavity is illuminated with a retroreflected standing wave transverse laser pump (TP). The system undergoes a self-organization phase transition for a sufficiently high intensity of the TP. In the self-organized phase, the atomic density is modulated and an intra-cavity light field is built up. The intra-cavity light field leaks through the cavity mirrors and is recorded by a heterodyne detector (green box). The signal processing unit calculates the present mean intra-cavity photon number and the feedback signal required to stabilize it to a set value (orange box). The feedback signal is sent to the intensity control of the TP (blue box). When the feedback is not enabled, an external control signal is sent to the intensity control of the TP.
Download figure:
Standard image High-resolution imageThe system can be described by the Hamiltonian of the Dicke model [7]:

The first term in equation (1) is the Hamiltonian of the cavity light field in the rotating frame of the TP. Here, Δc is the detuning between pump and cavity frequency neglecting the dispersive shift and (
) is the annihilation (creation) operator of the cavity mode. The second term describes the bare energy of the atomic part. It is based on a two-mode description of the atomic cloud in momentum space. The ground state
corresponds to the BEC. Here, px (pz) refers to the momentum along the x(z)-direction. The excited state
is a symmetric superposition of the momentum states
. The energy difference between the two modes is denoted by ω0, which is given by
in the limit of vanishing TP lattice depths. More details can be found in the supplementary information, which is available online at stacks.iop.org/NJP/22/033020/mmedia. Here, the recoil energy is defined as
, where
is the reduced Planck's constant,
is the wavenumber of the TP, and m refers to the atomic mass of 87Rb. By considering the zero-momentum state
and the excited momentum state
, the atomic system is cast into a pseudo spin-1/2 description with corresponding angular momentum operators
and
=
, where cj (
) is the annihilation (creation) operator in mode j. The pseudo spin has length
, with N being the number of atoms. The last term in the Hamiltonian describes the coherent coupling between the two motional states
and
induced by the TP. The microscopic processes related to this coupling can be identified: the atoms can scatter a photon from the TP into the cavity and vice versa, with a net change in momentum by
in both x- and z-direction. The coupling strength λ is proportional to the Rabi frequency of the Raman transition driven by the TP and the cavity (see supplementary information).
For small TP lattice depths (and therefore small coupling strengths), the ground state of the system is the BEC and no photons populate the cavity. For a critical coupling strength of [13, 24, 25]

the system undergoes a phase transition to a self-organized state, where the mode is macroscopically populated. As a spatial counterpart, a density modulation of the atomic cloud in the form of a chequerboard pattern emerges in the self-organized phase. Photons from the TP are coherently scattered off this density modulation into the cavity mode via Bragg scattering. In this way, a coherent intra-cavity light field is built up. The finite cavity decay rate amounts to
in the experiments described in this work.
The finite critical pump lattice depth required for this phase transition results from the competition between kinetic and potential energy. The kinetic energy cost, related to ω0, increases as the density modulation of the gas becomes more pronounced. At the same time, the potential energy is lowered due to the light shift experienced by the atoms in the interference lattice between TP and the red-detuned intra-cavity light field. Smaller cavity detunings Δc favor the build-up of a stronger intra-cavity light field, with the consequence of lowering the critical coupling.
The intra-cavity light field continuously leaks through the cavity mirrors due to the finite cavity decay rate. We detect the photons leaking out of the cavity with a heterodyne detector, which enables the frequency-resolved reconstruction of both electric field amplitude and phase within its detection bandwidth.
We can infer the phase transition from the normal to the self-organized state from two signatures [7]. The first signature is observed in absorption images of the atoms, which are taken after suddenly switching off all light fields and subsequent ballistic expansion of the atomic cloud. For a sufficiently long expansion time, the atomic image reveals the momentum distribution of the atomic cloud. The occurence of momentum peaks corresponding to a chequerboard modulation of the atomic density signals that the atoms were self-organized. This measurement is destructive. The second signature is the onset of a coherent intra-cavity light field, which is continuously recorded via the heterodyne detector. An exemplary trace of the mean intra-cavity photon number with no feedback employed is shown in gray in figure 3(a). The corresponding ramp of the transverse pump power is shown in figure 3(b).
The self-organization phase transition has previously been observed not only in a BEC [7], but also in thermal atomic clouds coupled to an optical cavity [26, 27]. In the latter case, the two-mode description introduced above is not valid, and the system cannot be mapped to the Dicke model. For self-organization to occur, the coupling strength needs to overcome the thermal energy.
When the TP lattice depth is ramped up to a constant value above the critical lattice depth, the mean intra-cavity photon number decreases over time. This decline is caused by the presence of the light fields which leads to heating of the atomic cloud through spontaneous emission and subsequent atom loss. A decrease in atom number results in a lower collective atom-field coupling strength and leads therefore to less intra-cavity photons for a constant TP lattice depth [7, see supplementary information]. This effect is especially pronounced for high mean intra-cavity photon numbers, where the system is deep in the self-organized phase. Furthermore, when the system just crossed the phase transition point, the atomic modulation is weak and the mean intra-cavity photon number low. As the system behaves nonlinearly close to criticality [13], it is highly susceptible to any drifts. This makes it very challenging to study the system at low intra-cavity photon numbers in the vicinity of the self-organization phase transition.
To overcome these challenges, we designed a microcontroller based feedback architecture which stabilizes the mean intra-cavity photon number. The idea is schematically shown in figure 1. We make use of the non-destructive and continuous readout of the system's state provided by the cavity. The signal processing unit calculates the current mean photon number from the signal recorded by the heterodyne detector. It determines the deviation from the desired mean photon number and calculates the feedback signal required to stabilize the system. The feedback signal is sent to the intensity control of the TP laser field.
3. Heterodyne detection and feedback circuit
The intra-cavity light field is detected in a balanced heterodyne measurement [28]. Figure 2 depicts the details of the optical heterodyne measurement, the signal processing stages and the feedback circuit.
The frequency of the intra-cavity light field in the self-organized phase equals the frequency of the TP because photons from the TP are scattered off a static atomic Bragg grating into the cavity mode. This Bragg grating is created in a self-consistent way as a result of the chequerboard potential provided by the interfering light fields. The position of the atomic pattern relatively to the TP lattice determines the phase of the cavity light field. The phase therefore stays constant during the lifetime of the self-organized state [28]. The light field leaking through one of the cavity mirrors is interfered with a local oscillator laser beam (LO) on a beam splitter and sent to two balanced photodiodes. The frequency of the LO is offset by from the TP frequency. The difference signal of the two photodiodes is split, partly phase-shifted and subsequently down-mixed in an analog mixer with a frequency of
. These two quadratures are further digitally down-mixed with a variable frequency δωD in the digital controller. The resulting quadratures are referred to as I(t) and Q(t). We use
to obtain information about the light field stemming from the self-organization process. We calculate the time-dependent mean intra-cavity photon number nph(t) of the intra-cavity light field from the quadratures I(t) and Q(t), see supplementary information. The error signal for the feedback algorithm is calculated as the difference between the set mean photon number n0(t) and the measured mean photon number nph(t):

The feedback signal c(t) is determined through the formula

where Δt is the temporal distance to the previously sent feedback signal and p is an adjustable gain.
Figure 2. Implementation of the feedback scheme. The intra-cavity light field leaking through the mirror is interfered with a local oscillator and detected by a balanced heterodyne detector (green box). The quadratures are first down-mixed in an analog circuit and then further processed in a digital controller to obtain the time-dependent quadratures I(t), Q(t) and the mean intra-cavity photon number nph(t) of the intra-cavity light field. The phase φ(t) of the intra-cavity light field can also be accessed. The determined photon number is compared with the set value n0(t) of the mean intra-cavity photon number (orange box). The calculated feedback signal c(t) is sent to the intensity regulation of an acousto-optic modulator (AOM) in the TP path, changing the value of the pump lattice depth VTP(t) (blue box). Frequencies are adjustable and set to ,
and
for the measurements presented.
Download figure:
Standard image High-resolution imageThe power of the TP is regulated via the transmission through an acousto-optic modulator controlled by analog electronics. This change in intensity is equivalent to a change in the TP lattice depth VTP(t). The setpoint for the TP intensity regulation is either steered externally or dynamically adjusted to the calculated feedback signal c(t). The bandwidth of the feedback circuit is . Technical details on the microcontroller and the implemented algorithm can be found in the supplementary material.
4. Measurements
We prepare a nearly pure BEC of 87Rb with N = 63(5) × 103 atoms in the Zeeman state , where F and mF refer to the total angular momentum manifold and the associated magnetic quantum number. The BEC is confined in a crossed-beam dipole trap with wavelength
centered inside a high-finesse optical cavity. The trap frequencies amount to
,
and
. The wavelength of the TP is set to
. The frequency of the TP is red-detuned by
=
from the cavity resonance frequency. This detuning is kept constant in the experiments presented in this work.
4.1. Comparison of self-organization with and without feedback
We compare measurements with and without the active feedback scheme in figure 3. For the measurement without feedback, the ramp of the TP lattice depth follows a predefined protocol: the lattice depth is first increased via an s-shaped ramp (see supplementary information) within 20 ms to its final value of 5.09(6)Erec. The lattice depth is subsequently held at a constant value for 1 s and finally decreased via an s-shaped ramp within 20 ms. The TP lattice depth for this protocol is depicted as the gray curve in figure 3(b). The mean intra-cavity photon number is shown in figure 3(a). It initially reaches nph = 17.76(8) (calculated using a moving average window of 5 ms), but decays to half its value within 160 ms. There is no signal from self-organization left after about 600 ms.
Figure 3. Comparison of self-organization with (red) and without (gray) feedback applied. (a) Mean intra-cavity photon number nph(t) recorded with the heterodyne detector for both scenarios. (b) Lattice depth VTP(t) of the transverse pump for the stabilized and non-stabilized case. Data is displayed using a moving average window of 5 ms. (c) Absorption images of the atomic cloud after a sudden switch off of the TP power for different hold times when the system is subject to feedback. The images are saturated to an optical density of OD = 0.8 to allow for better visibility of the momentum peaks. The dashed lines in (a) and (b) indicate the hold times for the images shown in (c).
Download figure:
Standard image High-resolution imageWe contrast these results with the measurements with feedback. The TP lattice depth is again increased via an s-shaped ramp within 20 ms to 4.45(5)Erec. The feedback scheme however takes over control as soon as half of the set mean photon number is reached and then steers the ramp of the TP lattice depth. The mean photon number is stabilized to nph = 12.2(3) for more than 1.8 s. The sudden decrease in the mean photon number results from switching off the TP power at a maximum lattice depth of 44.5(5)Erec. This value is the technical upper limit of the TP lattice depth in our feedback scheme due to the dynamic range of the microcontroller. With feedback, the TP lattice depth is increased in a nonlinear fashion.
In addition, we show in figure 3(c) absorption images of the atomic cloud for the case with feedback. The images were taken after suddenly switching off the TP after variable hold times. The hold time is defined as the time we let the system evolve after the initial ramp-up stage of 20 ms. We obtain the momentum distribution via the absorption images after 8 ms of ballistic expansion [8].
For short hold times, sharp momentum peaks are visible in the absorption images. We can understand these momentum peaks by invoking the two-mode description of the atomic cloud in momentum space as presented in section 2. The momentum peaks at px = 0, pz = 0 and at signal the occupation of the zero momentum mode
and of the excited momentum mode
. The momentum peaks at px = 0,
stem from the mere presence of the TP lattice.
The visibility of the higher momentum peaks decreases for increasing hold times, signaling a loss of atomic coherence (see supplementary information). Simultaneously, the thermal fraction increases. This behavior is due to the TP heating the atomic cloud, leading to the population of higher momentum states not captured by the two-mode description. Eventually we observe self-organization of a density-modulated but completely thermal atomic cloud [26]. The exact theoretical description of the evolution from self-organization in a two-mode system to self-organization in a multi-mode system is the subject of ongoing theoretical efforts [29, 30] and beyond the scope of this work. The measured intra-cavity light field did not exhibit any visible change within our detection sensitivity during the stabilization (see supplementary information).
As described above, the TP lattice depth for the case with feedback is increasing in a nonlinear and convex fashion. Qualitatively we comprehend the curve's shape by considering atom loss. The atom loss is compensated with a higher TP lattice depth, which leads to further heating and atom loss, requiring an even higher TP lattice depth. We performed measurements of the atom number for the case with feedback after adiabatically ramping down the TP lattice depth for different hold times up to 1 s, as presented in supplementary information. The atom number decay exhibits a linear trend. If we assume such a linear atom number decay during the complete time of the experiment and apply a relation between TP lattice depth and atom number valid in the two-mode description of the atomic cloud (see supplementary information), we obtain a similarly convex shaped curve as in figure 3(b).
4.2. Stabilization on different mean intra-cavity photon numbers
In an additional set of experiments, we explore the range of mean intra-cavity photon numbers the feedback scheme can stabilize. The data in figures 4(a) and (c) demonstrate that stabilization is possible in a wide range from nph = 0.17(4) to nph = 27.6(5). Due to the increased nonlinearity close to the phase transition, the stabilization at low photon numbers nph < 1 requires adjustments: the feedback algorithm takes over control only after the set photon number n0 is reached, and the gain settings both in the feedback software and TP intensity regulation need to be adjusted.
Figure 4. The feedback scheme can stabilize the mean intra-cavity photon number nph in a wide range. This implies that we can stabilize the system both close to the phase transition and deep in the self-organized phase. (a) Constant mean intra-cavity photon number traces of nph = 27.6(5), 20.6(4), 13.7(3), 6.8(2), 3.28(18), and (zoom in (c)) nph = 0.58(7), 0.30(6), 0.17(4). (b) Corresponding TP lattice depths, as set by the feedback circuit. Data is displayed using a moving average window of 5 ms.
Download figure:
Standard image High-resolution imageWe quantify where the stabilization is engaged relative to the phase transition. For the highest photon number nph = 27.6(5), the stabilization starts for being 50.7(20)% larger than the critical lattice depth of the self-organization phase transition. Here,
is the value of the TP lattice depth for which the desired photon number n0 is reached for the first time. The critical lattice depth is extracted from the experimental data (see supplementary information). For the lowest achieved photon number nph = 0.17(4), we start the stabilization at a
which is only 1.6(10)% above the critical TP lattice depth, well within the critical regime of this phase transition [20, 21]. The respective values for all stabilized nph can be found in the supplementary information.
5. Conclusion and outlook
We developed a microprocessor-based active feedback scheme which stabilizes the number of mean intra-cavity photons during self-organization of a BEC in an optical cavity. The feedback is acting on the intensity of the TP laser beam. The feedback scheme can easily be adjusted to act on other parameters, as for example the frequency detuning between TP and cavity. As the balanced heterodyne detection system records the electric field within the bandwidth of the detection system, the feedback software can be readily modified to calculate other physical quantities, like the phase or information in the spectrum of the light, and control the evolution of these observables. The detection of the cavity field and subsequent regulation can also be extended to the polarization degree of freedom.
The feedback scheme expands the possibilities for experimental studies of self-organization in hybrid atom-cavity systems. The current capabilities of the feedback scheme enable us to approach the phase transition starting from the self-organized phase in a highly controlled manner. Stabilizing the system close to the phase transition allows to make use of its increased sensitivity with respect to external perturbations for sensing applications [31–33]. Another interesting prospect is to study fluctuations over long times in the direct vicinity of the phase transition [21].
Instead of stabilizing on a constant mean intra-cavity photon number, the feedback scheme can be modified to modulate the atom-light coupling strength according to more complex schemes. Such control is important to engineer non-equilibrium phases and phase transitions [22, 23, 34, 35].
Acknowledgments
We thank Tobias Delbrück and Chang Gao from the Sensors Group in the Institute of Neuroinformatics, University of Zurich and ETH Zurich, for stimulating discussions and helpful advice. We thank Alexander Frank for electronic support and Fabian Finger for careful reading of the manuscript. We acknowledge funding from SNF, project numbers 182650 and 175329 (NAQUAS QuantERA), and NCCR QSIT, and funding from EU Horizon2020, ERCadvanced grant TransQ (Project Number 742579).