Search

We present an implementation for the use of Cholesky decomposition (CD) of two-electron integrals within the spin-free Dirac-Coulomb (SFDC) scheme that enables to perform high-accuracy coupled-cluster (CC) calculations at costs almost comparable to those of their non-relativistic counterparts. While for non-relativistic CC calculations atomic-orbital (AO) based algorithms, due to their significantly reduced disk-space requirements, are the key to efficient large-scale computations, such algorithms are less advantageous in the SFDC case due to their increased computational cost on that case. Here, molecular-orbital (MO) based algorithms exploiting the CD of the two-electron integrals allow to reduce disk-space requirements, and lead to computational cost in the CC step that are more or less the same as in the non-relativistic case. The only remaining overhead in a CD-SFDC-CC calculation are due to the need to compute additional two-electron integrals, the somewhat higher cost of the Hartree-Fock calculation in the SFDC case, and additional cost in the transformation of the Cholesky vectors from the AO to the MO representation. However, these additional costs typically amount to less than 5-15 % of the total wall time and are thus acceptable. We illustrate the efficiency of our CD scheme for SFDC-CC calculations on a series of illustrative calculations for the X(CO)$_4$ molecules with X = Ni, Pd, Pt.

We explore Riemannian optimization methods for Restricted-Open-shell Hartree-Fock (ROHF) and Complete Active Space Self-Consistent Field (CASSCF) methods. After showing that ROHF and CASSCF can be reformulated as optimization problems on so-called flag manifolds, we review Riemannian optimization basics and their application to these specific problems. We compare these methods to traditional ones and find robust convergence properties without fine-tuning of numerical parameters. Our study suggests Riemannian optimization as a valuable addition to orbital optimization for ROHF and CASSCF, warranting further investigation., Comment: 21 pages

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals., Comment: 27+8 pages, 6 figures. SI as an ancillary file

We present a new library designed to provide a simple and straightforward way to implement QM/AMOEBA and other polarizable QM/MM methods based on induced point dipoles. The library, herein referred to as OpenMMPol, is free and open-sourced and is engineered to address the increasing demand for accurate and efficient QM/MM simulations. OpenMMPol is specifically designed to allow polarizable QM/MM calculations of ground state energies and gradients, and excitation properties. Key features of OpenMMPol include a modular architecture facilitating extensibility, parallel computing capabilities for enhanced performance on modern cluster architectures, and a user-friendly interface for intuitive implementation and a simple and flexible structure for providing input data. To show the capabilities of fered by the library we present an interface with PySCF to perform QM/AMOEBA molecular dynamics, geometry optimization and excited state calculation based on (TD)DFT.