Femtojoule-scale all-optical latching and modulation via cavity nonlinear optics

Femtojoule-scale all-optical latching and modulation via cavity nonlinear optics

Yeong-Dae Kwon    Michael A. Armen    Hideo Mabuchi hmabuchi@stanford.edu Edward L. Ginzton Laboratory, Stanford University, Stanford, CA 94305
July 19, 2019
Abstract

We experimentally characterize Hopf bifurcation phenomena at femtojoule energy scales in a multi-atom cavity quantum electrodynamical (cavity QED) system, and demonstrate how such behaviors can be exploited in the design of all-optical memory and modulation devices. The data are analyzed using a semiclassical model that explicitly treats heterogeneous coupling of atoms to the cavity mode. Our results highlight the interest of cavity QED systems for ultra-low power photonic signal processing as well as for fundamental studies of mesoscopic nonlinear dynamics.

pacs:
37.30.+i, 42.50.Nn, 42.65.Pc, 42.65.Sf

The diverse phenomena of cavity nonlinear optics Mandel (2005) provide a rich basis for fundamental studies of dissipative nonlinear dynamics Lugiato and Wolf (1984); Otsuka (2000) and for the design of photonic signal processing devices Gibbs (1985). Recent experiments exploiting resonant atomic nonlinearities Rempe et al. (1991); Mabuchi et al. (1996); Sauer et al. (2004); Gupta et al. (2007); Mücke et al. (2010); Kerckhoff et al. (2011) and/or nanophotonic cavities Alton et al. (2011); Srinivasan and Painter (2007); Faraon et al. (2008); Nozaki et al. (2012); Groblacher et al. (2009); Weis et al. (2010); Venkataraman et al. (2013) have demonstrated that such phenomena can occur at very low (femtojoule–attojoule) energy scales in prototypical systems, reaching down to the quantum-physical few photon regime Carmichael (2007); Armen and Mabuchi (2006) and raising intriguing prospects for corresponding ultra-low power photonic information technology Vlasov (2012); Santori et al. (2010); Miller (2010). One of the less-explored dimensions of this developing scenario is the surprising complexity of dynamical behaviors, beyond simple thresholding and bistability, that can be achieved at low energy scales in cavities incorporating two-level atoms or comparable solid-state emitters. In this Letter we present experimental data on sub-critical Hopf bifurcation phenomena in an atom-cavity system and demonstrate a novel device configuration that latches a logic level and provides rf modulation of the output signal at the same time.

The general paradigm of two-level emitters coupled to an isolated optical mode can be realized in a wide range of gas, liquid, and solid-state systems with various types of optical cavities; our experiment uses a 12 mm Fabry-Perot resonator in ultrahigh vacuum with an ensemble of laser-cooled Cs atoms placed between the mirrors. We utilize only the F=m=4 to F=m=5 transition of the D2 line at 852 nm where the atoms are well modeled as two-level systems. The following semiclassical equations of motion is used to model the system composed of a cavity and N weakly-coupled, noninteracting atoms Drummond (1981):

(1a)
(1b)
(1c)

Here, is the intracavity field annihilation operator, are the rasing/lowering operators for the -th atom, and . The rates are the atom-cavity coupling strength for the -th atom, the atomic polarization decay rate, and the cavity field decay rate, respectively. and are the cavity and the atom detunings relative to the drive (laser) frequency, defined as and , and is the drive field amplitude.

The above Maxwell-Bloch-type equations are derived from an unconditional Jaynes-Cummings master equation in the same manner as in Armen and Mabuchi (2006) (see also Mabuchi (2008)), but with a slight generalization that incorporates distinct coupling strengths for the individual atoms. We note that the ‘factorization’ approximations made in the derivation, such as , mean that entanglement between the atoms and the field is ignored. This is justified because the coupling strengths for individual atoms are small, while the combined coupling strength as an ensemble is still large enough to cause a strong nonlinearity.

