Exploring the Accuracy Limits of Local Pair Natural Orbital Coupled-Cluster Theory

Chemical Affinity of Ag-Exchanged Zeolites for Efficient Hydrogen Isotope Separation

We report an ion-exchanged zeolite as an excellent candidate for large-scale application in hydrogen isotope separation. Ag(I)-exchanged zeolite Y has been synthesized through a standard ion-exchange procedure. The D₂/H₂ separation performance has been systematically investigated via thermal desorption spectroscopy (TDS). Undercoordinated Ag⁺ in zeolite AgY acts as a strong adsorption site and adorbs preferentially the heavier isotopologue even above liquid nitrogen temperature. The highest D₂/H₂ selectivity of 10 is found at an exposure temperature of 90 K. Furthermore, the high Al content of the zeolite structure leads to a high density of Ag sites, resulting in a high gas uptake. In the framework, approximately one-third of the total physisorbed hydrogen isotopes are adsorbed on the Ag sites, corresponding to 3 mmol/g. A density functional theory (DFT) calculation reveals that the isotopologue-selective adsorption of hydrogen at Ag sites contributes to the outstanding hydrogen isotope separation, which has been directly observed through cryogenic thermal desorption spectroscopy. The overall performance of zeolite AgY, showing good selectivity combined with high gas uptake, is very promising for future technical applications.

Silver-exchanged zeolite Y was synthesized for efficient and effective separation of hydrogen isotopes above liquid nitrogen temperature via chemical affinity sieving. The highest D₂/H₂ selectivity of 10 is achieved at 90 K combined with a high gas uptake, making zeolite AgY a promising candidate for large-scale deuterium enrichment.

Strong Bases Design: Predicted Limits of Basicity

Brønsted superbases have wide applications in organic chemistry due to their ability to activate C-H bonds. The strongest neutral bases to date are substituted aminophosphazenes developed in the late 1980s by Reinhard Schwesinger. Since then, much effort has been expended to create even stronger neutral bases. In this article, the reasons for the instability of very basic compounds are investigated by means of high-level quantum-chemical calculations. Theoretical basicity limits are suggested for solutions as well as for the gas phase. A record-breaking superbase most likely to be synthesizable and stable at ambient conditions is proposed. Hexamethylphosphoramide is considered a reliable ionizing solvent for superbases.

Can domain-based local pair natural orbitals approaches accurately predict phosphorescence energies?

Since the discovery of the peculiar conducting and optical properties of aromatics, many efforts have been made to characterize and predict their phosphorescence. This physical process is exploited in modern Organic Emitting Light Diodes (OLEDs), and it is also one of the processes decreasing the efficiency of Dye-sensitized solar cells (DSSCs). Herein, we propose a computational strategy for the accurate calculation of singlet-triplet gaps of aromatic compounds, which provides results that are in excellent agreement with available experimental data. Our approach relies on the domain-based local pair natural orbital (DLPNO) variant of the "gold standard" CCSD(T) method. The convergence of our results with respect to the key technical parameters of the calculation, such as the basis set used, the approximations employed in the perturbative triples correction, and the dimension of the PNOs space, was thoroughly discussed.

Collective Synthesis of Illudalane Sesquiterpenes via Cascade Inverse Electron Demand (4 + 2) Cycloadditions of Thiophene S,S-Dioxides

Thiophene S,S-dioxides are underutilized tools for the de novo construction of benzene rings in organic synthesis. We report a collective synthesis of nine illudalane sesquiterpenes using bicyclic thiophene S,S-dioxides as generalized precursors to the indane core of the natural products. Exploiting furans as unusual dienophiles in this inverse electron demand Diels-Alder cascade, this concise and convergent approach enables the synthesis of these targets in as little as five steps. Theoretical studies rationalize the reactivity of thiophene S,S-dioxides with both electron-poor and electron-rich dienophiles and reveal reaction pathways involving either nonpolar pericyclic or bifurcating ambimodal cycloadditions. Overall, this work demonstrates the wider potential of thiophene S,S-dioxides as convenient and flexible precursors to polysubstituted arenes.

Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics

Spin-state energetics of transition metal complexes remain one of the most challenging targets for electronic structure methods. Among single-reference wave function approaches, local correlation approximations to coupled cluster theory, most notably the domain-based local pair natural orbital (DLPNO) approach, hold the promise of bringing the accuracy of coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T), to molecular systems of realistic size with acceptable computational cost. However, recent studies on spin-state energetics of iron-containing systems raised doubts about the ability of the DLPNO approach to adequately and systematically approximate energetics obtained by the reference-quality complete active space second-order perturbation theory with coupled-cluster semicore correlation, CASPT2/CC. Here, we revisit this problem using a diverse set of iron complexes and examine several aspects of the application of the DLPNO approach. We show that DLPNO-CCSD(T) can accurately reproduce both CASPT2/CC and canonical CCSD(T) results if two basic principles are followed. These include the consistent use of the improved iterative (T₁) versus the semicanonical perturbative triple corrections and, most importantly, a simple two-point extrapolation to the PNO space limit. The latter practically eliminates errors arising from the default truncation of electron-pair correlation spaces and should be viewed as standard practice in applications of the method to transition metal spin-state energetics. Our results show that reference-quality results can be readily achieved with DLPNO-CCSD(T) if these principles are followed. This is important also in view of the applicability of the method to larger single-reference systems and multinuclear clusters, whose treatment of dynamic correlation would be challenging for multireference-based approaches.

