{"id":"91d08fc5-1439-432c-8184-796433be6eee","arxiv_id":"2506.18647","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A semiclassical cavity-QED simulation shows that coupling a cavity mode to a vibronic progression or to a hydrogen tunneling transition modifies molecular spectra and tunneling dynamics beyond the Born-Oppenheimer approximation.","lead":"This paper simulates a molecule in an optical cavity, where the cavity light field interacts with both electrons and a quantum proton. The simulations show that the cavity can couple to nuclear motions usually associated with electronic transitions and can modify hydrogen tunneling dynamics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tunneling simulations likely run in the ultrastrong-coupling regime, where the paper's neglect of the self-dipole term and its classical-cavity mean-field are both quantitatively unjustified.","rationale":"The reader's weakest assumption—mean-field treatment of the cavity—is correct in spirit, but the most acute version is in the tunneling section, where the paper's own parameters enter a regime that undermines its stated justification for two nested approximations. First, the dimensionless coupling g/omega estimated from the reported frequency and coupling strengths is not small; for d ~ 1-2 a.u., g/omega is ~0.3-1. This is a known ultrastrong-coupling range, where rotating-wave approximations fail and the mu^2 diamagnetic term is required for gauge invariance. The paper explicitly drops this term (Section II A) and claims strong coupling, but that claim is inconsistent with the tunneling parameters. Second, the classical-cavity mean-field in Eqs. 8a-8b is a self-consistent Ehrenfest approximation. In the ultrastrong regime, the Pauli-Fierz ground state contains virtual photons and is entangled; a classical trajectory cannot represent photon-number fluctuations. Both effects affect the tunneling dynamics, which is the paper's most significant and externally relevant claim. The self-dipole term, when included, modifies the proton potential because mu^2 differs between the delocalized bridge and the localized wells, so it can directly change the tunneling splitting and the energy exchange with the cavity. Without a rerun including this term, or a benchmark against a quantized-photon method, the claim that coupling a molecule to a cavity mode can alter hydrogen tunneling dynamics remains conditional. This does not invalidate the paper as a whole: the vibronic strong coupling results (Sec. IV A) use omega ~ 19 eV, where g/omega ~ 0.01-0.04, so the truncations are well controlled there. It means the tunneling claim should be viewed as unresolved. Since the reader already issued CONDITIONAL, the verdict is unchanged. Agreement with the reader: partial—the reader identified the mean-field as the weak point; we add that the parameter regime for claim 3 makes both the mean-field and the self-dipole truncation quantitatively doubtful, and we provide a concrete test.","tokens_in":17498,"tokens_out":15263,"duration_ms":156575,"concrete_test":"Repeat the Sec. IV B tunneling simulations (epsilon = 5e-5 and 10e-5 a.u., omega = 552 cm^-1) with the self-dipole term (epsilon^2/(2 * omega^2)) * <mu^2> added to the semiclassical Hamiltonian, recomputing the NEO-CASCI states and propagating the same initial superposition and cavity conditions. If the characteristic damping/phase-reversal pattern in Figs. 6-7 changes qualitatively, or if the tunneling splitting shifts by more than about 20%, the central claim depends on the neglected term and is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most consequential claim—that coupling to a cavity can alter hydrogen tunneling dynamics (Sec. IV B)—is the least protected by the paper's stated approximations. The paper neglects the mu^2 self-dipole term in Eq. 2, justifying this by saying only strong, not ultrastrong, coupling is considered. For the tunneling simulations, omega = 552 cm^-1 (0.0025 a.u.) and epsilon = 5e-5 and 10e-5 a.u. With a proton-transfer dipole matrix element d of order 1-2 a.u., the dimensionless coupling g/omega = (epsilon * d) / sqrt(2 * omega^3) is roughly 0.3-1.1, squarely in the ultrastrong regime (g/omega > 0.1). There, the self-dipole coefficient epsilon^2/(2 * omega^2) is about 8e-4 a.u., and multiplied by <mu^2> ~ 4 a.u.^2 it yields an energy around 660 cm^-1, comparable to the tunneling splitting itself; this term modifies the proton potential and can change the splitting. The classical-cavity mean-field (Eqs. 8a-8b) is also least tenable here, since photon-number fluctuations and light-matter entanglement are significant when g/omega is order unity. Thus the tunneling result may reflect the truncated Hamiltonian rather than genuine cavity modification. The paper explicitly acknowledges treating the cavity classically but gives no benchmark or error estimate for these parameters.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the nuclear-electronic orbital time-dependent configuration interaction (NEO-TDCI) method, in which electrons and a specified proton are quantized on equal footing without a Born-Oppenheimer separation, to include a cavity mode treated as a classical harmonic oscillator. The cavity is driven by the expectation value of the total molecular dipole (Eqs. 8a-8b) and feeds back on the molecule through a time-dependent dipole coupling term in the semiclassical Hamiltonian (Eq. 4). The method is applied to HeHHe+ (with HCN in the Supplementary Material). In Section IV A, the cavity is tuned to three vibronic transitions of the S0->S1 progression (omega_0, omega_2, omega_4) at three coupling strengths, and polaritonic spectra are decomposed into molecular and cavity components using projected coefficients and Fourier transforms of the cavity displacement. In Section IV B, the cavity is tuned to the computed 552 cm^-1 tunneling splitting, and two coupling strengths produce damped, phase-reversed proton dipole oscillations attributed to reversible energy exchange between the tunneling doublet and the cavity. The central claims are that this ab initio approach treats electronic and vibrational strong coupling in a unified way, that the cavity couples to nuclear motion even at ESC-scale cavity frequencies, and that coupling to a cavity can alter hydrogen tunneling dynamics.","tokens_in":17757,"tokens_out":9946,"duration_ms":93753,"significance":"The paper has several genuine strengths. The cavity-free tunneling propagation is validated against the analytic two-level solution (Fig. 6c), the Rabi-splitting patterns and their asymmetry are plausibly interpreted, and the projected-coefficient analysis provides a practical diagnostic for decomposing polaritonic states when the cavity is classical. The ability to treat a quantized proton and electrons on the same level, and to capture vibronic progressions with multiple nuclear quanta at first-principles level, goes beyond RT-NEO-TDDFT and is a useful methodological contribution. If the tunneling result survives a proper treatment of the cavity, the demonstration that a cavity can modify tunneling dynamics in an ab initio non-Born-Oppenheimer framework would be significant for polaritonic chemistry. However, the two most consequential claims, especially the cavity-modified tunneling, rest on the semiclassical mean-field treatment of the photon and on the neglect of the self-dipole term, and the parameter regime used for the tunneling simulations is exactly where those approximations are most vulnerable.","major_comments":[{"comment":"The tunneling simulations are run in the ultrastrong-coupling regime, which contradicts the paper's own justification in Section II A for neglecting the self-dipole term. With omega = 552 cm^-1 (0.0025 a.u.) and epsilon = 5e-5 and 1e-4 a.u., and a proton-transfer dipole d of order 1-2 a.u., the dimensionless ratio epsilon*d/sqrt(2*omega^3) is roughly 0.3-1.1, i.e., g/omega well above the ultrastrong threshold of 0.1. The neglected self-dipole coefficient epsilon^2/(2*omega^2) is 2e-4 to 8e-4 a.u.; multiplied by <mu^2> ~ 4 a.u.^2 this contributes 0.02-0.09 eV, comparable to the 552 cm^-1 (0.068 eV) tunneling splitting that is the target energy scale. Since this term modifies the proton potential, it can change the tunneling splitting and the cavity-modified dynamics at the same order as the effect being reported. The authors should include or otherwise quantitatively bound this term, or restrict the tunneling simulations to genuinely strong-coupling parameters.","section":"Sec. IV B; Eqs. (2) and (4)"},{"comment":"The classical-cavity mean-field treatment is least defensible in exactly the tunneling regime where g/omega is of order 0.3-1.1: photon-number fluctuations and light-matter entanglement are then significant, and the reversible damping and phase reversal seen in Figs. 6a,b could be an artifact of replacing the photon operator by its expectation value. Reference [58] (full-quantum RT-NEO-TDDFT) provides a natural benchmark, but no comparison of the semiclassical and fully quantized descriptions is reported for any system or parameter set. I would require, at minimum, a model-system benchmark (e.g., a two-level system plus a single quantized mode with the same g/omega) demonstrating that the mean-field error is small at the parameters used, or a full-quantum NEO calculation for the tunneling case; without this, the load-bearing claim that a cavity alters hydrogen tunneling is not established.","section":"Sec. II C, Eqs. (8a)-(8b); Sec. IV B"},{"comment":"No convergence tests are reported for the NEO-CASCI active spaces (4 electrons in 8 orbitals, 1 proton in 18 orbitals), for the protonic basis, or for the real-time propagation timesteps (0.005/0.01 a.u.), although quantitative outputs, including Rabi splittings, peak positions, and the 552 cm^-1 tunneling splitting used to set the cavity resonance, are central to the paper's claims. The manuscript should include at least a check that the relevant vibronic energies and the tunneling splitting are converged with respect to the active space and that the reported spectra are converged with respect to timestep and trajectory length.","section":"Sec. III"}],"minor_comments":[{"comment":"The text 'the hybrid light-matter nature of the the polaritonic states' contains a duplicated 'the'.","section":"Sec. IV A 2"},{"comment":"References 16 and 29 contain typographical or bibliographic errors ('optical eavities' in Ref. 16; a bare DOI-like string in the page field of Ref. 29) and should be corrected.","section":"References"},{"comment":"In Eq. (10), the adjoint symbol on the scalar coefficient C_{i,j}^{(n)} should be a complex conjugate; please clarify the notation.","section":"Sec. II D, Eq. (10)"},{"comment":"The statement that 'significantly smaller light-matter couplings than those used in the vibronic strong coupling case are required to observe cavity effects' would benefit from the quantitative dipole argument (the tunneling transition dipole is large), since the choice of epsilon = 5e-5 to 1e-4 a.u. currently appears ad hoc.","section":"Sec. IV B"},{"comment":"The caption states 'All plots in this row are identical'; it would be clearer to say that the four panels in row (d) are the same cavity-free spectrum shown for reference.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The main gap between the paper's claims and its evidence is the regime of validity of the semiclassical approximations in the tunneling section, which I flagged as a major concern. The paper's own justification for dropping the self-dipole term ('only strong coupling') is inconsistent with the parameters used in Section IV B; I treat this as a correctness risk rather than a presentation issue. The natural remedy, a benchmark against the group's own full-quantum RT-NEO-TDDFT [58] or a model two-level-plus-mode system, is within scope and should be requested. The manuscript is otherwise a good fit for a chemical-physics journal, and the extended NEO-TDCI framework with the polaritonic decomposition analysis adds genuine value beyond the existing RT-NEO-TDDFT cavity work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: this is a legitimate step forward, not a revolution. The genuinely new content is the combination of NEO-TDCI with a classical cavity mode to simulate vibronic strong coupling across multiple transitions and cavity-modified hydrogen tunneling in HeHHe+. The underlying pieces are the authors' own prior work, but the joint treatment of electrons and quantized protons without the Born-Oppenheimer separation lets them describe eV-scale vibronic progressions and cm^-1-scale tunneling in one framework. That is useful.\n\nThe vibronic section is the stronger half. The cavity-free dynamics is checked against the analytic two-level solution, the Rabi splittings are sensible, and the projected-coefficient decomposition into molecular and cavity character is a nice diagnostic. I found the explanation of the splitting asymmetries using nearby vibronic states convincing. The simulations are not fitting to target spectra; they are first-principles outputs, and the citation pattern looks fair.\n\nThe soft spot is real, and the stress-test note is right. In the tunneling simulations, omega = 552 cm^-1 and epsilon = 5e-5 or 1e-4 a.u. With a proton-transfer dipole matrix element of order 1–2 a.u., the dimensionless coupling g/omega is order 0.4–1.1. That is ultrastrong coupling by the usual definition. The paper explicitly justifies dropping the mu^2 self-dipole term because only strong coupling is considered, and then runs the tunneling calculation in precisely the regime where that term can shift energies by several hundred cm^-1, comparable to the tunneling splitting itself. The classical mean-field for the cavity is also on its weakest footing there, since photon-number fluctuations and light-matter entanglement should matter. The paper acknowledges the classical treatment but offers no benchmark against a quantized-cavity calculation or a calculation that includes the self-dipole term. So the claim that the cavity alters tunneling dynamics is plausible, and the damped/revival behavior is interesting, but it is not quantitatively secured. That is the load-bearing issue for the headline claim.\n\nMinor points: no convergence tests for basis/active space/timestep, and no public code or data yet, only a promise of Zenodo at publication. These are fixable.\n\nThe paper deserves a serious referee. I would not desk reject it. But I would ask the authors to either include the self-dipole term in the tunneling section or explicitly qualify those results as semiclassical/ultrastrong-coupling results with unknown DSE corrections. The vibronic section stands on its own.","headline":"The vibronic-strong-coupling part is solid; the tunneling part runs in the ultrastrong-coupling regime where the paper's own approximations are unbenchmarked.","tokens_in":18347,"tokens_out":4368,"would_cite":true,"duration_ms":46841,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Semiclassical NEO-TDCI treats electrons and a quantized proton on the same footing, letting one ab initio theory describe vibronic strong coupling and cavity-modified hydrogen tunneling.","keywords":["vibronic strong coupling","nuclear-electronic orbital","time-dependent configuration interaction","polaritonic chemistry","hydrogen tunneling","semiclassical cavity","non-Born-Oppenheimer dynamics","cavity-modified reactions"],"falsifier":"Perform the same HeHHe+ simulations with a fully quantized cavity mode at the same frequencies and coupling strengths; if the Rabi splittings, peak positions, or the timing and phase of the tunneling dipole oscillations differ by more than the spectral line widths, the semiclassical mean-field feedback is missing quantum correlations.","tokens_in":17245,"feed_emoji":"⚛️","tokens_out":7010,"duration_ms":66988,"temperature":0.7,"pith_summary":"This paper implements the semiclassical nuclear-electronic orbital time-dependent configuration interaction (NEO-TDCI) approach, in which electrons and a selected proton are quantized at the same level while the cavity mode is propagated as a classical harmonic oscillator. The aim is to show that one ab initio framework can describe electronic strong coupling and vibrational strong coupling without a Born-Oppenheimer separation between electrons and nuclei. Simulating the HeHHe+ cation, the authors find that a cavity tuned to different vibronic subpeaks produces polaritonic spectra with Rabi splittings, and that the cavity mode couples to nuclear motion even at frequencies normally associated with electronic transitions. They also find that coupling a cavity to the tunneling splitting near $\\sim 552\\ \\mathrm{cm}^{-1}$ changes the proton tunneling dynamics through reversible energy exchange between molecule and cavity. If correct, the method offers a first-principles route to polariton-modified chemistry where tunneling and nonadiabatic effects matter.","feed_headline":"One simulation method spans electronic and vibrational strong coupling","feed_subtitle":"Electrons and a tunneling proton are quantized together while the cavity stays classical.","key_machinery":"The central object is the semiclassical NEO-TDCI wavefunction coupled to a classical cavity mode. NEO-TDCI is a configuration interaction expansion over products of electronic Slater determinants and protonic Slater determinants, so a vibronic excitation is a joint electronic-nuclear excitation with no Born-Oppenheimer separation. The cavity mode follows classical oscillator equations of motion driven by the change in the total molecular dipole expectation value, and the wavefunction is propagated under the resulting time-dependent semiclassical Hamiltonian. A projected-coefficient analysis maps the time-dependent wavefunction onto cavity-free molecular vibronic states, while the Fourier transform of the cavity displacement identifies which spectral peaks carry field character.","core_discovery":"The central claim is that the same semiclassical NEO-TDCI simulation can capture vibronic strong coupling and cavity-modified hydrogen tunneling. In the vibronic case, with the cavity tuned to $\\omega_0$, $\\omega_2$, or $\\omega_4$ (the $S_0\\nu_0\\to S_1\\nu_0$, $S_0\\nu_0\\to S_1\\nu_2$, and $S_0\\nu_0\\to S_1\\nu_4$ transitions), the spectra show Rabi splittings that grow with coupling strength, and the projected-coefficient and field-displacement analysis shows that the polaritonic states mix many molecular vibronic states. At larger coupling nearly all transitions acquire cavity character, so a two-level Jaynes-Cummings model is insufficient. In the tunneling simulations, the cavity mode exchanges energy with the symmetric and antisymmetric proton vibronic states, damping the proton oscillation and reversing its phase before returning the energy. On this basis the paper claims that the cavity couples to nuclear motion even at ESC-associated frequencies and that a non-Born-Oppenheimer electron-proton treatment is necessary for polaritonic dynamics.","pith_inferences":["If the semiclassical feedback is reliable, the phase reversals in the tunneling dipole identify the cavity as a temporary energy reservoir; adding cavity loss or a donor-acceptor distance coordinate could convert this into net reaction-rate suppression or enhancement.","Applying the same approach to deuterated isotopologues would shift the tunneling splitting, so the cavity frequency and coupling strength needed to alter tunneling would shift in a testable way.","A fully quantized photon calculation at the same couplings would quantify how much photon-number correlations contribute beyond the mean-field feedback; the paper does not provide that benchmark."],"forward_implications":["Vibronic strong coupling can be simulated from first principles for realistic molecules, avoiding model Hamiltonians and diabatization.","The same method covers ESC- and VSC-scale physics without choosing an adiabatic representation for the nuclei.","A cavity mode tuned to an electronic transition can still couple to nuclear motion through the vibronic progression.","Cavity-modified hydrogen tunneling appears as reversible energy exchange, implying that coupling strength and cavity frequency control the tunneling dynamics.","At stronger coupling, polaritonic spectra require many molecular vibronic states, not just the resonant two-level transition."],"supporting_citations":[{"why":"Establishes the semiclassical cavity coupling prescription, treating the cavity as a classical oscillator driven by the time-dependent molecular dipole.","marker":"[57]"},{"why":"Supplies the NEO-TDCI method that produces vibronic progressions and real-time hydrogen tunneling dynamics.","marker":"[59]"},{"why":"Provides the nuclear-electronic orbital framework in which electrons and quantized nuclei are treated on equal footing.","marker":"[54]"},{"why":"Provides the full-quantum RT-NEO-TDDFT counterpart that the semiclassical approach is contrasted with.","marker":"[58]"},{"why":"Gives the length-gauge QED Hamiltonian that defines the cavity-matter coupling used here.","marker":"[31]"},{"why":"Offers the NEO multireference configuration interaction route cited for quantitatively accurate tunneling splittings.","marker":"[77]"},{"why":"Documents experimental vibronic strong coupling that motivates the modeled regime.","marker":"[61]"}],"fun_headline_variants":["Cavity tweaks proton tunneling in unified strong-coupling simulation","Same simulation captures electronic and vibrational cavity coupling","Hydrogen tunneling shifts when cavity mode joins the simulation","Semiclassical method unites electronic and vibrational strong coupling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing approximation is that the cavity mode can be treated as a classical oscillator driven by the expectation value of the molecular dipole, so quantum light-matter correlations and photon number fluctuations are assumed negligible at the coupling strengths used.","fun_headline_variants_meta":{"raw":{"variants":["Cavity tweaks proton tunneling in unified strong-coupling simulation","Same simulation captures electronic and vibrational cavity coupling","Hydrogen tunneling shifts when cavity mode joins the simulation","Semiclassical method unites electronic and vibrational strong coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000462,"raw_usage":{"total_tokens":2337,"prompt_tokens":995,"completion_tokens":1342,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":1276}},"tokens_in":611,"tokens_out":1342,"duration_ms":10499,"temperature":1.0,"reasoning_tokens":1276,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:45:14.578747+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the same HeHHe+ simulations with a fully quantized cavity mode at the same frequencies and coupling strengths; if the Rabi splittings, peak positions, or the timing and phase of the tunneling dipole oscillations differ by more than the spectral line widths, the semiclassical mean-field feedback is missing quantum correlations.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the semiclassical cavity coupling prescription, treating the cavity as a classical oscillator driven by the time-dependent molecular dipole."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the NEO-TDCI method that produces vibronic progressions and real-time hydrogen tunneling dynamics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the nuclear-electronic orbital framework in which electrons and quantized nuclei are treated on equal footing."},{"cited_title":"Light-Matter Entanglement in Real-Time Nuclear-Electronic Orbital Polariton Dynamics","cited_arxiv_id":"2506.06490","evidence_quote":"Provides the full-quantum RT-NEO-TDDFT counterpart that the semiclassical approach is contrasted with."},{"cited_title":"Flick , author M","cited_arxiv_id":null,"evidence_quote":"Gives the length-gauge QED Hamiltonian that defines the cavity-matter coupling used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Offers the NEO multireference configuration interaction route cited for quantitatively accurate tunneling splittings."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents experimental vibronic strong coupling that motivates the modeled regime."}],"review_version":2}