Improved Correlation Energy Extrapolation Schemes Based on Local Pair Natural Orbital Methods

Cyclo-N9 -: a Novel 5/6 Fused Polynitrogen Anion for High Energy Density Materials

After cyclo-pentazolate anion, a 5/6 fused structure of N9 - is constructed, and four novel nitrogen-rich ionic compounds are assembled on its basis. The results of the quantum calculations revealed an uneven distribution of electrons on cyclo-N9 -, with significant charge density near the N5/N9 atoms and an ADCH charge of -0.425. The relative strength of chemical bonds is assessed through bond order analysis, which is further supplemented by transition state theory and ab initio molecular dynamics, ultimately leading to the identification of the decomposition pathways of cyclo-N9 -. The aromaticity of cyclo-N9 - and its individual sub-rings is cleverly validated through a combination of NICS_ZZ and ICSS methods. Among the eight systems, cyclo-N9 - forms hydrogen bonds with other cations, and IGMH analysis revealed significant LP-π and π-π stacking interactions between [N7H4 +] and cyclo-N9 -, both of which enhance system stability. The theoretical energy densities in all systems are at the forefront in the currently emerging nitrogen-rich compounds. Attributed to its extraordinarily high enthalpy of formation, the detonation performance of [N7H4 +] [N9 -] is particularly excellent. However, [NH3OH+] [N9 -] exhibits better stability and most exciting performance, making it a highly promising candidate with application potential.

Impact of Intracellular Proteins on μ-Opioid Receptor Structure and Ligand Binding

Chronic pain is a prevalent problem affecting approximately one out of every five adults in the U.S. The most effective way to treat chronic pain is with opioids, but they cause dangerous side effects such as tolerance, addiction, and respiratory depression, which makes them quite deadly. Opioids, such as fentanyl, target the μ-opioid receptor (MOR), which can then bind to the intracellular Gi protein or the β-arrestin protein. The Gi pathway is primarily responsible for pain relief and potential side effects, but the β-arrestin pathway is chiefly responsible for the unwanted side effects. Ideally, an effective pain medication without side effects would bind to MOR, which would bias signaling solely through the Gi pathway. We used the Bio3D library to conduct principal component analysis to compare the cryo-electron microscopy MOR structure in complex with the Gi versus an X-ray crystallography MOR structure with a nanobody acting as a Gi mimic. Our results agree with a previous study by Munro, which concluded that nanobody-bound MOR is structurally different than Gi-bound MOR. Furthermore, we investigated the structural diversity of opioids that can bind to MOR. Quantum mechanical calculations show that the low energy solution structures of fentanyl differ from the one bound to MOR in the experimental structure, and pKa calculations reveal that fentanyl is protonated in aqueous solution. Glide docking studies show that higher energy structures of fentanyl in solution form favorable docking complexes with MOR. Our calculations show the relative abundance of each fentanyl conformation in solution as well as the energetic barriers that need to be overcome to bind to MOR. Docking studies confirm that multiple fentanyl conformations can bind to the receptor. Perhaps a variety of conformations of fentanyl can stabilize multiple conformations of the MOR, which can explain why fentanyl can induce different intracellular signaling and multiple physiological effects.

Au-Catalyzed 5C Reaction of Type II Diene-Ynenes toward Dihydrosemibullvalenes: Reaction Development and Mechanistic Study

We report an unexpected gold-catalyzed 5C reaction of type II diene-ynenes to synthesize dihydrosemibullvalenes, which are potential bioisosteres for drug discovery. This 5C reaction occurs through a sequence of elementary reactions of cyclopropanation/Cope rearrangement/carbon shift/cyclopropanation/C-H insertion (shortened here as the 5C reaction), supported by quantum chemistry calculations. Mechanistic studies have also been applied to answer why type-II diene-ynenes cannot access seven-membered carbocycles-embedded bridged molecules under the gold catalysis, finding that the chair-like Cope rearrangement transition state (not the traditional boat-like transition state) is the key to the change of regiochemistry.

Electrophilic Intermediates in the Nef and Meyer Reactions: A Computational Study

