1
|
Gray M, Herbert JM. Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes. J Chem Phys 2024; 161:054114. [PMID: 39105555 PMCID: PMC11305816 DOI: 10.1063/5.0206533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2024] [Accepted: 07/18/2024] [Indexed: 08/07/2024] Open
Abstract
The titular domain-based local pair natural orbital (DLPNO) approximation is the most widely used method for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies in large supramolecular complexes has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density and involve nonlocal dispersion that is discarded or approximated if the screening of pair natural orbitals (PNOs) is too aggressive. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases. Here, we test the DLPNO approximation at the level of second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] for a variety of large supramolecular complexes. DLPNO-MP2 interaction energies are within 3% of canonical values for small dimers with ≲10 heavy atoms, but for larger systems, the DLPNO approximation is often quite poor unless the results are extrapolated to the canonical limit where the threshold for discarding PNOs is taken to zero. Counterpoise correction proves to be essential in reducing errors with respect to canonical results. For a sequence of nanoscale graphene dimers up to (C96H24)2, extrapolated DLPNO-MP2 interaction energies agree with canonical values to within 1%, independent of system size, provided that the basis set does not contain diffuse functions; these cause the DLPNO approximation to behave erratically, such that results cannot be extrapolated in a meaningful way. DLPNO-CCSD(T) calculations are typically performed using looser PNO thresholds as compared to DLPNO-MP2, but this significantly impacts accuracy for large supramolecular complexes. Standard DLPNO-CCSD(T) settings afford errors of 2-6 kcal/mol for dimers involving coronene (C24H12) and circumcoronene (C54H18), even at the DLPNO-CCSD(T1) level.
Collapse
Affiliation(s)
- Montgomery Gray
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
2
|
Yang MY, Wu XP. Level-Shifted Embedded Cluster Method for Modeling the Chemistry of Metal Oxides. J Chem Theory Comput 2024. [PMID: 38300767 DOI: 10.1021/acs.jctc.3c01123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2024]
Abstract
The embedded cluster method has been used extensively in the study of the chemical and physical properties of metal oxides. This method has been a popular tool due to its relatively high accuracy and low computational cost. An even more promising option may entail integrating the embedded cluster method with the combined quantum mechanical and molecular mechanical (QM/MM) approach, thereby enabling further consideration of interactions within the entire system for superior results. We aim to accurately model the chemistry of metal oxides using this combined scheme. Here, using the prototypical MgO(100) surface as a test system, with Mg9O14 as the cluster in the quantum mechanical region, we show that the embedded cluster with untailored boundary effective core potentials (ECPs) can have frontier orbital energy levels that substantially deviate from the quantum mechanical reference results. This occurs even when Mg9O9, which retains the stoichiometry of MgO, is used as the cluster in the quantum mechanical region. As a result, the chemical properties of the embedded cluster models differ from those of the quantum mechanical reference model. To address this issue, we propose a new variant of the embedded cluster method called the level-shifted embedded cluster (LSEC) method, which allows the energy levels to be shifted to match the reference levels by tuning the boundary ECPs. Our validation calculations on the adsorption of various adsorbates with different properties on the MgO(100) surface show that the overall performance of QM/MM with the LSEC method is excellent for the adsorption energies, geometries, and charge properties. The excellent performance holds for both the nonstoichiometric and stoichiometric clusters (i.e., Mg9O14 and Mg9O9, respectively), demonstrating the robustness of the LSEC method. We expect that the LSEC method can be combined with QM/MM or used separately for future chemical studies of metal oxides and other ionically bonded systems.
Collapse
Affiliation(s)
- Ming-Yu Yang
- State Key Laboratory of Green Chemical Engineering and Industrial Catalysis, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, P.R. China
| | - Xin-Ping Wu
- State Key Laboratory of Green Chemical Engineering and Industrial Catalysis, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, P.R. China
| |
Collapse
|
3
|
Aghababaei F, Nemati-Kande E. PES and transport properties of the He⋯HBr complex from kinetic theory and molecular dynamics simulations. Phys Chem Chem Phys 2024; 26:4724-4735. [PMID: 38251971 DOI: 10.1039/d3cp05145e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2024]
Abstract
The ab initio intermolecular potential of the He⋯HBr van der Waals (vdW) complex was calculated at the CCSD(T)/a5zBF level of theory and expanded in terms of the orthogonal Legendre polynomials. The PES then was implemented to calculate the interaction viscosity (η12) and diffusion (D12) coefficients through classical Mason-Monchick approximation (MMA), quantum mechanical close-coupling (CC), and molecular dynamics (MD) simulations. Energy-dependent Senftleben-Beenakker (SB) cross-sections were calculated using rotationally averaged cross-sections, and Boltzmann averaging was used to reveal the temperature dependence of the SB cross-sections over the temperature range of T = 50-1000 K. The calculated transport properties from MMA are in close agreement with the CC results, especially for temperatures lower than T = 900 K. The ab initio potential data then were used to derive LJ (12,6) and Vashishta MD force fields, and the equilibrium MD simulation methods were implemented to extract η12 and D12 coefficients using Einstein formulas. It was found that the Vashishta 3-body interaction potential model shows better accuracy than the LJ (12,6) model in MD simulations of D12. For η12, however, both MD potential models are successful, and an average absolute deviation of lower than 1% was obtained when compared to the quantum mechanical CC method.
Collapse
Affiliation(s)
- Fatemeh Aghababaei
- Department of Physical Chemistry, Chemistry Faculty, Urmia University, Urmia, Iran.
| | - Ebrahim Nemati-Kande
- Department of Physical Chemistry, Chemistry Faculty, Urmia University, Urmia, Iran.
| |
Collapse
|
4
|
Jensen F. Basis Set Superposition Errors Are Partly Basis Set Imbalances. J Chem Theory Comput 2024; 20:767-774. [PMID: 38174405 DOI: 10.1021/acs.jctc.3c01156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
When calculating fragment interaction energies by electronic structure methods employing medium-sized atom-centered basis sets, it is often observed that the effect is systematically overestimated. The common interpretation is that the systematic error arises because the basis set for the complex is more complete than for the isolated fragments, and this is denoted basis set superposition errors. It has been observed, however, that the interaction energy in some cases is underestimated, which defies the interpretation in terms of basis set completeness, and instead suggests that the effect partly is due to basis set imbalance. The imbalance can be removed by explicit optimization of the basis sets for each structure, and it is shown that this to a significant extent reduces the systematic overestimation attributed to basis set superposition error.
Collapse
Affiliation(s)
- Frank Jensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus, Denmark
| |
Collapse
|
5
|
Ma T, Huang H, Liu Y, Peng Y. Theoretical investigation on a simple turn on fluorescent probe for detection of biothiols based on coumarin unit. Front Chem 2023; 11:1290745. [PMID: 38025079 PMCID: PMC10663294 DOI: 10.3389/fchem.2023.1290745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Accepted: 10/16/2023] [Indexed: 12/01/2023] Open
Abstract
The discovery of a simple and efficient detection method for biothiols would be scientifically significant due to the crucial role of them in various physiological processes. Recently, a simple fluorescent probe, DEMCA-NBSC, based on coumarin fragments, was developed by Ding et al., and provided an efficient way for real-time sensing of biothiols both in vivo and vitro. Theoretical insights to the fluorescence sensing mechanism of the probe were provided in this work. Details of the electron transfer process in the probe under optical excitation and the fluorescent character of the probe were analyzed using a quantum mechanical method. All these theoretical results could inspire the development of a highly convenient and efficient fluorescent probe to sense biothiols both in vivo and vitro.
Collapse
Affiliation(s)
- Tianhao Ma
- Affiliated 3rd Hospital, Jinzhou Medical University, Jinzhou, China
| | - He Huang
- College of Bio-Informational Engineering, Jinzhou Medical University, Jinzhou, China
| | - Yuling Liu
- College of Bio-Informational Engineering, Jinzhou Medical University, Jinzhou, China
| | - Yongjin Peng
- College of Bio-Informational Engineering, Jinzhou Medical University, Jinzhou, China
| |
Collapse
|
6
|
Pitteloud Q, Wind P, Jensen SR, Frediani L, Jensen F. Quantifying Intramolecular Basis Set Superposition Errors. J Chem Theory Comput 2023; 19:5863-5871. [PMID: 37595013 DOI: 10.1021/acs.jctc.3c00693] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/20/2023]
Abstract
We show that medium-sized Gaussian basis sets lead to significant intramolecular basis set superposition errors at Hartree-Fock and density functional levels of theory, with artificial stabilization of compact over extended conformations for a 186 atom deca-peptide. Errors of ∼80 and ∼10 kJ/mol are observed, with polarized double zeta and polarized triple zeta quality basis sets, respectively. Two different procedures for taking the basis set superposition error into account are tested. While both reduce the error, it appears that polarized quadruple zeta basis sets are required to reduce the error below a few kJ/mol. Alternatively, the basis set superposition error can be eliminated using multiresolution methods based on Multiwavelets.
Collapse
Affiliation(s)
- Quentin Pitteloud
- Hylleraas Centre, Department of Chemistry, UiT the Arctic University of Norway, Tromsø N-9037, Norway
| | - Peter Wind
- Hylleraas Centre, Department of Chemistry, UiT the Arctic University of Norway, Tromsø N-9037, Norway
| | - Stig Rune Jensen
- Hylleraas Centre, Department of Chemistry, UiT the Arctic University of Norway, Tromsø N-9037, Norway
| | - Luca Frediani
- Hylleraas Centre, Department of Chemistry, UiT the Arctic University of Norway, Tromsø N-9037, Norway
| | - Frank Jensen
- Department of Chemistry, Aarhus University, Aarhus DK-8000, Denmark
| |
Collapse
|
7
|
Novikov AS, Bolotin DS. Xenon Derivatives as Aerogen Bond-Donating Catalysts for Organic Transformations: A Theoretical Study on the Metaphorical "Spherical Cow in a Vacuum" Provides Insights into Noncovalent Organocatalysis. J Org Chem 2023; 88:1936-1944. [PMID: 35679603 DOI: 10.1021/acs.joc.2c00680] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Computations indicate that cationic and noncharged xenon derivatives should exhibit higher catalytic activity than their iodine-based noncovalent organocatalytic congeners. Perfluorophenyl xenonium(II) is expected to demonstrate the best balance between catalytic activity and chemical stability for use in organocatalysis. Comparing its catalytic activity with that of isoelectronic perfluoroiodobenzene indicates that the high catalytic activity of cationic noncovalent organocatalysts is predominantly attributed to the electrostatic interactions with the reaction substrates, which cause the polarization of ligated species during the reaction progress. In contrast, the electron transfer and covalent contributions to the bonding between the catalyst and substrate have negligible effects. The dominant effect of electrostatic interactions results in a strong negative correlation between the calculated Gibbs free energies of activation for the modeled reactions and the highest potentials of the σ-holes on the central atoms of the catalysts. No such correlation is observed for noncharged catalysts.
Collapse
Affiliation(s)
- Alexander S Novikov
- Institute of Chemistry, Saint Petersburg State University, Universitetskaya Nab. 7/9, Saint Petersburg 199034, Russian Federation.,Peoples' Friendship University of Russia (RUDN University), Moscow 117198, Russian Federation
| | - Dmitrii S Bolotin
- Institute of Chemistry, Saint Petersburg State University, Universitetskaya Nab. 7/9, Saint Petersburg 199034, Russian Federation
| |
Collapse
|
8
|
Rather IA, Riaz U, Ali R. Experimental and Computational Anion Binding Studies of meso-Substituted One-Walled Phthalimide-based Calix[4]pyrrole. J Mol Struct 2023. [DOI: 10.1016/j.molstruc.2023.135065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
9
|
Gray M, Bowling PE, Herbert JM. Systematic Evaluation of Counterpoise Correction in Density Functional Theory. J Chem Theory Comput 2022; 18:6742-6756. [PMID: 36251499 DOI: 10.1021/acs.jctc.2c00883] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A widespread belief persists that the Boys-Bernardi function counterpoise (CP) procedure "overcorrects" supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ∼150 atoms and ligand-protein models with ∼300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. This is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP "overcorrects" for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.
Collapse
Affiliation(s)
- Montgomery Gray
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| | - Paige E Bowling
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States.,Biophysics Graduate Program, The Ohio State University, Columbus, Ohio43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States.,Biophysics Graduate Program, The Ohio State University, Columbus, Ohio43210, United States
| |
Collapse
|
10
|
Mehta N, Martin JML. Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12. J Phys Chem Lett 2022; 13:9332-9338. [PMID: 36178852 PMCID: PMC9575149 DOI: 10.1021/acs.jpclett.2c02620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 09/26/2022] [Indexed: 06/16/2023]
Abstract
Following earlier work [Mehta, N.; Martin, J. M. L. J. Chem. Theory Comput.2022, 10.1021/acs.jctc.2c00426] that showed how the slow basis set convergence of the double hybrid density functional theory can be obviated by the use of F12 explicit correlation in the GLPT2 step (second order Görling-Levy perturbation theory), we demonstrate here for the very large and chemically diverse GMTKN55 benchmark suite that the CPU time scaling of this step can be reduced (asymptotically linearized) using the localized pair natural orbital (PNO-L) approximation at negligible cost in accuracy.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, Reḥovot7610001, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, Reḥovot7610001, Israel
| |
Collapse
|
11
|
Mehta N, Martin JML. Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database. J Chem Theory Comput 2022; 18:5978-5991. [PMID: 36099641 PMCID: PMC9558368 DOI: 10.1021/acs.jctc.2c00426] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Indexed: 11/28/2022]
Abstract
Double-hybrid density functional theory (DHDFT) offers a pathway to accuracy approaching composite wavefunction approaches such as G4 theory. However, the Görling-Levy second-order perturbation theory (GLPT2) term causes them to partially inherit the slow ∝L-3 (with L the maximum angular momentum) basis set convergence of correlated wavefunction methods. This could potentially be remedied by introducing F12 explicit correlation: we investigate the basis set convergence of both DHDFT and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large and chemically diverse general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) benchmark suite. The B2GP-PLYP-D3(BJ) and revDSD-PBEP86-D4 DHDFs are investigated as test cases, together with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets as large as cc-pVQZ-F12. We show that F12 greatly accelerates basis set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12 basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPPD in orbital-based approaches, and in fact comparable in quality to cc-pV(5+d)Z. Somewhat surprisingly, aug-cc-pVDZ-F12 is not required even for the anionic subsets. In conclusion, DHDF-F12/VDZ-F12 eliminates concerns about basis set convergence in both the development and applications of double-hybrid functionals. Mass storage and I/O bottlenecks for larger systems can be circumvented by localized pair natural orbital approximations, which also exhibit much gentler system size scaling.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, 7610001 Reḥovot, Israel
| |
Collapse
|
12
|
Karton A. π–π interactions between benzene and graphene by means of large-scale DFT-D4 calculations. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111606] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
13
|
Nagy PR, Gyevi-Nagy L, Lőrincz BD, Kállay M. Pursuing the basis set limit of CCSD(T) non-covalent interaction energies for medium-sized complexes: case study on the S66 compilation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2109526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Péter R. Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - László Gyevi-Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Balázs D. Lőrincz
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Mihály Kállay
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| |
Collapse
|
14
|
Brémond E, Li H, Sancho-García JC, Adamo C. Double Hybrids and Noncovalent Interactions: How Far Can We Go? J Phys Chem A 2022; 126:2590-2599. [PMID: 35438491 DOI: 10.1021/acs.jpca.2c01193] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate evaluation of weak noncovalent interactions in large, that is those containing up to thousand atoms, molecular systems represents a difficult challenge for any quantum chemical method. Indeed, some approximations are often introduced to render affordable these calculations. Here, we consider the PBE-QIDH/DH-SVPD protocol, combining a nonempirical double hybrid functional (PBE-QIDH) with a small basis set (DH-SVPD) tailored for noncovalent interactions with a double aim: (i) explore the robustness and accuracy of this protocol with respect to other Density Functional Approximations; (ii) illustrate how its performances are affected by the computational parameters underlying the calculation of the exact exchange and the Coulomb contribution, as well as the perturbative term. To this end, we consider three data sets, namely S66, L7, and CiM13, incorporating molecules of increasing size. On the bright side, our results suggest that the PBE-QIDH/DH-SVPD protocol is particularly accurate for large systems such as those contained in the CiM13 set (up to more than 1000 atoms and 14 000 basis functions), for which the DLPNO approximation leads to a significant speed-up for the evaluation of the perturbative correlation term. However, our analysis also points out the limit of this statistical exercise, when the quality of the reference data cannot be easily assessed, due to the size of the molecular complexes involved, and when the number of molecules is limited.
Collapse
Affiliation(s)
- Eric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | - Hanwei Li
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Health and Life Sciences, F-75005 Paris, France
| | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Health and Life Sciences, F-75005 Paris, France.,Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
15
|
Matveeva R, Falck Erichsen M, Koch H, Høyvik IM. The effect of midbond functions on interaction energies computed using MP2 and CCSD(T). J Comput Chem 2022; 43:121-131. [PMID: 34738658 DOI: 10.1002/jcc.26777] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2021] [Revised: 10/08/2021] [Accepted: 10/12/2021] [Indexed: 11/06/2022]
Abstract
In this article we use MP2 and CCSD(T) calculations for the A24 and S66 data sets to explore how midbond functions can be used to generate cost effective counterpoise corrected supramolecular interaction energies of noncovalent complexes. We use the A24 data set to show that the primary role of midbond functions is not to approach the complete basis set limit, but rather to ensure a balanced description of the molecules and the interaction region (unrelated to the basis set superposition error). The need for balance is a consequence of using atom centered basis sets. In the complete basis set limit, the error will disappear, but reaching the complete basis set limit is not feasible beyond small systems. For S66 we investigate the need for increasing the number of midbond centers. Results show that adding a second midbond center increases the accuracy, but the effect is secondary to changing the atom centered basis set. Further, by comparing calculations using the 3s3p2d1f1g midbond set with using aug-cc-pVDZ and aug-cc-pVTZ as midbond sets, we see that the requirements for the midbond set to be effective, is not just that it contains diffuse functions, but also that high angular momentum functions are included. By comparing two approaches for placing midbond centers we show that results are not particularly sensitive to placement as long as the placement is reasonable.
Collapse
Affiliation(s)
- Regina Matveeva
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Trondheim, Norway
| | - Merete Falck Erichsen
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Trondheim, Norway.,Scuola Normale Superiore, Pisa, Italy
| | - Ida-Marie Høyvik
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Trondheim, Norway
| |
Collapse
|
16
|
Santra G, Semidalas E, Mehta N, Karton A, Martin JML. S66x8 noncovalent interactions revisited: new benchmark and performance of composite localized coupled-cluster methods. Phys Chem Chem Phys 2022; 24:25555-25570. [DOI: 10.1039/d2cp03938a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The S66x8 noncovalent interactions benchmark has been re-evaluated at the “sterling silver” level. Against this, a selection of computationally more economical alternatives has been assayed, ranging from localized CC to double hybrids and SAPT(DFT).
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
- School of Science and Technology, University of New England, Armidale, NSW 2351, Australia
| | - Jan M. L. Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
17
|
Scheiner S. Maximal occupation by bases of π-hole bands surrounding linear molecules. J Comput Chem 2021; 43:319-330. [PMID: 34859910 DOI: 10.1002/jcc.26792] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Revised: 11/22/2021] [Accepted: 11/23/2021] [Indexed: 12/19/2022]
Abstract
Linear molecules such as CO2 contain a positive π-hole ring that surrounds C on the molecule's equator. Quantum calculations examine the question as to how many bases can simultaneously bind to this ring. Linear molecules examined are TO2 , where T = C, Si, Ge, Sn; bases are NCH and NH3 . CO2 engages in the weakest of the tetrel bonds, and can bind up to three NCH and two NH3 . Unlike σ-hole tetrel bonds, Si forms the strongest tetrel bonds, with interaction energies as high as 43 kcal/mol with NH3 . But like GeO2 , SiO2 can sustain only two bases in its equatorial ring. The π-hole ring of SnO2 can engage in up to four tetrel bonds with either NCH or NH3 , even though these bonds are weaker than those with GeO2 or SiO2 . As all of these complexes cast TO2 in the role of multiple electron acceptor, the resulting negative cooperativity makes each successive bond weaker than its predecessor as bases are added, as well as reducing the magnitude of the central molecule's π-hole.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah, USA
| |
Collapse
|
18
|
Hill JG, Shaw RA. Correlation consistent basis sets for explicitly correlated wavefunctions: Pseudopotential-based basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements. J Chem Phys 2021; 155:174113. [PMID: 34742216 DOI: 10.1063/5.0070638] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
New correlation consistent basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements have been developed specifically for use in explicitly correlated F12 calculations. This includes orbital basis sets for valence only (cc-pVnZ-PP-F12, n = D, T, Q) and outer core-valence (cc-pCVnZ-PP-F12) correlation, along with both of these augmented with additional high angular momentum diffuse functions. Matching auxiliary basis sets required for density fitting and resolution-of-the-identity approaches to conventional and F12 integrals have also been optimized. All of the basis sets are to be used in conjunction with small-core relativistic pseudopotentials [Figgen et al., Chem. Phys. 311, 227 (2005)]. The accuracy of the basis sets is determined through benchmark calculation at the explicitly correlated coupled-cluster level of theory for various properties of atoms and diatomic molecules. The convergence of the properties with respect to the basis set is dramatically improved compared to conventional coupled-cluster calculations, with cc-pVTZ-PP-F12 results close to conventional estimates of the complete basis set limit. The patterns of convergence are also greatly improved compared to those observed from the use of conventional correlation consistent basis sets in F12 calculations.
Collapse
Affiliation(s)
- J Grant Hill
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, United Kingdom
| | - Robert A Shaw
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, United Kingdom
| |
Collapse
|
19
|
Abstract
The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322-0300, USA
| |
Collapse
|
20
|
McKee NA, McKee ML. Evaluation of packing single and multiple atoms and molecules in the porous organic cage CC3- R. Phys Chem Chem Phys 2021; 23:19255-19268. [PMID: 34524296 DOI: 10.1039/d1cp01934a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The absorption of multiple atoms and molecules, including Kr, Xe, CH4, CO2, C2H2, H2O, and SF6, within CC3-R, a Porous Organic Cage (POC), was calculated and analyzed. The CC3-R molecule has one central cavity and four window sites. Most adsorbents were modeled with either one unit in the central cavity, four units in the window sites, or with five units in both sites. For Xe, the most favorable site was the central one. The CO2 molecule binds about 3 kcal mol-1 in free energy more strongly than CH4 in the central cavity of CC3-R at 300 K which may be enough to allow useful discrimination. Four C2H2 units and four CO2 units are calculated to bind similarly inside CC3-R (ΔH(298 K) = -8.6 and -7.7 kcal mol-1 per unit, respectively). Since H2O is smaller, more waters can easily fit inside. For twelve water molecules, the binding enthalpy per water is ΔH(298 K) = -16.4 kcal mol-1. For comparison, the binding enthalpy of (H2O)12 at the same level of theory (B3LYP/6-31G(d,p)-D3BJ//M06-2X/6-31G(d)) is predicted to be -12.3 kcal mol-1 per water. Finally, the dimerization of CC3-R and the association of CC3-R with CC3-S was studied as well as 3 to 9 iodine atoms enclosed in CC3-R.
Collapse
Affiliation(s)
- Nida A McKee
- Department of Chemistry and Biochemistry, Auburn, AL 36849, USA.
| | - Michael L McKee
- Department of Chemistry and Biochemistry, Auburn, AL 36849, USA.
| |
Collapse
|
21
|
Scheiner S. Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage. J Phys Chem A 2021; 125:6514-6528. [PMID: 34310147 DOI: 10.1021/acs.jpca.1c05431] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322-0300, United States
| |
Collapse
|
22
|
Brakestad A, Wind P, Jensen SR, Frediani L, Hopmann KH. Multiwavelets applied to metal-ligand interactions: Energies free from basis set errors. J Chem Phys 2021; 154:214302. [PMID: 34240981 DOI: 10.1063/5.0046023] [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
Transition metal-catalyzed reactions invariably include steps where ligands associate or dissociate. In order to obtain reliable energies for such reactions, sufficiently large basis sets need to be employed. In this paper, we have used high-precision multiwavelet calculations to compute the metal-ligand association energies for 27 transition metal complexes with common ligands, such as H2, CO, olefins, and solvent molecules. By comparing our multiwavelet results to a variety of frequently used Gaussian-type basis sets, we show that counterpoise corrections, which are widely employed to correct for basis set superposition errors, often lead to underbinding. Additionally, counterpoise corrections are difficult to employ when the association step also involves a chemical transformation. Multiwavelets, which can be conveniently applied to all types of reactions, provide a promising alternative for computing electronic interaction energies free from any basis set errors.
Collapse
Affiliation(s)
- Anders Brakestad
- Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Peter Wind
- Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Stig Rune Jensen
- Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Luca Frediani
- Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | | |
Collapse
|
23
|
Mandal SK, Munshi P. Predicting Accurate Lead Structures for Screening Molecular Libraries: A Quantum Crystallographic Approach. Molecules 2021; 26:molecules26092605. [PMID: 33946965 PMCID: PMC8124947 DOI: 10.3390/molecules26092605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 04/22/2021] [Accepted: 04/24/2021] [Indexed: 11/17/2022] Open
Abstract
Optimization of lead structures is crucial for drug discovery. However, the accuracy of such a prediction using the traditional molecular docking approach remains a major concern. Our study demonstrates that the employment of quantum crystallographic approach-counterpoise corrected kernel energy method (KEM-CP) can improve the accuracy by and large. We select human aldose reductase at 0.66 Å, cyclin dependent kinase 2 at 2.0 Å and estrogen receptor β at 2.7 Å resolutions with active site environment ranging from highly hydrophilic to moderate to highly hydrophobic and several of their known ligands. Overall, the use of KEM-CP alongside the GoldScore resulted superior prediction than the GoldScore alone. Unlike GoldScore, the KEM-CP approach is neither environment-specific nor structural resolution dependent, which highlights its versatility. Further, the ranking of the ligands based on the KEM-CP results correlated well with that of the experimental IC50 values. This computationally inexpensive yet simple approach is expected to ease the process of virtual screening of potent ligands, and it would advance the drug discovery research.
Collapse
|
24
|
Affiliation(s)
- Wesley Michaels
- Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
| | - Yan Zhao
- State Key Laboratory of Silicate Materials for Architectures, Wuhan University of Technology, Wuhan 430070, China
| | - Jian Qin
- Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
25
|
Scheiner S. Relative Strengths of a Pnicogen and a Tetrel Bond and Their Mutual Effects upon One Another. J Phys Chem A 2021; 125:2631-2641. [PMID: 33734698 DOI: 10.1021/acs.jpca.1c01211] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The ability of the T and Z atoms of TR3ZR2 to engage in a noncovalent interaction with NH3 is assessed by DFT calculations, where the T atom refers to C, Si, and Ge; Z = As, Sb, and P; and substituents R = H and F. In most instances, the tetrel bond (TB) is both stronger and shorter than the pnicogen bond (ZB). These two bond strengths can be equalized, or preference shifted to the ZB, if F substituents are placed on the Z and H on the T atoms. Employing C as the T atom results in a very weak TB, with the ZB clearly favored energetically. The simultaneous formation of both TB and ZB weakens both, particularly the latter, but both bonds survive intact. Geometric and spectroscopic perturbations of the subunits reflect the two types of noncovalent bonds.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University†Logan, Utah 84322-0300, United States
| |
Collapse
|
26
|
Karton A, Martin JML. Prototypical π-π dimers re-examined by means of high-level CCSDT(Q) composite ab initio methods. J Chem Phys 2021; 154:124117. [PMID: 33810692 DOI: 10.1063/5.0043046] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The benzene-ethene and parallel-displaced (PD) benzene-benzene dimers are the most fundamental systems involving π-π stacking interactions. Several high-level ab initio investigations calculated the binding energies of these dimers using the coupled-cluster with singles, doubles, and quasi-perturbative triple excitations [CCSD(T)] method at the complete basis set [CBS] limit using various approaches such as reduced virtual orbital spaces and/or MP2-based basis set corrections. Here, we obtain CCSDT(Q) binding energies using a Weizmann-3-type approach. In particular, we extrapolate the self-consistent field (SCF), CCSD, and (T) components using large heavy-atom augmented Gaussian basis sets [namely, SCF/jul-cc-pV{5,6}Z, CCSD/jul-cc-pV{Q,5}Z, and (T)/jul-cc-pV{T,Q}Z]. We consider post-CCSD(T) contributions up to CCSDT(Q), inner-shell, scalar-relativistic, and Born-Oppenheimer corrections. Overall, our best relativistic, all-electron CCSDT(Q) binding energies are ∆Ee,all,rel = 1.234 (benzene-ethene) and 2.550 (benzene-benzene PD), ∆H0 = 0.949 (benzene-ethene) and 2.310 (benzene-benzene PD), and ∆H298 = 0.130 (benzene-ethene) and 1.461 (benzene-benzene PD) kcal mol-1. Important conclusions are reached regarding the basis set convergence of the SCF, CCSD, (T), and post-CCSD(T) components. Explicitly correlated calculations are used as a sanity check on the conventional binding energies. Overall, post-CCSD(T) contributions are destabilizing by 0.028 (benzene-ethene) and 0.058 (benzene-benzene) kcal mol-1, and thus, they cannot be neglected if sub-chemical accuracy is sought (i.e., errors below 0.1 kcal mol-1). CCSD(T)/aug-cc-pwCVTZ core-valence corrections increase the binding energies by 0.018 (benzene-ethene) and 0.027 (benzene-benzene PD) kcal mol-1. Scalar-relativistic and diagonal Born-Oppenheimer corrections are negligibly small. We use our best CCSDT(Q) binding energies to evaluate the performance of MP2-based, CCSD-based, and lower-cost composite ab initio procedures for obtaining these challenging π-π stacking binding energies.
Collapse
Affiliation(s)
- Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
27
|
Scheiner S. Comparison of Bifurcated Halogen with Hydrogen Bonds. Molecules 2021; 26:molecules26020350. [PMID: 33445461 PMCID: PMC7827642 DOI: 10.3390/molecules26020350] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Revised: 01/08/2021] [Accepted: 01/09/2021] [Indexed: 01/30/2023] Open
Abstract
Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH3 and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, UT 84322-0300, USA
| |
Collapse
|
28
|
Scheiner S. Competition between a Tetrel and Halogen Bond to a Common Lewis Acid. J Phys Chem A 2020; 125:308-316. [DOI: 10.1021/acs.jpca.0c10060] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322-0300, United States
| |
Collapse
|
29
|
Semidalas E, Martin JML. Canonical and DLPNO-Based Composite Wavefunction Methods Parametrized against Large and Chemically Diverse Training Sets. 2: Correlation-Consistent Basis Sets, Core-Valence Correlation, and F12 Alternatives. J Chem Theory Comput 2020; 16:7507-7524. [PMID: 33200931 PMCID: PMC7735707 DOI: 10.1021/acs.jctc.0c01106] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A hierarchy
of wavefunction composite methods (cWFT), based on
G4-type cWFT methods available for elements H through Rn, was recently
reported by the present authors [2020, 16, 4238]. We extend this hierarchy
by considering the inner-shell correlation energy in the second-order
Møller–Plesset correction and replacing the Weigend–Ahlrichs
def2-mZVPP(D) basis sets used with complete basis
set extrapolation from augmented correlation-consistent core–valence
triple-ζ, aug-cc-pwCVTZ(-PP), and quadruple-ζ, aug-cc-pwCVQZ(-PP),
basis sets, thus creating cc-G4-type methods. For the large and chemically
diverse GMTKN55 benchmark suite, they represent a substantial further
improvement and bring WTMAD2 (weighted mean absolute deviation) down
below 1 kcal/mol. Intriguingly, the lion’s share of the improvement
comes from better capture of valence correlation; the inclusion of
core–valence correlation is almost an order of magnitude less
important. These robust correlation-consistent cWFT methods approach
the CCSD(T) complete basis limit with just one or a few fitted parameters.
Particularly, the DLPNO variants such as cc-G4-T-DLPNO are applicable
to fairly large molecules at a modest computational cost, as is (for
a reduced range of elements) a different variant using MP2-F12/cc-pVTZ-F12
for the MP2 component.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Rehovot, Israel
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Rehovot, Israel
| |
Collapse
|
30
|
Santos GFN, Carvalho LC, Oliveira DAS, Rego DG, Bueno MA, Oliveira BG. The definitive challenge of forming uncommon pseudo‐π···H–F and C···H–F hydrogen bonds on cyclic and cubic nonpolar hydrocarbons. J PHYS ORG CHEM 2020. [DOI: 10.1002/poc.4098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
| | - Leila Cardoso Carvalho
- Centro das Ciências Exatas e das Tecnologias Universidade Federal do Oeste da Bahia Barreiras Brazil
| | | | - Danilo Guimarães Rego
- Centro das Ciências Exatas e das Tecnologias Universidade Federal do Oeste da Bahia Barreiras Brazil
| | - Mauro Alves Bueno
- Centro das Ciências Exatas e das Tecnologias Universidade Federal do Oeste da Bahia Barreiras Brazil
| | - Boaz Galdino Oliveira
- Centro das Ciências Exatas e das Tecnologias Universidade Federal do Oeste da Bahia Barreiras Brazil
| |
Collapse
|
31
|
Lu J, Scheiner S. Relationships between Bond Strength and Spectroscopic Quantities in H-Bonds and Related Halogen, Chalcogen, and Pnicogen Bonds. J Phys Chem A 2020; 124:7716-7725. [DOI: 10.1021/acs.jpca.0c05936] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Jia Lu
- Department of Chemistry and Biochemistry, Utah State University , Logan, Utah 84322-0300, United States
| | - Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University , Logan, Utah 84322-0300, United States
| |
Collapse
|
32
|
Peng C, Lewis CA, Wang X, Clement MC, Pierce K, Rishi V, Pavošević F, Slattery S, Zhang J, Teke N, Kumar A, Masteran C, Asadchev A, Calvin JA, Valeev EF. Massively Parallel Quantum Chemistry: A high-performance research platform for electronic structure. J Chem Phys 2020; 153:044120. [DOI: 10.1063/5.0005889] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Cannada A. Lewis
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xiao Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | | | - Karl Pierce
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Fabijan Pavošević
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Samuel Slattery
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Jinmei Zhang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Nakul Teke
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Conner Masteran
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Andrey Asadchev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
33
|
de Oliveira BG, Zabardasti A, do Rego DG, Pour MM. The formation of H···X hydrogen bond, C···X carbon-halide or Si···X tetrel bonds on the silylene-halogen dimers (X = F or Cl): intermolecular strength, molecular orbital interactions and prediction of covalency. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-02644-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
34
|
Otero-de-la-Roza A, DiLabio GA. Improved Basis-Set Incompleteness Potentials for Accurate Density-Functional Theory Calculations in Large Systems. J Chem Theory Comput 2020; 16:4176-4191. [PMID: 32470304 DOI: 10.1021/acs.jctc.0c00102] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate calculation of chemical properties using density-functional theory (DFT) requires the use of a nearly complete basis set. In chemical systems involving hundreds to thousands of atoms, the cost of the calculations place practical limitations on the number of basis functions that can be used. Therefore, in most practical applications of DFT to large systems, there exists a basis-set incompleteness error (BSIE). In this article, we present the next iteration of the basis-set incompleteness potentials (BSIPs), one-electron potentials designed to correct for basis-set incompleteness error. The ultimate goal associated with the development of BSIPs is to allow the calculation of molecular properties using DFT with near-complete-basis-set results at a computational cost that is similar to a small basis set calculation. In this work, we develop BSIPs for 10 atoms in the first and second rows (H, B-F, Si-Cl) and 15 common basis sets of the Pople, Dunning, Karlsruhe, and Huzinaga types. Our new BSIPs are constructed to minimize BSIE in the calculation of reaction energies, barrier heights, noncovalent binding energies, and intermolecular distances. The BSIPs were obtained using a training set of 15 944 data points. The fitting approach employed a regularized linear least-squares method with variable selection (the LASSO method), which results in a much better fit to the training data than our previous BSIPs while, at the same time, reducing the computational cost of BSIP development. The proposed BSIPs are tested on various benchmark sets and demonstrate excellent performance in practice. Our new BSIPs are also transferable; i.e., they can be used to correct BSIE in calculations that employ density functionals other than the one used in the BSIP development (B3LYP). Finally, BSIPs can be used in any quantum chemistry program that have implemented effective-core potentials without changes to the software.
Collapse
Affiliation(s)
- A Otero-de-la-Roza
- Departamento de Quı́mica Fı́sica y Analítica and MALTA-Consolider Team, Facultad de Quı́mica, Universidad de Oviedo, 33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada.,Faculty of Management, University of British Columbia, Okanagan, 1137 Alumni Avenue, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
35
|
Semidalas E, Martin JML. Canonical and DLPNO-Based G4(MP2)XK-Inspired Composite Wave Function Methods Parametrized against Large and Chemically Diverse Training Sets: Are They More Accurate and/or Robust than Double-Hybrid DFT? J Chem Theory Comput 2020; 16:4238-4255. [PMID: 32456427 PMCID: PMC7366511 DOI: 10.1021/acs.jctc.0c00189] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The
large and chemically diverse GMTKN55 benchmark was used as
a training set for parametrizing composite wave function thermochemistry
protocols akin to G4(MP2)XK theory (Chan, B.; Karton, A.; Raghavachari,
K. J. Chem. Theory Comput. 2019, 15, 4478–4484). On account of their availability
for elements H through Rn, Karlsruhe def2 basis sets were employed.
Even after reparametrization, the GMTKN55 WTMAD2 (weighted mean absolute
deviation, type 2) for G4(MP2)-XK is actually inferior to that of
the best rung-4 DFT functional, ωB97M-V. By increasing the basis
set for the MP2 part to def2-QZVPPD, we were able to substantially
improve performance at modest cost (if an RI-MP2 approximation is
made), with WTMAD2 for this G4(MP2)-XK-D method now comparable to
the better rung-5 functionals (albeit at greater cost). A three-tier
approach with a scaled MP3/def2-TZVPP intermediate step, however,
leads to a G4(MP3)-D method that is markedly superior to even the
best double hybrids ωB97M(2) and revDSD-PBEP86-D4. Evaluating
the CCSD(T) component with a triple-ζ, rather than split-valence,
basis set yields only a modest further improvement that is incommensurate
with the drastic increase in computational cost. G4(MP3)-D and G4(MP2)-XK-D
have about 40% better WTMAD2, at similar or lower computational cost,
than their counterparts G4 and G4(MP2), respectively: detailed comparison
reveals that the difference lies in larger molecules due to basis
set incompleteness error. An E2/{T,Q} extrapolation and a CCSD(T)/def2-TZVP
step provided the G4-T method of high accuracy and with just three
fitted parameters. Using KS orbitals in MP2 leads to the G4(MP3|KS)-D
method, which entirely eliminates the CCSD(T) step and has no steps
costlier than scaled MP3; this shows a path forward to further improvements
in double-hybrid density functional methods. None of our final selections
require an empirical HLC correction; this cuts the number of empirical
parameters in half and avoids discontinuities on potential energy
surfaces. G4-T-DLPNO, a variant in which post-MP2 corrections are
evaluated at the DLPNO-CCSD(T) level, achieves nearly the accuracy
of G4-T but is applicable to much larger systems.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḩovot, Israel
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḩovot, Israel
| |
Collapse
|
36
|
Socha O, Dračínský M. Dimerization of Acetic Acid in the Gas Phase-NMR Experiments and Quantum-Chemical Calculations. Molecules 2020; 25:molecules25092150. [PMID: 32375390 PMCID: PMC7248931 DOI: 10.3390/molecules25092150] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2020] [Revised: 04/28/2020] [Accepted: 04/30/2020] [Indexed: 11/24/2022] Open
Abstract
Due to the nature of the carboxylic group, acetic acid can serve as both a donor and acceptor of a hydrogen bond. Gaseous acetic acid is known to form cyclic dimers with two strong hydrogen bonds. However, trimeric and various oligomeric structures have also been hypothesized to exist in both the gas and liquid phases of acetic acid. In this work, a combination of gas-phase NMR experiments and advanced computational approaches were employed in order to validate the basic dimerization model of gaseous acetic acid. The gas-phase experiments performed in a glass tube revealed interactions of acetic acid with the glass surface. On the other hand, variable-temperature and variable-pressure NMR parameters obtained for acetic acid in a polymer insert provided thermodynamic parameters that were in excellent agreement with the MP2 (the second order Møller–Plesset perturbation theory) and CCSD(T) (coupled cluster with single, double and perturbative triple excitation) calculations based on the basic dimerization model. A slight disparity between the theoretical dimerization model and the experimental data was revealed only at low temperatures. This observation might indicate the presence of other, entropically disfavored, supramolecular structures at low temperatures.
Collapse
Affiliation(s)
- Ondřej Socha
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague, Czech Republic;
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 166 10 Prague, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague, Czech Republic;
- Correspondence: ; Tel./Fax: +42-02-2018-3139
| |
Collapse
|
37
|
Zhao S, Wu Q, Pi J, Liu J, Zheng J, Hou S, Wei J, Li R, Sadeghi H, Yang Y, Shi J, Chen Z, Xiao Z, Lambert C, Hong W. Cross-plane transport in a single-molecule two-dimensional van der Waals heterojunction. SCIENCE ADVANCES 2020; 6:eaba6714. [PMID: 32524003 PMCID: PMC7259930 DOI: 10.1126/sciadv.aba6714] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2019] [Accepted: 03/31/2020] [Indexed: 06/01/2023]
Abstract
Two-dimensional van der Waals heterojunctions (2D-vdWHs) stacked from atomically thick 2D materials are predicted to be a diverse class of electronic materials with unique electronic properties. These properties can be further tuned by sandwiching monolayers of planar organic molecules between 2D materials to form molecular 2D-vdWHs (M-2D-vdWHs), in which electricity flows in a cross-plane way from one 2D layer to the other via a single molecular layer. Using a newly developed cross-plane break junction technique, combined with density functional theory calculations, we show that M-2D-vdWHs can be created and that cross-plane charge transport can be tuned by incorporating guest molecules. The M-2D-vdWHs exhibit distinct cross-plane charge transport signatures, which differ from those of molecules undergoing in-plane charge transport.
Collapse
Affiliation(s)
- Shiqiang Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Qingqing Wu
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK
| | - Jiuchan Pi
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Junyang Liu
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Jueting Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Songjun Hou
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK
| | - Junying Wei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Ruihao Li
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Hatef Sadeghi
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK
| | - Yang Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, China
| | - Jia Shi
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Zhaobin Chen
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Zongyuan Xiao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, China
| | - Colin Lambert
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK
| | - Wenjing Hong
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering & Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, China
| |
Collapse
|
38
|
Yu F, Wang Y. Dual‐hybrid direct random phase approximation and second‐order screened exchange with nonlocal van der Waals correlations for noncovalent interactions. J Comput Chem 2020; 41:1018-1025. [DOI: 10.1002/jcc.26149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 01/05/2020] [Accepted: 01/06/2020] [Indexed: 11/09/2022]
Affiliation(s)
- Feng Yu
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| | - Yaoting Wang
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| |
Collapse
|
39
|
Schmitz G, Elm J. Assessment of the DLPNO Binding Energies of Strongly Noncovalent Bonded Atmospheric Molecular Clusters. ACS OMEGA 2020; 5:7601-7612. [PMID: 32280904 PMCID: PMC7144154 DOI: 10.1021/acsomega.0c00436] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Accepted: 03/11/2020] [Indexed: 05/03/2023]
Abstract
This work assesses the performance of DLPNO-CCSD(T0), DLPNO-MP2, and density functional theory methods in calculating the binding energies of a representative test set of 45 atmospheric acid-acid, acid-base, and acid-water dimer clusters. The performance of the approximate methods is compared to high level explicitly correlated CCSD(F12*)(T)/complete basis set (CBS) reference calculations. Out of the tested density functionals, ωB97X-D3(BJ) shows the best performance with a mean deviation of 0.09 kcal/mol and a maximum deviation of 0.83 kcal/mol. The RI-CC2/aug-cc-pV(T+d)Z level of theory severely overpredicts the cluster binding energies with a mean deviation of -1.31 kcal/mol and a maximum deviation up to -3.00 kcal/mol. Hence, RI-CC2/aug-cc-pV(T+d)Z should not be utilized for studying atmospheric molecular clusters. The DLPNO variants are tested both with and without the inclusion of explicit correlation (F12) in the wavefunction, with different pair natural orbital (PNO) settings (loosePNO, normalPNO, and tightPNO) and using both double and triple zeta basis sets. The performance of the DLPNO-MP2 methods is found to be independent of PNO settings and yield low mean deviations of -0.84 kcal/mol or below. However, DLPNO-MP2 requires explicitly correlated wavefunctions to yield maximum deviations below 1.40 kcal/mol. For obtaining high accuracy, with maximum deviation below ∼1.0 kcal/mol, either DLPNO-CCSD(T0)/aug-cc-pVTZ (normalPNO) calculations or DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 (normalPNO) calculations are required. The most accurate level of theory is found to be DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 using a tightPNO criterion which yields a mean deviation of 0.10 kcal/mol, with a maximum deviation of 0.20 kcal/mol, compared to the CCSD(F12*)(T)/CBS reference.
Collapse
Affiliation(s)
- Gunnar Schmitz
- Department
of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Jonas Elm
- Department
of Chemistry and iClimate, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| |
Collapse
|
40
|
Szczepaniak M, Moc J. Elusive Cyanoform: Computational Probing Its Stability and Reactivity with Accurate Ab Initio Methods. J Phys Chem A 2020; 124:2634-2648. [DOI: 10.1021/acs.jpca.0c00540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Marek Szczepaniak
- Faculty of Chemistry, Wroclaw University, F. Joliot-Curie 14, 50-383 Wroclaw, Poland
| | - Jerzy Moc
- Faculty of Chemistry, Wroclaw University, F. Joliot-Curie 14, 50-383 Wroclaw, Poland
| |
Collapse
|
41
|
Nieman R, Aquino AJA, Lischka H. Benchmark ab initio calculations on intermolecular structures and the exciton character of poly(p-phenylenevinylene) dimers. J Chem Phys 2020; 152:044306. [DOI: 10.1063/1.5139411] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Reed Nieman
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Adelia J. A. Aquino
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, People’s Republic of China
| | - Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, People’s Republic of China
| |
Collapse
|
42
|
Martin JML, Santra G. Empirical Double‐Hybrid Density Functional Theory: A ‘Third Way’ in Between WFT and DFT. Isr J Chem 2019. [DOI: 10.1002/ijch.201900114] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Jan M. L. Martin
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| | - Golokesh Santra
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| |
Collapse
|
43
|
Abstract
Many-body dispersion has gained considerable attention over the past decade, particularly for condensed phase systems. However, quantitatively accurate studies of many-body dispersion have only recently become feasible due to challenges in reliability and accuracy. Currently available methodologies for calculating many-body dispersion have been challenged, with recent evidence suggesting, for example, that dispersion-corrected density functional theory (DFT) schemes cannot consistently predict many-body dispersion accurately. This study evaluates many-body dispersion energies using a composite approach that employs singles and doubles coupled cluster theory with perturbative/noniterative triples, CCSD(T), combined with an extrapolation to the complete basis set (CBS) limit. The combined CCSD(T)/CBS approach is applied to Arn and (H2O)n, n = 3-10, clusters, and a new data set called S22(3), which includes trimers generated based on the S22 data set. In these systems, the many-body dispersion provides a very small contribution to the total interaction energy of all of the systems studied, generally 3% or less of the total interaction energy. Two-body dispersion is the dominant dispersion contribution and many-body dispersion contributes no more than 5.7% of the total dispersion energy, generally staying below 2%.
Collapse
Affiliation(s)
- Melisa Alkan
- Department of Chemistry , Iowa State University , Ames , Iowa 50011 , United States.,Ames Laboratory , Ames , Iowa 50011 , United States
| | - Peng Xu
- Department of Chemistry , Iowa State University , Ames , Iowa 50011 , United States.,Ames Laboratory , Ames , Iowa 50011 , United States
| | - Mark S Gordon
- Department of Chemistry , Iowa State University , Ames , Iowa 50011 , United States.,Ames Laboratory , Ames , Iowa 50011 , United States
| |
Collapse
|
44
|
Wang C, Mo Y. Classical Electrostatic Interaction Is the Origin for Blue-Shifting Halogen Bonds. Inorg Chem 2019; 58:8577-8586. [DOI: 10.1021/acs.inorgchem.9b00875] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Changwei Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi’an 710119, China
| | - Yirong Mo
- Department of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, United States
| |
Collapse
|
45
|
Seifert NA, Hazrah AS, Jäger W. The 1-Naphthol Dimer and Its Surprising Preference for π-π Stacking over Hydrogen Bonding. J Phys Chem Lett 2019; 10:2836-2841. [PMID: 31002249 DOI: 10.1021/acs.jpclett.9b00646] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Using chirped-pulse Fourier transform microwave spectroscopy, we have experimentally identified the most abundant dimer of an alcohol analogue of naphthalene, 1-naphthol. The 1-naphthol dimer features a V-shaped, partially overlapping π-π stacked structure with no canonical hydrogen bonds between the subunits. This structural assignment is in contradiction to an earlier study of the 1-naphthol dimer using UV-IR dip double resonance spectroscopy in the O-H stretch region, which assigns a π-stacked but also canonically hydrogen-bonded structure. We use an improved theoretical analysis to resolve this discrepancy and show that the new, V-shaped structure is also consistent with the previously measured UV-IR data. These new computational and spectroscopic results shed light on the complicated nature of evaluating energetics and structures for larger, dispersion-bound systems.
Collapse
Affiliation(s)
- Nathan A Seifert
- Department of Chemistry University of Alberta 11227 Saskatchewan Drive , Edmonton , Alberta Canada T6G 2G2
| | - Arsh S Hazrah
- Department of Chemistry University of Alberta 11227 Saskatchewan Drive , Edmonton , Alberta Canada T6G 2G2
| | - Wolfgang Jäger
- Department of Chemistry University of Alberta 11227 Saskatchewan Drive , Edmonton , Alberta Canada T6G 2G2
| |
Collapse
|
46
|
Complexation of free-base and 3d transition metal(II) phthalocyanines with endohedral fullerene Sc3N@C80. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.03.019] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
47
|
Carvalho LC, Bueno MA, de Oliveira BG. The interplay and strength of the π⋯HF, C⋯HF, F⋯HF and F⋯HC hydrogen bonds upon the formation of multimolecular complexes based on C 2H 2⋯HF and C 2H 4⋯HF small dimers. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2019; 213:438-455. [PMID: 30738351 DOI: 10.1016/j.saa.2019.01.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2018] [Revised: 12/20/2018] [Accepted: 01/01/2019] [Indexed: 06/09/2023]
Abstract
The conception of this theoretical research was idealized aiming to unveil the intermolecular structures of complexes formed by acetylene or ethylene and hydrofluoric acid. At light of computational calculations by using the B3LYP/6-311++G(d,p) method, the geometries of the C2H2⋯(HF), C2H2⋯2(HF), C2H2⋯4(HF), C2H4⋯(HF), C2H4⋯2(HF) and C2H4⋯4(HF) hydrogen-bonded complexes were fully optimized. Moreover, the Post-Hartree-Fock calculations MP2/6-311++G(d,p), MP2/aug-cc-pVTZ, MP4(SDQ)/6-311++G(d,p) and CCSD/6-311++G(d,p) also were also used. The infrared spectra were analyzed in order to identify the new vibrational modes and frequencies of the proton donors shifted to red region. Through the modeling of charge-fluxes on the basis of the Quantum Theory of Atoms In Molecules (QTAIM) and, by contradicting the expectation of the hydrofluorination mechanisms of acetylene or ethylene, C⋯HF was recognized as a new type of hydrogen bond instead of the already well known π⋯H. The calculations of the Natural Bonding Orbital (NBO) and Charges derived from the Electrostatic Potential Grid-based (ChElPG) were also applied to interpret the shifting frequencies as well as measuring of the punctual charge-transfer after the formation of the complexes. Finally, the determination of the stabilization energy was carried out through the arguments of the Fock matrix in NBO basis and through the supermolecule approach. Also it is worthwhile to notice that some algebraic formulations were used for determining the electronic cooperative effect (CE).
Collapse
|
48
|
Ma Q, Werner HJ. Accurate Intermolecular Interaction Energies Using Explicitly Correlated Local Coupled Cluster Methods [PNO-LCCSD(T)-F12]. J Chem Theory Comput 2019; 15:1044-1052. [DOI: 10.1021/acs.jctc.8b01098] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
49
|
Kruse H, Banáš P, Šponer J. Investigations of Stacked DNA Base-Pair Steps: Highly Accurate Stacking Interaction Energies, Energy Decomposition, and Many-Body Stacking Effects. J Chem Theory Comput 2018; 15:95-115. [DOI: 10.1021/acs.jctc.8b00643] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 Listopadu 12, 77146 Olomouc, Czech Republic
| | - Jiřı́ Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 Listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
50
|
Maleki F, Salehzadeh S. Extension of the atom by atom scheme of counterpoise method and presentation of its new advantages. J Chem Phys 2018; 149:064116. [DOI: 10.1063/1.5037576] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Farahnaz Maleki
- Department of Chemistry, Bu-Ali Sina University, Hamedan, Iran
| | | |
Collapse
|