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

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.

Addressing the System-Size Dependence of the Local Approximation Error in Coupled-Cluster Calculations

Over the last two decades, the local approximation has been successfully used to extend the range of applicability of the "gold standard" singles and doubles coupled-cluster method with perturbative triples CCSD(T) to systems with hundreds of atoms. The local approximation error grows in absolute value with the increasing system size, i.e., by increasing the number of electron pairs in the system. In this study, we demonstrate that the recently introduced two-point extrapolation scheme for approaching the complete pair natural orbital (PNOs) space limit in domain-based pair natural orbital CCSD(T) calculations drastically reduces the dependence of the error on the system size, thus opening up unprecedented opportunities for the calculation of benchmark quality relative energies for large systems.

Another look at energetically quasi-degenerate structures of the gold cluster Au27 q with q = 1, 0, -1

Quantum chemical computations were used to reinvestigate the geometries, spectroscopic, and energetic properties of the gold clusters Au₂₇ ^q in three charge states (q = 1, 0, -1). Density functional theory (DFT) and the domain-based local pair natural orbital modification of the coupled-cluster theory DLPNO-CCSD(T) calculations revealed that, at variance with earlier reports in the literature, while the anion Au₂₇ ^- tends to exist in a tube-like form, both the lowest-energy Au₂₇ and Au₂₇ ⁺ isomers exhibit a pyramidal shape. However, several isomers were found to lie very close in energy, thus rendering a structural transition and their coexistence easy to occur. More specifically, the equilibrium geometry of the neutral Au₂₇ is a core-shell pyramid-like structure with one gold atom located inside. We also identified a novel ground state for the anion Au₂₇ ^- and located for the first time the global minimum of the cation Au₂₇ ⁺ . The vertical detachment energies of the neutral and anionic states were also computed and used to assign the available experimental photoelectron spectra. Although many Au₂₇ isomers were predicted to be energetically quasi-degenerate, the corresponding distinctive vibrational signatures can be used as fingerprints for the identification of cluster geometrical features.

The Catalytic Mechanism of Pdx2 Glutaminase Driven by a Cys-His-Glu Triad: A Computational Study

The catalytic mechanism of Pdx2 was studied with atomic detail employing the computational ONIOM hybrid QM/MM methodology. Pdx2 employs a Cys-His-Glu catalytic triad to deaminate glutamine to glutamate and ammonia - the source of the nitrogen of pyridoxal 5'-phosphate (PLP). This enzyme is, therefore, a rate-limiting step in the PLP biosynthetic pathway of Malaria and Tuberculosis pathogens that rely on this mechanism to obtain PLP. For this reason, Pdx2 is considered a novel and promising drug target to treat these diseases. The results obtained show that the catalytic mechanism of Pdx2 occurs in six steps that can be divided into four stages: (i) activation of Cys₈₇ , (ii) deamination of glutamine with the formation of the glutamyl-thioester intermediate, (iii) hydrolysis of the formed intermediate, and (iv) enzymatic turnover. The kinetic data available in the literature (19.1-19.5 kcal mol^-1 ) agree very well with the calculated free energy barrier of the hydrolytic step (18.2 kcal.mol^-11 ), which is the rate-limiting step of the catalytic process when substrate is readily available in the active site. This catalytic mechanism differs from other known amidases in three main points: i) it requires the activation of the nucleophile Cys₈₇ to a thiolate; ii) the hydrolysis occurs in a single step and therefore does not require the formation of a second tetrahedral reaction intermediate, as it is proposed, and iii) Glu₁₉₈ does not have a direct role in the catalytic process. Together, these results can be used for the synthesis of new transition state analogue inhibitors capable of inhibiting Pdx2 and impair diseases like Malaria and Tuberculosis.

A Phosphorus-Based Pacman Dication Generated by Cooperative Self-Activation of a Pacman Phosphane

Formal coordination of phosphorus(III) by a calix[4]pyrrole Schiff base ligand was achieved through the reaction of this ligand with PCl₃ under basic conditions. The reaction product adopts a Pacman conformation with two P-Cl moieties, one in exo and one in endo position. It represents the first non-metal compound of calix[4]pyrrole Schiff base ligands and of Pacman ligands in general. The spatial neighborhood of the two phosphorus atoms enables cooperative reactions. As a first example, the chloride abstraction with AgOTf is presented, yielding a macrocyclic dication with two embedded phosphorus(III) monocations, which both undergo a cooperative, internal activation reaction with an adjacent C=N double bond. This intramolecular redox process affords two pentacoordinated phosphorus(V) centers within the Pacman dication. All reaction products were fully characterized and all results are supported by computations.

Assessment of the van der Waals, Hubbard U parameter and spin-orbit coupling corrections on the 2D/3D structures from metal gold congeners clusters