The generation, interconversion, and reactivity of electrophilic species generated upon activation of nitroalkanes with protic acids (the Nef and Meyer reactions) were investigated by quantum-chemical calculations. N,N-Bis(hydroxy)iminium (R2C═N+(OH)2) and N-oxoiminium (R2C═N+═O) cations were shown to be produced independently from aci-nitroalkanes, while N-hydroxynitrilium cations (RC≡N+-OH) were formed via nearly barrierless C-H bond cleavage in N-oxoiminium cations. The N-oxoiminium and N-hydroxynitrilium cations whose formation is favored under highly acidic anhydrous conditions are strong electrophiles capable of reacting even with nonactivated arenes under ambient conditions. The N-oxoiminium cations R2C═N+═O are highly unusual ambident species containing three contiguous electrophilic centers (C, N, and O atoms). Nucleophilic addition at the oxygen atom is less preferred than the C- and N-attack yet possible in an intramolecular variant. These computational results shed light on some key aspects of the mechanisms of the Nef and Meyer reactions and predict the possibility of numerous interrupted versions of these reactions.

Assessment of the applicability of DFT methods to [Cp*Rh]-catalyzed hydrogen evolution processes

The present computational study provides a benchmark of density functional theory (DFT) methods in describing hydrogen evolution processes catalyzed by [Cp*Rh]-containing organometallic complexes. A test set was composed of 26 elementary reactions featuring chemical transformations and bonding situations essential for the field, including the emerging concept of non-innocent Cp* behavior. Reference values were obtained from a highly accurate 3/4 complete basis set and 6/7 complete PNO space extrapolated DLPNO-CCSD(T) energies. The performance of lower-level extrapolation procedures was also assessed. We considered 84 density functionals (DF) (including 13 generalized gradient approximations (GGA), nine meta-GGAs, 33 hybrids, and 29 double-hybrids) and three composite methods (HF-3c, PBEh-3c, and r2SCAN-3c), combined with different types of dispersion corrections (D3(0), D3BJ, D4, and VV10). The most accurate approach is the PBE0-DH-D3BJ (MAD of 1.36 kcal mol-1) followed by TPSS0-D3BJ (MAD of 1.60 kcal mol-1). Low-cost r2SCAN-3c composite provides a less accurate but much faster alternative (MAD of 2.39 kcal mol-1). The widely used Minnesota-family M06-L, M06, and M06-2X DFs should be avoided (MADs of 3.70, 3.94, and 4.01 kcal mol-1, respectively).

State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms

In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.

Experimental detection of the diamino-pentazolium cation and theoretical exploration of derived high energy materials

In this work, we realized the detection of diamino-pentazolium cation (DAPZ+) in the reaction solution experimentally and proved it to be meta-diamino-pentazole based on the transition state theory. Quantum chemical methods were used to predict its spectral properties, charge distribution, stability and aromaticity. Considering that DAPZ+ has excellent detonation properties, it was further explored by assembling it with N5-, N3- and C(NO2)3- anions, respectively. The results show a strong interaction between DAPZ+ and the three anions, which will have a positive effect on its stability. Thanks to the high enthalpy of formation and density, the calculated detonation properties of the three systems are exciting, especially [DAPZ+][N5-] (D: 10,016 m·s-1; P: 37.94 GPa), whose actual detonation velocity may very likely exceed CL-20 (D: 9773 m·s-1). There is no doubt that this work will become the precursor for the theoretical exploration of new polynitrogen ionic compounds.

Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes

Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.

Reduced Scaling Real-Time Coupled Cluster Theory

Real-time coupled cluster (CC) methods have several advantages over their frequency-domain counterparts, namely, response and equation of motion CC theories. Broadband spectra, strong fields, and pulse manipulation allow for the simulation of complex spectroscopies that are unreachable using frequency-domain approaches. Due to the high-order polynomial scaling, the required numerical time propagation of the CC residual expressions is a computationally demanding process. This scaling may be reduced by local correlation schemes, which aim to reduce the size of the (virtual) orbital space by truncation according to user-defined parameters. We present the first application of local correlation to real-time CC. As in previous studies of locally correlated frequency-domain CC, traditional local correlation schemes are of limited utility for field-dependent properties; however, a perturbation-aware scheme proves promising. A detailed analysis of the amplitude dynamics suggests that the main challenge is a strong time dependence of the wave function sparsity.

Further investigations into a Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals: Multipole correction, OSV extrapolation, and critical assessment

