Imaging using quantum noise properties of light

We show that it is possible to estimate the shape of an object by measuring only the fluctuations of a probing field, allowing us to expose the object to a minimal light intensity. This scheme, based on noise measurements through homodyne detection, is useful in the regime where the number of photons is low enough that direct detection with a photodiode is difficult but high enough such that photon counting is not an option. We generate a few-photon state of multi-spatial-mode vacuum-squeezed twin beams using four-wave mixing and direct one of these twin fields through a binary intensity mask whose shape is to be imaged. Exploiting either the classical fluctuations in a single beam or quantum correlations between the twin beams, we demonstrate that under some conditions quantum correlations can provide an enhancement in sensitivity when estimating the shape of the object.


Introduction
The use of quantum correlated resources to enhance sensing and measurement has been an active field of research in recent years [1]. Multi-spatial-mode quantum states of light have become available [2][3][4][5] and have allowed for the development of imaging techniques that can surpass classical limits on resolution [6][7][8] or enhance signal-to-noise [9][10][11]. Here we present a technique that exploits the noise properties of multi-spatial-mode vacuum-squeezed twin beams to detect the shape of an object being probed by one of the beams. While we are able to estimate the shape of the object by monitoring only the fluctuations of the beam that probes it, we obtain an enhancement in this estimation sensitivity under certain conditions by exploiting the quantum correlations between the twin beams. It is possible to estimate the shape of a binary (hard-apertured) intensity mask by probing it with a light field and performing a homodyne detection of the light that passes through it. The local oscillator (LO) spatially selects which portions of the probing field are detected, so the homodyne detection effectively acts as an imaging device [12,13]. By optimizing the shape of the LO to maximize the mode matching with the probing field after the mask (the LO-mask overlap), the object's shape can be inferred. If the mask is probed with a state of light consisting of a small number of photons, determining the mean value of the field becomes difficult due to unavoidable quantum noise. We can, however, use the same homodyne technique to measure the noise of the field rather than its mean value in order to extract information about the object's shape. Homodyne detection provides the additional advantage of allowing the discrimination of a signal against a background at other wavelengths.
A coherent state would not be a suitable choice to probe the mask since it is a displaced vacuum state. Its noise properties would thus remain unchanged regardless of how the probe field is altered or attenuated by the mask [14]. On the other hand, a state whose quadrature fluctuations are not at the shot noise limit (SNL) can provide information about the shape of the mask. For example, if a thermal state is used, the mask's shape can be estimated by manipulating the shape of the LO to optimize the geometric overlap between the field that passes through the mask and the LO to maximize the detected level of excess noise relative to the SNL. Since thermal states are classically obtainable, we will refer to such a single-beam method as a classical noise imaging technique.
In many situations a reference beam is used to eliminate correlated sources of noise and to improve the signal-to-noise ratio for a given measurement [15,16]. We extend this idea to noise imaging by exploiting the quantum correlations between vacuum-squeezed twin beams to improve the estimation of the shape of the mask over the classical noise imaging technique described above. We refer to this two-beam technique as a quantum noise imaging technique.
In this paper, we use twin beams generated by four-wave mixing (4WM) to directly compare these two techniques and find that the quantum technique provides an enhanced sensitivity to changes in the LO-mask overlap and thus a greater confidence in estimating the shape of the mask.