The coinage-metal clusters possess a natural complexity in their theoretical treatment that may be accompanied by inherent shortcomings in the methodological approach. Herein, we performed a scalar-relativistic density functional theory study, considering Perdew, Burke, and Ernzerhof (PBE) with (empirical and semi empirical) van der Waals (vdW), spin-orbit coupling (SOC), +U (Hubbard term), and their combinations, to treat the Cu 13 , Ag 13 , and Au 13 clusters in different structural motifs. The energetic scenario is given by the confirmation of the 3D lowest energy configurations for Cu 13 and Ag 13 within all approaches, while for Au 13 there is a 2D/3D competition, depending on the applied correction. The 2D geometry is 0.43 eV more stable with plain PBE than the 3D one, the SOC, +U, and/or vdW inclusion decreases the overestimated stability of the planar configurations, where the most surprising result is found by the D3 and D3BJ vdW corrections, for which the 3D configuration is 0.29 and 0.11 eV, respectively, more stable than the 2D geometry (with even higher values when SOC and/or +U are added). The D3 dispersion correction represents 7.9% (4.4%) of the total binding energy for the 3D (2D) configuration, (not) being enough to change the sd hybridization and the position of the occupied d -states. Our predictions are in agreement with experimental results and in line with the best results obtained for bulk systems, as well as with hybrid functionals within D3 corrections. The properties description undergoes small corrections with the different approaches, where general trends are maintained, that is, the average bond length is smaller (larger) for lower (higher)-coordinated structures, since a same number of electrons are shared by a smaller (larger) number of bonds, consequently, the bonds are stronger (weaker) and shorter (longer) and the sd hybridization index is larger (smaller). Thus, Au has a distinct behavior in relation to its lighter congeners, with a complex potential energy surface, where in addition to the relevant relativistic effects, correlation and dispersion effects must also be considered.

Reaction mechanism and dynamics for C8-hydroxylation of 9-methylguanine radical cation by water molecules

In contrast to their spontaneous deprotonation in aqueous solution, reactions of guanine and guanosine radical cations with water in the gas phase are exclusively initiated by hydration of the radical cations as reported in recent work (Y. Sun et al., Phys. Chem. Chem. Phys., 2018, 20, 27510). As gas-phase hydration reactions closely mimic the actual scenario for guanine radical cations in double-stranded DNA, exploration of subsequent reactions within their water complexes can provide an insight into the resulting oxidative damage to nucleosides. Herein guided-ion beam mass spectrometry experiment and direct dynamics trajectory simulations were carried out to examine prototype complexes of the 9-methylguanine radical cation with one and two water ligands (i.e., 9MG˙⁺·(H₂O)_1-2) in the gas phase, wherein the complexes were activated by collisional activation in the experiment and by thermal excitation at high temperatures in the simulations. Guided by mass spectroscopic measurements, trajectory results and reaction potential energy surface, three reaction pathways were identified. The first two reaction pathways start with H-atom abstraction from water by the O6 and N7 atoms in 9MG˙⁺ and are referred to as HA_O6 and HA_N7, respectively. The primary products of HA_O6 and HA_N7 reactions, including [9MG + H_O6]⁺/[9MG + H_N7]⁺ and ˙OH, react further to either form [8OH-9MG + H_O6]˙⁺ and [8OH-9MG + H_N7]˙⁺via C8-hydroxylation or form radical cations of 6-enol-guanine (6-enol-G˙⁺) and 7H-guanine (7HG˙⁺) via S_N2-type methanol elimination. The third reaction pathway corresponds to the formation of 8OH-9MG⁺ by H elimination from the complex, referred to as HE. Among these product channels, [8OH-9MG + H_N7]˙⁺ has the most favorable formation probability, especially in the presence of additional water molecules. This product may serve as a preceding structure to the 8-oxo-7,8-dihydroguanine lesion in DNA and has implications for health effects of radiation exposure and radiation therapy.

Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls

We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency. Most tested (meta)-GGA and hybrid density functionals perform better for open-shell than for closed-shell complexes; this is not the case for the double hybrids considered.

Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls

We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency. Most tested (meta)-GGA and hybrid density functionals perform better for open-shell than for closed-shell complexes; this is not the case for the double hybrids considered.

Assessing Density Functional Theory for Chemically Relevant Open-Shell Transition Metal Reactions