We report further investigations to aid the development of a Laplace MP2 (second-order Møller Plesset) method with a range separated Coulomb potential partitioned into short- and long-range parts. The implementation of the method extensively uses sparse matrix algebra, density fitting techniques for the short-range part, and a Fourier transformation in spherical coordinates for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas virtual space is described by orbital specific virtual orbitals (OSVs) associated with localized molecular orbitals. The Fourier transform is deficient for very large distances between localized occupied orbitals, and a multipole expansion for widely separated pairs is introduced for the direct MP2 contribution, which is applicable also to non-Coulombic potentials that do not satisfy the Laplace equation. For the exchange contribution, an efficient screening of contributing localized occupied pairs is employed, which is discussed more completely here. To mitigate errors due to the truncation of OSVs, a simple and efficient extrapolation procedure is used to obtain results close to MP2 for the full basis set of atomic orbitals Using a suitable set of default parameters, the accuracy of the approach is demonstrated. The current implementation of the approach is not very efficient, and the aim of this paper is to introduce and critically discuss ideas that can have more general applicability beyond MP2 calculations for large molecules.

Molecular Dynamics and Structural Studies of Zinc Chloroquine Complexes

Chloroquine (CQ) is a first-choice drug against malaria and autoimmune diseases. It has been co-administered with zinc against SARS-CoV-2 and soon dismissed because of safety issues. The structural features of Zn-CQ complexes and the effect of CQ on zinc distribution in cells are poorly known. In this study, state-of-the-art computations combined with experiments were leveraged to solve the structural determinants of zinc-CQ interactions in solution and the solid state. NMR, ESI-MS, and X-ray absorption and diffraction methods were combined with ab initio molecular dynamics calculations to address the kinetic lability of this complex. Within the physiological pH range, CQ binds Zn²⁺ through the quinoline ring nitrogen, forming [Zn(CQH)Cl_x(H₂O)_3-x]^(3+)-x (x = 0, 1, 2, and 3) tetrahedral complexes. The Zn(CQH)Cl₃ species is stable at neutral pH and at high chloride concentrations typical of the extracellular medium, but metal coordination is lost at a moderately low pH as in the lysosomal lumen. The pentacoordinate complex [Zn(CQH)(H₂O)₄]³⁺ may exist in the absence of chloride. This in vitro/in silico approach can be extended to other metal-targeting drugs and bioinorganic systems.

The driving effects of common atmospheric molecules for formation of prenucleation clusters: the case of sulfuric acid, formic acid, nitric acid, ammonia, and dimethyl amine

How secondary aerosols form is critical as aerosols' impact on Earth's climate is one of the main sources of uncertainty for understanding global warming. The beginning stages for formation of prenucleation complexes, that lead to larger aerosols, are difficult to decipher experimentally. We present a computational chemistry study of the interactions between three different acid molecules and two different bases. By combining a comprehensive search routine covering many thousands of configurations at the semiempirical level with high level quantum chemical calculations of approximately 1000 clusters for every possible combination of clusters containing a sulfuric acid molecule, a formic acid molecule, a nitric acid molecule, an ammonia molecule, a dimethylamine molecule, and 0-5 water molecules, we have completed an exhaustive search of the DLPNO-CCSD(T)/CBS//ωB97X-D/6-31++G** Gibbs free energy surface for this system. We find that the detailed geometries of each minimum free energy cluster are often more important than traditional acid or base strength. Addition of a water molecule to a dry cluster can enhance stabilization, and we find that the (SA)(NA)(A)(DMA)(W) cluster has special stability. Equilibrium calculations of SA, FA, NA, A, DMA, and water using our quantum chemical ΔG° values for cluster formation and realistic estimates of the concentrations of these monomers in the atmosphere reveals that nitric acid can drive early stages of particle formation just as efficiently as sulfuric acid. Our results lead us to believe that particle formation in the atmosphere results from the combination of many different molecules that are able to form highly stable complexes with acid molecules such as SA, NA, and FA.

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.

Comparative Quantum Chemistry Study on the Unimolecular Decomposition Channels of Pyrazole and Imidazole Energetic Materials

