skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Domain‐based local pair natural orbital methods within the correlation consistent composite approach
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.  more » « less
Award ID(s):
1900086
PAR ID:
10458954
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Journal of Computational Chemistry
Volume:
41
Issue:
8
ISSN:
0192-8651
Page Range / eLocation ID:
p. 800-813
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Fluorocarbons have been shown experimentally by Baker and coworkers to combine with the cyclopentadienylcobalt (CpCo) moiety to form fluoroolefin and fluorocarbene complexes as well as fluorinated cobaltacyclic rings. In this connection density functional theory (DFT) studies on the cyclopentadienylcobalt fluorocarbon complexes CpCo(L)(C n F 2n ) (L = CO, PMe 3 ; n = 3 and 4) indicate structures with perfluoroolefin ligands to be the lowest energy structures followed by perfluorometallacycle structures and finally by structures with perfluorocarbene ligands. Thus, for the CpCo(L)(C 3 F 6 ) (L = CO, PMe 3 ) complexes, the perfluoropropene structure has the lowest energy, followed by the perfluorocobaltacyclobutane structure and the perfluoroisopropylidene structure less stable by 8 to 11 kcal mol −1 , and the highest energy perfluoropropylidene structure less stable by more than 12 kcal mol −1 . For the two metal carbene structures Cp(L)CoC(CF 3 ) 2 and Cp(L)CoCF(C 2 F 5 ), the former is more stable than the latter, even though the latter has Fischer carbene character. For the CpCo(L)(C 4 F 8 ) (L = CO, PMe 3 ) complexes, the perfluoroolefin complex structures have the lowest energies, followed by the perfluorometallacycle structures at 10 to 20 kcal mol −1 , and the structures with perfluorocarbene ligands at yet higher energies more than 20 kcal mol −1 above the lowest energy structure. This is consistent with the experimentally observed isomerization of the perfluorinated cobaltacyclobutane complexes CpCo(PPh 2 Me)(–CFR–CF 2 –CF 2 –) (R = F, CF 3 ) to the perfluoroolefin complexes CpCo(PPh 2 Me)(RCFCF 2 ) in the presence of catalytic quantities of HN(SO 2 CF 3 ) 2 . Further refinement of the relative energies by the state-of-the-art DLPNO-CCSD(T) method gives results essentially consistent with the DFT results summarized above. 
    more » « less
  3. Abstract Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO2binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug‐cc‐pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug‐cc‐pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug‐cc‐pwCVDZ‐PP pseudopotential. 27 HOX⋯SO2complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO2complexes in this study ranged from 1.35 to 3.81 kcal mol−1. The single‐interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol−1, while the single‐interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol−1, indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO2which may guide further research of related systems. 
    more » « less
  4. Abstract Halogen bonding (XB) has emerged as an important bonding motif in supramolecules and biological systems. Although regarded as a strong noncovalent interaction, benchmark measurements of the halogen bond energy are scarce. Here, a combined anion photoelectron spectroscopy and density functional theory (DFT) study of XB in solvated Branions is reported. The XB strength between the positively‐charged σ‐hole on the Br atom of the bromotrichloromethane (CCl3Br) molecule and the Branion was found to be 0.63 eV (14.5 kcal mol−1). In the neutral complexes, Br(CCl3Br)1,2, the attraction between the free Br atom and the negatively charged equatorial belt on the Br atom of CCl3Br, which is a second type of halogen bonding, was estimated to have interaction strengths of 0.15 eV (3.5 kcal mol−1) and 0.12 eV (2.8 kcal mol−1). 
    more » « less
  5. 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