Due to the principle lack of systematic improvement possibilities of density functional theory, careful assessment of the performance of density functional approximations (DFAs) on well-designed benchmark sets, for example, for reaction energies and barrier heights, is crucial. While main-group chemistry is well covered by several available sets, benchmark data for transition metal chemistry is sparse. This is especially the case for larger, chemically relevant molecules. Addressing this issue, we recently introduced the MOR41 benchmark which covers chemically relevant reactions of closed-shell complexes. In this work, we extend these efforts to single-reference open-shell systems and introduce the "reactions of open-shell single-reference transition metal complexes" (ROST61) benchmark set. ROST61 includes accurate coupled-cluster reference values for 61 reaction energies with a mean reaction energy of -42.8 kcal mol^-1. Complexes with 13-93 atoms covering 20 d-block elements are included, but due to the restriction to single-reference open-shell systems, important elements such as iron or platinum could not be taken into account, or only to a small extent. We assess the performance of 31 DFAs in combination with three London dispersion (LD) correction schemes. Further, DFT-based composite methods, MP2, and a few semiempirical quantum chemical methods are evaluated. Consistent with the results for the MOR41 closed-shell benchmark, we find that the ordering of DFAs according to Jacob's ladder is preserved and that adding an LD correction is crucial, clearly improving almost all tested methods. The recently introduced r²SCAN-3c composite method stands out with a remarkable mean absolute deviation (MAD) of only 2.9 kcal mol^-1, which is surpassed only by hybrid DFAs with low amounts of Fock exchange (e.g., 2.3 kcal mol^-1 for TPSS0-D4/def2-QZVPP) and double-hybrid (DH) DFAs but at a significantly higher computational cost. The lowest MAD of only 1.6 kcal mol^-1 is obtained with the DH DFA PWPB95-D4 in the def2-QZVPP basis set approaching the estimated accuracy of the reference method. Overall, the ROST61 set adds important reference data to a sparsely sampled but practically relevant area of chemistry. At this point, it provides valuable orientation for the application and development of new DFAs and electronic structure methods in general.

Contrasting the Mechanism of H2 Activation by Monomeric and Potassium-Stabilized Dimeric AlI Complexes: Do Potassium Atoms Exert any Cooperative Effect?

Aluminyl anions are low-valent, anionic, and carbenoid aluminum species commonly found stabilized with potassium cations from the reaction of Al-halogen precursors and alkali compounds. These systems are very reactive toward the activation of σ-bonds and in reactions with electrophiles. Various research groups have detected that the potassium atoms play a stabilization role via electrostatic and cation ⋯ π interactions with nearby (aromatic)-carbocyclic rings from both the ligand and from the reaction with unsaturated substrates. Since stabilizing K⋯H bonds are witnessed in the activation of this class of molecules, we aim to unveil the role of these metals in the activation of the smaller and less polarizable H₂ molecule, together with a comprehensive characterization of the reaction mechanism. In this work, the activation of H₂ utilizing a NON-xanthene-Al dimer, [K{Al(NON)}]₂ (D) and monomeric, [Al(NON)]^- (M) complexes are studied using density functional theory and high-level coupled-cluster theory to reveal the potential role of K⁺ atoms during the activation of this gas. Furthermore, we aim to reveal whether D is more reactive than M (or vice versa), or if complicity between the two monomer units exits within the D complex toward the activation of H₂ . The results suggest that activation energies using the dimeric and monomeric complexes were found to be very close (around 33 kcal mol^-1 ). However, a partition of activation energies unveiled that the nature of the energy barriers for the monomeric and dimeric complexes are inherently different. The former is dominated by a more substantial distortion of the reactants (and increased interaction energies between them). Interestingly, during the oxidative addition, the distortion of the Al complex is minimal, while H₂ distorts the most, usually over 0.77 Δ E d i s t ≠ . Overall, it is found here that electrostatic and induction energies between the complexes and H₂ are the main stabilizing components up to the respective transition states. The results suggest that the K⁺ atoms act as stabilizers of the dimeric structure, and their cooperative role on the reaction mechanism may be negligible, acting as mere spectators in the activation of H₂ . Cooperation between the two monomers in D is lacking, and therefore the subsequent activation of H₂ is wholly disengaged.

Local energy decomposition analysis and molecular properties of encapsulated methane in fullerene (CH4@C60)