The difference in the initial decomposition step of pyrazoles and imidazoles was explored using the M062X method for optimization and G4-MP2 and approximated CCSD(T) methods for energies. Laplacian bond order analysis was used to study the effect of the nitro group on the bond strength and predict the bond dissociation energy (BDE) of the ring. Thermochemistry results show that the most possible decay channel of 1H-pyrazole and 3-nitropyrazole is the N₂ elimination, while the preferred initial step of 1H-imidazole is the CHN elimination. However, the nitro-nitrite isomerization dominates the decomposition of other nitro derivatives of 1H-pyrazole and 1H-imidazole. As for the formation of HO and HONO, the high energy barrier makes it difficult to take place. Based on the analysis of the lowest energy barrier and the BDE of NO₂ loss, it can be concluded that imidazoles are more stable than pyrazoles. This work contributes to revealing the difference in the initial step of energetic isomers and the understanding of the decomposition mechanism of energetic azoles.

o-Semiquinone radical anion isolated as an amorphous porous solid

The use of metal cations is a commonly applied strategy to create S > 1/2 stable molecular systems containing semiquinone radicals. Persistent mono-semiquinonato complexes of diamagnetic metal ions (S = 1/2) have been hitherto less common and mostly limited to the complexes of heavy metal ions. In this work, a mono-semiquinonato complex of aluminum, derived from 1,2-dihydroxybenzene, is obtained using a surprisingly short and uncomplicated procedure. The isolated product is an amorphous and porous solid that exhibits very good stability under ambient conditions. To characterise its molecular and electronic structure, 9.7, 34 and 406 GHz EPR spectroscopy was used in concert with computational techniques (DFT and DLPNO-CCSD). It was revealed that the radical complex is composed of two chemically equivalent aluminum cations and two catechol-like ligands with the unpaired electron uniformly distributed between the two organic molecules. The good stability and porous structure make this complex applicable in heterogeneous aerobic reactions.

Revisiting stacking interactions in tetrathiafulvalene and selected derivatives using tight-binding quantum chemical calculations and local coupled-cluster method

The engineering of supramolecular architectures needs accurate descriptions of the intermolecular interactions in crystal structures. Tetrathiafulvalene (TTF) is an effective building block used in the construction of promising functional materials. The parallel packing of the neutral TTF-TTF system was studied previously using the high-level quantum chemical method, advancing it as a valuable model system. The recently developed tight-binding quantum chemical method GFN2-xTB and local coupled-cluster method DLPNO-CCSD(T) were used to investigate the stacking interactions of TTF and selected derivatives deposited in the Cambridge Structural Database. Using the interaction energy of the TTF-TTF dimer calculated at the CCSD(T)/CBS level as the reference, the accuracies of the two methods are investigated. The energy decomposition analysis within the DLPNO-CCSD(T) framework reveals the importance of dispersion interaction in the TTF-related stacking systems. The dispersion interaction density plot vividly shows the magnitude and distribution of the dispersion interaction, providing a revealing insight into the stacking interactions in crystal structures. The results show that the GFN2-xTB and DLPNO-CCSD(T) methods could achieve accuracy at an affordable computational cost, which would be valuable in understanding the nature of parallel stacking in supramolecular systems.

Mechanism for the Reaction of White Phosphorus with Cp2Cr2(CO)6 Leading Ultimately to the Triple-Decker Sandwich Cp2Cr2(μ-η5,η5-P5): A Theoretical Study

The experimentally known reaction of Cp₂Cr₂(CO)₆ with white phosphorus (P₄) to give CpCr(CO)₂(η³-P₃), Cp₂Cr₂(CO)₄(μ-η,²η²-P₂), and the triple-decker sandwich Cp₂Cr₂(μ-η,⁵η⁵-P₅) is of interest since the P₄ reactant having a tetrahedral cluster of four phosphorus atoms is converted to products having P₂, P₃, and P₅ ligands. The mechanism of this obviously complicated reaction can be dissected into three stages using a coupled cluster theoretical method that has been benchmarked with the P₂, Mn(CO)₅, and CpCr(CO)₃ dimerization processes. The first stage of the Cp₂Cr₂(CO)₆/P₄ reaction mechanism generates the unsaturated singlet intermediate Cp₂Cr₂(CO)₅ that combines with the P₄ reactant. Decarbonylation of the resulting Cp₂Cr₂(CO)₅(P₄) complex provides a singlet tetracarbonyl readily fragmenting into the stable triphosphacyclopropenyl complex CpCr(CO)₂(η³-P₃) and the chromium phosphide CpCr(CO)₂(P). The isomeric triplet tetracarbonyl Cp₂Cr₂(CO)₄(P₄), readily fragments into CpCr(CO)₂(η²-P₂), which can generate the stable diphosphaacetylene complex Cp₂Cr₂(CO)₄(η,²η²-P₂) as well as the pentamer [CpCr(CO)₂]₅(P₁₀). Combination of the coordinately unsaturated CpCr(CO)(η³-P₃) with CpCr(CO)₂(η²-P₂) can lead to a ring expansion. This generates the P₅ pentagonal ligand in a Cp₂Cr₂(CO)₃(P₅) precursor to the experimentally observed carbonyl-free triple-decker sandwich Cp₂Cr₂(μ-η,⁵η⁵-P₅) after three successive decarbonylations.

