{"id":"f06170dc-3ea6-461c-aa97-562b1f3d41e1","arxiv_id":"2506.00217","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Static QED-HF/QED-DFT response theory shows the first hyperpolarizability of p-nitroaniline drops by more than 20% at lambda=0.05 a.u. in a z-polarized cavity, while the polarizability changes by only a few percent.","lead":"The authors implemented static electric dipole polarizability and first hyperpolarizability tensors for QED-Hartree-Fock and QED-DFT, and applied them to p-nitroaniline in a single-mode optical cavity. They find the hyperpolarizability can change by over 20% under strong coupling while the polarizability changes by only a few percent, suggesting cavity control of nonlinear optical properties.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fixed-geometry and fixed-orientation assumption underlies the -22% beta_bar change; the cavity potential favors the y-polarized orientation with only ~2% change, so the headline number may not survive cavity-induced relaxation or orientational averaging.","rationale":"The paper's central contribution is a response-theory implementation for QED-HF and QED-DFT, with numerical demonstration that the QED-HF first hyperpolarizability of p-nitroaniline can change by more than 20% under a z-polarized cavity at strong coupling. The strongest condition for this claim to hold as a physical statement is that the computed quantity represents a molecule that can actually adopt the assumed orientation and geometry: the molecule is held rigid at its isolated-molecule optimized structure, and the orientation is fixed rather than determined by the cavity potential. The paper itself flags this limitation in Sec. IV, noting that the lowest-energy orientation is y-polarized and that orientational averaging or cavity-induced rotations/distortions are needed for meaningful predictions. This is load-bearing because the cavity potential is orientation-dependent through the DSE term, and at the large coupling lambda = 0.05 a.u. the energy differences between orientations could be significant. The reader's weakest_assumption identifies exactly the same issue, and the paper's own text supports that reading. The formal QED-HF response equations appear internally consistent and are checked against finite-field calculations, and the basis-set convergence is documented, so I do not see a more serious internal flaw. The appropriate verdict remains conditional: the 20% hyperpolarizability change should be presented as a mean-field, rigid-geometry, fixed-orientation estimate until geometry relaxation and orientational averaging are assessed.","tokens_in":19396,"tokens_out":8275,"duration_ms":95354,"concrete_test":"Re-optimize the p-nitroaniline geometry at the QED-HF level with lambda = 0.05 a.u. and omega_cav = 0.1 Eh for each cavity polarization direction (x, y, z), using analytic or numerical gradients of the QED-HF energy. Then recompute beta_bar at the relaxed z-polarized geometry and compare it with the rigid-geometry value. In addition, compute a Boltzmann-weighted orientational average at room temperature using the QED-HF ground-state energies for the three orientations. If the relaxed and/or orientationally averaged change in beta_bar remains at least 20%, the fixed-geometry concern does not alter the headline; if the change drops below roughly 5-10%, the central claim should be reframed as conditional on rigid, artificially oriented molecules.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline numerical claim is the -22% change in the isotropically averaged first hyperpolarizability of p-nitroaniline at lambda = 0.05 a.u. with a z-polarized cavity mode (Sec. IV). This number is computed at the isolated-molecule optimized geometry and at a fixed molecular orientation, even though the QED-HF energy depends on molecular orientation through the dipole self-energy term. The authors explicitly note in Sec. IV that the lowest-energy orientation is y-polarized, which is the orientation giving the smallest change to beta_bar (~2%), and that meaningful predictions require orientational averaging or accounting for cavity-induced rotations and geometric distortions. Thus the reported 20% decrease is a rigid-frame, fixed-orientation result rather than a prediction for a molecule free to relax in the cavity. Because lambda = 0.05 a.u. corresponds to a picocavity-scale perturbation and the DSE term is quadratic in the dipole moment, cavity-induced changes in bond lengths, angles, and molecular orientation could substantially alter the beta tensor components entering Eq. 30, changing the magnitude or even the sign of the reported effect. The internal response-theory implementation is verified by finite-field calculations, so the concern is not about formal correctness of the QED-HF equations but about whether the computed quantity is the physically relevant one for a molecule in a cavity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives and implements static electric dipole response theory for mean-field cavity QED: QED-HF and QED-DFT for the polarizability, and QED-HF for the first hyperpolarizability. The implementation is verified against finite-field calculations and, in the non-QED limit, against Gaussian. The authors apply these methods to p-nitroaniline in a single-mode cavity and find that the isotropic polarizability changes by only a few percent up to |lambda| = 0.05 a.u., whereas the isotropic first hyperpolarizability can change by more than 20% when the cavity mode is polarized along the molecular z-axis. They also present basis-set convergence studies showing that diffuse functions are essential for qualitatively correct cavity-induced changes. The central formal result is a response-theory framework that extends standard HF/DFT response equations with photon contributions via the coherent-state representation.","tokens_in":19629,"tokens_out":5700,"duration_ms":59997,"significance":"If the numerical conclusions hold, this work provides the first ab initio response-theory treatment of static first hyperpolarizabilities in cavity QED environments, a capability that is timely given experimental interest in nonlinear optics under strong coupling. The derivation is self-contained and the implementation is carefully verified against finite-field references, which are explicit strengths. The paper also gives a useful and honest discussion of orientation dependence, including the observation that the energy is lowest for the orientation with the smallest hyperpolarizability change. However, the headline quantitative result is computed at a fixed geometry and fixed orientation, and the authors themselves state that orientational averaging or geometry relaxation is needed for meaningful predictions. This limitation, while acknowledged in Sec. IV, is not reflected in the abstract, and it is load-bearing for the physical interpretation of the reported 20% effect. The methodological contribution is sound, but the presentation of the numerical claim needs to be brought in line with the stated limitations.","major_comments":[{"comment":"The headline -22% change in the isotropically averaged first hyperpolarizability is computed at the isolated-molecule optimized geometry and at a fixed orientation with the cavity polarized along the z-axis. In Sec. IV the authors themselves note that the lowest-energy orientation is y-polarized, which gives only about a 2% change, and that 'meaningful predictions of cavity modified properties should either involve appropriate orientational averaging or account for cavity-induced rotations and geometric distortions.' Because the dipole self-energy term is quadratic in the dipole moment, cavity-induced reorientation or geometric distortion could alter the magnitude and even the sign of the hyperpolarizability modification. The abstract currently presents the -22% value without this caveat, so the manuscript should either provide relaxed or orientationally averaged results, or explicitly reframe the headline as a rigid-geometry, fixed-orientation model result that illustrates the potential sensitivity rather than a prediction for a molecule free to relax in the cavity.","section":"Abstract and Sec. IV (Fig. 2)"},{"comment":"In the first-order QED-HF response equation, the sum over jb is contracted with kappa^alpha_{e,ia}, but the response vector should carry the jb index: the term should read sum_{jb} (A_{ia,jb}+B_{ia,jb}) kappa^alpha_{e,jb}. As written, the equation is internally inconsistent because the summation index jb appears only in the A and B matrices and not in the response vector. Please correct this index and verify that the accompanying text and code documentation use the same convention, since this is a central equation of the formalism.","section":"Eq. (22)"}],"minor_comments":[{"comment":"The sentence 'we also use the first-order response function to evaluate the the first hyperpolarizability tensor' contains a duplicated article 'the'.","section":"Sec. I"},{"comment":"The phrase 'essentially converged at at any ζ-level' contains a duplicated 'at'; please fix this typo.","section":"Sec. IV"},{"comment":"The phrase 'oscillator strenghts' should read 'oscillator strengths'.","section":"Supporting Information"},{"comment":"The derivation would benefit from a brief statement that the photonic response contribution to the hyperpolarizability is entirely encoded in the kappa^alpha_p term within h^alpha_{ia,jb}, since a reader might otherwise expect explicit photon-number response terms in Eq. (25). This is a clarity suggestion, not a correctness concern.","section":"Sec. II B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid methodological contribution with careful numerical verification, and the authors are transparent about the fixed-geometry limitation in Sec. IV. The main issue for publication is that the abstract and conclusions present the -22% beta change as a headline result without prominently carrying the caveat that it is a rigid-orientation result; given the authors' own statement that the favored cavity orientation gives only ~2% change, this needs either additional calculations (e.g., cavity-optimized geometry or orientational averaging) or a clear reframing. The index typo in Eq. (22) is local but should be corrected before acceptance. I do not see grounds for rejection; the formal derivation and implementation appear sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First the punchline: the paper is a legitimate first. It implements static α and β tensors for QED-HF and QED-DFT response theory, verifies the code against finite-field and Gaussian in the non-QED limit, and reports a cavity-induced drop of more than 20% in the isotropic first hyperpolarizability of p-nitroaniline. That number, however, is a fixed-orientation, fixed-geometry mean-field estimate, not a free-molecule prediction. The authors know this and say it in Sec. IV, but the abstract still leads with the -22%.\n\nWhat is genuinely new: the QED-HF hyperpolarizability is the first ab initio response-theory treatment of a second-order static property in a cavity, and the derivation cleanly shows how the photon response enters through the first-order κ vector without any fitted parameters. The finite-field verification is real evidence that the implementation matches the equations. The basis-set study is also a useful service: it shows that non-augmented basis sets can reverse the sign of the cavity-induced α change, which is a practical warning for the field.\n\nSoft spots, in rough order of importance. First, the headline -22% depends on the z-polarized orientation and the isolated-molecule geometry. The DSE term is orientation-dependent, and the authors themselves note that the lowest-energy orientation is y-polarized, where the β change is only about 2%. Anyone citing this as 'a cavity can change β by 20%' will over-read the abstract; the paper should make the fixed-orientation caveat visible in the abstract or in the first paragraph of Sec. IV. Second, Eq. (22) has an index typo: the sum over jb multiplies κ^α_{e,ia} instead of κ^α_{e,jb}, so the contraction with (A+B) looks wrong on paper. The finite-field agreement says the code is right, but the typo should be fixed. Third, the unrestricted HF/DFT results come with no spin-contamination check, and β is computed only at QED-HF, which makes the title's 'mean-field approaches' slightly broader than the content for second-order properties. None of this undermines the method; it defines its scope.\n\nThe paper is honest about mean-field QED exaggerating cavity effects and cites the relevant concurrent work. I would send it to peer review. Revision should fix the typo, add or address the spin-contamination point, and rebalance the abstract so the -22% is explicitly labeled as an oriented, rigid-molecule mean-field estimate.","headline":"A genuinely new QED-HF/QED-DFT static response implementation, verified numerically, with a headline -22% β change that is a fixed-orientation mean-field result the authors themselves hedge in Sec. IV.","tokens_in":20164,"tokens_out":3794,"would_cite":true,"duration_ms":37482,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["33.15.Kr","42.50.Pq"],"model":"deepseek-v4-flash","headline":"The paper derives static response functions for mean-field cavity QED and finds, for p-nitroaniline, that a cavity mode polarized along the molecular axis cuts the isotropically averaged first hyperpolarizability by more than 20% while…","keywords":["cavity quantum electrodynamics","QED-HF","QED-DFT","electric dipole polarizability","first hyperpolarizability","response theory","p-nitroaniline","strong light-matter coupling"],"falsifier":"Recompute the QED-HF isotropically averaged first hyperpolarizability of $p$-nitroaniline at $|\\lambda| = 0.05$ a.u. with the molecular geometry re-optimized inside the cavity (or with an explicit orientational average); if the z-polarized cavity no longer gives a drop of more than 20% in $\\bar{\\beta}$, the fixed-geometry assumption is the load-bearing element of the prediction.","tokens_in":19182,"feed_emoji":"💡","tokens_out":20786,"duration_ms":190871,"temperature":0.7,"pith_summary":"This paper develops and implements static response theory for mean-field cavity quantum electrodynamics, targeting the electric dipole polarizability and first hyperpolarizability tensors of a molecule coupled to a single optical cavity mode. The central finding is an asymmetry: at a realistic single-molecule coupling strength ($|\\lambda| = 0.05$ a.u.), the isotropically averaged polarizability of $p$-nitroaniline changes by only a few percent, but the QED-HF first hyperpolarizability can be suppressed by more than 20% when the cavity polarization aligns with the molecular principal axis, while other polarizations slightly enhance it. The result matters because it identifies the first hyperpolarizability, not the polarizability, as the sensitive probe of strong light–matter coupling, and because it demonstrates, for the first time in an ab initio response framework, that a cavity can significantly renormalize higher-than-linear molecular response properties. The authors also establish practical computational requirements, showing that diffuse basis functions are essential for reliable cavity-induced changes.","feed_headline":"Cavity coupling cuts a molecule's nonlinear response by >20%","feed_subtitle":"Hyperpolarizability falls by >20% when the cavity is aligned with the molecular axis; polarizability shifts about 2-5%.","key_machinery":"The central object is the coherent-state representation of the Pauli–Fierz Hamiltonian, in which the photon displacement is absorbed into a unitary transformation so that the QED-HF and QED-DFT ground states reduce to standard self-consistent-field problems with dipole-self-energy-modified one- and two-electron integrals. The response machinery couples the usual RPA-like electronic Hessian $A+B$, augmented by dipole self-energy contributions, to a single photonic response amplitude $\\kappa_p$ via Eqs. (22)–(23); the static polarizability retains the standard HF form $\\alpha_{e,\\alpha\\beta} = 2\\sum_{ia}\\mu_{\\alpha,ia}\\kappa^{\\beta}_{e,ia}$, but the electronic response vector is determined together with the photon response. The QED-HF first hyperpolarizability $\\beta_{e,\\alpha\\beta\\gamma}$ is assembled from products of first-order response vectors and third energy derivatives, generalizing the standard HF expression. This structure lets the cavity renormalize $\\beta$ through the DSE-modified orbital Hessian and the photonic amplitude, which is how orientation-specific polarizations can produce a >20% suppression of $\\bar{\\beta}$ while leaving $\\bar{\\alpha}$ nearly intact.","core_discovery":"The authors work with the Pauli–Fierz Hamiltonian in the dipole approximation and length gauge, then apply the coherent-state transformation to define QED-HF and QED-DFT ground states. First-order response theory with respect to a static external electric field yields the electric dipole polarizability tensor $\\alpha_{e,\\alpha\\beta}$ in the same algebraic form as standard HF theory, but with the electronic response vector coupled to a photonic amplitude, and it yields the first hyperpolarizability tensor $\\beta_{e,\\alpha\\beta\\gamma}$ (the second-order response coefficient governing nonlinear optical effects such as second-harmonic generation) for QED-HF. For $p$-nitroaniline at $|\\lambda| = 0.05$ a.u. and $\\omega_{\\mathrm{cav}} = 0.1\\,E_h$, the isotropically averaged static polarizability $\\bar{\\alpha}$ changes by only $\\approx 2$ – $5\\%$, whereas the isotropically averaged QED-HF hyperpolarizability $\\bar{\\beta}$ decreases by more than $20\\%$ (from $133.5$ to $104.4$ a.u.) when the cavity mode is polarized along the molecular principal axis; x- and y-polarized modes instead increase $\\bar{\\beta}$ by as much as $7\\%$ and $2\\%$, respectively. The static response properties are independent of the cavity frequency, a direct consequence of the coherent-state formulation, and the cavity-induced changes in $\\bar{\\beta}$ are more sensitive to basis-set quality than those in $\\bar{\\alpha}$, with non-augmented basis sets able to yield qualitatively wrong signs for the polarizability shift.","pith_inferences":["A cautious extension: because mean-field QED approaches tend to exaggerate cavity-induced electronic structure changes (as the paper itself notes), the >20% hyperpolarizability suppression may represent an upper bound, with correlated treatments such as QED coupled cluster likely giving a smaller effect.","If cavity-induced geometry relaxation is allowed, the energy of $p$-nitroaniline is lowest when the cavity is y-polarized, exactly the orientation giving the smallest $\\bar{\\beta}$ change; freely rotating molecules should therefore show much smaller orientation-averaged modifications than the z-polarized rigid-geometry value.","The strong orientation dependence suggests a testable design rule: anchoring a molecule's orientation relative to the cavity polarization could toggle second-harmonic generation on and off, since rotating the polarization from x to z changes $\\bar{\\beta}$ from +7% to −22%.","The same first-order response formalism should extend to the second hyperpolarizability (third-order response) and to chiral or circularly polarized cavity modes, where orientation-dependent effects on $\\beta$ and higher tensors are likely to be even richer."],"forward_implications":["Cavity polarization direction becomes a control variable for second-order nonlinear response: for $p$-nitroaniline, rotating the cavity mode from x-polarized to z-polarized changes $\\bar{\\beta}$ from about a 7% enhancement to a >20% suppression.","Static polarizability is nearly insensitive to the cavity at these coupling strengths, so mean-field cavity-QED descriptions of electrostatically dominated properties need not be revised for first-order response.","Because static response properties are independent of $\\omega_{\\mathrm{cav}}$ in this coherent-state mean-field framework, any frequency dependence in polarizabilities or hyperpolarizabilities must enter through dynamic response theory or through beyond-mean-field correlations.","Diffuse basis functions are required for quantitatively reliable cavity-induced changes in both $\\bar{\\alpha}$ and $\\bar{\\beta}$; small non-augmented basis sets can give the wrong sign for the polarizability change."],"supporting_citations":[{"why":"supplies the coherent-state transformation operator used to define the QED-HF and QED-DFT ground states.","marker":"[66]"},{"why":"gives the prior polaritonic linear-response framework that this work extends to static second-order (hyperpolarizability) response.","marker":"[84]"},{"why":"supplies the dipole-self-energy-augmented RPA matrices that appear in the QED response equations.","marker":"[52]"},{"why":"is the source of the caveat that meaningful cavity-modified predictions require orientational averaging or cavity-induced geometry relaxation.","marker":"[55]"},{"why":"provides the standard Hartree–Fock hyperpolarizability expression that the QED-HF formula generalizes.","marker":"[117]"},{"why":"defines picocavity mode volumes, grounding the largest coupling strength considered as a realistic single-molecule value.","marker":"[134]"}],"fun_headline_variants":["Cavity coupling slashes hyperpolarizability by >20%","QED cavity: polarizability barely moves, hyperpolarizability dives 20%","Hyperpolarizability exposed to cavity effects; polarizability robust","Cavity-coupled QED cuts p-nitroaniline hyperpolarizability by >20%","Cavity mode polarization can swing hyperpolarizability by 20%+"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations use the isolated-molecule optimized geometry of $p$-nitroaniline for every cavity condition, so the reported changes assume the molecule neither reorients nor distorts inside the cavity; if it does either, the magnitude or even the sign of the hyperpolarizability shift could differ.","fun_headline_variants_meta":{"raw":{"variants":["Cavity coupling slashes hyperpolarizability by >20%","QED cavity: polarizability barely moves, hyperpolarizability dives 20%","Hyperpolarizability exposed to cavity effects; polarizability robust","Cavity-coupled QED cuts p-nitroaniline hyperpolarizability by >20%","Cavity mode polarization can swing hyperpolarizability by 20%+"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001118,"raw_usage":{"total_tokens":4750,"prompt_tokens":1141,"completion_tokens":3609,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":757,"completion_tokens_details":{"reasoning_tokens":3508}},"tokens_in":757,"tokens_out":3609,"duration_ms":30573,"temperature":1.0,"reasoning_tokens":3508,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:10:14.739100+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the QED-HF isotropically averaged first hyperpolarizability of $p$-nitroaniline at $|\\lambda| = 0.05$ a.u. with the molecular geometry re-optimized inside the cavity (or with an explicit orientational average); if the z-polarized cavity no longer gives a drop of more than 20% in $\\bar{\\beta}$, the fixed-geometry assumption is the load-bearing element of the prediction.","supporting_citations":[{"cited_title":"Castagnola , author R","cited_arxiv_id":null,"evidence_quote":"gives the prior polaritonic linear-response framework that this work extends to static second-order (hyperpolarizability) response."},{"cited_title":"Vu , author G","cited_arxiv_id":null,"evidence_quote":"supplies the dipole-self-energy-augmented RPA matrices that appear in the QED response equations."},{"cited_title":"Norman \\ and\\ author K","cited_arxiv_id":null,"evidence_quote":"provides the standard Hartree–Fock hyperpolarizability expression that the QED-HF formula generalizes."}],"review_version":1}