We demonstrate a novel all-optical quantum random number generator (RNG) based on above-threshold binary phase state selection in a degenerate optical parametric oscillator (OPO). Photodetection is not a part of the random process, and no post processing is required for the generated bit sequence. We show that the outcome is statistically random with 99% confidence, and verify that the randomness is due to the phase of initiating photons generated through spontaneous parametric down conversion of the pump, with negligible contribution of classical noise sources. With the use of micro- and nanoscale OPO resonators, this technique offers a promise for simple, robust, and high-speed on-chip all-optical quantum RNGs.
© 2012 OSA
True RNGs are desirable for applications ranging from cryptography to computer simulations. Quantum phenomena prove to be attractive for physical RNGs due to their fundamental randomness and immunity to attack [1–8]. It was recently demonstrated that quantum correlation of non-local states in two entangled atomic systems can violate Bell inequality principle and be used to generate random numbers . Entangled-photon and single-photon sources have been used for RNGs . Phase of vacuum fluctuations measured by a homodyne detection scheme , photon arrival time , and randomness in photon number distribution of a coherent source  have also been exploited as the engines of RNGs. Some other quantum RNGs are based on phase measurements of laser diode outputs [6, 7] and the Stokes light in Raman scattering . Most of these methods require massive post-processing to ensure randomness, in some, photodetection is an essential part of the random process, and several require a complicated experimental environment.
Optical parametric down conversion is an essential element in most quantum optical experiments including optical squeezing , and generation of entangled photons . In an OPO, photons generated through spontaneous down conversion of the pump initiate the oscillation in the absence of other inputs [11, 12]. This quantum process is the dominant effect during the oscillation build-up, leading to selection of one of the two possible phase states above threshold in a degenerate OPO . Building on this, we implement a synchronously pumped twin degenerate OPO, which comprises two identical independent OPOs in a single cavity, and measure the relative phase states of the OPO outputs above threshold as a bit value.
The quantum RNG presented in this paper does not require sophisticated detection techniques because the OPO operates above threshold. Our system is further simplified by using commercially available fiber lasers and well-developed quasi phase matching technology for nonlinear crystals . The optical beam at the output randomly toggles between two intensity levels with a well-defined clock signal. The randomness is intrinsically inherited from the quantum noise of the system, and we verify minimal contribution of classical noise sources by examining the dependence of randomness on the phase of initiating photons in the OPO. A sequence of 1 billion bits is generated at 10 kbps which has successfully passed all the statistical randomness tests provided by National Institute of Standards and Technology (NIST).
The operation principle is based on the phase response of a degenerate OPO. In a continuous wave (c.w.) type I degenerate OPO, where signal and idler are indistinguishable (we call it signal in this paper), stable oscillation above threshold can occur in one of the two possible phase states . This is because the gain mechanism in such an OPO, provided by degenerate parametric amplification, is phase sensitive and has a period of π for the signal phase .
Figure 1(a) schematically illustrates the incremental amplitude gain in such an amplifier as a function of the relative phase between the pump and the signal (θs). Maximum amplification occurs both at θs = 0 and θs = π, where the energy flows from the pump into the signal, and maximum deamplification at θs = π/2 where the energy flows from the signal to the pump. Equivalently, this period of π exists in the phase dependence of quantum fluctuations in a squeezed state produced by degenerate parametric down conversion . If the maximum gain exceeds the loss in the resonator (gmax > gth = δE, where δE is the electric field round-trip loss ), depending on the zero-point fluctuations of the signal modes, the OPO will stably oscillate with the phase of either θs = 0 or θs = π. This above-threshold phase state is inherited from the vacuum fluctuations , and the design parameters do not favor oscillation in one or another.
To make a quantum RNG based on this process, the relative phase states between two identical OPOs can be measured by interfering the outputs. This is demonstrated in the schematic in Fig. 1(b), where two degenerate OPOs are pumped with the same source, and an optical switch is used to restart the OPOs periodically. If both path lengths are matched, the interference of the output signals will randomly toggle between high and low intensity levels resulting in a sequence of “zeros” and “ones” at the photodetector output.
We realized the OPO-based quantum RNG using the synchronously pumped femtosecond degenerate OPO demonstrated in . For such a short-pulse OPO, similar to the c.w. case, the signal can take one of the two possible phase states  which destructively interfere with each other. Figure 1(c) depicts an example of the electric field of the pump and the signal in the two possible phase states. The existence of these phase states for a degenerate femtosecond OPO was experimentally verified in .
3. Experimental setup
To avoid implementation of two matched resonators, a twin OPO is used, in which two identical OPOs operate in the same ring resonator with the roundtrip path length twice the separation length of pump pulses, as depicted in Fig. 2(a). Figure 2(b) shows the pump and signal pulse trains for this twin OPO. The pump pulses generate two independent signal pulse trains, marked as OPO1 and OPO2 in the figure. These two temporally separated OPOs have half the repetition rate of the pump, the same polarization and spectral properties, and experience exactly the same optical paths in the oscillator. An unequal arm interferometer is used to measure the relative phase states of OPO1 and OPO2 by interfering their temporally separated pulse trains, as depicted in the schematic of Fig. 2(a).
The OPO is pumped by a 1560-nm mode-locked Er-fiber laser (Menlo Systems C-fiber, 100 MHz, 70 fs, 300 mW) where beam is conditioned by a mode-matching telescope for efficient pumping of the OPO. The resonator is a 6-m ring cavity corresponding to twice the spacing between the pump laser pulses. The cavity optics comprise one pair of concave mirrors with ROC=50 mm (M2 and M3), and six flat mirrors, five of which are gold coated with approximately 99% reflection (M4–M8). A single dielectric mirror (M1) is used to introduce the pump, which has 90% transmission for the pump and more than 99% reflection in the 2.8 – 4 μm. This mirror has a ‘chirped’ design of dielectric layers to compensate the dispersion of the nonlinear crystal.
Broadband gain centered around 3.1 μm is provided by 1-mm long MgO-doped periodically poled lithium niobate (MgO:PPLN) crystal. The poling period is 34.8 μm for broadband type-0 (e=e+e) phase matching at a temperature of 32°C. The crystal is cut such that the mid-IR beam propagates perpendicular to the poling domains when the beam enters at the Brewster angle. The beam waist for the signal in the crystal is ∼10 μm. The mirrors M2 and M3 are set to 5-degree angle of incidence to compensate the astigmatism caused by the Brewster angled crystal and allow stable resonances in the 6-m long cavity. The output is extracted with a pellicle beam splitter (OC) having ∼8% reflection over a broad bandwidth. The filters are AR coated Ge substrates to block the pump and transmit the mid-IR signal.
Oscillation occurs when signal and idler waves are brought into degenerate resonance by fine-tuning the cavity length with the piezo stage of M4. Three of these resonances occur separated by ∼1.5 μm of roundtrip cavity length, corresponding to half of the signal central wavelength . Continuous operation of the OPO is obtained by locking the cavity length to track the center of the strongest resonance using a dither-and-lock scheme . The twin OPO starts oscillating at a pump average power of about 120 mW, and the maximum mid-IR output power is 4 mW.
4. Results and discussions
The OPO stabilized with the servo loop could run continuously for several days. Figure 3(a) shows the signal spectrum centered at 3.1 μm and the pump centered at 1.56 μm. Complementary stable fringe patterns at the output of the interferometer were obtained as depicted in Fig. 3(b) when the beams in the arms are slightly angled vertically. Blocking and unblocking the pump resulted in random toggling between these two patterns (Media 1).
To capture a bit stream, the beam angles are well-aligned in the interferometer and a photodetector is used at the output, while an Acousto-Optic Modulator (AOM) causes periodic restarting of the twin OPO. Figure 3(c) depicts a sample of the signal at the output of the interferometer along with the clock signal at 10 Kbps applied to the AOM. The binary sequence is extracted from the interferometer output as depicted in the plot.
A sequence of 1 billion bits is taken with this method, and it proves to be random with an average of 0.5000. To verify the statistical randomness, a series of tests developed by NIST  are performed, and the summary of results are presented in Fig. 4. The 1-Gb sequence passed all the NIST statistical tests indicating it is random with 99% confidence.
The maximum bit-rate supported by this RNG depends on the turn-“on” and turn-“off” dynamics of the OPO. At the end of each clock cycle, the intracavity field must decay to the quantum noise level , or the residual field from the previous state will seed oscillation of the next state, and randomness will be lost. In the design presented here, the clock speed is slow enough (and OPO “on”-time long enough) that oscillation builds to a steady-state level. In this case, the clock rate is limited by the time it takes for the intensity to decay below noise, which is 10 to 20 times longer than the 1/e cavity decay time when the OPO is pumped well over threshold and allowed to reach steady state. Therefore, the RNG speed is limited by the cavity decay time, and faster bit rates, in the Gbps range are expected to be achievable using pumps with higher repetition rates, and shorter OPO cavities [20, 21].
However, one can operate the OPO closer to threshold, where the build-up time is longer than the decay time. In this case, the oscillation will not necessarily reach steady-state, but relative phase may still be measured with a sufficiently sensitive detector to establish the bit value, decoupling maximum RNG speed from the cavity decay time. A potential extreme case is eliminating the cavity and having a single pass parametric down conversion with a speed as high as the repetition rate but requiring either a sensitive detection system or a relatively high peak power. It is even possible to design a similar RNG with a c.w. pump to eliminate the repetition rate limit.
We exploit the OPO dynamics discussed above to verify that the randomness is due to the phase of photons initiating the oscillation. As mentioned before, our RNG operates in the regime where build-up is much faster than decay. This is due to the low (∼23%) modulation depth of the AOM that biases the pump only slightly below threshold during the “off” phase of the clock, resulting in a much longer decay time than if the OPO were un-pumped. This makes it easy to observe the transition to a fixed-state, ‘latched’ oscillation, when at faster clock rates the turn “on” edge arrives while many photons remain in the cavity. Figure 5(a) shows how the randomness – measured here by the rate of bit flip in a 100-kb long sequence – breaks as a function of clock speed. Figure 5(b), 5(c), 5(d) illustrate output samples for three points on this curve. This experiment shows that when an intracavity phase reference exists – in the form of residual photons from the previous clock cycle – the randomness of the sequence breaks; the influence of classical noise sources is negligible since they cannot change the phase state.
4.1. Estimated maximum speed
For quantum random outcome, the turn off time should be long enough to allow the intracavity power to decay from the steady state level, that is about 1 W, to the quantum noise level, i.e. one photon per mode (Pnoise = hνΔν ), which is about 1 μW. Here hν is the photon energy at the central signal wavelength of 3.1 μm, and Δν is the OPO bandwidth at 3-dB level, estimated to be ∼10 THz (Fig. 3(a)). The intensity decay time of the OPO can be estimated using:
We demonstrate a novel technique for all-optical quantum RNGs based on the randomness in spontaneous parametric down conversion by measuring the above-threshold phase of oscillation in a degenerate OPO. The twin system, made of two identical OPOs sharing the same cavity, facilitates the phase state measurement. Stable operation is achieved using a free running commercially available mode-locked fiber laser. Unlike many other quantum RNGs, this system enables quantum random number generation with no need for electronic or computer post processing on the generated bit sequence, and has the potential of extremely fast operation. Implementation of an OPO-based quantum RNGs is not limited to second-order nonlinear materials and free space elements as presented here; one can extend this approach to on-chip χ3 OPOs [22–24] where CMOS compatibility can open new possibilities for new generation of RNGs for future communication and computational systems. Combined with micro- and nano resonators, this technique paves the way for high-speed all-optical quantum RNGs with multi Gigabit per second rates.
The authors would like to thank S. E. Harris, Y. Yamamoto, M. M. Fejer, C. R. Phillips, J. S. Pelc, C. W. Rudy, and O. Crisafulli for invaluable discussions, and K. Urbanek for experimental support.
References and links
1. S. Pironio, A. Acin, S. Massar, A. B. de la Giroday, D. N. Matsukevich, P. Maunz, S. Olmschenk, D. Hayes, L. Luo, T. A. Manning, and C. Monroe, “Random numbers certified by Bell’s theorem,” Nature 464, 1021–1024 (2010). [CrossRef]
2. M. Fiorentino, C. Santori, S. M. Spillane, R. G. Beausoleil, and W. J. Munro, “Secure self-calibrating quantum random-bit generator,” Phys. Rev. A 75, 032334 (2007). [CrossRef]
3. C. Gabriel, C. Wittmann, D. Sych, R. Dong, W. Mauerer, U. L. Andersen, C. Marquardt, and G. Leuchs, “A generator for unique quantum random numbers based on vacuum states,” Nat. Photonics 4, 711–715 (2010). [CrossRef]
4. M. Wahl, M. Leifgen, M. Berlin, T. Rhlicke, H. Rahn, and O. Benson, “An ultrafast quantum random number generator with provably bounded output bias based on photon arrival time measurements,” Appl. Phys. Lett. 98, 171105 (2011). [CrossRef]
5. M. Ren, E. Wu, Y. Liang, Y. Jian, G. Wu, and H. Zeng, “Quantum random-number generator based on a photon-number-resolving detector,” Phys. Rev. A 83, 023820 (2011). [CrossRef]
6. M. Jofre, M. Curty, F. Steinlechner, G. Anzolin, J. P. Torres, M. W. Mitchell, and V. Pruneri, “True random numbers from amplified quantum vacuum,” Opt. Express 19, 20665–20672 (2011) [CrossRef]
8. P. J. Bustard, D. Moffatt, R. Lausten, G. Wu, I. A. Walmsley, and B. J. Sussman, “Quantum random bit generation using stimulated Raman scattering,” Opt. Express 19, 25173–25180 (2011). [CrossRef]
10. Y. H. Shih and C. O. Alley, “New Type of Einstein-Podolsky-Rosen-Bohm experiment using pairs of light quanta produced by optical parametric down conversion,” Phys. Rev. Lett. 61, 2921–2924 (1988). [CrossRef]
11. W. H. Louisell, A. Yariv, and A. E. Siegman, “Quantum fluctuations and noise in parametric processes I,” Phys. Rev. 124, 1646–1654 (1961). [CrossRef]
12. S. E. Harris, M. K. Oshman, and R. L. Byer, “Observation of tunable optical parametric fluorescence,” Phys. Rev. Lett. 18, 732–734 (1967). [CrossRef]
13. C. D. Nabors, S. T. Yang, T. Day, and R. L. Byer, “Coherence properties of a doubly-resonant monolithic optical parametric oscillator,” J. Opt. Soc. Am. B 7, 815–820 (1990). [CrossRef]
14. L. E. Myers, R. C. Eckardt, M. M. Fejer, R. L. Byer, W. R. Bosenberg, and J. W. Pierce, “Quasi-phase-matched optical parametric oscillators in bulk periodically poled LiNbO3,” J. Opt. Soc. Am. B 12, 2102–2116 (1995). [CrossRef]
15. A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University Press, 2007).
17. S. T. Wong, K. L. Vodopyanov, and R. L. Byer, “Self-phase-locked divide-by-2 optical parametric oscillator as a broadband frequency comb source,” J. Opt. Soc. Am. B 27, 876–882 (2010). [CrossRef]
18. A. Marandi, N. Leindecker, V. Pervak, R. L. Byer, and K. L. Vodopyanov, “Coherence properties of a broadband femtosecond mid-IR optical parametric oscillator operating at degeneracy,” Opt. Express 20, 7255–7262 (2012). [CrossRef]
19. A. Rukhin, J. Soto, J. Nechvatal, M. Smid, E. Barker, S. Leigh, M. Levenson, M. Van- gel, D. Banks, A. Heckert, J. Dray, and S. Vo, “A statistical test suite for random and pseudorandom number generators for cryptographic applications,” NIST special publication 800-22, Rev. 1-a, NIST, Gaithersburg, Maryland, USA, (2010).
20. S. Lecomte, R. Paschotta, S. Pawlik, B. Schmidt, K. Furusawa, A. Malinowski, D. J. Richardson, and U. Keller, “Synchronously pumped optical parametric oscillator with a repetition rate of 81.8 GHz,” IEEE Photon. Technol. Lett. 17, 483–485 (2005). [CrossRef]
21. J. U. Furst, D. V. Strekalov, D. Elser, A. Aiello, U. L. Andersen, Ch. Marquardt, and G. Leuchs, “Quantum light from a whispering-gallery-mode disk resonator,” Phys. Rev. Lett. 106, 113901 (2011). [CrossRef]
22. J. S. Levy, A. Gondarenko, M. A. Foster, A. C. Turner-Foster, A. L. Gaeta, and M. Lipson, “CMOS-compatible multiple-wavelength oscillator for on-chip optical interconnects,” Nat. Photonics 4, 37–40 (2010). [CrossRef]
23. L. Razzari, D. Duchesne, M. Ferrera, R. Morandotti, S. Chu, B. E. Little, and D. J. Moss, “CMOS-compatible integrated optical hyper-parametric oscillator,” Nat. Photonics 4, 41–45 (2010). [CrossRef]
24. X. Liu, R. M. Osgood, Y. A. Vlasov, and W. M. J. Green, “Mid-infrared optical parametric amplifier using silicon nanophotonic waveguides,” Nat. Photonics 4, 557–560 (2010). [CrossRef]