N-Heterocyclic Carbene Organocatalysis: With or Without Carbenes?

In this work the mechanism of the aldehyde umpolung reactions, catalyzed by azolium cations in the presence of bases, was studied through computational methods. Next to the mechanism established by Breslow in the 1950s that takes effect through the formation of a free carbene, we have suggested that these processes can follow a concerted asynchronous path, in which the azolium cation directly reacts with the substrate, avoiding the formation of the carbene intermediate. We hereby show that substituting the azolium cation, and varying the base or the substrate do not affect the preference for the concerted reaction mechanism. The concerted path was found to exhibit low barriers also for the reactions of thiamine with model substrates, showing that this path might have biological relevance. The dominance of the concerted mechanism can be explained through the specific structure of the key transition state, avoiding the liberation of the highly reactive, and thus unstable carbene lone pair, whereas activating the substrate through hydrogen-bonding interactions. Polar and hydrogen-bonding solvents, as well as the presence of the counterions of the azolium salts facilitate the reaction through carbenes, bringing the barriers of the two reaction mechanisms closer, in many cases making the concerted path less favorable. Thus, our data show that by choosing the exact components in a reaction, the mechanism can be switched to occur with or without carbenes.

Direct Kinetic Measurements and Master Equation Modelling of the Unimolecular Decomposition of Resonantly-Stabilized CH2CHCHC(O)OCH3 Radical and an Upper Limit Determination for CH2CHCHC(O)OCH3 + O2 Reaction

Methyl-Crotonate (MC, (E)-methylbut-2-enoate, CH3CHCHC(O)OCH3) is a potential component of surrogate fuels that aim to emulate the combustion of fatty acid methyl ester (FAME) biodiesels with significant unsaturated FAME content. MC has three allylic hydrogens that can be readily abstracted under autoignition and combustion conditions to form a resonantly-stabilized CH2CHCHC(O)OCH3 radical. In this study we have utilized photoionization mass spectrometry to investigate the O2 addition kinetics and thermal unimolecular decomposition of CH2CHCHC(O)OCH3 radical. First we determined an upper limit for the bimolecular rate coefficient of CH2CHCHC(O)OCH3 + O2 reaction at 600 K (k ≤ 7.5 × 10−17 cm3 molecule−1 s−1). Such a small rate coefficient suggest this reaction is unlikely to be important under combustion conditions and subsequent efforts were directed towards measuring thermal unimolecular decomposition kinetics of CH2CHCHC(O)OCH3 radical. These measurements were performed between 750 and 869 K temperatures at low pressures (<9 Torr) using both helium and nitrogen bath gases. The potential energy surface of the unimolecular decomposition reaction was probed at density functional (MN15/cc-pVTZ) level of theory and the electronic energies of the stationary points obtained were then refined using the DLPNO-CCSD(T) method with the cc-pVTZ and cc-pVQZ basis sets. Master equation simulations were subsequently carried out using MESMER code along the kinetically important reaction pathway. The master equation model was first optimized by fitting the zero-point energy corrected reaction barriers and the collisional energy transfer parameters Δ E down , ref $\Delta{E_{{\text{down}},\;{\text{ref}}}}$ and n to the measured rate coefficients data and then utilize the constrained model to extrapolate the decomposition kinetics to higher pressures and temperatures. Both the experimental results and the MESMER simulations show that the current experiments for the thermal unimolecular decomposition of CH2CHCHC(O)OCH3 radical are in the fall-off region. The experiments did not provide definite evidence about the primary decomposition products.

An Extensive Set of Accurate Fluoride Ion Affinities for p-Block Element Lewis Acids and Basic Design Principles for Strong Fluoride Ion Acceptors

