1
|
Ye HZ, Berkelbach TC. Periodic Local Coupled-Cluster Theory for Insulators and Metals. J Chem Theory Comput 2024. [PMID: 39376105 DOI: 10.1021/acs.jctc.4c00936] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/09/2024]
Abstract
We describe the implementation details of periodic local coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)] using local natural orbitals (LNOs) and k-point symmetry. We discuss and compare several choices for orbital localization, fragmentation, and LNO construction. By studying diamond and lithium, we demonstrate that periodic LNO-CC theory can be applied with equal success to both insulators and metals, achieving speedups of 2 to 3 orders of magnitude even for moderately sized k-point meshes. Our final predictions of the equilibrium cohesive energy, lattice constant, and bulk modulus for diamond and lithium are in good agreement with previous theoretical predictions and experimental results.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Initiative for Computational Catalysis, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
2
|
Nagy PR. State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms. Chem Sci 2024:d4sc04755a. [PMID: 39246365 PMCID: PMC11376132 DOI: 10.1039/d4sc04755a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2024] [Accepted: 07/30/2024] [Indexed: 09/10/2024] Open
Abstract
In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.
Collapse
Affiliation(s)
- Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics Műegyetem rkp. 3. H-1111 Budapest Hungary
- HUN-REN-BME Quantum Chemistry Research Group Műegyetem rkp. 3. H-1111 Budapest Hungary
- MTA-BME Lendület Quantum Chemistry Research Group Műegyetem rkp. 3. H-1111 Budapest Hungary
| |
Collapse
|
3
|
Zhang X, Li C, Ye HZ, Berkelbach TC, Chan GKL. Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics. J Chem Phys 2024; 161:014109. [PMID: 38949583 DOI: 10.1063/5.0212274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.
Collapse
Affiliation(s)
- Xing Zhang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chenghan Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
4
|
Zheng Y, Ni Z, Wang Y, Li W, Li S. Analytical Gradient Using Cluster-in-Molecule RI-MP2 Method for the Geometry Optimizations of Large Systems. J Chem Theory Comput 2024; 20:3626-3636. [PMID: 38626287 DOI: 10.1021/acs.jctc.4c00087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
We present an efficient analytical energy gradient algorithm for the cluster-in-molecule resolution-of-identity second-order Møller-Plesset perturbation (CIM-RI-MP2) method based on the Lagrange multiplier method. Our algorithm independently constructs the Lagrangian formalism within each cluster, avoiding the solution of the coupled-perturbed Hartree-Fock (CPHF) equation for the whole system. Due to this feature, the computational cost of the CIM-RI-MP2 gradients is much lower than that of other local MP2 algorithms. Benchmark calculations of several molecules containing up to 312 atoms demonstrate the general applicability of our CIM-RI-MP2 gradient algorithm. The optimized structure of a 244-atom molecule using the CIM-RI-MP2 method with the cc-pVDZ basis set is in good agreement with the corresponding crystal structure. A single-point gradient calculation conducted for a molecular cage containing 972 atoms and 9612 basis functions takes 48 h on 25 nodes, utilizing a total of 600 CPU cores. The present CIM-RI-MP2 gradient program is applicable for obtaining the optimized geometries of large systems with hundreds of atoms.
Collapse
Affiliation(s)
- Yang Zheng
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology of Ministry of Education, Hangzhou Normal University, Hangzhou 311121, P. R. China
| | - Yuqi Wang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| |
Collapse
|
5
|
Li W, Wang Y, Ni Z, Li S. Cluster-in-Molecule Local Correlation Method for Dispersion Interactions in Large Systems and Periodic Systems. Acc Chem Res 2023; 56:3462-3474. [PMID: 37991873 DOI: 10.1021/acs.accounts.3c00538] [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/2023]
Abstract
ConspectusThe noncovalent interactions, including dispersion interactions, control the structures and stabilities of complex chemical systems, including host-guest complexes and the adsorption process of molecules on the solid surfaces. The density functional theory (DFT) with empirical dispersion correction is now the working horse in many areas of applications. Post-Hartree-Fock (post-HF) methods have been well recognized to provide more accurate descriptions in a systematic way. However, traditional post-HF methods are mainly limited to small- or medium-sized systems, and their applications to periodic condensed phase systems are still very limited due to their expensive computational costs.To extend post-HF calculations to large molecules, the cluster-in-molecule (CIM) local correlation approach has been established, allowing highly accurate electron correlation calculations that are routinely available for very large systems. In the CIM approach, the electron correlation energy of a large molecule could be obtained from electron correlation calculations on a series of clusters, each of which contains a subset of occupied and virtual localized molecular orbitals. The CIM method could be massively and efficiently parallelized on general computer clusters. The CIM method has been implemented at various electron correlation levels, including second-order Mo̷ller-Plesset perturbation theory (MP2), coupled cluster singles and doubles (CCSD), CCSD with perturbative triples correction [CCSD(T)], etc. The CIM-MP2 energy gradient algorithm was developed and applied to the geometry optimizations of large systems. The CIM method has also been extended to condensed-phase systems under periodic boundary conditions (PBC-CIM). For periodic systems, the correlation energy per unit cell could be evaluated with correlation energy contributions from a series of clusters that are built with localized Wannier functions.CIM-based electron correlation calculations have been employed to investigate a number of chemical problems in which the dispersion interaction is important. CIM-based post-HF methods including CIM domain-based local pair natural orbital (DLPNO) CCSD(T) are applied to compute the relative or binding energies of biological systems or supramolecular complexes, the reaction barrier in a relatively complex chemical reaction. The CIM-MP2 method is used to obtain the optimized geometry of large systems. CIM-based post-HF calculations have also been used to compute the cohesive energies of molecular crystals and adsorption energies of molecules on the solid surfaces. The CIM and its PBC variant are expected to become a powerful theoretical tool for accurate calculations of the energies and structures for a broad range of large systems and condensed-phase systems with significant dispersion interactions.
Collapse
Affiliation(s)
- Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Yuqi Wang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Hangzhou Normal University, Hangzhou 311121, People's Republic of China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| |
Collapse
|
6
|
Wang Y, Guo Y, Neese F, Valeev EF, Li W, Li S. Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules. J Chem Theory Comput 2023; 19:8076-8089. [PMID: 37920973 DOI: 10.1021/acs.jctc.3c00627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
7
|
Hong B, Fang T, Li W, Li S. Predicting the structures and vibrational spectra of molecular crystals containing large molecules with the generalized energy-based fragmentation approach. J Chem Phys 2023; 158:044117. [PMID: 36725497 DOI: 10.1063/5.0137072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
The generalized energy-based fragmentation (GEBF) approach under periodic boundary conditions (PBCs) has been developed to facilitate calculations of molecular crystals containing large molecules. The PBC-GEBF approach can help predict structures and properties of molecular crystals at different theory levels by performing molecular quantum chemistry calculations on a series of non-periodic subsystems constructed from the studied systems. A more rigorous formula of the forces on translational vectors of molecular crystals was proposed and implemented, enabling more reliable predictions of crystal structures. Our benchmark results on several typical molecular crystals show that the PBC-GEBF approach could reproduce the forces on atoms and the translational vectors and the optimized crystal structures from the corresponding conventional periodic methods. The improved PBC-GEBF approach is then applied to predict the crystal structures and vibrational spectra of two molecular crystals containing large molecules. The PBC-GEBF approach can provide a satisfactory description on the crystal structure of a molecular crystal containing 312 atoms in a unit cell at density-fitting second-order Møller-Plesset perturbation theory and density functional theory (DFT) levels and the infrared vibrational spectra of another molecular crystal containing 864 atoms in a unit cell at the DFT level. The PBC-GEBF approach is expected to be a promising theoretical tool for electronic structure calculations on molecular crystals containing large molecules.
Collapse
Affiliation(s)
- Benkun Hong
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People's Republic of China
| | - Tao Fang
- Genesys Microelectronics (Shanghai) Co., Ltd., 6th Floor, 11th Building, No. 3000 LongDong Road, Pu Dong District, Shanghai, People's Republic of China
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People's Republic of China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People's Republic of China
| |
Collapse
|
8
|
Wang Y, Ni Z, Neese F, Li W, Guo Y, Li S. Cluster-in-Molecule Method Combined with the Domain-Based Local Pair Natural Orbital Approach for Electron Correlation Calculations of Periodic Systems. J Chem Theory Comput 2022; 18:6510-6521. [PMID: 36240189 DOI: 10.1021/acs.jctc.2c00412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The cluster-in-molecule (CIM) method was extended to systems with periodic boundary conditions (PBCs) in a previous work (PBC-CIM) [J. Chem. Theory Comput.2019, 15, 2933], which is able to compute the electronic structures of periodic systems at second-order Møller-Plesset perturbation theory (MP2) and coupled cluster singles and doubles (CCSD) levels. However, the high computational costs of CCSD with respect to the size of clusters limit the usage of PBC-CIM to crystals with small or medium unit cells. In this work, we further develop the PBC-CIM method by employing the domain-based local pair natural orbital (DLPNO) methods for the electron correlation calculations of clusters to reduce the computational costs. The combined approach allows CCSD with perturbative triples, denoted as CCSD(T), to be computationally available for accurate descriptions of periodic systems. The distant-pair correction is also implemented to improve the accuracy of PBC-CIM. As in the molecular cases, the distant pair correction significantly improves the accuracy of various PBC-CIM methods with few additional costs. The PBC-CIM-DLPNO-CCSD(T) approach has been applied to investigate the optimized lattice parameter of the cubic LiCl crystal and two adsorption problems (CO on the NaCl(100) surface and H2O on the h-BN surface). The results show that the CIM-DLPNO-CCSD(T) method offers accurate and efficient descriptions for the studied systems. Another application to the cohesive energy of the acetic acid crystal reveals that large basis sets are necessary for reliable calculations on the cohesive energies of molecular crystals.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology of Ministry of Education, Hangzhou Normal University, Hangzhou311121, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der RuhrD-45470, Germany
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong266237, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| |
Collapse
|
9
|
Røeggen I, Gao B. On the origin of bonding in metals: lithium as a case study. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2117662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
Affiliation(s)
- Inge Røeggen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, Tromsø, Norway
- Department of Physics and Technology, UiT The Arctic University of Norway, Tromsø, Norway
| | - Bin Gao
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, Tromsø, Norway
| |
Collapse
|
10
|
Li W, Ma H, Li S, Ma J. Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning. Chem Sci 2021; 12:14987-15006. [PMID: 34909141 PMCID: PMC8612375 DOI: 10.1039/d1sc02574k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2021] [Accepted: 10/12/2021] [Indexed: 12/11/2022] Open
Abstract
Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.
Collapse
Affiliation(s)
- Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Haibo Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| |
Collapse
|
11
|
Liang Q, Yang J. Third-Order Many-Body Expansion of OSV-MP2 Wave Function for Low-Order Scaling Analytical Gradient Computation. J Chem Theory Comput 2021; 17:6841-6860. [PMID: 34704757 DOI: 10.1021/acs.jctc.1c00581] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present a many-body expansion (MBE) formulation and implementation for efficient computation of analytical energy gradients from the orbital-specific-virtual second-order Møllet-Plesset perturbation theory (OSV-MP2) based on our earlier work (Zhou et al. J. Chem. Theory Comput. 2020, 16, 196-210). The third-order MBE(3) expansion of OSV-MP2 amplitudes and density matrices was developed to adopt the orbital-specific clustering and long-range termination schemes, which avoids term-by-term differentiations of the MBE energy bodies. We achieve better efficiency by exploiting the algorithmic sparsity that allows us to prune out insignificant fitting integrals and OSV relaxations. With these approximations, the present implementation is benchmarked on a range of molecules that show an economic scaling in the linear and quadratic regimes for computing MBE(3)-OSV-MP2 amplitude and gradient equations, respectively, and yields normal accuracy comparable to the original OSV-MP2 results. The MPI-3-based parallelism through shared memory one-sided communication is further developed for improving parallel scalability and memory accessibility by sorting the MBE(3) orbital clusters into independent tasks that are distributed on multiple processes across many nodes, supporting both global and local data locations in which selected MBE(3)-OSV-MP2 intermediates of different sizes are distinguished and accordingly placed. The accuracy and efficiency level of our MBE(3)-OSV-MP2 analytical gradient implementation is finally illustrated in two applications: we show that the subtle coordination structure differences of mechanically interlocked Cu-catenane complexes can be distinguished when tuning ligand lengths; and the porphycene molecular dynamics reveals the emergence of the vibrational signature arising from softened N-H stretching associated with hydrogen transfer, using an MP2 level of electron correlation and classical nuclei for the first time.
Collapse
Affiliation(s)
- Qiujiang Liang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| | - Jun Yang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| |
Collapse
|
12
|
Demel O, Lecours MJ, Habrovský R, Nooijen M. Toward Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals. J Chem Phys 2021; 155:154104. [PMID: 34686052 DOI: 10.1063/5.0060099] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report the development of a new Laplace MP2 (second-order Møller-Plesset) implementation using a range separated Coulomb potential, partitioned into short- and long-range parts. The implementation heavily relies on the use of sparse matrix algebra, density fitting techniques for the short-range Coulomb interactions, while a Fourier transformation in spherical coordinates is used for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas orbital specific virtual orbitals associated with localized molecular orbitals are obtained from the exchange matrix associated with specific localized occupied orbitals. The range separated potential is crucial to achieve efficient treatment of the direct term in the MP2, while extensive screening is employed to reduce the expense of the exchange contribution in MP2. The focus of this paper is on controllable accuracy and linear scaling of the data entering the algorithm.
Collapse
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michael J Lecours
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Richard Habrovský
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Marcel Nooijen
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
13
|
Liu J, Cheng L. Relativistic coupled‐cluster and equation‐of‐motion coupled‐cluster methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1536] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Junzi Liu
- Department of Chemistry The Johns Hopkins University Baltimore Maryland USA
| | - Lan Cheng
- Department of Chemistry The Johns Hopkins University Baltimore Maryland USA
| |
Collapse
|
14
|
Hansen AS, Aurbakken E, Pedersen TB. Smooth potential-energy surfaces in fragmentation-based local correlation methods for periodic systems. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1896046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- A. S. Hansen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - E. Aurbakken
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - T. B. Pedersen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| |
Collapse
|
15
|
Ni Z, Guo Y, Neese F, Li W, Li S. Cluster-in-Molecule Local Correlation Method with an Accurate Distant Pair Correction for Large Systems. J Chem Theory Comput 2021; 17:756-766. [PMID: 33410327 DOI: 10.1021/acs.jctc.0c00831] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The cluster-in-molecule (CIM) local correlation approach with an accurate distant pair correlation energy correction is presented. For large systems, the inclusion of distant pair correlation energies is essential for the accurate prediction of absolute correlation energies and relative energies. Here, we propose a simple and efficient scheme for evaluating the distant pair correlation energy correction for the CIM approaches. The corrections can be readily extracted from electron correlation calculations of clusters with almost no additional effort. Benchmark calculations show that the improved CIM approach can recover more than 99.94% of the correlation energy calculated by the parent method. By combining the CIM approach with the domain-based local pair natural orbital (DLPNO) local correlation approach, we have provided accurate binding energies at the CIM-DLPNO-CCSD(T) level for a test set consisting of eight weakly bound complexes ranging in size from 200 to 1027 atoms. With these results as the reference data, the accuracy and applicability of other electron correlation methods and a few density functional methods for large systems have been assessed.
Collapse
Affiliation(s)
- Zhigang Ni
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China.,College of Materials, Chemistry and Chemical Engineering, Hangzhou Normal University, Hangzhou 311121, China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.,FAccTs GmbH, Rolandstr. 67, 50677 Köln, Germany
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| |
Collapse
|
16
|
Hermes MR, Pandharkar R, Gagliardi L. Variational Localized Active Space Self-Consistent Field Method. J Chem Theory Comput 2020; 16:4923-4937. [DOI: 10.1021/acs.jctc.0c00222] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Matthew R. Hermes
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Riddhish Pandharkar
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
17
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|