Fig. 1. (a) Phase sensitive gain in a degenerate OPO – the amplitude incremental gain of the signal (g_|s|) versus the relative phase between the pump and the signal – (θ_s = ϕ_s − ϕ_p/2 +π/4, where ϕ_s and ϕ_p are the signal and pump phases, respectively), (b) Schematic of the all-optical quantum RNG. PD: Photodetector, BS: non-polarizing beam splitter, (c) Electric fields of the signal pulses compared to the pump in the two possible phase states of a short-pulse degenerate OPO.

Fig. 2. (a) Schematic of the experimental setup consisting of the twin degenerate OPO followed by an un-equal arm Michelson interferometer. The round-trip length difference of the arms in the interferometer is equal to the separation length of pump pulses. The output of the interferometer is filtered to eliminate the pump. An AOM is used to restart the twin OPO at a clock rate of 10 kbps, while the zeroth order output (non-diffracted beam) pumps the OPO. The detector PD1 and a piezoelectric transducer (PZT) on M₄ are parts of the “dither-and-lock” servo loop for stabilizing the cavity length to a resonance [16]. The interferometer arm length difference is stabilized using the interference of the undepleted pump and a PZT on M₁₁ (not shown in this schematic) resulting in a fringe visibility of greater than 60% throughout the experiment. OC: Output Coupler, (b) Pump and signal pulse trains for a twin OPO in which the cavity roundtrip time (T_cavity) is twice the repetition period of the pump (T_rep); this results in two independent OPOs with signal pulse trains at half the repetition frequency of the pump, marked by OPO1 and OPO2 on the pulse trains.

Fig. 3. (a) Optical spectrum of the pump and signal, (b) complementary fringe patterns at the output of the Michelson interferometer, i.e. maximum on one pattern is minimum on the other and vice versa, (c) (Media 1) sample time domain signal of the interferometer output, i.e. the random sequence, along with the clock signal at 10 kbps.

Fig. 4. Summary of the results of the NIST statistical tests [19] on a sequence of 1 billion bits. The sequence is chopped into 500 shorter sequences and each test is performed on all sequences. Horizontal axes shows the name of the tests, P – Value is the result of χ² distribution evaluation for all short sequences, and Proportion is the proportion of the short sequences passed each test with a significance level of 0.01. The final decision on if the whole sequence has passed a test is made based on these criteria: P – Value > 0.0001 and Proportion > 0.976 [19], which are satisfied for all the tests.

Fig. 5. (a) Dependence of the randomness on the clock speed; for each clock speed a sequence of 100 kb is recorded, and the rate of bit flip is calculated by counting the number of changes in the bit value and dividing it to the sequence length, for a random sequence the rate of bit flip is close to 0.5 and for a single-value sequence it is zero; The transition from random to non-random output occurred around 25 kbps which is close to the estimated value of 30 kbps taking into account the experiment parameters. This value can be increased by decreasing the overall pump power and approaching the threshold. To record the 1 Gb of data, the pump power was reduced so that the randomness would not break until 200 kbps (the measurement limit). (b, c, d) a sample of the output sequence when clock speed is (b) 20 kbps, (c) 30 kbps, and (d) 40 kbps.