By setting the time derivatives in Eqs. (1) to zero one can solve for the steady state values of , , and . The independent variable is then expressed as a function of to derive the input-output characteristic of the system, and the stability of any (hyperbolic) equilibrium point can be determined by local linearization. Alternatively, Eqs. (1) can be numerically integrated to produce sample trajectories that reflect nonlinear dynamical phenomena such as limit cycles, chaos (Gang et al. (1990); Yang et al. (2005)), etc. In this Letter we make use of both equilibrium and transient data to assess the accuracy of our theoretical modeling approach.

In our experiment, the atoms are first collected and cooled by a magneto-optical trap (MOT) formed on a gold-coated mirror surface Reichel et al. (1999). This surface MOT is located outside the cavity, 2 cm away from the cavity axis, where the cooling beams are not obstructed by the cavity mirrors. The atoms are then further cooled to sub-Doppler temperature (K) by polarization-gradient cooling, optically pumped to the m=4 state, and then loaded in the magnetic trap. The magnetic trap, formed by currents through copper wires (m diameter) buried underneath the mirror surface, can be moved by shifting the currents Hänsel et al. (2001), and is used as a conveyer to transport atoms (up to ) into the cavity mode. The number of atoms to be transported are adjusted by changing the initial MOT size. At the end of the transport process ( ms), the atoms are released from the magnetic trap and a large uniform magnetic field is applied in the direction of the cavity axis, which ensures that the quantization axes for the atoms and the field are aligned in the presence of possible stray fields near the surface. The strong bias field also lifts the degeneracy among the Zeeman states (Zeeman shift: MHz). Once the atoms are in position, the input laser is injected from one side of the Fabry-Perot cavity and drives the TEM mode. The beam that leaks out from the other side of the cavity serves as the output of the system. Due to the loss at the cavity mirrors and the existence of two possible directions for the leakage, we estimate , the probability of an intracavity photon making it to the output channel, to be 0.41. The output beam goes through an optical isolator and is then measured by a homodyne/heterodyne detection scheme Kerckhoff et al. (2011). The maximum atom-cavity coupling constant, the atomic polarization decay rate, and the measured cavity field decay rate are = (0.72, 2.6, 2.3) MHz, respectively.

When we apply the model to our experiment, we assume a large number of atoms distributed uniformly inside the cavity, so that the values of in Eqs. (1) follow directly from the gaussian standing-wave profile of the resonant field mode. Under this assumption the atomic ensemble is characterized by a single parameter , where is the maximum coupling constant at the anti-node center of the cavity. This simplification applies best for a cavity filled with a homogeneous dispersive/absorptive medium, a scenario for which extensive studies have already been performed using the Maxwell-Bloch Equations (MBEs) Rosenberger et al. (1991); Lugiato and Wolf (1984). The method used here differs from the conventional analytical treatment of the MBEs in that we do not simplify Eqs. (1a)-(1c) using the approximation of collective atomic variables, which assumes, among other things, low absorption rate by the medium Lugiato and Wolf (1984). We instead keep track of the evolution of individual atomic variables for better accuracy, at the cost of increased computational load.

The accuracy of our model in a regime of relatively simple dynamic behavior is illustrated in Figs. 1(a) and 1(b), which display the measured equilibrium input-output response of our atom-cavity system for two different parameter sets together with corresponding steady state solutions of Eqs. (1). For the parameters of Fig. 1(a) (, ) a soft thresholding behavior is observed with a corner near 6 nW input power. With increased atom number and a detuned probe (, ) as in Fig. 1(b) we observe classic bistability and hysteresis with coexisting stable equilibria from roughly 14 to 15.5 nW input power. For these measurements the input power was swept slowly (20 kHz) compared to the dynamical rates of the system () to probe the adiabatic response. Figures 1(a) and 1(b) each show an individual (but precisely repeatable) data trace of s duration that was recorded continuously as the input power was swept up and down.

Figure 1: (a) Threshold behavior of the system with , . The black line is the stable steady state solution derived from the model, and red (green) line is a typical data trace from the measurement when the input power is gradually increased (decreased). (b) Hysteretic behavior with , . Color scheme is the same as (a) except for the new black dotted line that indicates the unstable steady state solution. (c) Experimental demonstration of latching using bistable system parameters of (b). With the input power constant at 14.5nW the system maintains its current state. When a pulse is applied the system state is set high or low according to the pulse type (positive/negative).