Methane has been successfully encapsulated within cages of C₆₀ fullerene, which is an appropriate model system to study confinement effects. Its chemistry and physics are also relevant for theoretical model descriptions. Here we provide insights into intermolecular interactions and predicted spectroscopic responses of the CH₄@C₆₀ complex and compared them with results from other methods and with data from the literature. Local energy decomposition analysis (LED) within the domain-based local pair natural orbital coupled cluster singles, doubles, and perturbative triples (DLPNO-CCSD(T)) framework was used, and an efficient protocol for studies of endohedral complexes of fullerenes is proposed. This approach allowed us to assess energies in relation to electronic and geometric preparation, electrostatics, exchange, and London dispersion for the CH₄@C₆₀ endohedral complex. The calculated stabilization energy of CH₄ inside the C₆₀ fullerene was -13.5 kcal mol^-1 and its magnitude was significantly larger than the latent heat of evaporation of CH₄. Evaluation of vibrational frequencies and polarizabilities of the CH₄@C₆₀ complex revealed that the infrared (IR) and Raman bands of the endohedral CH₄ were essentially "silent" due to the dielectric screening effect of C₆₀, which acted as a molecular Faraday cage. Absorption spectra in the UV-vis domain and ionization potentials of C₆₀ and CH₄@C₆₀ were predicted. They were almost identical. The calculated ¹H/¹³C NMR shifts and spin-spin coupling constants were in very good agreement with experimental data. In addition, reference DLPNO-CCSD(T) interaction energies for complexes with noble gases (Ng@C₆₀; Ng = He, Ne, Ar, Kr) were calculated. The values were compared with those derived from supramolecular MP2/SCS-MP2 calculations and estimates with London-type formulas by Pyykkö and coworkers [Phys. Chem. Chem. Phys., 2010, 12, 6187-6203], and with values derived from DFT-based symmetry-adapted perturbation theory (DFT-SAPT) by Hesselmann & Korona [Phys. Chem. Chem. Phys., 2011, 13, 732-743]. Selected points at the potential energy surface of the endohedral He₂@C₆₀ trimer were considered. In contrast to previous theoretical attempts with the DFT/MP2/SCS-MP2/DFT-SAPT methods, our calculations at the DLPNO-CCSD(T) level of theory predicted the He₂@C₆₀ trimer to be thermodynamically stable, which is in agreement with experimental observations.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

Magnetic exchange coupling in Cu dimers studied with modern multireference methods and broken-symmetry coupled cluster theory

Spin-state energetics of exchange-coupled copper complexes pose a persistent challenge for applied quantum chemistry. Here, we provide a comprehensive comparison of all available theoretical approaches to the problem of exchange coupling in two antiferromagnetically coupled bis-μ-hydroxo Cu(II) dimers. The evaluated methods include multireference methods based on the density matrix renormalization group (DMRG), multireference methods that incorporate dynamic electron correlation either perturbatively, such as the N-electron valence state perturbation theory, or variationally, such as the difference-dedicated configuration interaction. In addition, we contrast the multireference results with those obtained using broken-symmetry approaches that utilize either density functional theory or, as demonstrated here for the first time in such systems, a local implementation of coupled cluster theory. The results show that the spin-state energetics of these copper dimers are dominated by dynamic electron correlation and represent an impossible challenge for multireference methods that rely on brute-force expansion of the active space to recover correlation energy. Therefore, DMRG-based methods even at the limit of their applicability cannot describe quantitatively the antiferromagnetic exchange coupling in these dimers, in contrast to dinuclear complexes of earlier transition metal ions. The convergence of the broken-symmetry coupled cluster approach is studied and shown to be a limiting factor for the practical application of the method. The advantages and disadvantages of all approaches are discussed, and recommendations are made for future developments.

Efficient Ab Initio Estimation of Formation Enthalpies for Organic Compounds: Extension to Sulfur and Critical Evaluation of Experimental Data

The efficient protocol for the estimation of gas-phase enthalpies of formation developed previously for C, H, O, N, and F elements was extended to sulfur. The protocol is based on a local coupled cluster with single, double, and perturbative triple excitation [CCSD(T)] approximation and allows rapid evaluation of compounds with sizes computationally prohibitive to canonical CCSD(T) using quadruple zeta basis sets. As a part of model development, a comprehensive review and critical evaluation of experimental data were performed for 87 sulfur-containing organic and inorganic compounds. A compact model with only three empirical parameters for sulfur introduced to address the effects beyond frozen core CCSD(T) was developed. The model exhibits approximately 2 kJ·mol^-1 standard deviation over a set of experimental values for a diverse collection of sulfur-containing compounds. The complete basis set version of the model demonstrates a similar performance and requires only one empirical parameter. Multiple problems with the existing experimental data were identified and discussed. In addition, a lack of reliable data for certain important classes of sulfur compounds was found to impede the model generalization and confident performance assessment.

Redox Isomerism in the S3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory

The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S3 state. Our results enabled us to evaluate conflicting models for the S3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems.

DFT Functionals for Modeling of Polyethylene Chains Cross-Linked by Metal Atoms. DLPNO-CCSD(T) Benchmark Calculations

