Pith. sign in

REVIEW 4 major objections 6 minor 16 references

Electronic Energy Singularities of Weakly H-bonded Ammonium Dimer

T0 review · 4 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2504.17107 v2 pith:FDMGFE5D submitted 2025-04-23 physics.chem-ph

Md Rejwan Ali This is my paper · ORCID
classification physics.chem-ph
keywords ammoniumdimerhydrogenbondelectronicenergysingularitytorsionanglescanabinitiomethodsmolecularmechanicsdipolemomentdiscontinuityquantumsignature
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

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.

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 (4)
  1. [Introduction and Computational Methodologies] 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.
  2. [Table 1, Figs. 4A/4B, 6A/6B] 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.
  3. [Conclusion] 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.
  4. [Fig. 5 and Results and Discussion] 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.
minor comments (6)
  1. [Abstract and Fig. 5] The word 'discreet' appears where 'discrete' is meant; this occurs in the abstract and in the figure captions.
  2. [References] The reference list heading is misspelled as 'REFRENCES', and reference [8] contains a trailing 'Szalewicz' that appears to be an artifact of citation formatting.
  3. [Results and Discussion] In the results text, 'B3YLP' should be 'B3LYP'.
  4. [Computational Methodologies and Table 1] 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.
  5. [Computational Methodologies] 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.
  6. [Table 1 and text] 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.

Circularity Check

2 steps flagged · score 6.0 of 10

Central singularity claim rests on the author's own prior papers and on relabeling noisy QM scans as physical singularities without convergence diagnostics.

  1. self citation load bearing [Introduction, first paragraph (refs [9,10])]
    "As per recent reported works, such weak and unstable H-bonded dimers’ pseudo-stable conformers, mostly generated via molecular mechanics method with a continuum of energy profile, are highly unlikely to have geometrical stability with torsion and must carry quantum break-point signature in electronic energy computation due to the intrinsic weak H-bond [9,10]."

    This sentence is the sole justification for the paper's central premise: a weak H-bond must carry a quantum 'break-point signature' in the electronic energy. Refs [9] and [10] are the author's own prior works (Ali & Mezei 2021; Ali 2022) that make the same interpretive claim for rivastigmine and water dimer. They are not independent derivations or machine-checked results. The paper then uses the noisy QM profiles as confirmation of that premise, so the conclusion 'weak H-bond ⇒ electronic-energy singularity' reduces to the earlier self-citation rather than to new, externally validated evidence.

  2. self definitional [Introduction, first paragraph; cf. Results and Discussion opening paragraph]
    "and such energy results are inherent in the system and not due to any numerical convergence issues in electronic energy computation [11]."

    By 'singularity feature' the paper means the noisy, discontinuous, hysteretic QM torsion energy and dipole profiles (Figs 2, 3, 5, 6). The only statement that separates the physical-singularity reading from a numerical-artifact reading is this sentence, which asserts that the results are 'inherent in the system' and cites a microwave-spectroscopy textbook [11] rather than any convergence study. The paper reports no SCF convergence-threshold sweeps, initial-guess variations, DIIS/damping changes, or algorithm comparisons. Therefore the central claim—weak H-bonding produces an electronic-energy singularity—is equivalent by construction to the observation that the QM scans are non-smooth, with the alternative explanation excluded by assertion rather than by evidence.

full rationale

The paper's reproduction of bond lengths, angles, and dipole moments against external references (Table 1 and refs [3,12-15]) is a legitimate external check for the optimized geometries, and those parts are not circular. The circularity is confined to the interpretation of the torsion scans: the 'singularity' is read off the same QM energy/dipole data whose interpretation is at issue, and the only cited authority for the claim that weak H-bonds must show such break-point signatures is the author's own prior work (refs [9,10]). No convergence study, no initial-guess variation, no state-crossing or conical-intersection analysis, and no comparison between SCF algorithms is reported, so the paper never separates a physical electronic-energy singularity from a scan hopping between self-consistent solutions. The MMFF/SYBYL control is also not a clean null because Table 1 gives MMFF an N1-H5 distance of 1.1 Å rather than an H-bond distance, so the QM/MM contrast may be a conformer difference rather than evidence of a quantum singularity. The Mathieu-equation paragraph is an analogy about torsional eigenstates, not a derivation that the electronic energy as a function of torsion must be singular. Because the central claim therefore rests on a self-citation chain plus an asserted rather than demonstrated exclusion of numerical artifacts, the paper is partially circular: score 6.

Assumptions & free parameters 0 free parameters · 3 assumptions · 1 invented entities

The paper's central claim rests on three unverified assumptions: that SCF convergence is reliable across the scan, that the observed discontinuities are physical rather than numerical, and that MM and QM are comparing the same hydrogen-bonded dimer. No free parameters are fitted, but the MMFF minimum itself (N-H = 1.1 Å) is an unphysical reference. The only invented entity is the 'anisotropic energy singularity', which has no independent evidence.

