A Supplementary Notes

Deterministic phase measurements exhibiting super-sensitivity and super-resolution


Phase super-sensitivity is obtained when the sensitivity in a phase measurement goes beyond the quantum shot noise limit, whereas super-resolution is obtained when the interference fringes in an interferometer are narrower than half the input wavelength. Here we show experimentally that these two features can be simultaneously achieved using a relatively simple setup based on Gaussian states and homodyne measurement. Using \num430 photons shared between a coherent- and a squeezed vacuum state, we demonstrate a \num22-fold improvement in the phase resolution while we observe a \num1.7-fold improvement in the sensitivity. In contrast to previous demonstrations of super-resolution and super-sensitivity, this approach is fully deterministic.

© Optica Vol. 5, Issue 1 (2018), pp. 60–64, Optical Society of America. One print or electronic copy may be made for personal use only. Systematic reproduction and distribution, duplication of any material in this paper for a fee or for commercial purposes, or modifications of the content of this paper are prohibited.


exponent-product = ⋅, tight-spacing = false

Current address: ]Centre for Engineered Quantum Systems, School of Mathematics and Physics, University of Queensland, St. Lucia, Queensland 4072, Australia

I Introduction

Quantum interference of light plays a pivotal role in high-precision quantum sensing 1, optical quantum computation 2 and quantum state tomography 3. It is typically understood as two-beam interference which can be observed, for instance, in a Mach–Zehnder interferometer or a double slit experiment. At the output, such interferometers create an oscillatory pattern with a periodicity given by half of the wavelength () of the radiation field, which may be referred to, in analogy to the resolution-benchmark in optical imaging, as the “Rayleigh criterion” for phase measurements. This limit can however be surpassed using different types of states or measurement schemes 4; 5; 6; 7; 8; 9; 10. In particular, measurement schemes which are based on parity detection 11; 9 or approximate parity detection via a phase-space relation 8 are utilised to beat this limit with classical states, i.e. they do not require quantum states 1. The arguably best-known quantum approach to observe a fringe narrowing uses NOON states, . Surpassing the Rayleigh criterion is referred to as super-resolution 13; 14 and is studied in the context of, e.g., optical lithography 5, matter-wave interferometry 15 and radar ranging 16.

In quantifying the performance for applications in quantum sensing and imaging, it is common to evaluate the Fisher information 17 or, equivalently, determine the sensitivity in the interferometric phase measurement. Using coherent states of light the optimal sensitivity is given by where is the mean number of photons of the state 18. This sensitivity constitutes the shot noise limit (SNL). Overcoming the SNL is commonly referred to as super-sensitivity and can be achieved by non-classical states 19; 20; 1; 21. Super-sensitivity based on squeezed states of light has proven to be a powerful and practical way to enhance the sensitivity of gravitational wave detectors 22; 23.

The effects of super-sensitivity and super-resolution can be obtained simultaneously. For example, optical NOON states offer a sensitivity with Heisenberg scaling, , and a phase resolution that scales as corresponding to fringes per half-wavelength. NOON states thus exhibit an equal scaling in the two effects. In contrast, this work shows how resolution and sensitivity are tuneable and can, in fact, compete with each other. Due to the high fragility of NOON states, the complexity in their generation and the commonly probabilistic way of generation, super-sensitivity and super-resolution have been only measured in the coincidence basis and in a highly probabilistic setting 6; 24; 19; 25. It has also been suggested to use two-mode squeezed vacuum states in combination with parity detection to attain the two “super-features” simultaneously 26. However, possibly due to the complications in implementing a parity detection scheme, it has so-far never been achieved experimentally. The complexity associated with the two schemes are due to the involved non-Gaussian states (NOON states) or the non-Gaussian measurements (parity detection). A natural question to ask is whether the same “super-features” can be realised using simple Gaussian operations. Here we answer this question in the affirmative.

We propose and experimentally demonstrate that, by using Gaussian states of light and Gaussian measurements, it is possible to realise a phase measurement which features super-resolution and super-sensitivity simultaneously. Using displaced squeezed states of light in conjunction with homodyne detection followed by a data-windowing technique, we show that the interferometric fringes can be made arbitrarily narrow while at the same time beating the shot noise limit. In stark contrast to the NOON state scheme which, in any practical setting, is highly probabilistic both in preparation and in detection, our approach provides a deterministic demonstration of super-resolution and super-sensitivity.

Ii Materials and methods

An illustration of the basic scheme is shown in Fig. 1. A vacuum squeezed state is combined with a coherent state of light at the entrance to a symmetric Mach–Zehnder interferometer. The Wigner function at the input is given by


where and are the amplitude- and phase-quadratures, is the amplitude of the coherent state, , where denotes the squeezing parameter, and represents the purity of the squeezed state. In the scheme, an amplitude modulated coherent state and a phase-squeezed vacuum state interfere on the first beam splitter of the interferometer. Then the resulting state acquires a relative phase shift , next interferes on the second beam splitter and finally one of the outputs is measured. As we used weak input signals, a homodyne readout scheme was employed. Fig. 1 illustrates the trajectory of the output state in phase space for different phase shifts.

If the interferometer is operated near a dark fringe, i.e. biasing the phase shift such that most of the light exits the second output of the interferometer, the phase-squeezed vacuum state will be detected. Thereby, the shot noise around the bias is suppressed and the phase sensitivity improved. The approach of feeding the commonly unused input mode with a vacuum squeezed state is equal to the proposal by Caves 20 to beat the shot noise limit in phase measurements. However, since the phase response for Caves’ scheme reads , which is an oscillating function with a period equal to , the resolution coincides with the mentioned “Rayleigh criterion” for phase measurements. In the following we show that by implementing a homodyne windowing scheme, the setup yields super-resolution and super-sensitivity.