Density functional theory (DFT) functionals for calculations of binding energies (BEs) of the polyethylene (PE) chains cross-linked by selected metal atoms (M) are benchmarked against DLPNO-CCSD(T) and DLPNO-CCSD(T1) data. PEX-M-PEX complexes as compared with plain parallel PEX···PEX chains with X = 3-9 carbon atoms are model species characterized by a cooperative effect of covalent C-M-C bonds and interchain dispersion interactions. The accuracy of DLPNO-CC methods was assessed by a comparison of BEs with the canonical CCSD(T) results for small PE3-M-PE3 complexes. Functionals for PEX···PEX and closed-shell PEX-M-PEX complexes (M = Be, Mg, Zn) were benchmarked against DLPNO-CCSD(T) BEs; open-shell complexes (M = Li, Ag, Au) were benchmarked against the DLPNO-CCSD(T1) method with iterative triples. Three dispersion corrections were combined with 25 DFT functionals for calculations of BEs with respect to PEX-M and PEX fragments employing def2-TZVPP and def2-QZVPP basis sets. Accuracy to within 5% for the closed-shell PEX-M-PEX complexes was achieved with five functionals. Less accurate are functionals for the open-shell PEX-M-PEX complexes; only two functionals deviate by less than 15% from DLPNO-CCSD(T1). Particularly problematic were PEX-Li-PEX complexes. A reasonable overall performance across all complexes in terms of the mean absolute percentage error is found for the range-separated hybrid functionals ωB97X-D3 and CAM-B3LYP/D3(BJ)-ABC.

3× Axial vs 3× Equatorial: The ΔGGA Value Is a Robust Computational Measure of Substituent Steric Effects

We define ΔG_GA as the free energy change for the formal equilibrium: [13]G-H + 1-X-adamantane → [13]G-X + adamantane, where [13]G-H is the C₁₃H₂₂ fragment of all-trans graphane with 3-fold symmetry. This compares with a situation where the group X is equatorial to three cyclohexane rings with one where it is axial to three rings. ΔG_GA values vary from 2.9 (CN) to 145.7 kJ mol^-1 (CCl₃), and this wide range means that ΔG can be calculated with confidence. ΔG_GA values for Me, Et, i-Pr, and t-Bu form a regular series, 34.9, 63.3, 101.6, and 142.0, and clearly reflect the steric size of the groups. We propose a model where the six axial hydrogens surrounding X on [13]G-X provide a nearly circular constriction on the substituent close to its point of attachment but which does not extend far above this. We compare these results with A values and with calculations on 2- and 7-substituted [1(2,3)4]pentamantanes. We show that electronic effects on ΔG_GA values are negligible but that they correlate well with computed cone and solid angles subtended by the substituent.

Rationalizing the Mechanism of Peroxyformate Decomposition: Computational Insights To Understand Solvent Influence

The heterolytic decomposition of tert‐butyl peroxyformate to tert‐butanol and carbon dioxide, catalyzed by pyridine, is a long‐known example of a reaction whose kinetics are strongly affected by solvent polarity. From DFT and ab initio methods together with the SMD implicit solvation model, an extension on the formerly accepted mechanism is proposed. This novel proposal involves the formation of a carbonic acid ester intermediate and its further decomposition, through an unreported pyridine‐mediated stepwise route. Computed barriers for this mechanism at DLPNO/CCSD(T)‐def2‐TZVP are in excellent agreement with experimental kinetic data across different solvents. Furthermore, the strong relationships between activation energies, geometric parameters in the transition state and the characteristics of the different solvents are also analyzed in depth.

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA-ellipticine and buckycatcher-fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA-ellipticine intercalation complex, and the buckycatcher-C₆₀ complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller-Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA-ellipticine, and buckycatcher-C₆₀, respectively. Such sizable discrepancy above "chemical accuracy" for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree-Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Tautomeric equilibrium, proton affinity and mass spectrometry fragmentation of flexible hydrogen-bonded precursors and rigid [Formula: see text] fluorescent dyes

The stability of two groups of conformationally locked molecules, similar in topology, but differing only by the type of the bridge rigidifying their structure, is studied. The series of the less-rigid 2-phenacylheterocyclic compounds and their stiff difluoroboranyl derivatives are investigated for the determination of the effect of [Formula: see text]/S/O replacement in a five-membered heterocyclic ring and the presence of a strong electron-donating group on the tautomeric equilibrium, protonation affinity, and fragmentation pattern observed in the structural elucidation by means of mass spectrometry technique. The results of the [Formula: see text]B97X-D/6-311++G(d,p) calculations, the topological analysis of electron density as well as the experimental MS measurements show the importance of the number of heteroatoms, their properties, and location in the molecule for the rational design of the systems of desired stable tautomers or the favorable protonation sites. The obtained data allow for the understanding of the fundamentals of the novel highly fluorescent difluoroborates fragmentation behavior, vital for their structural elucidation with the application of high-resolution tandem mass spectrometry methods.

Synthetic strategy toward ineleganolide: A cautionary tale