The theory curves in Figs. 1(a) and 1(b) were computed using manually-adjusted values of and to improve agreement with our data, but the fitted values differ from the independently measured ones by no more than 10%. is recorded for each experimental run by measuring the spectral response of the atom-cavity system using a very low intensity probe Thompson et al. (1992) right before and after the main measurement is made. Comparison of these two spectral curves also ensures that does not vary significantly during the measurement period. This procedure and statement of calibration accuracy apply to all remaining data-theory comparison plots in this Letter as well. Some discrepancy between the model and the measurement is expected, mainly because the magnetic trap from which the atoms are released just before the measurement has a similar size to the cavity mode (waist m), and therefore it is difficult to ensure that the atoms are truly uniformly distributed when the measurement is taken. Up to these inaccuracies however our model has proven to be a reliable predictor of the parameter values at which interesting nonlinear dynamical phenomena should be observed.

After confirming the parameters for simple bistable input-output behavior in our atom-cavity system, we performed an initial proof-of-principle experiment to demonstrate an all-optical set-reset (SR) latch, which could for example be used as an optical memory bit. Typical data is shown in Fig. 1(c). Here we envision a device configuration in which the set/reset control beam shares an optical input channel with the bias power beam (for example they could be combined by a beam splitter) so that they interfere constructively or destructively depending on the phase of the control beam. With this picture in mind, we demonstrated the set/reset control by superposing positive and negative pulses onto the laser beam that drives the atom-cavity system (upper blue trace of Fig. 1(c)). A positive pulse switches the response to the upper branch of the hysteresis curve, where it latches until a negative pulse resets it back to the lower branch. Our system required 15 nW optical bias power and the energy used for each control pulse was roughly .

We note that the width of the control pulses used in our experiment corresponds to the inverse cavity decay rate . Our switching time is thus longer than in nanoresonator-based devices, which generally have much larger due to their small size. However, in terms of the switching energy, ours appears to be the minimum energy demonstrated in an all-optical bistable device Nozaki et al. (2012). The switching energy potentially could be lowered even more as the energy scale in this experiment is still much larger than the optical shot noise (quantum fluctuation) limit Kerckhoff et al. (2011).

Figure 2: Experimental parameters , , are used. (a) Homodyne signal generated by the output beam as the input power is slowly swept up and down over the Hopf bifurcation point . The gray bars labeled with ‘b’ and ‘c’ indicate the regions of which the zoomed-in views are displayed in (b) and (c). (d) Theoretical prediction of the input-output characteristics of the system. The solid and the dotted black lines represent the stable and the unstable steady state solutions, respectively. Points and indicate Hopf bifurcations discussed in the text. The upper and lower blue lines trace the maximum and the minimum output power of the limit cycles observed in simulations. (e) On the bottom, the predicted output power of the system (d) near the subcritical bifurcation point is compared with the measurement. The measured data is plotted in red (up-sweep) and green (down-sweep), using the following rules: When no oscillation is detected, the mean output power is recorded by a cross-mark. When the oscillation is detected, a sinusoid is fitted to the measured signal and its max/min optical power is denoted by a pair of connected circles. On the top, we plot the frequencies of the both simulated and measured limit cycles. The simulation result is plotted in brown, while the measurements are denoted by red (up-sweep) and green (down-sweep) square marks.

With an appropriate choice of external parameters , and guided by theory, we also probed the atom-cavity dynamics in an unstable regime and observed “self-pulsing” behavior  Lugiato and Wolf (1984); Lugiato et al. (1984) in which the global attractor is a limit cycle. In Fig. 2(a) we show the system making transitions between a stable equilibrium point and a limit cycle as we sweep the driving input power. When the input power is large, the power of the output beam oscillates significantly at a frequency in the range of  MHz.