Figure 1: Scheme of the approach. A coherent state and a vacuum squeezed state are interfered on the first beam splitter. Insets show Wigner functions of the respective states, simulated for and . In one of the resulting modes, a variable phase shifter is placed. At the second beam splitter the modes interfere again, producing the depicted Wigner functions. Eight superimposed distributions illustrate the effect of the phase shift, that is each distribution is separated by . If at the squeezed vacuum state leaves the upper arm, the coherent state exits the lower one. Finally, the state is projected onto the quadrature eigenstate and partitioned by .

The quadrature measurement of the homodyne detector is divided into two bins, set by the ‘bin size’ : If the phase quadrature is measured, we categorise two different results which are associated with the intervals and . We describe such a measurement strategy by the projectors


The measurement observable can thus be written as , where and are the eigenvalues associated with the two measurement outcomes. Now the detector response is found by evaluating which, in the idealised case of , i.e. , and a pure squeezed vacuum state, yields


The full-width-half-maximum (FWHM) of this function follows for , thereby indicating that the interference fringes become narrower as is increasing and thus demonstrating super-resolution. It should be stressed that setting is an idealisation, as it means a projection on an infinitely squeezed state, i.e. even number state. However it points out that the operator is in some sense an approximation of the parity operator 11; 9. Considering instead a realistic setting where and the squeezed state is not pure (), the response function reads


where and


The scaling of the FWHM is preserved for a general value , i.e. FWHM . In Fig. 2a we plot the FWHM-improvement as a function of the squeezing parameter and the bin size . It is clear from this plot that the super-resolution feature only depends weakly on the degree of squeezing, and a similar conclusion is found for the purity of the state. The only critical parameter for attaining high resolution is the mean photon number of the input coherent state. More details and a derivation may be found in the Supplemental Document.

We now turn to the investigation of the sensitivity using the above scheme. The sensitivity can be found using the uncertainty propagation formula,


where , and for our measurement operator it follows


with the notation and . For a specific parameter regime defined by the purity , the bin size and the squeezing parameter , this sensitivity beats the shot noise limit. In Figs. 2b,c we plot the sensitivity relative to the shot noise limited sensitivity as a function of the bin size and the squeezing parameter for two different purities. It is shown in Fig. 2c that it is possible to achieve super-sensitivity in a setting where the squeezed state is impure. In conclusion, both super-sensitivity and super-resolution can be achieved in a practical setup for the parameter space shown in Fig. 2c. Furthermore, sensitivity and resolution features are neither independent nor fixed with respect to each other, but can be varied by the homodyne windowing technique. A discussion of the ultimate sensitivity may be found in the Supplemental Document.



Figure 2: Performance of the protocol for under variation of the bin size and squeezing parameter. a) Improvement of the FWHM compared to the Rayleigh criterion (). The achieved FWHM is extracted numerically from the response of a pure state. The improvement is monotonic in the sense that a smaller bin size always leads to a higher resolution. b) Maximum sensitivity compared to the SNL. The region with negative values describes the parameter space where the SNL is surpassed. The threshold is marked by the bold line. c) Unlike case b), we set the purity .

We proceed by discussing the experimental realisation depicted in Fig. 3. A squeezed vacuum state and a coherent state with a controllable photon number is injected into the input ports of a polarisation based Mach–Zehnder interferometer. The polarisation basis ensures high stability and quality of the interference. Furthermore it allows for simple control of the relative phase shift. The phase shift is varied by a half-wave plate mounted on a remote-controlled rotation stage. One output of the interferometer is measured with a high-efficiency homodyne detector exhibiting an overall quantum efficiency of \SI93\percent, given by \SI99\percent efficiency of the photo diodes and \SI97\percent visibility to the local oscillator (LO). The relative phase of the two input beams of the interferometer as well as the phase of the LO is actively stabilised via real-time feedback circuits, thereby recreating the scheme in Fig. 1 and projecting the output on the quadrature. A detailed description may be found in the Supplemental Document.

Figure 3: Experimental implementation. A vacuum squeezed state, created by parametric down-conversion, and a coherent state, generated via an electro-optic modulator (EOM), are sent into a polarisation-based Mach–Zehnder interferometer (MZI). A quarter-wave plate in combination with a motorised half-wave plate () forms the equivalent phase shift of a MZI where the two modes are spatially separated. The piezo transducers P1 and P2 stabilise the phase between the input states and the local oscillator (LO), respectively. A half-wave plate in front of the last polarising beam splitter is used to balance the photocurrent in the homodyne detector (HD). All cubes represent polarising beam splitters.

Squeezed vacuum is generated by parametric downconversion in a \SI10\milli\meter long periodically-poled KTP crystal embedded in a \SI23.5\milli\meter long cavity comprising a piezo-actuated curved mirror and a plane mirror integrated with end-facet of the crystal. A Pound–Drever–Hall (PDH) scheme is adopted to stabilise the cavity resonance. The downconversion process is pumped by a \SI45\milli\watt continuous-wave laser beam operating at \SI532\nano\meter, such that squeezed light is produced at \SI1064\nano\meter. To stabilise the pump phase, the radio-frequency signal used also for cavity stabilisation is down-mixed with a phase shift of \SI90\degree. Using a \SI5\milli\watt local oscillator, we observe \SI6.5 ±.1\decibel shot noise suppression at \SI5\mega\hertz sideband frequency, while the anti-squeezed quadrature is \SI11.3 ±.1\decibel above shot noise. The squeezed state parameters read, on average, and . A complete characterisation of the squeezed light source is presented in the Supplemental Document.