We present a case study to demonstrate how complex molecule synthesis can benefit from quantum mechanics (QM) calculations. Theory is applied in two contexts: testing the chemical intuition used in retrosynthetic planning, along with expediting the resolution of unexpected challenges encountered during the course of the synthesis. From a computational lens, we examine retrospectively the strategies employed and the decisions made during our synthetic efforts toward the diterpenoid natural product ineleganolide. Seemingly logical and robust hypotheses are found to be ill-fated after theoretical investigation. Prior knowledge of these issues may have potentially saved valuable time and resources during our synthetic efforts. This cautionary tale suggests that synthetic campaigns can benefit from computational evaluation of synthetic plans.

Antiradical Properties of N-Oxide Surfactants-Two in One

Surfactants are molecules that lower surface or interfacial tension, and thus they are broadly used as detergents, wetting agents, emulsifiers, foaming agents, or dispersants. However, for modern applications, substances that can perform more than one function are desired. In this study we evaluated antioxidant properties of two homological series of N-oxide surfactants: monocephalic 3-(alkanoylamino)propyldimethylamine-N-oxides and dicephalic N,N-bis[3,3′-(dimethylamino)propyl]alkylamide di-N-oxides. Their antiradical properties were tested against stable radicals using electron paramagnetic resonance (EPR) and UV-vis spectroscopy. The experimental investigation was supported by theoretical density functional theory (DFT) and ab initio modeling of the X–H bonds dissociation enthalpies, ionization potentials, and Gibbs free energies for radical scavenging reactions. The evaluation was supplemented with a study of biological activity. We found that the mono- and di-N-oxides are capable of scavenging reactive radicals; however, the dicephalic surfactants are more efficient than their linear analogues.

Pairing double hybrid functionals with a tailored basis set for an accurate thermochemistry of hydrocarbons

A collection of five challenging datasets, including noncovalent interactions, reaction barriers and electronic rearrangements of medium-sized hydrocarbons, has been selected to verify the robustness of double-hybrid functionals used in conjunction with the small DH-SVPD basis set, especially developed for noncovalent interactions. The analysis is completed by other, more standard functionals, for a total of 17 models, including also empirical corrections for dispersion. The obtained results show that the chemical accuracy threshold, that is an error lower than 1.0 kcal mol⁻¹, can be obtained by pairing the nonempirical PBE-QIDH functional with the DH-SVPD basis set, as well as by other semi-empirical functionals, such as DSD-PBEP86, using larger basis sets and empirical corrections. More in general, a significant improvement can be obtained using the DH-SVPD basis set with DHs, without resorting to any empirical corrections. This choice leads to a fast computational protocol that, avoiding any empirical potential, remains on a fully quantum ground.

The pairing of the PBE-QIDH double-hybrid functional with a tailored split-valence basis set leads to a fast computational protocol for the accurate evaluation of hydrocarbon thermochemistry, without resorting to any empirical correction.

Singlet Oxygen Oxidation of the Radical Cations of 8-Oxo-2'-deoxyguanosine and Its 9-Methyl Analogue: Dynamics, Potential Energy Surface, and Products Mediated by C5-O2 -Addition

8-Oxo-2'-deoxyguanosine (OG) is the most common DNA lesion. Notably, OG becomes more susceptible to oxidative damage than the undamaged nucleoside, forming mutagenic products in vivo. Herein the reactions of singlet O₂ with the radical cations of 8-oxo-2'-deoxyguanosine (OG^.+ ) and 9-methyl-8-oxoguanine (9MOG^.+ ) were investigated using ion-molecule scattering mass spectrometry, from which barrierless, exothermic O₂ -addition products were detected for both reaction systems. Corroborated by static reaction potential energy surface constructed using multi-reference CASPT2 theory and molecular dynamics simulated in the presence of the reactants' kinetic and internal energies, the C5-terminal O₂ -addition was pinpointed as the most probable reaction pathway. By elucidating the reaction mechanism, kinetics and dynamics, and reaction products and energetics, this work constitutes the first report unraveling the synergetic damage of OG by ionizing radiation and singlet O₂ .

Molecular properties affecting the hydration of acid-base clusters

In the atmosphere, water in all phases is ubiquitous and plays important roles in catalyzing atmospheric chemical reactions, participating in cluster formation and affecting the composition of aerosol particles. Direct measurements of water-containing clusters are limited because water is likely to evaporate before detection, and therefore, theoretical tools are needed to study hydration in the atmosphere. We have studied thermodynamics and population dynamics of the hydration of different atmospherically relevant base monomers as well as sulfuric acid-base pairs. The hydration ability of a base seems to follow in the order of gas-phase base strength whereas hydration ability of acid-base pairs, and thus clusters, is related to the number of hydrogen binding sites. Proton transfer reactions at water-air interfaces are important in many environmental and biological systems, but a deeper understanding of their mechanisms remain elusive. By studying thermodynamics of proton transfer reactions in clusters containing up to 20 water molecules and a base molecule, we found that that the ability of a base to accept a proton in a water cluster is related to the aqueous-phase basicity. We also studied the second deprotonation reaction of a sulfuric acid in hydrated acid-base clusters and found that sulfate formation is most favorable in the presence of dimethylamine. Molecular properties related to the proton transfer ability in water clusters are discussed.

