{"id":"dfc68418-51c0-4bd1-97c3-33e80df427a5","arxiv_id":"2506.11963","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Phase-space electronic structure theory predicts a three-dimensional branching plane and stable electronic momentum states at the BeH2 conical intersection, explaining complex Hartree-Fock instabilities.","lead":"Using a phase-space electronic Hamiltonian that couples nuclear momentum to electrons, the authors find that a conical intersection in BeH2 has a three-dimensional branching plane and stationary electronic states with nonzero electronic momentum. The results suggest that complex restricted Hartree-Fock instabilities near conical intersections have physical meaning.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central results hinge on the untested Γ'' angular-momentum coupling; the 3D branching plane and P_min double well could vanish if Γ'' is modified.","rationale":"The reader identified the general form of H_PS as the weakest assumption. I agree but make the concern more precise: the only source of a complex H12 near the CI is the cross term -2iℏP·Γ, and within the parametrization used this term is dominated by the ad hoc rotational operator Γ''. The paper's symmetry arguments establish that an off-diagonal matrix element can be nonzero, but not that it is of the right magnitude or even nonzero for the chosen Γ''. The double-well minimum position and the existence of the third branching direction are therefore not settled by the current data. The proposed computational test is a direct sensitivity check using the paper's own methods: removing Γ'' or changing the partition functions would reveal whether the qualitative claims are tied to the arbitrary parameters. The verdict remains CONDITIONAL because the paper's qualitative conclusions are plausible and the CRHF agreement provides some independent (if qualitative) support, but the model dependence must be demonstrated before the paper can be accepted.","tokens_in":741,"tokens_out":1791,"duration_ms":214742,"concrete_test":"Recompute the branching modes (Fig. 2) and the P-space scan (Fig. 3) for BeH2 with (i) Γ''=0 (only Γ' from Eq. 18) and (ii) two alternative choices of the partition functions θ_A (e.g., Becke fuzzy-cell weights vs. Voronoi) while keeping everything else fixed. If the f-branching mode vanishes, or if the minimum of the momentum scan moves to P=0, then the claimed 3D branching plane and P_min≠0 double well are artifacts of the arbitrary Γ'' parametrization, not robust predictions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Both headline results (three-dimensional branching plane, double-well minima at ±P_min ≠ 0) follow from the first-order term in P of H_PS, i.e., -2iℏ P·Γ. For the two degenerate CI states in BeH2, this term yields a nonzero Im H12 only through the matrix element ⟨Φ1|Γ|Φ2⟩. The Γ operator is not derived from the molecular Hamiltonian; Eq. 19's Γ'' term is an ad hoc angular-momentum operator containing position-dependent partition functions θ_A, locality parameters ζ_AB, and an atomically partitioned inertia tensor K_B. The paper provides no sensitivity analysis: it never reports whether the f-branching mode or the momentum double well persists when these arbitrary parameters are varied (or when Γ'' is removed). Since Γ'' is also the term that couples to rotation and produces the circulating currents of Fig. 4, the observed double well and the CAS-CI/CRHF agreement could be consequences of this particular parametrization rather than of the underlying phase-space physics. The branching analysis and momentum scan need to be repeated with different Γ'' choices to establish robustness.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates the conical intersection in BeH2 using a phase-space electronic Hamiltonian H_PS(R,P) that depends on both nuclear positions and momenta. The main claims are that (i) within this framework, the conical intersection's branching plane has dimension three, rather than the usual two-dimensional branching plane in Born-Oppenheimer theory, because the phase-space Hamiltonian is complex and requires an additional degeneracy condition Im H12 = 0; and (ii) at a geometry on the crossing seam, the ground-state energy as a function of nuclear momentum exhibits a double well with minima at ±P_min ≠ 0, indicating that the stationary electronic states carry finite electronic (angular) momentum. The authors further show that the electronic current densities obtained from phase-space CAS-CI calculations at P = ±P_min resemble those of complex restricted Hartree-Fock (CRHF) solutions at P = 0, and they propose a diabatization scheme based on diagonalizing L_z. The paper concludes with speculations about the dynamical and experimental consequences of these symmetry-broken states.","tokens_in":15459,"tokens_out":9537,"duration_ms":116142,"significance":"If the results hold, they would extend the concept of conical intersections into a phase-space framework and offer a physical interpretation of CRHF instabilities, which are often regarded as artifacts of the single-determinant approximation. The paper provides a concrete demonstration on a well-studied system (BeH2), with CAS-CI and CRHF calculations, and the qualitative agreement between the two levels is a valuable independent benchmark. The derivations of the branching-plane conditions and the double-well structure are internally consistent. However, the central predictions are contingent on the specific form of the phase-space Hamiltonian and its parameters; the paper does not establish robustness of these predictions to variations in the model, which tempers its immediate significance.","major_comments":[{"comment":"The central results, i.e., the three-dimensional branching plane and the double-well minima at ±P_min, are driven by the first-order term -2iℏ P·Γ in Eq. (17). The imaginary part of H12, which generates the f-branch and the momentum-space double well, is proportional to ⟨Φ1|Γ|Φ2⟩, and the operator Γ'' in Eq. (19) contains the ad hoc parameters θ_A, ζ_AB, and K_B. The paper provides no sensitivity analysis: it does not report whether these features persist when Γ'' is modified or removed. Because these predictions are the paper's principal claims, the authors should repeat the branching-plane scan and the momentum scan for at least a few different parameter choices (e.g., varying ζ_AB or setting Γ'' = 0) to demonstrate that the results are not artifacts of this particular parametrization.","section":"§III.A.2, Eqs. (17)-(19)"},{"comment":"The calculations use a modest 6-31G basis and a fixed active space (all orbitals except Be 1s). No basis-set or active-space convergence tests are reported. The quantitative values of P_min, the barrier height of the double well, and even the existence of the f-branch could depend on these choices. The authors should show that the qualitative features are stable when the basis is enlarged (e.g., to cc-pVDZ) or when the active space is expanded, to rule out basis-set artifacts.","section":"§III.A, Fig. 3"},{"comment":"The abstract states that the CAS-CI electronic momentum 'agrees with' the CRHF prediction, but the evidence presented is only qualitative. The paper does not provide numerical values for ⟨L_z⟩ from CAS-CI and CRHF, nor a quantitative comparison of P_min values or energy differences. To support the claim of agreement, the authors should report these numbers, along with the definitions used to compute P_min for the CRHF solutions. Without such numbers, the strength of the claimed agreement is unclear.","section":"§IV, Figs. 5-6"}],"minor_comments":[{"comment":"The abstract refers to 'full configuration interaction', but the calculations are CAS-CI with a frozen Be 1s orbital; this discrepancy should be corrected.","section":"Abstract and §III.A"},{"comment":"The definition of N_int^PS = 6N - 6 is not derived; the counting that leads to this expression and the role of translational/rotational momentum should be clarified.","section":"§III.A.2"},{"comment":"There are several typos and unclear phrases: 'we have still have a lot to learn' (abstract), 'are are parameterized' (Sec. I.C), 'is valid only valid' (Sec. I.B.1), and 'spin-broken symmetry states' (Sec. III.B) where 'time-reversal broken' or 'angular momentum broken' is meant.","section":"Throughout"},{"comment":"The sentence 'the minima is not always at P = 0' is grammatically incorrect and should be rephrased. Also, the reference to '2a direction' is ambiguous; figure panels should be cited explicitly.","section":"§III.B"},{"comment":"The claim that the phase-space Hamiltonian 'captures a non-trivial amount of the total non-adiabatic interaction' is not quantified; a reference to the specific earlier result would help.","section":"§II"}],"recommendation":"major_revision","confidential_remarks":"This is a well-written paper on a timely topic, but the central results are highly dependent on the assumed phase-space Hamiltonian, which is not derived within the manuscript. The lack of sensitivity and convergence analyses is a serious concern that should be addressed before publication. The paper might also benefit from a more direct quantitative comparison with the CRHF benchmark to substantiate the headline claim of agreement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one if you care about nonadiabatic dynamics or symmetry breaking in electronic structure. The group applies their phase-space Hamiltonian (built in earlier papers) to the BeH2 conical intersection and finds two things that are new: the branching plane becomes three-dimensional once nuclear momentum is included, and the electronic ground state at the seam develops a double well in P with minima at nonzero electronic momentum. They also show that CRHF solutions at P=0 carry currents resembling those of the phase-space CAS-CI states at P_min, and they interpret CRHF instabilities as physical manifestations of electronic momentum. That last step is the most valuable part of the paper, and it is a genuinely non-obvious connection.\n\nWhat the paper does well: the derivation of the branching conditions is transparent, the CAS-CI and CRHF calculations are consistent with each other, and the current-density plots make the symmetry breaking visually clear. The claim that a time-reversal-breaking term adds a third branching direction is mathematically correct for a complex Hamiltonian. There is no sleight of hand in the internal logic.\n\nThe soft spot is exactly what the stress-test note flags. The Γ operator, especially the Γ'' angular-momentum term in Eq. 19, is the load-bearing piece: it is what produces the imaginary coupling that generates the f-branch and the P double well. It is carried over from prior work, not derived in this manuscript, and the paper gives no sensitivity analysis. If that coupling matrix element were zero at the CI, the double well would vanish. The paper should show whether P_min and the f-branch survive when the parameters in Γ'' are varied or when Γ'' is dropped. That is a legitimate referee request, not a fatal flaw. The qualitative phenomenon is probably robust as long as some momentum-dependent coupling breaks time reversal, but the quantitative results are only as good as the parametrization.\n\nOther, smaller issues: the calculations are for one molecule, in a 6-31G basis, with a fixed active space and no convergence checks. That is standard for this kind of exploratory study, but it means the numerical P_min values should not be taken as predictions. The CRHF agreement is suggestive but not a validation of the exact Γ form; it is an independent benchmark that the symmetry breaking is not purely a single-determinant artifact.\n\nBottom line: yes, I would send this to a serious referee. It is a thought-provoking application with a testable interpretation, but the referee should require a robustness check on the Γ parametrization before publication. It deserves a full review in J. Chem. Phys. or similar.","headline":"A clean application of an existing phase-space Hamiltonian to a conical intersection, with a genuinely clever CRHF comparison, but the headline results lean on an untested parametrized operator that needs a sensitivity analysis.","tokens_in":15990,"tokens_out":3058,"would_cite":true,"duration_ms":44632,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"At a conical intersection, including nuclear momentum makes the branching plane three-dimensional and gives electrons nonzero momentum.","keywords":["conical intersection","phase space electronic structure","electronic momentum","branching plane","complex restricted Hartree-Fock","Born-Oppenheimer approximation","nonadiabatic dynamics","time-reversal symmetry breaking"],"falsifier":"Compute exact beyond-Born-Oppenheimer electronic momentum and current for a wavepacket traversing the BeH2 conical intersection at λ = 2.600 Bohr; if the emerging electronic state shows no net angular momentum and a momentum scan of the phase-space surface shows no stable double well at ±P_min, the central claim fails. A direct measurement showing Im(H12) = 0 throughout the seam would likewise contradict the third branching direction.","tokens_in":1596,"feed_emoji":"⚛️","tokens_out":2043,"duration_ms":82235,"temperature":0.7,"pith_summary":"This paper argues that the standard Born-Oppenheimer picture of conical intersections is incomplete because it sets electronic momentum to zero. Using a phase-space electronic Hamiltonian parameterized by both nuclear position and momentum, the authors show for BeH2 that the branching plane of the intersection seam becomes three-dimensional, not two. Fixing a geometry on the seam and scanning momentum reveals a double well whose minima sit at ±P_min, meaning the stationary electronic states carry real electronic angular momentum. They also show that this momentum, computed with full configuration interaction, matches what complex restricted Hartree-Fock predicts, suggesting that CRHF symmetry breaking near crossings has a physical basis. The larger point is that electronic momentum must be considered when studying photochemistry at conical intersections.","feed_headline":"Conical intersections get a third branching direction","feed_subtitle":"Phase-space theory gives electrons momentum at crossings and makes complex Hartree-Fock solutions physically meaningful.","key_machinery":"The load-bearing object is the phase-space electronic Hamiltonian $\\hat{H}_{PS}(\\mathbf{R},\\mathbf{P}) = \\sum_A \\frac{1}{2M_A}(\\mathbf{P}_A - i\\hbar \\hat{\\Gamma}_A(\\mathbf{R}))^2 + \\hat{H}_{el}(\\mathbf{R})$, where $\\hat{\\Gamma}_A$ is a one-electron operator that drags electronic density in the direction of nuclear momentum without diverging at crossings. $\\hat{\\Gamma}_A$ breaks time-reversal symmetry, allowing $\\mathrm{Im}(H_{12}) \\neq 0$ and producing the third branching direction; the double well in $\\mathbf{P}$ follows from mixing the two degenerate states with a complex phase. The paper uses BeH2 as a prototype, comparing CAS-CI (exact in the active space) with CRHF to connect the phase-space prediction to standard electronic-structure instabilities.","core_discovery":"The central discovery is that within phase-space electronic structure theory, the conical intersection between the ground and first excited singlet states of BeH2 is characterized by a three-dimensional branching plane generated by the gradients f, g, and h, corresponding to Im(H12), H11-H22, and Re(H12), whereas Born-Oppenheimer theory has only two. Because the P·Γ coupling breaks time-reversal symmetry, H12 can be complex, and Im(H12)=0 supplies the third degeneracy condition. At a fixed point on the CI seam, scanning nuclear momentum shows two degenerate minima at ±P_min, and the CAS-CI eigenstates display rotating electronic current density. The same currents appear in CRHF solutions at P=0, with ⟨L_z⟩ ≈ ±0.9ℏ, so the paper interprets complex Hartree-Fock instabilities as phase-space physical states carrying electronic angular momentum.","pith_inferences":["The three-dimensional branching plane suggests that momentum-space scans should be included when mapping conical intersection seams, and photochemical observables sensitive to electronic currents might distinguish the ±P_min states.","The CRHF-to-phase-space correspondence raises the possibility that other complex or triplet instabilities in linear-response theories also carry phase-space meaning rather than being mere artifacts.","If electronic angular momentum is generated near crossings, conserving total angular momentum will force nuclear rotation or spin polarization, potentially connecting these results to spin-selective chemistry, though the paper does not establish that link."],"forward_implications":["Nonadiabatic dynamics near conical intersections should include electronic momentum, since Born-Oppenheimer surfaces miss a physical torque on the electronic density.","For any even-electron molecule with a conical intersection, the phase-space branching plane has dimension three, shifting the seam dimension from N_int - 2 to N_int - 3 in phase-space internal coordinates.","Complex restricted Hartree-Fock instabilities near crossings can be viewed as capturing real electronic momentum, offering a cheap probe of these effects.","Diagonalizing L_z provides a phase-space route to diabatic states whose energies are nearly constant in R, giving a new way to build diabatic Hamiltonians around crossings.","If a wavepacket crosses a conical intersection on a single adiabatic surface, electronic angular momentum may be generated, requiring compensation by nuclear or spin angular momentum."],"supporting_citations":[{"why":"Introduces the original phase-space electronic Hamiltonian with derivative couplings, whose divergence near crossings motivates the use of the non-divergent Γ operator.","marker":"[42]"},{"why":"Defines the practical phase-space electronic Hamiltonian with a one-electron Γ operator used throughout this paper.","marker":"[43]"},{"why":"Provides the one-electron expression for electron rotational factors that underlies the angular-momentum part of Γ.","marker":"[44]"},{"why":"Gives the basis-free phase-space Hamiltonian that recovers electronic momentum and current density, validating the Γ form.","marker":"[45]"},{"why":"Shows that the Γ-based phase-space Hamiltonian captures nonadiabatic interaction in semiclassical dynamics, supporting its use near crossings.","marker":"[46]"},{"why":"Reports spin Coriolis symmetry breaking in phase-space theory, the spin analogue that motivates checking symmetry breaking without spin.","marker":"[48]"},{"why":"Supplies the BeH2 insertion model system with the conical intersection studied here.","marker":"[49]"},{"why":"Establishes complex restricted Hartree-Fock theory, the approximate method whose P=0 currents are compared with CAS-CI.","marker":"[50]"},{"why":"Proves that real electronic Hamiltonians make H12 real in Born-Oppenheimer theory, the fact that the phase-space Γ coupling overturns.","marker":"[34]"},{"why":"Derives perturbative electronic momentum beyond Born-Oppenheimer, establishing the baseline that BO electronic states have zero momentum.","marker":"[26]"}],"fun_headline_variants":["Phase-space theory gives conical intersections a third branch","Electrons gain momentum at conical intersections in phase space","Complex Hartree-Fock gets physical meaning via phase-space CI","BeH2 crossing reveals 3D branching and electronic angular momentum","Phase-space view: conical intersections carry electronic momentum"],"cache_read_input_tokens":18176,"weakest_assumption_plain":"The phase-space electronic Hamiltonian, specifically the form of the Γ operator that couples nuclear momentum to electronic density, is assumed to be a correct physical model; if that operator is wrong, the three-dimensional branching plane and the ±P_min double well need not exist.","fun_headline_variants_meta":{"raw":{"variants":["Phase-space theory gives conical intersections a third branch","Electrons gain momentum at conical intersections in phase space","Complex Hartree-Fock gets physical meaning via phase-space CI","BeH2 crossing reveals 3D branching and electronic angular momentum","Phase-space view: conical intersections carry electronic momentum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000164,"raw_usage":{"total_tokens":1283,"prompt_tokens":1017,"completion_tokens":266,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":633,"completion_tokens_details":{"reasoning_tokens":187}},"tokens_in":633,"tokens_out":266,"duration_ms":3134,"temperature":1.0,"reasoning_tokens":187,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:59:32.799537+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute exact beyond-Born-Oppenheimer electronic momentum and current for a wavepacket traversing the BeH2 conical intersection at λ = 2.600 Bohr; if the emerging electronic state shows no net angular momentum and a momentum scan of the phase-space surface shows no stable double well at ±P_min, the central claim fails. A direct measurement showing Im(H12) = 0 throughout the seam would likewise contradict the third branching direction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the original phase-space electronic Hamiltonian with derivative couplings, whose divergence near crossings motivates the use of the non-divergent Γ operator."},{"cited_title":"Tao , author T","cited_arxiv_id":null,"evidence_quote":"Gives the basis-free phase-space Hamiltonian that recovers electronic momentum and current density, validating the Γ form."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the BeH2 insertion model system with the conical intersection studied here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives perturbative electronic momentum beyond Born-Oppenheimer, establishing the baseline that BO electronic states have zero momentum."}],"review_version":1}