The coherent input state is produced by an electro-optical modulator (EOM) at a sideband-frequency of \SI5\mega\hertz. The chosen frequency ensures the creation of a coherent state far from low-frequency technical noise and with an amplitude that is conveniently controlled by the modulation depth of the EOM.

To measure the interferometer’s output state at \SI5\mega\hertz, the electronic output of the homodyne detector is down-mixed at this frequency, subsequently low-pass filtered at \SI100\kilo\hertz and then digitised with \SI14bit resolution. For each phase setting, \num[retain-unity-mantissa = false]1e6 samples are acquired at a sampling rate of \SI0.5\mega\hertz. The data is recorded on a computer for post-processing which includes the dichotomic windowing strategy given by (2) in which we set the bin size . After dividing the data according to , we calculate as well as the standard deviation for each phase setting from the data. Finally, is computed according to (6). The term in (6) is extracted directly from the data. Instead of calculating the derivative of also directly, it is estimated from the theoretical model of fitted to the data. This approach is chosen to increase the confidence in the computation of and a comparison between this and a direct evaluation is shown in the Supplemental Document. In the panels on the right of Fig. 4, is shown in comparison to the theoretical model given by (7).

Iii Results and discussion

The results for a mean photon number of \num33.6 and \num430, with a mean of \num2.8 photons contained in the squeezed state and an overall efficiency of circa \SI84\percent, are shown in Fig. 4 and compared with theoretical predictions. The latter is denoted by solid lines. It is clear from the plots that the scheme exhibits super-resolution as well as super-sensitivity for certain phase intervals. We expect that the resolution and sensitivity improves as the mean photon number is increased. This expectation is confirmed in Fig. 5 where the measurement of these two features for increasing photon numbers in the coherent state is depicted. Specifically, at we obtain a \num22-fold improvement in the phase resolution compared to a standard interferometer and a \num1.7-fold improvement in the sensitivity relative to the shot noise limit.



Figure 4: Results achieved for an input state with a) () and b) (). Left: The fringe after applying the dichotomy operator . A dashed line follows the fringe of a standard interferometer. Its FWHM is a) \num5.7 and b) \num22.2 times larger compared to our result. Right: The sensitivity derived from the experimental data. In a range of about , the SNL was surpassed by a factor of a) \num1.5 and b) \num1.7. The uncertainty of each data point is well within the ’’ symbol. We attribute the symmetric deviations at the wings to a systematic anomaly in the set phase-shift controlled by the HWP.


Figure 5: Summary of experimental results. The solid orange lines are theoretical predictions derived from the measured squeezing parameters and displacement amplitude. Each cross symbolises a measurement run. The uncertainties are much smaller than the symbol size. a) The FWHM under variation of the total average photon number of the input state. It always beats the Rayleigh criterion of . Comparing the theoretically predicted FWHM proves a stable performance of the setup. b) A comparison to four sensitivity limits. As for the resolution, the theoretical prediction affirms our experimental results. The SNL was outperformed throughout the experiment at a scaling of ; the Heisenberg scaling of is however not attainable. Using no windowing (), a better sensitivity and a scaling of can be achieved, however it comes at the cost of super-resolution. The ultimate bound of our protocol follows , assuming no losses and restrictions on the squeezing degree.

It is interesting to compare these results with a scheme exploiting pure optical NOON states which exhibit super-resolution and -sensitivity at the same photon-number scaling. Using such states, a similar improvement in resolution and sensitivity would require a \num23-photon and a \num3-photon NOON state, respectively. Importantly, this only holds for a lossless scenario. As of today, an optical \num5-photon NOON state has been produced which in principle will yield a \num5-fold improvement in resolution and a \num2.2-fold improvement in sensitivity 27. However, this realisation is intrinsically probabilistic and thus does not exhibit super-sensitivity in a deterministic setting. To the best of our knowledge, we found that the presented results constitute the first demonstration of super-resolution and super-sensitivity in a deterministic setting.

In summary, we proposed and experimentally demonstrated a simple approach to the simultaneous attainment of phase super-resolution and phase super-sensitivity. The approach is based on Gaussian squeezed states and Gaussian homodyne measurement followed by a windowing strategy, which is in stark contrast to previously proposed schemes realised with impractical and fragile NOON states, or high-efficiency parity detection. Our work is of fundamental interest as it highlights the fact that the observation of super-resolution is not a special quantum effect associated with non-Gaussian quantum states 6 or non-Gaussian measurements 7. In conclusion, we find that the actual quantum feature – that is super-sensitivity – may co-exist with the super-resolution feature without using advanced non-Gaussian states or non-Gaussian measurements. Assuming that the measurement’s figure of merit is phase sensitivity, we can not find an advantage in exploiting super-resolution in a Gaussian-noise governed context. Furthermore, we present the trade-off between resolution and sensitivity for the first time and show that significant super-resolution can be achieved at the cost of negligible increase of sensitivity at the scale of a fraction of SNL. This holds also in the presence of loss and classical Gaussian noise (discussed in the Supplemental Document). Our result sets a benchmark to evaluate super-resolving strategies, particularly under realistic imperfect conditions.

Funding Information

The work was supported by the Lundbeck Foundation and the Danish Council for Independent Research (Sapere Aude 4184-00338B and 0602-01686B). M.J. acknowledges support from the Czech Science Agency (project GB14-36681G).