The computed fluoride ion affinity (FIA) is a valuable descriptor to assess the Lewis acidity of a compound. Despite its widespread use, the varying accuracy of applied computational models hampers the broad comparability of literature data. Herein, we evaluate the performance of selected methods (like DLPNO‐CCSD(T)) in FIA computations against CCSD(T)/CBS data and guide for the choice of suitable density functionals that allow the treatment of larger Lewis acids. Based on the benchmarked methods, we computed an extensive set of gas‐phase and solvation corrected FIA, that is covering group 13–16 elements featuring moderate to strong electron‐withdrawing substituents (190 entries). It permits an unbiased comparison of FIA over a significant fraction of the periodic table, serves as a source of reference for future synthetic or theoretical studies, and allows to derive some simple design principles for strong fluoride ion acceptors. Finally, the manuscript includes a tutorial section for the computation of FIA with and without the consideration of solvation.

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pK_a values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

Structures and internal dynamics of diphenylether and its aggregates with water

We report on a detailed multi-spectroscopic analysis of the structures and internal dynamics of diphenylether and its aggregates with up to three water molecules by employing molecular beam experiments.

Accurate Band Gap Predictions of Semiconductors in the Framework of the Similarity Transformed Equation of Motion Coupled Cluster Theory

In this work, we present a detailed comparison between wave-function-based and particle/hole techniques for the prediction of band gap energies of semiconductors. We focus on the comparison of the back-transformed Pair Natural Orbital Similarity Transformed Equation of Motion Coupled-Cluster (bt-PNO-STEOM-CCSD) method with Time Dependent Density Functional Theory (TD-DFT) and Delta Self Consistent Field/DFT (Δ-SCF/DFT) that are employed to calculate the band gap energies in a test set of organic and inorganic semiconductors. Throughout, we have used cluster models for the calculations that were calibrated by comparing the results of the cluster calculations to periodic DFT calculations with the same functional. These calibrations were run with cluster models of increasing size until the results agreed closely with the periodic calculation. It is demonstrated that bt-PNO-STEOM-CC yields accurate results that are in better than 0.2 eV agreement with the experiment. This holds for both organic and inorganic semiconductors. The efficiency of the employed computational protocols is thoroughly discussed. Overall, we believe that this study is an important contribution that can aid future developments and applications of excited state coupled cluster methods in the field of solid-state chemistry and heterogeneous catalysis.

In this work, it is shown that a combination of the embedded cluster approach with wave-function-based ab initio methods in the framework of the Similarity Transformed Equation of Motion Coupled Cluster (bt-PNO STEOM-CC) provides an accurate protocol for band gap energy predictions in classes of organic and inorganic semiconductors.

Diagrammatic Coupled Cluster Monte Carlo

We propose a modified coupled cluster Monte Carlo algorithm that stochastically samples connected terms within the truncated Baker-Campbell-Hausdorff expansion of the similarity-transformed Hamiltonian by construction of coupled cluster diagrams on the fly. Our new approach-diagCCMC-allows propagation to be performed using only the connected components of the similarity-transformed Hamiltonian, greatly reducing the memory cost associated with the stochastic solution of the coupled cluster equations. We show that for perfectly local, noninteracting systems diagCCMC is able to represent the coupled cluster wavefunction with a memory cost that scales linearly with system size. The favorable memory cost is observed with the only assumption of fixed stochastic granularity and is valid for arbitrary levels of coupled cluster theory. Significant reduction in memory cost is also shown to smoothly appear with dissociation of a finite chain of helium atoms. This approach is also shown not to break down in the presence of strong correlation through the example of a stretched nitrogen molecule. Our novel methodology moves the theoretical basis of coupled cluster Monte Carlo closer to deterministic approaches.

Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods

We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).

Multistep Explicit Solvation Protocol for Calculation of Redox Potentials

