The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations

Extending the Reach of Accurate Wavefunction Methods

Multiconfigurational quantum chemistry methods, and especially the multiconfigurational self-consistent field (MCSCF) and multireference perturbation theory (MRPT2), are powerful tools, particularly suited to the accurate modeling of photochemical processes and transition metal catalysis. However, they are limited by their high computational cost compared to other methods, especially density functional theory. Moreover, there are areas where they would be expected to perform well, but where they are not applied due to lack of experience. This thesis addresses those issues. First, the efficiency of the Cholesky decomposition approximation to reduce the cost of MCSCF and MRPT2 without sacrificing their accuracy is demonstrated. This then motivates the extension of the Cholesky approximation to the computation of MCSCF nuclear gradients, thus strongly improving the ability to perform MCSCF non-adiabatic molecular dynamics. Typically, a tenfold speed-up is observed allowing dynamic simulation of larger systems or over longer times. Finally, multiconfigurational methods are applied to the computation of X-ray spectra of transition metal complexes. The importance of the different parameters in the calculation is systematically investigated, laying the base for wider applications of those accurate methods in the modeling of X-ray spectroscopy. A tool to analyze the resulting spectrum in terms of molecular orbitals is also presented, strengthening the interplay between theory and experiments. With these developments and other significant ones that have happened in recent years, multiconfigurational methods can now reach new grounds and contribute to important new discoveriesQC 20210216</p

Extending the Reach of Accurate Wavefunction Methods

Multiconfigurational quantum chemistry methods, and especially the multiconfigurational self-consistent field (MCSCF) and multireference perturbation theory (MRPT2), are powerful tools, particularly suited to the accurate modeling of photochemical processes and transition metal catalysis. However, they are limited by their high computational cost compared to other methods, especially density functional theory. Moreover, there are areas where they would be expected to perform well, but where they are not applied due to lack of experience. This thesis addresses those issues. First, the efficiency of the Cholesky decomposition approximation to reduce the cost of MCSCF and MRPT2 without sacrificing their accuracy is demonstrated. This then motivates the extension of the Cholesky approximation to the computation of MCSCF nuclear gradients, thus strongly improving the ability to perform MCSCF non-adiabatic molecular dynamics. Typically, a tenfold speed-up is observed allowing dynamic simulation of larger systems or over longer times. Finally, multiconfigurational methods are applied to the computation of X-ray spectra of transition metal complexes. The importance of the different parameters in the calculation is systematically investigated, laying the base for wider applications of those accurate methods in the modeling of X-ray spectroscopy. A tool to analyze the resulting spectrum in terms of molecular orbitals is also presented, strengthening the interplay between theory and experiments. With these developments and other significant ones that have happened in recent years, multiconfigurational methods can now reach new grounds and contribute to important new discoverie

Extending the Reach of Accurate Wavefunction Methods

Multiconfigurational quantum chemistry methods, and especially the multiconfigurational self-consistent field (MCSCF) and multireference perturbation theory (MRPT2), are powerful tools, particularly suited to the accurate modeling of photochemical processes and transition metal catalysis. However, they are limited by their high computational cost compared to other methods, especially density functional theory. Moreover, there are areas where they would be expected to perform well, but where they are not applied due to lack of experience. This thesis addresses those issues. First, the efficiency of the Cholesky decomposition approximation to reduce the cost of MCSCF and MRPT2 without sacrificing their accuracy is demonstrated. This then motivates the extension of the Cholesky approximation to the computation of MCSCF nuclear gradients, thus strongly improving the ability to perform MCSCF non-adiabatic molecular dynamics. Typically, a tenfold speed-up is observed allowing dynamic simulation of larger systems or over longer times. Finally, multiconfigurational methods are applied to the computation of X-ray spectra of transition metal complexes. The importance of the different parameters in the calculation is systematically investigated, laying the base for wider applications of those accurate methods in the modeling of X-ray spectroscopy. A tool to analyze the resulting spectrum in terms of molecular orbitals is also presented, strengthening the interplay between theory and experiments. With these developments and other significant ones that have happened in recent years, multiconfigurational methods can now reach new grounds and contribute to important new discoverie

MultiPsi : A python-driven MCSCF program for photochemistry and spectroscopy simulations on modern HPC environments