Fragment-Based Local Coupled Cluster Embedding Approach for the Quantification and Analysis of Noncovalent Interactions: Exploring the Many-Body Expansion of the Local Coupled Cluster Energy

Herein, we introduce a fragment-based local coupled cluster embedding approach for the accurate quantification and analysis of noncovalent interactions in molecular aggregates. Our scheme combines two different expansions of the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) energy: the many-body expansion (MBE) and the local energy decomposition (LED). The low-order terms in the MBE are initially computed in the presence of an environment that is treated at a low level of theory. Then, LED is used to decompose the energy of each term in the embedded MBE into additive fragment and fragment-pairwise contributions. This information is used to quantify the total energy of the system while providing at the same time in-depth insights into the nature and cooperativity of noncovalent interactions. Two different approaches are introduced and tested, in which the environment is treated at different levels of theory: the local coupled cluster in the Hartree-Fock (LCC-in-HF) method, in which the environment is treated at the HF level; and the electrostatically embedded local coupled cluster method (LCC-in-EE), in which the environment is replaced by point charges. Both schemes are designed to preserve as much as possible the accuracy of the parent local coupled cluster method for total energies, while being embarrassingly parallel and less memory intensive. These schemes appear to be particularly promising for the study of large and complex molecular aggregates at the coupled cluster level, such as condensed phase systems and protein-ligand interactions.

Photoredox Catalytic Pentafluorosulfanylative Domino Cyclization of α-Substituted Alkenes to Oxaheterocycles by Using SF6

Virtually inert sulfur hexafluoride becomes a precious pentafluorosulfanylation agent, if properly activated by photoredox catalysis, to access α-fluoro and α-alkoxy SF5 -compounds. This advanced protocol converts SF6 in the presence of alkynols as bifunctional C-C- and C-O-bond forming reagents directly into pentafluorosulfanylated oxygen-containing heterocycles in a single step from α-substituted alkenes. The proposed mechanism is supported by theoretical calculations and gives insights not only in the pentafluorosulfanylation step but also into formation of the carbon-carbon bond and is in full agreement with Baldwin's cyclization rules. The key step is a radical type 5-, 6- respectively 7-exo-dig-cyclization. The synthesized oxaheterocycles cannot be simply prepared by other synthetic methods, show a high level of structural complexity and significantly expand the scope of pentafluorosulfanylated building blocks valuable for medicinal and material chemistry.

Efficient Computation of Geometries for Gold Complexes

Computationally obtaining structural parameters along a reaction coordinate is commonly performed with Kohn‐Sham density functional theory which generally provides a good balance between speed and accuracy. However, CPU times still range from inconvenient to prohibitive, depending on the size of the system under study. Herein, the tight binding GFN2‐xTB method [C. Bannwarth, S. Ehlert, S. Grimme, J. Chem. Theory Comput. 2019, 15, 1652] is investigated as an alternative to produce reasonable geometries along a reaction path, that is, reactant, product and transition state structures for a series of transformations involving gold complexes. A small mean error (1 kcal/mol) was found, with respect to an efficient composite hybrid‐GGA exchange‐correlation functional (PBEh‐3c) paired with a double‐ζ basis set, which is 2–3 orders of magnitude slower. The outlined protocol may serve as a rapid tool to probe the viability of proposed mechanistic pathways in the field of gold catalysis.

Electron Attachment to Cytosine: The Role of Water

We present an EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) study on the electron attachment process to solvated cytosine. The electron attachment in the bulk solvated cytosine occurs through a doorway mechanism, where the initial electron is localized on water. The electron is subsequently transferred to cytosine by the mixing of electronic and nuclear degrees of freedom, which occurs on an ultrafast time scale. The bulk water environment stabilizes the cytosine-bound anion by an extensive hydrogen-bond network and drastically enhances the electron transfer rate from that observed in the gas phase. Microhydration studies cannot reproduce the effect of the bulk water environment on the electron attachment process, and one needs to include a large number of water molecules in the calculation to obtain converged results. The predicted adiabatic electron affinity and electron transfer rate obtained from our QM/MM calculations are consistent with the available experimental results.

Strong Binding of Noble Gases to [B12X11]-: A Theoretical Study