The calculation of molecular redox potentials in aqueous solution presents a challenge to quantum chemistry due to the need to calculate charged, open-shell species experiencing large solvent effects. Traditionally, redox potentials are calculated via the use of density functional theory and continuum solvation methods, but such protocols have been found to often suffer from large errors, particularly in the case of aqueous solution. While explicit solvation models hold promise of higher accuracy to describe solvent effects in general, their complicated use and lack of well-defined, reliable protocols has hindered their adoption. In this study, we present an explicit-solvation-based approach for the calculation of molecular redox potentials. We combine the use of affordable semiempirical QM/MM molecular dynamics (making use of the recently proposed GFN-xTB method by Grimme et al.) for both redox states and use the linear response approximation to relate vertical ionization energies to the adiabatic redox potential. Simulation length, averaging over snapshots, and accounting for bulk and polarization effects are systematically evaluated using phenol as a working example. We find that it is crucial to reliably account for bulk solvation effects in these calculations, as well as polarization effects which we divide up into short-range and long-range contributions. The short-range polarization contribution is accounted for via QM-region expansion, while the long-range contribution is accounted for via Drude-polarizable QM/MM. Our multistep protocol has been coded to be used in a fully automatic way in a local version of Chemshell. It has been evaluated on a test set of oxidation potentials of organic molecules and found to give gas-solution redox shifts with a mean absolute error of 0.13 eV with respect to experiment, compared to mean absolute errors of 0.26 and 0.21 eV with CPCM and SMD continuum models, respectively.

Performance of density functional theory for describing hetero-metallic active-site motifs for methane-to-methanol conversion in metal-exchanged zeolites

Methane-to-methanol conversion (MMC) can be facilitated with high methanol selectivities by copper-exchanged zeolites. There are however two open questions regarding the use of these zeolites to facilitate the MMC process. The first concerns the possibility of operating the three cycles in the stepwise MMC process by these zeolites in an isothermal fashion. The second concerns the possibility of improving the methanol yields by systematic substitution of some copper centers in these active sites with other earth-abundant transition metals. Quantum-mechanical computations can be used to compare methane activation by copper oxide species and analogous mixed-metal systems. To carry out such screening, it is important that we use theoretical methods that are accurate and computationally affordable for describing the properties of the hetero-metallic catalytic species. We have examined the performance of 47 exchange-correlation density functionals for predicting the relative spin-state energies and chemical reactivities of six hetero-metallic [M-O-Cu]²⁺ and [M-O₂ -Cu]²⁺ , (where MCo, Fe, and Ni), species by comparison with coupled cluster theory including iterative single, double excitations as well as perturbative treatment of triple excitations, CCSD(T). We also performed multireference calculations on some of these systems. We considered two types of reactions (hydrogen addition and oxygen addition) that are relevant to MMC. We recommend the use of τ-HCTH and OLYP to determine the spin-state energy splittings in the hetero-metallic motifs. ωB97, ωB97X, ωB97X-D3, and MN15 performed best for predicting the energies of the hydrogen and oxygen addition reactions. In contrast, local, and semilocal functionals do poorly for chemical reactivity. Using [Fe-O-Cu]²⁺ as a test, we see that the nonlocal functionals perform well for the methane CH activation barrier. In contrast, the semilocal functionals perform rather poorly. © 2018 Wiley Periodicals, Inc.

Systematic High-Accuracy Prediction of Electron Affinities for Biological Quinones

Quinones play vital roles as electron carriers in fundamental biological processes; therefore, the ability to accurately predict their electron affinities is crucial for understanding their properties and function. The increasing availability of cost-effective implementations of correlated wave function methods for both closed-shell and open-shell systems offers an alternative to density functional theory approaches that have traditionally dominated the field despite their shortcomings. Here, we define a benchmark set of quinones with experimentally available electron affinities and evaluate a range of electronic structure methods, setting a target accuracy of 0.1 eV. Among wave function methods, we test various implementations of coupled cluster (CC) theory, including local pair natural orbital (LPNO) approaches to canonical and parameterized CCSD, the domain-based DLPNO approximation, and the equations-of-motion approach for electron affinities, EA-EOM-CCSD. In addition, several variants of canonical, spin-component-scaled, orbital-optimized, and explicitly correlated (F12) Møller-Plesset perturbation theory are benchmarked. Achieving systematically the target level of accuracy is challenging and a composite scheme that combines canonical CCSD(T) with large basis set LPNO-based extrapolation of correlation energy proves to be the most accurate approach. Methods that offer comparable performance are the parameterized LPNO-pCCSD, the DLPNO-CCSD(T₀ ), and the orbital optimized OO-SCS-MP2. Among DFT methods, viable practical alternatives are only the M06 and the double hybrids, but the latter should be employed with caution because of significant basis set sensitivity. A highly accurate yet cost-effective DLPNO-based coupled cluster approach is used to investigate the methoxy conformation effect on the electron affinities of ubiquinones found in photosynthetic bacterial reaction centers. © 2018 Wiley Periodicals, Inc.