Appendix A Supplementary Notes

a.1 Squeezed state homodyne tomography

We start by considering the Wigner function of a mixed squeezed vacuum state as


where denotes the squeezing parameter and the purity. The squeezing parameter is related to  via .

Terming the axis as the phase quadrature, (8) describes a phase squeezed state. Its variances are


such that a homodyne tomography parametrised by the phase may be modelled by


In the experiment, we scanned the phase via a piezo-transducer driven by a triangle signal. To account for the nonlinearity of the piezo-transducer, the model function


was applied to describe the measured variance. The coefficients provide an approximate description of the piezo’s nonlinear behaviour and concatenated step functions recreate the phase intervals caused when the driving signal’s slope turns. From this model function,  and  have been extracted for a given optical pump power. Before doing so, the recordings were corrected for dark noise and biased with respect to shot noise. A conversion to a decibel scale was done by and for squeezing and antisqueezing, respectively.

a.2 Resolution and sensitivity expression of the protocol

To arrive at the expression for the response function , and in turn the resolution and the sensitivity , we start by propagating the Wigner function (Eq. (1) in the main text)


through a Mach–Zehnder interferometer. According to the main text, and denote the purity and squeezing parameter, respectively, while the coherent state amplitude . The field quadratures are written as and and their index represents a certain mode. The propagation may be described by a combination of a beam splitter, a phase shift and a second beam splitter transformation by


Doing so yields the output modes


where the output quadratures are labelled with a prime. Next, the Wigner function for the one of the output modes is recovered by tracing out the other one, e.g.


where , . To arrive at Eq. (3) from the main text which represents the idealised case of , i.e.  and , one performs the integral


and norms it with such that . To recover Eq. (4) which we used to process the experimental data, one solves the integral


where .

a.3 Ultimate sensitivity

The following arguments can be derived from Eq. (7) in the main text.

First, we note that the resolution of our protocol improves with increasing the coherent state amplitude . Thus the user aims for the highest laser power which does not corrupt the detection (or the sample). Then for each coherent state amplitude, the optimum squeezing , the optimum bin size and the phase with best sensitivity should be found. However, the optimum squeezing required is of about \SI13\decibel already for and increases quickly. The ultimate sensitivity, given by the maximal Fisher information, then scales as , i.e. , in the limit of large .

For practical purposes, the optimum can not be reached due to the high degree of squeezing. The best sensitivity using our resources is reached when which coincides with the limit for squeezed states and standard homodyne detection.

In Fig. 6a, various sensitivities are plotted. The figure shows that the protocol performs quite well with respect to the efficiency, discussed later, of our setup. Furthermore, it shows that the influence of the bin size , if set to , is small even when the goal is to reach the ultimate bound.

Fig. 6b illustrates the effect of losses for our protocol. The scenario for this simulation is: given the degree of squeezing ( or \SI6.5\decibel) used for the measurements, what is the effect of losses and the coherent state amplitude . In the Wigner function Eq. (1) (main text) we have treated losses in terms of purity which is very convenient when performing state tomography. For a treatment of losses, it is thus necessary to convert to . This can be achieved by comparing Eq. (1) to a Wigner function of a squeezed state which suffered losses via a beam splitter transformation. Once the ancillary modes introduced by the beam splitter are traced out, one may compare the variance of the purity versus the loss based approach and find the mapping


These expressions substitute the original and parameters in Eq. (7).



Figure 6: a) Comparison of for different states under variation of the coherent state amplitude . The SNL line marks the shot noise limit, given by the number of photons in the coherent state plus the (fixed) number of photons in the squeezed state we used in the experiment. Next, follows the sensitivity achieved in our experiment. To see how much better a pure squeezed state would have performed, was simulated. A pure state provides, in our configuration, a negligible enhancement. Next, shows the sensitivity under the assumption that and that the squeezing parameter is freely tunable. The bin size is fixed to . Finally, shows the ultimate sensitivity which is achievable when, in addition to and an optimised , also is optimised. b) A study of the sensitivity improvement when varying the efficiency and . For this simulation, we assume to use a fixed amount of photons in the squeezed state and a variable amount of photons in the coherent state. The “squeezing” photons are given by and in our measurement. The total number of photons is then used to calculate the shot noise limit. Negative-valued regions in the contour plot (above the thick contour line) indicate that the shot noise limit is surpassed. We find that losses of about \SI50\percent restrict the measurement from being super-sensitive. However, even at lower efficiencies, the feature of super-resolution is still maintained.

a.4 Impact of detector noise on sensitivity and resolution

So far we have studied the influence of losses and impure squeezing and considered them as effects that occur during the propagation through the interferometer.

Here we simulate the effect of technical detector noise which enters the mode after the output port of the interferometer. Hence, the thermal photons added by this process do not contribute to the shot noise limit as they do not bear any information about the phase shift. This study is especially important when the detection is the limiting factor in the system. To simulate the impact of detector noise, the interferometer’s output mode was “thermalised” by convoluted it with a Gaussian distribution of unity width. The results are shown in Fig. 7. The ratios of the sensitivity with versus without added noise are summarised in table 1.

First, and similar to the case of losses, our protocol is able maintain super-resolution despite realistic imperfections. Second, we note that super-resolution does not rely on squeezing, i.e. it is a classical effect. In fact, it can be seen that an increased resolution potentially decreases the sensitivity: to highlight this finding in this noise context, we simulated the response and sensitivity of our protocol with added noise, first with and second where was chosen to optimise the sensitivity (denoted in the figures). The case of yields a higher sensitivity at the cost of resolution.

