Twin-beam sub-shot-noise raster-scanning microscope with a hybrid detection scheme

By exploiting the quantised nature of light, we demonstrate a sub-shot-noise scanning optical transmittance microscope. Our microscope demonstrates, with micron scale resolution, a factor of improvement in precision of 1.76(9) in transmittance estimation gained per probe photon relative to an optimal classical version at the same optical power. This would allow us to observe photosensitive samples with nearly twice the precision,without sacrificing image resolution or increasing optical power to improve signal-to-noise ratio. Our setup uses correlated twin-beams produced by parametric down-conversion, and a hybrid detection scheme comprising photon-counting-based feed-forward and a highly efficient CCD camera.

Imaging that harnesses the quantum properties of light promises transformative new capabilities, including imaging an object's interaction with one wavelength of light but by detecting at another distinct wavelength [1], imaging below the diffraction limit [2], and reducing the uncertainty in the pointing of a laser beam [3].Quantum states of light are also known to enable higher precision absorption imaging [4][5][6] and phase imaging [7,8] relative to ideal classical light sources, which are ultimately limited by optical shot-noise.Sub-shot-noise measurements are particularly important for enhanced imaging and measurement of light sensitive photo-reactive biological samples, as higher precision measurements can provide better quality images for a given limited level of probe beam intensity [9].
The shot-noise of an ideal laser limits the minimum amplitude noise of an ideal classical light source, which is described by a Poisson distribution [10].Several reported approaches are able to surpass this limit, including the use of cooled and regularly pumped semiconductor lasers [11,12], generation of quadrature squeezed light [13], and using entangled Fock states [14].For the measurement of transmittance, the most precise sub-shot-noise approach per photon flux is the use of Fock states as a probe [15].These states can be readily generated by exploiting the photon-number correlations from twinbeam parametric processes such as spontaneous parametric down conversion (SPDC) [16] or spontaneous four-wave mixing (SFWM) [17].
For photon pair generation processes, the strong photonnumber correlations between each of the twin-beams enable intensity fluctuations of the probe light to be characterized or suppressed [18], which in turn enables a higher signal-tonoise ratio of the mean photon number of the beam passing through the sample for parameter estimation.However, optical loss degrades the intensity correlation of the generated light, such that this technique can only provide an advantage beyond classical techniques when the total loss is reduced be-low specific thresholds [4,19].Previous demonstrations of quantum sensing and imaging that did not satisfy these thresholds required postselection of successful detection events to observe the physics of sub-shot-noise performance [6,20], which underestimates the light exposure of the imaged sample.
Demonstrations of twin-beam sub-shot-noise single channel absorption estimation [21] showed how commercially available high quantum efficiency cooled CCD cameras enable precision beyond the ideal classical scenario where Poisson-distributed light is detected with a perfect 100% efficient detector.It has also been shown that fast optical gating, conditioned on single-photon detection events, can allow heralded photon counting twin-beam experiments to achieve subshot-noise absorption estimation without postselection [22], even when the components of the system have overall efficiencies below the required thresholds that would otherwise require post-selection.Here we combine optical gating of single-photons and detection using a low-noise, highefficiency CCD camera.
Combining these schemes increases the available quantum advantage offered by the individual constituent detection schemes and enables the use of photon-coincidence-detectionbased estimators with CCD camera intensity detection.The single-photon counting element allows for the removal of virtually all classically contributed fluctuations in the light illuminating the sample [22].This approach contrasts with methods based purely on intensity correlations where in practice the sample is illuminated with light comprising fluctuations that are partially of classical origin, due, for example, to losses in the reference photon channel.In the presence of such classical fluctuations, the best known intensity correlation strategy to extract maximum information about the sample is to use estimators that exploit both the quantum and the classical contributions to the light intensity fluctuations [21].However, compared to photon counting, this strategy is sub-arXiv:1906.05331v1[quant-ph] 12 Jun 2019 optimal because firstly, the classical contribution will be fundamentally limited by the shot-noise, and secondly it causes the estimator to be sensitive to technical fluctuations (that can be super-Poissonian) in the source intensity.In contrast, the new method reported here harnesses both the high efficiency of linear (non Geiger-mode) detectors and the classical noise suppression capability of optical gating with fast switching, conditional on single-photon detection events.
Many metrology applications seek to maximise measurement precision, which is statistically defined as the inverse of the variance of an experimentally estimated parameter.The variance is proportional to the uncertainty of experimental measurements used in the estimator.Here we focus on minimising the noise introduced into transmittance estimation caused by intensity fluctuations of an optical probe.
Fock states can be prepared from correlated twin-beams of either SPDC or SFWM, by separating signal and idler light into two modes deterministically, most commonly via wavelength or polarization selection.By detecting a photonnumber state (such as a single-photon) in one of the beams, the presence of its companion is heralded in the other beam and thus then be used as an optical probe.Any loss in the output Fock state's path degrades the heralded state into a mixture.However, this mixture still has a smaller photon-number variance than that of Poisson-distributed light with the same average intensity-it therefore can still be used to outperform classical light in transmittance estimation.
If the detector can resolve the time of arrival of singlephotons, then single-photon coincidence rates can be used to estimate transmittance of a sample using the ratio between the coincidence counts and the single counts of the reference photons ( N C /N R ).This parameter is also know as the heralding or Klyshko efficiency [23].To estimate the transmittance of a sample, η S , it would only be necessary to obtain the ratio of the of Klyshko efficiency of the probe path with (η P' ) and without (η P ) the sample in place [19]: Here N C is the number of coincidence counts between the reference and probe beam and N R is the number of single counts of the reference beam.Performance is then compared to that of a direct transmittance measurement with shot-noise limited light by estimating the ratio, Γ, of the precision obtained with each strategy.Sub-shot-noise (SSN) performance is indicated by Γ > 1.When the transmittance of the sample and the Klyshko efficiency of both the probe and reference channel are known, Γ can be estimated using From Eq. ( 2), we can see that SSN performance can only be achieved when η P + η R > 1, where η R is the Klyshko efficiency of the reference arm [19].
In order to use the estimator in Eq. ( 1), we use a singlephoton avalanche detector (SPAD) to implement feed-forward FIG.1: Simulation of precision ratio using twin-beams and feedforward.The plot compares Poisson-limited light (horizontal blue line at ratio of 1) to transmittance estimation with twin-beam strategies.The key refers to the different values of overall source Klyshko efficiency, which are symmetric for direct twin-beam exposure (dotted lines) and asymmetric for feed-forward (solid lines) due to the 15% loss in the switch and 10% leakage of unheralded photons.
with an optical delay and a fast optical switch so that probe photons are optically gated.This exposes the sample only when a herald photon was successfully detected, effectively reducing the losses of the reference channel and increasing the value of its Klyshko efficiency, η R , closer to 1, making it easier to satisfy the aforementioned condition for obtaining an advantage, as demonstrated in [22,24].We then use a CCD camera to detect the heralded probe photons.Since probe photons only reach the sample after being heralded, the number of photons that reach the CCD in the ideal case equals the number of photon coincidence events N C .Furthermore, using a modern CCD camera as a detector allows us to reach detection efficiencies that rival (albeit having more noise) cryogenically-cooled superconducting nanowire single-photon detectors [25], which have recently been used to achieve post-selection-free sub-shot-noise phase estimation with NOON states [14].
When considering realistic, imperfect components, any leakage of unheralded photons through the switch due to a poor extinction ratio will reduce the Klyshko efficiency of the reference channel η R < 1, while extra loss introduced in optical delays and optical switches will reduce η P .Therefore, when introducing an optical switch, there is a trade-off between increasing the Klyshko efficiency of the reference channel and introducing loss on the probe channel.Fig. 1 shows the precision ratio that can be obtained by implementing feed-forward on twin-beam sources with different levels of Klyshko efficiencies.For illustration, we have assumed that the loss on the probe beam, due to the feed-forward optics, is the same as that measured in our setup (15%) and has the same level of leakage of unheralded photons (10%).In Fig. 1 we observe, that feed-forward allows SSN performance to be achieved at lower transmittance than without feed-forward.It is also apparent that for Klyshko efficiencies below 70%, feed-FIG.2: Experimental setup and its characterisation.a) Photon-pairs are generated by SPDC in a PPKTP crystal, feed-forward of the idler photons is implemented using a polarisation independent switch, as described in Ref. [22].Heralded photons are focused first through the sample and subsequently refocused to a point on a CCD camera for detection.The sample is raster scanned in two dimensions, allowing an image to be reconstructed.A mechanical shutter before the crystal is triggered by activation of the camera sensor, initializing the experiment and allowing the pump beam into the crystal.The shutter prevents undesired exposure of the sample during the camera's dead time.The acquisition time of the SPAD is synchronized to the CCD using a software interface that communicates with a time tagging electronics system which only records counts during the acquisition time of the camera, minimising uncorrelated reference channel noise.b) Precision ratio of transmittance between our setup and a theoretically ideal Poisson-limited light source measured with the same detection efficiency, (horizontal blue line at a ratio of 1, which corresponds to the shot-noise limit).Data points are the ratio between the variance of our transmittance estimate and the calculated variance of the same measurement for an ideal coherent state, at different levels of sample transmittance (red circles).In each case, the data was normalized relative to the number of input photons incident on the sample.Points above the blue line exhibit subshot-noise performance.Each point corresponds to 13 series of 40 measurements with an integration time of 1 s at a rate of 40000 counts per second.Error bars correspond to the experimental standard deviation of the mean from the different data series.c) Electron beam microscope image of the platinum deposited target used for resolution characterisation.d) Results of imaging the resolution target using correlated photon pairs.Two distinct strips are apparent in the optical image for 3 µm spacing, but for smaller spacings the gap between parallel strips is not discernible.From this, we estimate that the spatial resolution of the SSN microscope to be 3 µm.
forward enables a higher precision ratio than directly using the correlated twin-beam source without optical gating.
The SSN microscope setup is shown in Fig. 2, which is adapted from the photon-pair source and optical switching circuit used in [21,22].Correlated photon pairs are generated by type II collinear SPDC in a 30 mm periodically poled potassium titanyl phosphate crystal (PPKTP) pumped with a continuous wave (CW) 404 nm laser (Toptica, TOPMode 405 HP † ).The wavelengths of the down-converted photons are tuned by adjusting the temperature of the crystal such that the correlated beams are emitted at 818 nm (probe photons) and 798 nm (reference photons).Generated photon pairs are filtered using a long pass and a band pass filter and then deter- † Disclaimer: Commercial equipment, instruments or materials are identified in this report to foster understanding.Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the materials or equipment are necessarily the best available for the purpose.
ministically separated using a polarizing beam splitter (PBS).After being separated, the photons are each coupled into separate single-mode fibers.The probe photon is then detected with a SPAD, heralding the presence of its companion and triggering an optical switch in the reference photon channel.
The optical switch is formed by a free-space Pockels cell modulator (Thorlabs, EO-AM-NR-C1 † ) inside a Sagnac loop as described in [22].When the optical switch has been triggered by detection of a probe photon, the corresponding reference photon is transmitted through the switch, focused onto the sample and then recollimated after the sample using a pair of 10x microscope objectives.Once the reference beam has traversed the sample and has been collimated it is then refocused onto the sensor of a CCD camera (Andor, iDus 416 † ) for detection.The sample is raster scanned in the plane perpendicular to the path of the beam using a pair of electronically driven linear stages.As the intensity of the light passing through the sample is measured, we reconstruct an image of the sample transmittance on a point by point basis.
The sample we imaged is a 3 mm thick AR coated N-BK7 Bright Laser (≈1 mW) SPDC (≈20 fW) Attenuated Laser (≈20 fW) window, in which features were engraved using ion-beammilling.By selectively thinning the anti-reflection coating of the window a low contrast transmittance figure was created.
The ion-beam-milled figures have high transmittance (≈95%) and low contrast (≈2%) which make them ideal to show the practicality of the twin-beam imaging scheme.In other sections of the sample, platinum deposition was used to create higher contrast markers for locating the low contrast sample and to provide targets with which to benchmark resolution.A high power reference image from the same imaging setup (Fig. 3 a), was taken by raster scanning a bright laser set at 1 mW of power (at the target), and measuring the transmitted power with a photodiode power meter.We characterized the resolution of our imaging system using a resolution target consisting of a series of parallel lines with widths ranging from 5 µm to 1 µm (Fig. 2 b).Using the correlated photon pairs to image, we resolve features with widths of down to 3 µm, (Fig. 2 c).
The leakage of unheralded photons through the optical switch was characterized by measuring the Klyshko efficiency η R of the reference beam using SPADs in both channels.We obtained η R = N C /N P = 90(3)%, which corresponds to the proportion of photons exposed to the sample that have been successfully heralded, and therefore the percentage of the photons detected on the camera can be treated as coincidences for our transmittance estimator.
We next characterized the system's performance for precision estimation using a variable neutral-density filter to act as an absorptive sample and estimated the transmittance for different values of attenuation.Using this procedure, we compared our system's precision in estimating transmittance 1/∆ 2 η Exp , (normalised by the mean intensity of the probe beam N PIn ) to that of an ideal shot-noise limited scheme, 1/∆ 2 η Coh , measured with a detector which has the same efficiency as the one used in this experiment (90%).
Here the number of probe photons is estimated from the mean number of detected probe photons, N PDet , corrected by the CCDs dark counts, N DC , the transmittance of the optics after the sample, η Opt , and the transmittance of the sample, η, The results of the system performance characterization are shown in Fig. 2 b.The maximum precision ratio obtained while characterizing our setup without a sample when compared with a direct measurement classical scheme using the same detector efficiency (90%) was Γ =1.76 (9), the reported uncertainty is given by the experimental standard deviation of the mean from the different series of data that were acquired.The precision ratio rises to Γ Dif =2.52(9) when compared to an ideal classical differential scheme.The system achieves sub-shot-noise performance for a transmittance higher than 0.4.Notably this system shows an absolute precision ratio of Γ Abs = 1.58 (9) over an ideal 100% detection efficiency classical experiment.
In Fig. 3 we compare images obtained using correlated photon pairs in our microscopy setup to those obtained using classical illumination.For the classical images (Fig. 3 c), we used a laser beam, attenuated to the same level of intensity as the probe photons input to the microscope from the SPDC source, together with a 50% reflective beamsplitter inserted after the optical switch and used the reflection to monitor power fluctuations on an additional SPAD, implementing a low power differential measurement.Fig. 3 a shows a reference image taken with the same method, but with the laser turned to 1 mW, and using two power meters as the detectors.
The time taken to acquire the low light intensity images in Fig. 3 (b,c) (≈ 1 s per pixel) in our current setup, precluded obtaining a statistically valid number of copies of one image to allow estimation of the variance between them.Instead we imaged eighty copies of a subset of the sample, at the same resolution, and computed the variance across them in a pixel-by-pixel analysis-the data for these images is shown in Fig. 4. When the variance across this sample is compared to the variance that would be obtained with a Poisson distributed illumination source we find a mean precision ratio of Γ =1.54 (26), which lies close to the expected value of Γ =1.73 calculated from the calibration curve in Fig. 2  the a 100% efficient detector, we find Γ Abs = 1.39 (26).As the transmittance is not homogeneous across the sample, the expected advantage in precision varies across the image.Fig. 4 shows a histogram of the transmittance of the eighty images used to translate transmittance into the expected precision ratio using Fig. 2 d, indicated by the top x-axis.We have implemented a 3 µm resolution scanning transmittance microscope that operates with precision that is up to 76% better than the shot-noise limit obtained using a classical (shot noise limited) light source in an equivalent apparatus.We implemented this microscope with a SPDC photon pair source together with single-photon feed-forward in a hybrid CCD-SPAD detection scheme, which requires no active stabilisation or locking of cavity optics.The hybrid detection method we use allows data gathered with a CCD-that by itself is incapable of acquiring photon arrival information-to be used in a coincidence-based estimator of transmittance, at the efficiency achievable in commercially available cooled-CCD sensors.Similar performance might be achieved using a high efficiency superconducting photon-counting detector [25], but this would come at the financial, space and power costs of current implementations of superconducting technology.Our analysis accounts for all photons passing through the sample, for both our measured photon pair experiment and the comparison to an idealised classical experiment.
In comparison with the wide field microscope presented in [5], our imaging system is slow due to the need to raster scan the sample.However, because we emit the probe beam from a single-mode fibre, it is straightforward for the microscope to achieve a high spatial resolution (down to 3 µm) close to the diffraction limit, whilst simultaneously maintaining a high quantum advantage.The resolution could be improved by increasing the numerical aperture of the confocal imaging lenses, in principle to about 0.4µm (half a wavelength) but would come at the cost of increased alignment sensitivity due to reduced depth of focus.
We believe that our scheme could be readily applied to scenarios where it is desired to image a sample with high resolution at the single-photon intensity level.For example, it could be used to obtain a precise spatial characterization of a singlephoton detector like the CCD sensor, or equivalently in singlephoton electrophysiology scenarios [26], where high precision targeted delivery of single-photons, and high resolution spatial information of responses may be required.However, a common critique of twin-beam noise reduction schemes operating within the photon counting regime is that their inherent low intensity is well below any damage threshold of realistic samples-we speculate our approach could be extended to higher intensities by combining optical delay and fast switching with the technique reported in [27], to post-select low noise high intensity beams, increasing the intensity on the sample, whilst maintaining sub-shot-noise performance.

FIG. 3 :
FIG.3: Experimental Results.a) Reference image of the sample, acquired with differential imaging and 1 mW of laser light passing through the sample.b) Noise-reduced image using correlated photon pairs at a rate of 40000 per second passing through the sample, which is equivalent to 20 fW of light passing through the sample.c) A differential image taken with a laser attenuated to the same rate of photons as used in b).Each image is 150 by 75 pixels.The step size of each pixel is 2 µm and the integration time per pixel for a) is 0.1 s and for b) and c) is 1 s.White regions correspond to low transmittance values that are below the scale.

FIG. 4 :
FIG. 4:Multiple scanned image analysis a) A 14x26 pixels subsection of the object imaged in Fig.3in the vicinity of letter "a".The step size between pixels is 2 µm.Each pixel corresponds to the mean of 80 recorded measurements with an integration time of 1 s.The gradient between the left and right side of the image is due to slow drift in the source brightness.b) Histogram of the transmittance and expected precision ratio.The red vertical dashed lines corresponds to the mean transmittance η =0.911 and the bin width is 20.