{"id":"569d2824-9ac1-45ff-9a90-d950359a1cf7","arxiv_id":"2411.13879","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A phase-space electronic Hamiltonian with GIAOs is implemented in Q-Chem, yielding gauge-, translation-, and rotation-invariant energies and predicting nonzero ground-state electronic momentum at finite magnetic fields.","lead":"This paper implements a phase-space electronic Hamiltonian for molecules in magnetic fields using gauge-invariant atomic orbitals (GIAOs) in the Q-Chem package. It shows that, at finite magnetic fields, the minimum energy occurs at nonzero nuclear kinetic momentum, implying electronic motion in the ground state.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central prediction Pi_eff_min != 0 rests on a dropped quadratic term and an untested Gamma ansatz; it needs an exact-limit or independent check.","rationale":"The reader identified the Gamma ansatz as the weakest assumption, and I agree that essentially the concern is that the central prediction depends on an unvalidated, non-unique ansatz. I partially disagree with the framing: the more precise, load-bearing issue is that the prediction Pi_eff_min != 0 is obtained only after dropping a negative-definite quadratic term in Eq. 60, so even within the paper's own model the location of the minimum is not fixed. The paper's conservation-law and GIAO-invariance proofs appear internally consistent and are a genuine strength, but they do not test the physical accuracy of Gamma. The numerical section is explicitly preliminary and the paper itself flags the Gamma^2 truncation and basis incompleteness, so the CONDITIONAL verdict is appropriate, with conditions being inclusion/estimation of the Gamma^2 term and validation against an exact limit or independent calculation.","tokens_in":24569,"tokens_out":2152,"duration_ms":21756,"concrete_test":"Recompute Fig. 2(d) at one finite field (e.g., B_y = 0.1 a.u.) with the previously dropped (Gamma)^2 term in Eq. 60 included. If Pi_eff_min shifts by a significant fraction of its reported value or returns to zero, the headline prediction is an artifact of truncation. In parallel, run the same phase-space GIAO calculation for a hydrogen atom and compare the predicted Pi_eff_min and electronic momentum against the known exact solution; a mismatch would show the Gamma ansatz is not quantitatively validated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that ground-state H2O2 in a static magnetic field has minimum phase-space energy at Pi_eff_min != 0, implying ground-state electronic current. That claim is computed from the truncated Eq. 60, which drops the -hbar^2 (Gamma^2) term. The retained Gamma-linear term is what makes the surface shift; the dropped term is negative and grows with Gamma, so the true minimum of the full E_PS could differ substantially. The Gamma operator itself is a phenomenological ansatz from Paper I, constrained by conservation laws (Eqs. 61-64) but not uniquely determined, and its weight functions theta_I use electronegativity-fitted sigma_I values. No calculation validates the resulting Pi_eff_min against an exact limit (e.g., hydrogen atom, where Paper I claims exactness) or an independent treatment. The symmetry proofs are consistency checks, not accuracy checks. Thus the quantitative prediction Pi_eff_min != 0, and even its sign or existence, is not yet robust.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript implements a phase-space electronic Hamiltonian for molecules in a static magnetic field in a finite atomic-orbital basis, using gauge-including atomic orbitals (GIAOs). It derives explicit GIAO matrix elements for the electron translation/rotation operator Γ, proves gauge-origin independence and translational/rotational invariance of the phase-space energy, and proves conservation of total pseudomomentum and canonical angular momentum in the GIAO basis. Preliminary Q-Chem calculations on H2O2 show that, for magnetic fields along y or z, the phase-space energy is minimized at a nonzero effective nuclear momentum Πeff_min, which the authors interpret as evidence of nonzero electronic motion in the ground state. The paper is a companion to a theoretical paper (Paper I) and focuses on the practical GIAO formulation and a first numerical demonstration.","tokens_in":24864,"tokens_out":10454,"duration_ms":119320,"significance":"If correct, the paper provides a practical route to beyond-Born-Oppenheimer electronic structure in magnetic fields: the new Γ terms are one-electron, the GIAO formulation removes gauge-origin artifacts, and the conservation laws are proven in a finite basis. The explicit working equations and the implementation in a standard package are valuable, and the prediction of a momentum-dependent potential surface with Πeff_min ≠ 0 is falsifiable. The main weaknesses are that the central numerical prediction is controlled by an ansatz for Γ that is constrained but not uniquely determined or independently validated, and that the momentum observables used for the physical interpretation are shown to suffer from basis-set artifacts. The specific concern that the dropped Γ² term shifts Πeff_min does not survive inspection, because that term is independent of Πeff and therefore does not change the stationarity condition; it only offsets the absolute energy.","major_comments":[{"comment":"The Γ operator is introduced as a 'meaningful form' satisfying the conservation constraints (61)-(64), but it is not derived uniquely, and its weight functions Θ_I depend on σ_I fitted to electronegativity. Since the minimum condition Πeff_min = iℏ⟨Γ⟩ and the electronic-momentum expressions in Eqs. (94)-(95) are linear in Γ, the central numerical prediction is entirely controlled by this ansatz. The paper does not validate Γ against an exact limit (e.g., the hydrogen atom, for which Paper I claims exactness) or against an independent finite-basis calculation. Please add such a benchmark or at least a sensitivity analysis with respect to the σ_I parameters.","section":"Sec. 3, Eqs. (52)-(57)"},{"comment":"Translational invariance of E_PS is asserted for Γ'' with the statement that the phase factor arises from transforming J and that K_J and (X_I - X0_J) are translationally invariant. This is a load-bearing identity for the pseudomomentum conservation proof in Sec. 4.2 and Appendix B. The paper should supply the derivation of Eq. (63) for both Γ' and Γ'', or an explicit appendix lemma, rather than leaving the central step as 'straightforward'.","section":"Sec. 3.2 and Eqs. (61)-(64)"},{"comment":"The paper acknowledges that the computed ⟨π_e⟩ at Πeff = 0 deviates from the exact zero due to an incomplete basis, and that GIAOs do not cure this. Because the central interpretation ('non-zero electronic motion in the ground state') relies on momentum observables, the paper should demonstrate basis-set convergence of ⟨π_e⟩ and ⟨r × π_e⟩ (e.g., cc-pVTZ or cc-pVQZ) or otherwise quantify the artifact at the field strengths where Πeff_min is reported in Fig. 2(d). Without this, the quantitative momentum signal supporting the main claim is not yet established.","section":"Sec. 5, Fig. 3"},{"comment":"The Γ² term is dropped as 'small' without quantification. Since Γ is independent of Πeff, dropping it does not change the stationarity condition for Πeff_min, so the core prediction of Fig. 2(d) is not affected by this truncation as far as the Πeff scan is concerned. However, the term does change absolute energies and would affect any geometry optimization; the paper should state this distinction explicitly and report the magnitude of the omitted term, especially for the energy curve in Fig. 2(e) and for the phrase 'minimum energy structures' in the abstract.","section":"Eq. (60) and Sec. 5"}],"minor_comments":[{"comment":"The symbol π appears to be used for what should be the two-electron integral g in the rotational-invariance identity; please correct the notation.","section":"Sec. 2.4, Eq. (45)"},{"comment":"The vertical shifts of the energy curves in Figs. 2(a)-(c) make it impossible to judge the magnitude of the energy variation with Πeff; please report at least one unshifted energy scale or the relevant energy differences.","section":"Sec. 5, Fig. 2"},{"comment":"The text states values 1/σ² = 2.42 and 3.79 for O and H without specifying units or the fitting procedure; a reference or a short description of the electronegativity parametrization would improve reproducibility.","section":"Sec. 5, Eq. (48)"},{"comment":"The phrase 'minimum energy structures' is stronger than what is computed, since only a minimum with respect to Πeff at a fixed nuclear geometry is demonstrated; please clarify that geometry optimization is not performed here.","section":"Abstract and Sec. 6"},{"comment":"The convention for factors of 1/(iℏ) in the relation between ⟨π_e⟩ and ⟨Γ⟩ is easy to misread; a short dimensional or notational comment would help.","section":"Sec. 5 near Eq. (94)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically interesting and the overall program is promising. My recommendation is driven by the need to validate the Γ ansatz and the momentum observables, not by any detected logical error in the symmetry derivations. The Γ² truncation, which might appear to undermine Πeff_min, is indeed not a problem for the stationarity condition because Γ is P-independent; this should be clarified to the authors so they do not overcorrect. The main risk is that the quantitative prediction in Fig. 2(d) is model-dependent until benchmarked or subjected to a parameter-sensitivity study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere is my read on Bhati et al. (2411.13879). The paper's real contribution is the GIAO machinery: explicit matrix elements for the Gamma' and Gamma'' operators, plus proofs that the phase-space energy is translationally and rotationally invariant in a finite basis and that pseudomomentum and angular momentum are conserved. Those derivations look careful and self-consistent. The authors are also honest about known blemishes: they flag the basis-set error in the linear momentum at Pi_eff = 0 (Fig 3a) and explicitly note they dropped the -hbar^2 Gamma^2 term in Eq 60.\n\nThe soft spot is the central physical prediction. The claim that ground-state H2O2 in a magnetic field has a minimum at Pi_eff_min != 0, implying circulating electronic current, is computed with Eq 60 truncated by that dropped term. The retained Gamma-linear term is what produces the surface shift; the dropped term is negative and grows with Gamma, so the true location of the minimum is not established until that term is quantified. The Gamma operator itself is a phenomenological ansatz from Paper I, constructed to satisfy conservation constraints rather than derived; the sigma_I parameters come from electronegativity fits. The symmetry proofs are consistency checks, not accuracy checks. The stress-test note got the right chord: no exact-limit or independent validation of the Pi_eff_min prediction appears in the paper.\n\nThat said, I do not think this is a fatal flaw. This is a technical methods paper, and the preliminary H2O2 numbers are labeled as such. The derivations in Secs 2-4 are the payload, and they look solid. The reader's \"conditional\" verdict is about right: accept if the authors can (1) bound or include the Gamma^2 term, (2) test sensitivity to the Gamma ansatz and sigma_I, and (3) run a check where an exact limit is known, e.g., the hydrogen atom from Paper I. The paper deserves a serious referee; it is not a desk reject. I would send it out with those conditions in mind. I would bring it to a reading group as a useful example of how to get symmetry properties in a finite basis, but I would not build on the Pi_eff_min result until it is confirmed.","headline":"The GIAO implementation and symmetry proofs are the real contribution; the headline Pi_eff_min != 0 prediction is suggestive but rests on a dropped term and an unvalidated ansatz, so treat the numerics as preliminary.","tokens_in":25334,"tokens_out":2828,"would_cite":true,"duration_ms":28686,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.15.ae","31.15.-p","33.15.-e"],"model":"deepseek-v4-flash","headline":"For a molecule in a static magnetic field, computing with gauge-invariant atomic orbitals makes the phase-space electronic Hamiltonian translationally and rotationally invariant, and the ground state can have nonzero electronic motion.","keywords":["phase-space electronic Hamiltonian","gauge-including atomic orbitals","magnetic field","electron translation factors","electron rotational factors","pseudomomentum conservation","angular momentum conservation","beyond Born-Oppenheimer"],"falsifier":"For the hydrogen atom in a uniform magnetic field, compute the phase-space Hamiltonian's ground-state energy and electronic current density in a large GIAO basis and compare with the exact analytical solution; the paper claims exactness for hydrogen, so any discrepancy would invalidate the $\\Gamma$ ansatz or its GIAO implementation.","tokens_in":24335,"feed_emoji":"🧲","tokens_out":10331,"duration_ms":90224,"temperature":0.7,"pith_summary":"This paper puts the phase-space electronic Hamiltonian from Paper I into a practical atomic-orbital form by dressing the basis with gauge-including atomic orbitals (GIAOs). It proves that, even in a truncated basis, the resulting energies are independent of the magnetic gauge origin and invariant under translation and rotation. It then shows that Hamilton's equations along an eigenstate conserve total pseudomomentum and total canonical angular momentum in the field direction. The central physical result is that, for a molecule like H2O2 in a finite magnetic field, the phase-space potential energy is minimized at a nonzero effective nuclear momentum, which the authors interpret as nonzero electronic motion in the ground state even when the nuclei are stationary.","feed_headline":"Magnetic fields set ground-state electrons in motion","feed_subtitle":"A phase-space Hamiltonian with gauge-invariant orbitals predicts nonzero electronic momentum even for stationary nuclei.","key_machinery":"The central object is the phase-space electronic Hamiltonian $\\hat{H}_{PS}(\\mathbf{X},\\mathbf{\\Pi}_{\\mathrm{eff}}) = \\sum_I (\\mathbf{\\Pi}_{\\mathrm{eff},I})^2/(2M_I) + V_{PS}(\\mathbf{X},\\mathbf{\\Pi}_{\\mathrm{eff}},G,\\mathbf{B})$, whose eigenstates are parameterized by both nuclear positions and an effective nuclear momentum. The coupling between nuclear momentum and electronic motion enters through the operator $\\hat{\\Gamma}_I$, a sum of electron translation and rotation factors built from the electronic pseudomomentum $\\hat{k}^I$ relative to each nucleus and an electronegativity-weighted projector $\\hat{\\Theta}_I$. Gauge-including atomic orbitals, which multiply each atomic function $\\phi_\\mu$ by the phase $\\exp(-ie \\mathbf{A}_M\\cdot\\mathbf{r})/\\hbar$, remove the dependence on the gauge origin $G$ and make all matrix elements (including those of $\\hat{\\Gamma}_I$) well defined in a finite basis; the conservation proofs rest on how these phases transform under translation and rotation.","core_discovery":"The paper establishes that the phase-space electronic Hamiltonian $\\hat{H}_{PS}(\\mathbf{X},\\mathbf{\\Pi}_{\\mathrm{eff}})$, evaluated in a GIAO basis, has all one- and two-electron integrals independent of the gauge origin $G$, and that the resulting energy $E_{PS}$ is translationally and rotationally invariant. From these invariances and the explicit matrix elements of the electron translation and rotation factors $\\hat{\\Gamma}_I$, the authors prove that the total pseudomomentum $\\mathbf{K}_{\\mathrm{mol}}$ and the $z$-component of the total canonical angular momentum $L^z_{\\mathrm{mol}}$ are conserved during dynamics. For $\\mathrm{H}_2\\mathrm{O}_2$, the phase-space potential surface at nonzero magnetic field has its minimum at $\\mathbf{\\Pi}_{\\mathrm{eff}}^{\\mathrm{min}} \\neq 0$; because the true kinetic momentum vanishes at that point, the nonzero $\\mathbf{\\Pi}_{\\mathrm{eff}}$ is not nuclear motion but rather a signature of circulating electronic current in the ground state.","pith_inferences":["If the nonzero $\\mathbf{\\Pi}_{\\mathrm{eff}}^{\\mathrm{min}}$ persists in large basis sets and with the $\\Gamma^2$ term included, standard Born-Oppenheimer magnetic-field calculations are missing a genuine electronic polarization that affects structure and dynamics.","The same GIAO phase-space machinery could be adapted to compute magnetic-field-induced currents and electronic momentum densities directly, giving a practical route to visualize electrons circulating in a molecular ground state.","A direct numerical test of the conservation laws (e.g., a trajectory of a diatomic in a uniform field) would either confirm the proofs or reveal a truncation artifact in the $\\Gamma$ ansatz; no such test is reported here.","The electron rotation factors in $\\Gamma''$ are constructed from the geometry of nuclear positions; applying this to nonrigid or floppy molecules may require a more flexible definition."],"forward_implications":["Quantum-chemistry packages can compute magnetic-field energy surfaces with the phase-space Hamiltonian using only one-electron additions to a standard Born-Oppenheimer code.","Molecular dynamics in magnetic fields can be run on an eigenstate of this Hamiltonian without tracking Berry forces, and the conserved pseudomomentum and angular momentum will be preserved.","Geometry optimizations in a magnetic field should allow the effective nuclear momentum to relax to a nonzero value, revealing electronic currents that Born-Oppenheimer calculations miss.","Spectroscopic observables that depend on electronic angular momentum, such as magnetic circular dichroism, can be connected to the $\\mathbf{\\Pi}_{\\mathrm{eff}}$-dependent terms in the phase-space potential."],"supporting_citations":[{"why":"The companion Paper I that defines the phase-space Hamiltonian and its basis-free conservation laws; this paper supplies the formalism being implemented.","marker":"[24]"},{"why":"The standard GIAO method (magnetic-field-dependent atomic orbital phases) that the paper adapts to the phase-space operator.","marker":"[30]"},{"why":"The origin-independent one-electron Hamiltonian construction with GIAOs that underlies the matrix-element proofs in Secs. 2 and 3.","marker":"[35]"},{"why":"Shows the field-free phase-space Hamiltonian recovers vibrational circular dichroism and electronic momenta, motivating the $\\Gamma$ terms.","marker":"[11]"},{"why":"A basis-free phase-space Hamiltonian that recovers beyond-Born-Oppenheimer electronic momentum and current densities; the basis for the $\\Gamma$ ansatz.","marker":"[29]"},{"why":"Provides the one-electron expression for electron rotational factors used in the $\\Gamma''$ term.","marker":"[27]"},{"why":"Textbook Wigner rule used to evaluate nuclear gradients of the density matrices without solving orbital response.","marker":"[41]"},{"why":"The electronic structure package in which the GIAO phase-space implementation is built.","marker":"[42]"}],"fun_headline_variants":["Ground-state electrons circulate in magnetic fields","Nonzero electronic momentum in ground state","Ground-state electrons move even with stationary nuclei","Magnetic fields drive ground-state electron motion","Phase-space Hamiltonian reveals circulating currents"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the electron translation and rotation factors, whose form is fixed by symmetry constraints and electronegativity-based parameters, correctly describe how nuclear momentum and a magnetic field affect the electrons; no exact or experimental benchmark is provided for this ansatz.","fun_headline_variants_meta":{"raw":{"variants":["Ground-state electrons circulate in magnetic fields","Nonzero electronic momentum in ground state","Ground-state electrons move even with stationary nuclei","Magnetic fields drive ground-state electron motion","Phase-space Hamiltonian reveals circulating currents"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000498,"raw_usage":{"total_tokens":2454,"prompt_tokens":977,"completion_tokens":1477,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":1415}},"tokens_in":593,"tokens_out":1477,"duration_ms":10621,"temperature":1.0,"reasoning_tokens":1415,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:45:50.358486+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the hydrogen atom in a uniform magnetic field, compute the phase-space Hamiltonian's ground-state energy and electronic current density in a large GIAO basis and compare with the exact analytical solution; the paper claims exactness for hydrogen, so any discrepancy would invalidate the $\\Gamma$ ansatz or its GIAO implementation.","supporting_citations":[{"cited_title":"A Phase-Space Electronic Hamiltonian for Molecules in a Static Magnetic Field I: Conservation of Total Pseudomomentum and Angular Momentum","cited_arxiv_id":"2411.13866","evidence_quote":"The companion Paper I that defines the phase-space Hamiltonian and its basis-free conservation laws; this paper supplies the formalism being implemented."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The standard GIAO method (magnetic-field-dependent atomic orbital phases) that the paper adapts to the phase-space operator."},{"cited_title":"An electronic Hamiltonian for origin independent calculations of magnetic properties","cited_arxiv_id":null,"evidence_quote":"The origin-independent one-electron Hamiltonian construction with GIAOs that underlies the matrix-element proofs in Secs. 2 and 3."},{"cited_title":"G.; Pei, Z.; Shao, Y.; Subotnik, J","cited_arxiv_id":null,"evidence_quote":"Shows the field-free phase-space Hamiltonian recovers vibrational circular dichroism and electronic momenta, motivating the $\\Gamma$ terms."},{"cited_title":"G.; Subotnik, J","cited_arxiv_id":null,"evidence_quote":"Provides the one-electron expression for electron rotational factors used in the $\\Gamma''$ term."}],"review_version":1}