Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods.

Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.

Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory.

With a growing demand for time-domain simulations of correlated many-body systems, the development of efficient and stable integration schemes for the time-dependent Schrödinger equation is of keen interest in modern electronic structure theory. In this work, we present two approaches for the formation of the quantum propagator for time-dependent equation-of-motion coupled cluster theory based on the Chebyshev and Arnoldi expansions of the complex, nonhermitian matrix exponential, respectively. The proposed algorithms are compared with the short-iterative Lanczos method of Cooper et al. [J. Phys. Chem. A 2021 125, 5438-5447], the fourth-order Runge-Kutta method, and exact dynamics for a set of small but challenging test problems. For each of the cases studied, both of the proposed integration schemes demonstrate superior accuracy and efficiency relative to the reference simulations.

A Perspective on Sustainable Computational Chemistry Software Development and Integration.

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Using Diffusion Maps to Analyze Reaction Dynamics for a Hydrogen Combustion Benchmark Dataset.

We use local diffusion maps to assess the quality of two types of collective variables (CVs) for a recently published hydrogen combustion benchmark dataset1 that contains ab initio molecular dynamics (MD) trajectories and normal modes along minimum energy paths. This approach was recently advocated in2 for assessing CVs and analyzing reactions modeled by classical MD simulations. We report the effectiveness of this approach to molecular systems modeled by quantum ab initio MD. In addition to assessing the quality of CVs, we also use global diffusion maps to perform committor analysis as proposed in.2 We show that the committor function obtained from the global diffusion map allows us to identify transition regions of interest in several hydrogen combustion reaction channels.

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory.

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. In the present work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn-Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

A parallel, distributed memory implementation of the adaptive sampling configuration interaction method.

The many-body simulation of quantum systems is an active field of research that involves several different methods targeting various computing platforms. Many methods commonly employed, particularly coupled cluster methods, have been adapted to leverage the latest advances in modern high-performance computing. Selected configuration interaction (sCI) methods have seen extensive usage and development in recent years. However, the development of sCI methods targeting massively parallel resources has been explored only in a few research works. Here, we present a parallel, distributed memory implementation of the adaptive sampling configuration interaction approach (ASCI) for sCI. In particular, we will address the key concerns pertaining to the parallelization of the determinant search and selection, Hamiltonian formation, and the variational eigenvalue calculation for the ASCI method. Load balancing in the search step is achieved through the application of memory-efficient determinant constraints originally developed for the ASCI-PT2 method. The presented benchmarks demonstrate near optimal speedup for ASCI calculations of Cr2 (24e, 30o) with 106, 107, and 3 × 108 variational determinants on up to 16 384 CPUs. To the best of the authors' knowledge, this is the largest variational ASCI calculation to date.

Computing the Many-Body Green's Function with Adaptive Variational Quantum Dynamics.

We present a method to compute the many-body real-time Green's function using an adaptive variational quantum dynamics simulation approach. The real-time Green's function involves the time evolution of a quantum state with one additional electron with respect to the ground state wave function that is first expressed as a linear-linear combination of state vectors. The real-time evolution and the Green's function are obtained by combining the dynamics of the individual state vectors in a linear combination. The use of the adaptive protocol enables us to generate compact ansatzes on-the-fly while running the simulation. In order to improve the convergence of spectral features, Padé approximants are applied to obtain the Fourier transform of the Green's function. We demonstrate the evaluation of the Green's function on an IBM Q quantum computer. As a part of our error mitigation strategy, we develop a resolution-enhancing method that we successfully apply on the noisy data from the real-quantum hardware.

PluginPlay: Enabling exascale scientific software one module at a time.

For many computational chemistry packages, being able to efficiently and effectively scale across an exascale cluster is a heroic feat. Collective experience from the Department of Energy's Exascale Computing Project suggests that achieving exascale performance requires far more planning, design, and optimization than scaling to petascale. In many cases, entire rewrites of software are necessary to address fundamental algorithmic bottlenecks. This in turn requires a tremendous amount of resources and development time, resources that cannot reasonably be afforded by every computational science project. It thus becomes imperative that computational science transition to a more sustainable paradigm. Key to such a paradigm is modular software. While the importance of modular software is widely recognized, what is perhaps not so widely appreciated is the effort still required to leverage modular software in a sustainable manner. The present manuscript introduces PluginPlay, https://github.com/NWChemEx-Project/PluginPlay, an inversion-of-control framework designed to facilitate developing, maintaining, and sustaining modular scientific software packages. This manuscript focuses on the design aspects of PluginPlay and how they specifically influence the performance of the resulting package. Although, PluginPlay serves as the framework for the NWChemEx package, PluginPlay is not tied to NWChemEx or even computational chemistry. We thus anticipate PluginPlay to prove to be a generally useful tool for a number of computational science packages looking to transition to the exascale.

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron-Sulfur Cubanes.

