Interferometric control of the photon-number distribution

We demonstrate deterministic control over the photon-number distribution by interfering two coherent beams within a disordered photonic lattice. By sweeping a relative phase between two equal-amplitude coherent fields with Poissonian statistics that excite adjacent sites in a lattice endowed with disorder-immune chiral symmetry, we measure an output photon-number distribution that changes periodically between super-thermal and sub-thermal photon statistics upon ensemble averaging. Thus, the photon-bunching level is controlled interferometrically at a fixed mean photon-number by gradually activating the excitation symmetry of the chiral-mode pairs with structured coherent illumination and without modifying the disorder level of the random system itself.

Optical interferometers implement deterministic field transformations that trace an interferogram by sweeping a phasewhich enables applications across all of optics and photonics 1 . The detected intensity is modulated, but the underlying photon number distribution P n does not change: the Poisson statistics associated with coherent light remain Poissonian, and the Bose-Einstein statistics that are the hallmark of thermal light remain Bose-Einstein. Introducing dynamical randomness in the interferometer can change the photon statistics. An early example in neutron interferometry implemented a chopper in one arm of a two-path interferometer 2,3 , while an optical 'stochastic interferometer' devised by De Martini et al. changes the photon statistics by inserting a randomly varying optical element [4][5][6][7] . Along a different vein, coherent light traversing a random medium may gradually acquire chaotic behavior and under certain conditions may even attain Bose-Einstein statistics. This phenomenon has been observed in several systems, including weakly transmitting barriers embedded in a waveguide 8 , quasi one-dimensional disordered samples 9 , and optically dense slabs containing multiple scatterers [10][11][12][13][14][15] . In these approaches, changing the photon statistics is predicated on varying the system disorder level -with the amount of fluctuations (photon bunching) typically proportional to the disorder level.
We present here a different strategy for tuning photon statistics in a random structure that does not require modifying its disorder level. Instead, the ensemble-averaged P n is varied determinstically by an external element in the form of a relative phase between two input channels of a random network -analogously to a traditional two-path interferometer. Reconstructing P n reveals deterministic interferometric tuning of photon-bunching across the super-thermal and sub-thermal regimes with fixed mean photon-numbern. This effect has its origin in the smooth and deterministic breaking of the excitation symmetry in certain random lattices, which can be achieved by sculpting the lattice excitation. Our experiment makes use of a one-dimensional (1D) disordered photonic lattice of evanescently coupled waveguides with offdiagonal disorder 16 . Additionally, we measure two statistical quantities: the normalized second-order photon correlation g (2) , which does not depend onn, and Mandel's Q-parameter, which does 17 . In previous work 18 , we measured the statistical parameter g (2) . Although g (2) can be computed from the photon-number distribution, the converse is not true. The counting distribution has features that are not captured by a single parameter, its normalized second-order moment g (2) . The novelty in our work lies in revealing the tunability of the shape of the distribution itself, by demonstrating the transition between a Poisson distribution, to a Bose-Einstein-like distribution, and the emergence of other intermediate photon counting distributions.
The concept of interferometric control over P n is depicted in Fig. 1. We start by considering coherent light with a Poissonian distribution P n =n n e −n /n!, n=0,1,2,···, that is fed into a single channel of a disordered lattice consisting of an array of coupled optical elements (waveguides [19][20][21] , resonators [22][23][24] , or fiber loops 25,26 ). The output P n can be changed by varying the lattice disorder level ∆ (to be defined below). Surprisingly, high-order coherences do not decline while increasing ∆ [Ref. 27 ]. If the lattice is endowed with disorder-immune 'chiral symmetry' [28][29][30] , a photonic thermalization gap 31 emerges upon ensemble averaging: the regime of sub-thermal photon statistics is forbidden at any disorder level, while super-thermal statistics are inaccessible to lattices lacking chiral symmetry 31 . We make use of the normalized second-order photon correlation g (2) = n(n−1) / n 2 ( . denotes ensemble averaging over both the quantum state and the lattice disorder) as a scalar measure of randomness to delineate the sub-thermal 1<g (2) <2 and super-thermal g (2) >2 regimes 32 . Instead of a monotonic trend towards 'thermalization' with increasing ∆ in lattices characterized by chiral symmetry (such as those with off-diagonal disorder 16,33 ), P n exhibits super-thermal statistics with a gradual decline towards thermal statistics upon increasing ∆ [ Fig. 1(a)]. Such a lattice can thus tune the photon statistics in the superthermal regime. Alternatively, in lattices lacking chiral symmetry (such as those with diagonal disorder 34,35 ), P n exhibits sub-thermal statistics with a gradual decline towards Poisson statistics upon increasing ∆ [ Fig. 1(b)]. This lattice can thus tune the photon statistics in the sub-thermal regime.
To tune the photon statistics without altering ∆, we sculpt the excitation over multiple lattice sites; e.g., by varying the relative phase between two sites. In traditional interferometry, a relative phase θ is introduced between two fields before they are superposed to trace an intensity interferogram 1 . In the interferometric scheme introduced here, a relative phase θ is introduced between adjacent channels of a random network with chiral symmetry fed with coherent light. The two fields superpose within the network and P n (θ) measured in a single output channel reveals a deterministic tuning of P n (θ) between super-thermal and sub-thermal regimes while sweeping θ [ Fig. 1(c)].
The above-described effect stems from the properties of the eigenmodes and eigenvalues of a lattice endowed with chiral symmetry. To appreciate the underlying physics, consider a generic tight-binding lattice model with nearest-neighbor-only coupling. The complex field amplitude E x (z) at site x after traveling along z is described by a set of coupled differential equations, where β x is the wave number for site x and C x,x±1 the coupling coefficient between site x and site x±1 -all of which may be arXiv:1706.08243v1 [physics.optics] 26 Jun 2017 x , whereβ is the average wave number. Since H is real and symmetric, {ψ (j) x } are all real. If H can be recast into a block-diagonal form after settinḡ β=0 in the interaction picture, this indicates that the lattice is endowed with chiral symmetry. Lattices with off-diagonal disorder (β x =β and C x,x±1 are random) satisfy chiral symmetry, whereas lattices featuring diagonal disorder (C x,x±1 =C and β x are random) do not. Henceforth, we focus our attention on lattices with offdiagonal disorder. A consequence of chiral symmetry is that x . However, for the impact of this skewsymmetry to be manifested, the members of the skew-symmetric paired modes with indices ±j must be excited with equal weights, which we refer to as 'activating' the chiral symmetry 18,31 . To visualize the activation of chiral symmetry, we associate a phasor with each mode at each lattice site. Then, in the rotating frame that is common to all modes, the eigenmodes appear in pairs whose members are equal magnitude counter-rotating phasors. This does not occur when the chiral symmetry is not activated or absent [ Fig. 1(d)].
This can be understood as follows. High-order coherences depend on the expected values of cross-correlations of the excited modal coefficients. In a disordered lattice, these typically all vanish except for the auto-correlation terms. However, in the presence of skew-symmetric modes in each realization of the random ensemble, the cross-correlation of their excitation coefficients survive statistical averaging, thereby adding a contribution to photon bunching and resulting in super-thermal statistics. This can be seen at both low ∆ via a perturbation analysis around the periodic lattice solutions, and high ∆ where localization reduces the number of lattice modes coupled to the excitation and thus also reduces the photon bunching 31,36 . One can also understand this behavior by examining the statistics of the field quadratures 37 . Under symmetric excitation configurations, the average field is constrained to a single quadrature in the complex plane. As a result, the intensity distribution approaches to chi-squared distri-bution with one degree of freedom 38 corresponding to modified Bose-Einstein photon number distribution 39 .
An input optical field E x (z=0) can be analyzed in a basis of lattice eigenmodes, E x . The field subsequently evolves after a distance z into E x (z)=∑ j c j ψ (j) x e ib j z . In the special case of a single-site excitation at the input E x (0)=δ x,0 , then |c j |=|c −j |=|ψ (±j) 0 | for all j, so that chiral symmetry is activated. This is not necessarily the case for more general field excitations. For example, when two adjacent sites are excited equally E x (0)=δ x,0 +δ x,1 , the modal coefficients are c ±j =ψ (j) 0 ±ψ (j) 1 , and chiral symmetry is activated |c j |=|c −j | only when the relative phase is θ=±π/2. Gradually varying the phase θ for fixed relative amplitudes (|A 0 |=|A 1 | in our experiment) tunes the chiralsymmetry breaking: maximal symmetry breaking at θ=0 or π, and symmetry activation at θ=±π/2.
The photonic lattice we utilized consists of an array of 101 identical 35-mm-long waveguides with nearest-neighbor-only evanescent coupling 40 . The average separation between the waveguides is 17¯m, resulting in an average coupling coefficientC≈1.71 cm −1 at a wavelength of 633 nm. The coupling coefficients are selected independently from a uniform probability distribution function with meanC and half-width ∆ in units ofC; our sample has ∆≈0.6.
The optical arrangement used in demonstrating deterministic interferometric control over P n is illustrated in Fig. 2. A singlemode coherent beam from a He-Ne laser is attenuated by a neutral density filter and split into two paths via a beam splitter. A relative phase shift θ is introduced via a delay in one path varied in 20-nm steps. The two beams are then brought together by a second beam splitter in parallel but closely spaced paths, which are imaged to two neighboring waveguides. The output facet of the array is imaged with a magnification of 8× via a lens (focal length f =3.5 cm) to a plane in which we scan a multimode fiber whose core optimally couples light from an individual waveguide. The multimode fiber delivers the collected light to a single-photoncounting module, and the output photon-number distribution P n is reconstructed using three photon-detection time windows: 20, 40, and 60 µs. The input intensity level is reduced to low levels The input laser is attenuated, split into two paths (beam splitter 1), a relative phase θ is introduced, the two beams are brought together in parallel paths (beam splitter 2), and then imaged into two adjacent waveguides in the array sample. The intensity distribution along the photonic lattice is plotted (in logarithmic scale) for a single realization of lattice disorder. The inset shows the coupling scheme, where E 0 and E 1 =E 0 e iθ are the coherent field amplitudes coupled to waveguides 0 and 1, respectively. Furthermore, a beam splitter at the output directs light that is imaged to a CCD camera, which helps ensure coupling into two neighboring waveguides for each disorder realizations (the camera and beam splitter are omitted for clarity).
so that only a few photons are detected within these windows while minimizing the accidental arrivals of multiple photons. We generate different disorder realizations by moving the array along the x-direction and coupling into a new pair of waveguides 16 . We measure single realizations of P n (θ) (averaged over 10 4 shots of the detection window) at the central lattice site x=0, and then average P n (θ) over an ensemble of 15 disorder realization for each value of θ by shifting the input excitation site.
We present in Fig. 3 our measurements confirming the deterministic interferometric tuning of P n (θ) in the excitation waveguide (x=0) while varying θ. As θ is swept, P n (θ) varies periodically (period π) between sub-thermal to super-thermal statistics [ Fig. 3(a)]. The measured mean photon-numbern(θ) does not vary with θ (n≈5.5 with a standard deviation ≈0.5). Althoughn remains constant, the photon-number distribution itself varies with θ, achieving maximal bunching when the chiral symmetry is fully activated θ=π/2,3π/2, and minimal bunching at θ=0,π when chiral symmetry is dormant. To better examine the salient changes in the photon statistics with θ, we plot in Fig. 3(b,c) P n (θ) for the extrema at θ=0 and π/2. The increase in the probability of higher photon numbers at the distribution tail when θ=π/2 as compared to that when θ=0 is a clear signature of photon bunching. In comparison, we also plot the Bose-Einstein distribution for thermal light with the experimentally determinedn≈5.5. At high n, the symmetric-excitation statistics exhibit higher probabilities than the Bose-Einstein distribution, whereas the broken-symmetric excitation has lower. Our experimental data [ Fig. 3(c)] is in excellent agreement with the simulations [ Fig. 3(b)].
Further analysis of the measured distributions P n (θ) helps bring about the changes that take place in the photon statistics. First, we examine a quantity extracted from P n (θ) that does not depend onn for the field considered here: the normalized second-order photon correlation function g (2) (θ) [Ref. 17 ]. We plot in Fig. 4(a) g (2) (θ) and note clearly that it varies sinusoidally with θ, between the sub-thermal and super-thermal regimes. The statistics are tuned from super-thermal (g (2) >2) to sub-thermal (g (2) <2), all while maintainingn fixed with θ [ Fig. 4(a), inset]. Althoughn changes according to the photon-counting window (n≈5. 5,11,16.5 for counting windows of 20,40,60 µs, respectively), the three interferometric traces of g (2) (θ) are indistinguishable.
Next we consider a quantity that does indeed depend onn: Mandel's Q-parameter, Q=Var(n)/ n −1=n(g (2) −1), which is thus linear inn for the field considered here 17 . Varying θ at a fixed n modulates Q(n;θ) between two limits identified by the dashed inclined lines in Fig. 4(b); we identify in Fig. 4(b) only the values (a) The measured P n (θ) while varying the phase θ between two coherent beams fed into adjacent sites of a disordered photonic lattice. The mean photon-numbern(θ)≈5.5 is independent of θ (with a standard deviation ≈0.5). (b) Simulated P n corresponding to θ=0 and π/2 obtained for n=5.5 and utilizing the physical parameters of the lattice. (c) Measured P n , extracted from (a) and corresponding to the simulations in (b). (a,b) Black-dotted line is Bose-Einstein distribution forn=5.5 corresponding to θ=0 and π/2, which are the extrema of this oscillation between the super-thermal and sub-thermal regimes. Increasingn (longer counting windows) leads to a linear growth in the two limits of Q modulation.
In conclusion, we have demonstrated that the photon-number distribution P n -and hence any photon statistic such as g (2) or Mandel's Q-parameter -can be tuned deterministically by varying a relative phase between two equal-amplitude beams launched into adjacent sites of a disordered photonic lattice. Such interferometric control over P n is possible by judicious excitation-  (2) (θ) and Mandel's Q-parameter. (a) g (2) as a function of θ obtained from the measured P n (θ) for photon-counting windows of 20, 40, and 60 µs. The g (2) -interferograms show no dependence onn. The inset showsn(θ) for the different photon-counting windows:n increases with the detection window but is independent of θ. (b) Mandel's Q-parameter as a function ofn at θ=0 (minimally activated chiral mode pairs) and θ=π/2 (maximally activated chiral mode pairs) obtained from the measured P n . The data points correspond to photon-counting windows in (a). The blue dashed line and the red dot-dashed lines are linear fits that fall in the super-thermal (blue-black color scheme) and sub-thermal (white-red color scheme) regimes, respectively. The border between the two color schemes demarcated by a solid gray line corresponds to Bose-Einstein statistics (thermal light). The dashed inclined lines correspond to horizontal dashed lines in (a). symmetry breaking of the skew-symmetric chiral mode pairs in disordered lattices that exhibit such symmetries.
Coupling light into a single site is a highly symmetric configuration. By sculpting the complex spatial distribution of the excitation over multiple sites, one may gradually break the excitation symmetry, which facilitates a smooth transition across the edge of the thermalization gap. This interferometric control strategy can be generalized in multiple ways. One may vary the relative amplitude of the fields exciting the two lattice sites while maintaining a fixed relative phase, or one may excite an extended section of the lattice while alternating the relative phases between adjacent sites. The latter approach has the advantage of producing the same phase-tunable P n with fixedn across the lattice.
Finally, we note that multiple experiments have been reported in which non-classical light is coupled into ordered or disordered photonic lattices [41][42][43][44] . In these reports, only small photon numbers are involved (usually two) in either Fock or entangled photon states, and while the output states different from those at the input, there have been no studies in the nature of the underlying photon statistics. Our work extends these studies into photon distributions with higher mean photon-number and demonstrates unambiguously interferometric control over P n itself.
A.F.A. was supported in part by the U.S. Office of Naval Research (ONR) under contract N00014-14-1-0260.