2.3 4.0 3.2 2.0 3.3
2.2 3.8 3.0 1.8
Table 1: The relative sensitivity-degradation for the schemes shown in Fig. 7. The higher the number, the stronger the impact of detector noise.

In quantum enhanced sensing where sensitivity is the figure of merit and in face of Gaussian noise, we find no advantage in using super-resolution. The treatment of non-Gaussian detector effects such as thresholds or the effect of analogue-to-digital converters in the signal chain are question left open for future projects.



Figure 7: Interferometer response (top panels) and phase sensitivity (bottom panels) for different detection schemes and input states. Comparison for a coherent state amplitude of a_) and b_) . Three scenarios are simulated, each with and without added Gaussian detector noise: no squeezing (), squeezing () and direct intensity detection (). In case of , the bin size has been set to and to which optimises the sensitivity. Furthermore, we assumed a pure squeezed state. Panels a_) also show the scenario of using intensity detection and squeezing, which is a technique suggested by Caves (28).

Appendix B Supplementary Methods

b.1 Basic experimental setup

A schematic representation of the basic experimental setup for the generation of squeezed light is shown in Fig. 8. The setup was powered by a continuous-wave Nd:YAG laser (Innolight GmbH Diabolo) providing \SI1064\nano\meter and \SI532\nano\meter radiation, whereas the latter was produced by second harmonic generation from the fundamental. The second harmonic field was employed as pump for squeezed light generation. For increased mode matching efficiency, the pump light was filtered by means of a triangular-shaped travelling-wave mode cleaning cavity (MCC) prior to coupling into the squeezing cavity. The MCC was stabilised using a standard Pound–Drever–Hall (PDH) scheme exploiting the internal phase modulation of the laser at \SI12\mega\hertz for error signal generation.

The fundamental beam was split into two, one serving as steering beam in the generation of squeezed vacuum state while the second and most intense part was used as local oscillator for homodyne detection and coherent state generation.

Figure 8: Simplified setup to generate the local oscillator, coherent- and squeezed state. MCC: Mode cleaning cavity. DM: Dichroic mirror. EOM: Electro-optical modulator. ESA: Electronic spectrum analyser. PID: Servomechanism.

In the path of the steering beam, an electro-optic modulator was placed to generate phase modulation sidebands for PDH stabilisation of the squeezer cavity. A second feedback loop actuating the phase of the pump beam was used for locking the squeezed light source to deamplification or amplification, i.e. amplitude- or phase-squeezing, respectively.

b.2 Characterisation of squeezed light source

The squeezed-light source consisted of a linear Fabry–Pérot resonator enclosing a \SI10\milli\meter periodically poled potassium titanyl phosphate (ppKTP) crystal with two flat end-facets. One end-facet was high-reflective coated for both wavelengths, the fundamental at \SI1064\nano\meter and the pump at \SI532\nano\meter, serving as end mirror for the resonator. The other end-facet was anti-reflective coated for both wavelengths. The coupling mirror was attached to a piezo transducer and had a reflectivity of \SI90\percent for \SI1064\nano\meter and \SI20\percent for \SI532\nano\meter. The mirror was polished to a radius of curvature of \SI20\milli\meter and was placed \SI13\milli\meter from the crystal. This yielded a full-width-half-maximum cavity bandwidth of about \SI80\mega\hertz. To achieve phase matching, the non-linear crystal was attached to a Peltier element. A phase matching temperature was reached at \SI36.30\celsius.

The cavity was locked by a PDH phase modulation-demodulation technique at \SI37.22\mega\hertz using a \SI550\micro\watt steering beam launched into the cavity from the high-reflective mirror. The pump phase was locked using the same phase modulation as for the cavity lock, but with a demodulation phase of \SI90\degree with respect to the cavity lock error signal. To achieve amplitude squeezing the pump phase was locked to deamplification of the steering beam.

Prior to sending the squeezed beam into the actual interferometric setup, a characterisation of the squeezing degree was performed. A homodyne detector was set up close to the output of the squeezed light source to lower optical losses. For exemplification, Fig. 9 shows the variance at a pump power of \SI2\milli\watt and \SI92\milli\watt.

Figure 9: Squeezing spectra at \SI2\milli\watt (left) and \SI92\milli\watt optical pump power. \num20000 samples were recorded over \SI200\milli\second at a resolution bandwidth of \SI270\kilo\hertz and a video bandwidth of \SI680\hertz with a Agilent N9000A CXA set to \SI5\mega\hertz. The local oscillator (LO) phase was varied by means of a piezo actuated mirror. The light yellow line on top of the experimental data in green was fitted to (11). This data was taken downstream the squeezed light source, i.e. it bypassed the interferometer used for phase sensing.

Theoretically, the variance of (anti-)squeezed light versus pump power follows (33)


where denotes the threshold power, the total efficiency of the system, the cavity decay rate and the probing frequency. The upper sign shall be chosen to describe the antisqueezing variance. (19) was used to fit the recorded variance of the homodyne tomography and estimate , and . Inserting the standard deviations of the fit parameters provided the prediction bounds shown in Fig. 10. Estimated values are summarised in table 2.