Four-wave mixing and squeezed light detection
To generate the light states needed for this experiment, we have used 4WM in a double-Λ configuration [17] in a hot 85 Rb vapor ( Fig. 1(a)-(b)). The 4WM process converts two photons from a single pump beam into a pair of photons emitted into twin fields referred to as the probe and conjugate. These twin beams exhibit strong amplitude correlations such that the amplitude difference noise is below the SNL [18]. Taken individually, these fields are thermal states exhibiting uniformly distributed excess noise. Quantum correlations between the probe and the conjugate can be measured by performing a balanced homodyne detection of each field and subtracting the homodyne signals for appropriately chosen phases of the LOs. The noise of a single mode's generalized quadraturê X θ = 1 √ θ p and θ c ,X J θ is squeezed under the condition that θ p + θ c = π. We therefore only need to control the phase of one of the LOs (using a PZT-mounted mirror) to detect the level of squeezing between the probe and conjugate.
To implement the 4WM process, a linearly polarized intense pump beam (300 mW) is focused down to a 1 e 2 diameter of 1.2 mm inside a 1.25 cm long 85 Rb cell heated to 110 • C. The probe and conjugate frequencies are unseeded, so the resulting probe and conjugate fields are generated from spontaneous emission. The frequency of the pump beam is detuned 800 MHz to the blue of the |5 2 S 1/2 , F = 2 → |5 2 P 1/2 , F = 3 transition at 795 nm. After separating the pump beam from the probe and conjugate beams with a polarizing beam splitter, each beam is sent to a balanced homodyne detector using a pair of matched photodiodes with quantum efficiencies of approximately 95%. The total optical path losses from the 4WM process to the homodyne detection are approximately 4%.

Experimental procedure
To characterize each technique's sensitivity we reduced the LO search to a one-dimensional problem by choosing our mask to be in the shape of a bow tie that can be oriented at an arbitrary angle around its center of symmetry ( Fig. 1(c)). The position of the center of symmetry was fixed, and the only free parameter to be estimated was the mask's angular orientation. We generated probe and conjugate LOs with congruent bow tie shapes but rotated at an arbitrary common angle relative to the mask. The estimated parameter is the overlap between the conjugate after the mask and its LO, which is a function of the angle. We then recorded the resultant quadrature noise powers with a spectrum analyzer as a function of the overlap for the classical and quantum noise imaging techniques. In order to reliably measure the amount of squeezing in the quantum case, the phases of the LOs were quantum noise locked [19] such that the twin beam quadrature noise difference was always minimized for each overlap value.
An important feature of the correlated fields generated by our 4WM set-up is that they are multi-spatial-mode [20], which can in general provide greater flexibility in obtaining an enhancement with these states [14]. Quantum correlations between the probe and conjugate are localized to pairwise correlated regions within the beams referred to as coherence areas [21,22] ( Fig. 1(c)). In order to measure and compare only correlated coherence areas we symmetrically generate LOs for the probe and conjugate by shaping the probe beam using a spatial light modulator (SLM) to shape the spatial profile of a coherent state at the probe frequency. We then use this beam to seed a second 4WM process (Fig. 2). This technique [23] ensures that only the spatial modes supported by the process will be present in the LOs and that any mode distortions due to Kerr lensing [20] will be accounted for as the beams propagate to the far field. Additionally, the use of a second 4WM process guarantees that the spatial profiles of the probe and conjugate LOs are manipulated synchronously so that only correlated regions of the probe and conjugate fields are detected as the shapes of the LOs are manipulated. Without any mask present, we align the LOs to the twin fields in order to maximize the squeezing. The mask is then inserted into the conjugate field, and we begin the search for the mask's shape by manipulating the transverse profile of the LOs.
We calibrated the LO-mask overlap by splitting the seed beam used to generate the LOs and sending it to seed the 4WM process formerly used to generate the vacuum-squeezed twin beams, so that bright twin beams are obtained in their place. These bright twin beams can then be interfered with their LOs in order to optimize mode matching and ensure proper initial alignment for the homodyne detection. As we rotated the angular orientation of the bow tie shapes on the SLM, we observed that the homodyne mode-matching efficiencies between the bright beams and their LOs stayed near 97% for both the probe and conjugate. We then inserted the mask into the bright seeded conjugate beam and directly detected the fraction of transmitted power as the angle of the bow tie was changed using the SLM. This procedure allowed us to factor any non-uniformity in the beam's intensity profile into the evaluation of the overlap at any LO orientation.
The noise generated by the homodyne difference signal was detected at 750 kHz using a spectrum analyzer in zero-span mode with a resolution bandwidth of 30 kHz, a video bandwidth of 100 Hz, and a sweep time of 1 s. The detected noise was digitized into 460 points, and the average noise power of these 460 points was taken to be a single quadrature noise measurement. To characterize each measurement's statistical uncertainty, we divided each trace into 46 segments consisting of 10 points each. We adopted this procedure after verifying that clustering the data into 10 point segments maintained statistical independence among tthese segments. The average value of each segment was tabulated, and the standard deviation of these 46 values was calculated to characterize the measurement's statistical uncertainty. Finally, 10 series of measurements were taken at 15 different local oscillator bow tie angles for both the quantum and classical noise imaging techniques. This data is presented for the classical and quantum techniques in Fig. 3(a). Since each series of measurements was separated by approximately 20 minutes, the data in Fig. 3(a) indicates an experimental stability over several hours.  We first consider the quantum noise imaging technique (blue triangles in Fig. 3). If the conjugate LO does not spatially overlap with the conjugate field that passes through the mask its measured fluctuations will be near the SNL. On the other hand, since the probe field passes unobstructed to its homodyne detector, its homodyne detection will produce excess noise. Therefore the homodyne difference signal will yield a net excess noise relative to the twin beam SNL. As the overlap between the conjugate and its LO increases, however, the quantum correlations between the twin beams are recovered. Assuming that the noise properties are identical for all the coherence areas, the search for the mask's shape is complete when the amount of squeezing is maximized for the largest LO area possible.
The classical noise imaging technique is implemented by simply blocking the probe beam (red squares in Fig. 3). The only difference is that the conjugate beam, taken alone, is a thermal state with excess noise. In this case, the search for the shape of the mask consists of maximizing the excess noise by manipulating the LO.

