ECCE-QM IRC: both directions end in the same valley (SN2 backward branch reaches the product complex); false 'minimum reached'

Offen
#764 2 Kommentare 0 Reaktionen 1 zugewiesene Person Auf GitHub ansehen

@ohlincha arbeitet bereits daran.

Seit 01.11.2026.

Bewertung

Dieses Issue wurde noch nicht bewertet.

Beschreibung

SN2 Cl⁻ + CH₃Br, B3LYP/aug-cc-pVDZ, charge −1, from the ORCA TS geometry (radium, ecce-qm from main, 2026-11-01).

  • irc_step 0.1 (default), irc_max_points 40: both directions run to the 40-point limit; backward at s = −1.98 E = −3074.14836, s = −2.98 E = −3074.15365, s = −3.88 E = −3074.15546 Eh, still descending.
  • irc_step 0.2, irc_max_points 80: ends with 'forward 32 points, gradient below threshold: minimum reached; backward 21 points, gradient below threshold: minimum reached'. Backward at s = −1.99 E = −3074.14576, s = −2.98 E = −3074.14774, last point s = −4.08 E = −3074.14833 Eh, still descending by ~5e-5 Eh per step.

So with the larger step the backward branch is already 2.6 mEh above the step-0.1 path at s = −2, flattens, and is reported as a minimum about 7 mEh above where the step-0.1 run already was. Forward looks consistent.

To find: whether the corrector step loses the path at 0.2 (step-size control, Hessian update), and why the minimum test fires while the energy still decreases (gradient threshold alone; add an energy-change criterion or an optimisation of the end point). A false 'minimum reached' must not be reported.

Vorherrschende Sprache
C
Sterne
10
Forks
8
Ø Merge
21 Min.
Gemergte PRs (30 T.)
1

Entwicklungsumgebung

Dieses Projekt bietet weder Dev-Container noch Dockerfile noch Beitragsleitfaden – die Einrichtung liegt bei Ihnen. Beginnen Sie mit der README; die allgemeinen Schritte stehen in unserem Leitfaden für den ersten Beitrag.

Erste Schritte

  1. Lesen Sie das ganze Issue und danach den Beitragsleitfaden des Projekts.
  2. Schreiben Sie ins Issue, dass Sie es übernehmen — das erspart doppelte Arbeit.
  3. Forken Sie das Repository und arbeiten Sie in einem Branch.
  4. Öffnen Sie einen Pull Request, der die Issue-Nummer nennt.

Mehr aus FriendsofECCE/ECCE

Alle Issues in FriendsofECCE/ECCE

Ähnliche Issues

Weitere Issues zu C

Neue Issues direkt in Ihr Postfach

Eine kurze Übersicht über anfängerfreundliche GitHub-Issues.