pith. sign in

arxiv: 1907.09572 · v2 · pith:22E6X2QOnew · submitted 2019-07-22 · 🪐 quant-ph

Steady states, squeezing, and entanglement in intracavity triplet down conversion

Mathew D. E. Denys , Murray K. Olsen , Luke S. Trainor , Harald G. L. Schwefel , Ashton S. Bradley This is my paper

Pith reviewed 2026-05-24 17:46 UTC · model grok-4.3

classification 🪐 quant-ph
keywords triplet down conversionintracavitysqueezingentanglementpositive-P representationsteady statesquantum opticsnonlinear optics
0
0 comments X

The pith

Quantum effects like squeezing and entanglement peak just above the semi-classical threshold in intracavity triplet down conversion.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines degenerate intracavity triplet down conversion, in which one high-energy photon splits into three lower-energy photons inside an optical resonator. Stochastic differential equations derived from the truncated positive-P representation are solved to track the time evolution of the intracavity mode populations and to extract their steady-state values as a function of pump strength. The resulting output-field spectra show that quantum corrections become largest immediately above the classical pumping threshold, where they deviate from semi-classical expectations. From these spectra the authors locate parameter regions in which measurable squeezing and bipartite entanglement should appear. The method is checked against Monte Carlo wave-function simulations and matches well when photon numbers remain low.

Core claim

In the steady state of intracavity triplet down conversion, quantum effects are strongest immediately above the semi-classical pumping threshold; the output spectra there exhibit squeezing and bipartite entanglement that semi-classical theory does not predict.

What carries the argument

The truncated positive-P representation, which converts the master equation into stochastic differential equations whose ensemble averages yield the intracavity populations and output spectra.

If this is right

  • Steady-state populations of the pump and signal modes can be obtained directly from the stochastic equations as pump intensity is varied.
  • Regimes exist in which the output fields display measurable quadrature squeezing.
  • Bipartite entanglement between the signal modes is present in identifiable operating regions.
  • Numerical results differ significantly from semi-classical predictions near threshold.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Experimental searches for triplet down conversion could target pump powers just above threshold to maximize observable quantum correlations.
  • The same stochastic approach may be applied to non-degenerate or multimode versions of the process without major reformulation.
  • If the threshold behavior generalizes, similar deviations from classical predictions could appear in other three-photon nonlinear optical systems.

Load-bearing premise

The truncated positive-P representation remains accurate for the quantum dynamics even when its validation is restricted to low mode populations.

What would settle it

Observation of cavity output spectra that match semi-classical predictions rather than the reported quantum deviations immediately above threshold would falsify the central claim.

Figures

Figures reproduced from arXiv: 1907.09572 by Ashton S. Bradley, Harald G. L. Schwefel, Luke S. Trainor, Mathew D. E. Denys, Murray K. Olsen.

