1
|
Barnes TA, Ellis S, Chen J, Plimpton SJ, Nash JA. Plugin-based interoperability and ecosystem management for the MolSSI Driver Interface Project. J Chem Phys 2024; 160:214114. [PMID: 38832733 DOI: 10.1063/5.0214279] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2024] [Accepted: 05/15/2024] [Indexed: 06/05/2024] Open
Abstract
The MolSSI Driver Interface (MDI) Project is an effort to simplify and standardize the process of enabling tight interoperability between independently developed code bases and is supported by numerous software packages across the domain of chemical physics. It enables a wide variety of use cases, including quantum mechanics/molecular mechanics, advanced sampling, path integral molecular dynamics, machine learning, ab initio molecular dynamics, etc. We describe two major developments within the MDI Project that provide novel solutions to key interoperability challenges. The first of these is the development of the MDI Plugin System, which allows MDI-supporting libraries to be used as highly modular plugins, with MDI enforcing a standardized application programming interface across plugins. Codes can use these plugins without linking against them during their build process, and end-users can select which plugin(s) they wish to use at runtime. The MDI Plugin System features a sophisticated callback system that allows codes to interact with plugins on a highly granular level and represents a significant advancement toward increased modularity among scientific codes. The second major development is MDI Mechanic, an ecosystem management tool that utilizes Docker containerization to simplify the process of developing, validating, maintaining, and deploying MDI-supporting codes. Additionally, MDI Mechanic provides a framework for launching MDI simulations in which each interoperating code is executed within a separate computational environment. This eliminates the need to compile multiple production codes within a single computational environment, reducing opportunities for dependency conflicts and lowering the barrier to entry for users of MDI-enabled codes.
Collapse
Affiliation(s)
- T A Barnes
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - S Ellis
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - J Chen
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - S J Plimpton
- Temple University, Philadelphia, Pennsylvania 19122, USA
| | - J A Nash
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| |
Collapse
|
2
|
Carrasco-Busturia D, Ippoliti E, Meloni S, Rothlisberger U, Olsen JMH. Multiscale biomolecular simulations in the exascale era. Curr Opin Struct Biol 2024; 86:102821. [PMID: 38688076 DOI: 10.1016/j.sbi.2024.102821] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 04/03/2024] [Accepted: 04/04/2024] [Indexed: 05/02/2024]
Abstract
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.
Collapse
Affiliation(s)
- David Carrasco-Busturia
- DTU Chemistry, Technical University of Denmark (DTU), Kongens Lyngby, DK-2800, Denmark. https://twitter.com/@DavidCdeB
| | - Emiliano Ippoliti
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, DE-52428, Germany
| | - Simone Meloni
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie (DOCPAS), Università degli Studi di Ferrara (Unife), Ferrara, I-44121, Italy. https://twitter.com/@smeloni99
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, CH-1015, Switzerland. https://twitter.com/@lcbc_epfl
| | | |
Collapse
|
3
|
Kriebel MH, Tecmer P, Gałyńska M, Leszczyk A, Boguslawski K. Accelerating Pythonic Coupled-Cluster Implementations: A Comparison Between CPUs and GPUs. J Chem Theory Comput 2024; 20:1130-1142. [PMID: 38306601 PMCID: PMC10867805 DOI: 10.1021/acs.jctc.3c01110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 01/12/2024] [Accepted: 01/13/2024] [Indexed: 02/04/2024]
Abstract
In this work, we benchmark several Python routines for time and memory requirements to identify the optimal choice of the tensor contraction operations available. We scrutinize how to accelerate the bottleneck tensor operations of Pythonic coupled-cluster implementations in the Cholesky linear algebra domain, utilizing a NVIDIA Tesla V100S PCIe 32GB (rev 1a) graphics processing unit (GPU). The NVIDIA compute unified device architecture API interacts with CuPy, an open-source library for Python, designed as a NumPy drop-in replacement for GPUs. Due to the limitations of video memory, the GPU calculations must be performed batch-wise. Timing results of some contractions containing large tensors are presented. The CuPy implementation leads to a factor of 10-16 speed-up of the bottleneck tensor contractions compared to computations on 36 central processing unit (CPU) cores. Finally, we compare example CCSD and pCCD-LCCSD calculations performed solely on CPUs to their CPU-GPU hybrid implementation, which leads to a speed-up of a factor of 3-4 compared to the CPU-only variant.
Collapse
Affiliation(s)
- Maximilian H. Kriebel
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Paweł Tecmer
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Marta Gałyńska
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Aleksandra Leszczyk
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| |
Collapse
|
4
|
Mejia-Rodriguez D, Aprà E, Autschbach J, Bauman NP, Bylaska EJ, Govind N, Hammond JR, Kowalski K, Kunitsa A, Panyala A, Peng B, Rehr JJ, Song H, Tretiak S, Valiev M, Vila FD. NWChem: Recent and Ongoing Developments. J Chem Theory Comput 2023; 19:7077-7096. [PMID: 37458314 DOI: 10.1021/acs.jctc.3c00421] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Eric J Bylaska
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff R Hammond
- Accelerated Computing, NVIDIA Helsinki Oy, Porkkalankatu 1, 00180 Helsinki, Finland
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Zapata Computing, Inc., 100 Federal Street, Boston, Massachusetts 02110, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Marat Valiev
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
5
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
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.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
6
|
Abstract
Drug development is a wide scientific field that faces many challenges these days. Among them are extremely high development costs, long development times, and a small number of new drugs that are approved each year. New and innovative technologies are needed to solve these problems that make the drug discovery process of small molecules more time and cost efficient, and that allow previously undruggable receptor classes to be targeted, such as protein-protein interactions. Structure-based virtual screenings (SBVSs) have become a leading contender in this context. In this review, we give an introduction to the foundations of SBVSs and survey their progress in the past few years with a focus on ultralarge virtual screenings (ULVSs). We outline key principles of SBVSs, recent success stories, new screening techniques, available deep learning-based docking methods, and promising future research directions. ULVSs have an enormous potential for the development of new small-molecule drugs and are already starting to transform early-stage drug discovery.
Collapse
Affiliation(s)
- Christoph Gorgulla
- Harvard Medical School and Physics Department, Harvard University, Boston, Massachusetts, USA;
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
- Current affiliation: Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee, USA
| |
Collapse
|
7
|
Williams-Young DB, Asadchev A, Popovici DT, Clark D, Waldrop J, Windus TL, Valeev EF, de Jong WA. Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory. J Chem Phys 2023; 158:234104. [PMID: 37326157 DOI: 10.1063/5.0151070] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2023] [Accepted: 05/26/2023] [Indexed: 06/17/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- David B Williams-Young
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Andrey Asadchev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Doru Thom Popovici
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - David Clark
- NVIDIA Corporation, Santa Clara, California 95051, USA
| | - Jonathan Waldrop
- Chemical and Biological Sciences Division, Ames National Laboratory, Ames, Iowa 50011, USA
| | - Theresa L Windus
- Chemical and Biological Sciences Division, Ames National Laboratory, Ames, Iowa 50011, USA
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Wibe A de Jong
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
8
|
Raghavan B, Paulikat M, Ahmad K, Callea L, Rizzi A, Ippoliti E, Mandelli D, Bonati L, De Vivo M, Carloni P. Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations. J Chem Inf Model 2023. [PMID: 37319347 DOI: 10.1021/acs.jcim.3c00557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The initial phases of drug discovery - in silico drug design - could benefit from first principle Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics (MD) simulations in explicit solvent, yet many applications are currently limited by the short time scales that this approach can cover. Developing scalable first principle QM/MM MD interfaces fully exploiting current exascale machines - so far an unmet and crucial goal - will help overcome this problem, opening the way to the study of the thermodynamics and kinetics of ligand binding to protein with first principle accuracy. Here, taking two relevant case studies involving the interactions of ligands with rather large enzymes, we showcase the use of our recently developed massively scalable Multiscale Modeling in Computational Chemistry (MiMiC) QM/MM framework (currently using DFT to describe the QM region) to investigate reactions and ligand binding in enzymes of pharmacological relevance. We also demonstrate for the first time strong scaling of MiMiC-QM/MM MD simulations with parallel efficiency of ∼70% up to >80,000 cores. Thus, among many others, the MiMiC interface represents a promising candidate toward exascale applications by combining machine learning with statistical mechanics based algorithms tailored for exascale supercomputers.
Collapse
Affiliation(s)
- Bharath Raghavan
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Department of Physics, RWTH Aachen University, Aachen 52074, Germany
| | - Mirko Paulikat
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| | - Katya Ahmad
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| | - Lara Callea
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 1, 20126 Milan, Italy
| | - Andrea Rizzi
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Atomistic Simulations, Italian Institute of Technology, Genova 16163, Italy
| | - Emiliano Ippoliti
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| | - Davide Mandelli
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| | - Laura Bonati
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 1, 20126 Milan, Italy
| | - Marco De Vivo
- Molecular Modelling and Drug Discovery, Italian Institute of Technology, Genova 16163, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Department of Physics and Universitätsklinikum, RWTH Aachen University, Aachen 52074, Germany
| |
Collapse
|
9
|
Richard RM, Keipert K, Waldrop J, Keçeli M, Williams-Young D, Bair R, Boschen J, Crandall Z, Gasperich K, Mahmud QI, Panyala A, Valeev E, van Dam H, de Jong WA, Windus TL. PluginPlay: Enabling exascale scientific software one module at a time. J Chem Phys 2023; 158:2890211. [PMID: 37171197 DOI: 10.1063/5.0147903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Accepted: 04/18/2023] [Indexed: 05/13/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Ryan M Richard
- Ames National Laboratory, Ames, Iowa 50011, USA
- Iowa State University, Ames, Iowa 50011, USA
| | | | | | - Murat Keçeli
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | | | - Raymond Bair
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Jeffery Boschen
- Ames National Laboratory, Ames, Iowa 50011, USA
- Iowa State University, Ames, Iowa 50011, USA
| | - Zachery Crandall
- Ames National Laboratory, Ames, Iowa 50011, USA
- Iowa State University, Ames, Iowa 50011, USA
| | | | | | - Ajay Panyala
- Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | | | | | - Wibe A de Jong
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Theresa L Windus
- Ames National Laboratory, Ames, Iowa 50011, USA
- Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
10
|
Ma Y, Li Z, Chen X, Ding B, Li N, Lu T, Zhang B, Suo B, Jin Z. Machine-learning assisted scheduling optimization and its application in quantum chemical calculations. J Comput Chem 2023; 44:1174-1188. [PMID: 36648254 DOI: 10.1002/jcc.27075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/16/2022] [Accepted: 12/23/2022] [Indexed: 01/18/2023]
Abstract
Easy and effective usage of computational resources is crucial for scientific calculations, both from the perspectives of timeliness and economic efficiency. This work proposes a bi-level optimization framework to optimize the computational sequences. Machine-learning (ML) assisted static load-balancing, and different dynamic load-balancing algorithms can be integrated. Consequently, the computational and scheduling engine of the ParaEngine is developed to invoke optimized quantum chemical (QC) calculations. Illustrated benchmark calculations include high-throughput drug suit, solvent model, P38 protein, and SARS-CoV-2 systems. The results show that the usage rate of given computational resources for high throughput and large-scale fragmentation QC calculations can primarily profit, and faster accomplishing computational tasks can be expected when employing high-performance computing (HPC) clusters.
Collapse
Affiliation(s)
- Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - ZhiYing Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Xin Chen
- ShenZhen Bay Laboratory, Shenzhen, China
| | - Bowen Ding
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
| | - Ning Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
- College of Chemistry and Materials Engineering, Wenzhou University, Wen Zhou, China
| | - Teng Lu
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Baohua Zhang
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - BingBing Suo
- Department of Physics, Northwest University, Xi'an, China
| | - Zhong Jin
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
11
|
Bhat V, Callaway CP, Risko C. Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials. Chem Rev 2023. [PMID: 37141497 DOI: 10.1021/acs.chemrev.2c00704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Connor P Callaway
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| |
Collapse
|
12
|
Pathak H, Panyala A, Peng B, Bauman NP, Mutlu E, Rehr JJ, Vila FD, Kowalski K. Real-Time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-Body Methods Infrastructure. J Chem Theory Comput 2023; 19:2248-2257. [PMID: 37096369 DOI: 10.1021/acs.jctc.3c00045] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ J. Chem. Phys. 2020, 152, 174113] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams-Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photoemission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.
Collapse
Affiliation(s)
- Himadri Pathak
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
13
|
Abstract
Combustion is a reactive oxidation process that releases energy bound in chemical compounds used as fuels─energy that is needed for power generation, transportation, heating, and industrial purposes. Because of greenhouse gas and local pollutant emissions associated with fossil fuels, combustion science and applications are challenged to abandon conventional pathways and to adapt toward the demand of future carbon neutrality. For the design of efficient, low-emission processes, understanding the details of the relevant chemical transformations is essential. Comprehensive knowledge gained from decades of fossil-fuel combustion research includes general principles for establishing and validating reaction mechanisms and process models, relying on both theory and experiments with a suite of analytic monitoring and sensing techniques. Such knowledge can be advantageously applied and extended to configure, analyze, and control new systems using different, nonfossil, potentially zero-carbon fuels. Understanding the impact of combustion and its links with chemistry needs some background. The introduction therefore combines information on exemplary cultural and technological achievements using combustion and on nature and effects of combustion emissions. Subsequently, the methodology of combustion chemistry research is described. A major part is devoted to fuels, followed by a discussion of selected combustion applications, illustrating the chemical information needed for the future.
Collapse
|
14
|
Jiang A, Turney JM, Schaefer HF. Tensor Hypercontraction Form of the Perturbative Triples Energy in Coupled-Cluster Theory. J Chem Theory Comput 2023; 19:1476-1486. [PMID: 36802552 PMCID: PMC10018738 DOI: 10.1021/acs.jctc.2c00996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
Abstract
We present the working equations for a reduced-scaling method of evaluating the perturbative triples (T) energy in coupled-cluster theory, through the tensor hypercontraction (THC) of the triples amplitudes (tijkabc). Through our method, we can reduce the scaling of the (T) energy from the traditional O(N7) to a more modest O(N5). We also discuss implementation details to aid future research, development, and software realization of this method. Additionally, we show that this method yields submillihartree (mEh) differences from CCSD(T) when evaluating absolute energies and sub-0.1 kcal/mol energy differences when evaluating relative energies. Finally, we demonstrate that this method converges to the true CCSD(T) energy through the systematic increasing of the rank or eigenvalue tolerance of the orthogonal projector, as well as exhibiting sublinear to linear error growth with respect to system size.
Collapse
Affiliation(s)
- Andy Jiang
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Justin M Turney
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Henry F Schaefer
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
15
|
Unsleber JP, Liu H, Talirz L, Weymuth T, Mörchen M, Grofe A, Wecker D, Stein CJ, Panyala A, Peng B, Kowalski K, Troyer M, Reiher M. High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation. J Chem Phys 2023; 158:084803. [PMID: 36859110 DOI: 10.1063/5.0136526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.
Collapse
Affiliation(s)
- Jan P Unsleber
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Hongbin Liu
- Microsoft Quantum, Redmond, Washington 98052, USA
| | | | - Thomas Weymuth
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Adam Grofe
- Microsoft Quantum, Redmond, Washington 98052, USA
| | - Dave Wecker
- Microsoft Quantum, Redmond, Washington 98052, USA
| | - Christopher J Stein
- Department of Chemistry, TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, D-85748 Garching, Germany
| | - Ajay Panyala
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Bo Peng
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Karol Kowalski
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | | | - Markus Reiher
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
16
|
Manathunga M, Aktulga HM, Götz AW, Merz KM. Quantum Mechanics/Molecular Mechanics Simulations on NVIDIA and AMD Graphics Processing Units. J Chem Inf Model 2023; 63:711-717. [PMID: 36720086 DOI: 10.1021/acs.jcim.2c01505] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
We have ported and optimized the graphics processing unit (GPU)-accelerated QUICK and AMBER-based ab initio quantum mechanics/molecular mechanics (QM/MM) implementation on AMD GPUs. This encompasses the entire Fock matrix build and force calculation in QUICK including one-electron integrals, two-electron repulsion integrals, exchange-correlation quadrature, and linear algebra operations. General performance improvements to the QUICK GPU code are also presented. Benchmarks carried out on NVIDIA V100 and AMD MI100 cards display similar performance on both hardware for standalone HF/DFT calculations with QUICK and QM/MM molecular dynamics simulations with QUICK/AMBER. Furthermore, with respect to the QUICK/AMBER release version 21, significant speedups are observed for QM/MM molecular dynamics simulations. This significantly increases the range of scientific problems that can be addressed with open-source QM/MM software on state-of-the-art computer hardware.
Collapse
Affiliation(s)
- Madushanka Manathunga
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan48824-1322, United States
| | - Hasan Metin Aktulga
- Department of Computer Science and Engineering, Michigan State University, East Lansing, Michigan48824-1322, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California92093-0505, United States
| | - Kenneth M Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan48824-1322, United States
| |
Collapse
|
17
|
Bylaska EJ, Tratnyek PG, Torralba-Sanchez TL, Edwards KC, Dixon DA, Pignatello JJ, Xu W. Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN). J Phys Chem A 2022; 126:9059-9075. [PMID: 36417759 DOI: 10.1021/acs.jpca.2c06014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web-based tools. For each mechanism, Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH- at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than in previous studies, but still found to be ΔG‡ ≈ 18 kcal/mol using PBE0 AIMD/MM free energy simulations, suggesting that the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater (ΔG‡ > 29 kcal/mol PBE0 AIMD/MM). The most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH3)OH-, and subsequently, the methoxy anion (-OCH3) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.
Collapse
Affiliation(s)
- Eric J Bylaska
- Fundamental Sciences, Pacific Northwest National Laboratory, Richland, Washington99354, United States
| | - Paul G Tratnyek
- OHSU-PSU School of Public Health, Oregon Health & Science University, Portland, Oregon97239, United States
| | - Tifany L Torralba-Sanchez
- OHSU-PSU School of Public Health, Oregon Health & Science University, Portland, Oregon97239, United States
| | - Kyle C Edwards
- Department of Chemistry & Biochemistry, The University of Alabama, Tuscaloosa, Alabama35487-0336, United States
| | - David A Dixon
- Department of Chemistry & Biochemistry, The University of Alabama, Tuscaloosa, Alabama35487-0336, United States
| | - Joseph J Pignatello
- Department of Environmental Sciences, The Connecticut Agricultural Experiment Station, New Haven, Connecticut06511, United States
| | - Wenqing Xu
- Civil and Environmental Engineering, Villanova University, Villanova, Pennsylvania19085, United States
| |
Collapse
|
18
|
Molecular Modeling in Anion Exchange Membrane Research: A Brief Review of Recent Applications. Molecules 2022; 27:molecules27113574. [PMID: 35684512 PMCID: PMC9182285 DOI: 10.3390/molecules27113574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 05/24/2022] [Accepted: 05/30/2022] [Indexed: 12/04/2022] Open
Abstract
Anion Exchange Membrane (AEM) fuel cells have attracted growing interest, due to their encouraging advantages, including high power density and relatively low cost. AEM is a polymer matrix, which conducts hydroxide (OH−) ions, prevents physical contact of electrodes, and has positively charged head groups (mainly quaternary ammonium (QA) groups), covalently bound to the polymer backbone. The chemical instability of the quaternary ammonium (QA)-based head groups, at alkaline pH and elevated temperature, is a significant threshold in AEMFC technology. This review work aims to introduce recent studies on the chemical stability of various QA-based head groups and transportation of OH− ions in AEMFC, via modeling and simulation techniques, at different scales. It starts by introducing the fundamental theories behind AEM-based fuel-cell technology. In the main body of this review, we present selected computational studies that deal with the effects of various parameters on AEMs, via a variety of multi-length and multi-time-scale modeling and simulation methods. Such methods include electronic structure calculations via the quantum Density Functional Theory (DFT), ab initio, classical all-atom Molecular Dynamics (MD) simulations, and coarse-grained MD simulations. The explored processing and structural parameters include temperature, hydration levels, several QA-based head groups, various types of QA-based head groups and backbones, etc. Nowadays, many methods and software packages for molecular and materials modeling are available. Applications of such methods may help to understand the transportation mechanisms of OH− ions, the chemical stability of functional head groups, and many other relevant properties, leading to a performance-based molecular and structure design as well as, ultimately, improved AEM-based fuel cell performances. This contribution aims to introduce those molecular modeling methods and their recent applications to the AEM-based fuel cells research community.
Collapse
|
19
|
|
20
|
Lesiuk M. Quintic-scaling rank-reduced coupled cluster theory with single and double excitations. J Chem Phys 2022; 156:064103. [DOI: 10.1063/5.0071916] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Waldrop JM, Windus TL, Govind N. Projector-Based Quantum Embedding for Molecular Systems: An Investigation of Three Partitioning Approaches. J Phys Chem A 2021; 125:6384-6393. [PMID: 34260852 DOI: 10.1021/acs.jpca.1c03821] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Projector-based embedding is a relatively recent addition to the collection of methods that seek to utilize chemical locality to provide improved computational efficiency. This work considers the interactions between the different proposed procedures for this method and their effects on the accuracy of the results. The interplay between the embedded background, projector type, partitioning scheme, and level of atomic orbital (AO) truncation are investigated on a selection of reactions from the literature. The Huzinaga projection approach proves to be more reliable than the level-shift projection when paired with other procedural options. Active subsystem partitioning from the subsystem projected AO decomposition (SPADE) procedure proves slightly better than the combination of Pipek-Mezey localization and Mulliken population screening (PMM). Along with these two options, a new partitioning criteria is proposed based on subsystem von Neumann entropy and the related subsystem orbital occupancy. This new method overlaps with the previous PMM method, but the screening process is computationally simpler. Finally, AO truncation proves to be a robust option for the tested systems when paired with the Huzinaga projection, with satisfactory results being acquired at even the most severe truncation level.
Collapse
Affiliation(s)
| | - Theresa L Windus
- Ames Laboratory, Ames, Iowa 50011, United States.,Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
22
|
Manathunga M, Jin C, Cruzeiro VWD, Miao Y, Mu D, Arumugam K, Keipert K, Aktulga HM, Merz KM, Götz AW. Harnessing the Power of Multi-GPU Acceleration into the Quantum Interaction Computational Kernel Program. J Chem Theory Comput 2021; 17:3955-3966. [PMID: 34062061 DOI: 10.1021/acs.jctc.1c00145] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report a new multi-GPU capable ab initio Hartree-Fock/density functional theory implementation integrated into the open source QUantum Interaction Computational Kernel (QUICK) program. Details on the load balancing algorithms for electron repulsion integrals and exchange correlation quadrature across multiple GPUs are described. Benchmarking studies carried out on up to four GPU nodes, each containing four NVIDIA V100-SXM2 type GPUs demonstrate that our implementation is capable of achieving excellent load balancing and high parallel efficiency. For representative medium to large size protein/organic molecular systems, the observed parallel efficiencies remained above 82% for the Kohn-Sham matrix formation and above 90% for nuclear gradient calculations. The accelerations on NVIDIA A100, P100, and K80 platforms also have realized parallel efficiencies higher than 68% in all tested cases, paving the way for large-scale ab initio electronic structure calculations with QUICK.
Collapse
Affiliation(s)
- Madushanka Manathunga
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Chi Jin
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Vinícius Wilian D Cruzeiro
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093-0505, United States.,Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Yipu Miao
- Facebook, 1 Hacker Way, Menlo Park, California 94025, United States
| | - Dawei Mu
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, 1205 W Clark Street, Urbana, Illinois 61801, United States
| | - Kamesh Arumugam
- NVIDIA Corporation, Santa Clara, California 95051, United States
| | | | - Hasan Metin Aktulga
- Department of Computer Science and Engineering, Michigan State University, 428 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Kenneth M Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093-0505, United States
| |
Collapse
|
23
|
Hu W, Chen M. Editorial: Advances in Density Functional Theory and Beyond for Computational Chemistry. Front Chem 2021; 9:705762. [PMID: 34322476 PMCID: PMC8311291 DOI: 10.3389/fchem.2021.705762] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 06/25/2021] [Indexed: 11/30/2022] Open
Affiliation(s)
- Wei Hu
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, China
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
| | - Mohan Chen
- HEDPS, CAPT, College of Engineering, Peking University, Beijing, China
| |
Collapse
|