Results
We wish to quantify the relative sensitivity between the quantum and classical noise imaging techniques in detecting changes in the LO-mask overlap for different overlaps. To do this we introduce the relative uncertainty in the overlap estimation ∆O est : ∆N represents the measured standard deviation of a given noise power N and ∂ N ∂ O is the slope of the noise power as a function of the overlap. In other words, ∆N quantifies the "noise on the noise," which incorporates both sources of statistical uncertainty and technical noise. ∂ N ∂ O is set by both the spatial mode composition of the light illuminating the mask as well as the mask's shape [14]. To estimate ∂ N ∂ O , we computed a linear fit to the noise data at overlaps greater than 0.8 and extrapolated to an overlap of unity. We then included this extrapolated point in a 3rd order polynomial fit and plotted the result in Fig. 3(a). Fig. 3(a) confirms that for overlaps near unity the excess noise and squeezing are maximized for the classical and quantum techniques, respectively. In Fig. 3(b), we present curves comparing ∆O est for the classical and quantum schemes. As Fig. 3(b) suggests, the quantum noise imaging technique provides a higher sensitivity to changes in the overlap than the classical technique. Although the comparable magnitudes of Fig. 3(a) might suggest that the quantum and classical cases should exhibit comparable values for ∆O est , the difference between the curves in Fig. 3(b) can be explained by the variation of ∆N with overlap. Specifically, the value of ∆N is expected to scale with the noise power N. For a given overlap the magnitudes of | ∂ N ∂ O | are similar for each technique, but the smaller value of ∆N for the quantum technique will yield a smaller value of ∆O est .
Since the objective of these techniques is to estimate the shape of the mask, we wish to operate in the regime where the overlap is close to unity. For overlaps of 0.9 and above, the quantum case provides an enhancement over the classical case by a factor of 6.3 ± 0.4, corresponding to enhancement in the estimation of the angle by a factor of 3.8 for small angles. The enhancement factor of 6.3 was calculated by averaging ∆O est for overlaps above 0.9 in Fig. 3(b) for the quantum and classical techniques. The angular enhancement factor of 3.8 was obtained by separately calibrating the measured LO-mask overlap as a function of the angle between the LO bow tie and the mask. The indicated uncertainties here and in the figures represent one standard deviation, combined statistical and systematic uncertainties.
We have verified that the classical technique is limited by fluctuations of the conjugate field. On the other hand, we are able to observe excess technical noise introduced by the quantum noise lock used for the measurement of the squeezed noise power for the quantum technique. This leads to a degradation in sensitivity when estimating the angle of the mask. Since this technical noise is not present on the excess noise measurements of the classical technique it decreases the enhancement provided by quantum noise imaging technique over the classical technique. Accordingly, even larger enhancements should be possible. The advantage that can be obtained with the quantum technique over the single-beam classical technique is dependent on the amount of squeezing present [14]. Losses after the twin beam source will erode this advantage depending on how they are distributed between the probe and conjugate. In particular, unbalanced losses between the beams would alter the sensitivity and could eliminate any advantage.