We systematically explore the stability and properties of [B₁₂X₁₁NG]^- adducts resulting from the binding of noble gas atoms to anionic [B₁₂X₁₁]^- clusters in the gas phase of mass spectrometers. [B₁₂X₁₁]^- can be obtained by stripping one X^- off the icosahedral closo-dodecaborate dianion [B₁₂X₁₂]^2-. We study the binding of the noble gas atoms He, Ne, Ar, Kr, and Xe to [B₁₂X₁₁]^- with substituents X = F, Cl, Br, I, and CN. While He cannot be captured by these clusters and Ne only binds at low temperatures, the complexes with the heavier noble gas atoms Ar, Kr, and Xe show appreciable complexation energies and exceed 1 eV at room temperature in the case of [B₁₂(CN)₁₁Xe]^-. The predicted B-NG equilibrium distance in the complexes with Ar, Kr, and Xe is only 0.10-0.25 Å longer than the sum of the covalent radii of the two corresponding atoms, and a significant charge transfer from the noble gas atom to the icosahedral B₁₂ cage is observed.

DFT study of H2 adsorption at a Cu-SSZ-13 zeolite: a cluster approach

In this work the H2 adsorption at a Cu(i)-SSZ-13 exchanged zeolite was theoretically investigated. A systematic cluster approach was used and different density functionals (B3LYP, B3LYP-D3(BJ), M06L, PBE, PBE-D3(BJ) and ωB97XD) and a def2-SVP basis set were benchmarked. In order to select the best approach to the H2 adsorption over a Cu(i)-SSZ-13 cluster with 78 atoms (16 T-sites), two main tasks were performed: (1) a comparison between theoretical and experimental structures and (2) a comparison between theoretical and experimental adsorption enthalpies. By employing the most suitable functional - the ωB97X-D - the H2 interaction with the zeolite structure was studied by means of NBO, NCI, AIM and DLPNO-CCSD(T)/LED analyses.

Experimental and theoretical assessment of protonated Hoogsteen 9-methylguanine-1-methylcytosine base-pair dissociation: kinetics within a statistical reaction framework

We investigated the collision-induced dissociation (CID) reactions of a protonated Hoogsteen 9-methylguanine-1-methylcytosine base pair (HG-[9MG·1MC + H]+), which aims to address the mystery of the literature reported "anomaly" in product ion distributions and compare the kinetics of a Hoogsteen base pair with its Watson-Crick isomer WC-[9MG·1MC + H]+ (reported recently by Sun et al.; Phys. Chem. Chem. Phys., 2020, 22, 24986). Product ion cross sections and branching ratios were measured as a function of center-of-mass collision energy using guided-ion beam tandem mass spectrometry, from which base-pair dissociation energies were determined. Product structures and energetics were assessed using various theories, of which the composite DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97XD/6-311++G(d,p) was adopted as the best-performing method for constructing a reaction potential energy surface. The statistical Rice-Ramsperger-Kassel-Marcus theory was found to provide a useful framework for rationalizing the dominating abundance of [1MC + H]+ over [9MG + H]+ in the fragment ions of HG-[9MG·1MC + H]+. The kinetics analysis proved the necessity for incorporating into kinetics modeling not only the static properties of reaction minima and transition states but more importantly, the kinetics of individual base-pair conformers that have formed in collisional activation. The analysis also pinpointed the origin of the statistical kinetics of HG-[9MG·1MC + H]+vs. the non-statistical behavior of WC-[9MG·1MC + H]+ in terms of their distinctively different intra-base-pair hydrogen-bonds and consequently the absence of proton transfer between the N1 position of 9MG and the N3' of 1MC in the Hoogsteen base pair. Finally, the Hoogsteen base pair was examined in the presence of a water ligand, i.e., HG-[9MG·1MC + H]+·H2O. Besides the same type of base-pair dissociation as detected in dry HG-[9MG·1MC + H]+, secondary methanol elimination was observed via the SN2 reaction of water with nucleobase methyl groups.

Calculating Reliable Gibbs Free Energies for Formation of Gas-Phase Clusters that Are Critical for Atmospheric Chemistry: (H2SO4)3

The effects of atmospheric aerosols on our climate are one of the biggest uncertainties in global climate models. Calculating the pathway for the formation of pre-nucleation clusters that become aerosols is challenging, requiring a comprehensive analysis of configurational space and highly accurate Gibbs free energy calculations. We identified a large set of minimum energy configurations of (H₂SO₄)₃ using a sampling technique based on a genetic algorithm and a stepwise density functional theory (DFT) approach and computed the thermodynamics of formation of these configurations with more accurate wavefunction-based electronic energies computed on the DFT geometries. The DLPNO-CCSD(T) methods always return more positive energies compared to the DFT energies. Within the DLPNO-CCSD(T) methods, extrapolating to the complete basis set limit gives more positive free energies compared to explicitly correlated single-point energies. The CBS extrapolation was shown to be robust as both the 4-5 inverse polynomial and Riemann zeta function schemes were within chemical accuracy of one another.

CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions?

We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.