Figure 1
Figure 1. Figure 1: FIG. 1: Our model involves an optical cavity with two occupied res [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: The steady-states of [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: The time evolution of the expected mode populations for a [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The steady-states of [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The steady-state Duan–Simon spectra of the output fields for [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The spectrum of quadrature variances for the output fields [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: For each mode, the steady-state standard deviations of the [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: A comparison between the time evolution of the expected [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
read the original abstract

Triplet down conversion, the process of converting one high-energy photon into three low-energy photons, may soon be experimentally feasible due to advances in optical resonator technology. We use quantum phase-space techniques to analyse the process of degenerate intracavity triplet down conversion by solving stochastic differential equations within the truncated positive-P representation. The time evolution of both intracavity mode populations are simulated, and the resulting steady-states are examined as a function of the pump intensity. Quantum effects are most pronounced in the region immediately above the semi-classical pumping threshold, where our numerical results differ significantly from semi-classical predictions. Regimes of measurable squeezing and bipartite entanglement are identified from steady-state spectra of the cavity output fields. We validate the truncated positive-P description against Monte Carlo wave function simulations, finding good agreement for low mode populations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper analyzes degenerate intracavity triplet down conversion using the truncated positive-P representation to solve stochastic differential equations. It reports time evolution and steady-state populations as a function of pump intensity, identifies a regime immediately above the semi-classical threshold where quantum effects (squeezing and bipartite entanglement) are pronounced and differ from semi-classical predictions, and extracts measurable signatures from steady-state output spectra. Validation against Monte Carlo wave-function simulations is reported only for low mode populations.

Significance. If the numerical results hold, the work identifies experimentally accessible parameter regimes for observable quantum correlations in a process that may become feasible with current resonator technology, extending phase-space methods to a three-photon nonlinear process.

major comments (2)
  1. [Abstract] Abstract and validation section: good agreement with MCWF is stated only for low mode populations, yet the central claim concerns the above-threshold regime where intracavity populations increase; without truncation-error bounds or convergence checks in that regime, the reported spectra and entanglement measures rest on an unverified extrapolation of the approximation.
  2. [Results] Results on squeezing and entanglement: the identification of measurable regimes from steady-state spectra assumes the truncated positive-P SDEs remain accurate above threshold; the manuscript provides no parameter details, error bars, or population values at the reported operating points to allow assessment of truncation validity.
minor comments (2)
  1. [Abstract] The abstract and main text should state the truncation order and any cutoff parameters used in the positive-P representation.
  2. [Figures] Figure captions and text should report the specific pump intensities and steady-state populations at which squeezing and entanglement are claimed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive report. The comments highlight important points regarding the validation of the truncated positive-P method in the above-threshold regime. We address each major comment below and will revise the manuscript to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Abstract] Abstract and validation section: good agreement with MCWF is stated only for low mode populations, yet the central claim concerns the above-threshold regime where intracavity populations increase; without truncation-error bounds or convergence checks in that regime, the reported spectra and entanglement measures rest on an unverified extrapolation of the approximation.

    Authors: We agree that the validation is explicitly limited to low populations and that additional checks would strengthen the claims for the above-threshold regime. In the revised manuscript we will report the specific intracavity population values at the operating points used for the spectra and entanglement calculations. We will also add a discussion of the expected range of validity of the truncation, drawing on the known properties of the positive-P representation for similar nonlinear processes, and include any feasible additional MCWF comparisons at moderate populations. revision: yes

  2. Referee: [Results] Results on squeezing and entanglement: the identification of measurable regimes from steady-state spectra assumes the truncated positive-P SDEs remain accurate above threshold; the manuscript provides no parameter details, error bars, or population values at the reported operating points to allow assessment of truncation validity.

    Authors: The referee is correct that the manuscript would benefit from explicit parameter values and population numbers to allow readers to assess truncation validity. We will revise the results section to include the pump intensities, steady-state populations, and any available statistical error estimates from the stochastic trajectories at the reported operating points. This will enable direct evaluation of the regime in which the reported squeezing and entanglement signatures are obtained. revision: yes

Circularity Check

0 steps flagged

No circularity: direct numerical integration validated against independent Monte Carlo method

full rationale

The paper solves stochastic differential equations in the truncated positive-P representation and compares steady-state spectra to Monte Carlo wave function simulations, reporting agreement only for low populations. No load-bearing step reduces a prediction to a fitted input, self-definition, or self-citation chain; the central results are outputs of the numerical integration rather than redefinitions of its inputs. The limited validation regime is a scope limitation, not a circularity in the derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the domain assumption that the truncated positive-P representation remains accurate for the studied pump intensities and that the Monte Carlo comparison validates it at low populations. No free parameters or invented entities are mentioned.

axioms (1)
  • domain assumption Truncated positive-P representation is a valid approximation for the intracavity triplet down conversion dynamics
    Invoked to justify solving the stochastic differential equations and extracting steady-state spectra.

pith-pipeline@v0.9.0 · 5680 in / 1165 out tokens · 68700 ms · 2026-05-24T17:46:31.170951+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

  1. [1]

    to simulate the truncated positive-P SDE (6). Expectation values of the mode populations⟨ˆna⟩ =⟨ˆa† ˆa⟩ and⟨ˆnb⟩ =⟨ˆb† ˆb⟩ were calculated by averaging over 106 trajectories, which was found to be sufficient for convergence of the second-order mo- ments used to characterize the system throughout this work. For comparison with a semi-classical treatment, mod...

    work page

  2. [2]

    2: The time evolution of the expected mode populations for a system with the parameters κ = 0.001γa,γ b = 2γa,ϵ a = 5γa, and ϵb = 200γa

    5 3 ·104 (a) γat population |α s|2 (analytic) |β s|2 (analytic) ⟨na⟩ (ODE) ⟨nb⟩ (ODE) ⟨na⟩ (SDE) ⟨nb⟩ (SDE) 12 14 16 18 20 22 24 26 28 30 6, 400 6, 500 6, 600 6, 700 (b) γat ⟨na⟩ ⟨na⟩ (SDE, ϵ a turned off) 12 14 16 18 20 22 24 26 28 30120 140 160 180 (c) γat ⟨nb⟩ ⟨nb⟩ (SDE, ϵ a turned off) FIG. 2: The time evolution of the expected mode populations for a sy...

    work page

  3. [3]

    6: The steady-state Duan–Simon spectra of the output fields for the cases that an injected signal is (orange) and is not (blue) present in the steady-state

    5 5 ω/γ a DS± DS+ (ϵa = 5γa) DS− (ϵa = 5γa) DS+ (ϵa = 0) DS− (ϵa = 0) FIG. 6: The steady-state Duan–Simon spectra of the output fields for the cases that an injected signal is (orange) and is not (blue) present in the steady-state. The Duan–Simon inequality (21) is violated by DS−, indicating that the two fields are entangled. The parameters used are κ = 0....

    work page

  4. [4]

    8 1 C[Ya, Y a] 0 50 100 150 200 250 3000 5 10 15 20 pumping, ϵb/γ a frequency, ω/γ a 0

    work page

  5. [5]

    7: The spectrum of quadrature variances for the output fields as a function of the pump intensity ϵb

    8 1 C[Yb, Y b] FIG. 7: The spectrum of quadrature variances for the output fields as a function of the pump intensity ϵb. The parameters used are κ = 0.001γa,γ b = 2γa and ϵa = 0, and simulations were averaged over 106 trajectories. The region enclosed by the dashed lines indicates where the linearized fluctuation analysis does not apply due to the magnitud...

    work page

  6. [6]

    M. K. Olsen, Phys. Rev. A 96, 063839 (2017)

    work page 2017

  7. [7]

    D. B. Kolker, N. Y . Kostyukova, A. A. Boyko, V . V . Badikov, D. V . Badikov, A. G. Shadrintseva, N. N. Tretyakova, K. G. Zenov, A. A. Karapuzikov, and J.-J. Zondy, J. Phys. Comm. 2, 035039 (2018)

    work page 2018

  8. [8]

    M. K. Olsen, Opt. Comm. 410, 966 (2018). 9

    work page 2018

  9. [9]

    Li and M

    J. Li and M. K. Olsen, Phys. Rev. A 97, 043856 (2018)

    work page 2018

  10. [10]

    J. U. Frst, D. V . Strekalov, D. Elser, A. Aiello, U. L. Andersen, C. Marquardt, and G. Leuchs, Phys. Rev. Lett. 106, 113901 (2011)

    work page 2011

  11. [11]

    Peano, H

    V . Peano, H. G. L. Schwefel, C. Marquardt, and F. Marquardt, Phys. Rev. Lett. 115, 243603 (2015)

    work page 2015

  12. [12]

    U. L. Andersen, T. Gehring, C. Marquardt, and G. Leuchs, Phys. Scr. 91, 053001 (2016)

    work page 2016

  13. [13]

    Vahlbruch, M

    H. Vahlbruch, M. Mehmet, K. Danzmann, and R. Schnabel, Phys. Rev. Lett. 117, 110801 (2016)

    work page 2016

  14. [14]

    Takei, H

    N. Takei, H. Yonezawa, T. Aoki, and A. Furusawa, Phys. Rev. Lett. 94, 220502 (2005)

    work page 2005

  15. [15]

    J. Jing, J. Zhang, Y . Yan, F. Zhao, C. Xie, and K. Peng, Phys. Rev. Lett. 90, 167903 (2003)

    work page 2003

  16. [16]

    Cavanna, F

    A. Cavanna, F. Just, X. Jiang, G. Leuchs, M. V . Chekhova, P. S. Russell, and N. Y . Joly, Optica3, 952 (2016)

    work page 2016

  17. [17]

    P. D. Drummond and C. W. Gardiner, J. Phys. A 13, 2353 (1980)

    work page 1980

  18. [18]

    Bajer, J

    J. Bajer, J. Mod. Opt. 38, 1085 (1991)

    work page 1991

  19. [19]

    Felbinger, S

    T. Felbinger, S. Schiller, and J. Mlynek, Phys. Rev. Lett. 80, 492 (1998)

    work page 1998

  20. [20]

    Lindblad, Comm

    G. Lindblad, Comm. Math. Phys. 48, 119 (1976)

    work page 1976

  21. [21]

    A. J. Daley, Adv. Phys. 63, 77 (2014)

    work page 2014

  22. [22]

    C. W. Gardiner and P. Zoller, F oundations of Quantum Optics, 2nd ed. (Imperial College Press, 2014)

    work page 2014

  23. [23]

    L. I. Plimak, M. Fleischhauer, M. K. Olsen, and M. J. Collett, Phys. Rev. A 67, 013812 (2003)

    work page 2003

  24. [24]

    M. K. Olsen, L. I. Plimak, and M. Fleischhauer, Phys. Rev. A 65, 053806 (2002)

    work page 2002

  25. [25]

    P. D. Drummond and M. Hillery, The Quantum Theory of Non- linear Optics (Cambridge University Press, 2014)

    work page 2014

  26. [26]

    J. B. Surya, X. Guo, C.-L. Zou, and H. X. Tang, Optica 5, 103 (2018)

    work page 2018

  27. [27]

    J. B. Surya and H. X. Tang, (private communication) (2019)

    work page 2019

  28. [28]

    T. Ning, O. Hyvrinen, H. Pietarinen, T. Kaplas, M. Kauranen, and G. Genty, Opt. Express 21, 2012 (2013)

    work page 2012

  29. [29]

    Mateen, J

    F. Mateen, J. Boales, S. Erramilli, and P. Mohanty, Microsys. Nanoeng. 4, 14 (2018)

    work page 2018

  30. [30]

    Kischkat, S

    J. Kischkat, S. Peters, B. Gruska, M. Semtsiv, M. Chashnikova, M. Klinkmller, O. Fedosenko, S. Machulik, A. Aleksandrova, G. Monastyrskyi, Y . Flores, and W. T. Masselink, Applied Op- tics 51, 6789 (2012)

    work page 2012

  31. [31]

    Rackauckas and Q

    C. Rackauckas and Q. Nie, Discrete Continuous Dyn. Syst. Ser. B 22, 2731 (2016)

    work page 2016

  32. [32]

    Bezanson, A

    J. Bezanson, A. Edelman, S. Karpinski, and V . B. Shah, SIAM Review 59, 65 (2017)

    work page 2017

  33. [33]

    Rackauckas and Q

    C. Rackauckas and Q. Nie, J. Open Source Softw. 5, 15 (2017)

    work page 2017

  34. [34]

    M. K. Olsen and A. S. Bradley, Opt. Commun. 282, 3924 (2009)

    work page 2009

  35. [35]

    M. K. Olsen, Phys. Rev. A 70, 035801 (2004)

    work page 2004

  36. [36]

    C. W. Gardiner, Handbook of Stochasitc Methods , 4th ed. (Springer, 2009)

    work page 2009

  37. [37]

    Chaturvedi and C

    S. Chaturvedi and C. W. Gardiner, J. Stat. Phys.17, 429 (1977)

    work page 1977

  38. [38]

    D. F. Walls and G. J. Milburn, Quantum Optics (Springer, 1994)

    work page 1994

  39. [39]

    C. W. Gardiner and M. J. Collett, Phys. Rev. A31, 3761 (1985)

    work page 1985

  40. [40]

    C. M. Caves, Phys. Rev. D 23, 1693 (1981)

    work page 1981

  41. [41]

    D. F. Walls, Nature 306, 141 (1983)

    work page 1983

  42. [42]

    L.-M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 84, 2722 (2000)

    work page 2000

  43. [43]

    Simon, Phys

    R. Simon, Phys. Rev. Lett. 84, 2726 (2000)

    work page 2000

  44. [44]

    A. S. Bradley, M. K. Olsen, O. Pfister, and R. C. Pooser, Phys. Rev. A 72, 053805 (2005)

    work page 2005

  45. [45]

    Dalibard, Y

    J. Dalibard, Y . Castin, and K. Mølmer, Phys. Rev. Lett.68, 580 (1992)

    work page 1992

  46. [46]

    R. Dum, P. Zoller, and H. Ritsch, Phys. Rev. A45, 4879 (1992)

    work page 1992

  47. [47]

    Krmer, D

    S. Krmer, D. Plankensteiner, L. Ostermann, and H. Ritsch, Comput. Phys. Commun. 227, 109 (2018)

    work page 2018