{"id":"929a7791-5215-4921-a179-87344af5efce","arxiv_id":"2507.19180","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An unpolarized Fabry-Perot cavity is modeled in QED-CC by two degenerate perpendicular modes, and point-group symmetry is exploited to assign and compute polaritonic excited states.","lead":"A new way to compute how molecules and light interact in an unpolarized optical cavity treats two perpendicular light directions at once. It uses molecular symmetry to predict and classify the mixed light-matter states, with test calculations on benzene, fluorobenzene, azulene, and H2.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The two-mode equal-coupling truncation is the load-bearing idealization: Hamiltonian-level C∞ symmetry is proven, but the step from that model to a real unpolarized Fabry-Pérot cavity is unquantified and untested against mode-continuum effects.","rationale":"The central claim is that two degenerate perpendicular modes with equal coupling represent an unpolarized Fabry-Pérot cavity and that point-group symmetry then determines which polaritonic states mix and can be computed selectively. What must be true for that claim to hold is, first, that the two-mode model is the right physical description of the cavity and, second, that the CC/EOM implementation actually realizes the symmetry in its computed states. The paper proves the second point at the Hamiltonian level and, encouragingly, the H2 k∥ tables show exact B2u/B3u degeneracy, which is consistent with the predicted C∞ symmetry for that orientation. The first point, however, is assumed rather than demonstrated: the physical cavity is replaced by a single frequency with two polarizations, and the contributions of other modes, mode continua, and mass renormalization are acknowledged but not quantified. This is not an internal inconsistency, and it is a standard single-mode idealization, but it is load-bearing because every symmetry-derived selection rule and numerical prediction in Sections 3.1-3.5 is a statement about that truncated model. The reader's weakest assumption identifies exactly this issue, and the recommended conditional verdict remains appropriate: the methodology and symmetry analysis are credible, but the domain of validity of the two-mode approximation should be established before the unpolarized-cavity results are used as quantitative predictions. A mode-convergence test is the most direct way to settle whether the concern lands.","tokens_in":47026,"tokens_out":26937,"duration_ms":290443,"concrete_test":"Perform a mode-convergence test: repeat the H2 k∥ spectrum and the benzene polarized-versus-unpolarized density-shift calculations with a growing hierarchy of cavity modes, starting with the present two-mode pair and adding 4, 8, and 16 modes built from degenerate perpendicular pairs at the same frequency, then pairs at nearby frequencies within the light-matter coupling window, including the full dipole self-energy sum rather than a single pair. If the B2u/B3u degeneracy, the point-group state labels, or the benzene Δρ change by more than about 10% as the mode set grows, the two-mode truncation is the decisive approximation and the central claim must be restated as applying only to the two-mode model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The algebraic core of the paper is sound: for a fixed pair of degenerate perpendicular modes, Eqs. (9)-(18) correctly show that the two-mode model Hamiltonian is C∞-invariant, and the direct-product point-group analysis follows from that. The load-bearing step is the modeling assumption that precedes the symmetry proof: that a real unpolarized Fabry-Pérot cavity is faithfully represented by one frequency with exactly two orthonormal polarization vectors and equal coupling strength λα = λ̄α. This enters through Eqs. (8), (14)-(15), and Section 2.4. A planar cavity supports a continuum of transverse modes for each longitudinal frequency, and the dipole self-energy contains contributions from all of them; the associated mass renormalization is explicitly deferred to future work in Sections 2.1 and 4. If additional near-resonant modes or the continuum renormalize the levels or induce couplings between symmetry blocks that the two-mode group analysis keeps separate, then the computed state assignments, degeneracies, and density shifts (including the reported approximately doubled benzene density shift) are properties of the truncated model rather than of the physical unpolarized cavity. The eT cross-check in Section 2.7 validates only the single-polarization special case, not the two-mode unpolarized construction. The central claim is therefore established for an idealized two-mode model, but its transfer to experiment rests on an unquantified convergence assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends polaritonic coupled-cluster theory (QED-CC/EOM-QED-CC) to unpolarized optical cavities by explicitly including two degenerate cavity modes with perpendicular linear polarizations and equal coupling. Sections 2.2 and 2.3 derive the rotational and reflection symmetries of the two-mode Hamiltonian at the operator level and argue that the coupled light-matter system is form-invariant under rotations of the polarization basis. Section 2.6 adapts direct-product point-group decomposition to the photonic indices, giving selection rules and block-structured tensors for the coupled-cluster implementation. Numerical results obtained with the Qcumbre/CFOUR implementation cover ground-state densities of benzene, fluorobenzene, and azulene in linearly polarized versus unpolarized cavities, and excited-state energy scans of H2 in both cavity types, including assignments of polaritonic states by irreducible representation and a discussion of truncation-induced crossings. The paper concludes that unpolarized cavities can be viewed as effectively combining two orthogonal linear polarizations and that point-group symmetry enables targeted excited-state calculations.","tokens_in":47366,"tokens_out":18368,"duration_ms":183312,"significance":"If the model assumption is accepted, the paper delivers a practically useful and internally consistent symmetry framework: the two-mode Hamiltonian in Eqs. (8)-(18) is explicitly rotationally invariant, the direct-product decomposition is a parameter-free algebraic argument that yields concrete selection rules in Eqs. (41)-(44), and the H2 state assignments follow the stated symmetry rules. The comparison with an independent QEDFT calculation for azulene is a welcome external check. Strengths include the explicit operator-level derivation, the use of the eT benchmark for the single-polarization limit, and the substantial numerical tables. The main caveats are that the mapping from a real Fabry-Pérot cavity to the two-mode equal-coupling model is not quantified, and the new two-mode numerical implementation is not benchmarked against any independent method; these limit the strength of the physical conclusions but do not undermine the internal consistency of the symmetry derivation itself.","major_comments":[{"comment":"As written, the operator in Eq. (11) does not implement the reflection/C2 transformation whose action is stated in Eq. (12). For n = (1,0) (n1 = 1, n2 = 0), Eq. (11) reduces to exp(iπ N̄/2), which maps ᾱ to -i ᾱ, not to -ᾱ as required by Eq. (12); it also does not square to the identity on odd photon-number states. The correct generator for a reflection across the plane containing n is iπ times the projector onto the perpendicular direction, i.e., iπ[n2^2 α†α + n1^2 ᾱ†ᾱ - n1n2(ᾱ†α + α†ᾱ)], without the factors 1/2 on the diagonal terms. Please correct Eq. (11) and confirm that it satisfies σv^2 = 1, since this is part of the central symmetry derivation.","section":"§2.2, Eq. (11)"},{"comment":"The central claim that a real unpolarized Fabry-Pérot cavity is faithfully represented by two discrete modes with identical frequency and identical coupling strength is asserted rather than justified. A planar cavity supports a continuum of transverse modes for each longitudinal frequency, and the dipole self-energy in Eq. (6) receives contributions from all of them, with mass renormalization growing as more modes are included (as the authors note in Section 2.1 and defer to future work in Section 4). All symmetry statements and numerical results, including the H2 level structure and the density shifts for benzene and azulene, are therefore properties of the truncated two-mode Hamiltonian. Please provide a quantitative estimate of the neglected-mode and continuum contribution, for example a convergence study in the number of explicitly included modes for a small system or an analytic estimate of the effect on Δρ and polariton splittings, or explicitly restrict the physical claims to the two-mode model.","section":"§2.1 and §2.4; Eqs. (8), (15), (20)"},{"comment":"The implementation is verified only in the single-polarization limit against the eT program. The new two-mode unpolarized results, including density differences, degeneracies, and avoided crossings, are not checked against any independent implementation. Given that this is the paper's main methodological novelty, please add at least one independent validation for the two-mode case, for example (i) a QED-FCI or QED-CI comparison for H2 with the same two-mode Hamiltonian, or (ii) a numerical test that energies and Δρ are invariant under an arbitrary polarization-basis rotation angle θ in Eq. (9). Such a check would confirm that the reported state assignments and crossings are not artifacts of the implementation.","section":"§2.7 and Results, Figs. 3-12"}],"minor_comments":[{"comment":"The coupling strength is stated as λ = 0.1 in the caption of Fig. 3, whereas Section 3 and Table 4 use λ = 0.05; please harmonize these values and state explicitly which value was used for each figure.","section":"§3.1, Fig. 3 vs. Table 4"},{"comment":"The caption states a cavity frequency of 12.68 eV, while the text and Tables 5-12 use 12.48 eV; please correct the discrepancy.","section":"Fig. 12 caption"},{"comment":"The last two columns are both labelled ϵ∥(B1u); the second one should likely be ϵ∥(B3u) to match the other tables and the symmetry discussion.","section":"Table 5 header"},{"comment":"Please clarify that Cθ in Eq. (9) is a photon-basis rotation and that Eq. (18) is a form-invariance statement: the Hamiltonian is invariant under simultaneous rotation of the polarization vectors, and for the full symmetry operation also of the electronic coordinates. As written, a photon-only Cθ acting on D·d gives D'·d with rotated polarization vectors, not literally D·d.","section":"§2.3, Eq. (18)"},{"comment":"The comparison with Flick et al. (Ref. 15) is qualitative because the molecular geometry is said to differ; please state the geometry difference explicitly or provide a quantitative metric (for example the integrated density difference) so that the comparison is more than visual.","section":"§3.3, azulene comparison"},{"comment":"The symmetrization notation around Eqs. (25) and (26) is confusing; please define once whether γαβ, γ̄ᾱβ̄, and γ̄αβ are symmetric in their photonic indices before presenting the unrestricted-sum formula.","section":"§2.5, Eqs. (25)-(26)"},{"comment":"The integral in Eq. (46) is written with both limits as ∞; it should be ∫ from -∞ to ∞ over d3r.","section":"§3, Eq. (46)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and presents a useful extension of polaritonic coupled-cluster theory. The identified issues, including the incorrect reflection operator in Eq. (11) and the absence of an independent two-mode validation, are addressable; I do not recommend rejection. Please require the authors to re-check all symmetry operator formulas and to add the two-mode benchmark or an explicit model re-scoping before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuine, modest advance in QED-CC methodology, not a breakthrough. The authors extend single-polarization QED-CC to two degenerate perpendicular modes representing an unpolarized Fabry-Perot cavity, and they show how direct-product point-group symmetry applies to the photonic indices. The C∞ invariance argument for the two-mode Hamiltonian is correct, and the H2 avoided-crossing assignments follow the stated selection rules. I checked the key equations; they hold together. Credit where due: they also compare against an independent QEDFT calculation for azulene, and they explicitly flag mass renormalization and multi-mode effects as future work. That is honest.\n\nThe soft spots are mostly about verification and presentation. The two-mode implementation is only cross-checked in the single-polarization limit against eT. There is no independent FCI or perturbative benchmark for the genuine two-mode couplings. The claim that the unpolarized cavity produces roughly twice the density shift is plausible from the self-energy structure but is presented without a derivation; it needs either a short proof or a qualifying statement. The tables contain several typos (duplicate column labels in Table 5, stray 'e' formatting in Table 11) that should be cleaned. These are minor but they make the numerics harder to trust than they should be.\n\nThe bigger conceptual point: the entire calculation is for a two-mode model, and the step to a real cavity with a continuum of transverse modes is not quantified. That is exactly the limitation the authors acknowledge, so I do not count it as a hidden flaw. But it does mean the results are properties of the truncated model, not of the physical unpolarized cavity. Readers should not over-interpret the benzene/azulene densities as experimental predictions.\n\nWho benefits: people working on ab initio polaritonic methods, especially QED-CC and EOM-CC variants. The symmetry exploitation is genuinely useful for targeting states. This paper deserves a serious referee: it is clearly within the scope of a methods journal, the algebra is checkable, and the limitations are stated. My recommendation is acceptance after minor-to-moderate revisions: add an independent two-mode benchmark, fix the tables, and either derive or soften the 'twice as high' density claim. I would not desk-reject this, but I would also not accept it in its current form without those changes.","headline":"Real, modest extension of QED-CC to two-mode unpolarized cavities with a sound symmetry argument; lacks an independent two-mode benchmark but deserves peer review.","tokens_in":47821,"tokens_out":3898,"would_cite":true,"duration_ms":40229,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.50.Pq","31.15.-p"],"model":"deepseek-v4-flash","headline":"Two perpendicular, equally coupled cavity modes preserve the rotational symmetry of an unpolarized Fabry-Pérot cavity, which lets point-group symmetry target individual polaritonic states.","keywords":["quantum electrodynamics coupled cluster","unpolarized Fabry-Pérot cavity","point-group symmetry","direct product decomposition","polaritonic excited states","equation-of-motion coupled cluster","avoided crossings","dipole self-energy"],"falsifier":"Compute the two-mode QED-CC ground-state energy and first polaritonic gap for a small molecule such as H$_2$, then repeat with a larger set of modes at the same frequency and coupling; if the two-mode results shift by more than chemical accuracy, the discrete two-mode representation of the unpolarized cavity fails. A complementary test is to measure the coupling strengths for the two orthogonal polarizations of a real Fabry-Pérot cavity and check whether they are equal, since unequal $\\lambda_\\epsilon$ and $\\lambda_{\\bar\\epsilon}$ would lift the predicted rotational degeneracy of molecular orientations.","tokens_in":46833,"feed_emoji":"⚛️","tokens_out":9086,"duration_ms":83272,"temperature":0.7,"pith_summary":"This paper extends quantum-electrodynamic coupled cluster theory to molecules inside an unpolarized optical Fabry-Pérot cavity, a common experimental setting. The central move is to describe the cavity field by two modes with perpendicular polarizations and identical frequency and coupling strength, which keeps the Hamiltonian invariant under rotations about the cavity axis. The authors then adapt point-group symmetry, via direct-product decomposition, to the polaritonic coupled cluster equations so that excited states carry irreducible representation labels and can be computed individually. They apply the method to benzene, fluorobenzene, azulene and H$_2$, finding dense avoided-crossing landscapes and polarization-dependent electron density shifts, including removal of density from a carbon-fluorine bond.","feed_headline":"Two perpendicular modes bring polariton theory to unpolarized cavities","feed_subtitle":"Equal-strength perpendicular photon modes preserve rotational symmetry and let symmetry labels target polaritonic states.","key_machinery":"The central object is the two-mode rotation operator $\\hat C_\\theta = \\exp[\\theta(\\hat\\alpha^\\dagger\\bar{\\hat\\alpha} - \\bar{\\hat\\alpha}^\\dagger\\hat\\alpha)]$ for the perpendicular mode pair, together with the equal-coupling condition $\\lambda_\\alpha=\\bar\\lambda_\\alpha$. This operator generates the $\\mathrm{C}_\\infty$ rotations of the displacement field, and combined with the equal-coupling condition it makes the bilinear term $\\hat D\\cdot\\hat d$ (and hence the full Hamiltonian) rotationally invariant. The second piece of machinery is the direct-product decomposition of the point group applied to the QED-CC equations, which blocks the bilinear coupling integrals and cluster amplitudes by irreducible representation so that the non-linear coupled cluster equations factor into symmetry sectors and the $\\hat R$ excitation operator of EOM-CC can be built for a chosen irrep.","core_discovery":"The paper argues that an unpolarized Fabry-Pérot cavity is faithfully represented, in the dipole approximation, by explicitly including two modes per cavity frequency with perpendicular polarization vectors $\\epsilon$ and $\\bar\\epsilon$ having equal frequency and equal coupling strength $\\lambda_\\alpha=\\bar\\lambda_\\alpha$. With this choice, the bare cavity Hamiltonian $\\hat H_{\\rm bare} = \\sum_\\alpha \\omega_\\alpha(\\hat\\alpha^\\dagger\\hat\\alpha + \\bar{\\hat\\alpha}^\\dagger\\bar{\\hat\\alpha})$ is invariant under the continuous rotation operator $\\hat C_\\theta = \\exp[\\theta \\sum_\\alpha(\\hat\\alpha^\\dagger\\bar{\\hat\\alpha} - \\bar{\\hat\\alpha}^\\dagger\\hat\\alpha)]$; the displacement field $\\hat D$ itself is not invariant, but the bilinear coupling $\\hat D\\cdot\\hat d$ is, so the full polaritonic Hamiltonian retains the $\\mathrm{D}_{\\infty h}$ symmetry of the bare cavity. Treating the two polarizations in the coupled cluster ansatz (QED-CCSD-12-SD with the $\\hat\\Gamma_1$ and $\\hat\\Gamma_2$ photonic operators) and exploiting point-group symmetry through direct-product decomposition reduces the floating-point cost by a factor of $h^2$ (the group order) and, crucially, allows polaritonic excited states to be assigned to irreducible representations and to be targeted individually. The result is a symmetry-exact description of unpolarized cavities, which the paper uses to show that such cavities produce dense avoided-crossing landscapes, allowed crossings between states differing by more than one photon, and electron density shifts roughly twice as large as in a linearly polarized cavity.","pith_inferences":["If the two-mode model is correct, the same point-group logic extends to any number of degenerate mode pairs, so a multi-mode unpolarized cavity could be handled by the same symmetry-blocked machinery; the practical obstacle would then be the choice of relevant frequencies and the growing mass renormalization, which the paper leaves to future work.","The basis-independence shown in the appendix (real versus complex polarization vectors related by a unitary transformation) implies that the unpolarized-cavity description contains both linear and circular polarization descriptions as special bases, so computed polaritonic spectra should be invariant under that change—a direct check of the implementation.","The truncation-sensitive crossings suggest that predictions about photochemistry in the strong-coupling regime should be tested against calculations that include higher photonic excitations, since the present scheme may artificially allow crossings that a more complete treatment would convert into avoided crossings, changing the dynamical picture.","The equal-coupling condition $\\lambda_\\epsilon=\\lambda_{\\bar\\epsilon}$ is a testable experimental assumption: if a carefully characterized Fabry-Pérot cavity shows different coupling strengths for the two orthogonal polarizations, the rotational symmetry and the resulting degeneracies would be lifted, and the model would need to be relaxed."],"forward_implications":["In an unpolarized cavity, any orientation of an asymmetric molecule rotated about the cavity wave vector $\\mathbf k$ is degenerate in energy; for benzene this means the unpolarized cavity stabilizes the orientations with $\\mathbf k$ in the molecular plane as an infinite degenerate family, whereas a linearly polarized cavity selects a single orientation.","Point-group symmetry turns polaritonic excited states into labeled objects: two states mix (forming upper/lower polaritons) only if the direct product of their electronic and photonic irreducible representations contains the totally symmetric representation, so one can predict which transitions will show Rabi splitting.","States of the same irreducible representation can still cross if they differ by more than one photon, because the CCSD-12-SD similarity-transformed Hamiltonian has no two-photon creation with a single electronic de-excitation; this produces allowed crossings that a less truncated treatment (adding $\\hat\\Gamma_3$ or $\\hat S^2_1$) would turn into avoided crossings.","For the H$_2$ molecule, the unpolarized cavity is effectively a superposition of two orthogonal linearly polarized cavities, so the excited-state landscape contains degenerate pairs ($\\Pi_g$, $\\Delta_g$) that carry direct information about the rotational symmetry.","The ground-state electron density shifts in an unpolarized cavity are roughly twice as large as in a linearly polarized one at the same coupling strength, because two polarization components contribute to the interaction."],"supporting_citations":[{"why":"Defines the base polaritonic coupled cluster ansatz that this work generalizes to two perpendicular modes.","marker":"[11]"},{"why":"Previous treatment of polarization symmetry in chiral cavities that motivates the two-mode picture and the sigma_h transformation.","marker":"[26]"},{"why":"Provides the Pauli-Fierz / dipole-approximation Hamiltonian in the coherent-state basis from which the two-mode model is derived.","marker":"[7]"},{"why":"Identifies the Fabry-Pérot cavity as the experimental system being modeled, motivating the unpolarized two-mode description.","marker":"[28]"},{"why":"Supplies the direct-product decomposition technique for exploiting point-group symmetry in the coupled cluster equations.","marker":"[31]"},{"why":"Provides the QEDFT results for azulene in a linearly polarized cavity used as the comparison baseline for the density differences.","marker":"[15]"},{"why":"Describes the QED-CC amplitude and Lambda-equation formalism that the present implementation solves.","marker":"[39]"}],"fun_headline_variants":["Symmetry-exact polaritonic theory for unpolarized cavities","Two perpendicular modes preserve cavity symmetry in polariton CCSD","Point-group symmetry targets polaritonic states in unpolarized cavities","Unpolarized cavities get exact symmetry via two-mode coupled cluster","Two perpendicular polarizations make cavity theory symmetry-exact"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that an unpolarized Fabry-Pérot cavity is faithfully captured by two discrete modes with identical frequency and identical coupling strength $\\lambda$ to the molecule, with all other cavity modes, mode continua, and mass renormalization effects neglected — an approximation the paper states explicitly and flags as future work.","fun_headline_variants_meta":{"raw":{"variants":["Symmetry-exact polaritonic theory for unpolarized cavities","Two perpendicular modes preserve cavity symmetry in polariton CCSD","Point-group symmetry targets polaritonic states in unpolarized cavities","Unpolarized cavities get exact symmetry via two-mode coupled cluster","Two perpendicular polarizations make cavity theory symmetry-exact"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000699,"raw_usage":{"total_tokens":3207,"prompt_tokens":1045,"completion_tokens":2162,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":661,"completion_tokens_details":{"reasoning_tokens":2091}},"tokens_in":661,"tokens_out":2162,"duration_ms":16318,"temperature":1.0,"reasoning_tokens":2091,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:59:23.219588+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the two-mode QED-CC ground-state energy and first polaritonic gap for a small molecule such as H$_2$, then repeat with a larger set of modes at the same frequency and coupling; if the two-mode results shift by more than chemical accuracy, the discrete two-mode representation of the unpolarized cavity fails. A complementary test is to measure the coupling strengths for the two orthogonal polarizations of a real Fabry-Pérot cavity and check whether they are equal, since unequal $\\lambda_\\epsilon$ and $\\lambda_{\\bar\\epsilon}$ would lift the predicted rotational degeneracy of molecular orientations.","supporting_citations":[{"cited_title":"S.; Ronca, E.; Kj nstad, E","cited_arxiv_id":null,"evidence_quote":"Defines the base polaritonic coupled cluster ansatz that this work generalizes to two perpendicular modes."},{"cited_title":"T.; Angelico, S.; Kjønstad, E","cited_arxiv_id":null,"evidence_quote":"Previous treatment of polarization symmetry in chiral cavities that motivates the two-mode picture and the sigma_h transformation."},{"cited_title":"Understanding polaritonic chemistry from ab initio quantum electrodynamics","cited_arxiv_id":null,"evidence_quote":"Provides the Pauli-Fierz / dipole-approximation Hamiltonian in the coherent-state basis from which the two-mode model is derived."},{"cited_title":"Resolution of gauge ambiguities in ultrastrong-coupling cavity quantum electrodynamics","cited_arxiv_id":null,"evidence_quote":"Identifies the Fabry-Pérot cavity as the experimental system being modeled, motivating the unpolarized two-mode description."},{"cited_title":"Ab Initio Optimized Effective Potentials for Real Molecules in Optical Cavities: Photon Contributions to the Molecular Ground State","cited_arxiv_id":null,"evidence_quote":"Provides the QEDFT results for azulene in a linearly polarized cavity used as the comparison baseline for the density differences."},{"cited_title":"R.; Barlini, A.; Ronca, E.; Koch, H","cited_arxiv_id":null,"evidence_quote":"Describes the QED-CC amplitude and Lambda-equation formalism that the present implementation solves."}],"review_version":2}