Optical absorption and shape transition in neutral SnN clusters with N ≤ 40: a photodissociation spectroscopy and electric beam deflection study

Neutral SnN clusters with N = 6-20, 25, 30, 40 are investigated in a joint experimental and quantum chemical study with the aim to reveal their optical absorption in conjunction with their structural evolution. Electric beam deflection and photodissociation spectroscopy are applied as molecular beam techniques at nozzle temperatures of 16 K, 32 K and 300 K. The dielectric response is probed following the approach in S. Schäfer et al., J. Phys. Chem A, 2008, 112, 12312-12319. It is improved on those findings and the cluster size range is extended in order to cover the prolate growth regime. The impact of the electric dipole moment, rotational temperature and vibrational excitation on the deflection profiles is discussed thoroughly. Photodissociation spectra of tin clusters are recorded for the first time, show similarities to spectra of silicon clusters and are demonstrated to be significantly complicated by the presence of multiphoton absorption in the low-energy region and large excess energies upon dissociation which is modelled by the RRKM theory. In both experiments two isomers for the clusters with N = 8, 11, 12, 19 need to be considered to explain the experimental results. Triple-capped trigonal prisms and double-capped square antiprisms are confirmed to be the driving building units for almost the entire size range. Three dominating fragmentation channels are observed, i.e. the loss of a tin atom for N < 12, a Sn7 fragment for N < 19 and a Sn10 fragment for N ≥ 19 with Sn15 subunits constituting recurring geometric motifs for N > 20. The prolate-to-quasispherical structural transition is found to occur at 30 < N ≤ 40 and is analyzed with respect to the observed optical behavior taking quantum chemical calculations and the Mie-Gans theory into account. Limitations of the experimental approach to study the geometric and electronic structure of the clusters at elevated temperatures due to vibrational excitation is also thoroughly discussed.

Probing the Transition State-to-Intermediate Continuum: Mechanistic Distinction between a Dry versus Wet Perepoxide in the Singlet Oxygen "Ene" Reaction at the Air-Water Interface

A mechanistic study is reported for the reactions of singlet oxygen (¹O₂) with alkene surfactants of tunable properties. Singlet oxygen was generated either top-down (photochemically) by delivery as a gas to an air-water interface or bottom-up (chemically) by transport to the air-water interface as a solvated species. In both cases, reactions were carried out in the presence of 7-carbon (7C), 9-carbon (9C), or 11-carbon (11C) prenylsurfactants [(CH₃)₂C═CH(CH₂)_nSO₃^- Na⁺ (n = 4, 6, 8)]. Higher "ene" hydroperoxide regioselectivities (secondary ROOH 2 to tertiary ROOH 3) were reached in delivering ¹O₂ top-down through air as compared to bottom-up via aqueous solution. In the photochemical reaction, ratios of 2:3 increased from 2.5:1 for 7C, to 2.8:1 for 9C, and to 3.2:1 for 11C. In contrast, in the bubbling system that generated ¹O₂ chemically, the selectivity was all but lost, ranging only from 1.3:1 to 1:1. The phase-dependent regioselectivities appear to be correlated with the "ene" reaction with photochemically generated, drier ¹O₂ at the air-water interface vs those with wetter ¹O₂ from the bubbling reactor. Density functional theory-calculated reaction potential energy surfaces (PESs) were used to help rationalize the reaction phase dependence. The reactions in the gas phase are mediated by perepoxide transition states with 32-41 kJ/mol binding energy for C═C(π)···¹O₂. The perepoxide species, however, evolve to well-defined stationary structures in the aqueous phase, with covalent C-O bonds and 85-88 kJ/mol binding energy. The combined experimental and computational evidence points to a unique mechanism for ¹O₂ "ene" tunability in a perepoxide continuum from a transition state to an intermediate.

Mechanisms, Challenges, and Opportunities of Dual Ni/Photoredox-Catalyzed C(sp2)-C(sp3) Cross-Couplings

The merging of photoredox and nickel catalysis has revolutionized the field of C-C cross-coupling. However, in comparison to the development of synthetic methods, detailed mechanistic investigations of these catalytic systems are lagging. To improve the mechanistic understanding, computational tools have emerged as powerful tools to elucidate the factors controlling reactivity and selectivity in these complex catalytic transformations. Based on the reported computational studies, it appears that the mechanistic picture of catalytic systems is not generally applicable, but is rather dependent on the specific choice of substrate, ligands, photocatalysts, etc. Given the complexity of these systems, the need for more accurate computational methods, readily available and user-friendly dynamics simulation tools, and data-driven approaches is clear in order to understand at the molecular level the mechanisms of these transformations. In particular, we anticipate that such improvement of theoretical methods will become crucial to advance the understanding of excited-state properties and dynamics of key species, as well as to enable faster and unbiased exploration of reaction pathways. Further, with greater collaboration between computational, experimental, and spectroscopic communities, the mechanistic investigation of photoredox/Ni dual-catalytic reactions is expected to thrive quickly, facilitating the design of novel catalytic systems and promoting our understanding of the reaction selectivity.

