REVIEW 2 major objections 4 minor 61 references
Schr{\"o}dinger cat state formation in small bosonic Josephson junctions at finite temperatures and dissipation
T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Two-mode Bose condensates with attractive interactions form Schrödinger cat states only below a critical temperature set by the junction's Rabi frequency and interaction strength.
desk verdict The zero-temperature analysis is solid but the finite-temperature crossover claim is built on a misapplied metastable-decay formalism: in a symmetric double well the temperature that destroys cat coherence is set by the exponentially small tunneling splitting, not by the barrier-top frequency. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is an effective-potential mapping plus the imaginary-time bounce technique. The two-mode condensate is replaced by a fictitious macroscopic particle with effective mass $m_{\rm eff}=s\hbar\sqrt{1-x^2}/\omega_R$ moving in the potential $V(x)=-(s\hbar\omega_R/2)(\lambda x^2+2\sqrt{1-x^2})$, where $x=(n_b-n_a)/N$ is the normalized particle-number imbalance. For $\lambda>1$ this potential is a double well whose minima and barrier height are analytic functions of $\lambda$. The finite-temperature analysis inverts the potential to $U(x)=-V(x)$, introduces a 'thermon' quasiparticle that oscillates near the saddle point with frequency $\omega=\omega_R\sqrt{\lambda-1}$, and identifies the crossover temperature from the standard bounce relation $T_c=\hbar\omega/(2\pi k_B)$; Ohmic dissipation is folded in through the factor $\alpha$. That conversion of a many-body superposition problem into single-particle metastable decay is what carries the quantitative claims.
What would settle it
One could measure the escape or relaxation rate of a small attractive two-mode condensate as a function of temperature and check for a crossover from Arrhenius behavior $\Gamma\simeq(\omega_0/2\pi)e^{-V_0/k_BT}$ at $T>T_c$ to temperature-independent quantum tunneling at $T<T_c$, with the crossover temperature scaling as $\sqrt{\lambda-1}$ and vanishing at $\lambda=1$; a direct probe would be to record the Wigner function of the ground state and observe whether the two peaks and interference fringes survive below $T_c$ and wash out above it.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that the ground state of a two-mode BJJ with attractive interactions is governed by the parameter $\lambda$, with a quantum phase transition at $\lambda=1$: for $\lambda<1$ the state is a coherent Gaussian distribution over particle-number imbalance, while for $\lambda>1$ it becomes a superposition of two wave packets centered at $x_{\pm}=\pm\sqrt{1-\lambda^{-2}}$, and in the limit $\lambda\to\infty$ it approaches the N00N state $(|N,0\rangle+|0,N\rangle)/\sqrt{2}$. The finite-temperature extension shows that quantum tunneling, which is responsible for the cat state, dominates only for $T\le T_c=\hbar\omega/(2\pi k_B)$ with $\omega=\omega_R\sqrt{\lambda-1}$, and weak Ohmic dissipation renormalizes this to $T_c=\hbar\omega\alpha/(2\pi k_B)$ with $\alpha=\sqrt{1+(\gamma/2\omega)^2}-\gamma/2\omega$. The paper also derives the minimum bounce exponent $B_c=s\pi(\lambda-1)^{3/2}/(\alpha\lambda)$ for observing the cat state, and concludes that the required temperatures are low but experimentally reachable.
Load-bearing premise
The argument assumes that the single-particle imaginary-time bounce theory of metastable decay in an inverted potential correctly says when the many-body two-mode condensate actually forms or keeps its Schrödinger cat superposition, rather than merely describing escape from a metastable well.
Editorial extensions
If this is right
- For a given junction, cooling below $T_c=\hbar\omega_R\sqrt{\lambda-1}/(2\pi k_B)$ should suffice for cat-state formation, while above $T_c$ thermal activation destroys the superposition.
- The fidelity of the N00N-like state grows with $\lambda$, and the ideal N00N state is reached only as $\lambda\to\infty$, so experiments should target moderate $\lambda$ with a fidelity threshold rather than exact N00N states.
- Because $T_c$ vanishes at $\lambda=1$, the quantum phase transition itself is purely quantum; approaching it from the cat side requires ever-lower temperatures.
- Weak Ohmic dissipation lowers the crossover temperature by the factor $\alpha$, so stronger damping pushes cat-state formation to lower temperatures.
- For the $^7$Li parameters used, the temperature scale is $T_0=\hbar\omega_R/(2\pi k_B)=1.6$ nK, so $T_c=T_0\sqrt{\lambda-1}$ is far below the $\sim100$ nK condensation temperature and within reach of current cooling.
Reading between the lines
- A natural extension the paper does not develop is to turn the static crossover into a dynamical prediction: below $T_c$, the same bounce formalism could estimate the lifetime of the cat superposition under Ohmic noise, giving a direct decoherence time for practical preparation protocols.
- The Gaussian and inverted-parabola approximations break down near $\lambda=1$ and for $\lambda\to\infty$; a numerical solution of the two-mode Schrödinger equation with a thermal bath would test whether $T_c$ and $B_c$ survive beyond the analytic approximations.
- Because the transition at $T_c$ is described as second-order, one could look for critical slowing down or enhanced fluctuations in the particle-number imbalance near the crossover, similar to equilibrium phase transitions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies a two-mode bosonic Josephson junction (BJJ) with attractive interactions. It maps the Fock-space dynamics to a continuous Schrödinger equation for a fictitious particle in an effective potential, identifies a zero-temperature structural change at lambda = 1, uses Wigner functions to characterize cat and N00N-like states, and then invokes the imaginary-time bounce formalism to compute a crossover temperature Tc = hbar*omega_R*sqrt(lambda-1)/(2*pi*k_B) below which quantum tunneling is claimed to dominate and Schrödinger cat states can be observed in the presence of weak dissipation. The paper closes with numerical estimates for a 7Li condensate, finding Tc = 1.6 nK, well below the BEC temperature.
Significance. If the central finite-temperature claim were established, the paper would provide a simple quantitative temperature scale for preparing cat states in small BJJs, which would be useful for experiment planning. The paper has explicit analytic strengths: the effective-potential mapping, closed-form expressions for the barrier height and oscillation frequencies, and a Wigner-function representation of the cat-state interference are all given without fitted free parameters, and the lambda = 1 structural change is taken from prior literature rather than introduced ad hoc. However, the finite-temperature interpretation does not follow from the metastable-decay machinery used, and the quoted temperature is the wrong physical scale for coherence of the symmetric double-well ground state. As it stands, the paper's advertised main conclusion is not supported, although the zero-temperature analysis may have value as a separate contribution.
major comments (2)
- [Sec. 5, Eqs. (29)-(37)] The manuscript identifies Tc = hbar*omega/(2*pi*k_B) as the threshold below which Schrödinger cat states can be observed, but the bounce formalism (Refs. [53-57]) computes the crossover between thermal activation and quantum tunneling for metastable decay, not the temperature at which the equilibrium state of a symmetric double well is a coherent cat. For the effective potential (13) with lambda > 1, the two minima are degenerate, so the low-energy eigenstates are symmetric and antisymmetric superpositions split by an exponentially small energy Delta ~ hbar*omega_0*exp(-S0/hbar). The equilibrium off-diagonal coherence between the two localized configurations is proportional to tanh(Delta/(2*k_B*T)), so the coherence is already negligible once k_B*T >> Delta. The weak-metastability condition in Eq. (32), V0/(hbar*omega_0) >> 1, makes Delta exponentially suppressed, so Delta/k_B can be orders of magnitude smaller than Tc. For the Section 5 parameters (N=100, lambda approximately 2, omega_R=2*pi*208 Hz), the two-level estimate gives Delta/k_B far below the quoted Tc = 1.6 nK. The statement in Eq. (37) that SC state observation 'may be obtained for barriers possessing B >= Bc' therefore rests on an identification between an incoherent escape-rate crossover and the existence of a coherent superposition that the paper never derives. A direct calculation of the thermal-state fidelity or coherence in the two-mode model is needed before this claim can be accepted.
- [Sec. 5, numerical example after Eq. (37)] The numerical parameters are internally inconsistent. With uN/k_B = 10 nK and omega_R = 2*pi*208 Hz, the parameter lambda = uN/(hbar*omega_R) = (uN/k_B)/(hbar*omega_R/k_B) = 10/1.6 is about 6.25 at N = N0 = 100, not 'nearby the critical value lambda = 1' as stated in the text. To realize lambda approximately 1 with the same omega_R, one would need uN/k_B approximately 1.6 nK, which changes the quoted critical temperature and the associated experimental estimates. The numerical example should be reworked so that the stated values of N0, u, and omega_R are mutually consistent.
minor comments (4)
- [Eq. (19)] The printed formula sigma = 2*m_eff/hbar^2*c_2^{-1/4} is dimensionally inconsistent; the harmonic-oscillator ground-state width should be sigma = (hbar^2/(2*m_eff*c_2))^{1/4}. The subsequent explicit expressions in Eqs. (20)-(21) are consistent with the corrected form, so this appears to be a typographical error that should be fixed.
- [Sec. 2 and Sec. 6] The transition at lambda = 1 is repeatedly called a quantum phase transition, but for a finite-N system with N = 100 the ground-state energy is analytic in lambda; this is a finite-size crossover rather than a true phase transition. The terminology should be qualified.
- [Fig. 3 caption and Abstract] There are several typos: Fig. 3 caption contains 'other other parameters', and the Abstract contains 'in the the zero temperature limit'. These should be corrected.
- [Eq. (12) and surrounding text] The statement that the sqrt(1-x^2) factor in the effective mass 'plays no essential role' is not substantiated; at the minima for lambda > 1 this factor changes the effective mass by a factor lambda, which is explicitly used in Eq. (22). Please clarify the intended approximation.
Circularity Check
No circular reduction: Tc is obtained from the standard instanton crossover formula with the barrier frequency computed from U(x)=-V(x); the T≤Tc-to-SC identification is an unsupported physical assumption, not an input-output equivalence.
full rationale
The derivation chain is self-contained rather than circular. Starting from the two-mode Hamiltonian (1), the continuum limit Eq. (10) and the effective potential V(x) in Eq. (13) are algebraic; the λ=1 double-well onset follows from V''(0), and the Gaussian widths, ω0, and Wigner function are computed from the model parameters. The central finite-temperature result Tc=ℏω/(2πkB) (Eq. 29) is not fitted and is not defined to be the SC condition: it is the standard crossover temperature of the bounce/instanton formalism attributed to Refs. [53–57], and ω=ω_R√(λ-1) is obtained from U''_xx(0)/m of the inverted potential U(x)=-V(x) (Eqs. 26–30). The dissipation-modified Tc (Eq. 33) comes from the same external formalism. Numerics insert reported physical parameters (a_sc, ω⊥, ω_R) and do not reverse-engineer Tc. Self-citations to Refs. [21] and [42] provide background on the spectrum and the effective-potential/inverted-potential method; the load-bearing instanton formulas are external. The genuine weakness is physical, not circular: Eq. (37) asserts that 'the SC state observation, which is based on the essentially quantum tunneling processes, may be obtained for the barriers possessing B ≥ Bc,' but the paper never derives the relation between the metastable-decay crossover and the exponentially small tunneling splitting that controls equilibrium cat coherence in a symmetric double well. That is a missing justification for the central claim, not a reduction of the prediction to its inputs. Score 1 reflects minor, non-load-bearing self-citation only.
Assumptions & free parameters
assumptions (5)
- domain assumption Two-mode approximation is valid for the considered particle numbers (up to N ≈ 1000).
- domain assumption The discrete three-term recurrence (3) can be accurately replaced by the continuous Schrödinger equation (10) with position-dependent mass.
- ad hoc to paper The imaginary-time bounce formalism for a single quantum particle in a double-well potential applies to the many-body BJJ system and its effective potential.
- domain assumption Dissipation is Ohmic and weak (γ ≪ ω), entering only through the factor α in Eq. (33).
- domain assumption Gaussian wave-packet approximation (15) and harmonic expansion (17) are valid away from λ = 1.
invented entities (2)
-
Fictitious macroscopic (FM) quantum particle
-
Thermon quasiparticle
Cite this review
Pith. "Pith review of Schr{\"o}dinger cat state formation in small bosonic Josephson junctions at finite temperatures and dissipation." pith.science (2026). https://pith.science/paper/STQJXR3U
@misc{pith2026250716032,
author = {Pith},
title = {Pith review of: Schr\"odinger cat state formation in small bosonic Josephson junctions at finite temperatures and dissipation},
year = {2026},
howpublished = {\url{https://pith.science/paper/STQJXR3U}},
note = {Machine review of arXiv:2507.16032}
}
abstract
In this work, we consider the feasibility of Schr{\"o}dinger cat (SC) and $N00N$ states formation by a convenient bosonic Josephson junction (BJJ) system in two-mode approximation. Starting with purely quantum description of two-mode Bose-Einstein condensate we investigate the effective potential approach that provides an accurate analytical description for the system with a large number of particles. We show that in the zero temperature limit SC states result from a quantum phase transition that occurs when the nonlinear strength becomes comparable with the Josephson coupling parameter. The Wigner function approach demonstrates the growth of the SC state halves separation and formation of $N00N$-like states (a Fock state superposition) with the particle number increase. We examine the possibility to attain the SC state at finite temperatures and a weak dissipation leading to appearing of some critical temperature; it defines the second-order phase transition from classical activation process to the SC state formation through the quantum tunneling phenomenon. Numerical estimations demonstrate that the critical temperature is sufficiently below the temperature of atomic condensation. The results obtained may be useful for experimental observation of SC states with small condensate Josephson junctions.
Reference graph
Works this paper leans on
-
[42]
Lebedev M E et al 2017 Scientific reports 7 1
work page 2017
-
[1]
Schr¨ odinger E 1935 Naturwissenschaften, Bd.23 844
work page 1935
-
[2]
Wheeler J A and Zurek W 1983 Princeton University Press 260
work page 1983
-
[3]
Leggett A J 1984 Contemp. Phys. 25 583
work page 1984
-
[4]
Dodonov V V, Malkin I A and Man’ko V I 1974 Physica 72 597
work page 1974
-
[5]
Gerry C C and Knight P L 1997 American Journal of Physics 65 964
work page 1997
-
[6]
Haroche S and Raimond J 2006 Oxford University Press 82 86
work page 2006
-
[7]
Ourjoumtsev A et al 2007 Nature 448 784
work page 2007
Show all 61 references
-
[8]
Sychev D et al 2017 Nature Photon 11 382
2017
-
[9]
Grimm A et al 2020 Nature 584 205
2020
-
[10]
Cosacchi M et al M 2021 Phys. Rev. Research 3 023088
2021
-
[11]
Johnson K G et al 2017 Nat Commun 8 697
2017
-
[12]
Xu K et al 2020 Science Advances, 6 eaba4935
2020
-
[13]
Hacker B et al 2019 Nature Photon 13 110
2019
-
[14]
Eaton M, Nehra R and Pfister O 2019 New J. Phys. 21 113034
2019
-
[15]
Tatsuta M, Matsuzaki Y and Shimizu A 2019 Phys. Rev. A 100 032318
2019
-
[16]
Tsarev D V et al 2018 Optics express 26 19583
2018
-
[17]
Tsarev D V et al 2019 New Journal of Physics 21 083041
2019
-
[18]
Alodjants A et al 2022 Physical Review A 105 012606
2022
-
[19]
Mikheev E V, Pugin A S, Kuts D A, Podoshvedov S A and Nguyen B A 2019 Sci Rep 9 14301
2019
-
[20]
Takase K, Yoshikawa J, Asavanant W, Endo M and Furusawa A 2021 Phys. Rev. A 103 013710
2021
-
[21]
Tsarev D et al 2020 New Journal of Physics 22 113016 14
2020
-
[22]
Paw lowski K, Fadel M, Treutlein P, Castin Y and Sinatra A 2017 Phys. Rev. A 95 063609
2017
-
[23]
Jin R et al 1988 Appl. Phys. Lett. 53 1791
1988
-
[24]
Albiez M et al 2005 Phys. Rev. Lett. 95 010402
2005
-
[25]
Vretenar M, Kassenberg B, Bissesar S, Toebes C and Klaers J 2021 Phys. Rev. Research 3 023167
2021
-
[26]
Devoret M H, Wallraff A, Martinis J M 2004 arXiv preprint cond-mat/0411174
2004 arXiv
-
[27]
Cirac J I, Lewenstein M, Mølmer K and Zoller P 2018 Phys. Rev. A 57 1208
2018
-
[28]
Mazzarella G, Salasnich L, Parola A and Toigo F 2011 Phys. Rev. A 83 053607
2011
-
[29]
Piazza F, Pezze’ L and Smerzi A 2008 Phys. Rev. A 78 051601
2008
-
[30]
Haigh T J, Ferris A J and Olsen M K 2010 Optics communications 283 3540
2010
-
[31]
Bollinger J J, Itano W M, Wineland D J and Heinzen D J 1996 Phys. Rev. A 54 R4649
1996
-
[32]
Boto A N, Kok P, Abrams D S, Braunstein S L, Williams C P and Dowling J P 2000 Phys. Rev. Lett. 85 2733
2000
-
[33]
B: Quantum Semiclass
Kok P, Braunstein S L and Dowling J P 2004 Opt. B: Quantum Semiclass. Opt. 6 S811
2004
-
[34]
Dowling J P 2008 Contemporary physics 49 125
2008
-
[35]
Burlakov A V, Chekhova M V, Karabutova O A, Klyshko D N and Kulik S P 1999 Phys. Rev. A 60 R4209(R)
1999
-
[36]
Yurke B and Stoler D 1986 Phys. Rev. Lett. 57 13
1986
-
[37]
Chudnovsky E M and Garanin D A 1997 Phys. Rev. Lett. 79 4469
1997
-
[38]
Anglin J R and A Vardi 2001 Phys. Rev. A 64 013605
2001
-
[39]
Zaslavskii O 1990 Phys. Lett. A 145 471
1990
-
[40]
Heiss W D 1988 Zeitschrift f¨ ur Physik A Atomic Nuclei 329 133
1988
-
[41]
Glick A J, Lipkin H J and Meshkov N 1965 Nuclear Physics 62 211
1965
-
[43]
Caldeira A O and Leggett A J 1981 Phys. Rev. Letts 46 211 15
1981
-
[44]
Grabert H, Weiss U and Hanggi P 1984 Phys. Rev. Letts 52 2193
1984
-
[45]
Larkin A I and Ovchinnikov Y 1984 Zh. Eksp. Teor. Fiz. 86 726
1984
-
[46]
Leggett A J 2001 Rev. Mod. Phys. 73 307
2001
-
[47]
Pezze L et al 2018 Rev. Mod. Phys. 90 035005
2018
-
[48]
Smerzi A et al 1997 Phys. Rev. Letts 79 4950
1997
-
[49]
Kohler S and Sols F 2002 Phys. Rev. Letts 89 060403
2002
-
[50]
Rubeni D, Foerster A, Mattei E and Roditi I 2012 Nucl. Phys. B 856 698
2012
-
[51]
Hillery M et al 1984 Phys. Rep. 106 121
1984
-
[52]
Mukunda N 1978 Pramana 11 1
1978
-
[53]
Langer S 1967 Ann. Phys. (N.Y.) 41 108
1967
-
[54]
Callan C G and Coleman S 1977 Physical Review D 16 1762
1977
-
[55]
Mel’nikov V I and Meshkov S V 1983 ZhETF Pisma Redaktsiiu 38 111
1983
-
[56]
Hanggi P, Grabert H, Ingold G L and Weiss U 1985 Phys. Rev. Letts. 55 761
1985
-
[57]
Grabert H, Olschowski P and Weiss U 1987 Phys. Rev. B. 36 1931
1987
-
[58]
Khaykovich L et al 2002 Science 296 1290
2002
-
[59]
Parkins A S and Walls D F 1998 Phys. Rep. 303 1
1998
-
[60]
Bradley C C et al 1955 Phys. Rev. letts 75 1687
1955
-
[61]
Bradley C C, Sackett C A, Hulet R G 1997 Phys. Rev. Letts. 78 985 16
1997
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.