Iron-sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. To provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III)4S4(SMe)4]-2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape.

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters.

The predominance of Kohn-Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn-Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

Solving Coupled Cluster Equations by the Newton Krylov Method.

We describe using the Newton Krylov method to solve the coupled cluster equation. The method uses a Krylov iterative method to compute the Newton correction to the approximate coupled cluster amplitude. The multiplication of the Jacobian with a vector, which is required in each step of a Krylov iterative method such as the Generalized Minimum Residual (GMRES) method, is carried out through a finite difference approximation, and requires an additional residual evaluation. The overall cost of the method is determined by the sum of the inner Krylov and outer Newton iterations. We discuss the termination criterion used for the inner iteration and show how to apply pre-conditioners to accelerate convergence. We will also examine the use of regularization technique to improve the stability of convergence and compare the method with the widely used direct inversion of iterative subspace (DIIS) methods through numerical examples.

Relativistic Real-Time Time-Dependent Equation-of-Motion Coupled-Cluster.

We present a relativistic time-dependent equation-of-motion coupled-cluster with single and double excitations (TD-EOM-CCSD) formalism. Unlike other explicitly time-dependent quantum chemical methods, the present approach considers the time correlation function of the dipole operator, as opposed to the expectation value of the time-dependent dipole moment. We include both scalar relativistic effects and spin-orbit coupling variationally in this scheme via the use of the exact two-component (X2C) wave function as the reference that enters the coupled-cluster formalism. In order to evaluate the accuracy of X2C-TD-EOM-CCSD, we compare zero-field splitting in atomic absorption spectra of open-shell systems (Na, K, Mg+, and Ca+) with values obtained from experiment. In closed-shell species (Na+, K+, Mg2+, and Ca2+), we observe singlet-triplet mixing in the X2C-TD-EOM-CC calculations, which results from the use of the X2C reference. The effects of the X2C reference are evaluated by comparing spectra derived from X2C-TD-EOM-CC calculations to those from TD-EOM-CC calculations using a complex generalized Hartree-Fock ([Formula: see text]-GHF) reference.

Simulating Magnetic Circular Dichroism Spectra with Real-Time Time-Dependent Density Functional Theory in Gauge Including Atomic Orbitals.

Magnetic circular dichroism (MCD) spectra are able to provide insights into the geometric, electronic, and magnetic properties of chemical systems. However, they can be challenging to understand and simulate given the need to simultaneously treat both the finite magnetic and optical fields. Thus, efficient simulations are desired to understand the spectra and resolve the molecular electronic states. Real-time dynamics are used widely in the simulation of electronic spectroscopies such as absorption as well as electronic circular dichroism, but simulating MCD with real-time dynamics is technically and theoretically challenging. In this work, we introduce a real-time dynamics-based ab initio method with a nonperturbative treatment of a static magnetic field with London orbitals for simulating the MCD spectra of closed shell systems. Effects of a magnetic field are included variationally in the spin-free nonrelativistic Hamiltonian. Real-time time-dependent density functional theory dynamics are then performed, from which we compute the response function in the presence of the external magnetic field, giving the MCD spectrum. The method developed in this paper is applied to simulate the MCD spectra for pyrimidine, pyrazine, and 1,4-naphthoquinone. Results are discussed and compared to the experiment.

Embedding non-collinear two-component electronic structure in a collinear quantum environment.

Spin-containing materials are important for spintronic applications. In this work, we present a computational framework to embed nonrelativistic, two-component calculations in a one-component environment. In this framework, both embedding scalar potential and magnetic field can be included to describe the interaction between quantum subsystems. In this current development, a generalized Kohn-Sham density functional theory electronic structure is embedded in unrestricted Kohn-Sham density functional theory. Two test systems are studied: a Li3 on a closed-shell He lattice and a Li3 on a He lattice containing a Li atom defect. The noncollinearity of Li3 is unaffected upon embedding in a closed-shell environment through the scalar potential embedding. However, the open-shell nature of the Li atom defect introduces an effective magnetic field that couples to the magnetic components of the generalized Kohn-Sham Hamiltonian. These results show that noncollinear quantum embedding in an open-shell collinear environment may modify the spin structure of the embedded system. The formalism developed herein may serve as a useful tool in the modeling of inhomogeneous magnetic fields in two-component calculations.

Approximate Green's Function Coupled Cluster Method Employing Effective Dimension Reduction.