Double Hybrids and Noncovalent Interactions: How Far Can We Go?

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.

Ab initio development of generalized Lennard-Jones (Mie) force fields for predictions of thermodynamic properties in advanced molecular-based SAFT equations of state

A methodology for obtaining molecular parameters of a modified statistical associating fluid theory for variable-range interactions of Mie form (SAFT-VR Mie) equation of state (EoS) from ab initio calculations is proposed for non-associative species that can be modeled as single spherical segments. The methodology provides a strategy to map interatomic or intermolecular potentials obtained from ab initio quantum-chemistry calculations to the corresponding Mie potentials that can be used within the SAFT-VR Mie EoS. The inclusion of corrections for quantum and many-body effects allows for an excellent, fully predictive description of the vapor-liquid envelope and other bulk thermodynamic properties of noble gases; this description is of similar or superior quality to that obtained using SAFT-VR Mie with parameters regressed in the traditional way using experimental thermodynamic-property data. The methodology is extended to an anisotropic species, methane, where similar levels of accuracy are obtained. The efficacy of using less-accurate quantum-chemistry methods in this methodology is explored, showing that these methods do not provide satisfactory results, although we note that the description is nevertheless substantially better than those obtained using the conductor-like screening model for describing real solvents (COSMO-RS), the only other fully predictive ab initio method currently available. Overall, the reliance on thermophysical data is completely dispensed with, providing the first extensible, wholly predictive SAFT-type EoSs.

Collision-induced dissociation of homodimeric and heterodimeric radical cations of 9-methylguanine and 9-methyl-8-oxoguanine: correlation between intra-base pair proton transfer originating from the N1-H at a Watson-Crick edge and non-statistical dissociation

It has been shown previously in protonated, deprotonated and ionized guanine-cytosine base pairs that intra-base pair proton transfer from the N1-H at the Watson-Crick edge of guanine to the complementary nucleobase prompts non-statistical dissociation of the base-pair system, and the dissociation of a proton-transferred base-pair structure is kinetically more favored than that of the starting, conventional base-pair structure. However, the fundamental chemistry underlying this anomalous and intriguing kinetics has not been completely revealed, which warrants the examination of more base-pair systems in different structural contexts in order to derive a generalized base-pair structure-kinetics correlation. The purpose of the present work is to expand the investigation to the non-canonical homodimeric and heterodimeric radical cations of 9-methylguanine (9MG) and 9-methyl-8-oxoguanine (9MOG), i.e., [9MG·9MG]˙⁺, [9MOG·9MG]˙⁺ and [9MOG·9MOG]˙⁺. Experimentally, collision-induced dissociation tandem mass spectrometry coupled with an electrospray ionization (ESI) source was used for the formation of base-pair radical cations, followed by detection of dissociation product ions and cross sections in the collisions with Xe gas under single ion-molecule collision conditions and as a function of the center-of-mass collision energy. Computationally, density functional theory and coupled cluster theory were used to calculate and identify probable base-pair structures and intra-base pair proton transfer and hydrogen transfer reactions, followed by kinetics modeling to explore the properties of dissociation transition states and kinetic factors. The significance of this work is twofold: it provides insight into base-pair opening kinetics in three biologically-important, non-canonical systems upon oxidative and ionization damage; and it links non-statistical dissociation to intra-base pair proton-transfer originating from the N1-H at the Watson-Crick edge of 8-oxoguanine, enhancing understanding towards the base-pair fragmentation assisted by proton transfer.

Open-Shell Variant of the London Dispersion-Corrected Hartree-Fock Method (HFLD) for the Quantification and Analysis of Noncovalent Interaction Energies

The London dispersion (LD)-corrected Hartree–Fock (HF) method (HFLD) is an ab initio approach for the quantification and analysis of noncovalent interactions (NCIs) in large systems that is based on the domain-based local pair natural orbital coupled-cluster (DLPNO-CC) theory. In the original HFLD paper, we discussed the implementation, accuracy, and efficiency of its closed-shell variant. Herein, an extension of this method to open-shell molecular systems is presented. Its accuracy is tested on challenging benchmark sets for NCIs, using CCSD(T) energies at the estimated complete basis set limit as reference. The HFLD scheme was found to be as accurate as the best-performing dispersion-corrected exchange-correlation functionals, while being nonempirical and equally efficient. In addition, it can be combined with the well-established local energy decomposition (LED) for the analysis of NCIs, thus yielding additional physical insights.

Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals

The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T₀) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T₀) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T₀). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.

Hydrogen-Bond Topology Is More Important Than Acid/Base Strength in Atmospheric Prenucleation Clusters

We explored the hypothesis that on the nanoscale level, acids and bases might exhibit different behavior than in bulk solution. Our study system consisted of sulfuric acid, formic acid, ammonia, and water. We calculated highly accurate Domain-based Local pair-Natural Orbital- Coupled-Cluster/Complete Basis Set (DLPNO-CCSD(T)/CBS) energies on DFT geometries and used the resulting Gibbs free energies for cluster formation to compute the overall equilibrium constants for every possible cluster. The equilibrium constants combined with the initial monomer concentrations were used to predict the formation of clusters at the top and the bottom of the troposphere. Our results show that formic acid is as effective as ammonia at forming clusters with sulfuric acid and water. The structure of formic acid is uniquely suited to form hydrogen bonds with sulfuric acid. Additionally, it can partner with water to form bridges from one side of sulfuric acid to the other, hence demonstrating that hydrogen bonding topology is more important than acid/base strength in these atmospheric prenucleation clusters.