assumptions (3)
  • domain assumption SCF convergence to the same electronic state is achieved at every torsion step.
    The claim of a physical singularity presumes the QM calculator is correctly converged at each geometry. The forward/reverse hysteresis and many discrete dipole levels contradict this, but the paper does not test it.
  • ad hoc to paper Observed discontinuities in QM energy are intrinsic quantum features, not numerical artifacts.
    Stated in the introduction ('inherent in the system and not due to any numerical convergence issues') and conclusion ('purely due to the quantum nature of bond break-up'), with no supporting test or independent evidence.
  • domain assumption MMFF/SYBYL optimized geometry is a valid reference for the same H-bonded dimer.
    Table 1 shows MMFF N1H5 distance of 1.1 Å, which is a covalent N-H bond distance, not a hydrogen bond. The comparison of MM and QM energy surfaces assumes both describe the same physical species, which is false.
invented entities (1)
  • Anisotropic energy singularities
    purpose: Explains the non-smooth, discontinuous QM torsion energy profiles as an inherent quantum feature of weak H-bond breaking.
    The 'singularity' is asserted to be physical, but no falsifiable prediction is offered beyond the same QM outputs that motivated the claim. No mechanism, no symmetry analysis, no direct observable is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Electronic Energy Singularities of Weakly H-bonded Ammonium Dimer." pith.science (2026). https://pith.science/paper/FDMGFE5D

@misc{pith2026250417107,
  author       = {Pith},
  title        = {Pith review of: Electronic Energy Singularities of Weakly H-bonded Ammonium Dimer},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FDMGFE5D}},
  note         = {Machine review of arXiv:2504.17107}
}
read the original abstract

Quantum and molecular mechanics based electronic energy studies of weak H-bonded ammonium dimer show distinctive feature in energy profile when computed by different QM methods contrast to MM methods. MM based MMFF and SYBYL methods show smoothly varying dihedral energy profile for torsion angle variation around weak N1-H5 held by H-bond strength of around 13 KJ/mol. All the QM based methods HF, B3LYP and MP2 show noisy and unstable torsion dependent electronic energy profile for H-bonded ammonium dimer. Exploring energy surface beyond bond length shows singularities and discontinuities. QM-based computation of dipole moment shows several discreet values with jumps and discontinuities with torsion angle variation for ammonium dimer. Also repeated computations and reverse torsion energy profile show persistent singularity feature observed in all standard QM techniques. KEY WORDS: ammonium dimer, H-bond, quantum signature, anisotropic energy singularities

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 15 canonical work pages

  1. [1]

    Ammonia dimer: A surprising structure

    DD Nelson, GT Fraser, W Klemperer: “Ammonia dimer: A surprising structure” Journal of Chemical Physics, 83, pp 6201–6208 1985 https://doi.org/10.1063/1.449566

    work page doi:10.1063/1.449566 1985

  2. [2]

    Ammonia Dimer, Linear or Cyclic?

    J Sadlej and L Lapinski, “Ammonia Dimer, Linear or Cyclic?” Journal of Molecular Structure: THEOCHEM 139, pp 233-240 1986 https://doi.org/10.1016/0166-1280(86)87040-3

    work page doi:10.1016/0166-1280(86)87040-3 1986

  3. [3]

    Extensive Theoretical Studies of the Hydrogen‐bonded Complexes (H2O)2, (H2O)2H+, (HF)2, (HF)2H+, F2H−, and (NH3)2

    MJ Frisch, JED Bene, JS Binkley, HF Schaefer: “Extensive Theoretical Studies of the Hydrogen‐bonded Complexes (H2O)2, (H2O)2H+, (HF)2, (HF)2H+, F2H−, and (NH3)2.” Journal of Chemical Physics 84 pp2279–2289 1986 https://doi.org/10.1063/1.450390

    work page doi:10.1063/1.450390 1986

  4. [4]

    The Ammonia Dimer Spectrum in Cold Helium Clusters

    M Behrens, U Buck, R Fröchtenicht, M Hartmann; M Havenith: “The Ammonia Dimer Spectrum in Cold Helium Clusters.”Journal of Chemical Physics 107 pp7179–7186 1997 https://doi.org/10.1063/1.474957

    work page doi:10.1063/1.474957 1997

  5. [5]

    Measurement of the Intermolecular Vibration–Rotation Tunneling Spectrum of the Ammonia Dimer by Tunable Far Infrared Laser Spectroscopy

    M. Havenith, RC Cohen, KL Busarow, DH Gwo, YT Lee, RJ Saykally: “Measurement of the Intermolecular Vibration–Rotation Tunneling Spectrum of the Ammonia Dimer by Tunable Far Infrared Laser Spectroscopy.” Journal of Chemical Physics 94 pp4776–4789 1991 https://doi.org/10.1063/1.460562

    work page doi:10.1063/1.460562 1991

  6. [6]

    Ammonia Dimer: Further Structural Studies

    DD Nelson, W. Klemperer, GT. Fraser; FJ Lovas, RD Suenram “Ammonia Dimer: Further Structural Studies.” Journal of Chemical Physics 87 pp 6364–6372 1987 https://doi.org/10.1063/1.453466

    work page doi:10.1063/1.453466 1987

  7. [7]

    Hydrogen Bonding in Liquid Ammonia

    A Krishnamoorthy K Nomura, N Baradwaj, K Shimamura, et al., “Hydrogen Bonding in Liquid Ammonia.” The Journal of Physical Chemistry Letters, 13 , pp 7051–7057 2022 https://doi.org/10.1021/acs.jpclett.2c01608 [8]A Jing, K Szalewicz, van der Avoird: "Ammonia Dimer: Extremely Fluxional but still hydrogen bonded. Nature Communications, 13 1470 2022 https://do...

    work page doi:10.1021/acs.jpclett.2c01608 2022

  8. [9]

    Observation of Quantum Signature in Rivastigmine Chemical Bond Break-up and Quantum Energetics, Spectral Studies of Anti-Alzheimer Inhibitors

    M. Rejwan Ali and Mihaly Mezei, “Observation of Quantum Signature in Rivastigmine Chemical Bond Break-up and Quantum Energetics, Spectral Studies of Anti-Alzheimer Inhibitors.” Journal of Biomolecular Structure and Dynamics, 39 pp 118–128. 2021 https://doi.org/10.1080/07391102.2019.1708462 [10]MR Ali, “Quantum Signature of Anisotropic Singularities in Hyd...

    arXiv 2021

