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
|
Smith GJ, Koutsoukos S, Lancaster B, Becker J, Welton T, Hunt PA. Unravelling ionic liquid solvent effects for a non-polar Cope rearrangement reaction. Phys Chem Chem Phys 2024; 26:12453-12466. [PMID: 38625536 DOI: 10.1039/d4cp00156g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
The impact of ionic liquids (ILs) on polar reactions is well recognised, however the impact of ILs on non-polar reactions is less well understood or explored. Pericyclic Cope rearrangements are highly concerted, exhibit minimal charge localisation and pass through an uncharged but well-defined transition state, and thus provide a good mechanism for exploring the impact of IL polarizability on chemical reactivity. Recently, a 10× rate enhancement has been observed for the Cope rearrangement of 3-phenyl-1,5-hexadiene in the IL 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C4C1im][NTf2] compared to benzene. In this work we undertake a DFT based computational study (B3LYP-D3BJ/6-311+G(d,p) and M06-2X-D3/6-311+G(d,p)) of the Cope rearrangement of 3-phenyl-1,5-hexadiene and 3-propyl-hexa-1,5-diene in molecular solvents (acetonitrile, benzene and ethanol) and the IL [C4C1im][NTf2] using the SMD solvation model. The impact of benzene and [C4C1im][NTf2] on the Cope rearrangement of 3-phenyl-1,5-hexadiene is studied in more detail and we provide insight into the reason for the rate enhancement in an IL. The volume of activation is evaluated and the potential impact of 'solvent pressure' is discussed. We identify two potential mechanisms for volume effects to contribute to the rate enhancement. Solvent association energies are evaluated at the DLNPO-CCSD(T) level. Specific solvent interactions are explored through atomic partial charge, molecular orbital and bond critical point analysis, as well as via non-colvalent interaction (NCI) plots, electrostatic potential (ESP) differences and density difference Δρ(r) plots. We find that the cation and anion together form an extensive van der Waals pocket in-which the transition state (TS) sits. Electron density within the TS is anisotropically polarised via a 'push-pull' effect due to the dual cation-anion nature of the IL, stabilising the TS relative to benzene. We also provide experimental evidence that these effects are generalisable to other ILs. Overall, our aim is to provide a deeper moleuclar level understanding of the impact of ILs on non-polar reactions.
Collapse
Affiliation(s)
- Gavin J Smith
- Department of Chemistry, Imperial College London, UK
| | | | - Ben Lancaster
- Department of Chemistry, Imperial College London, UK
| | - Julian Becker
- Department of Chemistry, Imperial College London, UK
| | - Tom Welton
- Department of Chemistry, Imperial College London, UK
| | - Patricia A Hunt
- Department of Chemistry, Imperial College London, UK
- School of Chemical and Physical Sciences, Victoria University of Wellington, New Zealand.
| |
Collapse
|
3
|
Wu X, Yao X, Xie B, Wang P, Huo W, Zhu Y, Hou Q, Wu M, Wu Y, Zhang F. Unraveling the atmospheric oxidation mechanism and kinetics of naphthalene: Insights from theoretical exploration. CHEMOSPHERE 2024; 352:141356. [PMID: 38309603 DOI: 10.1016/j.chemosphere.2024.141356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/24/2023] [Revised: 01/28/2024] [Accepted: 01/31/2024] [Indexed: 02/05/2024]
Abstract
Naphthalene, the most abundant polycyclic aromatic hydrocarbon in the atmosphere, significantly influences OH consumption and secondary organic aerosol (SOA) formation. Naphthoquinone (NQ) is a significant contributor to ring-retaining SOA from naphthalene degradation, impacting the redox properties and toxicity of ambient particles. However, inconsistencies persist regarding concentrations of its isomers, 1,2-NQ and 1,4-NQ. In present work, our theoretical investigation into naphthalene's reaction with OH and subsequent oxygenation unveils their role in SOA formation. The reaction kinetics of initial OH and subsequent O2 oxidation was extensively studied using high-level quantum chemical methods (DLPNO-CCSD(T)/aug-ccpVQZ//M052x-D3/6-311++G(d,p)) combined with RRKM/master equation simulations. The reactions mainly proceed through electrophilic addition and abstraction from the aromatic ring. The total rate coefficient of naphthalene + OH at 300 K and 1 atm from our calculation (7.2 × 10-12 cm3 molecule-1 s-1) agrees well with previous measurements (∼1 × 10-11 cm3 molecule-1 s-1). The computed branching ratios facilitate accurate product yield determination. The largest yield of 1-hydroxynaphthalen-1-yl radical (add1) producing the major precursor of RO2 is computed to be 93.8 % in the ambient environment. Our calculated total rate coefficient (5.2 × 10-16 cm3 molecule-1 s-1) for add1 + O2 closely matches that of limited experimental data (8.0 × 10-16 cm3 molecule-1 s-1). Peroxy radicals (RO2) generated from add1 + O2 include 4-cis/trans-(1-hydroxynaphthalen-1-yl)-peroxy radical (add1-4OOadd-cis/trans, 66.0 %/17.5 %), 2-cis/trans-(1-hydroxynaphthalen-1-yl)-peroxy radical (add1-2OOadd-cis/trans, 10.3 %/6.3 %). Regarding the debated predominance of 1,4-NQ (corresponding to the parent RO2, i.e., add1-4OOadd-cis/trans) and 1,2-NQ (corresponding to the parent RO2, i.e., add1-2OOadd-cis/trans) in the atmosphere, our findings substantiate the dominance of 1,4-NQ. This study also indicates potential weakening of 1,4-NQ's dominance due to competition from decomposition reactions of add1-4OOadd-cis/trans and add1-2OOadd-cis/trans. Precise reaction kinetics data are essential for characterizing SOA transformation derived from naphthalene and assessing their climatic impacts within modeling frameworks.
Collapse
Affiliation(s)
- Xiaoqing Wu
- College of Information Engineering, China Jiliang University, Hangzhou, 310018, PR China; Science and Technology on Plasma Dymamics Lab, Aeronautics and Astronautics Engineering College, Air Force Engineering University, Xian, 710038, PR China.
| | - Xiaoxia Yao
- Science and Technology on Plasma Dymamics Lab, Aeronautics and Astronautics Engineering College, Air Force Engineering University, Xian, 710038, PR China.
| | - Binbin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou, 311231, PR China.
| | - Pengfei Wang
- National Institute of Extremely-Weak Magnetic Field Infrastructure, Hangzhou, 310051, PR China.
| | - Wanli Huo
- College of Information Engineering, China Jiliang University, Hangzhou, 310018, PR China.
| | - Yifei Zhu
- Institute of Aero-engine, School of Mechanical Engineering, Xi'an Jiaotong University, Xi'an, PR China.
| | - Qifeng Hou
- National Institute of Extremely-Weak Magnetic Field Infrastructure, Hangzhou, 310051, PR China.
| | - Mengqi Wu
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, PR China; Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui, 230088, PR China.
| | - Yun Wu
- Science and Technology on Plasma Dymamics Lab, Aeronautics and Astronautics Engineering College, Air Force Engineering University, Xian, 710038, PR China.
| | - Feng Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, PR China; Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui, 230088, PR China.
| |
Collapse
|
4
|
Pogrebetsky J, Siklitskaya A, Kubas A. MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations. J Chem Theory Comput 2023. [PMID: 37338422 DOI: 10.1021/acs.jctc.3c00444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method has been proven to provide accurate single-point energies at a fraction of the cost of canonical CCSD(T) calculations. However, the desired "chemical accuracy" can only be obtained with a large PNO space and extended basis set. We present a simple yet accurate and efficient correction scheme based on a perturbative approach. Here, in addition to DLPNO-CCSD(T) energy, one calculates DLPNO-MP2 correlation energy with the same settings as in the preceding coupled-cluster calculation. In the next step, the canonical MP2 correlation energy is obtained in the same orbital basis. This can be efficiently performed for essentially all molecule sizes accessible with the DLPNO-CCSD(T) method. By taking the difference between the canonical MP2 and DLPNO-MP2 energies, we obtain a correction term that can be added to the DLPNO-CCSD(T) correlation energy. This way, one can obtain the total correlation energy close to the limit of the complete PNO space (cPNO). The presented approach allows us to significantly increase the accuracy of the DLPNO-CCSD(T) method for both closed- and open-shell systems. The latter are known to be especially challenging for locally correlated methods. Unlike the previously developed PNO extrapolation procedure by Altun, Neese, and Bistoni ( J. Chem. Theory Comput. 2020, 16, 6142-6149), this strategy allows us to get the DLPNO-CCSD(T) correlation energy at the cPNO limit in a cost-efficient way, resulting in a minimal overall increase in calculation time as compared to the uncorrected method.
Collapse
Affiliation(s)
- James Pogrebetsky
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| | - Alexandra Siklitskaya
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| |
Collapse
|
5
|
Rublev P, Tkachenko NV, Dub PA, Boldyrev AI. On the existence of CO 32- microsolvated clusters: a theoretical study. Phys Chem Chem Phys 2023; 25:14046-14055. [PMID: 37161655 DOI: 10.1039/d3cp00955f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Microsolvated clusters of multiply charged anions play a crucial role in atmospheric chemistry and some of them were previously registered experimentally. At the same time, there are no experimental observations of [CO3·(H2O)n]2-. The reasons for this may be related to the thermodynamical or kinetical instability of microsolvated CO32- toward autoionization or autoprotonation processes. In this study we theoretically investigate the potential stability of the [CO3·(H2O)n]2- microsolvated clusters from both perspectives - thermodynamic and kinetic - and we claim they are stable toward autoionization and kinetically semi-stable toward autoprotonation. In addition, the behaviour of CO32- anions in bulk water solvent was analysed to highlight important precautions for synthetic purposes.
Collapse
Affiliation(s)
- Pavel Rublev
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, USA.
| | - Nikolay V Tkachenko
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, USA.
| | - Pavel A Dub
- Schrödinger Inc., San Diego, California 92121, USA
| | - Alexander I Boldyrev
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, USA.
| |
Collapse
|
6
|
Kee CW. Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities. Molecules 2023; 28:molecules28041715. [PMID: 36838703 PMCID: PMC9966076 DOI: 10.3390/molecules28041715] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2022] [Revised: 02/03/2023] [Accepted: 02/05/2023] [Indexed: 02/25/2023] Open
Abstract
Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.
Collapse
Affiliation(s)
- Choon Wee Kee
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR), 1 Pesek Road, Jurong Island, Singapore 627833, Republic of Singapore
| |
Collapse
|