{"id":"c638a6a1-007d-41c7-a23b-c5fcee46638f","arxiv_id":"2501.08051","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new Møller-Plesset perturbation theory built on strong-coupling QED Hartree-Fock orbitals accurately captures cavity-induced electron-photon correlation and avoids the long-range artifacts of alternative QED-MP2 approaches.","lead":"This paper develops SC-QED-MP2, a new perturbative method for molecules inside optical cavities that starts from orbitals already dressed by the cavity field. It shows this method tracks high-level QED coupled cluster benchmarks more closely than earlier MP2 variants and avoids spurious long-range effects in molecular interactions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The long-range part of the central claim is not directly benchmarked: the intermolecular dissociation plots (Figs. 3-6) contain no QED-CCSD reference, so SC-QED-MP2's well-behaved plateau is not shown to be the physically correct ground-state curve.","rationale":"The reader's weakest assumption concerns QED-CCSD adequacy. I agree this is the right area but would sharpen it: the adequacy issue matters most where QED-CCSD is absent, i.e., the intermolecular long-range curves. For the single-molecule tests, QED-HF-based QED-CCSD is a defensible benchmark because the systems are neutral and isolated; the paper's critique of QED-HF (origin dependence, non-size-intensivity) does not invalidate those comparisons. The genuinely load-bearing gap is that the one claim that distinguishes SC-QED-MP2 from the simpler QED(np-HF)-MP2—the behavior at long range—is presented without any reference data. The theoretical size-intensivity analysis is plausible and well argued, but numerical verification against an independent method is needed to show that the SC-QED-MP2 plateau is the physical curve rather than just a non-divergent one. I also note the paper's honesty about qualitative trends and the non-variational nature of the methods, which lowers the risk of overclaiming. Overall, the derivation is careful and the method is promising; the missing reference curve is an addressable gap, hence CONDITIONAL remains the right verdict.","tokens_in":22703,"tokens_out":31438,"duration_ms":306714,"concrete_test":"Compute QED-CCSD (with the same QED-HF reference used in the paper) dissociation curves for the H2 dimer (Figure 3), water dimer (Figure 4), and benzene-water (Figure 5) systems at the same geometries, cavity parameters, and aug-cc-pVDZ basis, and overlay them with the SC-QED-MP2 curves. Settling criterion: if QED-CCSD matches SC-QED-MP2's plateau and does not exhibit the QED-MP2/LF-MP2 divergence, the concern is resolved; if QED-CCSD diverges or disagrees materially, the central claim that SC-QED-MP2 provides the physically correct long-range behavior is not established by the data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central differentiator is that SC-QED-MP2 avoids the unphysical long-range behavior of QED-MP2 and LF-MP2. This is demonstrated in Figures 3-6 by showing that SC-QED-MP2 (and QED(np-HF)-MP2) reaches a plateau while QED-MP2/LF-MP2 diverge. However, no QED-CCSD reference is shown for any intermolecular curve. The only benchmark used in the paper, QED-CCSD, is built on QED-HF orbitals, which the same paper argues are non-size-intensive and, for charged systems, origin dependent (Section 2, Eqs. 14-16). For the single-molecule tests (coupling/frequency dispersions, orientational effects) this is not a serious issue because the molecules are neutral and isolated. But for the long-range intermolecular regime, the absence of a reference means the central claim rests solely on the theoretical size-intensivity argument. That argument (Eq. 24) explains why QED-MP2 diverges as 1/R^2, and the SC-QED-HF Fock size-intensivity (Ref 45) suggests SC-QED-MP2 should not. This is credible, but it does not establish that the SC-QED-MP2 plateau is the exact or QED-CCSD-level long-range energy. The paper itself notes that QED-CCSD built on QED-HF 'is expected to capture more correlation energy than SC-QED-MP2' and that an SC version of QED-CC is under development. Without a long-range reference, the ranking of SC-QED-MP2 versus its competitors in the regime that most distinguishes them is not numerically verified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a second-order Møller-Plesset perturbation theory built on the strong-coupling QED Hartree-Fock reference (SC-QED-MP2). The authors partition the Pauli-Fierz Hamiltonian after an orbital-dependent coherent-state transformation, define a zeroth-order Fock Hamiltonian in the dipole basis, and derive the second-order energy correction as sums over double electronic excitations with arbitrary photon number, single excitations with at least one photon, and purely photonic excitations with at least two photons. They compare SC-QED-MP2 with QED-MP2, QED(np-HF)-MP2, LF-MP2, and QED-CCSD on coupling and frequency dispersions for ammonia, intermolecular dissociation curves for the hydrogen dimer, the water dimer, and benzene-water, and polarization orientation scans for chloroethylene and water. The central claims are that SC-QED-MP2 accurately reproduces QED-CCSD electron-photon correlation while remaining affordable, and that, unlike QED-MP2 and LF-MP2, it does not display unphysical long-range intermolecular behavior because its Fock operator and orbitals are size-intensive.","tokens_in":23078,"tokens_out":11020,"duration_ms":108889,"significance":"The proposed method is a natural and potentially useful extension of SC-QED-HF: it adds perturbative correlation while preserving the variational orbital-specific coherent-state reference. The derivation in the Supporting Information is systematic and does not rely on fitted parameters; the {eta_p} parameters are variationally optimized. The authors provide a transparent scaling argument (Eq. 24) for the QED-MP2 long-range artifact and identify the basis-dependent origin of the LF-MP2 problem (Eq. 57). The data are deposited at a persistent DOI, and the calculations are reproducible in principle. These are real strengths. However, the numerical validation is incomplete exactly in the regime that distinguishes the method: the intermolecular long-range curves are not benchmarked against an independent reliable reference, and the only coupled-cluster benchmark used is built on the QED-HF reference that the paper itself criticizes. The central accuracy claim is therefore plausible but not yet fully established.","major_comments":[{"comment":"The paper's headline differentiator is that SC-QED-MP2 avoids the unphysical long-range behavior of QED-MP2 and LF-MP2, but this is never checked against an accurate reference. The dissociation curves contain only the perturbative methods, and the text explicitly notes that a SC version of QED-CC is under development. A plateau relative to QED-MP2 is not enough to show that the plateau is the correct ground-state energy; it could be a wrong but well-behaved limit. I request at least one intermolecular curve with a QED-CCSD (or QED-FCI for a small model) reference, or an equivalent independent benchmark, to support the claim.","section":"Section 3, Figures 3-6"},{"comment":"The benchmark QED-CCSD is built on QED-HF orbitals, and Section 2 (Eqs. 13-16) argues that the QED-HF Fock operator is non-size-intensive and origin-dependent for charged systems. For neutral single molecules this is a reasonable benchmark, but for the long-range intermolecular regime the reference itself may inherit the same artifact. The statement 'the comparison is justified as we focus on electron-photon correlation effects' is qualitative; a numerical demonstration that QED-CCSD's long-range interaction energy is stable is needed.","section":"Section 3, first paragraph"},{"comment":"The size-intensivity claim for SC-QED-MP2 is carried over from the SC-QED-HF Fock matrix (Ref. 45), but the second-order energy in Eq. (49) is an infinite sum over photonic excitations, and its size-intensivity is not demonstrated analytically or numerically. An explicit argument, or a numerical check that the truncated energy is additive for separated subsystems, would close this gap and directly support the long-range claim.","section":"Section 2, Eqs. (41)-(49)"}],"minor_comments":[{"comment":"'Pauli-Fiertz' should be 'Pauli-Fierz'.","section":"Section 2, Eq. (1)"},{"comment":"The sentence 'Specifically, the method are built starting from two possible reference states' contains a subject-verb agreement error; it should be 'the methods are built'.","section":"Section 1, paragraph 3"},{"comment":"The sentence 'QED-HF is unable to account for the cavity-induced non size-extensive effects' is confusing because QED-HF was just called size-extensive; the intended term is likely 'non-size-intensive effects'.","section":"Section 2, QED-HF discussion"},{"comment":"The offset procedure for the frequency dispersion curves is described too tersely; the shifts should be specified explicitly in an equation or table so that the comparison can be reproduced.","section":"Section 3, Figure 2"},{"comment":"'Sytem' is a typo for 'system'.","section":"Section 3, Figure 5"},{"comment":"The overline notation for the Löwdin-orthogonalized basis is not defined in the main text; please define it before Eq. (57).","section":"Section 2, Eq. (57)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for the journal. The self-citations (Refs. 45, 46) are directly relevant and not inappropriate. The main gap is benchmark coverage: if the authors can supply long-range QED-CCSD or equivalent reference data, the central claim would be substantially strengthened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new piece here is real: SC-QED-MP2 is not QED-MP2 with a different orbital guess. It uses the SC-QED-HF reference, which dresses the electrons with cavity photons and captures some electron-photon correlation at mean-field level, and the second-order energy expression with infinite photon sums is new. The size-intensivity analysis is also useful, especially the scaling argument for the QED-MP2 long-range divergence. The derivation in the SI is parameter-free and coherent, and the derivation of the general multi-boson Rayleigh-Schrodinger framework is a nice bonus. Credit is due for that. Self-citation of SC-QED-HF is fine; that method is published and independently tested.\n\nNow the soft spots, in proportion. The stress-test note is right: in the intermolecular plots, the central claim is that SC-QED-MP2 avoids the unphysical divergence of QED-MP2 and LF-MP2, but no QED-CCSD reference is shown in those curves. So the reader cannot see whether the SC-QED-MP2 plateau is the physically correct long-range energy or just a different approximation that happens not to blow up. The theoretical argument for size-intensivity is credible, and I do not think this is a load-bearing flaw, but it is a real gap in the numerical evidence. Relatedly, the benchmark used throughout is QED-CCSD built on QED-HF, a reference the paper itself criticizes as non-size-intensive and origin-dependent for charged systems. The authors justify this by saying they focus on electron-photon correlation; for neutral isolated molecules that is reasonable, but for long-range intermolecular interactions it is exactly the regime where the criticism matters most. So the ranking of methods in the regime that most distinguishes them is not numerically verified.\n\nMinor: the implementation is not released, only available on request, though the data are on Zenodo. That is tolerable for a methods paper but worth flagging.\n\nWho gets value: method developers in ab initio polaritonics, and people choosing a correlated method for strongly coupled molecular cavities. The paper deserves a serious referee; I would send it to review with a request for a long-range reference curve in a small system where a more trustworthy method (or an SC-based coupled cluster, once available) can be computed, and for clarification of the LF-MP2 size-intensivity argument. The derivation is solid, and the central idea is worth engaging with.","headline":"A genuinely new MP2 variant for strongly coupled polaritons with a credible derivation; the long-range benchmark gap is real but does not sink the paper.","tokens_in":23645,"tokens_out":1578,"would_cite":true,"duration_ms":19596,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.15.Md","31.15.xq","42.50.Pq"],"model":"deepseek-v4-flash","headline":"SC-QED-MP2, a perturbation theory built on cavity-consistent orbitals, claims to capture field-induced electron-photon correlation at mean-field level and avoid unphysical long-range behavior seen in QED-MP2 and LF-MP2.","keywords":["polaritonic chemistry","cavity QED","Møller-Plesset perturbation theory","strong coupling","QED Hartree-Fock","molecular orbitals","electron-photon correlation","size-intensivity"],"falsifier":"Take two hydrogen molecules far apart inside a cavity with the polarization along the displacement direction and compute the SC-QED-MP2 dissociation curve: if the curve diverges rather than reaching a plateau, the claimed size-intensivity of the zeroth-order Hamiltonian fails. Alternatively, compute SC-QED-MP2 energies for a charged molecule after translating the origin: any change would contradict the claimed origin invariance.","tokens_in":22523,"feed_emoji":"⚛️","tokens_out":7611,"duration_ms":67983,"temperature":0.7,"pith_summary":"Inside an optical cavity, molecules and quantized fields form polaritons, and the cheapest reliable route to their correlated ground states is Møller-Plesset perturbation theory. This paper argues that previous QED-MP2 schemes fail because their zeroth-order Hamiltonian uses molecular orbitals that are not consistent with the cavity: the QED-HF orbitals are not size-intensive and, for charged systems, not origin-invariant. The authors develop SC-QED-MP2, built on strong-coupling QED Hartree-Fock orbitals that already dress electrons with cavity photons. They show that this method reproduces QED-coupled-cluster reference trends for coupling and frequency dispersions, and it avoids the unphysical long-range intermolecular curves that QED-MP2 and its Lang-Firsov variant produce. The paper thereby states that a fully consistent orbital framework for the zeroth-order Hamiltonian—not just a well-chosen wave function parametrization—is what carries the accuracy.","feed_headline":"SC-QED-MP2 reproduces polariton energies without long-range blowups","feed_subtitle":"Consistent cavity orbitals capture electron–photon correlation at mean-field level, fixing dissociation curves.","key_machinery":"The machinery is the strong-coupling QED Hartree-Fock reference: $|\\psi_{SC}\\rangle = \\exp(-\\lambda/\\sqrt{2\\omega}\\sum_p \\eta_p \\tilde{E}_{pp}(b-b^\\dagger))|HF,0\\rangle$, an orbital-specific coherent-state dressing in the basis that diagonalizes the dipole operator $(d\\cdot\\epsilon)$. The parameters $\\eta_p$ are variationally optimized, and Gaussian factors $Q_{pq}=\\exp(-\\lambda^2/(4\\omega)(\\eta_p-\\eta_q)^2)$ built into the transformed Hamiltonian carry cavity-induced correlation into the mean-field Fock operator, making it origin-invariant and size-intensive. The second-order energy correction sums double electronic excitations with arbitrary photon number $n$, single excitations with $n\\geq 1$, and purely photonic excitations with $n\\geq 2$, with denominators $n\\omega$ plus orbital energy differences.","core_discovery":"The central claim is that SC-QED-MP2 accurately reproduces field-induced electron-photon correlation effects because those effects are already present at the mean-field level, in the strong-coupling QED Hartree-Fock reference. The reference is built in the dipole basis, the basis that diagonalizes the dipole operator, with an orbital-specific coherent-state transformation; the resulting Fock operator is origin-invariant and size-intensive, unlike the QED-HF Fock operator. On top of this reference, the second-order correction captures single, double, and purely photonic excitations across photon numbers. In benchmark comparisons against QED-CCSD, SC-QED-MP2 matches the reference trends for cavity-coupling and frequency dispersions of ammonia and gives physical dissociation curves for hydrogen, water, and benzene-water complexes, while QED-MP2 and LF-MP2 show unphysical long-range behavior when the cavity polarization has a component along the molecular displacement.","pith_inferences":["A natural next test is ultrastrong coupling, beyond lambda around 0.05 a.u.; the paper's infinite-coupling exactness argument predicts SC-QED-MP2 should continue to improve relative to QED-CCSD as lambda grows, while QED(np-HF)-MP2 should degrade, an ordering that is directly measurable.","If the size-intensivity result transfers, SC-QED-MP2 should become the default affordable method for cavity-modified intermolecular interactions, including cases such as the benzene-water metastable complex where QED-MP2 incorrectly turns an unbounded interaction into a bound one.","The paper's emphasis on the dipole basis suggests that multi-mode cavities cannot be handled by simply diagonalizing each mode; an orbital framework that simultaneously treats multiple non-commuting dipole directions will be needed for realistic cavities."],"forward_implications":["SC-QED-MP2 reproduces the QED-CCSD coupling and frequency dispersions for ammonia across the tested range, and becomes the most accurate perturbative method at large coupling because its reference becomes exact in the infinite-coupling limit.","QED-MP2 and LF-MP2 produce unphysical, diverging dissociation curves for two far-apart molecules when the polarization has a component along the displacement direction; SC-QED-MP2 and QED(np-HF)-MP2 remain well behaved, identifying the orbital basis as the source of the failure.","QED(np-HF)-MP2 is well behaved but is expected to lose accuracy at very strong coupling, since its zeroth-order Hamiltonian contains no cavity effects on the orbitals; SC-QED-MP2 improves exactly in that regime.","Because SC-QED-MP2 is size-intensive and based on a mean-field reference that already includes electron-photon correlation, it offers an affordable MP2-level route to strongly coupled polaritonic ground states, and the same reference should support QED versions of CC2, CC3, and active-space methods."],"supporting_citations":[{"why":"Supplies the original QED-MP2 and QED(np-HF)-MP2 formulations that this paper extends and compares against.","marker":"[44]"},{"why":"Introduces the SC-QED-HF theory whose orbitals form the zeroth-order Hamiltonian of SC-QED-MP2.","marker":"[45]"},{"why":"Demonstrates the unphysical properties of QED-HF orbitals and provides the QED-CCSD method used as the benchmark.","marker":"[28]"},{"why":"Presents the Lang-Firsov MP2 scheme, the alternative with a similar wave function parametrization but different basis, whose long-range failure is contrasted with SC-QED-MP2.","marker":"[47]"},{"why":"Provides the second-order convergence algorithm for SC-QED-HF that makes the MP2 implementation practical.","marker":"[46]"},{"why":"Supplies the eT program in which the new perturbative methods are implemented.","marker":"[58]"}],"fun_headline_variants":["Strong-coupling QED-MP2 tames long-range polariton blowups","Orbital choice fixes cavity QED perturbation theory","Cavity orbitals make QED-MP2 physical at strong coupling","SC-QED-MP2: mean-field captures electron-photon correlation","Dipole-basis orbitals cure QED perturbation theory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical ranking of SC-QED-MP2 against its competitors assumes that QED-CCSD built on QED-HF is an accurate reference for strongly coupled ground states, even though the paper itself argues that QED-HF orbitals are ill-defined, non-size-intensive, and origin-dependent for charged systems.","fun_headline_variants_meta":{"raw":{"variants":["Strong-coupling QED-MP2 tames long-range polariton blowups","Orbital choice fixes cavity QED perturbation theory","Cavity orbitals make QED-MP2 physical at strong coupling","SC-QED-MP2: mean-field captures electron-photon correlation","Dipole-basis orbitals cure QED perturbation theory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000844,"raw_usage":{"total_tokens":3634,"prompt_tokens":864,"completion_tokens":2770,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":480,"completion_tokens_details":{"reasoning_tokens":2680}},"tokens_in":480,"tokens_out":2770,"duration_ms":19941,"temperature":1.0,"reasoning_tokens":2680,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:29:30.162673+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take two hydrogen molecules far apart inside a cavity with the polarization along the displacement direction and compute the SC-QED-MP2 dissociation curve: if the curve diverges rather than reaching a plateau, the claimed size-intensivity of the zeroth-order Hamiltonian fails. Alternatively, compute SC-QED-MP2 energies for a charged molecule after translating the origin: any change would contradict the claimed origin invariance.","supporting_citations":[{"cited_title":"Perturbation theoretical approaches to strong light--matter coupling in ground and excited electronic states for the description of molecular polaritons","cited_arxiv_id":null,"evidence_quote":"Supplies the original QED-MP2 and QED(np-HF)-MP2 formulations that this paper extends and compares against."},{"cited_title":"R.; Haugland, T","cited_arxiv_id":null,"evidence_quote":"Introduces the SC-QED-HF theory whose orbitals form the zeroth-order Hamiltonian of SC-QED-MP2."},{"cited_title":"S.; Ronca, E.; Kj nstad, E","cited_arxiv_id":null,"evidence_quote":"Demonstrates the unphysical properties of QED-HF orbitals and provides the QED-CCSD method used as the benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the Lang-Firsov MP2 scheme, the alternative with a similar wave function parametrization but different basis, whose long-range failure is contrasted with SC-QED-MP2."},{"cited_title":"R.; Castagnola, M.; Koch, H","cited_arxiv_id":null,"evidence_quote":"Provides the second-order convergence algorithm for SC-QED-HF that makes the MP2 implementation practical."},{"cited_title":"D.; Kj nstad, E","cited_arxiv_id":null,"evidence_quote":"Supplies the eT program in which the new perturbative methods are implemented."}],"review_version":1}