This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Fundamental limits to the generation of highly displaced bright squeezed light using linear optics and parametric amplifiers

Steve M. Young [email protected] Sandia National Laboratories, Albuquerque, New Mexico 87123    Daniel Soh [email protected] Wyant College of Optical Sciences, The University of Arizona, Tucson, Arizona 85721, USA
Abstract

High quality squeezed light is an important resource for a variety of applications. Multiple methods for generating squeezed light are known, having been demonstrated theoretically and experimentally. However, the effectiveness of these methods – in particular, the inherent limitations to the signals that can be produced – has received little consideration. Here we present a comparative theoretical analysis for generating a highly-displaced high-brightness squeezed light from a linear optical method – a beam-splitter mixing a squeezed vacuum and a strong coherent state – and parametric amplification methods including an optical parametric oscillator, an optical parametric amplifier, and a dissipative optomechanical squeezer seeded with coherent states. We show that the quality of highly-displaced high-brightness squeeze states that can be generated using these methods is limited on a fundamental level by the physical mechanism utilized; across all methods there are significant tradeoffs between brightness, squeezing, and overall uncertainty. We explore the nature and extent of these tradeoffs specific to each mechanism and identify the optimal operation modes for each, and provide an argument for why this type of tradeoff is unavoidable for parametric amplifier type squeezers.

I Introduction

Vacuum squeezed light has long been studied and used in various quantum applications [1]. Particular interest lies in a squeezed light with high displacement and high brightness. This kind of quantum state becomes useful for various applications: for example, it can improve the measurement accuracy dramatically when the highly-displaced high-brightness phase-squeezed light is used for measuring optical phase delays [2, 3]. The phase measurement uncertainty is dictated by the phase variance of the state of the light, and highly-displaced high-brightness squeezed light is tremendously advantageous in reducing the phase uncertainty. It has been widely discussed that any metrology and sensing application would require both sensitivity and precision where the former can be improved by squeezed noise below the shot noise while the latter could only be achieved by high power probe light [4]. Therefore, a highly-displaced high-brightness squeezed light will accomplish both the enhanced sensitivity and improved precision. A recent result of a super-resolution imaging using such highly-displaced high-brightness squeezed light attests the critical usage of such a quantum state of light [5]. Alternatively, one can use a high-brightness amplitude-squeezed light to control an optically driven energy transition process [6, 7]. This scenario becomes particularly relevant if the uncertainty of the final energy state after transition depends upon the uncertainty in the driving field’s amplitude [8].

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 1: Three distinct methods to create high-brightness squeezed light: 1 beam-splitter mixing, 1 seeded optical parametric oscillator, 1 seeded optical parametric amplifier,and 1 seeded dissipative optomechanical light squeezing.

Broadly, there are two possible approaches to generating highly-displaced high-brightness squeezed light. The first is to generate a squeezed vacuum state using a given process and displace it through a beam splitter seeded with a strong coherent state of light. The second is to seed the process for generating squeezed light so that it directly produces a highly-displaced squeezed signal rather than a squeezed vacuum. However, in addition to engineering constraints, methods of squeezing face fundamental limits to the quality of the squeezed light that can be generated even for the rudimentary squeezed vacuum states. One can expect that obtaining highly-displaced high-brightness squeezed states will come with additional restrictions.

In this work we compare methods to generate squeezed light with high displacement and high brightness and examine the fundamental limits to the brightness and squeezing that can be obtained. First, we consider the case of beam-splitter mixing between a highly squeezed vacuum, provided by a method of choice, and high brightness coherent light. Second, we consider the direct generation of highly-displaced bright squeezed light through seeding for the cases of optical parametric oscillators, optical parametric amplifiers, and dissipative optomechanical squeezers. We analyze the performance and the limitation of each method and evaluate the trade-offs among the brightness, the level of squeezing, and overall uncertainty of quadrature measurements.

II Beam-splitter mixing method

A straightforward method for creating a light state that has a reduced noise below the standard quantum limit (i.e.,shot noise) is to use a beam splitter and mix a squeezed vacuum and a high-brightness coherent state light (see Fig. 1(a)). The two output modes from a beam splitter are generally entangled, unless the two input modes are coherent states. Therefore, the output that carries a reduced noise is indeed a mixed state after tracing out the other nearly coherent state light. Still, the expected variance of the noise in the main squeezed output would be below the shot noise provided that the mixing ratio is sufficiently small.

The initial state of the two input modes in the tensor Hilbert space is

|ψi=Ss(B)|0sDp()|0p,\left|\psi_{i}\right\rangle=S_{s}(B)\left|0\right\rangle_{s}\otimes D_{p}(\mathcal{E})\left|0\right\rangle_{p}, (1)

where the squeezing operator is given as [9]

Ss(B)=exp(12(Bas2Bas2)),S_{s}(B)=\exp\left(\frac{1}{2}(B^{*}a_{s}^{2}-Ba_{s}^{\dagger 2})\right), (2)

a,aa,a^{\dagger} are the annihilation and the creation operators for the first input mode, and the displacement operator is

Db()=exp(apap).D_{b}(\mathcal{E})=\exp\left(\mathcal{E}a_{p}^{\dagger}-\mathcal{E}^{*}a_{p}\right). (3)

Here, \mathcal{E} is the complex amplitude of the displacement. Additionally, |0s,p\left|0\right\rangle_{s,p} represent a vacuum state in mode s,ps,p, respectively.

The action of a beam splitter is best explained by a beam-splitter Hamiltonian given as

HBS=η(asap+asap),H_{\rm BS}=-\hbar\eta(a_{s}^{\dagger}a_{p}+a_{s}a_{p}^{\dagger}), (4)

where η\eta is the interaction (photon-exchange) strength. As shown in Appendix A, for XX displaced seed states in the limit of high pump displacement, the variance of quadrature measurements can be written as