We present MultiPsi, an open-source MCSCF program for the calculation of ground and excited states properties of strongly correlated systems. The program currently implements a general MCSCF code with excited states available using either state-averaging or linear response. It is written in a highly modular fashion using Python/C++ which makes it well suited as a development platform, enabling easy prototyping of novel methods, and as a teaching tool using interactive notebooks. The code is also very efficient and designed for modern high-performance computing environments using hybrid OpenMP/MPI parallelization. This efficiency is demonstrated with the calculation of the CASSCF energy and linear response of a molecule with more than 700 atoms as well as a fully optimized conventional CI calculation on more than 400 billion determinants. This article is categorized under: Software > Quantum Chemistry Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy

Extending the Reach of Accurate Wavefunction Methods

Multiconfigurational quantum chemistry methods, and especially the multiconfigurational self-consistent field (MCSCF) and multireference perturbation theory (MRPT2), are powerful tools, particularly suited to the accurate modeling of photochemical processes and transition metal catalysis. However, they are limited by their high computational cost compared to other methods, especially density functional theory. Moreover, there are areas where they would be expected to perform well, but where they are not applied due to lack of experience. This thesis addresses those issues. First, the efficiency of the Cholesky decomposition approximation to reduce the cost of MCSCF and MRPT2 without sacrificing their accuracy is demonstrated. This then motivates the extension of the Cholesky approximation to the computation of MCSCF nuclear gradients, thus strongly improving the ability to perform MCSCF non-adiabatic molecular dynamics. Typically, a tenfold speed-up is observed allowing dynamic simulation of larger systems or over longer times. Finally, multiconfigurational methods are applied to the computation of X-ray spectra of transition metal complexes. The importance of the different parameters in the calculation is systematically investigated, laying the base for wider applications of those accurate methods in the modeling of X-ray spectroscopy. A tool to analyze the resulting spectrum in terms of molecular orbitals is also presented, strengthening the interplay between theory and experiments. With these developments and other significant ones that have happened in recent years, multiconfigurational methods can now reach new grounds and contribute to important new discoveriesQC 20210216</p

Extending the Reach of Accurate Wavefunction Methods

Multiconfigurational quantum chemistry methods, and especially the multiconfigurational self-consistent field (MCSCF) and multireference perturbation theory (MRPT2), are powerful tools, particularly suited to the accurate modeling of photochemical processes and transition metal catalysis. However, they are limited by their high computational cost compared to other methods, especially density functional theory. Moreover, there are areas where they would be expected to perform well, but where they are not applied due to lack of experience. This thesis addresses those issues. First, the efficiency of the Cholesky decomposition approximation to reduce the cost of MCSCF and MRPT2 without sacrificing their accuracy is demonstrated. This then motivates the extension of the Cholesky approximation to the computation of MCSCF nuclear gradients, thus strongly improving the ability to perform MCSCF non-adiabatic molecular dynamics. Typically, a tenfold speed-up is observed allowing dynamic simulation of larger systems or over longer times. Finally, multiconfigurational methods are applied to the computation of X-ray spectra of transition metal complexes. The importance of the different parameters in the calculation is systematically investigated, laying the base for wider applications of those accurate methods in the modeling of X-ray spectroscopy. A tool to analyze the resulting spectrum in terms of molecular orbitals is also presented, strengthening the interplay between theory and experiments. With these developments and other significant ones that have happened in recent years, multiconfigurational methods can now reach new grounds and contribute to important new discoverie

Towards a Computational Ecotoxicity Assay

Thousands of anthropogenic chemicals are released into the environment each year, posing potential hazards to human and environmental health. Toxic chemicals may cause a variety of adverse health effects, triggering immediate symptoms or delayed effects over longer periods of time. It is thus crucial to develop methods that can rapidly screen and predict the toxicity of chemicals, to limit the potential harmful impacts of chemical pollutants. Computational methods are being increasingly used in toxicity predictions. Here, the method of molecular docking is assessed for screening potential toxicity of a variety of xenobiotic compounds, including pesticides, pharmaceuticals, pollutants and toxins deriving from the chemical industry. The method predicts the binding energy of the pollutants to a set of carefully selected receptors, under the assumption that toxicity in many cases is related to interference with biochemical pathways. The strength of the applied method lies in its rapid generation of interaction maps between potential toxins and the targeted enzymes, which could quickly yield molecularlevel information and insight into potential perturbation pathways, aiding in the prioritisation of chemicals for further tests. Two scoring functions are compared, Autodock Vina and the machine-learning scoring function RF-Score-VS. The results are promising, though hampered by the accuracy of the scoring functions. The strengths and weaknesses of the docking protocol are discussed, as well as future directions for improving the accuracy for the purpose of toxicity predictions.<br /