Sie haben nicht angegeben, bei welcher Geometrie Konvergenzprobleme auftreten. Daher dachte ich zunächst, dass dies passiert, wenn alle drei $ \ ce {N} $ -Atome weit voneinander entfernt sind. Aber dann habe ich die Berechnung selbst durchgeführt und festgestellt, dass dies nicht der Fall ist: Der SCF-Konvergenzfehler tritt bei R1 = 1.2
und R2 = 1.1
auf. Bei dieser Geometrie beginnt SCF nach 20 Iterationen zurück zu schwingen und zu erzwingen, wie in der folgenden Grafik gezeigt. Eine einfache Erhöhung der Anzahl der SCF-Iterationen würde also nicht helfen.
Für Aus diesem Grund habe ich schnell einige andere bekannte Tricks überprüft, um die Probleme mit der SCF-Konvergenz zu vermeiden (siehe z. B. hier), und ich empfehle die Verwendung eines quadratisch konvergenten SCF (QC) -SCF) durch Hinzufügen von scf = qc
zum Routenabschnitt:
# HF / 6-31G Testscan scf = qc
Beachten Sie jedoch, dass QC-SCF zwar zuverlässiger als der Standard-DIIS-Algorithmus ist, aber auch langsamer . Erwarten Sie daher eine Verlängerung der Rechenzeit. Aus diesem Grund gibt es auch eine interessante Option scf = xqc
, die zuerst den Standardalgorithmus ausprobiert und nur dann zu einem quadratisch konvergenten Algorithmus wechselt, wenn der Standard nicht konvergiert hat.
# HF / 6-31G-Testscan scf = xqc
Das Umschalten auf einen quadratisch konvergenten Algorithmus hat jedoch nicht geholfen, da SCF immer noch nicht konvergiert, sondern bei einigen anderen Geometrien. Daher habe ich die Angelegenheit genauer untersucht und zusätzlich zu scf = xqc
empfehle ich, sowohl die radiale als auch die Winkelflexibilität Ihres Basissatzes durch Hinzufügen von Diffus- und Polarisationsfunktionen zu verbessern hinein
# HF / 6-31 + G (d, p) Testscan scf = xqc
Nach dem Ausführen der Berechnung für a Stunde erreichte es 2656 Punkte (was ungefähr 40% aller Punkte entspricht) ohne SCF-Konvergenzprobleme, aber falls diese Variante immer noch fehlschlägt:
- Erwägen Sie, die Konvergenzkriterien zu senken, indem Sie
scf = xqc
in scf = (xqc, conver = n)
ändern, wobei n
beispielsweise ist , 6 (standardmäßig 8). - Wenn Fehler auftreten, wenn mindestens ein $ \ ce {N} $ weit von zwei anderen entfernt ist, sollten Sie die maximale Entfernung für Ihren PES-Scan verringern. Meiner Meinung nach sind 8,5 Å viel zu groß, 5 Å wären wahrscheinlich mehr als genug.
Endlich gelang mir der folgende Routenabschnitt
# HF / 6-31 + G (d, p) SCF scannen = (XQC, Konverter = 6)
Hier sind die Ergebnisse und unten ist die Darstellung der SPE. : D