"Alphabet gun" test
While we have shown that the quantum noise imaging technique can improve the sensitivity in determining the maximum LO-mask overlap, we also show that this technique can be used with a simple search algorithm to recognize a mask among a given set of choices. We place a mask cut into the shape of the letter Z into the path of the conjugate field and then shape the probe and conjugate LOs with the SLM to form the various letters of the Latin alphabet ( Fig. 4(a)). All letters are positioned such that their glyph widths exactly overlap (Fig. 5). Each letter essentially functions as a guess for the shape of the mask. With no mask inserted into the path of the conjugate beam, each choice of LO letter yields slightly different levels of detected single beam excess noise and twin beam squeezing due to the spatial mode structure of the twin beams. These noise levels serve as baselines for comparison when the mask is inserted. We then look at the deviation from these baselines upon insertion of the mask: D i denotes the deviation in the noise power for LO letter i, N i masked is the noise power recorded by the homodyne detection with the mask present, and N i baseline signifies the unmasked baseline noise power.
For the classical technique, the presence of the mask leads to varying degrees of reduction in the measured excess noise (Fig. 4(b)). For the quantum technique, the mask causes the measured quadrature-difference noise to increase and, for most choices of LO letter, for the squeezing to be lost. In fact, only the correct choice of the LO letter (Z) continued to yield any squeezing. It should be noted that the letter I in the chosen font did not produce a strong enough LO to elevate the single beam shot noise significantly above the electronic noise floor, so it is not included in Fig. 4(b)-(c).
Although both the classical and quantum estimation techniques correctly suggest that the letter Z is the best choice of LO, the quantum technique provides a higher degree of confidence. Uncertainties in the deviation incorporate contributions from both the baseline and masked noise measurements. These uncertainties are plotted along with their associated changes in squeezing and excess noise in Fig. 4. For the quantum case (Fig. 4(c)), the correct choice of LO is separated from the next closest choice (the letter T) by over 7 standard deviations. For the classical case, however, the quadrature noises for the letters T and Z lie within their respective uncertainties.

Conclusion
In conclusion, we have demonstrated the ability to estimate the shape of a binary intensity mask using homodyne detection of the noise of vacuum-squeezed twin beams used to illuminate the object. In the classical and quantum techniques described, the mask is exposed to exactly the same state of light, which contains too few photons for direct detection of the amplitude quadrature and too many photons for current photon counters. These techniques might therefore prove useful for imaging applications where a minimal exposure to light is critical (e.g. when dealing with low optical damage thresholds or when it is desirable that the detection go unnoticed). A promising application is related to imaging of biological samples under controlled conditions where squeezing could be maintained and the exposure of the sample minimized [24]. The classical noise imaging technique exploits the fact that, taken individually, each of the twin beams behaves like a thermal state of a light characterized by uniform excess quadrature noise. Although the excess noise properties of a single beam can be used, we have demonstrated that the quantum correlations of the twin beams offer an enhancement in sensitivity for the regime of high LO-mask overlap. As explained in [14], these correlations must be quantum for an enhancement to be possible. Additionally, we showed that a mask can be recognized among a given alphabet of choices with increased confidence using the quantum noise imaging technique. As one can notice from Fig. 5, the resolution with which we generate the LOs is limited. The constraint on the resolution is the limited pump power available for the 4WM process, which places a practical upper bound on the pump's available transverse size. Increasing the pump's transverse size at a given intensity would increase the resolution of the LOs [11] and should improve our ability to estimate the shape of the mask.
Finally, we note that one could in principle identify a phase mask separate from, or in addition to, an intensity mask with a similar technique. The search for a phase mask would require a scan over the LO phase in only the beam that passes through the mask. Such a search would not be possible, however, using the "classical technique" since for a thermal beam all phases have the same noise level.