The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures

Analysis of bonding motifs in unusual molecules I: planar hexacoordinated carbon atoms

The bonding structures of CO3Li3+ and CS3Li3+ are studied by means of oriented quasi-atomic orbitals (QUAOs) to assess the possibility of these molecules being planar hexacoordinated carbon (phC) systems. CH3Li and CO32- are employed as reference molecules. It is found that the introduction of Li+ ions into the molecular environment of carbonate has a greater effect on the orbital structure of the O atoms than it does on the C atom. Partial charges computed from QUAO populations imply repulsion between the positively charged C and Li atoms in CO3Li3+. Upon the transition from CO3Li3+ to CS3Li3+, the analysis reveals that the substitution of O atoms by S atoms inverts the polarity of the carbon-chalcogen σ bond. This is linked to the difference in s- and p-fractions of the QUAOs of C and S, as element electronegativities do not explain the observed polarity of the CSσ bond. Partial charges indicate that the larger electron population on the C atom in CS3Li3+ makes C-Li attraction possible. Upon comparison with the C-Li bond in methyllithium, it is found that the C-Li covalent interactions in CO3Li3+ and CS3Li3+ have about 14% and 6% of the strength of the C-Li covalent interaction in CH3Li, respectively. Consequently, it is concluded that only CS3Li3+ may be considered to be a phC system.

Relativistic Coupled Cluster with Completely Renormalized and Perturbative Triples Corrections

We have implemented noniterative triples corrections to the energy from coupled-cluster with single and double excitations (CCSD) within the 1-electron exact two-component (1eX2C) relativistic framework. The effectiveness of both the CCSD(T) and the completely renormalized (CR) CC(2,3) approaches are demonstrated by performing all-electron computations of the potential energy curves and spectroscopic constants of copper, silver, and gold dimers in their ground electronic states. Spin-orbit coupling effects captured via the 1eX2C framework are shown to be crucial for recovering the correct shape of the potential energy curves, and the correlation effects due to triples in these systems change the dissociation energies by about 0.1-0.2 eV or about 4-7%. We also demonstrate that relativistic effects and basis set size and contraction scheme are significantly more important in Au2 than in Ag2 or Cu2.

Nonadiabatic Exciton Dynamics and Energy Gradients in the Framework of FMO-LC-TDDFTB

We introduce a novel methodology for simulating the excited-state dynamics of extensive molecular aggregates in the framework of the long-range corrected time-dependent density-functional tight-binding fragment molecular orbital method (FMO-LC-TDDFTB) combined with the mean-field Ehrenfest method. The electronic structure of the system is described in a quasi-diabatic basis composed of locally excited and charge-transfer states of all fragments. In order to carry out nonadiabatic molecular dynamics simulations, we derive and implement the excited-state gradients of the locally excited and charge-transfer states. Subsequently, the accuracy of the analytical excited-state gradients is evaluated. The applicability to the simulation of exciton transport in organic semiconductors is illustrated on a large cluster of anthracene molecules. Additionally, nonadiabatic molecular dynamics simulations of a model system of benzothieno-benzothiophene molecules highlight the method's utility in studying charge-transfer dynamics in organic materials. Our new methodology will facilitate the investigation of excitonic transfer in extensive biological systems, nanomaterials, and other complex molecular systems consisting of thousands of atoms.

A Merger of the Spin-Flip ORMAS Approach and the MC-PDFT Method

The SF-ORMAS-PDFT (spin-flip occupation restricted multiple active space-pair density functional theory) approach combines the SF-ORMAS-CI method with the MC-PDFT method to treat both static and dynamic correlation in multiconfigurational systems. The static correlation description is generated via the spin-flip approach, which uses a high-spin single reference determinant to treat excited states with multiconfigurational characters. The on-top pair density functional theory uses a translation scheme applied to GGA density functionals. The SF-ORMAS-PDFT scheme has also been combined with virtual valence orbitals (VVO), a well-defined subspace of the virtual molecular orbitals, giving rise to significant speedups relative to the use of the full virtual space. The accuracy of the SF-ORMAS-PDFT method is tested by calculating 65 vertical excitation energies of 12 small- and medium-sized organic molecules. The SF-ORMAS-PDFT vertical excitation energies calculated with VVOs are comparable to those calculated with the full virtual space. The SF-ORMAS-PDFT/6-31G(d) level of theory predicts the rotational barrier of ethylene to be 65.5 and 65.9 kcal/mol, with full virtual space and VVOs, respectively. These predicted barrier heights compare well with the experimental value of 65 kcal/mol.

Multiscale biomolecular simulations in the exascale era

The complexity of biological systems and processes, spanning molecular to macroscopic scales, necessitates the use of multiscale simulations to get a comprehensive understanding. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations are crucial for capturing processes beyond the reach of classical MD simulations. The advent of exascale computing offers unprecedented opportunities for scientific exploration, not least within life sciences, where simulations are essential to unravel intricate molecular mechanisms underlying biological processes. However, leveraging the immense computational power of exascale computing requires innovative algorithms and software designs. In this context, we discuss the current status and future prospects of multiscale biomolecular simulations on exascale supercomputers with a focus on QM/MM MD. We highlight our own efforts in developing a versatile and high-performance multiscale simulation framework with the aim of efficient utilization of state-of-the-art supercomputers. We showcase its application in uncovering complex biological mechanisms and its potential for leveraging exascale computing.

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree-Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. The accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, Rcut. Rcut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small Rcut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase

The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).

Accelerating Coupled-Cluster Calculations with GPUs: An Implementation of the Density-Fitted CCSD(T) Approach for Heterogeneous Computing Architectures Using OpenMP Directives

An algorithm is presented for the coupled-cluster singles, doubles, and perturbative triples correction [CCSD(T)] method based on the density fitting or the resolution-of-the-identity (RI) approximation for performing calculations on heterogeneous computing platforms composed of multicore CPUs and graphics processing units (GPUs). The directive-based approach to GPU offloading offered by the OpenMP application programming interface has been employed to adapt the most compute-intensive terms in the RI-CCSD amplitude equations with computational costs scaling as O ( N O 2 N V 4 ) , O ( N O 3 N V 3 ) , and O ( N O 4 N V 2 ) (where NO and NV denote the numbers of correlated occupied and virtual orbitals, respectively) and the perturbative triples correction to execute on GPU architectures. The pertinent tensor contractions are performed using an accelerated math library such as cuBLAS or hipBLAS. Optimal strategies are discussed for splitting large data arrays into tiles to fit them into the relatively small memory space of the GPUs, while also minimizing the low-bandwidth CPU-GPU data transfers. The performance of the hybrid CPU-GPU RI-CCSD(T) code is demonstrated on pre-exascale supercomputers composed of heterogeneous nodes equipped with NVIDIA Tesla V100 and A100 GPUs and on the world's first exascale supercomputer named "Frontier", the nodes of which consist of AMD MI250X GPUs. Speedups within the range 4-8× relative to the recently reported CPU-only algorithm are obtained for the GPU-offloaded terms in the RI-CCSD amplitude equations. Applications to polycyclic aromatic hydrocarbons containing 16-66 carbon atoms demonstrate that the acceleration of the hybrid CPU-GPU code for the perturbative triples correction relative to the CPU-only code increases with the molecule size, attaining a speedup of 5.7× for the largest circumovalene molecule (C66H20). The GPU-offloaded code enables the computation of the perturbative triples correction for the C60 molecule using the cc-pVDZ/aug-cc-pVTZ-RI basis sets in 7 min on Frontier when using 12,288 AMD GPUs with a parallel efficiency of 83.1%.