Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, an...
Saved in:
| Published in: | Journal of chemical theory and computation Vol. 16; no. 6; p. 3757 |
|---|---|
| Main Authors: | , , , , , |
| Format: | Journal Article |
| Language: | English |
| Published: |
United States
09.06.2020
|
| ISSN: | 1549-9626, 1549-9626 |
| Online Access: | Get more information |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Abstract | By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based
molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (
) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H
O)
) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory. |
|---|---|
| AbstractList | By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory. By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical ( ) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H O) ) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory. |
| Author | Ko, Hsin-Yu Santra, Biswajit Wu, Xifan Jia, Junteng DiStasio, Jr, Robert A Car, Roberto |
| Author_xml | – sequence: 1 givenname: Hsin-Yu orcidid: 0000-0003-1619-6514 surname: Ko fullname: Ko, Hsin-Yu organization: Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States – sequence: 2 givenname: Junteng surname: Jia fullname: Jia, Junteng organization: Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States – sequence: 3 givenname: Biswajit orcidid: 0000-0003-3609-2106 surname: Santra fullname: Santra, Biswajit organization: Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States – sequence: 4 givenname: Xifan surname: Wu fullname: Wu, Xifan organization: Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States – sequence: 5 givenname: Roberto surname: Car fullname: Car, Roberto organization: Department of Physics, Princeton University, Princeton, New Jersey 08544, United States – sequence: 6 givenname: Robert A orcidid: 0000-0003-2732-194X surname: DiStasio, Jr fullname: DiStasio, Jr, Robert A organization: Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States |
| BackLink | https://www.ncbi.nlm.nih.gov/pubmed/32045232$$D View this record in MEDLINE/PubMed |
| BookMark | eNpNkE9PGzEQxS1ERQjl3lPlYw9sao_3T3xMAylIQSAVztGsdzZx5LWpvXvYT8FXZiVSqacZzftp3tObs3MfPDH2TYqFFCB_okmLo-nNQtdCyrI6Y5eyyHWmSyjP_9tnbJ7SUQilclAXbKZA5AUouGTvdx5rZ_2ebzHuKftj0BFfB9-QT9RkzwdMxO_HOtqG30432498M3jT2-DR8ZcDhTjyXxPV8FXNH7ydFP4YHJnBYeS3o8fOTjm5XJzoG75y-xBtf-huOPqGP1NsQ-zQG_rKvrToEl2f5hV73dy9rO-z7dPvh_VqmxlV5H2GslKm1tCIMl8KRCnqujVSm7JUYLRYEuR5UaFRAmDSqkYZyAtDutXQAsEV-_H59y2GvwOlftfZZMg59BSGtIPJRlZaLeWEfj-hQ91Rs3uLtsM47v6VCB8E_XcZ |
| CitedBy_id | crossref_primary_10_1039_D2CP02978B crossref_primary_10_1021_jacs_3c08226 crossref_primary_10_1063_5_0262451 crossref_primary_10_1063_5_0260892 crossref_primary_10_1016_j_scitotenv_2024_169969 crossref_primary_10_3390_molecules25133084 crossref_primary_10_1021_acs_jctc_4c00019 crossref_primary_10_1002_jcc_26816 crossref_primary_10_1557_s43579_022_00273_7 crossref_primary_10_1016_j_jpcs_2024_112020 crossref_primary_10_1088_2515_7639_adeaec crossref_primary_10_3389_fchem_2023_1232425 crossref_primary_10_1002_ejic_202100742 crossref_primary_10_1016_j_jhazmat_2024_135420 crossref_primary_10_1016_j_ssc_2025_115861 crossref_primary_10_1063_5_0225396 crossref_primary_10_1103_PhysRevB_108_165142 crossref_primary_10_1039_D5CC00230C crossref_primary_10_1063_5_0211238 crossref_primary_10_1063_5_0067781 crossref_primary_10_1063_5_0208103 crossref_primary_10_1063_5_0239984 crossref_primary_10_1021_acs_inorgchem_4c05218 crossref_primary_10_1063_5_0201698 crossref_primary_10_1016_j_jpcs_2023_111855 |
| ContentType | Journal Article |
| DBID | NPM 7X8 |
| DOI | 10.1021/acs.jctc.9b01167 |
| DatabaseName | PubMed MEDLINE - Academic |
| DatabaseTitle | PubMed MEDLINE - Academic |
| DatabaseTitleList | MEDLINE - Academic PubMed |
| Database_xml | – sequence: 1 dbid: NPM name: PubMed url: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed sourceTypes: Index Database – sequence: 2 dbid: 7X8 name: MEDLINE - Academic url: https://search.proquest.com/medline sourceTypes: Aggregation Database |
| DeliveryMethod | no_fulltext_linktorsrc |
| Discipline | Chemistry |
| EISSN | 1549-9626 |
| ExternalDocumentID | 32045232 |
| Genre | Journal Article |
| GroupedDBID | 4.4 53G 55A 5GY 5VS 7~N AABXI ABJNI ABMVS ABQRX ABUCX ACGFS ACIWK ACS ADHLV AEESW AENEX AFEFF AHGAQ ALMA_UNASSIGNED_HOLDINGS AQSVZ BAANH CS3 CUPRZ D0L DU5 EBS ED~ F5P GGK GNL IH9 J9A JG~ NPM P2P RNS ROL UI2 VF5 VG9 W1F 7X8 ABBLG ABLBI |
| ID | FETCH-LOGICAL-c354t-a173cb92d06480aa10bbfc19c6632c908e24457ac30220bb7d3c245ce9f92f2e2 |
| IEDL.DBID | 7X8 |
| ISICitedReferencesCount | 49 |
| ISICitedReferencesURI | http://www.webofscience.com/api/gateway?GWVersion=2&SrcApp=Summon&SrcAuth=ProQuest&DestLinkType=CitingArticles&DestApp=WOS_CPL&KeyUT=000541503600027&url=https%3A%2F%2Fcvtisr.summon.serialssolutions.com%2F%23%21%2Fsearch%3Fho%3Df%26include.ft.matches%3Dt%26l%3Dnull%26q%3D |
| ISSN | 1549-9626 |
| IngestDate | Sun Nov 09 11:02:17 EST 2025 Thu Jan 02 22:58:39 EST 2025 |
| IsPeerReviewed | true |
| IsScholarly | true |
| Issue | 6 |
| Language | English |
| LinkModel | DirectLink |
| MergedId | FETCHMERGED-LOGICAL-c354t-a173cb92d06480aa10bbfc19c6632c908e24457ac30220bb7d3c245ce9f92f2e2 |
| Notes | ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 23 |
| ORCID | 0000-0003-3609-2106 0000-0003-2732-194X 0000-0003-1619-6514 |
| PMID | 32045232 |
| PQID | 2354179381 |
| PQPubID | 23479 |
| ParticipantIDs | proquest_miscellaneous_2354179381 pubmed_primary_32045232 |
| PublicationCentury | 2000 |
| PublicationDate | 2020-06-09 |
| PublicationDateYYYYMMDD | 2020-06-09 |
| PublicationDate_xml | – month: 06 year: 2020 text: 2020-06-09 day: 09 |
| PublicationDecade | 2020 |
| PublicationPlace | United States |
| PublicationPlace_xml | – name: United States |
| PublicationTitle | Journal of chemical theory and computation |
| PublicationTitleAlternate | J Chem Theory Comput |
| PublicationYear | 2020 |
| SSID | ssj0033423 |
| Score | 2.525613 |
| Snippet | By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby... |
| SourceID | proquest pubmed |
| SourceType | Aggregation Database Index Database |
| StartPage | 3757 |
| Title | Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance |
| URI | https://www.ncbi.nlm.nih.gov/pubmed/32045232 https://www.proquest.com/docview/2354179381 |
| Volume | 16 |
| WOSCitedRecordID | wos000541503600027&url=https%3A%2F%2Fcvtisr.summon.serialssolutions.com%2F%23%21%2Fsearch%3Fho%3Df%26include.ft.matches%3Dt%26l%3Dnull%26q%3D |
| hasFullText | |
| inHoldings | 1 |
| isFullTextHit | |
| isPrint | |
| link | http://cvtisr.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwpV1baxQxFA7VCvqi9dJateUIPjbbnWR2s3mSddulBbssqLBvS3KStBU7U51R6K_wL3vOXOyTIPRlXiaBTJLJ9-XcPiHe4XBML8ZKBqO0zD2idC5PMtdoAnptg2nFJsxiMVmt7LIzuFVdWGV_JjYHdSiRbeSHSo9YLIsA5v31d8mqUexd7SQ07olNTVSGQ7rM6q8XQXN1u6Zeas5VKFXvpiRYO3RYDb5ijQPrG1fEvwlmAzTzJ3cd4pZ43FFMmLZ74qnYiMUz8XDWK7s9F7-POWGKMAs-chy4_ETrFGFWshhuFYNcXhC0wckNJ3PBEUe41zcwJwBs7YbQ5vPDB2oVYOrhlAOQSjjrlXbhqJW5rwaQDbrWBzD9dk6DrS-uDsAVAZa3GQsvxJf58efZieyEGSTSJ9bSZUajtyoQn5kMncuG3ifMLBJ9UWiHk0ikYWQcas7j9d4EjSofYbTJqqSi2hb3i7KILwUkYjhjp1WK1uU22Emi-06Wu2Scti5lu-JtP9drmiX2Zrgilj-r9e1s74qddsHW122FjrXmIvvEFV_9R-_X4pHiOzRbVuwbsZnot4974gH-qi-rH_vNjqLnYnn2B2PD15o |
| linkProvider | ProQuest |
| openUrl | ctx_ver=Z39.88-2004&ctx_enc=info%3Aofi%2Fenc%3AUTF-8&rfr_id=info%3Asid%2Fsummon.serialssolutions.com&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=Enabling+Large-Scale+Condensed-Phase+Hybrid+Density+Functional+Theory+Based+Ab+Initio+Molecular+Dynamics.+1.+Theory%2C+Algorithm%2C+and+Performance&rft.jtitle=Journal+of+chemical+theory+and+computation&rft.au=Ko%2C+Hsin-Yu&rft.au=Jia%2C+Junteng&rft.au=Santra%2C+Biswajit&rft.au=Wu%2C+Xifan&rft.date=2020-06-09&rft.issn=1549-9626&rft.eissn=1549-9626&rft.volume=16&rft.issue=6&rft.spage=3757&rft_id=info:doi/10.1021%2Facs.jctc.9b01167&rft.externalDBID=NO_FULL_TEXT |
| thumbnail_l | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/lc.gif&issn=1549-9626&client=summon |
| thumbnail_m | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/mc.gif&issn=1549-9626&client=summon |
| thumbnail_s | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/sc.gif&issn=1549-9626&client=summon |