This dissertation is focused around enabling large-scale quantum chemistry calculations by developing novel algorithms on new heterogeneous hardware architectures. Section 2 contains a study on assessing the importance of many-body dispersion phenomena in molecular clusters and understanding the need to incorporate the phenomena into existing quantum chemistry methods to accurately capture long-range interactions. It is known that highly accurate quantum chemistry methods are computationally expensive; therefore, there is a need for new computer algorithms that will be able to vastly improve the time to solution and enable correlated quantum chemistry calculations on large molecular systems. Sections 3 and 4 focus on the development and implementation of various algorithms that enable quantum chemistry calculations on molecules of unprecedent sizes. Specifically, Section 3 describes a novel CUDA-based algorithm that enables GPU-acceleration of the main parts of the Self-Consistent Field (SCF) cycle, resulting in the fastest SCF code to date. Section 4 introduces a novel OpenMP-based implementation of the Fock build formation, which enables GPU-accelerated calculations on all existing GPU and CPU hardware. Section 5 expands upon the OpenMP work and describes a new GPU-accelerated library named LibERI for evaluating molecular integrals – often the main computational bottleneck in many quantum chemistry methods. LibERI is intended to serve as a library for fast quantum chemistry calculations that can be interfaced with any existing quantum chemistry software packages through the MolSSI Driver Interface (MDI). MDI is a standard interface between scientific computational codes. One of LibERI’s objectives is to seamlessly integrate with a Molecular Mechanics (MM) package using MDI to enable fast QM/MM calculations. This integration would allow for efficient GPU-accelerated treatment of large QM regions via LibERI. Lastly, Section 6 concludes this dissertation with a study that applies existing quantum chemistry methods to unravel the mechanism of a Pd-catalyzed carboacylation reaction involving the activation of ester C-O bonds.
This dissertation is focused around enabling large-scale quantum chemistry calculations by developing novel algorithms on new heterogeneous hardware architectures. Section 2 contains a study on assessing the importance of many-body dispersion phenomena in molecular clusters and understanding the need to incorporate the phenomena into existing quantum chemistry methods to accurately capture long-range interactions. It is known that highly accurate quantum chemistry methods are computationally expensive; therefore, there is a need for new computer algorithms that will be able to vastly improve the time to solution and enable correlated quantum chemistry calculations on large molecular systems. Sections 3 and 4 focus on the development and implementation of various algorithms that enable quantum chemistry calculations on molecules of unprecedent sizes. Specifically, Section 3 describes a novel CUDA-based algorithm that enables GPU-acceleration of the main parts of the Self-Consistent Field (SCF) cycle, resulting in the fastest SCF code to date. Section 4 introduces a novel OpenMP-based implementation of the Fock build formation, which enables GPU-accelerated calculations on all existing GPU and CPU hardware. Section 5 expands upon the OpenMP work and describes a new GPU-accelerated library named LibERI for evaluating molecular integrals – often the main computational bottleneck in many quantum chemistry methods. LibERI is intended to serve as a library for fast quantum chemistry calculations that can be interfaced with any existing quantum chemistry software packages through the MolSSI Driver Interface (MDI). MDI is a standard interface between scientific computational codes. One of LibERI’s objectives is to seamlessly integrate with a Molecular Mechanics (MM) package using MDI to enable fast QM/MM calculations. This integration would allow for efficient GPU-accelerated treatment of large QM regions via LibERI. Lastly, Section 6 concludes this dissertation with a study that applies existing quantum chemistry methods to unravel the mechanism of a Pd-catalyzed carboacylation reaction involving the activation of ester C-O bonds.
Graphical Processing Units (GPUs) are highly parallel, programmable accelerators boasting high peak floating point performance. Over the last couple of years the use of GPUs for general purpose computing have revolutionized quantum chemistry. The computational bottleneck in an ab-initio quantum method is the calculation of a large number of twoelectron integrals. To date, a number of GPU accelerated two-electron integral implementations have been developed significantly improving the performance of a static quantum mechanical (QM) calculation. However, when performing an ab-initio QM gradient calculation, optimization, QM or Hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) dynamics simulation the twoelectron integral derivatives arise as an additional bottleneck. Hybrid QM/MM methods particularly dynamics methods are commonly used to study large chemical/biological systems. These methods are a popular choice used for studying reaction mechanisms, conformational and configurational structures important in glycobiology. Usually a semiempirical QM method is used however, these have shown variable accuracy for the study of carbohydrate conformation and prevents an analytical investigation of electronic structure. The use of a higher level of theory, such as an ab initio method, is desirable however this comes at a much greater computational cost. Together with the bottlenecks above this cost results from the polarization of the QM region within an electrostatic embedding scheme, which requires the calculation of a large number of one-electron integral derivatives. Thus for QM/MM calculations the one-electron integral derivatives becomes a third bottleneck together with the two mentioned above. Recently, the above bottlenecks have become popular GPU acceleration targets. This thesis describes an extension of the GPU based Quantum Supercharger Library (QSL) to perform the above calculations/simulations. In contrast to GPU packages developed from the ground up, the QSL is a library of routines aimed at accelerating legacy codes, such as GAMESS-UK, GAMESS-US and NWChem, used in electronic structure calculations. Algorithms are presented for accelerating the one- and two-electron integral derivatives on a GPU. In addition to the derivatives, the one-electron integral calculation was ported to the GPU in order to remove a small additional cost arising for large QM/MM systems. Furthermore, the use of automatic code generation for generating GPU kernels was explored and compared to the original approach. Using the QSL library implemented in GAMESS-UK several benchmark calculations and simulations were performed. These were performed in double precision on a single GPU (Kepler K20) and compared to a single CPU using the 6-31G basis set. A speedup of up to 9.3X is achieved for an ab initio gradient calculation compared to the CPU running an optimized serial version of GAMESS-UK using the Schlegel method. For a single point QM/MM calculation of cellobiose in different sized water spheres (3267-24843 point charges) speedups of between 13X and 34X is achieved. QSL/GAMESS-UK coupled to the CHARMM molecular dynamics package was then used in order to perform accelerated molecular dynamics simulations. Benchmark QM and QM/MM molecular dynamics simulations were performed on cellobiose in vacuo and in a water sphere (45 QM atoms and 24843 point charges respectively). The QSL is able to perform 9.7 ps/day of ab initio QM dynamics and 6.4 ps/day of QM/MM dynamics. Testing of the auto-generated version of the QSL showed better performance for lower angular momentum classes but reduced performance for higher angular momentum classes. The efficiency of the integral and derivative routines within the QSL library was tested on a computationally intense realistic glycobiological condensed phase free energy computation. Ab-initio pucker free energy surfaces/volumes of !-Ribofuranose and !-Glucopyranose in vacuum and in water were computed. These are the first converged Hartree-Fock/MM free energy simulations for these carbohydrates performed in solution. The value of the ab initio Free energy surfaces/volumes was demonstrated through an analysis of solvent polarization effects that are evident through a comparison of the vacuum and solution minimum free energy pathways. In particular the water structure around the pucker conformers, analysis of the primary alcohol distribution as well as analysis of the electronic structure of these conformers reveals that water significantly affects the free energy pathways, primary alcohol distribution and the barriers for inter-conversion between pucker conformers.
Predicting the chemico-physical properties of large molecular systems is a formidable challenge in chemistry and materials science. Traditional quantum mechanical methods, while accurate, have impractical scaling for large molecules with thousands of atoms, which are crucial in the development of novel therapeutics, catalysts, and nanomaterials. To address this, molecular fragmentation algorithms have been proposed to improve scalability and enable extensive parallelism. In this article, we introduce a significant enhancement to the Fragment Molecular Orbital (FMO) method, termed the Coulomb-Perturbed Fragmentation (CPF) method. CPF incorporates algorithmic improvements and implementation enhancements to optimize performance on heterogeneous computing systems equipped with a large number of GPUs. Key developments include a significant simplification of iteratitve self-consistent field (SCF) algorithm, advanced data management through a one-sided communication model, topology-aware optimizations, and a hybrid communication strategy for intra-group exchanges. Moreover, CPF integrates a distributed dynamic multi-layer load balancing scheme to optimise fragment distribution and workload management across nodes and GPUs. Performance evaluations on a 420-atom benzene molecule system comprising 35 fragments reveal that CPF outperforms existing GPU/CPU-based FMO algorithms in both efficiency and accuracy. When deployed on the Gadi supercomputer, CPF achieves over 97% parallel efficiency on 20 nodes, with scalability maintaining above 98% and 90% efficiency in weak-scaling tests for smaller and larger systems, respectively. Notably, CPF matches or exceeds the computational accuracy of conventional FMO methods, marking a substantial progress in the field of computational chemistry for fragmentation-based large-scale molecular modelling.
Heterogeneous nodes composed of a multicore CPU and at least one graphics processing unit (GPU) are increasingly common in high-performance scientific computing, and significant programming effort is currently being undertaken to port existing scientific algorithms to these unique architectures. We present implementations for two many-body quantum chemistry methods on heterogeneous nodes: the coupled-cluster with single and double excitations (CCSD) and time-dependent configuration interaction with single and double excitations (TD-CISD) methods. Both methods can be implemented on a computer as a series of dense matrix-matrix multiplications, operations that GPUs are particularly adept at performing. The GPU-accelerated CCSD algorithm is as much as 4.3 times faster than the corresponding CPU algorithm and 9.7 times faster than the algorithm in the Molpro package. The TD-CISD algorithm is accelerated by as much as a factor of 3.9 by GPUs. Enhanced performance is achieved by overlapping CPU and GPU computations.
The dynamic rearrangement of a molecule at finite temperature is not captured by a single static input geometry to a quantum chemical calculation. Multilevel workflows have established themselves as standard procedures to capture the dependence of the molecule’s configuration on computed properties. Each step within this workflow uses an electronic structure theory method that features the optimal computational time to accuracy ratio. For exploratory steps the computational expense of a method needs to be minimal as several thousand energy and nuclear gradient evaluations are needed. Semiempirical tight binding methods, such as GFN2-xTB are often used for such purposes. However, these methods usually lack the accuracy to rank conformers relative to each other within chemical accuracy (2 kcal/mol). Therefore, so-called “low-cost” DFT methods are commonly used for geometry optimizations and conformer ranking. The combination of (semi)empirical correction models with small basis sets leads to good accuracy at moderate computational cost. However, even these highly efficient methods reach the limits of feasibility for very large systems (> 200 atoms). This work proposes solutions to partially alleviate those restrictions and allow routine calculations of these system sizes outside of large high-performance computing facilities. The key to unlock those calculations is to combine tailored algorithmic developments with consumer-grade GPUs. A GPU-enabled linear algebra library is developed to accelerate further development and make the entry to GPU-accelerated software development easier. To demonstrate the ease of implementation, the addition of Mulliken-approximated exchange to GFN2-xTB will be shown. In most semiempirical methods, the general eigenvalue solver algorithm becomes the computational bottleneck, dominating the walltime with increasing system size. A density matrix purification algorithm is tailored to give up to six-fold acceleration of GFN2-xTB calculations using a mixed-precision approach. Lastly, 3c-methods are added to the GPU-accelerated TeraChem software package. In combination with the “mixed-precision” scheme for computing the electron repulsion integrals, this addition decreases the runtime of DFT calculations by up to an order of magnitude on mid-level consumer hardware. The overall goal of these projects is to make GPU-accelerated calculations in quantum chemistry more available. This achieve by exposing algorithmic developments that are performant on reasonably priced hardware and making GPU computing more accessible to non-experts.
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 N O and N V 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 (C 66 H 20 ). The GPU-offloaded code enables the computation of the perturbative triples correction for the C 60 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%.
The advancements in classical computational components such as CPUs and GPUs have had a reverberating impact on the democratization of computational power available, with the field of computational chemistry reaping significant benefits from the propulsion of the limits of systems studied. The accurate prediction of chemical properties and reaction dynamics through computational simulation is now a centerpiece of modern chemical research, significantly accelerating discoveries research areas of materials science and drug development, as well as industrially relevant processes such as the emission of aldehydes from soft foams, flame retardants and pyrolysis and the carbonylation of epoxides. These achievements of computational chemistry still are overshadowed in the face of fundamental bottlenecks where classical methods remain constrained by the exponential scaling when treating large and/or strongly correlated quantum many-body systems. This innate computational complexity often constitutes chemically accurate calculations based on first principles impractical, or plain impossible, for systems beyond modest sizes. Quantum computing offers a potential pathway around these limitations, originally formally formulated by David Deutsch as a tool for simulating quantum systems, with the first quantum algorithms indicating speed-ups compared to classical counterparts arising from Peter Shor and Lov Grover for the prime factorization of an integer number and unstructured search respectively. Relying on fundamental theorems by Claude Shannon, it has been demonstrated that through the use of error-correcting codes fault-tolerant quantum computers are able to safeguard information against noise and decoherence, strongly indicating towards exponential speedups in solving central computational chemistry problems, such as determining ground state energies or reaction pathways of chemical systems. Substantial strides in experimental quantum computing over the past decade, including proof-of-principle demonstrations of quantum error correction, mean to bring full-scale fault-tolerant quantum computing within the grasp of humanity. The breakthroughs needed to attain scalable, error-corrected quantum computation with capabilities exceeding those of classical methods in real-world scenarios is likely still some years away. While striving to address the engineering challenges standing in the way of achieving fault tolerant quantum computing, substantial effort has been focused toward developing and benchmarking hybrid algorithms tailored for the Noisy Intermediate- Scale Quantum (NISQ) devices, available now to researchers. These hybrid quantum-classical algorithms, with the most notable being the Variational Quantum Eigensolver (VQE), allow the probing of quantum computing’s utility in computational chemistry even in the face of significant hardware noise and limited coherence times. The work presented in this thesis introduces several key additions the VQE framework, geared especially towards to NISQ hardware, where the mindful use of available resources is imperative. First, the implementation of the Conjugate Model Search optimizer is described, an optimization method that integrates simple models that are surrogate to the original function, combined with conjugate gradient optimization to navigate effectively parameter landscapes typical of variational algorithms, prone to noise. Second, this thesis contains the proposal and extensive study of a hyperparameter tuning procedure employing genetic algorithms, allowing for comprehensive simulation benchmarks that augment the robustness and efficiency of VQE implementations with respect to convergence time through real wall-time models of quantum processors. Continuing, a bridge is presented between quantum simulation and quantum chemistry, through qubit-to-fermion mappings of existing quantum gates, further enabling the direct and accurate modeling of chemical systems on trapped atom quantum hardware with the added arsenal provided by techniques designed for qubit architectures. Additionally, experimental data contained in this thesis include detailed benchmarking data collected from state-of-the-art trapped ion quantum processors available at the Johannes Gutenberg Universität Mainz, demonstrating the actual viability and performance of various methodologies for extracting results. Lastly, the state of VQE is examined, as it served as a springboard for a plethora of advancements in the development of quantum computing hardware, as well as the direction the field of simulation of computational quantum chemistry via quantum devices is headed towards in the coming years.