# Multimode Quantum State Tomography Using Unbalanced Array Detection

###### Abstract

We measure the joint Q-function of a multi-spatial-mode field using a charge-coupled device array detector in an unbalanced heterodyne configuration. The intensity pattern formed by interference between a weak signal field and a strong local oscillator is resolved using Fourier analysis and used to reconstruct quadrature amplitude statistics for 22 spatial modes simultaneously. The local oscillator and signal propagate at an angle of thus shifting the classical noise to modes that do not overlap with the signal. In this configuration, balanced detection is not necessary.

###### pacs:

42.50.Ar, 03.65.Ta, 03.65.Wj]www.amcdawes.com

Multimode quantum optical systems are increasingly important in applications such as slow and stored light devices Grodecka-Grad et al. (2012). A complete understanding of these devices requires a fidelity measurement of the stored quantum state, but there has not yet been a simple way to measure the quantum state of a multi-spatial mode field. Balanced homodyne detection (BHD) is the standard measurement technique used to construct the complete quantum mechanical state of light Leonhardt (1997). This technique has one primary weakness: multimode fields or fields in an unknown spatial mode suffer significant losses due to mode mismatch between the local oscillator (LO) and the signal field Leonhardt (1997); Beck (2000).

Balanced array detection overcomes mode-matching losses and has been used to perform measurements of more than two modes at once Raymer et al. (1993). It is possible to obtain simultaneous (but not joint) measurements of the Wigner functions of many modes by using array detectors Dawes and Beck (2001); Dawes et al. (2003). It has also been shown theoretically that array detectors may be used to measure the joint Q-function of a multimode field Iaconis et al. (2000). One of the difficulties of these methods is the need to balance two array detectors and perform pixel-by-pixel subtraction to eliminate the classical intensity noise of the LO Yuen and Chan (1983). These alignment requirements become prohibitive when trying to measure two transverse directions. This difficulty was removed by Beck et al. who achieved the same effect with a single array, using unbalanced detection, by arranging the LO and signal fields to isolate the classical LO noise Beck et al. (2001).

In this Letter, we describe a multi-spatial-mode detection method that is analogous to the multi-temporal-mode detection method of Beck et al. We simultaneously measure many spatial modes without having to vary the propagation direction or the phase of the LO; this reduces the amount of data that needs to be acquired and enables simultaneous joint measurements of two modes. This measurement technique can measure joint statistics between two independent spatial modes. Knowledge of these inter-mode correlations could provide a new route to information storage in a wide variety of systems.

To develop the theory of unbalanced array detection, we consider the spatial intensity incident on the array detector due to the combination of an LO field and the signal field . For simplicity, we describe only one transverse coordinate (i.e. horizontal); generalization to both transverse coordinates is straightforward. If the LO is normally incident on the detector and the signal field propagates at a small angle to the LO, the spatial intensity at the detector is given by

(1) |

The Fourier transform of this measurement yields

(2) |

where and denotes the convolution. The first of the three terms peaks at and contains the second-order classical LO noise that would be removed in balanced detection. The second term is negligible if the signal field is weak. If is large enough, the function contains all information on the heterodyned signal. This function has a peak value at and is therefore separate from the classical noise at . Just as Beck et al. Beck et al. (2001) eliminate classical LO noise by temporally separating the signal and LO fields, we eliminate classical LO noise by separating the LO and signal fields by their propagation direction—they have different transverse components of their propagation wave vector Dawes et al. (2013).

To describe the detection process from the perspective of quantum mechanics, we consider the spatial pattern detected at an array, as in Fig. 1, where is the operator corresponding to the number of photons incident on pixel of the array. Each pixel measures a different spatial position and there are pixels labeled by . We can also express the measured field in terms of plane-wave modes , using the Fourier relation

(3) |

We assume the LO occupies the plane-wave modes near the center of the spectral window. The signal occupies the plane-wave modes with positive and the plane-wave modes with negative contain only the vacuum field. This is explicitly written as

(4) |

The operator we measure corresponds to the inverse Fourier transform of

(5) |

Combining Eqns. 3 and 4 yields an expression for the per-pixel number operator in terms of plane-wave mode operators . Terms that combine two weak fields and are discarded as second-order. From here, Eqn. 5 is used to relate in terms of the ’s. Finally, the large-amplitude of the LO field means the LO operators can be replaced with their coherent state amplitudes . For we find

