We have developed a new composite model chemistry method called WMS (Wuhan–Minnesota scaling method) with three characteristics: (1) a composite scheme to approximate the complete configuration interaction valence energy with the affordability condition of requiring no calculation more expensive than CCSD(T)/jul-cc-pV(T+d)Z, (2) low-cost methods for the inner-shell correlation contribution and scalar relativistic correction, and (3) accuracy comparable to methods with post-CCSD(T) components. The new method is shown to be accurate for the W4-17 database of 200 atomization energies with an average mean unsigned error (averaged with equal weight over strongly correlated and weakly correlated subsets of the data) of 0.45 kcal mol −1 , and the performance/cost ratio of these results compares very favorably to previously available methods. We also assess the WMS method against the DBH24-W4 database of diverse barrier heights and the energetics of the reactions of three strongly correlated Criegee intermediates with water. These results demonstrate that higher-order correlation contributions necessary to obtain high accuracy for molecular thermochemistry may be successfully extrapolated from the lower-order components of CCSD(T) calculations, and chemical accuracy can now be obtained for larger and more complex molecules and reactions. 
                        more » 
                        « less   
                    
                            
                            Scalable Ab Initio Electronic Structure Methods with Near Chemical Accuracy for Main Group Chemistry
                        
                    
    
            This study evaluates the precision of widely recognized quantum chemical methodologies, CCSD(T), DLPNO−CCSD(T), and localized ph-AFQMC, for determining the thermochemistry of main group elements. DLPNO−CCSD- (T) and localized ph-AFQMC, which o(er greater scalability compared to canonical CCSD(T), have emerged over the past decade as pivotal in producing precise benchmark chemical data. Our investigation includes closed-shell, neutral molecules, focusing on their heat of formation and atomization energy sourced from four specific small molecule data sets. First, we selected molecules from the G2 and G3 data sets, noted for their reliable experimental heat of formation data. Additionally, we incorporate molecules from the W4−11 and W4−17 sets, which provide high-level theoretical reference values for atomization energy at 0 K. Our findings reveal that both DLPNO−CCSD(T) and ph-AFQMC methods are capable of achieving a root-mean-square deviation of less than 1 kcal/mol across the combined data set, aligning with the threshold for chemical accuracy. Moreover, we make e(orts to confine the maximum deviations within 2 kcal/mol, a degree of precision that significantly broadens the applicability of these methods in fields such as biology and materials science. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2245592
- PAR ID:
- 10581704
- Publisher / Repository:
- The of Physical Chemistry A
- Date Published:
- Journal Name:
- The Journal of Physical Chemistry A
- Volume:
- 128
- Issue:
- 28
- ISSN:
- 1089-5639
- Page Range / eLocation ID:
- 5796 to 5807
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol−1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol−1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms).more » « less
- 
            Obtaining sub-chemical accuracy (1 kJ mol−1) for reaction energies of medium-sized gas-phase molecules is a longstanding challenge in the field of thermochemical modeling. The perturbative triples correction to coupled-cluster single double triple [CCSD(T)] constitutes an important component of all high-accuracy composite model chemistries that obtain this accuracy but can be a roadblock in the calculation of medium to large systems due to its O(N7) scaling, particularly in HEAT-like model chemistries that eschew separation of core and valence correlation. This study extends the work of Lesiuk [J. Chem. Phys. 156, 064103 (2022)] with new approximate methods and assesses the accuracy of five different approximations of (T) in the context of a subset of molecules selected from the W4-17 dataset. It is demonstrated that all of these approximate methods can achieve sub-0.1 kJ mol−1 accuracy with respect to canonical, density-fitted (T) contributions with a modest number of projectors. The approximation labeled Z̃T appears to offer the best trade-off between cost and accuracy and shows significant promise in an order-of-magnitude reduction in the computational cost of the CCSD(T) component of high-accuracy model chemistries.more » « less
- 
            ABSTRACT In this work, by means of quantum chemistry (Density Functional Theory (DFT), PW6B95/def2-TZVPP; DLPNO-CCSD(T)/CBS), HCN polymerization [(HCN)1 − 4] initiated and catalysed by a siloxyl radical (Si-O•) on a model silica surface is analysed. Linear HCN polymers (pHCN) are obtained by a radical initiated mechanism at a SiO• site and are characterized by a -(HC-N)- skeleton due to radical localization on the terminal N atom and radical attack on the C centre. NC heterocycles are formed by cyclization of the linear SiO-(HCN)3 − 4 and are always thermodynamically preferred over their linear counterparts, acting as thermodynamic sinks. Of particular interest to the astrochemistry community is the formation of the N-heterocycle 1,3,5-triazine that can be released into the gas phase at relatively low T (ΔG† = 23.3 kcal/mol). Full hydrogenation of SiO-(HCN•) follows two reaction channels with products: (a) SiO-CH3 + •NH2 or (b) amino-methanol + Si•, though characterized by slow kinetics. Nucleophilic addition of H2O to the electron-rich SiO-(HCN•) shows an unfavourable thermodynamics as well as a high-activation energy. The cleavage of the linear (HCN)1−4 from the SiO• site also shows a high thermodynamic energy penalty (ΔG≥82.0 kcal/mol). As a consequence, the silicate surface will be passivated by a chemically active ‘pHCN brush’ modifying the surface physico-chemical properties. The prospect of surface-catalysed HCN polymers exhibiting a high degree of chemical reactivity and proposed avenues for the formation of 1,3,5-triazine and amino-methanol opens exciting new chemical pathways to Complex Organic Matter formation in astrochemistry.more » « less
- 
            Domain‐based local pair natural orbital methods within the correlation consistent composite approachAb 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.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    