Our model predicts a range of input powers ( nW) for which a stable equilibrium point and a stable limit cycle coexist. Starting from low input, as the system is driven past the sub-critical Hopf bifurcation point (1200 nW, shown in Figs. 2(d) and 2(e) the equilibrium point becomes unstable and the system jumps to the the limit cycle. With the system oscillating, the input power can be lowered back down and the limit cycle remains until another bifurcation occurs at a lower input power (1150 nW) than . As shown in Fig. 2(e), such hysteresis can be observed clearly in our experiment, and we find that the measured amplitude and frequency of the limit cycle closely match predictions. Moreover, combined measurements of the amplitude and phase quadratures of the output field confirm that the oscillation of the intracavity field follows a trajectory on the optical phase plane that is well predicted by theory (See Supplementary Fig. S1).

Near the predicted super-critical Hopf bifurcation point ( 2100nW), on the other hand, the match between our measurements and theory was not as good. Simulation predicts that the amplitude of the limit cycle should gradually converge to zero as the system approaches from below, at which point a stable equilibrium point defines the steady-state. However, the experimental system continued to show large-amplitude oscillatory behavior at random intervals in time even when the input power was far above the predicted location of . Simulations with added technical noise suggest that such behavior is quite reasonable, however, as the equilibrium points near are highly under-damped.

Figure 3: The system shown in Fig. 2 is used to demonstrate latching between a stable equilibrium point and a limit cycle. The spectrogram in the middle is generated using s segments of the homodyne data, which contain roughly 5 oscillation cycles. 0 dB power level was arbitrary chosen. Note that more bandwidth is retained in the homodyne signal displayed here than in Fig. 1(c), in order to preserve the limit cycle oscillations. In the bottom frame we show a zoomed-in view of the homodyne signal from s to s.

The limit cycle behavior of our system converts the constant input power into an rf-modulated output power. We are thus able to demonstrate an elementary physical mechanism that potentially could be adopted to make self-contained inline optical oscillators in a nanophotonic setting, for which the use of conventional optical-oscillator-type devices such as Q-switched lasers might not be practical. Our system furthermore is capable of transferring baseband power modulation to rf modulation of the output via the correlation between the input power and the frequency of the limit cycle (Fig. 2(e), upper plot). This suggests the possibility of an all-optical baseband-to-rf up-converter. And finally, hysteresis around the sub-critical Hopf bifurcation can be used to make a latch of rather unique characteristics, which switches between dc and modulated signal formats. Optical dc-rf latching devices may find use in optical communication systems where signal up-conversion is a common practice to avoid the high level of the background noise at low frequencies. Fig. 3 shows a demonstration of such latching performance. In this context we note that solid-state implementations of this type of cavity nonlinear optical device, for example using ensembles of quantum dots Hennessy et al. (2007) or vacancy centers Barclay et al. (2009), could achieve much higher values of and therefore much higher limit cycle frequencies, potentially in the GHz range.

We find that the energy stored in the optical dc-rf latch in Fig. 3 is comparable to that of the elementary latch in Fig. 1, as in both cases the cavity contains about 500-1000 photons (100-200 attojoule) on average in the high energy state and much fewer when in the low energy state. However, the energy required for switching between the two states was more than an order of magnitude higher for the dc-rf latch (30 fJ). This is mostly due to the mismatch between the input laser frequency and the resonances of the system. In Supplementary Section Fig. S2 we show a demonstration of sub-femtojoule dc-to-rf switching operation of the same Hopf bifurcation latch using a secondary near-resonant beam. We also note that the use of pulse shaping on the control pulse Sandhu et al. (2010) could potentially reduce the switching energy of both the elementary and the dc-rf latch further.

In conclusion, we have demonstrated that a cavity-based optical device that incorporates a resonant nonlinear response can exhibit complex dynamic behaviors such as hysteretic limit cycle formation at femtojoule energy scales. We have suggested novel signal processing functions based on such phenomena, and have demonstrated an all-optical configuration that simultaneously latches a logic level and performs rf modulation of the output power. Intrinsic self-oscillatory behavior of the atom-cavity system, once a subject of active research using thermal gases and mW of laser power Orozco et al. (1984, 1989); Segard and Macke (1988), here has been probed with a laser-cooled atoms and W power levels. Analogous phenomena could be realized in the context of nanophotonics, where there is already an emerging interest in observing and functionalizing self-oscillations Grigoriev and Biancalana (2011); Malaguti et al. (2011); Gu et al. (2012).

This work has been supported by the ARO (W911NF-08-1-0427) and DARPA (N66001-11-1-4106).

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
316616
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description