(6) |

If the LO occupies only a single () plane-wave mode, then Eq. 6 simplifies to

(7) |

For an LO in a single plane-wave mode, a measurement of (the Fourier transform of the photocount data) returns a complex number. Using Eq. 7, this complex number can be interpreted as a measurement of , the signal mode annihilation operator, plus a vacuum contribution (the second term in Eq. 7). The annihilation operator is itself the sum of the two field quadrature amplitudes

(8) |

Therefore, the real and imaginary components of each correspond to simultaneous measurement of the quadrature amplitudes and . Of course these observables are noncommuting so the ability to measure them simultaneously comes at the price of additional vacuum noise Arthurs and Kelly Jr (1965); Leonhardt (1997). It is precisely this additional noise that prevents reconstruction of the Wigner function, instead allowing reconstruction of the Q-function.

It is important to note that each Fourier-transformed exposure of the array returns a set of complex numbers. Each of these, indexed by , has a real and imaginary part that corresponds to the signal mode field quadrature . If one value of is selected, and the corresponding pairs are histogrammed, the result tends toward the Q distribution for the field quadratures of mode . Because the data for all values of are collected simultaneously, any mode can be analyzed from a single set of data. Additionally, joint Q-distributions can be computed for any pair or combination of modes.

Our implementation of this system, shown in Fig. 2, begins with a
external-cavity diode laser. The laser output is sent
through an acousto-optic modulator (80 MHz). The first order diffracted beam is
then spatially filtered using an optical fiber (780-HP). After the fiber output,
a telescope, starting with a 50-mm focal length lens, expands the beam diameter
to . A 50/50 beam splitter separates the LO from the signal
and a 250-mm focal length lens collimates the LO beam after expansion. The
signal beam remains slightly diverging after passing through a 300-mm focal
length lens. The mode mismatch between the signal and the LO provides a simple
multimode signal for our demonstration of this method. The signal is attenuated
by a factor of using neutral density filters. The signal and LO interfere
on a charge coupled device (CCD) camera ^{1}^{1}1Princeton Instruments, PyLoN
400BR eXcelon.. The CCD is a 1340400 array of square pixels with quantum efficiency of 98% at
, cooled to to achieve a dark current
rate of 2-3 electrons per pixel per hour.

The LO power incident on the CCD is and the LO and signal beams interfere on the CCD at an angle of . The AOM signal is modulated by a function generator (Rigol DG4062) to create a rectangular pulse. The pulse is triggered after the start of the exposure to keep the CCD dark during readout.

The largest noise source in our measurements is the electronic readout noise. We characterize this noise by measuring the variance in photocounts with the LO illuminating the CCD and without illumination. The variance with illumination is an average of above the variance without illumination so our signal to noise ratio is more than adequate.

If the signal is blocked, the signal mode entering the detector is then the vacuum mode, and the quadrature amplitudes measured should be zero. Despite our spatial filtering, the residual spatial components in the LO beam leave a background signal in the FFT output that should be subtracted before computing the Q-function. To carry out this correction, we collect 500 exposures with the signal blocked (i.e. the signal in the vacuum state) and average the FFT output across these 500 exposures. This vacuum average is then subtracted from subsequent data prior to further calculations.

We collect 1000-8000 shots of CCD data from an ROI at the center of the CCD (600
pixels wide by 10 pixels tall), sum vertically the ten rows to obtain data points and compute the FFT of the resulting 600-element array. The
output of this calculation is an array of 600 complex values of where the
first 300 values correspond to unique modes indexed by ^{2}^{2}2Given that
the raw CCD data is real, the FFT output is symmetric about the midpoint..

We compute the quadratures by scaling the data by the total number of photons detected per shot in the CCD region of interest. For weak signal fields, is essentially equal to the average number of photons in the LO and therefore related to the classical LO amplitude as

Therefore, the quadratures in Eq. 8 are calculated from experimental values of as

(9) | |||

(10) |

From the quadratures, we compute the Q-function for a specific mode . This calculation is done in two ways. First, we simply histogram the and pairs using 30 bins in each dimension. For large data sets, this histogram will approximate . Additionally, we calculate by applying a kernel density estimate (KDE) to the and pairs. This method essentially places a narrow Gaussian function at each point and computes the sum. The width of the Gaussian kernel is determined using Scott’s rule and results in some smoothing relative to our raw histogram Scott (2009).

