REVIEW 4 major objections 5 minor 24 references
Fractional squeezing: spectra and dynamics from generalized squeezing Hamiltonian with fractional orders
T0 review · 4 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read By treating the generalized-squeezing order n as a continuous variable, this paper argues that the squeezing Hamiltonian's spectrum switches from continuous to discrete exactly at n=2, and that the oscillation amplitude switches from diverg
desk verdict A clever fractional-order interpolation of the squeezing Hamiltonian with finite-size scaling identifies n=2 and n=4 as critical points, but the conclusions may depend on the interpolation, which is not tested. 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 central object is the fractional-order squeezing Hamiltonian matrix H(n) whose off-diagonal elements are sqrt(Gamma((k+1)n+1)/Gamma(kn+1)), obtained by replacing factorials with Gamma functions so that n becomes a continuous parameter. The extrapolation machinery consists of fitting the truncation-size dependence of the smallest positive eigenvalue to E_min = E_min,∞ + δE_min N^{-α} and the renormalized photon number to ⟨m⟩ = A N^{-α} + B, which lets the paper read off the infinite-N limits and the scaling exponents. A secondary explanatory device is a toy 'hierarchical' Hamiltonian with rapidly increasing off-diagonal elements, which shows why the middle-spectrum eigenstates become loca
What would settle it
Compute the smallest positive eigenvalue E_min for n=2 using much larger truncation sizes than those reported; if the extrapolated E_min tends to a positive finite value rather than zero, the claimed continuous spectrum at n=2 is wrong. Likewise, at n=4, compute ⟨m⟩ for truncation sizes beyond the range used here: if the growth is faster than logarithmic (e.g., a small power law), the claimed logarithmic scaling and the n=4 critical boundary would be refuted.
Extended reading notes
Core claim
The author generalizes the generalized-squeezing Hamiltonian to fractional order by replacing the factorial ratios in the matrix elements with Gamma-function ratios, so that the squeezing order n can be any positive real number. Using finite-N truncations and power-law extrapolations of the smallest positive eigenvalue and the renormalized photon number, the paper identifies two critical points: n=2, where the spectrum changes from continuous to discrete, and n=4, where the photon-number oscillation amplitude changes from asymptotically infinite to finite. At n=4 the data follows a logarithmic scaling law, which is why the power-law fit fails there; for n>4 the extrapolated photon number app
Load-bearing premise
The argument rests on trusting finite-N truncations and power-law fits to extrapolate to the infinite-N limit near n=2 and n=4, even though fractional n lacks a physical multiphoton interpretation and the fit at n=4 fails, as the paper itself acknowledges.
Editorial extensions
If this is right
- The spectrum of two-photon squeezing, previously debated in finite-N simulations, is continuous, consistent with the standard continuous-spectrum treatment of two-photon squeezing dynamics.
- At n=4 the photon-number oscillation amplitude diverges only logarithmically, so it is infinite in the infinite-size limit but weaker than any power-law divergence.
- For n>4 the photon-number oscillation amplitude is finite, and for n≥5 it stabilizes close to the value 0.5 predicted by a two-level truncation.
- The fractional-order extrapolation method can predict behavior at computationally hard integer points, offering a template for studying other parameter-dependent quantum-optical Hamiltonians such as the multiphoton quantum Rabi model.
- In the large-n limit, the low-energy physics of generalized squeezing is governed almost entirely by the two lowest Fock states, giving a simple asymptotic description.
Reading between the lines
- If the Gamma-function interpolation is more than a numerical device, it may be possible to define fractional powers of the creation and annihilation operators algebraically, and to check whether the resulting spectra match physical implementations of nonlinear interactions.
- The linear scaling of the exponent α on both sides of n=2 hints at a closed-form scaling law for the small-n discrete-to-continuous transition that the paper does not derive; extracting that law from the fit data could be a direct follow-up.
- The same finite-N fitting strategy could be applied to the multiphoton quantum Rabi model to locate its critical coupling thresholds, which the paper mentions only as a suggestion.
- The parity-sensitivity of dynamics for n>4, explained by the hierarchical toy model, suggests that experiments probing high-order squeezing may need to control the truncation parity of the effective Hilbert space, an implication the paper leaves implicit.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript generalizes the generalized-squeezing Hamiltonian H_n = i[(a^†)^n - a^n] to real (fractional) orders n by replacing the factorial matrix elements in Eq. (4) with Gamma functions, giving Eq. (5). The author performs finite-N truncation diagonalizations for n ∈ [0,6] and fits the lowest positive eigenvalue E_min and the renormalized photon-number expectation value ⟨m⟩ to power-law-plus-offset forms (Eqs. (7) and (9)). From these fits, the paper infers that the spectrum is continuous for n ≤ 2 and discrete for n > 2, that the generalized-squeezing oscillation amplitude diverges for n < 4 and becomes finite for n ≥ 4, with n = 4 exhibiting logarithmic scaling, and that large-n behavior is captured by a hierarchical toy Hamiltonian. The central claim is that fractional-order calculations provide a reliable way to locate and characterize the critical points n = 2 and n = 4, which are hard to treat by direct integer-order simulation.
Significance. If the inferred critical points are robust, this would provide a clean global picture of generalized squeezing and resolve ambiguities left by integer-order studies, e.g., the continuity of the spectrum at n = 2 and the divergence behavior at n = 4. The paper is honest in emphasizing that fractional n is a mathematical tool rather than a physical parameter (Sec. II) and in acknowledging the fit failure at n = 4 (Sec. III.B). The Gamma-function interpolation is explicit and reproducible from the equations, and the hierarchical toy model gives a simple, testable explanation of the large-n trend. However, the central inference rests on a non-unique interpolation and on three-parameter fits without error bars at exactly the points where the fits are least reliable. The manuscript does not yet establish that the critical values n = 2 and n = 4 are independent of the interpolation choice or that the finite-size extrapolations are converged.
major comments (4)
- [Sec. II, Eq. (5)] The continuous family Ĥ(n) is defined by replacing the factorial matrix elements of Eq. (4) with Gamma functions. This is one of infinitely many smooth extensions that agree at integer n, and the paper itself states that fractional n has no physical interpretation. The inferred critical points n = 2 and n = 4 and the associated scaling exponents could, in principle, depend on the specific interpolation: the leading large-k growth h_k ~ (kn)^{n/2} is interpolation-independent in a broad class of extensions, but sub-leading terms can shift the thresholds. Because the paper's conclusions about the physical integer cases are drawn from the fractional-n data, please provide a robustness check under alternative factorial interpolations (e.g., Γ(n+1) times a slowly varying or subexponential factor) or an analytic derivation of the thresholds from an interpolation-independent quantity. Without
- [Sec. III.A, Eq. (7)] The conclusion that E_min,∞ = 0 for n ≤ 2 is obtained by fitting E_min(N) to E_min,∞ + δE_min N^{-α}. In the critical region near n = 2, the fitted α deviates from the extrapolated linear behavior, and the paper attributes this to 'likely' numerical errors without supporting statistics. Since E_min,∞ is itself the fitted intercept, the zero/nonzero status of the asymptotic eigenvalue at n = 2 is read off from the same fitting function used to define it. Please provide convergence evidence that is not based solely on this fit: for example, include odd-N data, Richardson extrapolation, or a rigorous upper/lower bound on E_min(N), and report confidence intervals for the fitted parameters. The current graphical evidence is suggestive but not conclusive at the critical point.
- [Sec. III.B, Eq. (9)] The fit ⟨m⟩ = A N^{-α} + B is explicitly stated to fail at n = 4 ('the fitting fails and generates an error in our numerical fitting calculations'). This is precisely the value of n at which the paper locates the transition from divergent to finite ⟨m⟩. The logarithmic scaling at n = 4 is taken from Ref. [17] rather than demonstrated by the present fitting procedure. Consequently, the boundary n = 4 and the claim that α crosses zero at n = 4 are inferred indirectly from fits that are unreliable exactly at the point of interest. Please add a dedicated analysis at n = 4 (e.g., a direct logarithmic fit with residuals), and quantify the uncertainty in the inferred transition point, for instance by bootstrapping over truncation-size ranges or by comparing fits with and without the n = 4 data.
- [Sec. III, Figs. 2 and 5] All load-bearing conclusions rely on three-parameter nonlinear fits (Eqs. (7) and (9)) to truncation-size data. No error bars, residual analyses, or goodness-of-fit metrics are reported, and only even values of N are used. This is especially concerning because the fit models are used to locate critical points at which the models themselves break down. Even if the qualitative picture is correct, the paper would be substantially strengthened by reporting uncertainties on E_min,∞, α, A, and B and by showing that the results are stable under removing the smallest few N data points.
minor comments (5)
- [Sec. III.B] Typo: 'intuitively undestood' should be 'intuitively understood'; 'logarithim scaling law' should be 'logarithmic scaling law'.
- [Figure 2 caption] The inset and axis labels contain apparent rendering issues: '10□5' and '10□3' should presumably be 10^{-5} and 10^{-3} or similar. Please check the figure output.
- [Sec. IV] The sentence 'We can therefore restrict ourselves to the state space extending from 1 to N−2' after removing the two highest states is slightly confusing about indexing; clarify whether states are numbered 1,...,N and which indices remain.
- [Sec. II] The paragraph on fractional calculus and fractional spatial dimensions is useful context, but it is not connected to the rest of the paper. A sentence explaining why these analogies are not pursued further would help set expectations.
- [Sec. III.A, Eq. (7)] In Eq. (7), the notation δE_min is used as a fitting parameter, which is nonstandard and could be confused with a small error. Consider renaming it c or E_1 to avoid ambiguity.
Circularity Check
No significant circularity: the fractional-n construction is an independent interpolation, and the extrapolated critical points are inferred from numerical fits rather than inserted by definition.
full rationale
The paper's central move is to replace factorial matrix elements in Eq. (4) with Gamma functions in Eq. (5), producing a continuous family of truncated Hamiltonians. The integer-n results are not used to define this family: Eq. (5) is an independent extension that agrees with the integer cases but is not constrained by the conclusions (continuous spectrum for n<=2, finite oscillation amplitude for n>=4). The subsequent fits, Eqs. (7) and (9), estimate asymptotic intercepts and scaling exponents from finite-N data; the inferred critical values n=2 and n=4 are extrapolations of those fitted curves, not quantities already contained in the fit inputs. The paper explicitly flags the conceptual limitation that fractional n has no physical interpretation (Sec. II after Eq. (5)); this is a robustness concern about interpolation nonuniqueness, not circularity. The citations to Refs. [17] and [19] are to prior published work by the same group, but they are used as context and benchmarks (e.g., the n=4 logarithmic scaling in Ref. [17]) rather than as unverified uniqueness theorems that force the conclusions; the present paper's own finite-N data independently generates the fitted trends. No equation in the paper reduces to another by construction, and no fitted parameter is renamed as a prediction. The derivation chain is self-contained given the Gamma-interpolation ansatz; whether that ansatz is physically privileged is a validity question, not a circularity.
Assumptions & free parameters
free parameters (4)
- E_min,∞ (per-n asymptotic smallest eigenvalue) =
≈0 for n < 2; rises rapidly for n > 2
- α (per-n convergence exponent for E_min) =
Decreases linearly to 0 at n=2 for n<2; increases linearly for n>2
- A, B, α in the ⟨m⟩ scaling fit =
B→0.5 as n→∞; α changes sign at n=4
- Toy-model hierarchy ratio h_{j+1}/h_j =
10 or 100
assumptions (5)
- ad hoc to paper Gamma-function interpolation of the squeezing-order Hamiltonian (replacing factorials by Γ) defines a meaningful continuous family interpolating integer squeezing orders.
- domain assumption Finite-size extrapolation using Eq. (7) correctly gives the infinite-N asymptotic E_min, with a power-law ansatz.
- domain assumption The two eigenstates closest to zero characterize the squeezed-vacuum dynamics and spectrum for n>2.
- domain assumption Even truncation sizes N are sufficient; parity effects do not change the inferred boundaries.
- standard math Standard finite-dimensional linear algebra and numerical diagonalization are valid for the truncated matrix.
invented entities (1)
-
Fractional-order squeezing Hamiltonian Ĥ(n) with Gamma-function matrix elements
Cite this review
Pith. "Pith review of Fractional squeezing: spectra and dynamics from generalized squeezing Hamiltonian with fractional orders." pith.science (2026). https://pith.science/paper/YHOYXEV5
@misc{pith2026260115693,
author = {Pith},
title = {Pith review of: Fractional squeezing: spectra and dynamics from generalized squeezing Hamiltonian with fractional orders},
year = {2026},
howpublished = {\url{https://pith.science/paper/YHOYXEV5}},
note = {Machine review of arXiv:2601.15693}
}
abstract
We generalize the generalized-squeezing problem to include fractional values of the squeezing order $n$. This approach allows us to determine the locations of critical points at which qualitative changes in behaviour occur and accurately predict the behaviour at these critical points, which are challenging for conventional computational methods. Based on our numerical calculations, we identify with a high degree of confidence the point at which the spectrum turns from continuous to discrete and the point at which oscillations turn from having asymptotically infinite amplitudes to having finite amplitudes. Furthermore, we numerically investigate the behaviour in the large $n$ regime and provide an intuitive explanation for the numerical results.
Figures
Reference graph
Works this paper leans on
-
[17]
Ashhab and M
S. Ashhab and M. Ayyash, Properties and dynamics of generalized squeezed states, New J. Phys.27, 054104 (2025)
2025
-
[1]
D. F. Walls and G. J. Milburn,Quantum Optics(Springer, Berlin, 1994)
1994
-
[2]
M. O. Scully and M. S. Zubairy,Quantum Optics(Cambridge University Press, 1997). 13
1997
-
[3]
P. D. Drummond and M. Hillery,The Quantum Theory of Nonlinear Optics(Cambridge University Press, 2014)
2014
-
[4]
Felicetti, D
S. Felicetti, D. Z. Rossatto, E. Rico, E. Solano, and P. Forn-D ´ ıaz, Two-photon quantum Rabi model with superconducting circuits, Phys. Rev. A97, 013851 (2018)
2018
-
[5]
C. W. S. Chang, C. Sab ´ ın, P. Forn-D ´ ıaz, F. Quijandr ´ ıa, A. M. Vadiraj, I. Nsanzineza, G. Johansson, and C. M. Wilson, Observation of three-photon spontaneous parametric down- conversion in a superconducting parametric cavity, Phys. Rev. X10, 011011 (2020)
2020
-
[6]
A. M. Eriksson, T. S´ epulcre, M. Kervinen, T. Hillmann, M. Kudra, S. Dupouy, Y. Lu, M. Khanahmadi, J. Yang, C. Castillo-Moreno, P. Delsing, and S. Gasparinetti, Universal control of a bosonic mode via drive-activated native cubic interactions, Nature Commun.15, 2512 (2024)
2024
-
[7]
Ayyash, X
M. Ayyash, X. Xu, S. Ashhab, and M. Mariantoni, Driven multiphoton qubit-resonator inter- actions, Phys. Rev. A110, 053711 (2024)
2024
Show all 24 references
-
[8]
B˘ az˘ avan, S
O. B˘ az˘ avan, S. Saner, D. J. Webb, E. M. Ainley, P. Drmota, D. P. Nadlinger, G. Araneda, D. M. Lucas, C. J. Ballance, and R. Srinivas, Squeezing, trisqueezing, and quadsqueezing in a spin-oscillator system, arXiv:2403.05471
-
[9]
Saner, O
S. Saner, O. B˘ az˘ avan, D. J. Webb, G. Araneda, D. M. Lucas, C. J. Ballance, and R. Srini- vas, Generating arbitrary superpositions of nonclassical quantum harmonic oscillator states, arXiv:2409.03482
-
[10]
R. A. Fisher, M. M. Nieto and V. D. Sandberg, Impossibility of naively generalizing squeezed coherent states, Phys. Rev. D29, 1107 (1984)
1984
-
[11]
Hillery, M
M. Hillery, M. S. Zubairy, and K. W´ odkiewicz, Squeezing in higher order nonlinear optical processes, Phys. Lett. A103, 259 (1984)
1984
-
[12]
C. K. Hong and L. Mandel, Higher-order squeezing of a quantum field, Phys. Rev. Lett.54, 323 (1985)
1985
-
[13]
S. L. Braunstein and R. I. McLachlan, Generalized squeezing, Phys. Rev. A35, 1659 (1987)
1987
-
[14]
S. L. Braunstein and C. M. Caves, Phase and homodyne statistics of generalized squeezed states, Phys. Rev. A42, 4115 (1990)
1990
-
[15]
Hillery, Photon number divergence in the quantum theory of n-photon down conversion, Phys
M. Hillery, Photon number divergence in the quantum theory of n-photon down conversion, Phys. Rev. A42, 498 (1990)
1990
-
[16]
Braak, Thek-photon quantum Rabi model, to appear in T
D. Braak, Thek-photon quantum Rabi model, to appear in T. Takagiet al.(Eds.), 14 Mathematical Foundations for Post-Quantum Cryptography(Springer, Singapore, 2026); arXiv:2401.02370
2026 arXiv
-
[18]
Properties and dynamics of generalized squeezed states
R. Gordillo and R. Puebla, Comment on “Properties and dynamics of generalized squeezed states”, arXiv:2507.12250
-
[19]
Ashhab, F
S. Ashhab, F. Fischer, D. Lonigro, D. Braak, and D. Burgarth, Finite-dimensional approxi- mations of generalized squeezing, Phys. Rev. A113, 013703 (2026)
2026
-
[20]
Fischer, D
F. Fischer, D. Burgarth, and D. Lonigro, Self-adjoint realizations of higher-order squeezing operators, arXiv:2508.09044
-
[21]
Properties and dynamics of generalized squeezed states
S. Ashhab and M. Ayyash, Reply to Comment on “Properties and dynamics of generalized squeezed states”, arXiv:2601.07227
-
[22]
Laskin, Fractional Quantum Mechanics, Phys
N. Laskin, Fractional Quantum Mechanics, Phys. Rev. E62, 3135 (2000)
2000
-
[23]
Herrmann,Fractional Calculus: An Introduction for Physicists(World Scientific, Singa- pore, 2011)
R. Herrmann,Fractional Calculus: An Introduction for Physicists(World Scientific, Singa- pore, 2011)
2011
-
[24]
Goldenfeld,Lectures On Phase Transitions And The Renormalization Group(Taylor & Francis, Boca Raton, 1992)
N. Goldenfeld,Lectures On Phase Transitions And The Renormalization Group(Taylor & Francis, Boca Raton, 1992). 15
1992
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.