The Green's function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structure of molecular and periodic systems, especially when electrons of the system are strongly correlated. However, in order for the GFCC to become a method that may be routinely used in the electronic structure calculations, robust numerical techniques and approximations must be employed to reduce its extremely high computational overhead. In our recent studies, it has been demonstrated that the GFCC equations can be solved directly in the frequency domain using iterative linear solvers, which can be easily distributed in a massively parallel environment. In the present work, we demonstrate a successful application of model-order-reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires high-resolution descriptions of spectral function, instead of solving the GFCC linear equation of full dimension for every single frequency point of interest, an efficiently solvable linear system model of a reduced dimension may be built upon projecting the original GFCC linear system onto a subspace. From this reduced order model is obtained a reasonable approximation to the full dimensional GFCC linear equations in both interpolative and extrapolative spectral regions. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some selected frequencies within the spectral region of interest. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The method is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems such as carbon monoxide, 1,3-butadiene, benzene, and adenine. To reach the same level of accuracy as that of the original GFCC method, the application of MOR in the GFCC method is able to significantly lower the original computational cost for the aforementioned molecules in designated frequency regimes. As a byproduct, the reduced order model obtained by this method is found to provide a high-quality initial guess, which improves the convergence rate for the existing iterative linear solver.

An ab Initio Linear Response Method for Computing Magnetic Circular Dichroism Spectra with Nonperturbative Treatment of Magnetic Field.

Magnetic circular dichroism (MCD) experiments provide a sensitive tool for exploring geometric, magnetic, and electronic properties of chemical complexes and condensed matter systems. They are also challenging to simulate because of the need to simultaneously treat the perturbations of a finite magnetic field as well as an optical field. In this work, we introduce an ab initio approach that treats the external magnetic field nonperturbatively with London orbitals for simulating the MCD spectra of closed-shell systems. Effects of a magnetic field are included variationally in the spin-free nonrelativistic Hamiltonian, followed by a linear response formalism to directly calculate the difference in absorption between the left and right circularly polarized light. In addition to the presentation of underlying mathematical formalism and implementation, the method developed in this paper has been applied to simulations of MCD spectra of the sodium anion, 2,2,6,6-tetramethylcyclohexanone, and 3-methyl-2-hexanone. Results are discussed and compared to experiments.

Generalized Hartree-Fock with Nonperturbative Treatment of Strong Magnetic Fields: Application to Molecular Spin Phase Transitions.

In this work, we present a framework of an ab initio variational approach to effectively explore electronic spin phase transitions in molecular systems inside of a homogeneous magnetic field. In order to capture this phenomenon, the complex generalized Hartree-Fock ([Formula: see text]) method is used in the spinor formalism with London orbitals. Recursive algorithms for computing the one- and two-electron integrals of London orbitals are also provided. A Pauli matrix representation of the [Formula: see text] method is introduced to separate spin contributions from the scalar part of the Fock matrix. Next, spin phase transitions in two different molecular systems are investigated in the presence of a strong magnetic field. Noncollinear spin configurations are observed during the spin phase transitions in H2 and a dichromium complex, [(H3N)4Cr(OH)2Cr(NH3)4]4+, with an increase in magnetic field strength. The competing driving forces of exchange coupling and the spin Zeeman effect have been shown to govern the spin phase transition and its transition rate. Additionally, the energetic contributions of the spin Zeeman, orbital Zeeman, and diamagnetic terms to the potential energy surface are also analyzed.

A Well-Tempered Hybrid Method for Solving Challenging Time-Dependent Density Functional Theory (TDDFT) Systems.

The time-dependent Hartree-Fock (TDHF) and time-dependent density functional theory (TDDFT) equations allow one to probe electronic resonances of a system quickly and inexpensively. However, the iterative solution of the eigenvalue problem can be challenging or impossible to converge, using standard methods such as the Davidson algorithm for spectrally dense regions in the interior of the spectrum, as are common in X-ray absorption spectroscopy (XAS). More robust solvers, such as the generalized preconditioned locally harmonic residual (GPLHR) method, can alleviate this problem, but at the expense of higher average computational cost. A hybrid method is proposed which adapts to the problem in order to maximize computational performance while providing the superior convergence of GPLHR. In addition, a modification to the GPLHR algorithm is proposed to adaptively choose the shift parameter to enforce a convergence of states above a predefined energy threshold.

Efficient Implementation of Variation after Projection Generalized Hartree-Fock.

Projected Hartree-Fock (PHF) theory can restore important symmetries to broken symmetry wave functions. Variation after projection (VAP) implementations make it possible to deliberately break and then restore a given symmetry by directly minimizing the projected energy expression. This technique can be applied to any symmetry that can be broken from relaxing constraints on single Slater determinant wave functions. For instance, generalized Hartree-Fock (GHF) wave functions are eigenfunctions of neither Sz nor S2. By relaxing these constraints, the wave function can explore a larger variational space and can reach lower energies than more constrained HF solutions. We have implemented spin-projected GHF (SGHF), which retains many of the advantages of breaking symmetry while also being a spin eigenfunction, with some notable improvements over previous implementations. Our new algorithm involves the formation of new intermediate matrices not previously discussed in the literature. Discretization of the necessary integration over the rotation group SO(3) is also accomplished much more efficiently using Lebedev grids. A novel scheme to incrementally build rotated Fock matrices is also introduced and compared with more standard approaches.