Combining Accuracy and Efficiency: An Incremental Focal-Point Method Based on Pair Natural Orbitals

In this work, we present a new pair natural orbitals (PNO)-based incremental scheme to calculate CCSD(T) and CCSD(T0) reaction, interaction, and binding energies. We perform an extensive analysis, which shows small incremental errors similar to previous non-PNO calculations. Furthermore, slight PNO errors are obtained by using T_PNO = T_TNO with appropriate values of 10^-7 to 10^-8 for reactions and 10^-8 for interaction or binding energies. The combination with the efficient MP2 focal-point approach yields chemical accuracy relative to the complete basis-set (CBS) limit. In this method, small basis sets (cc-pVDZ, def2-TZVP) for the CCSD(T) part are sufficient in case of reactions or interactions, while some larger ones (e.g., (aug)-cc-pVTZ) are necessary for molecular clusters. For these larger basis sets, we show the very high efficiency of our scheme. We obtain not only tremendous decreases of the wall times (i.e., factors >10²) due to the parallelization of the increment calculations as well as of the total times due to the application of PNOs (i.e., compared to the normal incremental scheme) but also smaller total times with respect to the standard PNO method. That way, our new method features a perfect applicability by combining an excellent accuracy with a very high efficiency as well as the accessibility to larger systems due to the separation of the full computation into several small increments.

Basic Phosphonium Ionic Liquids as Wittig Reagents

The possibility of designing a solvent/reagent for Wittig reactions from basic phosphonium salts is explored theoretically. In the suggested R4P+PhO- and Ph3PR+PhO- ionic liquids (ILs), the phenolate anion is prone to remove the α-proton from the alkyl chains, forming a phosphorous ylide. Significant hydrogen bonding between the oxygen atoms of the anions and α-hydrogen atoms of the cations were found by molecular dynamics simulations of these substances; therefore, proton transfer between the two ions is inherently supported by the structure of the liquid as well. The subsequent steps of the Wittig reaction from the phosphorous ylide were also found to be energetically possible. The mesoscopic structure of these materials exhibits a significant segregation into polar and nonpolar domains, which may also allow an easy dissolution of the substrates. The formation of a pentacoordinated phosphorous derivative through P-O bond formation was found to be also possible in the gas phase for both kind of compounds. Accordingly, having such basic anions in phosphonium-based ILs may produce such a neutral and therefore volatile species, which may hold further significant applications for these solvents in ion-exchange and separation techniques and in synthesis.

DLPNO-CCSD(T) scaled methods for the accurate treatment of large supramolecular complexes

In this work, we present scaled variants of the DLPNO-CCSD(T) method, dubbed as (LS)DLPNO-CCSD(T) and (NS)DLPNO-CCSD(T), to obtain accurate interaction energies in supramolecular complexes governed by noncovalent interactions. The novel scaled schemes are based on the linear combination of the DLPNO-CCSD(T) correlation energies calculated with the standard (LoosePNO and NormalPNO) and modified (Loose2PNO and Normal2PNO) DLPNO-CCSD(T) accuracy levels. The scaled DLPNO-CCSD(T) variants provide nearly TightPNO accuracy, which is essential for the quantification of weak noncovalent interactions, with a noticeable saving in computational cost. Importantly, the accuracy of the proposed schemes is preserved irrespective of the nature and strength of the supramolecular interaction. The (LS)DLPNO-CCSD(T) and (NS)DLPNO-CCSD(T) protocols have been used to study in depth the role of the CH-π versus π-π interactions in the supramolecular complex formed by the electron-donor truxene-tetrathiafulvalene (truxTTF) and the electron-acceptor hemifullerene (C₃₀ H₁₂ ). (NS)DLPNO-CCSD(T)/CBS calculations clearly reveal the higher stability of staggered (dominated by CH-π interactions) versus bowl-in-bowl (dominated by π-π interactions) arrangements in the truxTTF•C₃₀ H₁₂ heterodimer. Hemifullerene and similar carbon-based buckybowls are therefore expected to self-assemble with donor compounds in a richer way other than the typical concave-convex π-π arrangement found in fullerene-based aggregates. © 2017 Wiley Periodicals, Inc.