Analysis of Conformational Preferences in Caffeine

High level DLPNO−CCSD(T) electronic structure calculations with extended basis sets over B3LYP−D3 optimized geometries indicate that the three methyl groups in caffeine overcome steric hindrance to adopt uncommon conformations, each one placing a C−H bond on the same plane of the aromatic system, leading to the C−H bonds eclipsing one carbonyl group, one heavily delocalized C−N bond constituent of the fused double ring aromatic system, and one C−H bond from the imidazole ring. Deletion of indiscriminate and selective non-Lewis orbitals unequivocally show that hyperconjugation in the form of a bidirectional −CH3 ⇆ aromatic system charge transfer is responsible for these puzzling conformations. The structural preferences in caffeine are exclusively determined by orbital interactions, ruling out electrostatics, induction, bond critical points, and density redistribution because the steric effect, the allylic effect, the Quantum Theory of Atoms in Molecules (QTAIM), and the non-covalent interactions (NCI), all predict wrong energetic orderings. Tiny rotational barriers, not exceeding 1.3 kcal/mol suggest that at room conditions, each methyl group either acts as a free rotor or adopts fluxional behavior, thus preventing accurate determination of their conformations. In this context, our results supersede current experimental ambiguity in the assignation of methyl conformation in caffeine and, more generally, in methylated xanthines and their derivatives.

Gauging the Steric Effects of Silyl Groups with a Molecular Balance

We present an experimental and computational study of a cyclooctatetraene (COT)-based molecular balance disubstituted with commonly used silyl groups. Such groups often serve as protecting groups and are typically considered innocent bystanders. Our motivation here is to determine the actual steric effects of such groups by employing a molecular balance. While in the unfolded 1,4-valence isomer the silyl groups are far apart (d_σ-σ ≥ 5.15 Å), the folded 1,6-isomer is affected greatly by noncovalent interactions due to close σ-σ contacts (d_σ-σ ≤ 2.58 Å). In order to investigate the thermodynamic equilibrium between the 1,6- and 1,4-valence isomers, we employed temperature-dependent nuclear magnetic resonance measurements. Additionally, we assessed the nature of attractive and repulsive interactions in 1,6-disilyl-COT derivatives via a combination of local energy decomposition analysis (LED) and symmetry-adapted perturbation theory (SAPT) at the DLPNO-CCSD(T)/def2-TZVP and sSAPT0/aug-cc-pVDZ levels of theory. We identified London dispersion interactions as the main contributor to the molecular stability of the folded states, whereas Pauli exchange repulsion and a resulting internal strain favor the unfolded diastereomer.

A Hierarchical Theoretical Study of the Hydrogen Abstraction Reactions of H2/C1-C4 Molecules by the Methyl Peroxy Radical and Implications for Kinetic Modeling

The hydrogen atom abstraction by the methyl peroxy radical (CH₃O₂) is an important reaction class in detailed chemical kinetic modeling of the autoignition properties of hydrocarbon fuels. Systematic theoretical studies are performed on this reaction class for H₂/C₁-C₄ fuels, which is critical in the development of a base model for large fuels. The molecules include hydrogen, alkanes, alkenes, and alkynes with a carbon number from 1 to 4. The B2PLYP-D3/cc-pVTZ level of theory is employed to optimize the geometries of all of the reactants, transition states, and products and also the treatments of hindered rotation for lower frequency modes. Accurate benchmark calculations for abstraction reactions of hydrogen, methane, and ethylene with CH₃O₂ are performed by using the coupled cluster method with explicit inclusion of single and double electron excitations and perturbative inclusion of triple electron excitations (CCSD(T)), the domain-based local pair-natural orbital coupled cluster method (DLPNO-CCSD(T)), and the explicitly correlated CCSD(T)-F12 method with large basis sets. Reaction rate constants are computed via conventional transition state theory with quantum tunneling corrections. The computed rate constants are compared with literature values and those employed in detailed chemical kinetic mechanisms. The calculated rate constants are implemented into the recently developed NUIGMECH1.1 base model for kinetic modeling of ignition properties.

The effect of alkylation on the micro-solvation of ethers revealed by highly localized water librational motion

The specific far-infrared spectral signatures associated with highly localized large-amplitude out-of-plane librational motion of water molecules have recently been demonstrated to provide sensitive spectroscopic probes for the micro-solvation of organic molecules [Mihrin et al., Phys. Chem. Chem. Phys. 21(4), 1717 (2019)]. The present work employs this direct far-infrared spectroscopic approach to investigate the non-covalent intermolecular forces involved in the micro-solvation of a selection of seven ether molecules with systematically varied alkyl substituents: dimethyl ether, diethyl ether, diisopropyl ether, ethyl methyl ether, t-butyl methyl ether, and t-butyl ethyl ether. The ranking of the observed out-of-plane water librational band signatures for this selected series of ether–water complexes embedded in inert neon matrices at 4 K reveals information about the interplay of directional intermolecular hydrogen bond motifs and non-directional and long-range dispersion interactions for the micro-solvated structures. These far-infrared observables differentiate minor subtle effects introduced by specific alkyl substituents and serve as rigorous experimental benchmarks for modern quantum chemical methodologies of various levels of scalability, which often fail to accurately predict the structural variations and corresponding vibrational signatures of the closely related systems. The accurate interaction energies of the series of ether–water complexes have been predicted by the domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations, followed by a local energy decomposition analysis of the energy components. In some cases, the secondary dispersion forces are in direct competition with the primary intermolecular hydrogen bonds as witnessed by the specific out-of-plane librational signatures.

Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches

