-
Spin parameter optimization for spin-polarized extended tight-binding methods
Authors:
Siyavash Moradi,
Rebecca Tomann,
Josie Hendrix,
Martin Head-Gordon,
Christopher J. Stein
Abstract:
We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential par…
▽ More
We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential parameters for a given reference dataset, allowing for a nuanced understanding of their impact on diverse molecular properties. In the case of the W4-11 dataset, substantial error reduction was achieved, demonstrating the potential of the optimization. Transferability of the optimized spin-polarization constants over different properties, however, is limited, as we demonstrate by applying the optimized parameters on a set of singlet-triplet gaps in carbenes. Further studies on ionization potentials and electron affinities highlight some inherent limitations of current extended tight-binding methods that can not be resolved by simple parameter optimization. We conclude that the significantly improved accuracy strongly encourages the present re-optimization of the spin-polarization constants, whereas the limited transferability motivates a property-specific optimization strategy.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Benchmarking DFT-based excited-state methods for intermolecular charge-transfer excitations
Authors:
Nicola Bogo,
Christopher J. Stein
Abstract:
Intermolecular charge-transfer is a highly important process in biology and energy-conversion applications where generated charges need to be transported over several moieties. However, its theoretical description is challenging since the high accuracy required to describe these excited states must be accessible for calculations on large molecular systems. In this benchmark study, we identify reli…
▽ More
Intermolecular charge-transfer is a highly important process in biology and energy-conversion applications where generated charges need to be transported over several moieties. However, its theoretical description is challenging since the high accuracy required to describe these excited states must be accessible for calculations on large molecular systems. In this benchmark study, we identify reliable low-scaling computational methods for this task. Our reference results were obtained from highly accurate wavefunction calculations that restrict the size of the benchmark systems. However, the density-functional theory based methods that we identify as accurate can be applied to much larger systems. Since targeting charge-transfer states requires the unambiguous classification of an excited state, we first analyze several charge-transfer descriptors for their reliability concerning intermolecular charge-transfer and single out DCT as an optimal choice for our purposes. In general, best results are obtained for orbital-optimized methods - and among those, IMOM proved to be the most numerically stable variant - but optimally-tuned range-separated hybrid functionals combined with rather small basis sets proved to yield surprisingly good results. This makes these fast calculations attractive for high-throughput screening applications.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Automated, Consistent, and Even-handed Selection of Active Orbital Spaces for Quantum Embedding
Authors:
Elena Kolodzeiski,
Christopher J. Stein
Abstract:
A widely used strategy to reduce the computational cost in quantum-chemical calculations is to partition the system into an active subsystem, which is the focus of the computational efforts and an environment that is treated at a lower computational level. The system partitioning is mostly based on localized molecular orbitals. When reaction paths or energy differences are to be calculated, it is…
▽ More
A widely used strategy to reduce the computational cost in quantum-chemical calculations is to partition the system into an active subsystem, which is the focus of the computational efforts and an environment that is treated at a lower computational level. The system partitioning is mostly based on localized molecular orbitals. When reaction paths or energy differences are to be calculated, it is crucial to keep the orbital space consistent for all structures. Inconsistencies in the orbital space can lead to unpredictable errors in the potential energy surface. While successful strategies to ensure this consistency have been established for organic and even metal-organic systems, these methods often fail for metal clusters or nanoparticles with a high density of near-degenerate and delocalized molecular orbitals. However, such systems are highly relevant for catalysis. Accurate yet feasible quantum-mechanical ab initio calculations are therefore highly desired. In this work, we present an approach based on the SPADE algorithm that allows us to ensure an automated and consistent partitioning even for systems with delocalized and near-degenerate molecular orbitals and demonstrate the validity of this method for the binding energies of small molecules on transition-metal clusters.
△ Less
Submitted 15 June, 2023;
originally announced June 2023.
-
A statistical perspective on microsolvation
Authors:
Mohammad Rahbar,
Christopher J. Stein
Abstract:
The lack of a procedure to determine equilibrium thermodynamic properties of a small system interacting with a bath is frequently seen as a weakness of conventional statistical mechanics. A typical example for such a small system is a solute surrounded by an explicit solvation shell. One way to approach this problem is to enclose the small system of interest in a large bath of explicit solvent mol…
▽ More
The lack of a procedure to determine equilibrium thermodynamic properties of a small system interacting with a bath is frequently seen as a weakness of conventional statistical mechanics. A typical example for such a small system is a solute surrounded by an explicit solvation shell. One way to approach this problem is to enclose the small system of interest in a large bath of explicit solvent molecules, considerably larger than the system itself. The explicit inclusion of the solvent degrees of freedom is obviously limited by the available computational resources. A potential remedy to this problem is a microsolvation approach where only a few explicit solvent molecules are considered and surrounded by an implicit solvent bath. Still, the sampling of the solvent degrees of freedom is challenging with conventional grand canonical Monte Carlo methods, since no single chemical potential for the solvent molecules can be defined in the realm of small-system thermodynamics. In this work, a statistical thermodynamic model based on the grand canonical ensemble is proposed that avoids the conventional system size limitations and accurately characterizes the properties of the system of interest subject to the thermodynamic constraints of the bath. We extend an existing microsolvation approach to a generalized multi-bath "micro-statistical" model and show that the previously derived approaches result as a limit of our model. The framework described here is universal and we validate our method numerically for a Lennard-Jones model fluid.
△ Less
Submitted 14 December, 2022;
originally announced December 2022.
-
High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation
Authors:
Jan P. Unsleber,
Hongbin Liu,
Leopold Talirz,
Thomas Weymuth,
Maximilian Mörchen,
Adam Grofe,
Dave Wecker,
Christopher J. Stein,
Ajay Panyala,
Bo Peng,
Karol Kowalski,
Matthias Troyer,
Markus Reiher
Abstract:
Quantum chemical calculations on atomistic systems have evolved into a standard approach to study molecular matter. These calculations often involve a significant amount of manual input and expertise although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for e…
▽ More
Quantum chemical calculations on atomistic systems have evolved into a standard approach to study molecular matter. These calculations often involve a significant amount of manual input and expertise although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput lectronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject hem to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all omponents of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow at the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.
△ Less
Submitted 13 April, 2023; v1 submitted 26 November, 2022;
originally announced November 2022.
-
NewtonNet: A Newtonian message passing network for deep learning of interatomic potentials and forces
Authors:
Mojtaba Haghighatlari,
Jie Li,
Xingyi Guan,
Oufan Zhang,
Akshaya Das,
Christopher J. Stein,
Farnaz Heidar-Zadeh,
Meili Liu,
Martin Head-Gordon,
Luke Bertels,
Hongxia Hao,
Itai Leven,
Teresa Head-Gordon
Abstract:
We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from trainable latent force vectors, and physics-infused operators that are inspired by the Newtonian physics, the entire model remains rotationally equivariant, and many-body interactions are inferr…
▽ More
We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from trainable latent force vectors, and physics-infused operators that are inspired by the Newtonian physics, the entire model remains rotationally equivariant, and many-body interactions are inferred by more interpretable physical features. We test NewtonNet on the prediction of several reactive and non-reactive high quality ab initio data sets including single small molecule dynamics, a large set of chemically diverse molecules, and methane and hydrogen combustion reactions, achieving state-of-the-art test performance on energies and forces with far greater data and computational efficiency than other deep learning models.
△ Less
Submitted 5 August, 2021;
originally announced August 2021.
-
The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients
Authors:
Christopher J. Stein,
John M. Herbert,
Martin Head-Gordon
Abstract:
We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions and…
▽ More
We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions and a Stern-layer correction. Recasting the Poisson-Boltzmann equations into Euler-Lagrange equations then significantly simplifies the derivation of the free energy of solvation for these approximate models. The parameters of the model are then either fit directly to experimental observables, e.g. the finite ion size, or optimized for agreement with experimental results. Experimental data for this optimization is available in the form of Sechenov coefficients that describe the linear dependence of the salting-out effect of solutes with respect to the electrolyte concentration. In the final part we rationalize the qualitative disagreement of the finite ion size modification to the Poisson-Boltzmann model with experimental observations by taking into account the electrolyte concentration dependence of the Stern layer. A route towards a revised model that captures the experimental observations while including the finite ion size effects is then outlined. This implementation paves the way for the study of electrochemical and electrocatalytic processes of molecules and cluster models with accurate electronic structure methods.
△ Less
Submitted 8 October, 2019;
originally announced October 2019.
-
Semiclassical Dispersion Corrections efficiently improve Multi-Configurational Theory with Short-Range Density-Functional Dynamic Correlation
Authors:
Christopher J. Stein,
Markus Reiher
Abstract:
Multi-configurational wave functions are known to describe electronic structure across a Born-Oppenheimer surface qualitatively correct. However, for quantitative reaction energies, dynamical correlation originating from the many configurations involving excitations out of the restricted orbital space, the active space, must be considered. Standard procedures involve approximations that eventually…
▽ More
Multi-configurational wave functions are known to describe electronic structure across a Born-Oppenheimer surface qualitatively correct. However, for quantitative reaction energies, dynamical correlation originating from the many configurations involving excitations out of the restricted orbital space, the active space, must be considered. Standard procedures involve approximations that eventually limit the ultimate accuracy achievable (most prominently, multi-reference perturbation theory). At the same time, the computational cost increase dramatically due to the necessity to obtain higher-order reduced density matrices. It is this disproportion that leads us here to propose a MC-srDFT-D hybrid approach of semiclassical dispersion (D) corrections to cover long-range dynamical correlation in a multi-configurational (MC) wave function theory which includes short-range (sr) dynamical correlation by density functional theory (DFT) without double counting. We demonstrate that the reliability of this approach is very good (at negligible cost), especially when considering that standard second-order multi-reference perturbation theory usually overestimates dispersion interactions.
△ Less
Submitted 17 March, 2020; v1 submitted 22 August, 2019;
originally announced August 2019.
-
autoCAS: a program for fully automated multi-configurational calculations
Authors:
Christopher J. Stein,
Markus Reiher
Abstract:
We present our implementation autoCAS for fully automated multi-configurational calculations, which we also make available free of charge on our webpages. The graphical user interface of autoCAS connects a general electronic structure program with a density matrix renormalization group program to carry out our recently introduced automated active space selection protocol for multi-configurational…
▽ More
We present our implementation autoCAS for fully automated multi-configurational calculations, which we also make available free of charge on our webpages. The graphical user interface of autoCAS connects a general electronic structure program with a density matrix renormalization group program to carry out our recently introduced automated active space selection protocol for multi-configurational calculations [J. Chem. Theory Comput., 2016, 12, 1760]. Next to this active space selection, autoCAS carries out several steps of multi-configurational calculations so that only a minimal input is required to start them, comparable to that of a standard Kohn-Sham density functional theory calculation, so that black-box multi-configurational calculations become feasible. Furthermore, we introduce a new extension to the selection algorithm that facilitates automated selections for molecules with large valence orbital spaces consisting of several hundred orbitals.
△ Less
Submitted 26 May, 2019; v1 submitted 29 March, 2019;
originally announced April 2019.
-
Optimization of highly excited matrix product states with an application to vibrational spectroscopy
Authors:
Alberto Baiardi,
Christopher J. Stein,
Vincenzo Barone,
Markus Reiher
Abstract:
Configuration-interaction-type calculations on electronic and vibrational structure are often the method of choice for the reliable approximation of many-particle wave functions and energies. The exponential scaling, however, limits their application range. An efficient approximation to the full configuration interaction solution can be obtained with the density matrix renormalization group (DMRG)…
▽ More
Configuration-interaction-type calculations on electronic and vibrational structure are often the method of choice for the reliable approximation of many-particle wave functions and energies. The exponential scaling, however, limits their application range. An efficient approximation to the full configuration interaction solution can be obtained with the density matrix renormalization group (DMRG) algorithm without a restriction to a predefined excitation level. In a standard DMRG implementation, however, excited states are calculated with a ground-state optimization in the space orthogonal to all lower lying wave function solutions. A trivial parallelization is therefore not possible and the calculation of highly excited states becomes prohibitively expensive, especially in regions with a high density of states. Here, we introduce two variants of the density matrix renormalization group algorithm that allow us to target directly specific energy regions and therefore highly excited states. The first one, based on shift-and-invert techniques, is particularly efficient for low-lying states, but is not stable in regions with a high density of states. The second one, based on the folded auxiliary operator, is less efficient, but more accurate in targeting high-energy states. We apply the algorithm to the solution of the nuclear Schroedinger equation, but emphasize that it can be applied to the diagonalization of general Hamiltonians as well, such as the electronic Coulomb Hamiltonian to address X-ray spectra. In combination with several root-homing algorithms and a stochastic sampling of the determinant space, excited states of interest can be adequately tracked and analyzed during the optimization. We demonstrate that we can accurately calculate prominent spectral features of large molecules such as the sarcosyn-glycine dipeptide.
△ Less
Submitted 16 January, 2019; v1 submitted 18 October, 2018;
originally announced October 2018.
-
Vibrational Density Matrix Renormalization Group
Authors:
Alberto Baiardi,
Christopher J. Stein,
Vincenzo Barone,
Markus Reiher
Abstract:
Variational approaches for the calculation of vibrational wave functions and energies are a natural route to obtain highly accurate results with controllable errors. However, the unfavorable scaling and the resulting high computational cost of standard variational approaches limit their application to small molecules with only few vibrational modes. Here, we demonstrate how the density matrix reno…
▽ More
Variational approaches for the calculation of vibrational wave functions and energies are a natural route to obtain highly accurate results with controllable errors. However, the unfavorable scaling and the resulting high computational cost of standard variational approaches limit their application to small molecules with only few vibrational modes. Here, we demonstrate how the density matrix renormalization group (DMRG) can be exploited to optimize vibrational wave functions (vDMRG) expressed as matrix product states. We study the convergence of these calculations with respect to the size of the local basis of each mode, the number of renormalized block states, and the number of DMRG sweeps required. We demonstrate the high accuracy achieved by vDMRG for small molecules that were intensively studied in the literature. We then proceed to show that the complete fingerprint region of the sarcosyn-glycin dipeptide can be calculated with vDMRG.
△ Less
Submitted 3 July, 2017; v1 submitted 27 March, 2017;
originally announced March 2017.
-
Automated Identification of Relevant Frontier Orbitals for Chemical Compounds and Processes
Authors:
Christopher J. Stein,
Markus Reiher
Abstract:
Quantum-chemical multi-configurational methods are required for a proper description of static electron correlation, a phenomenon inherent to the electronic structure of molecules with multiple (near-)degenerate frontier orbitals. Here, we review how a property of these frontier orbitals, namely the entanglement entropy is related to static electron correlation. A subset of orbitals, the so-called…
▽ More
Quantum-chemical multi-configurational methods are required for a proper description of static electron correlation, a phenomenon inherent to the electronic structure of molecules with multiple (near-)degenerate frontier orbitals. Here, we review how a property of these frontier orbitals, namely the entanglement entropy is related to static electron correlation. A subset of orbitals, the so-called active orbital space, is an essential ingredient for all multi-configurational methods. We proposed an automated selection of this active orbital space, that would otherwise be a tedious and error prone manual procedure, based on entanglement measures. Here, we extend this scheme to demonstrate its capability for the selection of consistent active spaces for several excited states and along reaction coordinates.
△ Less
Submitted 1 February, 2017;
originally announced February 2017.
-
Measuring Multi-Configurational Character by Orbital Entanglement
Authors:
Christopher J. Stein,
Markus Reiher
Abstract:
One of the most critical tasks at the very beginning of a quantum chemical investigation is the choice of either a multi- or single-configurational method. Naturally, many proposals exist to define a suitable diagnostic of the multi-configurational character for various types of wave functions in order to assist this crucial decision. Here, we present a new orbital-entanglement based multi-configu…
▽ More
One of the most critical tasks at the very beginning of a quantum chemical investigation is the choice of either a multi- or single-configurational method. Naturally, many proposals exist to define a suitable diagnostic of the multi-configurational character for various types of wave functions in order to assist this crucial decision. Here, we present a new orbital-entanglement based multi-configurational diagnostic termed $Z_{s(1)}$. The correspondence of orbital entanglement and static (or nondynamic) electron correlation permits the definition of such a diagnostic. We chose our diagnostic to meet important requirements such as well-defined limits for pure single-configurational and multi-configurational wave functions. The $Z_{s(1)}$ diagnostic can be evaluated from a partially converged, but qualitatively correct, and therefore inexpensive density matrix renormalization group wave function as in our recently presented automated active orbital selection protocol. Its robustness and the fact that it can be evaluated at low cost make this diagnostic a practical tool for routine applications.
△ Less
Submitted 16 November, 2016; v1 submitted 8 September, 2016;
originally announced September 2016.
-
The Delicate Balance of Static and Dynamic Electron Correlation
Authors:
Christopher J. Stein,
Vera von Burg,
Markus Reiher
Abstract:
Multi-configurational approaches yield universal wave function parameterizations that can qualitatively well describe electronic structures along reaction pathways. For quantitative results, multi-reference perturbation theory is required to capture dynamic electron correlation from the otherwise neglected virtual orbitals. Still, the overall accuracy suffers from the finite size and choice of the…
▽ More
Multi-configurational approaches yield universal wave function parameterizations that can qualitatively well describe electronic structures along reaction pathways. For quantitative results, multi-reference perturbation theory is required to capture dynamic electron correlation from the otherwise neglected virtual orbitals. Still, the overall accuracy suffers from the finite size and choice of the active orbital space and peculiarities of the perturbation theory. Fortunately, the electronic wave functions at equilibrium structures of reactants and products can often be well described by single-reference methods and hence are accessible to accurate coupled cluster calculations. Here, we calculate the heterolytic double dissociation energy of four 3d-metallocenes with the complete active space self-consistent field method and compare to highly accurate coupled cluster data. Our coupled cluster data are well within the experimental error bars. This accuracy can also be approached by complete active space calculations with an orbital selection based on information entropy measures. The entropy based active space selection is discussed in detail. We find a very subtle balance between static and dynamic electron correlation effects that emphasizes the need for algorithmic active space selection and that differs significantly from restricted active space results for identical active spaces reported in the literature.
△ Less
Submitted 23 May, 2016;
originally announced May 2016.
-
Automated Selection of Active Orbital Spaces
Authors:
Christopher J. Stein,
Markus Reiher
Abstract:
One of the key challenges of quantum-chemical multi-configuration methods is the necessity to manually select orbitals for the active space. This selection requires both expertise and experience and can therefore impose severe limitations on the applicability of this most general class of ab initio methods. A poor choice of the active orbital space may yield even qualitatively wrong results. This…
▽ More
One of the key challenges of quantum-chemical multi-configuration methods is the necessity to manually select orbitals for the active space. This selection requires both expertise and experience and can therefore impose severe limitations on the applicability of this most general class of ab initio methods. A poor choice of the active orbital space may yield even qualitatively wrong results. This is obviously a severe problem, especially for wave function methods that are designed to be systematically improvable. Here, we show how the iterative nature of the density matrix renormalization group combined with its capability to include up to about one hundred orbitals in the active space can be exploited for a systematic assessment and selection of active orbitals. These benefits allow us to implement an automated approach for active orbital space selection, which can turn multi-configuration models into black box approaches.
△ Less
Submitted 11 February, 2016;
originally announced February 2016.
-
New Approaches for ab initio Calculations of Molecules with Strong Electron Correlation
Authors:
Stefan Knecht,
Erik Donovan Hedegård,
Sebastian Keller,
Arseny Kovyrshin,
Yingjin Ma,
Andrea Muolo,
Christopher J. Stein,
Markus Reiher
Abstract:
Reliable quantum chemical methods for the description of molecules with dense-lying frontier orbitals are needed in the context of many chemical compounds and reactions. Here, we review developments that led to our newcomputational toolbo x which implements the quantum chemical density matrix renormalization group in a second-generation algorithm. We present an overview of the different components…
▽ More
Reliable quantum chemical methods for the description of molecules with dense-lying frontier orbitals are needed in the context of many chemical compounds and reactions. Here, we review developments that led to our newcomputational toolbo x which implements the quantum chemical density matrix renormalization group in a second-generation algorithm. We present an overview of the different components of this toolbox.
△ Less
Submitted 21 April, 2016; v1 submitted 31 December, 2015;
originally announced December 2015.