Figure 10: Scaling of (anti-)squeezing variance with respect to pump power. Two measurement sets are shown. Open circles denote anti-squeezing. The filled area represents the lower and upper prediction band of the fit with a confidence level of \SI99\percent. Parameter uncertainties are smaller than the denoting symbols. At \SI92\milli\watt pump power, a shot noise suppression of \SI12.45 ±.05\decibel was achieved.
\SI113.73 ±2.89\milli\watt
\SI94.230 ±0.251\percent
\SI519.61 ±76.87\radian\mega\hertz
Table 2: Estimates with standard uncertainties of the model parameters in (19). Values are based on the two recorded sets shown in Fig. 10.

b.3 Use and efficiency of squeezed light in the interferometer

In the actual interferometer used for this experiment, the settings in table 3 apply. As stated, the pump power was set to \SI45\milli\watt, i.e. only a fraction of the available degree of squeezing was employed. This choice was made to avoid saturation of the homodyne detector, as to resolve a higher degree of squeezing a stronger local oscillator is required. When the phase in the interferometer is set such that the squeezed state is detected, this causes no issues. However as the phase of the interferometer is swept, a much stronger coherent state will be detected. At high local oscillator power, this easily saturates the electronics.

To characterise the efficiency of squeezed light detection downstream the interferometer, a homodyne tomography similar to the one discussed above was performed. Theoretically, the loss-reduced degree of squeezing is given by


such that we estimated using the values from table 3.

Set Measured
Parameter Value Parameter Value
Crystal temperature \SI36.6\celsius Dark noise clearance \SI18\decibel
Pump power \SI45\milli\watt \SI99\percent
Steering beam power \SI400\micro\watt \SI97\percent
Local oscillator power \SI5\milli\watt \SI-6.5 ±.1\decibel
Down-mix frequency \SI5\mega\hertz \SI11.3 ±.1\decibel
EOM driving frequency \SI5\mega\hertz \num.582 ±.001
Table 3: Parameter values set or measured during the experiment. The visibility was measured for the inference between the local oscillator (LO) and beam for the coherent state () and squeezed state (). The degree of (anti-)squeezing and purity were characterised via a fit to a homodyne tomography.

b.4 Interferometer setup

To perform the quantum enhanced phase measurement scheme presented in the main text, a Mach–Zehnder interferometer was built. Instead of using the standard configuration, we implemented a polarisation based Mach–Zehnder interferometer where the two spatial modes are substituted by orthogonal polarisation modes. This means that only one spatial mode comprises the to-be-interfered beams and inherently increases the mechanical stability, as demonstrated experimentally in similar setups (31).

Fig. 11 illustrates a scaled version of the optical setup. The electronic components involved to stabilise the measurement are summarised in Fig. 12.

Figure 11: To-scale schematic of the interferometer. All cubes represent polarising beam splitters. The wave plate labelled by was mounted in a remote-controlled rotation stage. Detector D1 monitored the power and the polarisation of the local oscillator. The interference between a fraction of the coherent- and the squeezed beam was detected at D2. Its signal was part of a feedback loop to establish a stable phase stabilisation between the the squeezed- and the coherent state. PM-SMF: Polarisation maintaining single-mode fibre. HD: Homodyne detector. P: Mirror mounted on piezo-transducer. EOM: Electro-optical modulator. D: Photodetector. Red path: Squeezed state. Green path: Coherent state. Blue path: Local oscillator.
Figure 12: Diagram of the employed electronics. All shown devices were built in-house. External frequencies were generated by a digital synthesiser based on a Analog Devices AD9959 chip. The photocurrent from D2 was mixed with an electronic local oscillator at the same frequency as used for the phase modulation of the steering beam. Thereby we could reference the coherent state and the steering beam, carrying the squeezed state, to the beam providing the local oscillator. This holds as the steering beam and local oscillator were derived from the same source. The AC-coupled output of the homodyne photocurrent was split in two: The first half was mixed with a sinusoidal signal at \SI5\mega\hertz and low-pass filtered. This scheme mimics a spectrum analyser set to \SI5\mega\hertz. After down-mixing the signal, it was acquired by a 14-bit data acquisition (DAQ) card. To match the bandwidth of the card, a low-pass filter of \SI100\kilo\hertz was connected between the output of the homodyne detector and the DAQ card. The second half served for locking the homodyne detection on the amplitude quadrature. A stable measurement of the amplitude quadrature was ensured by a servomechanism. HV amp: High voltage amplifier.

It is convenient to separate the description into three parts: The beam representing the squeezed state (drawn in red in Fig. 11), the one representing the coherent state (green) and finally the local oscillator (blue).

The local oscillator and the coherent state were derived from the infrared beam which was coupled into a polarisation maintaining single-mode fibre (Thorlabs P3-1064PM-FC-2). A fibre-based beam splitter sent \SI90\percent of the power directly to a fibre collimator (Thorlabs TC12APC-1064). To control the optical power and polarisation properties, photodetector D1 was placed after a removable mirror and a half wave plate / polarising beam splitter (HWP / PBS) combination. The interference visibility between local oscillator and signal beam, which combined the squeezed- and the coherent state, was optimised by means of a telescope. For the measurement run, the visibility read \SI97\percent, mainly limited by the large beam waist of \SI3.4\milli\meter of the signal beam due to a propagation length of circa \SI7\meter measured from the squeezer cavity (32). The propagation length of \SI7\meter is due to the fact that the squeezing source and the actual experiment were built on two different optical tables. Piezo P2 actuated a mirror to control the local oscillator’s phase. A HWP set the local oscillator power before combining it spatially with the signal beam. The next HWP in combination with a PBS rotated both the local oscillator’s and the signal beam’s polarisation state, thereby splitting up the two beams into equal halves. Finally, a homodyne detector converted the interference signal into a photocurrent.

