{"id":"b632b8e7-a3ab-475c-947a-4cf45cc96073","arxiv_id":"2512.14174","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A Heisenberg-picture perturbative expansion adds beyond-semiclassical emitter-dynamics corrections to HHG theory, predicting squeezing that grows with emitter number while g(2) approaches 1.","lead":"This paper derives a perturbative Heisenberg-picture framework for high-order harmonic generation that includes quantum corrections to the emitting electron's motion, beyond the usual semiclassical treatment. It predicts that squeezing of the emitted light grows with the number of emitters while photon statistics become Poissonian, and that the new emitter-dynamics corrections significantly boost squeezing.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Perturbative control is not established: the expansion parameter N g0 |p~(omega)| can be order one at N=10^7, and the first O(g0^4 N^3) correction is admitted in App. D to be numerically uncontrolled.","rationale":"The paper develops a genuinely new Heisenberg-picture perturbative approach and derives closed-form expressions. The central claim that the PHD is 'accurately controlled' in experimental regimes depends on the smallness of epsilon = N g0 |p~(omega)|. The reader correctly identified this as the weakest assumption. In good faith, I do not see an internal inconsistency in the leading-order algebra; the concern is about the regime of validity. The manuscript itself flags the unresolved mode sum in App. D, which is exactly the missing truncation-error estimate. For N=10^7 and the stated bounds on g0 and |p~|, epsilon can exceed 1, so the formal condition (23) is not satisfied. The selection-rule argument in Sec. II B bounds only the spectrum and does not control the quadrature variance or g^(2), which are central to the headline claims. Therefore the appropriate verdict remains CONDITIONAL: the framework is promising and the leading-order results may be correct, but the paper should provide a concrete evaluation of the App. D mode sum (or an alternative rigorous bound) and state for which N the truncation is valid before the 'controlled/accurate' wording is accepted. I agree with the reader's weakest_assumption.","tokens_in":37315,"tokens_out":11745,"duration_ms":99030,"concrete_test":"Using the same Crank-Nicolson TDSE data as in Fig. 3, evaluate the App. D correction Eq. (D6) to the harmonic spectrum for a representative harmonic (e.g., the 9th) with a concrete mode discretization: mode spacing Delta_omega = 2*pi/T_pulse and cutoff at the highest harmonic plotted, then double the cutoff to test convergence. Compare the magnitude of the O(g0^4 N^3) correction to S_coh of Eq. (30) for N=10^5, 10^6, and 10^7. If the correction is not << S_coh or depends strongly on the cutoff, the perturbative truncation underlying the large-N squeezing predictions is not controlled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim that PHD is 'accurately controlled' rests on the inequality N g0 |p~(omega)| << 1 (Eq. (23)) justifying truncation of U' at second order. This is the load-bearing condition. It is not met in the regime where the headline N-scaling is shown: with g0=4e-8 and N=10^7, |p~(omega)| < 5 a.u. gives epsilon = N g0 |p~| up to ~2, not << 1; even at N=10^6, epsilon can be ~0.2. The authors' own App. D shows the first non-vanishing correction to the HHG spectrum is O(g0^4 N^3) and involves a sum over all field modes 'whose resolution and termination point are not well defined nor understood, making the evaluation uncontrolled.' The selection-rule argument bounding this term constrains only the spectrum; it provides no bound on corrections to the quadrature variance (Eq. (32)) or g^(2) (Eq. (33)) at the same order. Thus the perturbative control that underlies Eqs. (25)-(27) and the claimed beyond-semiclassical enhancement of squeezing is not established for the large-N results, even though the leading-order expressions may be correct.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Heisenberg-picture perturbative framework, termed PHD, for quantum-optical high-order harmonic generation. After displacing the coherent driving field and transforming to the semiclassical and free-field frames, the authors expand the interaction-picture evolution operator to second order in the effective coupling N g0 |p~(ω)|, obtaining closed-form expressions for the harmonic spectrum, quadrature squeezing, and g^(2)(0). The central new ingredient is the quantum correction to the semiclassical emitter dynamics, Eq. (25)-(27), which is claimed to significantly enhance the predicted squeezing. The framework is applied to an atomic ensemble and to an extended Fermi-Hubbard model, yielding N-scaling results: the coherent spectrum grows as N^2, the squeezing (variance) grows linearly in N, and g^(2)(0) tends to 1 for large N. The authors argue that the PHD is accurately controlled and preferable to Schrödinger-picture approaches, which they benchmark against.","tokens_in":37699,"tokens_out":3601,"duration_ms":37247,"significance":"If the central control claim held, the paper would provide a useful and computationally light tool for strong-field quantum optics, with clear N-scaling relations and an explicit link between photon statistics and emitter correlations. The derivation is largely self-contained and the appendices provide substantial detail, including exact factorizations of many-emitter expectation values and a comparison with the Schrödinger-picture product ansatz. The authors also explicitly identify a limitation in App. D, which is commendable. However, as detailed in the major comments, the load-bearing assertion that the expansion is accurately controlled in the experimentally relevant large-N regime is not established, and the paper's own App. D contains an admission that the leading correction cannot be evaluated. The qualitative predictions — increasing squeezing with N and Poissonian statistics for large N — may well be correct, but the quantitative 'controlled and accurate' claim requires either a rigorous remainder bound or a restriction to smaller N.","major_comments":[{"comment":"The perturbative control condition is Eq. (23), N g0 |p~(omega)| << 1. With g0 = 4e-8 a.u. and |p~(omega)| < 5 a.u., the product is about 2 at N = 1e7 and about 0.2 at N = 1e6. Yet the headline results in Fig. 3 are shown for N = 1e5, 1e6, and 1e7, including the largest value where the expansion parameter is not small. The claim that higher-order terms are negligible is therefore not supported in the very regime used to demonstrate the large-N scaling and the g^(2) -> 1 limit. The authors should either restrict the validity claim to N satisfying Eq. (23) with a numerical margin, or provide a quantitative estimate of the remainder beyond second order.","section":"Sec. II B 1, Eq. (23) and Sec. III B, Fig. 3"},{"comment":"The first non-vanishing correction to the harmonic spectrum is shown to be O(g0^4 N^3), and the authors write that the mode sum needed to evaluate it has a 'resolution and termination point [that] are not well defined nor understood, making the evaluation uncontrolled.' This is a direct admission that the leading correction to the central observable cannot be quantified. The main text (Sec. II B 2) nevertheless states that the truncation is justified by Eq. (23). The selection-rule argument that follows constrains only the spectral correction; it provides no bound on the O(g0^4 N^3) corrections to the quadrature variance in Eq. (32) or to g^(2) in Eq. (33), which are the quantities behind the claimed beyond-semiclassical enhancement of squeezing. The 'accurately controlled' characterization is therefore not established for the non-spectral observables.","section":"App. D, after Eq. (D7)"},{"comment":"The argument that the experimental observation of only odd harmonics 'proves that the coherent part of the spectrum is the dominating term, which puts both an upper and lower boundary on N' is an empirical consistency check, not a mathematical control of the expansion. It assumes that the O(g0^4 N^3) correction, if non-negligible, would necessarily produce detectable even-order contributions and that no other physical mechanism suppresses them. This does not replace a bound on the remainder, especially because the same argument cannot be applied to the quadrature variance or g^(2). The validity of the PHD in the large-N regime should be justified from the expansion itself, not from the experimental outcome it is meant to predict.","section":"Sec. II B 2"}],"minor_comments":[{"comment":"The quantity p~(omega) is used to state the expansion condition but is never defined precisely. Please specify whether it is the Fourier transform of the single-emitter dipole expectation value, its maximum over the pulse, or another quantity, and provide its value for the two model systems.","section":"Sec. II B 1, Eq. (23)"},{"comment":"There are several typographical slips: 'conferring with' should be 'cf.'; 'indeces' should be 'indices'; 'consitute' should be 'constitute'. These do not affect the results.","section":"App. A and App. B"},{"comment":"The estimate N ~ 6.3e6 for the Hubbard chains relies on several simplifying assumptions listed in the text. Since the paper identifies N as the key parameter controlling both the validity condition and the observable scalings, it would be helpful to give an uncertainty range for this estimate and to state explicitly whether the quoted N satisfies Eq. (23) for the Hubbard parameters.","section":"Sec. III C"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the uncontrolled higher-order correction admitted in App. D. I would ask the authors to either provide a genuine bound on the remainder or substantially weaken the 'accurately controlled' claim and restrict the quantitative conclusions to N values where Eq. (23) is safely satisfied. If no such bound can be provided, the paper may need to be reframed as a leading-order Heisenberg-picture calculation with an openly uncontrolled truncation, rather than a validated controlled expansion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper before citing it as a controlled Heisenberg-picture expansion for quantum-optical HHG. The core new thing is real: the authors derive beyond-semiclassical corrections to the emitter dynamics (their Q_q terms) that earlier Heisenberg treatments by Sundaram–Milonni and Stammer et al. did not include, and they work out closed-form expressions for the spectrum, squeezing, and g(2) with transparent N-scaling. The algebra in the appendices is self-contained and generally careful, and benchmarking against Schrödinger-picture results for the Hubbard model gives reasonable agreement for spectra. That is genuine progress.\n\nThe problem is the control claim. The whole expansion is justified by Eq. (23), N g0 |p~(ω)| << 1. With their values, g0 = 4e-8 and |p~| up to 5 a.u., at N=10^7 that product is about 2, and even at N=10^6 it is ~0.2. So the headline large-N results are outside the regime where the truncation is justified. Worse, Appendix D admits that the first O(g0^4 N^3) correction to the spectrum involves a sum over all modes \"whose resolution and termination point are not well defined nor understood, making the evaluation uncontrolled.\" The authors argue this term must vanish for inversion-symmetric targets because experiments see only odd harmonics, but that symmetry argument bounds the spectrum, not the quadrature variance or g(2), which get corrections at the same order. So the claimed significant enhancement of squeezing from beyond-semiclassical dynamics is not actually established — it may be an artifact of the truncation. The paper would need truncation-error estimates, or a better bound on the mode sum, before the word \"controlled\" is earned.\n\nOne more note: the appeal to experimental odd-harmonic observations as proof of validity is a consistency check, not an independent test. Still, the leading-order expressions for the spectrum are likely correct, and the formalism is a useful addition with advantages for deriving N-scaling and connecting photonic observables to emitter correlations.\n\nWho should read this: anyone working on quantum-optical HHG, especially people studying squeezing predictions. It deserves a serious referee — I would send it to review, but with a clear request to resolve the perturbative-control issue and provide a numerical estimate of the O(g0^4 N^3) correction to the observables.","headline":"Promising Heisenberg-picture treatment of quantum-optical HHG whose central control claim is not supported in the large-N regime; the squeezing enhancement may be an artifact of the truncation.","tokens_in":38134,"tokens_out":2835,"would_cite":true,"duration_ms":27119,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By expanding the dynamics in the Heisenberg picture, this paper derives closed-form quantum corrections to the emitters in high-order harmonic generation and shows these corrections can substantially increase the squeezing of the emitted li","keywords":["high-order harmonic generation","quantum optics","Heisenberg picture","squeezed light","photon statistics","strong-field physics","perturbative expansion","Fermi-Hubbard model"],"falsifier":"Compute the O(g0^4 N^3) correction to the harmonic spectrum of Appendix D using a convergent regularization of the mode sum; if for N=10^6–10^7 it is comparable to the coherent spectrum, the expansion is not controlled. Alternatively, measure quadrature squeezing and g(2)(0) from a gas jet as the emitter density is varied: the theory predicts squeezing increasing linearly with N and g(2) approaching 1, so observing g(2) clearly below 1 at large N, or squeezing not increasing with density, would contradict the central claim.","tokens_in":37248,"feed_emoji":"⚛️","tokens_out":5727,"duration_ms":46841,"temperature":0.7,"pith_summary":"This paper tries to establish a controlled, closed-form quantum-optical description of high-order harmonic generation (HHG) that goes beyond the usual semiclassical treatment of the emitting medium. The central claim is that the emitters' dynamics acquire small but important corrections from their coupling to the quantized electromagnetic field, and that these beyond-semiclassical corrections, captured by a perturbative Heisenberg-picture expansion, significantly change the predicted nonclassical properties of the emitted light. In particular, the paper claims the degree of squeezing grows with the number of independent emitters, while the photon statistics become Poissonian (classical) in the many-emitter limit—so squeezed light and classical photon counting can coexist. A sympathetic reader would care because this gives an experimentally accessible route to squeezing in the UV/XUV range and a direct link between a photonic observable and the underlying electron dynamics.","feed_headline":"Squeezing of high-harmonic light rises with emitter count","feed_subtitle":"Quantum back-action boosts harmonic squeezing as photon statistics turn classical","key_machinery":"The central object is the strong-field quantum-optical perturbative Heisenberg dynamics (PHD) expansion: after displacing the intense coherent driver into the Hamiltonian and moving to the rotating frames of both the semiclassical Hamiltonian and the free field, the time-evolution operator U'(t) is expanded to second order in the small effective coupling N g0 |p~(ω)|. This yields Eq. (25): Q'(t) = Q_sc(t) + Q_q(t), with the quantum correction Q_q(t) given by nested commutators of the interaction V(t) = A(t)·Q_sc(t) with the semiclassically driven dipole (Eq. (27)). The expansion turns the intractable joint light-matter dynamics into closed-form integrals over semiclassical two-time correlati","core_discovery":"The paper's key result is a set of equations—(25)–(27) of the paper—that express the exact emitter momentum Q'(t) in a displaced, rotating frame as the semiclassically driven response Q_sc(t) plus a quantum correction Q_q(t) built from commutators of the interaction with the driven emitter. Expanding the time-evolution operator to second order in the effective coupling N g0 |p~(ω)| yields closed expressions for the harmonic spectrum, the quadrature variance (squeezing), and the second-order correlation function g(2)(0). For the spectrum, the leading-order prediction coincides with the semiclassical result, so quantum corrections to the emitter do not alter the harmonic spectrum; for squeezin","pith_inferences":["One consequence the authors leave implicit: because squeezing grows with N while g(2) approaches 1, experiments that detect nonclassicality via photon statistics alone could miss the squeezing; quadrature measurements are the discriminating probe.","A testable extension would be to drive the same atomic model with a bright squeezed vacuum input; the PHD expansion could be adapted to check whether the predicted N-scaling of the nonclassical output survives when the input already carries squeezing.","The uncontrolled O(g0^4 N^3) mode sum identified in Appendix D suggests the formalism's accuracy claim is only as strong as a regularization of that sum; a future calculation supplying a convergent evaluation would either confirm or bound the perturbative regime.","The resonance-enhanced squeezing seen at the ninth harmonic in atoms and at the Mott exciton in the Hubbard model hints that resonantly tuned driving fields could amplify squeezing without raising N—an optimization the authors mention but do not pursue quantitatively."],"forward_implications":["To leading order, the HHG spectrum is unchanged by the beyond-semiclassical emitter dynamics, so standard semiclassical spectrum calculations remain valid in the regimes studied.","The degree of squeezing scales linearly with the number of independent emitters and is significantly enhanced by the quantum correction term, implying that larger phase-matched ensembles should produce more squeezed harmonic light.","In the many-emitter limit, the second-order correlation function g(2)(0) tends to 1, so the emitted light can be simultaneously squeezed and Poissonian in photon statistics.","The framework justifies the mode-decoupling product ansatz used in Schrödinger-picture treatments, since to leading order different harmonic modes do not couple.","The approach applies to both atomic gases and strongly correlated solids (Fermi-Hubbard chains), indicating that a strong electronic resonance boosts squeezing in both classes of emitters."],"fun_headline_variants":["Squeezing rises with emitter count in high-harmonic light","More emitters amplify squeezing of harmonic light","Quantum back-action boosts squeezing in HHG","Heisenberg approach reveals enhanced harmonic squeezing","Harmonic squeezing grows with number of emitters"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole calculation rests on the assumption that the effective light-matter coupling N g0 |p~(ω)| is much smaller than one; for the largest emitter numbers considered here (10^7) that product can approach one, and the paper itself notes the leading correction term cannot be evaluated because the required sum over all light modes is uncontrolled.","fun_headline_variants_meta":{"raw":{"variants":["Squeezing rises with emitter count in high-harmonic light","More emitters amplify squeezing of harmonic light","Quantum back-action boosts squeezing in HHG","Heisenberg approach reveals enhanced harmonic squeezing","Harmonic squeezing grows with number of emitters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000187,"raw_usage":{"total_tokens":1204,"prompt_tokens":824,"completion_tokens":380,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":568,"completion_tokens_details":{"reasoning_tokens":310}},"tokens_in":568,"tokens_out":380,"duration_ms":4416,"temperature":1.0,"reasoning_tokens":310,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T16:14:00.393510+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the O(g0^4 N^3) correction to the harmonic spectrum of Appendix D using a convergent regularization of the mode sum; if for N=10^6–10^7 it is comparable to the coherent spectrum, the expansion is not controlled. Alternatively, measure quadrature squeezing and g(2)(0) from a gas jet as the emitter density is varied: the theory predicts squeezing increasing linearly with N and g(2) approaching 1, so observing g(2) clearly below 1 at large N, or squeezing not increasing with density, would contradict the central claim.","supporting_citations":[],"review_version":1}