1
|
Drabik G, Radoń M. Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes. J Chem Theory Comput 2024; 20:3199-3217. [PMID: 38574194 PMCID: PMC11044276 DOI: 10.1021/acs.jctc.4c00092] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 03/19/2024] [Accepted: 03/20/2024] [Indexed: 04/06/2024]
Abstract
Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.
Collapse
Affiliation(s)
- Gabriela Drabik
- Jagiellonian
University, Doctoral School
of Exact and Natural Sciences, Łojasiewicza 11, 30-348 Kraków, Poland
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| | - Mariusz Radoń
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| |
Collapse
|
2
|
North SC, Wilson AK. Ab Initio Composite Approaches for Heavy Element Energetics: Ionization Potentials for the Actinide Series of Elements. J Phys Chem A 2022; 126:3027-3042. [PMID: 35427146 DOI: 10.1021/acs.jpca.2c01007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The first, second, and third gas-phase ionization potentials have been determined for the actinide series of elements using an ab initio composite scalar and fully relativistic approach, employing the coupled cluster with single, double, and perturbative triple excitations (CCSD(T)) and Dirac Hartree-Fock (DHF) methods, extrapolated to the complete basis set (CBS) limit. The impact of electron correlation and basis set choice within this framework are examined. Additionally, the first three ionization potentials were obtained using an ab initio heavy element correlation-consistent Composite Approach (here referred to as α-ccCA). This is the first utilization of a ccCA for actinide species.
Collapse
Affiliation(s)
- Sasha C North
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
3
|
Pham BQ, Datta D, Gordon MS. PDG: A Composite Method Based on the Resolution of the Identity. J Phys Chem A 2021; 125:9421-9429. [PMID: 34658243 DOI: 10.1021/acs.jpca.1c06186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller-Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles-doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the "standard" G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. The efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4-5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.
Collapse
|
4
|
Mosbach S, Menon A, Farazi F, Krdzavac N, Zhou X, Akroyd J, Kraft M. Multiscale Cross-Domain Thermochemical Knowledge-Graph. J Chem Inf Model 2020; 60:6155-6166. [PMID: 33242243 DOI: 10.1021/acs.jcim.0c01145] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In this paper, we develop a set of software agents which improve a knowledge-graph containing thermodynamic data of chemical species by means of quantum chemical calculations and error-canceling balanced reactions. The knowledge-graph represents species-associated information by making use of the principles of linked data, as employed in the Semantic Web, where concepts correspond to vertices and relationships between the concepts correspond to edges of the graph. We implement this representation by means of ontologies, which formalize the definition of concepts and their relationships, as a critical step to achieve interoperability between heterogeneous data formats and software. The agents, which conduct quantum chemical calculations and derive the estimates of standard enthalpies of formation, update the knowledge-graph with newly obtained results, improving data values, and adding nodes and connections between them. A key distinguishing feature of our approach is that it extends an existing, general-purpose knowledge-graph, called J-Park Simulator (http://theworldavatar.com), and its ecosystem of autonomous agents, thus enabling seamless cross-domain applications in wider contexts. To this end, we demonstrate how quantum calculations can directly affect the atmospheric dispersion of pollutants in an industrial emission use-case.
Collapse
Affiliation(s)
- Sebastian Mosbach
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom.,Cambridge Centre for Advanced Research and Education in Singapore (CARES), CREATE Tower #05-05, 1 Create Way, Singapore, 138602
| | - Angiras Menon
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - Feroz Farazi
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - Nenad Krdzavac
- Cambridge Centre for Advanced Research and Education in Singapore (CARES), CREATE Tower #05-05, 1 Create Way, Singapore, 138602
| | - Xiaochi Zhou
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - Jethro Akroyd
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom.,Cambridge Centre for Advanced Research and Education in Singapore (CARES), CREATE Tower #05-05, 1 Create Way, Singapore, 138602
| | - Markus Kraft
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom.,Cambridge Centre for Advanced Research and Education in Singapore (CARES), CREATE Tower #05-05, 1 Create Way, Singapore, 138602.,School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, Singapore, 637459
| |
Collapse
|
5
|
Patel P, Kuntz DM, Jones MR, Brooks BR, Wilson AK. SAMPL6 logP challenge: machine learning and quantum mechanical approaches. J Comput Aided Mol Des 2020; 34:495-510. [PMID: 32002780 PMCID: PMC10817701 DOI: 10.1007/s10822-020-00287-0] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 01/08/2020] [Indexed: 10/25/2022]
Abstract
Two different types of approaches: (a) approaches that combine quantitative structure activity relationships, quantum mechanical electronic structure methods, and machine-learning and, (b) electronic structure vertical solvation approaches, were used to predict the logP coefficients of 11 molecules as part of the SAMPL6 logP blind prediction challenge. Using electronic structures optimized with density functional theory (DFT), several molecular descriptors were calculated for each molecule, including van der Waals areas and volumes, HOMO/LUMO energies, dipole moments, polarizabilities, and electrophilic and nucleophilic superdelocalizabilities. A multilinear regression model and a partial least squares model were used to train a set of 97 molecules. As well, descriptors were generated using the molecular operating environment and used to create additional machine learning models. Electronic structure vertical solvation approaches considered include DFT and the domain-based local pair natural orbital methods combined with the solvated variant of the correlation consistent composite approach.
Collapse
Affiliation(s)
- Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824-1322, USA
| | - David M Kuntz
- Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, Denton, TX, 76203-5070, USA
| | - Michael R Jones
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, 20852-5690, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, 20852-5690, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824-1322, USA.
- Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, Denton, TX, 76203-5070, USA.
| |
Collapse
|
6
|
Patel P, Wilson AK. Domain-based local pair natural orbital methods within the correlation consistent composite approach. J Comput Chem 2019; 41:800-813. [PMID: 31891196 DOI: 10.1002/jcc.26129] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 10/29/2019] [Accepted: 12/01/2019] [Indexed: 01/15/2023]
Abstract
Ab initio composite approaches have been utilized to model and predict main group thermochemistry within 1 kcal mol-1 , on average, from well-established reliable experiments, primarily for molecules with less than 30 atoms. For molecules of increasing size and complexity, such as biomolecular complexes, composite methodologies have been limited in their application. Therefore, the domain-based local pair natural orbital (DLPNO) methods have been implemented within the correlation consistent composite approach (ccCA) framework, namely DLPNO-ccCA, to reduce the computational cost (disk space, CPU (central processing unit) time, memory) and predict energetic properties such as enthalpies of formation, noncovalent interactions, and conformation energies for organic biomolecular complexes including one of the largest molecules examined via composite strategies, within 1 kcal mol-1 , after calibration with 119 molecules and a set of linear alkanes. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, Michigan, 48824
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan, 48824
| |
Collapse
|
7
|
Feller D. Estimating the intrinsic limit of the Feller-Peterson-Dixon composite approach when applied to adiabatic ionization potentials in atoms and small molecules. J Chem Phys 2017; 147:034103. [DOI: 10.1063/1.4993625] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Minenkov Y, Wang H, Wang Z, Sarathy SM, Cavallo L. Heats of Formation of Medium-Sized Organic Compounds from Contemporary Electronic Structure Methods. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00335] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Yury Minenkov
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), KAUST
Catalysis Center (KCC), 23955-6900 Thuwal, Saudi Arabia
| | - Heng Wang
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), Clean
Combustion Research Center (CCRC), 23955-6900 Thuwal, Saudi Arabia
| | - Zhandong Wang
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), Clean
Combustion Research Center (CCRC), 23955-6900 Thuwal, Saudi Arabia
| | - S. Mani Sarathy
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), Clean
Combustion Research Center (CCRC), 23955-6900 Thuwal, Saudi Arabia
| | - Luigi Cavallo
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), KAUST
Catalysis Center (KCC), 23955-6900 Thuwal, Saudi Arabia
| |
Collapse
|
9
|
Sheng X, Qian S, Hu F. On the performance of an complete basis set extrapolation scheme tailored for the equilibrium distance applied to the helium dimer potential energy surface. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.02.016] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
10
|
Silva CDS, Pereira DH, Custodio R. G4CEP: A G4 theory modification by including pseudopotential for molecules containing first-, second- and third-row representative elements. J Chem Phys 2017; 144:204118. [PMID: 27250290 DOI: 10.1063/1.4952427] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The G4CEP composite method was developed from the respective G4 all-electron version by considering the implementation of compact effective pseudopotential (CEP). The G3/05 test set was used as reference to benchmark the adaptation by treating in this work atoms and compounds from the first and second periods of the periodic table, as well as representative elements of the third period, comprising 440 thermochemical data. G4CEP has not reached a so high level of accuracy as the G4 all-electron theory. G4CEP presented a mean absolute error around 1.09 kcal mol(-1), while the original method presents a deviation corresponding to 0.83 kcal mol(-1). The similarity of the optimized molecular geometries between G4 and G4CEP indicates that the core-electron effects and basis set adjustments may be pointed out as a significant factor responsible for the large discrepancies between the pseudopotential results and the experimental data, or even that the all-electron calculations are more efficient either in its formulation or in the cancellation of errors. When the G4CEP mean absolute error (1.09 kcal mol(-1)) is compared to 1.29 kcal mol(-1) from G3CEP, it does not seem so efficient. However, while the G3CEP uncertainty is ±4.06 kcal mol(-1), the G4CEP deviation is ±2.72 kcal mol(-1). Therefore, the G4CEP theory is considerably more reliable than any previous combination of composite theory and pseudopotential, particularly for enthalpies of formation and electron affinities.
Collapse
Affiliation(s)
- Cleuton de Souza Silva
- Instituto de Química, Universidade Estadual de Campinas, Barão Geraldo, P.O. Box 6154, 13083-970 Campinas, São Paulo, Brazil
| | - Douglas Henrique Pereira
- Departamento de Ciências Exatas e Biotecnológicas, Universidade Federal do Tocantins, Campus de Gurupi, 77410-530 Gurupi, Tocantins, Brazil
| | - Rogério Custodio
- Instituto de Química, Universidade Estadual de Campinas, Barão Geraldo, P.O. Box 6154, 13083-970 Campinas, São Paulo, Brazil
| |
Collapse
|
11
|
Vasiliu M, Peterson KA, Dixon DA. Benchmark-Quality Atomization Energies for BeH and BeH2. J Chem Theory Comput 2017; 13:649-653. [DOI: 10.1021/acs.jctc.6b01154] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Monica Vasiliu
- Department
of Chemistry, The University of Alabama, Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| | - Kirk A. Peterson
- Department
of Chemistry, Washington State University, Pullman, Washington 99164-4630 United States
| | - David A. Dixon
- Department
of Chemistry, The University of Alabama, Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| |
Collapse
|
12
|
|
13
|
The Impact of Larger Basis Sets and Explicitly Correlated Coupled Cluster Theory on the Feller–Peterson–Dixon Composite Method. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1016/bs.arcc.2016.02.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
14
|
Peterson C, Penchoff D, Wilson A. Prediction of Thermochemical Properties Across the Periodic Table. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1016/bs.arcc.2016.04.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|