A DLPNO-CCSD(T) benchmarking study of intermolecular interactions of ionic liquids

Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes

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.

Unravelling ionic liquid solvent effects for a non-polar Cope rearrangement reaction

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.

Unraveling the atmospheric oxidation mechanism and kinetics of naphthalene: Insights from theoretical exploration

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.

MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations

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.

On the existence of CO32- microsolvated clusters: a theoretical study

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 [CO₃·(H₂O)_n]^2-. The reasons for this may be related to the thermodynamical or kinetical instability of microsolvated CO₃^2- toward autoionization or autoprotonation processes. In this study we theoretically investigate the potential stability of the [CO₃·(H₂O)_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 CO₃^2- anions in bulk water solvent was analysed to highlight important precautions for synthetic purposes.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

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.