Hydrated transition metal ions are prototypical systems that can be used to model properties of transition metals in complex chemical environments. These seemingly simple systems present challenges for computational chemistry and are thus crucial in evaluations of quantum chemical methods for spin-state and redox energetics. In this work, we explore the applicability of the domain-based pair natural orbital implementation of coupled cluster (DLPNO-CC) theory to the calculation of ionization energies and redox potentials for hydrated ions of all first transition row (3d) metals in the 2+/3+ oxidation states, in connection with various solvation approaches. In terms of model definition, we investigate the construction of a minimally explicitly hydrated quantum cluster with a first and second hydration layer. We report on the convergence with respect to the coupled cluster expansion and the PNO space, as well as on the role of perturbative triple excitations. A recent implementation of the conductor-like polarizable continuum model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent redox potentials at the coupled cluster level. Our results establish conditions for the convergence of DLPNO-CCSD(T) energetics and stress the absolute necessity to explicitly consider the second solvation sphere even when CPCM is used. The achievable accuracy for redox potentials of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore, multilayer approaches that combine a higher-level DLPNO-CCSD(T) description of the first solvation sphere with a lower-level description of the second solvation layer are investigated. The present work establishes optimal and transferable methodological choices for employing DLPNO-based coupled cluster theory, the associated CPCM implementation, and cost-efficient multilayer derivatives of the approach for open-shell transition metal systems in complex environments.

Gem-Diaurated Gold(III) Complexes: Synthesis, Structure, Aurophilic Interaction, and Catalytic Activity

We present a protocol to synthesize air stable gem-diaurated gold(III) compounds from 1,3-diketones in a single cycloauration step with tetrachloroauric acid. So far related species were only accessible from phosphonium bis(ylide) ligands which hold the two gold atoms in close proximity. Lacking such a constraint, our compounds show the longest Au-Au distances of all gem-diaurated carbons, ranging from 3.26 to 3.32 Å. Modeling based on results of CCSD(T) calculations shows no stabilization by aurophilic interactions for our gold(III) systems, compared to 9.1 kcal/mol for gold(I) gem-diauration. This demonstrates no aurophilic interactions are needed for the isolation of air stable gem-diaurated gold(III) complexes. We show the new gem-diaurated gold(III) compounds are active in the gold-catalyzed phenol synthesis and highly active in the cycloisomerization of an N-propargylcarboxamide; here, we obtained the so far highest known TON of over 2500 per gold atom with respect to the oxazole formation.

Stability of the polyynic form of C18, C22, C26, and C30 nanorings: a challenge tackled by range-separated double-hybrid density functionals