In the following figures, both the raw-data histogram and the KDE versions of the Q-function are presented for comparison. Shown in Fig. 3 is for a mode with photons and . The Q-function can be computed for any mode using the same set of data. Because data is collected for all modes simultaneously, we can also compute a quantity of interest (such as average photon number and measure that quantity for all modes simultaneously. In Fig. 4 we plot average photon number vs. plane wave mode angle . This illustrates the mode spectrum of the signal field. In particular, with a multimode field there are several modes with comparable amplitudes. While Fig. 4 shows the 22 modes near the signal, we measure 600 modes simultaneously; highly multimode signals can be measured using this technique.

To demonstrate the detection of a signal with variable phase, we dither the position of a piezo-mounted mirror in the signal beam path and observe modulation through of phase. The Q-function calculated from such data is shown in Fig. 5. The characteristic donut shape corresponds to phase noise although the phase modulation is not completely random as evidenced by the peaks present in the Q-function.

It is important to note that we can also compute any quantity that is a function of the quadratures and . One such quantity is the joint Q-function, . The joint Q-function shown in Fig. 6 is for two nearby modes and . As expected, these modes exhibit strong positive correlation between and .

In conclusion, we point out that this method could be used to measure the quantum state of light stored in a multimode quantum memory such as that described in Grodecka-Grad et al. (2012). The peak sensitivity of our CCD system is near the and resonances of Rubidium. Very recently, quantum state tomography measurements have been performed on single-photon states retrieved from a stored-light system Bimbard et al. (2014). Multimode measurements of such systems could be made using methods presented in this Letter and would yield new information about correlations between modes retrieved from stored light systems.

###### Acknowledgements.

We thank M. Beck for helpful discussions. This material is based upon work supported by the National Science Foundation under Grant No. 1205828. Additional financial support was provided by the Research Corporation for Science Advancement and the Pacific Research Institute for Science and Mathematics.## References

- Grodecka-Grad et al. (2012) A. Grodecka-Grad, E. Zeuthen, and A. S. Sørensen, Phys. Rev. Lett. 109, 133601 (2012).
- Leonhardt (1997) U. Leonhardt, Measuring the quantum state of light, Vol. 22 (Cambridge Univ Pr, 1997).
- Beck (2000) M. Beck, Physical Review Letters 84, 5748 (2000).
- Raymer et al. (1993) M. G. Raymer, J. Cooper, and M. Beck, Physical Review A 48, 4617 (1993).
- Dawes and Beck (2001) A. M. Dawes and M. Beck, Physical Review A 63, 040101 (2001).
- Dawes et al. (2003) A. M. Dawes, M. Beck, and K. Banaszek, Physical Review A 67, 032102 (2003).
- Iaconis et al. (2000) C. Iaconis, E. Mukamel, and I. A. Walmsley, Journal of Optics B: Quantum and Semiclassical Optics 2, 510 (2000).
- Yuen and Chan (1983) H. P. Yuen and V. W. S. Chan, Opt. Lett. 8, 177 (1983).
- Beck et al. (2001) M. Beck, C. Dorrer, and I. A. Walmsley, Phys. Rev. Lett. 87, 253601 (2001).
- Dawes et al. (2013) A. M. C. Dawes, N. T. Holte, and H. A. Dassonville, Proc. SPIE 8636, 863610 (2013).
- Arthurs and Kelly Jr (1965) E. Arthurs and J. Kelly Jr, Bell Syst. Tech. J 44, 725 (1965).
- (12) Princeton Instruments, PyLoN 400BR eXcelon.
- (13) Given that the raw CCD data is real, the FFT output is symmetric about the midpoint.
- Scott (2009) D. W. Scott, Multivariate density estimation: theory, practice, and visualization, Vol. 383 (John Wiley & Sons, 2009).
- Bimbard et al. (2014) E. Bimbard, R. Boddeda, N. Vitrant, A. Grankin, V. Parigi, J. Stanojevic, A. Ourjoumtsev, and P. Grangier, Phys. Rev. Lett. 112, 033601 (2014).