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

オープン
  • 何か問題がありますか？
#764 コメント 2 件 リアクション 0 件 担当者 1 名 GitHub で見る

@ohlincha がすでに取り組んでいます。

2026年11月1日 から。

評価

この issue はまだ評価されていません。

説明

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.

主要言語
C
スター
10
フォーク
8
平均マージ
21分
マージ済み PR（30日）
1

環境構築

このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。

はじめの一歩

  1. issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
  2. 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
  3. リポジトリをフォークし、ブランチを切って変更します。
  4. issue 番号を参照したプルリクエストを送ります。

FriendsofECCE/ECCE のほかの issue

FriendsofECCE/ECCE の issue をすべて見る

似ている issue

C の issue をもっと見る

新しい issue をメールで受け取る

初心者向けの GitHub issue を短くまとめたダイジェスト。