We calculate the relative energy between the cumulene and polyyne structures of a set of C_4k+2 (k = 4-7) rings (C₁₈, C₂₂, C₂₆, and C₃₀ prompted by the recent synthesis of the cyclo[18]carbon (or simply C₁₈) compounds. Reference results were obtained by a costly Quantum Monte-Carlo (QMC) approach, providing thus very accurate values allowing to systematically compare the performance of a variety of wavefunction methods [(i.e., MP2, SCS-MP2, SOS-MP2, DLPNO-CCSD, and DLPNO-CCSD(T)] as well as DFT approaches, applying for the latter a diversity of density functionals covering global and range-separated hybrid and double-hybrid models. The influence of the use of a range-separation scheme for density functionals, for both hybrid and double-hybrid expressions, is discussed according to its key role. Overall, range-separated double-hybrid functionals (e.g., RSX-QIDH) behave very accurately and provide competitive results compared with DLPNO-CCSD(T), at a more reasonable computational cost.

The Mechanism and Kinetics Model of Degradation of Dicarboxylic Acids by Hydroxyl Radicals under Atmospheric Conditions

The atmospheric degradation mechanism of dicarboxylic acids (DCAs) initiated by hydroxyl radicals has been theoretically investigated at the DLPNO-CCSD(T)/def2-TZVP//BH&HLYP/6-311++G(d,p) level of theory. In the presence of O₂, the degradation of DCAs by hydroxyl radicals takes place through a two-step mechanism: the α-H elimination and the degradation of the peroxyl radical intermediate. The latter degradation mechanism is easy to proceed for the exothermic process of radical recombination. Therefore, the degradation rate of DCAs is determined by an α-H elimination step, which is accelerated in the case of long carbon-chain DCAs with a lower energy barrier. Canonical variational transition state theory has been employed to estimate the rate constants of the H-elimination step of the DCA degradation reaction by hydroxyl radicals over the temperature range of 220-1000 K.

Gas-phase thermochemistry of polycyclic aromatic hydrocarbons: an approach integrating the quantum chemistry composite scheme and reaction generator

We introduce a protocol aimed at predicting the accurate gas-phase enthalpies of formation of polycyclic aromatic hydrocarbons (PAHs). Automatic generation of a dataset of equilibrated chemical reactions preserving the number of carbon atoms in each hybridization state on each side of equations is at the core of our scheme. The performed tests suggest the recommended enthalpy of formation to be derived via a two-step scheme. First, we consider the reactions with a minimal sum of the total number of particles involved, N, and the absolute difference between the total number of products and reactants, |ΔN|. Second, among these reactions, we identify the one with the smallest absolute reaction enthalpy change, . This approach has been applied to predict the gas-phase enthalpies of formation of 113 PAHs via the Feller-Peterson-Dixon approach. Our calculated values provide the mean absolute deviations of 1.7, 1.9, 4.2, 8.1, and 18.5 kJ mol^-1 with respect to the literature group-based error corrected (GBEC) G3MP2B3, ATOMIC (HC), group equivalent M06-2X, GBEC B3LYP, and G4MP2 values. Our predicted values give the mean signed and mean absolute errors of -7.5 and 12.9 kJ mol^-1 with respect to the experimental enthalpies of formation. The combination of our predicted and the experimental values provide the solid-state enthalpies of formation, , which are not available for a few species. Approaching these values as well as , producing large discrepancies from the experimental side, would be indispensable for testing and further tuning of computational chemistry approaches.

Simple Synthesis of Complex Amines from the Diels–Alder Adducts of (–)-Cytisine

The Diels–Alder reaction of N-benzylcytisine with N-methyl- and N-benzylmaleimides is 100% endo-selective and gives the corresponding syn- and anti-diastereomers in 11–42% isolated yields. The studies of the reaction progress with LCMS and NMR along with detailed quantum chemical calculations revealed that some Diels–Alder adducts are kinetically and their isomers are thermodynamically controlled products. The Pd/C-catalyzed hydrogenation of benzyl-protected cytisine amine derivatives resulted in the removal of the benzyl group and the addition of hydrogen to the C=C double bond to give the corresponding secondary amines in 45–84% yield. The complete reduction of carbonyl groups in a cytisine derivative with LiAlH4 in THF under reflux afforded the respective tricyclic triamine. Quantum mechanical calculations for the mechanism of the Diels–Alder reaction between the simplest model compounds are presented.

The Adsorption of Small Molecules on the Copper Paddle-Wheel: Influence of the Multi-Reference Ground State

We report a theoretical study of the adsorption of a set of small molecules (C₂H₂, CO, CO₂, O₂, H₂O, CH₃OH, C₂H₅OH) on the metal centers of the “copper paddle-wheel”—a key structural motif of many MOFs. A systematic comparison between DFT of different rungs, single-reference post-HF methods (MP2, SOS–MP2, MP3, DLPNO–CCSD(T)), and multi-reference approaches (CASSCF, DCD–CAS(2), NEVPT2) is performed in order to find a methodology that correctly describes the complicated electronic structure of paddle-wheel structure together with a reasonable description of non-covalent interactions. Apart from comparison with literature data (experimental values wherever possible), benchmark calculations with DLPNO–MR–CCSD were also performed. Despite tested methods show qualitative agreement in the majority of cases, we showed and discussed reasons for quantitative differences as well as more fundamental problems of specific cases.

CASPT2 molecular geometries of Fe(II) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.

The MOBH35 Metal–Organic Barrier Heights Reconsidered: Performance of Local-Orbital Coupled Cluster Approaches in Different Static Correlation Regimes

We have revisited the MOBH35 (Metal–Organic Barrier Heights, 35 reactions) benchmark [Iron, Janes, J. Phys. Chem. A, 2019, 123 ( (17), ), 3761−378130973722; ibid. 2019, 123, 6379–6380] for realistic organometallic catalytic reactions, using both canonical CCSD(T) and localized orbital approximations to it. For low levels of static correlation, all of DLPNO-CCSD(T), PNO-LCCSD(T), and LNO-CCSD(T) perform well; for moderately strong levels of static correlation, DLPNO-CCSD(T) and (T₁) may break down catastrophically, and PNO-LCCSD(T) is vulnerable as well. In contrast, LNO-CCSD(T) converges smoothly to the canonical CCSD(T) answer with increasingly tight convergence settings. The only two reactions for which our revised MOBH35 reference values differ substantially from the original ones are reaction 9 and to a lesser extent 8, both involving iron. For the purpose of evaluating density functional theory (DFT) methods for MOBH35, it would be best to remove reaction 9 entirely as its severe level of static correlation makes it just too demanding for a test. The magnitude of the difference between DLPNO-CCSD(T) and DLPNO-CCSD(T₁) is a reasonably good predictor for errors in DLPNO-CCSD(T₁) compared to canonical CCSD(T); otherwise, monitoring all of T₁, D₁, max|t_i^A|, and 1/(ε_LUMO – ε_HOMO) should provide adequate warning for potential problems. Our conclusions are not specific to the def2-SVP basis set but are largely conserved for the larger def2-TZVPP, as they are for the smaller def2-SV(P): the latter may be an economical choice for calibrating against canonical CCSD(T). Finally, diagnostics for static correlation are statistically clustered into groups corresponding to (1) importance of single excitations in the wavefunction; (2a) the small band gap, weakly separated from (2b) correlation entropy; and (3) thermochemical importance of correlation energy, as well as the slope of the DFT reaction energy with respect to the percentage of HF exchange. Finally, a variable reduction analysis reveals that much information on the multireference character is provided by T₁, I_ND/I_tot, and the exchange-based diagnostic A₁₀₀[TPSS].

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

The accuracy of correlation energy recovered by coupled cluster single-, double-, and perturbative triple-excitations, CCSD(T), has led to the method being considered the gold standard of computational chemistry. The application of CCSD(T) has been limited to medium-sized molecular systems due to its steep scaling with molecular size. The recent development of alternative domain-based local pair natural orbital coupled-cluster method, DLPNO-CCSD(T), has significantly broadened the range of chemical systems to which CCSD(T) level calculations can be applied. Condensed systems such as ionic liquids (ILs) have a large contribution from London dispersion forces of up to 150 kJ mol^-1 in large-scale clusters. Ionic liquids show appreciable charge transfer effects that result in the increased valence orbital delocalization over the entire ionic network, raising the question whether the application of methods based on localized orbitals is reliable for these semi-Coulombic materials. Here the performance of DLPNO-CCSD(T) is validated for the prediction of correlation interaction energies of two data sets incorporating single-ion pairs of protic and aprotic ILs. DLPNO-CCSD(T) produced results within chemical accuracy with tight parameter settings and a non-iterative treatment of triple excitations. To achieve spectroscopic accuracy of 1 kJ mol^-1 , especially for hydrogen-bonded ILs and those containing halides, the DLPNO settings had to be increased by two orders of magnitude and include the iterative treatment of triple excitations, resulting in a 2.5-fold increase in computational cost. Two new sets of parameters are put forward to produce the performance of DLPNO-CCSD(T) within chemical and spectroscopic accuracy.

Modulating the reactivity of phosphanylidenephosphoranes towards water with Lewis acids

Phosphanylidenephosphoranes of the type R−P(PR’3), also known as phospha-Wittig reagents, can be utilized in a variety of bond activation reactions exploiting their phosphinidenoid reactivity. In here, we thoroughly show that...

S66x8 noncovalent interactions revisited: new benchmark and performance of composite localized coupled-cluster methods

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).

Gas Phase Protolysis of Trisarylzincate Anions

We have applied a combination of tandem-mass spectrometry, quantum-chemical calculations, and statistical rate theory computations to examine the gas phase reactions between the trisarylzincate anions Ar^XZnPh₂^- (Ar^X = p-X-C₆H₄, X = NMe₂, OMe, Me, H, F, and Cl) and 2,2,2-trifluoroethanol at T = 310 ± 20 K. The observed reactions bring about the protonation of one of the aryl anions, which is then released as the corresponding arene, while the formed alkoxide binds to the zinc center. The protonation is faster for the more electron-rich aryl groups and shows a linear Hammett plot if the rate constant for X = NMe₂ is discarded from the analysis. Although the reactions are highly exothermic, they proceed only with relatively low efficiencies (0.1% ≤ φ ≤ 1.3%). According to the quantum-chemical calculations, this behavior can be ascribed to the reactions proceeding through a double-well potential with a tight transition structure located at the central barrier. Based on these potential energy surfaces, the statistical rate theory computations can reproduce the measured rate constants within factors of 2 to 8. A comparison of the protolysis of the trisarylzincates with that of the corresponding free aryl anions demonstrates how the coordination to the metal center not only stabilizes the carbanions energetically but also moderates their reactivity. Thus, our gas phase study contributes to a better understanding of the fundamentals of organometallic reactivity.

Molecular Structure of Nickel Octamethylporphyrin—Rare Experimental Evidence of a Ruffling Effect in Gas Phase

The structure of a free nickel (II) octamethylporphyrin (NiOMP) molecule was determined for the first time through a combined gas-phase electron diffraction (GED) and mass spectrometry (MS) experiment, as well as through quantum chemical (QC) calculations. Density functional theory (DFT) calculations do not provide an unambiguous answer about the planarity or non-planar distortion of the NiOMP skeleton. The GED refinement in such cases is non-trivial. Several approaches to the inverse problem solution were used. The obtained results allow us to argue that the ruffling effect is manifested in the NiOMP molecule. The minimal critical distance between the central atom of the metal and nitrogen atoms of the coordination cavity that provokes ruffling distortion in metal porphyrins is about 1.96 Å.

Singlet O2 Reactions with Radical Cations of 8-Bromoguanine and 8-Bromoguanosine: Guided-Ion Beam Mass Spectrometric Measurements and Theoretical Treatments

8-Bromoguanosine is generated in vivo as a biomarker for early inflammation. Its formation and secondary reactions lead to a variety of biological sequelae at inflammation sites, most of which are mutagenic and linked to cancer. Herein, we report the formation of radical cations of 8-bromoguanine (8BrG^•+) and 8-bromoguanosine (8BrGuo^•+) and their reactions toward the lowest excited singlet molecular oxygen (¹O₂)─a common reactive oxygen species generated in biological systems. This work aims to investigate synergistic, oxidatively generated damage of 8-brominated guanine and guanosine that may occur upon ionizing radiation, one-electron oxidation, and ¹O₂ oxidation. Capitalizing on measurements of reaction product ions and cross sections of 8BrG^•+ and 8BrGuo^•+ with ¹O₂ using guided-ion beam tandem mass spectrometry and augmented by computational modeling of the prototype reaction system, 8BrG^•+ + ¹O₂, using the approximately spin-projected ωB97XD/6-31+G(d,p) density functional theory, the coupled cluster DLPNO-CCSD(T)/aug-cc-pVTZ and the multireference CASPT2(21,15)/6-31G**, probable reaction products, and potential energy surfaces (PESs) were mapped out. 8BrG^•+ and 8BrGuo^•+ present similar exothermic oxidation products, and their reaction efficiencies with ¹O₂ increase with decreasing collision energy. Both single- and multireference theories predicted that the two most energetically favorable reaction pathways correspond to ¹O₂-addition to the C8 and C5-positions of 8BrG^•+, respectively. The CASPT2-calculated PES represents the best quantitative agreement with the experimental benchmark, in that the oxidation exothermicity is close to the water hydration energy of product ions and, thus, is able to eliminate a water ligand in the product ions.

BH9, a New Comprehensive Benchmark Data Set for Barrier Heights and Reaction Energies: Assessment of Density Functional Approximations and Basis Set Incompleteness Potentials

The calculation of accurate reaction energies and barrier heights is essential in computational studies of reaction mechanisms and thermochemistry. To assess methods regarding their ability to predict these two properties, high-quality benchmark sets are required that comprise a reasonably large and diverse set of organic reactions. Due to the time-consuming nature of both locating transition states and computing accurate reference energies for reactions involving large molecules, previous benchmark sets have been limited in scope, the number of reactions considered, and the size of the reactant and product molecules. Recent advances in coupled-cluster theory, in particular local correlation methods like DLPNO-CCSD(T), now allow the calculation of reaction energies and barrier heights for relatively large systems. In this work, we present a comprehensive and diverse benchmark set of barrier heights and reaction energies based on DLPNO-CCSD(T)/CBS called BH9. BH9 comprises 449 chemical reactions belonging to nine types common in organic chemistry and biochemistry. We examine the accuracy of DLPNO-CCSD(T) vis-a-vis canonical CCSD(T) for a subset of BH9 and conclude that, although there is a penalty in using the DLPNO approximation, the reference data are accurate enough to serve as a benchmark for density functional theory (DFT) methods. We then present two applications of the BH9 set. First, we examine the performance of several density functional approximations commonly used in thermochemical and mechanistic studies. Second, we assess our basis set incompleteness potentials regarding their ability to mitigate basis set incompleteness errors. The number of data points, the diversity of the reactions considered, and the relatively large size of the reactant molecules make BH9 the most comprehensive thermochemical benchmark set to date and a useful tool for the development and assessment of computational methods.

Near-Exact CCSDT Energetics from Rank-Reduced Formalism Supplemented by Non-iterative Corrections

We introduce a non-iterative energy correction, added on top of the rank-reduced coupled-cluster method with single, double, and triple substitutions, that accounts for excitations excluded from the parent triple excitation subspace. The formula for the correction is derived by employing the coupled-cluster Lagrangian formalism, with an additional assumption that the parent excitation subspace is closed under the action of the Fock operator. Owing to the rank-reduced form of the triple excitation amplitudes tensor, the computational cost of evaluating the correction scales as N⁷, where N is the system size. The accuracy and computational efficiency of the proposed method is assessed for both total and relative correlation energies. We show that the non-iterative correction can fulfill two separate roles. If the accuracy level of a fraction of kJ/mol is sufficient for a given system, the correction significantly reduces the dimension of the parent triple excitation subspace needed in the iterative part of the calculations. Simultaneously, it enables reproducing the exact CCSDT results to an accuracy level below 0.1 kJ/mol, with a larger, yet still reasonable, dimension of the parent excitation subspace. This typically can be achieved at a computational cost only several times larger than required for the CCSD(T) method. The proposed method retains the black-box features of the single-reference coupled-cluster theory; the dimension of the parent excitation subspace remains the only additional parameter that has to be specified.

Direct orbital selection within the domain-based local pair natural orbital coupled-cluster method

Domain-based local pair natural orbital coupled cluster (DLPNO-CC) has become increasingly popular to calculate relative energies (e.g., reaction energies and reaction barriers). It can be applied within a multi-level DLPNO-CC-in-DLPNO-CC ansatz to reduce the computational cost and focus the available computational resources on a specific subset of the occupied orbitals. We demonstrate how this multi-level DLPNO-CC ansatz can be combined with our direct orbital selection (DOS) approach [M. Bensberg and J. Neugebauer, J. Chem. Phys. 150, 214106 (2019)] to automatically select orbital sets for any multi-level calculation. We find that the parameters for the DOS procedure can be chosen conservatively such that they are transferable between reactions. The resulting automatic multi-level DLPNO-CC method requires no user input and is extremely robust and accurate. The computational cost is easily reduced by a factor of 3 without sacrificing accuracy. We demonstrate the accuracy of the method for a total of 61 reactions containing up to 174 atoms and use it to predict the relative stability of conformers of a Ru-based catalyst.

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol^-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.