Following the proposed input scheme (Fig. 1 in the main text), the amplitude quadrature of the signal beam had to be detected. Hence, the local oscillator and signal beam had to be locked at a fringe maximum or minimum. For locking at this point, we implemented an electronic mixer to detected the interference at \SI37.22\mega\hertz, which was the modulation frequency of the steering beam’s phase. This technique is usually referred to as ‘AC lock’ (29).

The smaller fraction of light split by the fibre-based PBS (Thorlabs PBC1064SM-APC) travelled through an electro-optical modulator (EOM) driven at \SI5\mega\hertz before coupled out by a fibre collimator. The EOM (Photline NIR-MPX-LN-0.1-P-P-FA-FA) modulated the phase of the field, such that we prepared, at the given frequency, a coherent state in the phase quadrature. This preparation scheme contrasts the concept illustrated in the main text, where an amplitude-displaced coherent state enters the interferometer. In fact, as we prepared an amplitude squeezed state, also the squeezed state does not correspond to the configuration in the main text. Hence, the phase of both input states had to be rotated by \SI90\degree to match the illustrated configuration. As outlined in the previous paragraph, a rotation of the relative phase can be achieved by implementing an AC lock. Equal to the homodyne detection lock, the reference signal was provided by the phase modulation at \SI37.22\mega\hertz of the steering beam. The interference visibility with the local oscillator measured \SI99\percent at the homodyne detector. Depending on the amplitude of the applied modulation signal, the coherent state’s average photon number was controlled. To calibrate the voltage-to-photon-number conversion, a homodyne tomography was performed.

Next, we turn to the beam transferring the squeezed state, which we term “squeezed beam”. As mentioned above, the squeezer was operated at amplitude squeezing, i.e. de-amplification. This operation is beneficial to the noise features of the steering beam, as an operation at amplification leads to an increase of technical noise from the laser source. Furthermore, the decreased power helped to prevent saturation of the detector electronics. In front of the squeezer cavity, the optical power of the steering beam measured \SI400\micro\watt. The pump beam at \SI532\nano\meter had a power of \SI45\milli\watt. To split a small fraction from the squeezed beam for interference with the coherent beam at detector D2, a HWP was placed in front of the PBS which combined the two beams spatially. After this PBS, the squeezed beam was polarised, orthogonal to the polarisation of the coherent beam. However, both beams share the same spatial mode, which was guaranteed by a visibility of \SI99\percent at detector D2. In this situation, the two polarisation modes constitute the interferometer arms. To delay one mode with respect to the other, i.e. to create the phase shift, a HWP was used. Mounted in a remote-controlled rotation stage (Thorlabs K10CR1/M), we could control the experiment from a PC. To mimic a common Mach–Zehnder interferometer, the relative phase shifts caused by reflections and the employed locking technique had to be considered. Analysing the phase shifts by means of the Jones formalism (30) showed that an additional quarter wave plate was required to mimic the spatial Mach–Zehnder interferometer. Its polarisation axis should be oriented parallel to either of the beams.

The important experimental parameters are summarised in table 3.

b.5 Data acquisition and processing

The homodyne detector was a direct photocurrent subtraction design equipped with photodiodes with a quantum efficiency of . The detector circuit featured three outputs: A DC output with a bandwidth of \SI330\kilo\hertz, an AC output with a high-pass filter of \SI1\mega\hertz and an output of a signal created by mixing the AC signal with an electronic local oscillator. The latter signal was low-pass filtered at \SI500\kilo\hertz to provide a down-mixed signal for data acquisition. In this way, the data was recorded on a PC without using an electronic spectrum analyser. Since the coherent state was prepared by means of the EOM driven at a frequency of \SI5\mega\hertz, the electronic local oscillator was set to the same frequency. To minimise the required electrical power for the mixing process, the input for the electronic local oscillator was equipped with a resonant filter.

A digital oscilloscope with a PCIe interface (GaGe CSE8384) sampled the data from the homodyne detector. It featured a bandwidth of \SI100\kilo\hertz at a sampling rate of \SI500\kilo\hertz and resolution of \SI14bit. Per measurement point, \num[retain-unity-mantissa = false]1e6 samples were acquired. A measurement point was defined by the setting of the motorised HWP. Centred about the null phase, \num57 points were recorded every \SI.9\degree (\SI15.7\milli\radian). Beyond this central region, \num40 points were taken with an increment of \SI2\degree. These points were used to validate the response of the stage actuating the HWP. The overall procedure was repeated for seven different coherent state amplitudes.

The actual data processing was done as follows:

  1. Characterise the squeezed- and coherent state individually via homodyne tomography. Thus we account for all photons (potentially) interacting with the sample to faithfully calibrate the sensitivity . The photons contributed by the local oscillator are not taken into account, as they do not enter the interferometer.

  2. Null the data for a correct zero phase. This step is important for further processing, as it symmetrises the data such that the \SI0\degree setting implies detecting a squeezed state.

  3. Choose a value for and sort the data of each measurement point according to


    with and . As mentioned in the main text, was set to . This step yields .

  4. Fit the function to the experimental results.

  5. Calculate the variance at each data point. With an analytic expression of , the derivative in the sensitivity expression


    may be calculated. Given the variance and derivative at each data point, the sensitivity is computed. From the fit to from the previous step, a comparison to an analytic expression can be derived. To justify the validity of this approach, Fig. 13 compares the sensitivity found by evaluating (22) as discussed, i.e. with a algebraic evaluation of the derivative, and an entirely numerical evaluation.

    Figure 13: Comparison of  where the derivative in (22) is evaluated by inserting the experimental data into a fit (square symbols) and where the derivative is evaluated numerically. The phase points of the “all numerical” approach are shifted, because the derivative is approximated at . On the other hand, using an algebraic expression for allows for an evaluation at . The left panel shows the sensitivity for , the right one shows it for (Fig. 4b in the main text).