(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =e2B+(1e2B)α2,\displaystyle=e^{-2B}+\left(1-e^{-2B}\right)\alpha^{2},
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =e2B+(1e2B)α2.\displaystyle=e^{2B}+\left(1-e^{2B}\right)\alpha^{2}.

where α2(=sin2ηt)\alpha^{2}(=\sin^{2}\eta t), where tt is the time duration of beam splitter interaction, is the squared ratio of seed output displacement and pump input displacement. The squeezing, which is in amplitude for B>0B>0 and in phase for B<0B<0, is inversely related to the displacement squared of the output. This result is quite physically intuitive as the two terms on the right-hand side explain the proportional noise addition between a squeezed vacuum and a displaced vacuum according to the beam-splitter ratio ηt\eta t. Enhanced seed amplitude results from the mixing with the pump state, which contributes noise from the pump in addition to displacement from the pump. The overall uncertainty that we define as the product of the standard deviations in the measurements of both quadratures is

ΔXsΔPs=1+4α2(1α2)sinh2(B).\displaystyle\Delta X_{s}\Delta P_{s}=\sqrt{1+4\alpha^{2}\left(1-\alpha^{2}\right)\sinh^{2}\left(B\right)}. (5)

reflecting the mixed state produced by the beam splitter. Thus, in addition to accepting less squeezing as we increase the mixing ratio, we see that we also deviate increasingly from perfect squeezing that maintains the minimum overall uncertainty. Allowing the beams to mix for a longer interaction time increases the output displacement but at the expense of squeezing magnitude and quality.

The consequences of this are shown in Fig.2. First, we emphasize that there is a tradeoff between squeezing and output brightness, as the increase in noise variance is directly proportionate to α2\alpha^{2}. Additionally, the mixing of the pump signal impacts overall uncertainty. If a weak (i.e.,low-displacement), but highly-squeezed output is desired it is optimal to use a highly squeezed vacuum input with low mixing rate (i.e.,small ηt\eta t). However, for brighter output states, using highly-squeezed vacuum input significantly increases overall uncertainty (Eq. (5)); minimizing overall uncertainty requires using a less squeezed vacuum state, which is somewhat counterintuitive. In Fig.2 we plot the maximum squeezing possible for a given α2\alpha^{2}, obtained by varying ηt\eta t and BB constrained by constrained α2\alpha^{2}, allowing for the visualization of both tradeoffs together. Ultimately, there is an unavoidable tradeoff between brightness and squeezing of the output inherent in the beam splitter mechanism.

Refer to caption
Refer to caption
Figure 2: 2 The top panel shows the degradation of output squeezing factor as the relative squared output displacement α2\alpha^{2} is increased for two different squeezing levels of the seed input. The bottom panel shows the accompanying increase in overall uncertainty. 2 The best possible squeezing of the output obtainable with overall uncertainty below a given threshold. The plots are representative both phase and amplitude squeezing, as the system operates symmetrically with respect to the two modes.

III Optical Parametric Oscillator/Amplifier

We next analyze and predict the performance of the second method - a seeded degenerate optical parametric amplifier or oscillator (Fig. 1 (b,c)). In these systems a nonlinear medium couples pump signal to a half-frequency seed signal via the interaction Hamiltonian

HOP\displaystyle H_{\rm OP} =ωasas+2ωapap+i(g/2)(as2apas2ap),\displaystyle=\hbar\omega a_{s}^{\dagger}a_{s}+\hbar 2\omega a_{p}^{\dagger}a_{p}+i\hbar(g/2)(a_{s}^{\dagger 2}a_{p}-a_{s}^{2}a_{p}^{\dagger}), (6)

where apa_{p} and asa_{s} are the field operators for the pump and seed signals, respectively, gg is the field coupling, and ω\omega is the seed signal frequency. At present it is the method of choice for generating squeezed light, with greater than 1010 dB squeeze factors obtainable [10, 11, 12, 13].

Depending on the relative phase of the pump (pp) and seed (ss) signals, the modified seed signal is either amplified or deamplified by the interaction. Here we will take the pump and seed inputs to be XX displaced, with the operating mode depending on the relative sign of the pump displacement with respect to the seed displacement. Additionally, the output experiences phase or amplitude squeezing, respectively. In the case of a parametric amplifier the signals are fed directly through the medium, whereas in the case of a parametric oscillator the medium in placed inside of a cavity, into which the signals enter and interact.

We first consider the case of the parametric oscillator where input and output occur on a single side of the cavity. This was first analyzed in [14] ignoring pump field noise, where it was found that the squeezing depends on the ratio gapκ/2\frac{g\langle a_{p}\rangle}{\kappa/2} and that maximally phase-squeezed output is possible when gapκ/2=1\frac{g\langle a_{p}\rangle}{\kappa/2}=1. However, this result was for the non-seeded case so that the squeezed output state is a vacuum. It was shown later [15] that the presence of a seed signal requires explicit consideration of the noise present in the pump signal and significantly alters the squeezing that can be obtained. In particular, unlike the beam-splitter case the choice of squeezing axis dramatically impacts the squeezing performance, which is significantly impaired by high seed power.

The equations of motions and their solutions can be found in Appendix B. The results are easiest to understand in the limit of small sin\mathcal{E}_{s}^{\rm in} (which is the input field displacement.) Perturbatively expanding around the internal seed displacement sin=0\mathcal{E}_{s}^{\rm in}=0, the relative squared output displacement is

α2=(soutpin)2=(1+𝒞01𝒞0sinpin)2,\displaystyle\alpha^{2}=\left(\frac{\mathcal{E}_{s}^{\rm out}}{\mathcal{E}_{p}^{\rm in}}\right)^{2}=\left(\frac{1+\mathcal{C}_{0}}{1-\mathcal{C}_{0}}\frac{\mathcal{E}_{s}^{\rm in}}{\mathcal{E}_{p}^{\rm in}}\right)^{2},

where \mathcal{E} are the input and output field displacements for the pump and seed, as indicated by the subscripts and superscripts. Here, 𝒞0=gAp0κ/2\mathcal{C}_{0}=\frac{gA_{p}^{0}}{\kappa/2}, where Ap0A_{p}^{0} is the intracavity pump amplitude for sin=0\mathcal{E}_{s}^{\rm in}=0, and κ\kappa is the cavity dissipation rate. The variances of the output quadratures (Eq.(15)) are

(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =(3𝒞022(1+𝒞0)2α2(1+𝒞0)3𝒞022(1+𝒞0)2α2+(1𝒞0))2+4𝒞02(1+𝒞0)2(3𝒞022(1+𝒞0)2α2+(1𝒞0))2α2,\displaystyle=\left(\frac{\frac{3\mathcal{C}_{0}^{2}}{2\left(1+\mathcal{C}_{0}\right)^{2}}\alpha^{2}-\left(1+\mathcal{C}_{0}\right)}{\frac{3\mathcal{C}_{0}^{2}}{2\left(1+\mathcal{C}_{0}\right)^{2}}\alpha^{2}+\left(1-\mathcal{C}_{0}\right)}\right)^{2}+\frac{\frac{4\mathcal{C}_{0}^{2}}{\left(1+\mathcal{C}_{0}\right)^{2}}}{\left(\frac{3\mathcal{C}_{0}^{2}}{2\left(1+\mathcal{C}_{0}\right)^{2}}\alpha^{2}+\left(1-\mathcal{C}_{0}\right)\right)^{2}}\alpha^{2},
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =(𝒞022(1+𝒞0)2α2(1𝒞0)𝒞022(1+𝒞0)2α2+(1+𝒞0))2+4𝒞02(1+𝒞0)2(𝒞022(1+𝒞0)2α2+(1+𝒞0))2α2.\displaystyle=\left(\frac{\frac{\mathcal{C}_{0}^{2}}{2\left(1+\mathcal{C}_{0}\right)^{2}}\alpha^{2}-\left(1-\mathcal{C}_{0}\right)}{\frac{\mathcal{C}_{0}^{2}}{2\left(1+\mathcal{C}_{0}\right)^{2}}\alpha^{2}+\left(1+\mathcal{C}_{0}\right)}\right)^{2}+\frac{\frac{4\mathcal{C}_{0}^{2}}{\left(1+\mathcal{C}_{0}\right)^{2}}}{\left(\frac{\mathcal{C}_{0}^{2}}{2\left(1+\mathcal{C}_{0}\right)^{2}}\alpha^{2}+\left(1+\mathcal{C}_{0}\right)\right)^{2}}\alpha^{2}.

As with the beam splitter case, we have written the quadrature variances in terms of α2\alpha^{2} to highlight the connection to the noise.

The first terms on the RHS can be attributed to noise originating from the seed input, while the second terms are due to the noise from the pump input. We can see that for vacuum output we recover the original expressions in [14] that allow for arbitrarily strong squeezing, while increasing the desired output through a displaced seed input adds noise to both quadratures proportional to α2\alpha^{2}, necessarily compromising the squeezing and overall uncertainty. To see this quantitatively, results calculated using the exact expressions Eqs. (15) and (B) are shown in Fig. 3.

Refer to caption
(a) Phase squeezing
Refer to caption
(b) Amplitude squeezing
Figure 3: Characteristics of the output signal for 3(a) phase-squeezing and 3(b) amplitude-squeezing. In the top panels, the amplitude of the seed output relative to the pump input is plotted as a function of time for different values of 𝒞0=gAp0κ/2\mathcal{C}_{0}=\frac{gA_{p}^{0}}{\kappa/2}. In the center and bottom panels, the corresponding squeezing factor and overall uncertainty of the output are plotted against relative squared output displacement α2\alpha^{2}. For amplitude squeezing, the squeezing and uncertainty curves are cutoff where the relationship to target output power becomes nonmonotonic due to the high seed input required to overcome deamplification.

When p\mathcal{E}_{p} is negative and s\mathcal{E}_{s} is positive, the seed signal is amplified and the output signal is phase squeezed. This case is shown in Fig. 3(a) for different values of gApκ/2\frac{gA_{p}}{\kappa/2}. The amplification of the signal aids in the generation of a bright signal for a given level of squeezing; the parameter regime that results in the strongest amplification coincides with the greatest level squeezing. For any desired brightness the greatest level of squeezing is obtained for gApκ/2\frac{gA_{p}}{\kappa/2} very close to 1. However, this comes at a cost: the overall uncertainty is increased dramatically in this regime.

When p\mathcal{E}_{p} and s\mathcal{E}_{s} are both positive (Fig 3(b)), the seed signal is deamplified and the output signal is amplitude squeezed. In contrast to the previous case, the high values of gApκ/2\frac{gA_{p}}{\kappa/2} yielding greater squeezing suppress the brightness of the output. This leads to a more complicated picture of squeezing vs. power; depending on the desired power, the greatest squeezing will be associated with different values of gApκ/2\frac{gA_{p}}{\kappa/2}. The situation is also much worse in terms of overall uncertainty.

Refer to caption
(a) Phase squeezing
Refer to caption
(b) Amplitude squeezing
Figure 4: The optimal squeezing attainable for a parametric oscillator operating in the 4(a) amplifying, phase squeezing and 4(b) deamplifying, amplitude squeezing regimes. Each curve shows the best squeezing possible while maintaining total uncertainty below a given threshold.

As in the beam-splitting case, we can summarize these results by plotting the optimal squeezing vs. α2\alpha^{2} in Fig.4 for different thresholds of overall uncertainty. As anticipated, generation of phase-squeezed light can outperform the beam splitter while generation of amplitude-squeezed light underperforms. In this latter case, the optimal strategy would be to generate a squeezed vacuum state and displace it via a beam splitter.

Turning our attention to the case of the parametric amplifier, we must consider non-steady-state dynamics as the pump and seed fields evolve in the medium. This duration becomes the relevant timescale in place of 2/κ2/\kappa in the oscillator case. The equations of motion are presented in Appendix C. To compute the squeezing we integrate equations for the noise operators in Eq.(18) numerically. Fig. 5 shows example trajectories of relevant quantities for both the amplifying, phase squeezing regime and the deamplifying, amplitude squeezing regime. The top panel shows the field amplitudes. It is clear that the amplitudes of both fields change dramatically, even for minimal seeding. Importantly, at long enough times the seed amplitude always vanishes.

Refer to caption
(a) Phase squeezing
Refer to caption
(b) Amplitude squeezing
Figure 5: Dynamics of the squared field amplitudes (top), squeezing (middle), and overall uncertainty (bottom) for 5(a) amplification/phase-squeezing and 5(b) deamplification/amplitude-squeezing.

In the deamplifying, amplitude-squeezing case, the dynamics are straightforward. The seed is depleted while the squeezing increases until the impact of pump noise mixing causes squeezing to degrade along with increasing overall uncertainty. In the amplifying, amplitude-squeezing case the seed amplitude increases as the pump is depleted. Then, at the point where the pump can no longer be depleted, the pump amplitude changes its course to increase and the system enters the deamplifying regime, under which the seed is depleted and the similar dynamics to the deamplifying case prevail. In middle and bottom panels we show the corresponding trajectories of squeezing and overall uncertainty. As with amplitude, the squeezing changes dramatically over time. Initially the squeezing increases linearly. However, the rate of squeezing degrades as the pump depletes and the pump noise and seed noise mix. This reverses when the mixing-in of the pump noise, which includes a substantial component of mixed-in seed noise of opposite sign, exceeds the squeezing. Once the system enters the deamplification regime, the long term effect is phase squeezing and rapidly increasing overall uncertainty as the seed and pump noise continue to mix. This is distinctively different from the amplitude-squeezing only case where such noise mixing does occur only at the excessively squeezing stage. Importantly, the maximum squeezing occurs early in evolution, well before the timing of the local extrema of the amplitudes.

Refer to caption
Refer to caption
Figure 6: The optimal squeezing attainable for a parametric amplifier operating in the 6 amplifying, phase squeezing and 6 deamplifying, amplitude-squeezing regimes constrained to have overall uncertainty below several threshold values. See the details of the dashed s=p\mathcal{E}_{s}=\mathcal{E}_{p} curve in the main text.

These dynamics result in a complicated parameter dependence of squeezing. In particular, in both cases, amplitude maxima and squeezing maxima do not coincide, meaning that optimal operation will depend on choosing the interaction duration appropriately. As before, we plot (Fig. 6) the squeezing as a function of relative squared output displacement α2=As(t)2p2\alpha^{2}=\frac{A_{s}(t)^{2}}{\mathcal{E}_{p}^{2}}. In this case we choose the maximum squeezing subject to uncertainty constraints for a given As(t)2A_{s}(t)^{2} for all tt. The results are broadly consistent with the parametric oscillator performance; however, while phase-squeezing performance is nearly identical, amplitude squeezing performance is superior when overall uncertainty is allowed to be greater. In fact, highly squeezed states are possible at high relative brightness; at higher power levels the performance even exceeds that of beam splitter mixing in the previous section. However, this comes with a caveat: in this regime the required seed input brightness is substantial. At high relative output displacement the seed input must be more than an order of magnitude more powerful than the pump input. The gray dashed line in Fig. 6 shows the cutoff introduced if we restrict the seed input power to that of the pump input power, such that outputs to the right of the line are unavailable with this constraint. So while, in contrast to the parametric oscillator case, the amplifier benefits from significant increases to seed input brightness as a counter to deamplification, in practice the inputs required render the advantage in terms of the relative output power moot.

IV Optomechanical Squeezer

The third method that we consider relies on coupling of optical fields to a mechanical resonator. While squeezing has been demonstrated using a wide variety of implementations [16, 17, 18, 19], here we consider a dissipative optomechanical light squeezing scheme [20]. In this method, the light is squeezed via a strong optomechanical coupling where the mechanical degree of freedom is squeezed due to the interaction with the two-tone optical pumping into the optical cavity. While the authors in [20] considered only the case of vacuum optical seeding (i.e., zero-amplitude seeding into the optical cavity), we expand the analysis to include now the case where the optical cavity is seeded by a non-zero-amplitude coherent light whose frequency is resonant with the optical cavity.

The optomechanical system’s Hamiltonian is given by [21]:

HOM=ωcavaa+Ωbbg(b+b)aa,\displaystyle H_{\mathrm{OM}}=\hbar\omega_{\mathrm{cav}}a^{\dagger}a+\hbar\Omega b^{\dagger}b-\hbar g(b^{\dagger}+b)a^{\dagger}a, (7)

where a,aa,a^{\dagger} are the annihilation and the creation operators for the photonic field, b,bb,b^{\dagger} are those for the phononic field, and ωcav\omega_{\mathrm{cav}}, Ω\Omega are the resonant frequencies of the photonic and phononic cavities, respectively. Here, gg is the single-photon optomechanical coupling coefficient. The cavity is coupled to the external environment (including the input fields) given by the field operators

ain\displaystyle a_{\mathrm{in}} =ain0+eiωcavt(Ep(t)+Es(t)),\displaystyle=a_{\mathrm{in}}^{0}+{\rm e}^{i\omega_{\mathrm{cav}}t}\left(E_{p}(t)+E_{s}(t)\right),
Ep(t)\displaystyle E_{p}(t) =+eiΩt+eiΩt,\displaystyle=\mathcal{E}_{+}e^{i\Omega t}+\mathcal{E}_{-}e^{-i\Omega t},
Es(t)\displaystyle E_{s}(t) =0,\displaystyle=\mathcal{E}_{0}, (8)

where ain0a_{\mathrm{in}}^{0} is a vacuum noise operator and \mathcal{E} is the classical driving field [22]. We set the seed to be XX displaced, so that 0in\mathcal{E}_{0}^{\rm in} is real. The phononic cavity is coupled to thermal noise binb_{\mathrm{in}}.

Assuming that Γ<κΩ\Gamma<\kappa\ll\Omega, generating and solving the equations of motion (Appendix C) we obtain amplitude squeezing for ±\mathcal{E}_{\pm} imaginary given by

α2\displaystyle\alpha^{2} =(1𝒞𝒟)32C2(1+𝒟2),\displaystyle=\frac{\left(1-\mathcal{C}\mathcal{D}\right)^{3}}{2C^{2}\left(1+\mathcal{D}^{2}\right)},
(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =(1+𝒟2)(1𝒞𝒟)2+𝒞𝒟2(2n¯+1),\displaystyle=\frac{\left(1+\mathcal{D}^{2}\right)\left(1-\mathcal{C}\mathcal{D}\right)}{2}+\mathcal{C}\mathcal{D}^{2}\left(2\bar{n}+1\right),
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =(1𝒞𝒟)2(1+𝒞𝒟)2+4(1+1𝒞𝒟)21𝒞𝒟2(2n¯+1),\displaystyle=\frac{\left(1-\mathcal{C}\mathcal{D}\right)^{2}}{\left(1+\mathcal{C}\mathcal{D}\right)^{2}}+\frac{4}{\left(1+\frac{1}{\mathcal{C}\mathcal{D}}\right)^{2}}\frac{1}{\mathcal{C}\mathcal{D}^{2}}\left(2\bar{n}+1\right), (9)

where n¯\bar{n} is the thermal occupation of the mechanical oscillator, 𝒞=4g2|+|2ΓΩ\mathcal{C}=\frac{4g^{2}\left|\mathcal{E}_{+}-\mathcal{E}_{-}\right|^{2}}{\Gamma\Omega} and is determined by the system parameters, 𝒟=12|++||+|\mathcal{D}=\frac{1}{2}\frac{\left|\mathcal{E}_{+}+\mathcal{E}_{-}\right|}{\left|\mathcal{E}_{+}-\mathcal{E}_{-}\right|} and captures the difference in amplitude of the two probe inputs, and we have defined the output squared displacement ratio to be α2=0out||2+|+|2\alpha^{2}=\frac{\mathcal{E}_{0}^{\rm out}}{\left|\mathcal{E}_{-}\right|^{2}+\left|\mathcal{E}_{+}\right|^{2}}. When ±\mathcal{E}_{\pm} are real, the expressions for ΔXs\Delta X_{s} and ΔPs\Delta P_{s} are interchanged, yielding phase squeezing.

The first terms on RHS of the quadrature expressions can be traced to noise contributed by from the probe fields, while the second term originates from the mechanical oscillator. We note that, for a vacuum input leading to a squeezed vacuum output, 𝒞𝒟=1\mathcal{C}\mathcal{D}=1 so that

α2\displaystyle\alpha^{2} =0,\displaystyle=0,
(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =1𝒞(2n¯+1),\displaystyle=\frac{1}{\mathcal{C}}\left(2\bar{n}+1\right),
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =𝒞(2n¯+1),\displaystyle=\mathcal{C}\left(2\bar{n}+1\right),

for amplitude squeezing squeezing (phase squeezing of course yields the same result with the expressions for the quadratures interchanged), consistent with [20]. However, allowing even a low amplitude seed ouptut has a dramatic impact on squeezing. To better understand the relationship between the quadrature variances and the output power, we can write for n¯=0\bar{n}=0 the quadratures to leading order in α2\alpha^{2} as

(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} 𝒟[1+(1𝒟)2𝒟((1+1𝒟2)4α2)13],\displaystyle\approx\mathcal{D}\left[1+\frac{\left(1-\mathcal{D}\right)^{2}}{\mathcal{D}}\left(\frac{\left(1+\frac{1}{\mathcal{D}^{2}}\right)}{4}\alpha^{2}\right)^{\frac{1}{3}}\right],
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} 1𝒟[1(1𝒟)((1+1𝒟2)4α2)23].\displaystyle\approx\frac{1}{\mathcal{D}}\left[1-\left(1-\mathcal{D}\right)\left(\frac{\left(1+\frac{1}{\mathcal{D}^{2}}\right)}{4}\alpha^{2}\right)^{\frac{2}{3}}\right].
Refer to caption
Figure 7: The optimal squeezing attainable for an optomechanical dissipative squeezer assuming n¯=0\bar{n}=0. Each curve shows the best squeezing possible while maintaining total uncertainty below a given threshold.

It is clear that the squeezing is rapidly reduced as the output power increases, as does the deviation of the total uncertainty from unity. When 𝒟\mathcal{D} is small, which is the condition that maximizes squeezing in the vacuum output case, the additional amplitude noise increases rapidly with desired output power, an effect that is exacerbated by the cube root applied to the α2\alpha^{2}. While the phase noise is somewhat reduced, the scaling of this reduction with seed output power is much weaker. This is shown numerically (using Eq.(9)) in Fig. 7, which as before shows the squeezing as a function of relative squared output displacement under different overall uncertainty constraints, here generated for different values of 𝒞\mathcal{C} and 𝒟\mathcal{D}. For brighter outputs modest squeezing is possible regardless of uncertainty desired, but at low finite outputs the performance is quite poor as demonstrated above, with less squeezing actually possible for low desired uncertainty than at higher powers. The brightness vs. squeezing tradeoffs for an optomechanical squeezer are severe enough that the optimal strategy is to produce a squeezed vacuum and displace via a beam splitter.

V General Case

All of the methods explored here exhibit trade-offs between displacement and both squeezing and overall uncertainty. While the specific physics involved dictate the form and severity of these tradeoffs, their existence follows from the use of a pump field to both squeeze and displace the output, and we expect the performance of any parametric-amplification-based method to be similarly afflicted. To see this we note that these Hamiltonians contain cubic terms that yield terms of the form cscpc_{s}c_{p}, where csc_{s} and cpc_{p} may be creation or annihilation operators in the fields, in the equation of motion for the intended squeezed output operator. Partitioning the fields into amplitudes and noise components and linearizing, the contribution to the noise evolution will have the form Cscp+CpcsC_{s}c_{p}+C_{p}c_{s}, where CpC_{p} and CsC_{s} are the C-number amplitudes. The term CpcsC_{p}c_{s} generates the squeezing while CscpC_{s}c_{p} mixes in the noise from the pump field; when the seed amplitude is zero, only squeezing is performed and a pure state results, while nonzero seed amplitude – required to produce displaced output states – necessarily introduces and adds additional noise from the pump field. Therefore any system with this form of Hamiltonian will be subject to the tradeoffs discussed.

VI Conclusion

We presented the performance analysis of generating highly-displaced bright squeezed light from a beam-splitter mixing given a prepared squeezed vacuum as well as generating directly from a seeded optical parametric oscillator, a seeded optical parametric amplifier, and a seeded dissipative optomechanical light squeezer. For each method, we provided a full-quantum-model-based derivation of the expected squeezed (reduced) noise from the main output of the system as a function of mixing degree between the vacuum squeezed mode and a strong coherent-state input. Our results demonstrate that there are limits inherent to the physics of each method that result in tradeoffs between output power (brightness), squeezing, and overall uncertainty that are generic to use of pump field to cause squeezing. All squeezing methods suffer degradation of both squeezing and total uncertainty as the desired brightness is increased. Phase squeezing in either parametric oscillators or amplifiers, as well as amplitude squeezing in parametric amplifiers, are able to outperform the beam splitter mixing of a prepared squeezed vacuum with a bright coherent state, while the optomechanical dissipative squeezer and amplitude-squeezing parametric oscillator perform worse. Additionally, in all cases, it is important to note that parameters that produce the greatest squeezed vacuum do not necessarily provide the best results for bright squeezed states.

Acknowledgements

The article was supported by the Laboratory Directed Research and Development program at Sandia National Laboratories and has been authored by an employee of National Technology & Engineering Solutions of Sandia, LLC under Contract No. DE-NA0003525 with the U.S. Department of Energy (DOE). The employee owns all right, title and interest in and to the article and is solely responsible for its contents. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this article or allow others to do so, for United States Government purposes. The DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan https://www.energy.gov/downloads/doe-public-access-plan This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

Appendix A Beamsplitter

The Heisenberg equations of motion for the beam-splitter Hamiltonian (Eq.(4)) in the rotating frame are simply

as˙=iηas,\displaystyle\dot{a_{s}}=i\eta a_{s},
ap˙=iηap,\displaystyle\dot{a_{p}}=i\eta a_{p},

which can be straightforwardly solved to yield

ap(t)=ap(0)cosηtiap(0)sinηt.\displaystyle a_{p}(t)=a_{p}(0)\cos\eta t-ia_{p}(0)\sin\eta t.

Based on the initial state described above, we have

as(0)\displaystyle a_{s}(0) =Ss(B)asSs(B)=ascoshreiγassinhr,\displaystyle=S_{s}(B)^{\dagger}a_{s}S_{s}(B)=a_{s}\cosh r-e^{i\gamma}a_{s}^{\dagger}\sinh r,
ap(0)\displaystyle a_{p}(0) =Dp()apDp()=ap+,\displaystyle=D_{p}(\mathcal{E})^{\dagger}a_{p}D_{p}(\mathcal{E})=a_{p}+\mathcal{E},

so that we can compute the output variances (for BB and \mathcal{E} real) as

(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =ψo|aXs2|ψos(ψo|sXs|ψos)2\displaystyle=\left\langle\psi_{o}\right|_{a}X_{s}^{2}\left|\psi_{o}\right\rangle_{s}-\left(\left\langle\psi_{o}\right|_{s}X_{s}\left|\psi_{o}\right\rangle_{s}\right)^{2}
=e2Bcos2(ηt)+sin2(ηt),\displaystyle=e^{-2B}\cos^{2}(\eta t)+\sin^{2}(\eta t),
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =ψo|aPs2|ψos(ψo|sPs|ψos)2\displaystyle=\left\langle\psi_{o}\right|_{a}P_{s}^{2}\left|\psi_{o}\right\rangle_{s}-\left(\left\langle\psi_{o}\right|_{s}P_{s}\left|\psi_{o}\right\rangle_{s}\right)^{2}
=e2Bcos2(ηt)+sin2(ηt).\displaystyle=e^{2B}\cos^{2}(\eta t)+\sin^{2}(\eta t).

This gives the squeezing, which is in XX for B>0B>0 and in PP for B<0B<0. The overall uncertainty is then

ΔXsΔPs=1+2cos2(ηt)sin2(ηt)(cosh2B1).\displaystyle\Delta X_{s}\Delta P_{s}=\sqrt{1+2\cos^{2}(\eta t)\sin^{2}(\eta t)\left(\cosh 2B-1\right)}. (10)

The relative squared displacement of the output signal is computed as

α2\displaystyle\alpha^{2} =ψo|Xs|ψo2+ψo|Ps|ψo4||22\displaystyle=\frac{\left\langle\psi_{o}\right|X_{s}\left|\psi_{o}\right\rangle^{2}+\left\langle\psi_{o}\right|P_{s}\left|\psi_{o}\right\rangle}{4\left|\mathcal{E}\right|^{2}}^{2}
=||2sin2(ηt)||2\displaystyle=\frac{\left|\mathcal{E}\right|^{2}\sin^{2}(\eta t)}{\left|\mathcal{E}\right|^{2}}
=sin2(ηt)\displaystyle=\sin^{2}(\eta t) (11)

Appendix B Parametric Oscillator

The rotating frame equations of motion for a parametric oscillator where input and output occur on a single side of the cavity, using the Hamiltonian in equation (6), are given as

as˙\displaystyle\dot{a_{s}} =gasapκ2asκasin,\displaystyle=ga_{s}^{\dagger}a_{p}-\frac{\kappa}{2}a_{s}-\sqrt{\kappa}a_{s}^{\rm in},
ap˙\displaystyle\dot{a_{p}} =(g/2)as2κ2apκapin,\displaystyle=-(g/2)a_{s}^{2}-\frac{\kappa}{2}a_{p}-\sqrt{\kappa}a_{p}^{\rm in},
asout(t)\displaystyle a_{s}^{\rm out}(t) =κas(t)+asin,\displaystyle=\sqrt{\kappa}a_{s}(t)+a_{s}^{\rm in},
apout(t)\displaystyle a_{p}^{\rm out}(t) =κap(t)+apin,\displaystyle=\sqrt{\kappa}a_{p}(t)+a_{p}^{\rm in}, (12)

where asa_{s} and apa_{p} are intracavity operators, asin/outa_{s}^{\rm in/out} and apin/outa_{p}^{\rm in/out} are the input/output field operators, and sin\mathcal{E}_{s}^{\rm in} and pin\mathcal{E}_{p}^{\rm in} are the input field displacements. Here, κ\kappa is the cavity photon decay rate.

Assuming steady state operation, by linearizing the above equations by letting apAp+apa_{p}\rightarrow A_{p}+a_{p}, asAs+asa_{s}\rightarrow A_{s}+a_{s}, where Ap=apA_{p}=\langle a_{p}\rangle and As=asA_{s}=\langle a_{s}\rangle, and dropping terms of order greater than one in the new field operators, we obtain

As\displaystyle A_{s} =2gκAsAp2κsin,\displaystyle=\frac{2g}{\kappa}A_{s}^{*}A_{p}-\frac{2}{\sqrt{\kappa}}\mathcal{E}_{s}^{\rm in},
Ap\displaystyle A_{p} =gκAs22κpin,\displaystyle=-\frac{g}{\kappa}A_{s}^{2}-\frac{2}{\sqrt{\kappa}}\mathcal{E}_{p}^{\rm in},
sout\displaystyle\mathcal{E}_{s}^{\rm out} =κAs+sin,\displaystyle=\sqrt{\kappa}A_{s}+\mathcal{E}_{s}^{\rm in},
as˙\displaystyle\dot{a_{s}} =g(asAp+Asap)κ2asκasin,\displaystyle=g\left(a_{s}^{\dagger}A_{p}+A_{s}^{*}a_{p}\right)-\frac{\kappa}{2}a_{s}-\sqrt{\kappa}a_{s}^{\rm in},
ap˙\displaystyle\dot{a_{p}} =gAsasκ2apκapin,\displaystyle=-gA_{s}a_{s}-\frac{\kappa}{2}a_{p}-\sqrt{\kappa}a_{p}^{\rm in},
asout(t)\displaystyle a_{s}^{\rm out}(t) =κas(t)+asin,\displaystyle=\sqrt{\kappa}a_{s}(t)+a_{s}^{\rm in}, (13)

where pin\mathcal{E}_{p}^{\rm in} and sin\mathcal{E}_{s}^{\rm in} are input field amplitudes and sout\mathcal{E}_{s}^{\rm out} is the squeezed output amplitude. In frequency space the field equations of motion become

as(ω)\displaystyle a_{s}(\omega) =g(as(ω)Ap+Asap(ω))κasin(ω),\displaystyle=g\left(a_{s}(-\omega){{}^{\dagger}}A_{p}+A_{s}^{*}a_{p}(\omega)\right)-\sqrt{\kappa}a_{s}^{\rm in}(\omega),
ap(ω)\displaystyle a_{p}(\omega) =gAsas(ω)κapin(ω).\displaystyle=-gA_{s}a_{s}(\omega)-\sqrt{\kappa}a_{p}^{\rm in}(\omega).

Solving for AsA_{s},ApA_{p} and as(0)a_{s}(0) as in [15], we find

χ\displaystyle\chi =gsin(11+12g2(sin)2(4κgpin+κ23κ)3),\displaystyle=g\mathcal{E}_{s}^{\rm in}\left(1-\sqrt{1+\frac{1}{2g^{2}(\mathcal{E}_{s}^{\rm in})^{2}}\left(\frac{4\sqrt{\kappa}g\mathcal{E}_{p}^{\rm in}+\kappa^{2}}{3\kappa}\right)^{3}}\right),
As\displaystyle A_{s} =κχ1/3g+4gpin+κ3/26gχ1/3,\displaystyle=\frac{-\sqrt{\kappa}\chi^{1/3}}{g}+\frac{4g\mathcal{E}_{p}^{\rm in}+\kappa^{3/2}}{6g\chi^{1/3}},
Ap\displaystyle A_{p} =gAs2κ2κpin,\displaystyle=-\frac{gA_{s}^{2}}{\kappa}-\frac{2}{\sqrt{\kappa}}\mathcal{E}_{p}^{\rm in},
as(0)\displaystyle a_{s}(0) =κ(g2As2+κ24)[κ2asin(0)+gAsapin(0)]+κgAp2[κ2asin(0)+gAsapin(0)](g2|As|2+κ24)2(κgAp2)2\displaystyle=-\sqrt{\kappa}\frac{\left(g^{2}A_{s}^{2}+\frac{\kappa^{2}}{4}\right)\left[\frac{\kappa}{2}a_{s}^{\rm in}(0)+gA_{s}a_{p}^{\rm in}(0)\right]+\frac{\kappa gA_{p}}{2}\left[\frac{\kappa}{2}a_{s}^{\rm in\dagger}(0)+gA_{s}a_{p}^{\rm in\dagger}(0)\right]}{\left(g^{2}\left|A_{s}\right|^{2}+\frac{\kappa^{2}}{4}\right)^{2}-\left(\frac{\kappa gA_{p}}{2}\right)^{2}}

where we have taken ApA_{p} and AsA_{s} to be real. This yields the output quadratures

(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =(g2As2κgAp2κ24)2+(κgAs)2(g2As2κgAp2+κ24)2,\displaystyle=\frac{\left(g^{2}A_{s}^{2}-\frac{\kappa gA_{p}}{2}-\frac{\kappa^{2}}{4}\right)^{2}+\left(\kappa gA_{s}\right)^{2}}{\left(g^{2}A_{s}^{2}-\frac{\kappa gA_{p}}{2}+\frac{\kappa^{2}}{4}\right)^{2}}, (14)
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =(g2As2+κgAp2κ24)2+(κgAs)2(g2As2+κgAp2+κ24)2,\displaystyle=\frac{\left(g^{2}A_{s}^{2}+\frac{\kappa gA_{p}}{2}-\frac{\kappa^{2}}{4}\right)^{2}+\left(\kappa gA_{s}\right)^{2}}{\left(g^{2}A_{s}^{2}+\frac{\kappa gA_{p}}{2}+\frac{\kappa^{2}}{4}\right)^{2}}, (15)

and relative squared output displacement

α2=(soutpin)2=(sin+κAspin)2.\displaystyle\alpha^{2}=\left(\frac{\mathcal{E}_{s}^{\rm out}}{\mathcal{E}_{p}^{\rm in}}\right)^{2}=\left(\frac{\mathcal{E}_{s}^{\rm in}+\sqrt{\kappa}A_{s}}{\mathcal{E}_{p}^{\rm in}}\right)^{2}.

Appendix C Parametric Amplifier

The equations of motion for the parametric amplifier are

as˙\displaystyle\dot{a_{s}} =gasap,\displaystyle=ga_{s}^{\dagger}a_{p},
ap˙\displaystyle\dot{a_{p}} =(g/2)as2,\displaystyle=-(g/2)a_{s}^{2}, (16)

with an initial condition of

|ψ(0)=Ds(s)|0sDp(p)|0p,\left|\psi(0)\right\rangle=D_{s}(\mathcal{E}_{s})\left|0\right\rangle_{s}\otimes D_{p}(\mathcal{E}_{p})\left|0\right\rangle_{p}, (17)

where initially the signal is seeded with a coherent state with an amplitude s\mathcal{E}_{s} and the pump’s initial condition is also a coherent state with an amplitude p\mathcal{E}_{p}. Here, Ds,p()D_{s,p}(\mathcal{E}) is a displacement operator as in equation (3). Once again letting apAp+apa_{p}\rightarrow A_{p}+a_{p}, asAs+asa_{s}\rightarrow A_{s}+a_{s}, where Ap=apA_{p}=\langle a_{p}\rangle and As=asA_{s}=\langle a_{s}\rangle,

As˙\displaystyle\dot{A_{s}} =gAsAp,\displaystyle=gA_{s}^{*}A_{p},
Ap˙\displaystyle\dot{A_{p}} =g2As2,\displaystyle=-\frac{g}{2}A_{s}^{2},
as˙\displaystyle\dot{a_{s}} =g(asAp+Asap+asap),\displaystyle=g\left(a_{s}^{\dagger}A_{p}+A_{s}^{*}a_{p}+a_{s}{{}^{\dagger}}a_{p}\right),
ap˙\displaystyle\dot{a_{p}} =g2(2Asasas2),\displaystyle=-\frac{g}{2}\left(2A_{s}a_{s}-a_{s}^{2}\right),
As(0)\displaystyle A_{s}(0) =s,\displaystyle=\mathcal{E}_{s},
Ap(0)\displaystyle A_{p}(0) =p.\displaystyle=\mathcal{E}_{p}. (18)

It is important to emphasize that the pump that drives the parametric process cannot be assumed to be a constant as before and in previous analyses [23]. To see this we can analytically solve for As(t)A_{s}(t) and Ap(t)A_{p}(t), yielding

As(t)\displaystyle A_{s}(t) =(2c1(1tanh2(gtc1+c1c2))1/2,\displaystyle=\left(2c_{1}(1-\tanh^{2}(gt\sqrt{c_{1}}+\sqrt{c_{1}}c_{2})\right)^{1/2},
Ap(t)\displaystyle A_{p}(t) =c1tanh(gtc1+c1c2),\displaystyle=-\sqrt{c_{1}}\tanh(gt\sqrt{c_{1}}+\sqrt{c_{1}}c_{2}), (19)

where

c1\displaystyle c_{1} =p2+s2/2,\displaystyle=\mathcal{E}_{p}^{2}+\mathcal{E}_{s}^{2}/2,
c2\displaystyle c_{2} =tanh1(p/p2+s2/2)/p2+s2/2.\displaystyle=-\tanh^{-1}\left(\mathcal{E}_{p}/\sqrt{\mathcal{E}_{p}^{2}+\mathcal{E}_{s}^{2}/2}\right)/\sqrt{\mathcal{E}_{p}^{2}+\mathcal{E}_{s}^{2}/2}. (20)

Appendix D Optomechanical Squeezer

Starting from the optomechanical system’s Hamiltonian and input field descriptions (Eqs. (7) and (8)), the equations of motion in the rotating frame are

a˙\displaystyle\dot{a} =ig(beiΩt+beiΩt)aκ2aκ(ain+p(t)+s(t)),\displaystyle=ig\left(b{{}^{\dagger}}{\rm e}^{i\Omega t}+b{\rm e}^{-i\Omega t}\right)a-\frac{\kappa}{2}a-\sqrt{\kappa}\left(a^{in}+\mathcal{E}_{p}(t)+\mathcal{E}_{s}(t)\right),
b˙\displaystyle\dot{b} =igaaeiΩtΓ2Γbin.\displaystyle=iga{{}^{\dagger}}a{\rm e}^{i\Omega t}-\frac{\Gamma}{2}-\sqrt{\Gamma}b^{in}.

which in Fourier space become

a(ω)\displaystyle a(\omega) =igdω[a(ω)b(ω+ω+Ω)+a(ω)b(ωω+Ω)]κ(ain(ω)+p(ω)+s(ω))κ2iω,\displaystyle=\frac{ig\int d\omega^{\prime}\left[a(\omega^{\prime})b(-\omega+\omega^{\prime}+\Omega){{}^{\dagger}}+a(\omega^{\prime})b(\omega-\omega^{\prime}+\Omega)\right]-\sqrt{\kappa}\left(a^{in}(\omega)+\mathcal{E}_{p}(\omega)+\mathcal{E}_{s}(\omega)\right)}{\frac{\kappa}{2}-i\omega},
b(ω)\displaystyle b(\omega) =ig𝑑ωa(ω)a(ωωΩ)Γbin(ω)Γ2iω,\displaystyle=\frac{ig\int d\omega^{\prime}a(-\omega^{\prime}){{}^{\dagger}}a(\omega-\omega^{\prime}-\Omega)-\sqrt{\Gamma}b^{in}(\omega)}{\frac{\Gamma}{2}-i\omega},

we take a(ω)A(ω)+a(ω)a(\omega)\rightarrow A(\omega)+a(\omega) and b(ω)b(ω)+b(ω)b(\omega)\rightarrow b(\omega)+b(\omega) and discard the terms higher than linear order in the field operators to obtain

A(ω)\displaystyle A(\omega) =igdωA(ω)[B(ω+ω+Ω)+B(ωω+Ω)]κ(p(ω)+s(ω))κ2iω,\displaystyle=\frac{ig\int d\omega^{\prime}A(\omega^{\prime})\left[B(-\omega+\omega^{\prime}+\Omega){{}^{\dagger}}+B(\omega-\omega^{\prime}+\Omega)\right]-\sqrt{\kappa}\left(\mathcal{E}_{p}(\omega)+\mathcal{E}_{s}(\omega)\right)}{\frac{\kappa}{2}-i\omega},
B(ω)\displaystyle B(\omega) =ig𝑑ωA(ω)A(ωωΩ)Γ2iω,\displaystyle=\frac{ig\int d\omega^{\prime}A(-\omega^{\prime}){{}^{\dagger}}A(\omega-\omega^{\prime}-\Omega)}{\frac{\Gamma}{2}-i\omega},
a(ω)\displaystyle a(\omega) =1κ2iω{igdω[A(ω)(b(ω+ω+Ω)+b(ωω+Ω))\displaystyle=\frac{1}{\frac{\kappa}{2}-i\omega}\left\{ig\int d\omega^{\prime}\left[A(\omega^{\prime})\left(b(-\omega+\omega^{\prime}+\Omega){{}^{\dagger}}+b(\omega-\omega^{\prime}+\Omega)\right)\right.\right.
+a(ω)(B(ω+ω+Ω)+B(ωω+Ω))]κain(ω)},\displaystyle\left.\left.+a(\omega^{\prime})\left(B(-\omega+\omega^{\prime}+\Omega){{}^{\dagger}}+B(\omega-\omega^{\prime}+\Omega)\right)\right]-\sqrt{\kappa}a^{in}(\omega)\right\},
b(ω)\displaystyle b(\omega) =ig𝑑ω[A(ω)a(ωωΩ)+a(ω)A(ωωΩ)]Γbin(ω)Γ2iω.\displaystyle=\frac{ig\int d\omega^{\prime}\left[A(-\omega^{\prime}){{}^{\dagger}}a(\omega-\omega^{\prime}-\Omega)+a(-\omega^{\prime}){{}^{\dagger}}A(\omega-\omega^{\prime}-\Omega)\right]-\sqrt{\Gamma}b^{in}(\omega)}{\frac{\Gamma}{2}-i\omega}.

The driving fields will ensure that at steady state the average amplitudes will be delta functions at multiples of Ω\Omega; i.e.,we can write A(ω)=nAnδ(ωnΩ)A(\omega)=\sum_{n}A_{n}\delta(\omega-n\Omega) and B(ω)=nBnδ(ωnΩ)B(\omega)=\sum_{n}B_{n}\delta(\omega-n\Omega). We then have

A0\displaystyle A_{0} =igjAj[B1+j+B1j]κ0κ2,\displaystyle=\frac{ig\sum_{j}A_{j}\left[B_{1+j}{{}^{\dagger}}+B_{1-j}\right]-\sqrt{\kappa}\mathcal{E}_{0}}{\frac{\kappa}{2}},
A+1\displaystyle A_{+1} =igjAj[Bj+B2j]κ+κ2iΩ,\displaystyle=\frac{ig\sum_{j}A_{j}\left[B_{j}{{}^{\dagger}}+B_{2-j}\right]-\sqrt{\kappa}\mathcal{E}_{+}}{\frac{\kappa}{2}-i\Omega},
A1\displaystyle A_{-1} =igjAj[B2+j+Bj]κκ2+iΩ,\displaystyle=\frac{ig\sum_{j}A_{j}\left[B_{2+j}{{}^{\dagger}}+B_{-j}\right]-\sqrt{\kappa}\mathcal{E}_{-}}{\frac{\kappa}{2}+i\Omega},
B0\displaystyle B_{0} =igAjAj1Γ2,\displaystyle=\frac{igA_{-j}{{}^{\dagger}}A_{-j-1}}{\frac{\Gamma}{2}},
a(ω)\displaystyle a(\omega) =1κ2iω{ig[Aj(b(ω+(1+j)Ω)+b(ω+(1j)Ω))+\displaystyle=\frac{1}{\frac{\kappa}{2}-i\omega}\left\{ig\left[A_{j}\left(b(-\omega+(1+j)\Omega){{}^{\dagger}}+b(\omega+(1-j)\Omega)\right)+\right.\right.
a(ω+(j1)Ω)Bj+a(ω+(1j)Ω)Bj]κain(ω)},\displaystyle\left.\left.a(\omega+(j-1)\Omega)B_{j}{{}^{\dagger}}+a(\omega+(1-j)\Omega)B_{j}\right]-\sqrt{\kappa}a^{in}(\omega)\right\},
b(ω)\displaystyle b(\omega) =ig𝑑ω[Aja(ω(j+1)Ω)+a(ωΩ)Aj]Γbin(ω)Γ2iω.\displaystyle=\frac{ig\int d\omega^{\prime}\left[A_{-j}{{}^{\dagger}}a(\omega-(j+1)\Omega)+a(\omega-\Omega){{}^{\dagger}}A_{j}\right]-\sqrt{\Gamma}b^{in}(\omega)}{\frac{\Gamma}{2}-i\omega}.

We will assume that κΩ\kappa\ll\Omega so that we can restrict our consideration to A0A_{0} and A±ΩA_{\pm\Omega}, the latter of which we will write as A±A_{\pm} for brevity. As a result, since ultimately we desire a(0)a(0), upon linearizing we need only concern ourselves with the noise operator at ω=0\omega=0 and ω=±Ω\omega=\pm\Omega and write them as a0a_{0} and a±a_{\pm}, respectively. We will also assume ΓΩ\Gamma\ll\Omega so that we need only consider b0=b(0)b_{0}=b(0). We then have

A0\displaystyle A_{0} =1κ112+g2κΓ(|A+|2|A|2)0,\displaystyle=\frac{1}{\sqrt{\kappa}}\frac{1}{\frac{1}{2}+\frac{g^{2}}{\kappa\Gamma}\left(\left|A_{+}\right|^{2}-\left|A_{-}\right|^{2}\right)}\mathcal{E}_{0},
A±\displaystyle A_{\pm} =κ(κ2iΩ)2[2g2A02Γ(κ2iΩ±2g2|A0|2Γ)±],\displaystyle=\frac{\sqrt{\kappa}}{\left(\frac{\kappa}{2}\mp i\Omega\right)^{2}}\left[\mp\frac{2g^{2}A_{0}^{2}}{\Gamma}\mathcal{E}_{\mp}^{*}-\left(\frac{\kappa}{2}\mp i\Omega\pm\frac{2g^{2}\left|A_{0}\right|^{2}}{\Gamma}\right)\mathcal{E}_{\pm}\right],
a0\displaystyle a_{0} =g2A0[AaA+a++AaA+a+]igΓ[Ab0in,+A+b0in]κΓ2a0inκΓ4+g2(|A+|2|A|2),\displaystyle=\frac{g^{2}A_{0}\left[A_{-}{{}^{\dagger}}a_{-}-A_{+}{{}^{\dagger}}a_{+}+A_{-}a_{-}{{}^{\dagger}}-A_{+}a_{+}{{}^{\dagger}}\right]-ig\sqrt{\Gamma}\left[A_{-}b_{0}^{in,\dagger}+A_{+}b_{0}^{in}\right]-\frac{\sqrt{\kappa}\Gamma}{2}a^{in}_{0}}{\frac{\kappa\Gamma}{4}+g^{2}\left(\left|A_{+}\right|^{2}-\left|A_{-}\right|^{2}\right)},
a±\displaystyle a_{\pm} ±g2[(2AA0+A0A±)a0+A0A±a0+A0A0a]igΓA0b0in,±κΓ2a±inκΓ4iΓΩ2g2|A0|2.\displaystyle\approx\frac{\pm g^{2}\left[\left(2A_{\mp}{{}^{\dagger}}A_{0}+A_{0}{{}^{\dagger}}A_{\pm}\right)a_{0}+A_{0}A_{\pm}a_{0}{{}^{\dagger}}+A_{0}A_{0}a_{\mp}{{}^{\dagger}}\right]-ig\sqrt{\Gamma}A_{0}b_{0}^{in,\pm\dagger}-\frac{\sqrt{\kappa}\Gamma}{2}a^{in}_{\pm}}{\frac{\kappa\Gamma}{4}\mp i\frac{\Gamma\Omega}{2}\mp g^{2}\left|A_{0}\right|^{2}}.

Then, assuming A0A_{0} is real, solving these equations we find that the output quadrature

(ΔXs)2\displaystyle\left(\Delta X_{s}\right)^{2} =|κx0+1|2+κ|x+|2+κ|x|2+κ|xb|2,\displaystyle=\left|\sqrt{\kappa}x_{0}+1\right|^{2}+\kappa\left|x_{+}\right|^{2}+\kappa\left|x_{-}\right|^{2}+\kappa\left|x_{b}\right|^{2},
x0\displaystyle x_{0} =[i2Cg4A02(Cg2A02CΩ2[A+A])CΩ6|A+A+|4+12(1+Δ)(Δ+2Cg4A02CΩ4|A+A+|2)](1+Δ)(12(1+Δ)+Cg4A02CΩ4|A+A+|2),\displaystyle=\frac{\left[i\frac{2C_{g}^{4}A_{0}^{2}\left(C_{g}^{2}A_{0}^{2}-C_{\Omega}^{2}\Im\left[A_{+}A_{-}\right]\right)}{C_{\Omega}^{6}\left|A_{-}^{*}+A_{+}\right|^{4}}+\frac{1}{2}\left(1+\Delta\right)\left(\Delta+2\frac{C_{g}^{4}A_{0}^{2}}{C_{\Omega}^{4}\left|A_{-}^{*}+A_{+}\right|^{2}}\right)\right]}{\left(1+\Delta\right)\left(\frac{1}{2}\left(1+\Delta\right)+\frac{C_{g}^{4}A_{0}^{2}}{C_{\Omega}^{4}\left|A_{-}^{*}+A_{+}\right|^{2}}\right)},
x±\displaystyle x_{\pm} =1κCg2A0|A+A+|(Cg2A02|A+A+|2A+A±|A+A+|iA±|A+A+|CΩ2)CΩ4(12(1+Δ)+Cg4A02CΩ4|A+A+|2),\displaystyle=\frac{1}{\sqrt{\kappa}}\frac{-C_{g}^{2}\frac{A_{0}}{\left|A_{-}^{*}+A_{+}\right|}\left(C_{g}^{2}\frac{A_{0}^{2}}{\left|A_{-}^{*}+A_{+}\right|^{2}}\frac{A_{\mp}+A_{\pm}^{*}}{\left|A_{-}^{*}+A_{+}\right|}-i\frac{A_{\pm}^{*}}{\left|A_{-}^{*}+A_{+}\right|}C_{\Omega}^{2}\right)}{C_{\Omega}^{4}\left(\frac{1}{2}\left(1+\Delta\right)+\frac{C_{g}^{4}A_{0}^{2}}{C_{\Omega}^{4}\left|A_{-}^{*}+A_{+}\right|^{2}}\right)},
xb\displaystyle x_{b} =1Γi2Cg[i2Cg2CΩ2A02|A+A+|2A+A+|A+A+|(Cg2A02(Cg2A02CΩ2[A+A])CΩ4|A+A+|4+1+Δ2)+1+Δ2AA+|A+A+|](1+Δ)(12(1+Δ)+Cg4A02CΩ4|A+A+|2),\displaystyle=\frac{1}{\sqrt{\Gamma}}\frac{i2C_{g}\left[i2\frac{C_{g}^{2}}{C_{\Omega}^{2}}\frac{A_{0}^{2}}{\left|A_{-}^{*}+A_{+}\right|^{2}}\frac{A_{-}^{*}+A_{+}}{\left|A_{-}^{*}+A_{+}\right|}\left(\frac{C_{g}^{2}A_{0}^{2}\left(C_{g}^{2}A_{0}^{2}-C_{\Omega}^{2}\Im\left[A_{+}A_{-}\right]\right)}{C_{\Omega}^{4}\left|A_{-}^{*}+A_{+}\right|^{4}}+\frac{1+\Delta}{2}\right)+\frac{1+\Delta}{2}\frac{A_{-}^{*}-A_{+}}{\left|A_{-}^{*}+A_{+}\right|}\right]}{\left(1+\Delta\right)\left(\frac{1}{2}\left(1+\Delta\right)+\frac{C_{g}^{4}A_{0}^{2}}{C_{\Omega}^{4}\left|A_{-}^{*}+A_{+}\right|^{2}}\right)},
(ΔPs)2\displaystyle\left(\Delta P_{s}\right)^{2} =|κp0+1|2+κ|p+|2+κ|p|2+κ|pb|2,\displaystyle=\left|\sqrt{\kappa}p_{0}+1\right|^{2}+\kappa\left|p_{+}\right|^{2}+\kappa\left|p_{-}\right|^{2}+\kappa\left|p_{b}\right|^{2},
p0\displaystyle p_{0} =1κ(1+Δ),\displaystyle=-\frac{1}{\sqrt{\kappa}\left(1+\Delta\right)},
p±\displaystyle p_{\pm} =0,\displaystyle=0,
pb\displaystyle p_{b} =2CgΓ(1+Δ),\displaystyle=\frac{2C_{g}}{\sqrt{\Gamma}\left(1+\Delta\right)},
α2\displaystyle\alpha^{2} =Δ2(1+Δ)02|+|2+||2,\displaystyle=\frac{\Delta^{2}}{\left(1+\Delta\right)}\frac{\mathcal{E}_{0}^{2}}{\left|\mathcal{E}_{+}\right|^{2}+\left|\mathcal{E}_{-}\right|^{2}},

where we have defined the quantities

Cg\displaystyle C_{g} =g|A++A|κ,\displaystyle=\frac{g\left|A_{+}+A_{-}^{*}\right|}{\kappa},
CΩ2\displaystyle C_{\Omega}^{2} =ΓΩ2κ2,\displaystyle=\frac{\Gamma\Omega}{2\kappa^{2}},

for convenience. The variance can be nearly – though not quite exactly – minimized by setting Δ=2Cg4A02CΩ4|A+A+|2\Delta=-2\frac{C_{g}^{4}A_{0}^{2}}{C_{\Omega}^{4}\left|A_{-}^{*}+A_{+}\right|^{2}} and taking ±\mathcal{E}_{\pm} to be purely imaginary and of opposite sign, resulting in Eqs.(9)

References

  • Schnabel [2017] R. Schnabel, Squeezed states of light and their applications in laser interferometers, Physics Reports 684, 1 (2017), squeezed states of light and their applications in laser interferometers.
  • Zhuang et al. [2018] Q. Zhuang, Z. Zhang, and J. H. Shapiro, Distributed quantum sensing using continuous-variable multipartite entanglement, Physical Review A 97, 032329 (2018).
  • Guo et al. [2020] X. Guo, C. R. Breum, J. Borregaard, S. Izumi, M. V. Larsen, T. Gehring, M. Christandl, J. S. Neergaard-Nielsen, and U. L. Andersen, Distributed quantum sensing in a continuous-variable entangled network, Nature Physics 16, 281 (2020).
  • Atkinson et al. [2021] G. Atkinson, E. Allen, G. Ferranti, A. McMillan, and J. Matthews, Quantum enhanced precision estimation of transmission with bright squeezed light, Physical Review Applied 16, 044031 (2021).
  • Soh and Chatterjee [2023] D. Soh and E. Chatterjee, Label-free quantum super-resolution imaging using entangled multi-mode squeezed light, New Journal of Physics , 093001 (2023).
  • Kuzmich et al. [1997] A. Kuzmich, K. Mølmer, and E. Polzik, Spin squeezing in an ensemble of atoms illuminated with squeezed light, Physical review letters 79, 4782 (1997).
  • Milburn [1984] G. Milburn, Interaction of a two-level atom with squeezed light, Optica Acta: International Journal of Optics 31, 671 (1984).
  • Polzik et al. [1992] E. Polzik, J. Carri, and H. Kimble, Atomic spectroscopy with squeezed light for sensitivity beyond the vacuum-state limit, Applied Physics B 55, 279 (1992).
  • Caves [1981] C. M. Caves, Quantum-mechanical noise in an interferometer, Physical Review D 23, 1693 (1981).
  • Vahlbruch et al. [2008] H. Vahlbruch, M. Mehmet, S. Chelkowski, B. Hage, A. Franzen, N. Lastzka, S. Goßler, K. Danzmann, and R. Schnabel, Observation of squeezed light with 10-db quantum-noise reduction, Phys. Rev. Lett. 100, 033602 (2008).
  • Eberle et al. [2010] T. Eberle, S. Steinlechner, J. Bauchrowitz, V. Händchen, H. Vahlbruch, M. Mehmet, H. Müller-Ebhardt, and R. Schnabel, Quantum enhancement of the zero-area sagnac interferometer topology for gravitational wave detection, Phys. Rev. Lett. 104, 251102 (2010).
  • Stefszky et al. [2012] M. S. Stefszky, C. M. Mow-Lowry, S. S. Y. Chua, D. A. Shaddock, B. C. Buchler, H. Vahlbruch, A. Khalaidovski, R. Schnabel, P. K. Lam, and D. E. McClelland, Balanced homodyne detection of optical quantum states at audio-band frequencies and below, Classical and Quantum Gravity 29, 145015 (2012).
  • Vahlbruch et al. [2016] H. Vahlbruch, M. Mehmet, K. Danzmann, and R. Schnabel, Detection of 15 db squeezed states of light and their application for the absolute calibration of photoelectric quantum efficiency, Phys. Rev. Lett. 117, 110801 (2016).
  • Collett and Gardiner [1984] M. J. Collett and C. W. Gardiner, Squeezing of intracavity and traveling-wave light fields produced in parametric amplification, Phys. Rev. A 30, 1386 (1984).
  • Lariontsev and Zolotoverkh [2002] E. Lariontsev and I. Zolotoverkh, Characteristics of bright-squeezed light produced in a below-threshold optical parametric oscillator, Journal of Optics B: Quantum and Semiclassical Optics 4, 15 (2002).
  • Brooks et al. [2012] D. W. C. Brooks, T. Botter, S. Schreppler, T. P. Purdy, N. Brahms, and D. M. Stamper-Kurn, Non-classical light generated by quantum-noise-driven cavity optomechanics, Nature 488, 476 (2012).
  • Safavi-Naeini et al. [2013] A. H. Safavi-Naeini, S. Gröblacher, J. T. Hill, J. Chan, M. Aspelmeyer, and O. Painter, Squeezed light from a silicon micromechanical resonator, Nature 500, 185 (2013).
  • Purdy et al. [2013] T. P. Purdy, P.-L. Yu, R. W. Peterson, N. S. Kampel, and C. A. Regal, Strong optomechanical squeezing of light, Phys. Rev. X 3, 031012 (2013).
  • Xiong et al. [2022] B. Xiong, S. Chao, C. Shan, and J. Liu, Optomechanical squeezing with pulse modulation, Opt. Lett. 47, 5545 (2022).
  • Kronwald et al. [2014] A. Kronwald, F. Marquardt, and A. A. Clerk, Dissipative optomechanical squeezing of light, New Journal of Physics 16, 063058 (2014).
  • Aspelmeyer et al. [2014] M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, Cavity optomechanics, Reviews of Modern Physics 86, 1391 (2014).
  • Gardiner and Collett [1985] C. W. Gardiner and M. J. Collett, Input and output in damped quantum systems: Quantum stochastic differential equations and the master equation, Physical Review A 31, 3761 (1985).
  • Milburn and Walls [1981] G. Milburn and D. Walls, Production of squeezed states in a degenerate parametric amplifier, Optics Communications 39, 401 (1981).