Show all 16 references
  1. [11]

    Microwave Molecular Spectra

    W Gordy and RL Cook, "Microwave Molecular Spectra" Wiley, New York 1984

    1984

  2. [12]

    The Ammonia Dimer Equilibrium Dissociation Energy: Convergence to the Basis Set Limit at the Correlated Level

    L Stalring M Schutz R Lindh G Karlström P Widmar “The Ammonia Dimer Equilibrium Dissociation Energy: Convergence to the Basis Set Limit at the Correlated Level.” Molecular Physics 100 pp 3389-3399 2002 https://doi.org/10.1080/00268970210162718

    2002 doi

  3. [13]

    Ab Initio Search for the Equilibrium Structure of the Ammonia Dimer

    F Tao and W Klemperer, “Ab Initio Search for the Equilibrium Structure of the Ammonia Dimer.” Journal of Chemical Physics 99 pp 5976–5982 1993 https://doi.org/10.1063/1.465896

    1993 doi

  4. [14]

    Ab Initio and DFT Calculations of Some Weakly Bound Dimers and Complexes. I. The Dimers of Ammonia and Phosphine

    JA Altmann, MG Govender TA Ford “Ab Initio and DFT Calculations of Some Weakly Bound Dimers and Complexes. I. The Dimers of Ammonia and Phosphine.” Molecular Physics 103 pp 949-961 2005 https://doi.org/10.1080/00268470412331333555

    2005 doi

  5. [15]

    Structure, Internal Mobility, and Spectrum of the Ammonia Dimer: Calculation of the Vibration–Rotation‐tunneling States

    EHT Olthof; A van der Avoird; PES Wormer “Structure, Internal Mobility, and Spectrum of the Ammonia Dimer: Calculation of the Vibration–Rotation‐tunneling States.” Journal of Chemical. Physics 101 pp8430–8442 1994 https://doi.org/10.1063/1.468105

    1994 doi

  6. [16]

    Structure and Properties of Liquid Ammonia

    WL Jorgensen and M Ibrahim: “Structure and Properties of Liquid Ammonia.” Journal of the American Chemical Society, 102 pp3309–3315 1980 https://doi.org/10.1021/ja00530a001

    1980 doi

  7. [17]

    Intermolecular Potentials for Ammonia Based on the Test Particle Model and the Coupled Pair Functional Method

    KP Sagarik, R Ahlrichs, S Brode: “Intermolecular Potentials for Ammonia Based on the Test Particle Model and the Coupled Pair Functional Method”, Molecular Physics 57 pp 1247–1264 1986 https://doi.org/10.1080/00268978600100891

    1986 doi

  8. [18]

    Exploring Chemistry with Electronic Structure Methods: A Guide to Using Gaussian

    JB Foresman, AE Frisch: "Exploring Chemistry with Electronic Structure Methods: A Guide to Using Gaussian". 3rd ed. Wallingford CT 2015 Electronic Energy Singularities of Weakly H-bonded Ammonium Dimer Figure Slides Fig.1A & 1B: Optimized Geometry of Ammonium Dimer Computed by...

    2015

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.