{"id":"dea4e92f-2328-4d1e-b9f8-922fa906d2ea","arxiv_id":"2504.17107","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"Quantum-chemistry torsional energy scans of the ammonium dimer show non-smooth, history-dependent profiles that the author labels electronic energy singularities, but the analysis does not rule out SCF convergence artifacts.","lead":"The paper reports that quantum-chemistry energy scans of the ammonium dimer show discontinuities and noise under torsion, while molecular mechanics gives smooth curves. It interprets this as a quantum signature of weak hydrogen bond breaking, but the evidence points to numerical artifacts rather than physical singularities.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The evidence for physical electronic-energy singularities is indistinguishable from SCF convergence failure, and the paper provides no convergence tests to separate them.","rationale":"The reader's weakest assumption—that the discontinuities might be numerical SCF artifacts rather than physical singularities—is exactly the load-bearing vulnerability. The paper includes no convergence diagnostics, no derivation of why weak H-bonds should produce singular electronic energies, and no state-resolved analysis. The reported hysteresis and multiple dipole levels are particularly telling because they indicate different SCF solutions being selected, not a single well-converged surface. The concrete test is decisive and inexpensive: tighten SCF convergence and vary protocol/initial guesses. If the singularities persist under all protocols, the paper's central claim gains credibility; if not, it is an artifact. I therefore see no reason to change the reader's REJECT verdict, and I agree with the identification of this as the weakest assumption. Secondary problems (invalid MMFF control, missing methodological detail, reliance on self-cited prior work) reinforce but do not replace this concern. No exceptional circumstances such as machine-checked proofs or reproducible code are present to offset the evidentiary gap.","tokens_in":4655,"tokens_out":3114,"duration_ms":30891,"concrete_test":"Repeat one representative B3LYP/6-311G* forward and reverse torsion scan (same 1° resolution, same constrained N1-H5 distances from Table 1) at an SCF convergence threshold of 10^-12 instead of the default, with and without DIIS acceleration, and from two initial guesses: (i) a converged wavefunction from the neighboring torsion angle and (ii) a fresh guess from the optimized dimer at each angle. If the energy and dipole profiles remain identical at the same torsion grids—same jump locations, same branch values within 10^-10 Hartree—the singularity interpretation survives this test. If the jumps shift, vanish, or depend on the initial guess, the observed features are SCF artifacts and the central claim is falsified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the non-smooth, discontinuous, hysteretic QM torsion profiles are physical electronic-energy singularities inherent to the weak H-bond. The paper asserts this in the Introduction ('inherent in the system and not due to any numerical convergence issues') but offers no supporting evidence. All reported signatures—noise, isolated discontinuities, forward/reverse hysteresis, and multiple discrete dipole levels in Fig. 5—are the standard fingerprints of SCF convergence failure or of a scan hopping between different self-consistent solutions. No SCF convergence thresholds, initial-guess strategies, DIIS/damping settings, or comparisons between algorithms are reported. The MMFF/SYBYL control is also not a clean null: Table 1 shows MMFF's N1-H5 distance is 1.1 Å rather than an H-bond distance, so the QM/MM contrast may reflect different conformers rather than a quantum singularity. A physical singularity in a Born-Oppenheimer potential-energy surface as a function of a torsion angle is not a generic consequence of weak hydrogen bonding; such surfaces are smooth except at genuine degeneracies, and no conical intersection or state-crossing analysis is given. The Mathieu-equation paragraph in the Conclusion is an analogy about torsional eigenstates, not a derivation that electronic energy must be singular. If the QM jumps disappear when SCF convergence is tightened or when the initial guess is varied, the headline claim collapses.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports torsion-angle scans of the ammonium dimer computed with molecular mechanics (MMFF, SYBYL) and with quantum-chemical methods (HF, B3LYP, MP2). It finds that MM methods give smooth, continuous dihedral energy profiles, while QM methods give noisy, discontinuous profiles with apparent jumps and hysteresis. The author interprets these non-smooth QM profiles as physical electronic-energy singularities intrinsic to the weak N1-H5 hydrogen bond, and claims that the effect is persistent across all standard ab initio techniques. Additional torsion-dependent dipole calculations are presented as showing discrete dipole levels in QM but not MM, and a Mathieu-equation analogy is offered in the conclusion.","tokens_in":4858,"tokens_out":2867,"duration_ms":29299,"significance":"If the central claim were correct, it would be highly significant: the existence of genuine electronic-energy singularities in a closed-shell molecular potential-energy surface as a function of a torsion angle would challenge standard Born-Oppenheimer surface theory and would imply a new quantum signature of weak hydrogen bonds. The manuscript does use several electronic-structure methods and reports quantitative bond lengths and binding energies, which is a reasonable starting point. However, the paper provides no SCF convergence tests, no error analysis, no basis-set superposition or level-shift diagnostics, and no independent check distinguishing a physical singularity from a numerical artifact. The evidence presented is fully consistent with well-known SCF convergence failures or with scans hopping between different self-consistent solutions. As presented, the claim is not supported, and the paper's interpretative framework is circular in that the 'quantum signature' is inferred from the very QM outputs whose reliability is at issue.","major_comments":[{"comment":"The central assertion that the observed discontinuities are 'inherent in the system and not due to any numerical convergence issues' is never tested. The manuscript reports no SCF convergence thresholds, no initial-guess strategies, no DIIS or damping settings, no comparison of different SCF algorithms, and no check that the noisy energy profiles are converged with respect to tightening criteria. Noise, isolated discontinuities, forward/reverse hysteresis, and multiple discrete dipole levels are classic fingerprints of SCF convergence failure or of a scan switching between different local solutions. Because the paper offers no diagnostic to exclude these mundane explanations, the load-bearing distinction between a physical singularity and a numerical artifact is not established.","section":"Introduction and Computational Methodologies"},{"comment":"The MMFF and SYBYL controls are not a valid null model. Table 1 lists the MMFF N1-H5 distance as 1.1 Å and the SYBYL distance as 1.6 Å, whereas the QM methods give 2.2-2.4 Å. An N1-H5 distance of 1.1 Å is not a hydrogen bond; it is a chemically different structure, essentially a merged or rearranged geometry. Consequently, the smooth MM torsion profiles may simply reflect a different conformer or a different potential-energy region, not the absence of a quantum singularity in the hydrogen-bonded system. The contrast in Figs. 4 and 6 therefore cannot support the paper's interpretation without first establishing that the MM and QM scans sample the same physical region.","section":"Table 1, Figs. 4A/4B, 6A/6B"},{"comment":"The Mathieu-equation paragraph is an analogy about torsional eigenstates of a periodic potential, not a derivation that an electronic Born-Oppenheimer energy surface must be singular. For a closed-shell system with no symmetry-enforced degeneracy, electronic energy as a function of a torsion angle is generically smooth. The manuscript presents no conical-intersection analysis, no state-crossing search, and no argument explaining why weak hydrogen bonding would create a true electronic singularity. Therefore this paragraph does not provide theoretical support for the headline claim.","section":"Conclusion"},{"comment":"The interpretation of multiple discrete dipole-moment levels as a quantum signature is circular: the same QM outputs that are asserted to be singular are used as evidence for the singular nature. Multiple discrete dipole values with jumps are precisely what is observed when an SCF calculation converges to different local minima or to different charge-localized solutions along a scan. Without reporting the dipole values at each step, the convergence history, or a comparison with coupled-cluster or other robust methods, the dipole data cannot distinguish a physical quantization from a numerical artifact.","section":"Fig. 5 and Results and Discussion"}],"minor_comments":[{"comment":"The word 'discreet' appears where 'discrete' is meant; this occurs in the abstract and in the figure captions.","section":"Abstract and Fig. 5"},{"comment":"The reference list heading is misspelled as 'REFRENCES', and reference [8] contains a trailing 'Szalewicz' that appears to be an artifact of citation formatting.","section":"References"},{"comment":"In the results text, 'B3YLP' should be 'B3LYP'.","section":"Results and Discussion"},{"comment":"The basis-set labels are inconsistent ('HF/6-31G+', 'HF/6-311G+', 'B3LYP/6-311G+', 'MP2/aug-cc-PVDz') and do not specify polarization/diffuse functions consistently; '6-31G+' and '6-311G+' are not standard unambiguous basis-set names.","section":"Computational Methodologies and Table 1"},{"comment":"The scan setup is incompletely specified: the atoms defining the torsion angle are not named in the text, the convention for forward versus reverse scans is not defined, and the constrained-distance minimization protocol is described only qualitatively.","section":"Computational Methodologies"},{"comment":"The text refers to 'dissociation energy' estimated from well depth, while Table 1 labels the quantity 'Dimer Binding Energy', without specifying whether counterpoise correction or zero-point energy is included; the 13.0 kJ/mol value appears to conflate binding energy with the depth of a constrained one-dimensional scan.","section":"Table 1 and text"}],"recommendation":"reject","confidential_remarks":"The paper is a short report whose central claim is unsupported by the evidence it provides. The missing convergence tests are not a cosmetic deficiency; they are the only way to test whether the observed non-smooth QM profiles are physical. Given that the reported signatures are all standard SCF-convergence-failure indicators, and the MM control structures are not chemically comparable to the QM structures, the manuscript's conclusion is likely to be an artifact. Even with additional calculations, the paper would need major new content (e.g., coupled-cluster scans, variational stability checks, state symmetry analysis) to justify reconsideration. I would not encourage resubmission in its current form, though a much more thorough study of the underlying numerics could be a separate worthwhile project."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper makes a strong claim: that the Born-Oppenheimer electronic energy of the weakly H-bonded ammonium dimer is singular as a function of torsion. The evidence, however, is all numerical noise consistent with SCF convergence failure. The author never runs a single convergence test.\n\nWhat's actually new: the torsional scans for ammonium dimer at HF/B3LYP/MP2 are new data, and the raw data is on Zenodo. The observation that QM scans are noisy while MM scans are smooth is reproducible in principle. But that is a computational observation, not a discovery.\n\nThe soft spots are load-bearing. The claim that the singularities are physical is asserted in the Introduction ('inherent in the system and not due to any numerical convergence issues') and never tested. No SCF thresholds, initial guesses, DIIS settings, or algorithm comparisons are reported. The forward/reverse hysteresis and multiple dipole levels in Fig. 5 are textbook signs of SCF convergence failure. The MMFF control is invalid: Table 1 lists N1-H5 = 1.1 Å, not a hydrogen-bond distance, so the QM/MM contrast compares different structures. The Mathieu-equation paragraph in the Conclusion is an analogy, not a derivation. The only conceptual support is the author's own prior work (refs 9, 10), making the interpretation circular.\n\nTo the author's credit, the data repository is provided and the figures show the raw profiles. The paper is honest in reporting repeated and reverse scans, which is more than some papers do. But the interpretation is not supported.\n\nThis paper is for a reader studying SCF convergence artifacts in weak dimers—the data is a nice illustration of what can go wrong. As a claim about physical singularities, it is not worth referee time unless the author returns with convergence tests. Desk reject is appropriate.","headline":"The paper's central claim—physical electronic-energy singularities in the ammonium dimer—is indistinguishable from SCF convergence failure, and the author offers no tests to separate them.","tokens_in":5431,"tokens_out":3300,"would_cite":false,"duration_ms":29831,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that the weak hydrogen bond in the ammonium dimer produces genuine electronic energy singularities under torsion, visible as jumps and noise in HF, B3LYP, and MP2 profiles, rather than numerical artifacts.","keywords":["ammonium dimer","hydrogen bond","electronic energy singularity","torsion angle scan","ab initio methods","molecular mechanics","dipole moment discontinuity","quantum signature"],"falsifier":"Record the same 360-degree torsion scan for the dimer while progressively tightening the electronic-structure solver's convergence thresholds and switching initial guess orbitals; if the jumps move, vanish, or coincide with convergence warnings, they are numerical artifacts rather than electronic singularities. A multi-reference calculation along the same path that returns a smooth profile would also falsify the claim.","tokens_in":4382,"feed_emoji":"⚛️","tokens_out":6947,"duration_ms":62809,"temperature":0.7,"pith_summary":"This paper tries to establish that the weakly hydrogen-bonded ammonium dimer carries genuine electronic energy singularities: under a 360-degree torsion of the N1-H5 bond, ab initio methods (HF, B3LYP, MP2) produce noisy, discontinuous energy and dipole-moment profiles, while molecular mechanics (MMFF, SYBYL) produces smooth curves. The author argues these non-smooth features are inherent quantum-mechanical consequences of weak hydrogen-bond breaking under torsion, not numerical convergence artifacts. If true, the result would mean standard quantum-chemistry scans expose bond-breaking criticality that force-field models smooth over, and that other weak H-bonded and van der Waals dimers should show similar anisotropic singularities.","feed_headline":"Quantum chemistry finds energy jumps in weak H-bonded dimer","feed_subtitle":"If real, the jumps mark bond-breaking criticality that smooth force-field models miss.","key_machinery":"The central probe is a 1-degree-resolution torsion scan around the weak N1-H5 hydrogen bond of the ammonium dimer, performed with constrained N1-H5 distances in the 2.3 to 2.8 angstrom range and compared between ab initio methods (HF/6-311G*, B3LYP/6-311G*, MP2/aug-cc-pVDZ) and molecular mechanics (MMFF94, SYBYL). The comparison of energy and dipole-moment profiles is the mechanism that carries the argument: the repeated contrast of non-smooth quantum profiles against smooth molecular-mechanics profiles is offered as evidence that the singularities are electronic in origin. The conclusion additionally frames quantized torsional energy as the natural solution of a torsion-potential Schrödinger equation, via the Mathieu equation, the standard differential equation for a particle in a periodic potential.","core_discovery":"The paper's central discovery is that, for the ammonium dimer held together by a weak N1-H5 hydrogen bond, every ab initio method tested here (HF, B3LYP, MP2) yields torsion-dependent electronic energy profiles that are not smooth: they contain noise, jumps, and discontinuities, whereas the molecular-mechanics methods MMFF and SYBYL give smooth continuous profiles for the same rotation. The same contrast appears in computed dipole moments, with quantum methods showing several discrete levels and jumps while molecular mechanics shows two smooth levels. The author asserts this non-smoothness is an inherent electronic energy singularity of weak hydrogen-bond breaking under torsion, not a convergence artifact, and that repeated forward and reverse scans reproduce the feature.","pith_inferences":["If the claim holds, the same torsion-scan approach could be applied to other weak dimers such as water, HF, or mixed rare-gas-ammonia clusters to map where singularities appear and whether their angular position correlates with hydrogen-bond strength.","A sharper test than the paper performs is to check whether each energy jump coincides with an avoided crossing or a change in the character of the highest-occupied molecular orbital; that would connect the singularities to a concrete electronic-structure mechanism.","The claim implies that force-field-based sampling over torsion angles in weakly bonded clusters visits configurations that correspond to no stable electronic state, which would affect how conformational searches are interpreted."],"forward_implications":["If the singularities are real, molecular-mechanics torsion profiles misrepresent the electronic energy of weakly H-bonded dimers, so force-field scans should not be used to judge which conformers are energetically accessible.","The persistence of the jumps across repeated forward and reverse scans inside each method means the feature is reproducible, not a one-time convergence failure within a single run.","The same anisotropic singularity behavior is expected in other weak H-bonded and van der Waals dimers with dissociation energies near 10 to 25 kJ/mol.","Quantum-computed dipole moments along the torsion should show discrete stepped levels rather than smooth variation, offering a second observable signature of the same effect."],"supporting_citations":[{"why":"Establishes that the ammonium dimer is weakly hydrogen-bonded and extremely fluxional, the premise of the study.","marker":"[8]"},{"why":"Supplies the earlier finding of a quantum electronic signature at bond break-up that motivates interpreting the jumps as physical.","marker":"[9]"},{"why":"Shows the analogous anisotropic singularity in the water dimer, the pattern this paper extends to ammonium.","marker":"[10]"},{"why":"Benchmarks the ab initio dissociation energy and geometry values against converged basis-set results.","marker":"[12]"},{"why":"Provides the intermolecular potentials used to estimate the roughly 13 kJ/mol hydrogen-bond strength.","marker":"[17]"},{"why":"Gives the protocol for extracting dissociation energy from computed energy scans.","marker":"[18]"}],"fun_headline_variants":["Quantum methods show energy jumps in weak H-bonded dimer","Erratic QM energy profiles mark weak H-bonded ammonium dimer","Ab initio energy discontinuities in ammonium dimer under torsion","QM vs MM: energy singularities in weak H-bonded dimer","Torsion reveals electronic energy singularities in ammonium dimer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the discontinuous, noisy energy profiles from HF, B3LYP, and MP2 are physical electronic features of the weak hydrogen bond rather than failures of the self-consistent-field iteration to converge at certain geometries.","fun_headline_variants_meta":{"raw":{"variants":["Quantum methods show energy jumps in weak H-bonded dimer","Erratic QM energy profiles mark weak H-bonded ammonium dimer","Ab initio energy discontinuities in ammonium dimer under torsion","QM vs MM: energy singularities in weak H-bonded dimer","Torsion reveals electronic energy singularities in ammonium dimer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000214,"raw_usage":{"total_tokens":1369,"prompt_tokens":832,"completion_tokens":537,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":448,"completion_tokens_details":{"reasoning_tokens":452}},"tokens_in":448,"tokens_out":537,"duration_ms":4852,"temperature":1.0,"reasoning_tokens":452,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:48:25.706197+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Record the same 360-degree torsion scan for the dimer while progressively tightening the electronic-structure solver's convergence thresholds and switching initial guess orbitals; if the jumps move, vanish, or coincide with convergence warnings, they are numerical artifacts rather than electronic singularities. A multi-reference calculation along the same path that returns a smooth profile would also falsify the claim.","supporting_citations":[{"cited_title":"The Ammonia Dimer Equilibrium Dissociation Energy: Convergence to the Basis Set Limit at the Correlated Level","cited_arxiv_id":null,"evidence_quote":"Benchmarks the ab initio dissociation energy and geometry values against converged basis-set results."},{"cited_title":"Intermolecular Potentials for Ammonia Based on the Test Particle Model and the Coupled Pair Functional Method","cited_arxiv_id":null,"evidence_quote":"Provides the intermolecular potentials used to estimate the roughly 13 kJ/mol hydrogen-bond strength."},{"cited_title":"Exploring Chemistry with Electronic Structure Methods: A Guide to Using Gaussian","cited_arxiv_id":null,"evidence_quote":"Gives the protocol for extracting dissociation energy from computed energy scans."}],"review_version":1}