Supplementary References


  1. A simple approach for increasing the fringe resolution by classical means is a multi-pass interferometer. As the mentioned techniques can be mapped to multi-pass configurations, we focus on the single-pass configuration.


  1. V. Giovannetti, S. Lloyd,  and L. Maccone, Science 306, 1330 (2004).
  2. P. Kok, W. J. Munro, K. Nemoto, T. C. Ralph, J. P. Dowling,  and G. J. Milburn, Rev. Mod. Phys. 79, 135 (2007).
  3. A. I. Lvovsky and M. G. Raymer, Rev. Mod. Phys. 81, 299 (2009).
  4. J. G. Rarity, P. R. Tapster, E. Jakeman, T. Larchuk, R. A. Campos, M. C. Teich,  and B. E. A. Saleh, Phys. Rev. Lett. 65, 1348 (1990).
  5. A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, C. P. Williams,  and J. P. Dowling, Phys. Rev. Lett. 85, 2733 (2000).
  6. M. W. Mitchell, J. S. Lundeen,  and A. M. Steinberg, Nature 429, 161 (2004).
  7. K. J. Resch, K. L. Pregnell, R. Prevedel, A. Gilchrist, G. J. Pryde, J. L. O’Brien,  and A. G. White, Phys. Rev. Lett. 98, 223601 (2007).
  8. E. Distante, M. Ježek,  and U. L. Andersen, Phys. Rev. Lett. 111, 033603 (2013).
  9. Y. Gao, P. M. Anisimov, C. F. Wildfeuer, J. Luine, H. Lee,  and J. P. Dowling, J. Opt. Soc. Am. B 27, A170 (2010).
  10. L. Cohen, D. Istrati, L. Dovrat,  and H. S. Eisenberg, Opt. Express 22, 11945 (2014).
  11. C. C. Gerry, Phys. Rev. A 61, 043811 (2000).
  12. A simple approach for increasing the fringe resolution by classical means is a multi-pass interferometer. As the mentioned techniques can be mapped to multi-pass configurations, we focus on the single-pass configuration.
  13. J. Jacobson, G. Björk, I. Chuang,  and Y. Yamamoto, Phys. Rev. Lett. 74, 4835 (1995).
  14. E. J. S. Fonseca, C. H. Monken,  and S. Pádua, Phys. Rev. Lett. 82, 2868 (1999).
  15. J. P. Dowling, Phys. Rev. A 57, 4736 (1998).
  16. K. Jiang, H. Lee, C. C. Gerry,  and J. P. Dowling, J. Appl. Phys. 114, 193102 (2013).
  17. S. L. Braunstein and C. M. Caves, Phys. Rev. Lett. 72, 3439 (1994).
  18. A. Monras, Phys. Rev. A 73, 033821 (2006).
  19. T. Nagata, R. Okamoto, J. L. O’Brien, K. Sasaki,  and S. Takeuchi, Science 316, 726 (2007).
  20. C. M. Caves, Phys. Rev. D 23, 1693 (1981).
  21. J. C. F. Matthews, X.-Q. Zhou, H. Cable, P. J. Shadbolt, D. J. Saunders, G. A. Durkin, G. J. Pryde,  and J. L. O’Brien, npj Quantum Information 2, 16023 (2016).
  22. H. Grote, K. Danzmann, K. L. Dooley, R. Schnabel, J. Slutsky,  and H. Vahlbruch, Phys. Rev. Lett. 110, 181101 (2013).
  23. LIGO Scientific Collaboration, Nat. Photon. 7, 613 (2013).
  24. D. Bouwmeester, Nature 429, 139 (2004).
  25. Y. Israel, S. Rosen,  and Y. Silberberg, Phys. Rev. Lett. 112, 103604 (2014).
  26. P. M. Anisimov, G. M. Raterman, A. Chiruvelli, W. N. Plick, S. D. Huver, H. Lee,  and J. P. Dowling, Phys. Rev. Lett. 104, 103602 (2010).
  27. I. Afek, O. Ambar,  and Y. Silberberg, Science 328, 879 (2010).
  28. C. M. Caves. Quantum-mechanical noise in an interferometer. Phys. Rev. D, 23:1693–1708, 4 1981.
  29. A. Freise and K. A. Strain. Interferometer Techniques for Gravitational-Wave Detection. Living Rev. Relativity, 13, 2010.
  30. E. Hecht. Optics. Addison-Wesley Longman, 4 edition, 2002.
  31. M. Mičuda, E. Doláková, I. Straka, M. Miková, M. Dušek, J. Fiurášek, and M. Ježek. Highly stable polarization independent Mach-Zehnder interferometer. Rev. Sci. Instrum., 85(8), 2014.
  32. A large beam waist is beneficial to decrease the divergence of a beam, however it renders the localisation of the waist difficult. The less precise the waist and its localisation is known, the more adjustment is required to match such a beam to a second beam. In principle, the use of a telescope in the signal beam path could have decreased the waist to a smaller size, but to reduce the optical losses we tried to maintain a low count of optics.
  33. Ling-An Wu, M. Xiao, and H. J. Kimble. Squeezed states of light from an optical parametric oscillator. J. Opt. Soc. Am. B, 4(10):1465–1475, 10 1987.
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.