Eine zuverlässige und effiziente Berechnungsmethode zum Auffinden von Übergangszuständen in chemischen Reaktionen

Es wurde eine rechnerische Methode zum Auffinden von Übergangszuständen in chemischen Reaktionen entwickelt, die den Rechenaufwand bei hoher Zuverlässigkeit erheblich reduziert. Im Vergleich zur am weitesten verbreiteten bestehenden Methode reduziert die vorliegende Methode den Gesamtrechenaufwand um etwa 50 bis 70 %.

Die Entwicklung, verfügbar auf GitHubist bereit, Fortschritte in der Materialwissenschaft zu beschleunigen und die Erforschung chemischer Reaktionen zugänglicher und effizienter zu machen. Dies könnte zu schnelleren wissenschaftlichen Entdeckungen und technologischen Innovationen führen.

Bei chemischen Reaktionen gehen Stoffe von einem energetisch stabilen Zustand in einen anderen über und durchlaufen dabei einen instabilen Übergangszustand. Dieser Vorgang ähnelt dem Finden der Route mit der niedrigsten Höhe über einen Berg, wenn man von einer Seite zur anderen überquert. Das Verständnis des Übergangszustands – des Gipfels dieses metaphorischen Bergpfads – ist entscheidend für ein tiefgreifendes Verständnis der Reaktionsmechanismen.

Aufgrund der vorübergehenden und instabilen Natur dieser Zustände ist ihre experimentelle Beobachtung und Identifizierung jedoch eine Herausforderung und erfordert oft eine rechnerische Untersuchung.

Diese Studie konzentriert sich auf rechnerische Methoden zum Finden eines Übergangszustands zwischen einem bekannten Reaktanten und einem Produkt. Diese Art der Übergangszustandssuche optimiert den Pfad, der Produkt und Reaktant verbindet, sodass es den Übergangszustand durchquert. Da der Pfad normalerweise durch mehrere Punkte auf dem Pfad (oft als Bilder bezeichnet) dargestellt wird, wird der Pfad tatsächlich durch die schrittweise Aktualisierung der Bilder optimiert.

Die heute am häufigsten verwendete Methode ist die Nudged Elastic Band (NEB)-Methode. Eine der größten Herausforderungen dieser Methode besteht darin, dass sie rechenintensiv ist. Dafür gibt es zwei Hauptgründe. Zum einen ist eine große Anzahl von Bildern erforderlich, um die Auflösung der Suche zu erhöhen. Der andere Grund besteht darin, dass das Suchprinzip nicht variiert (also eine Zielfunktion minimiert) und daher tendenziell auch die Anzahl der Aktualisierungen pro Bild groß ist.

Die in dieser Studie neu implementierte Methode löst diese Probleme auf innovative Weise. Erstens kann die Anzahl der Bilder auf etwa 3 reduziert werden, da nur die Region um den Übergangszustand intensiv durchsucht wird. Darüber hinaus ist das Suchprinzip variierend und kann daher effizienter gelöst werden. Konkret ist die Zielfunktion als Linienintegral der Exponentialfunktion der Energie entlang des Pfades definiert.

Die Leistung der neuen Methode wurde anhand von 121 chemischen Reaktionen bewertet und die Ergebnisse mit der NEB-Methode und ihrer verbesserten Version verglichen. Erstens identifizierte die vorliegende Methode Übergangszustände in 98 % der Fälle korrekt. Diese Genauigkeit ist viel höher als bei der NEB-Methode und vergleichbar mit der verbesserten Version. Zweitens zeigte die vorliegende Methode eine deutliche Reduzierung der Gesamtrechenkosten – etwa 70 % weniger als die NEB-Methode und 50 % weniger als ihre verbesserte Version.

Um eine breitere Anwendung zu ermöglichen, haben die Forscher ihr Rechenprogramm auf GitHub verfügbar gemacht. Es ist in Python geschrieben und für die Verwendung mit der Atomic Simulation Environment (ASE) konzipiert und ermöglicht Forschern die einfache Untersuchung von Übergangszuständen durch die Angabe von Reaktanten und Produkten.

Mit Blick auf die Zukunft sind die Auswirkungen dieser Forschung enorm. Indem die Methode die Suche nach Übergangszuständen einfacher und schneller macht, ist sie bereit, die Forschung und Entwicklung in allen Bereichen der Naturwissenschaften mithilfe der Computerchemie zu beschleunigen.

Die Forschung ist veröffentlicht im Zeitschrift für chemische Theorie und Berechnung.