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

Uncertainty-Aware First-Principles Exploration of Chemical Reaction Networks

Exploring large chemical reaction networks with automated exploration approaches and accurate quantum chemical methods can require prohibitively large computational resources. Here, we present an automated exploration approach that focuses on the kinetically relevant part of the reaction network by interweaving (i) large-scale exploration of chemical reactions, (ii) identification of kinetically relevant parts of the reaction network through microkinetic modeling, (iii) quantification and propagation of uncertainties, and (iv) reaction network refinement. Such an uncertainty-aware exploration of kinetically relevant parts of a reaction network with automated accuracy improvement has not been demonstrated before in a fully quantum mechanical approach. Uncertainties are identified by local or global sensitivity analysis. The network is refined in a rolling fashion during the exploration. Moreover, the uncertainties are considered during kinetically steering of a rolling reaction network exploration. We demonstrate our approach for Eschenmoser-Claisen rearrangement reactions. The sensitivity analysis identifies that only a small number of reactions and compounds are essential for describing the kinetics reliably, resulting in efficient explorations without sacrificing accuracy and without requiring prior knowledge about the chemistry unfolding.

Encapsulation of charged halogens by the 512 water cage

This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.

Accurate Enthalpies of Formation for PFAS from First-Principles: Combining Different Levels of Theory in a Generalized Thermochemical Hierarchy

The enthalpies of formation are computed for a large number of per- and poly fluoroalkyl substances (PFAS) using a connectivity-based hierarchy (CBH) approach. A combination of different electronic structure methods are used to provide the reference data in a hierarchical manner. The ANL0 method, in conjunction with the active thermochemical tables, provides enthalpies of formation for smaller species with subchemical accuracy. Coupled-cluster theory with explicit correlations are used to compute enthalpies of formation for intermediate species, based upon the ANL0 results. For the largest PFAS, including perfluorooctanoic acid (PFOA) and heptafluoropropylene oxide dimer acid (GenX), coupled-cluster theory with local correlations is used. The sequence of homodesmotic reactions proposed by the CBH are determined automatically by a new open-source code, AutoCBH. The results are the first reported enthalpies of formation for the majority of the species. A convergence analysis and global uncertainty quantification confirm that the enthalpies of formation at 0 K should be accurate to within ±5 kJ/mol. This new approach is not limited to PFAS, but can be applied to many chemical systems.

Gas-Phase Oxidation of Atmospherically Relevant Unsaturated Hydrocarbons by Acyl Peroxy Radicals

This study assesses the atmospheric impact of reactions between unsaturated hydrocarbons such as isoprene and monoterpenes and peroxy radicals containing various functional groups. We find that reactions between alkenes and acyl peroxy radicals have reaction rates high enough to be feasible in the atmosphere and lead to high molar mass accretion products. Moreover, the reaction between unsaturated hydrocarbons and acyl peroxy radicals leads to an alkyl radical, to which molecular oxygen rapidly adds. This finding is confirmed by both theoretical calculations and experiments. The formed perester peroxy radical may either undergo further H-shift reactions or react bimolecularly. The multifunctional oxygenated compounds formed through acyl peroxy radical + alkene reactions are potentially important contributors to particle formation and growth. Thus, acyl peroxy radical-initiated oxidation chemistry may need to be included in atmospheric models.

Analytical Gradient Using Cluster-in-Molecule RI-MP2 Method for the Geometry Optimizations of Large Systems

We present an efficient analytical energy gradient algorithm for the cluster-in-molecule resolution-of-identity second-order Møller-Plesset perturbation (CIM-RI-MP2) method based on the Lagrange multiplier method. Our algorithm independently constructs the Lagrangian formalism within each cluster, avoiding the solution of the coupled-perturbed Hartree-Fock (CPHF) equation for the whole system. Due to this feature, the computational cost of the CIM-RI-MP2 gradients is much lower than that of other local MP2 algorithms. Benchmark calculations of several molecules containing up to 312 atoms demonstrate the general applicability of our CIM-RI-MP2 gradient algorithm. The optimized structure of a 244-atom molecule using the CIM-RI-MP2 method with the cc-pVDZ basis set is in good agreement with the corresponding crystal structure. A single-point gradient calculation conducted for a molecular cage containing 972 atoms and 9612 basis functions takes 48 h on 25 nodes, utilizing a total of 600 CPU cores. The present CIM-RI-MP2 gradient program is applicable for obtaining the optimized geometries of large systems with hundreds of atoms.

Backside versus Frontside SN2 Reactions of Alkyl Triflates and Alcohols

Nucleophilic substitution reactions are elementary reactions in organic chemistry that are used in many synthetic routes. By quantum chemical methods, we have investigated the intrinsic competition between the backside SN2 (SN2-b) and frontside SN2 (SN2-f) pathways using a set of simple alkyl triflates as the electrophile in combination with a systematic series of phenols and partially fluorinated ethanol nucleophiles. It is revealed how and why the well-established mechanistic preference for the SN2-b pathway slowly erodes and can even be overruled by the unusual SN2-f substitution mechanism going from strong to weak alcohol nucleophiles. Activation strain analyses disclose that the SN2-b pathway is favored for strong alcohol nucleophiles because of the well-known intrinsically more efficient approach to the electrophile resulting in a more stabilizing nucleophile-electrophile interaction. In contrast, the preference of weaker alcohol nucleophiles shifts to the SN2-f pathway, benefiting from a stabilizing hydrogen bond interaction between the incoming alcohol and the leaving group. This hydrogen bond interaction is strengthened by the increased acidity of the weaker alcohol nucleophiles, thereby steering the mechanistic preference toward the frontside SN2 pathway.

The impact of solvation on the structure and electric field strength in Li+GlyGly complexes

To scrutinise the impact of electric fields on the structure and vibrations of biomolecules in the presence of water, we study the sequential solvation of lithium diglycine up to three water molecules with cryogenic infrared action spectroscopy. Conformer-specific IR-IR spectroscopy and H2O/D2O isotopic substitution experiments provide most of the information required to decipher the structure of the observed conformers. Additional confirmation is provided by scaled harmonic vibrational frequency calculations using MP2 and DFT. The first water molecule is shown to bind to the Li+ ion, which weakens the electric field experienced by the peptide and as a consequence, also the strength of an internal NH⋯NH2 hydrogen bond in the diglycine backbone. The strength of this hydrogen bond decreases approximately linearly with the number of water molecules as a result of the decreasing electric field strength and coincides with an increase in the number of conformers observed in our spectra. The addition of two water molecules is already sufficient to change the preferred conformation of the peptide backbone, allowing for Li+ coordination to the lone pair of the terminal amine group.

Catalytic role of the enol ether intermediate in the intramolecular Stetter reaction: a computational perspective

The intramolecular Stetter reaction catalyzed by a carbene is investigated by density functional theory (DFT) calculations and kinetic simulations. Catalyst 1 first reacts with aldehyde 2 to give the primary adduct (PA). The PA undergoes the intramolecular oxa-Michael reaction to irreversibly generate enol ether intermediate 9. The conversion of the enol ether to the Breslow intermediate (BI) requires the assistance of a base such as the PA. The next step involves formation of a carbon-carbon bond through the Michael addition, and expulsion of the catalyst generates the Stetter product 7. Calculations show that the catalytic cycle is composed of two irreversible processes: the first one involves the exergonic formation of the enol ether intermediate, while the second one is the conversion of the enol ether to the final product. Kinetic simulations using initial concentrations of [1]0 = 0.005 M and [2]0 = 0.025 M demonstrate that under a steady-state condition, 35% of the catalyst rests on the state of the enol ether (0.0018 M). The catalyst resting state therefore consists of the unbound form (the free catalyst) and its bound form (the enol ether species). According to variable time normalization analysis, the reaction exhibits a second-order dependence (first order in catalyst and first order in substrate), which agrees with experiments. The oxa-Michael reaction to form the enol ether is identified to be turnover limiting in the intramolecular Stetter reaction, which rationalizes the observed electronic effect of the Michael acceptor on the reactivity, as well as the measured isotope effect with respect to the aldehydic proton/deuteron. The base that participates in the BI formation has a significant effect on the build-up of the resting state 9 and the active catalyst concentration. In addition, the thermodynamic stability of the enol ether is found to depend on the tether length between the aromatic aldehyde and the Michael acceptor, as well as the chemical nature of the carbene catalyst. The favorability for the oxa-Michael reaction is therefore suggested to govern the reactivity of the intramolecular Stetter transformation.

The mechanism and kinetics of the atmospheric oxidation of CF3(CF2)2CHCH2 (HFC-1447fz) by hydroxyl radicals: ab initio investigation

The oxidation of 3,3,4,4,5,5,5-heptafluoro-1-pentene (HFC-1447fz) by hydroxyl radicals plays a crucial role in atmospheric conditions. By employing the CCSD(T)/cc-pVTZ//M06-2X/6-311++G(d,p) level of theory, the detailed reaction mechanism, kinetics and atmospheric implications of the degradation of HFC-1447fz by hydroxyl radicals were investigated. Compared to H-abstraction channels, the OH addition reaction is determined to be more favorable initial pathways in the degradation processes of HFC-1447fz. The overall rate coefficient of the degradation of HFC-1447fz by OH radicals is estimated to be 1.66 × 10-12 cm3 molecule-1 s-1 and the lifetime of HFC-1447fz is found to be 7 days at 298 K, which are in good agreement with the reported experimental results. The global warming potential (GWP) for HFC-1447fz on the 50, 100 and 500-year time horizons is estimated using the calculated rate coefficient. Furthermore, the mechanisms of the subsequent reactions of two OH-addition adducts have also been investigated. By TD-DFT calculations, it was found that eleven species can undergo photodissociation, while ten other species are photolytically stable under sunlight.

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

Phosphaarsenes - Combining Phospha- and Arsa-Wittig-Reagents

Dipnictenes of the type RPn=PnR (Pn=P, As, Sb, Bi) can be viewed as dimers of the corresponding pnictinidenes R-Pn. Phosphanylidene- and arsanylidenephosphoranes (R-Pn(PMe3); Pn=P, As) have been shown to be versatile synthetic surrogates for the delivery of pnictinidene fragments. We now report that thermal treatment of 1 : 1 mixtures of R-P(PMe3) and R'-As(PMe3) gives access to arsaphosphenes of the type RP=AsR'. Three examples are presented and the properties and reactivity of Mes*P=AsDipTer (1) (Mes*=2,4,6-tBu3-C6H2; DipTer=2,6-(2,6-iPr2C6H3)2-C6H3) were studied in detail. Solid state 31P NMR spectroscopy revealed a large 31P NMR chemical shift anisotropy with a span of ca. 920 ppm for 1 while computational methods were employed to investigate this pronounced magnetic deshielding of the P atom in 1. In the presence of the carbene IMe4 (IMe4=:C(MeNCMe)2) 1 is shown to be split into the corresponding NHC adducts Mes*P(IMe4) and DipTerAs(IMe4), which is additionally shown for diarsenes.

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree-Fock exchange. These results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

Tensor Hypercontraction of Cluster Perturbation Theory: Quartic Scaling Perturbation Series for the Coupled Cluster Singles and Doubles Ground-State Energies

Even though cluster perturbation theory has been shown to be a robust noniterative alternative to coupled cluster theory, it is still plagued by high order polynomial computational scaling and the storage of higher order tensors. We present a proof-of-concept strategy for implementing a cluster perturbation theory ground-state energy series for the coupled cluster singles and doubles energy with N4 computational scaling using tensor hypercontraction (THC). The reduction in computational scaling by two orders is achieved by decomposing two electron repulsion integrals, doubles amplitudes and multipliers, as well as selected double intermediates to the THC format. Using the outlined strategy, we showcase that the THC pilot implementations retain numerical accuracy to within 1 kcal/mol relative to corresponding conventional and density fitting implementations, and we empirically verify the N4 scaling.

Role of curvature in acridone for 1 O2 oxidation of a natural product homoallylic alcohol: A novel iso-hydroperoxide intermediate

A density functional theoretical (DFT) study is presented, implicating a 1 O2 oxidation process to reach a dihydrobenzofuran from the reaction of the natural homoallylic alcohol, glycocitrine. Our results predict an interconversion between glycocitrine and an iso-hydroperoxide intermediate [R(H)O+ -O- ] that provides a key path in the chemistry which then follows. Formations of allylic hydroperoxides are unlikely from a 1 O2 'ene' reaction. Instead, the dihydrobenzofuran arises by 1 O2 oxidation facilitated by a 16° curvature of the glycocitrine ring imposed by a pyramidal N-methyl group. This curvature facilitates the formation of the iso-hydroperoxide, which is analogous to the iso species CH2 I+ -I- and CHI2 + -I- formed by UV photolysis of CH2 I2 and CHI3 . The iso-hydroperoxide is also structurally reminiscent of carbonyl oxides (R2 C=O+ -O- ) formed in the reaction of carbenes and oxygen. Our DFT results point to intermolecular process, in which the iso-hydroperoxide's fate relates to O-transfer and H2 O dehydration reactions for new insight into the biosynthesis of dihydrobenzofuran natural products.

Accidental triplet harvesting in donor-acceptor dyads with low spin-orbit coupling

We present an accidental mechanism for efficient intersystem crossing (ISC) between singlet and triplet states with low spin-orbit coupling (SOC) in molecules having donor-acceptor (D-A) moieties separated by a Sigma bond. Our study shows that SOC between the lowest singlet excited state and the higher-lying triplet states, together with nuclear motion-driven coupling of this triplet state with lower-lying triplet states during the free rotation about a Sigma bond, is one of the possible ways to achieve the experimentally observed ISC rate for a class of D-A type photoredox catalysts. This mechanism is found to be the dominant contributor to the ISC process with the corresponding rate reaching a maximum at a dihedral angle in the range of 72°-78° between the D-A moieties of 10-(naphthalen-1-yl)-3,7-diphenyl-10H-phenoxazine and other molecules included in the study. We have further demonstrated that the same mechanism is operative in a specific spirobis[anthracene]dione molecule, where the D and A moieties are interlocked near to the optimal dihedral angle, indicating the plausible effectiveness of the proposed mechanism. The present finding is expected to have implications in strategies for the synthesis of new generations of triplet-harvesting organic molecules.

Affinity of Telluronium Chalcogen Bond Donors for Lewis Bases in Solution: A Critical Experimental-Theoretical Joint Study

Telluronium salts [Ar2 MeTe]X were synthesized, and their Lewis acidic properties towards a number of Lewis bases were addressed in solution by physical and theoretical means. Structural X-ray diffraction analysis of 21 different salts revealed the electrophilicity of the Te centers in their interactions with anions. Telluroniums' propensity to form Lewis pairs was investigated with OPPh3 . Diffusion-ordered NMR spectroscopy suggested that telluroniums can bind up to three OPPh3 molecules. Isotherm titration calorimetry showed that the related heats of association in 1,2-dichloroethane depend on the electronic properties of the substituents of the aryl moiety and on the nature of the counterion. The enthalpies of first association of OPPh3 span -0.5 to -5 kcal mol-1 . Study of the affinity of telluroniums for OPPh3 by state-of-the-art DFT and ab-initio methods revealed the dominant Coulombic and dispersion interactions as well as an entropic effect favoring association in solution. Intermolecular orbital interactions between [Ar2 MeTe]+ cations and OPPh3 are deemed insufficient on their own to ensure the cohesion of [Ar2 MeTe ⋅ Bn ]+ complexes in solution (B=Lewis base). Comparison of Grimme's and Tkatchenko's DFT-D4/MBD-vdW thermodynamics of formation of higher [Ar2 MeTe ⋅ Bn ]+ complexes revealed significant molecular size-dependent divergence of the two methodologies, with MBD yielding better agreement with experiment.

Exploring the stability and aromaticity of rare earth doped tin cluster MSn16- (M = Sc, Y, La)

Rare earth elements have high chemical reactivity, and doping them into semiconductor clusters can induce novel physicochemical properties. The study of the physicochemical mechanisms of interactions between rare earth and tin atoms will enhance our understanding of rare earth functional materials from a microscopic perspective. Hence, the structure, electronic characteristics, stability, and aromaticity of endohedral cages MSn16- (M = Sc, Y, La) have been investigated using a combination of the hybrid PBE0 functional, stochastic kicking, and artificial bee colony global search technology. By comparing the simulated results with experimental photoelectron spectra, it is determined that the most stable structure of these clusters is the Frank-Kasper polyhedron. The doping of atoms has a minimal influence on density of states of the pure tin system, except for causing a widening of the energy gap. Various methods such as ab initio molecular dynamics simulations, the spherical jellium model, adaptive natural density partitioning, localized orbital locator, and electron density difference are employed to analyze the stability of these clusters. The aromaticity of the clusters is examined using iso-chemical shielding surfaces and the gauge-including magnetically induced currents. This study demonstrates that the stability and aromaticity of a tin cage can be systematically adjusted through doping.

Cost-effective approach for atmospheric accretion reactions: a case of peroxy radical addition to isoprene

We present an accurate and cost-effective method for investigating the accretion reactions between unsaturated hydrocarbons and oxidized organic radicals. We use accretion between isoprene and primary, secondary and tertiary alkyl peroxy radicals as model reactions. We show that a systematic semiempirical transition state search can lead to better transition state structures than relaxed scanning with density functional theory with a significant gain in computational efficiency. Additionally, we suggest accurate and effective quantum chemical methods to study accretion reactions between large unsaturated hydrocarbons and oxidized organic radicals. Furthermore, we examine the atmospheric relevance of these types of reactions by calculating the bimolecular reaction rate coefficients and formation rates under atmospheric conditions from the quantum chemical reaction energy barriers.

Taming nonclassical carbocations to control small ring reactivity

Prediction of the outcome of ring opening of small organic rings under cationic conditions can be challenging due to the intermediacy of nonclassical carbocations. For example, the solvolysis of cyclobutyl or cyclopropylmethyl derivatives generates up to four products on nucleophilic capture or elimination via cyclopropylcarbinyl and bicyclobutonium ions. Here, we show that such reaction outcomes can be controlled by subtle changes to the structure of nonclassical carbocation. Using bicyclo[1.1.0]butanes as cation precursors, the regio- and stereochemistry of ring opening is shown to depend on the degree and nature of the substituents on the cationic intermediates. Reaction outcomes are rationalized using computational models, resulting in a flowchart to predict product formation from a given cation precursor.

Exploring CPS-Extrapolated DLPNO-CCSD(T1) Reference Values for Benchmarking DFT Methods on Enzymatically Catalyzed Reactions

Domain-based local pair natural orbital coupled-cluster singles doubles with perturbative triples [DLPNO-CCSD(T)] is regularly used to calculate reliable benchmark reference values at a computational cost significantly lower than that of canonical CCSD(T). Recent work has shown that even greater accuracy can be obtained at only a small additional cost through extrapolation to the complete PNO space (CPS) limit. Herein, we test two levels of CPS extrapolation, CPS(5,6), which approximates the accuracy of standard TightPNO, and CPS(6,7), which surpasses it, as benchmark values to test density functional approximations (DFAs) on a small set of organic and transition-metal-dependent enzyme active site models. Between the different reference levels of theory, there are changes in the magnitudes of the absolute deviations for all functionals, but these are small and there is minimal impact on the relative rankings of the tested DFAs. The differences are more significant for the metalloenzymes than the organic enzymes, so we repeat the tests on our entire ENZYMES22 set of organic enzyme active site models [Wappett, D. A.; Goerigk, L. J. Phys. Chem. A 2019, 123, 7057-7074] to confirm that using the CPS extrapolations for the reference values has negligible impact on the benchmarking outcomes. This means that we can particularly recommend CPS(5,6) as an alternative to standard TightPNO settings for calculating reference values, increasing the applicability of DLPNO-CCSD(T) in benchmarking reaction energies and barrier heights of larger models of organic enzymes. DLPNO-CCSD(T1)/CPS(6,7) energies for ENZYMES22 are finally presented as updated reference values for the set, reflecting the recent improvements in the method.

Bimolecular versus Trimolecular Reaction Pathways for H2O2 with Hypochlorous Species and Implications for Wastewater Reclamation

The benchmark advanced oxidation technology (AOT) that uses UV/H2O2 integrated with hypochlorous species exhibits great potential in removing micropollutants and enhancing wastewater treatability for reclamation purposes. Although efforts have been made to study the reactions of H2O2 with hypochlorous species, there exist great discrepancies in the order of reaction kinetics, the rate constants, and the molecule-level mechanisms. This results in an excessive use of hypochlorous reagents and system underperformance during treatment processes. Herein, the titled reaction was investigated systematically through complementary experimental and theoretical approaches. Stopped-flow spectroscopic measurements revealed a combination of bi- and trimolecular reaction kinetics. The bimolecular pathway dominates at low H2O2 concentrations, while the trimolecular pathway dominates at high H2O2 concentrations. Both reactions were simulated using direct dynamics trajectories, and the pathways identified in the trajectories were further validated by high-level quantum chemistry calculations. The theoretical results not only supported the spectroscopic data but also elucidated the molecule-level mechanisms and helped to address the origin of the discrepancies. In addition, the impact of the environmental matrix was evaluated by using two waters with discrete characteristics, namely municipal wastewater and ammonium-rich wastewater. Municipal wastewater had a negligible matrix effect on the reaction kinetics of H2O2 and the hypochlorous species, making it a highly suitable candidate for this integration technique. The obtained in-depth reaction mechanistic insights will enable the development of a viable and economical technology for safe water reuse.

New Mechanistic Insights into the Primary Thermolysis Reactions of 1,3,4,6-Tetranitrooctahydroimidazo-[4,5-d]imidazole (BCHMX) from Predictive Local Coupled Cluster Calculations

Theoretical studies of the decomposition mechanism of energetic materials quite often scrutinize only the primary thermolysis reactions. However, the secondary reactions are crucial, inter alia, for proper building of the combustion models and understanding the autocatalytic processes. In the present study, we applied predictive DLPNO-CCSD(T) calculations to elucidate the kinetics and decomposition mechanism of a novel promising energetic material, 1,3,4,6-tetranitrooctahydroimidazo [4,5-d] imidazole (BCHMX). We identified eight previously unknown BCHMX conformers, both cis and trans in accordance to the spatial position of the H atoms bonded to a carbon bridge. Among them, the relative enthalpies of cis isomers lie in the narrow range ∼10 kJ mol-1 rendering them thermally accessible in the course of decomposition. The radical N-NO2 bond cleavage via one of the novel conformers is the dominant primary decomposition channel of BCHMX with the kinetic parameters Ea = 168.4 kJ mol-1 and log(A, s-1) = 18.5. We also resolved several contradictory assumptions on the mechanism and key intermediates of BCHMX thermolysis. To get a deeper understanding of the decomposition mechanism, we examined a series of unimolecular and bimolecular secondary channels of BCHMX. Among the former reactions, the C-C bond unzipping followed by another radical elimination of a nitro group is the most energetically favorable pathway with an activation barrier ∼113 kJ mol-1. However, contrary to the literature assumptions, the bimolecular H atom abstraction from a pristine BCHMX molecule by a primary nitramine radical product, not the nitro one, followed by another NO2 radical elimination, is the most important bimolecular secondary thermolysis reaction of BCHMX at lower temperatures. The isokinetic temperature of the bimolecular and unimolecular secondary reactions is ∼620 K. Unimolecular reactions might be important in dilute solutions, where bimolecular reactions are suppressed. The secondary reactions considered in the present work might be pertinent in the case of related energetic nitramines (e.g., RDX, HMX, and CL-20).

Gas-Phase Alkali-Metal Cation Affinities of Stabilized Enolates

The chemistry of alkali-metal enolates is dominated by ion pairing. To improve our understanding of the intrinsic interactions between the alkali-metal cations and the enolate anions, we have applied Cooks' kinetic method to determine relative M+ (M=Li, Na, K) affinities of the stabilized enolates derived from acetylacetone, ethyl acetoacetate, diethyl malonate, ethyl cyanoacetate, 2-cyanoacetamide, and methyl malonate monoamide in the gas phase. Quantum chemical calculations support the experimental results and moreover afford insight into the structures of the alkali-metal enolate complexes. The affinities decrease with increasing size of the alkali-metal cations, reflecting weaker electrostatic interactions and lower charge densities of the free M+ ions. For the different enolates, a comparison of their coordinating abilities is complicated by the fact that some of the free anions undergo conformational changes resulting in stabilizing intramolecular interactions. If these complicating effects are disregarded, the M+ affinities correlate with the electron density of the chelating functionalities, that is, the carbonyl and/or the nitrile groups of the enolates. A comparison with the known association constants of the corresponding alkali-metal enolates in solution points to the importance of solvation effects for these systems.

Shape resonance induced electron attachment to cytosine: The effect of aqueous media

We have investigated the impact of microsolvation on shape-type resonance states of nucleobases, taking cytosine as a case study. To characterize the resonance position and decay width of the metastable states, we employed the newly developed DLPNO-based EA-EOM-CCSD method in conjunction with the resonance via Padé (RVP) method. Our calculations show that the presence of water molecules causes a redshift in the resonance position and an increase in the lifetime for the three lowest-lying resonance states of cytosine. Furthermore, there are some indications that the lowest resonance state in isolated cytosine may get converted to a bound state in the presence of an aqueous environment. The obtained results are extremely sensitive to the basis set used for the calculations.

Effects of Intra-Base Pair Proton Transfer on Dissociation and Singlet Oxygenation of 9-Methyl-8-Oxoguanine-1-Methyl-Cytosine Base-Pair Radical Cations

8-Oxoguanosine is the most common oxidatively generated base damage and pairs with complementary cytidine within duplex DNA. The 8-oxoguanosine-cytidine lesion, if not recognized and removed, not only leads to G-to-T transversion mutations but renders the base pair being more vulnerable to the ionizing radiation and singlet oxygen (1 O2 ) damage. Herein, reaction dynamics of a prototype Watson-Crick base pair [9MOG ⋅ 1MC]⋅+ , consisting of 9-methyl-8-oxoguanine radical cation (9MOG⋅+ ) and 1-methylcystosine (1MC), was examined using mass spectrometry coupled with electrospray ionization. We first detected base-pair dissociation in collisions with the Xe gas, which provided insight into intra-base pair proton transfer of 9MOG⋅+  ⋅ 1MC← → ${{\stackrel{ {\rightarrow} } { {\leftarrow} } } }$ [9MOG - HN1 ]⋅ ⋅ [1MC+HN3' ]+ and subsequent non-statistical base-pair separation. We then measured the reaction of [9MOG ⋅ 1MC]⋅+ with 1 O2 , revealing the two most probable pathways, C5-O2 addition and HN7 -abstraction at 9MOG. Reactions were entangled with the two forms of 9MOG radicals and base-pair structures as well as multi-configurations between open-shell radicals and 1 O2 (that has a mixed singlet/triplet character). These were disentangled by utilizing approximately spin-projected density functional theory, coupled-cluster theory and multi-referential electronic structure modeling. The work delineated base-pair structural context effects and determined relative reactivity toward 1 O2 as [9MOG - H]⋅>9MOG⋅+ >[9MOG - HN1 ]⋅ ⋅ [1MC+HN3' ]+ ≥9MOG⋅+  ⋅ 1MC.

Improved Configurational Sampling Protocol for Large Atmospheric Molecular Clusters

The nucleation process leading to the formation of new atmospheric particles plays a crucial role in aerosol research. Quantum chemical (QC) calculations can be used to model the early stages of aerosol formation, where atmospheric vapor molecules interact and form stable molecular clusters. However, QC calculations heavily depend on the chosen computational method, and when dealing with large systems, striking a balance between accuracy and computational cost becomes essential. We benchmarked the binding energies and structures and found the B97-3c method to be a good compromise between the accuracy and computational cost for studying large cluster systems. Further, we carefully assessed configurational sampling procedures for targeting large atmospheric molecular clusters containing up to 30 molecules (approximately 2 nm in diameter) and proposed a funneling approach with highly improved accuracy. We find that several parallel ABCluster explorations lead to better guesses for the cluster global energy minimum structures than one long exploration. This methodology allows us to bridge computational studies of molecular clusters, which typically reach only around 1 nm, with experimental studies that often measure particles larger than 2 nm. By employing this workflow, we searched for low-energy configurations of large sulfuric acid-ammonia and sulfuric acid-dimethylamine clusters. We find that the binding free energies of clusters containing dimethylamine are unequivocally more stable than those of the ammonia-containing clusters. Our improved configurational sampling protocol can in the future be applied to study the growth and dynamics of large clusters of arbitrary compositions.

Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules

In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

Assessment of DLPNO-MP2 Approximations in Double-Hybrid DFT

The unfavorable scaling (N5) of the conventional second-order Møller-Plesset theory (MP2) typically prevents the application of double-hybrid (DH) density functionals to large systems with more than 100 atoms. A prominent approach to reduce the computational demand of electron correlation methods is the domain-based local pair natural orbital (DLPNO) approximation that is successfully used in the framework of DLPNO-CCSD(T). Its extension to MP2 [Pinski P.; Riplinger, C.; Valeev, E. F.; Neese, F. J. Chem. Phys. 2015, 143, 034108.] paved the way for DLPNO-based DH (DLPNO-DH) methods. In this work, we assess the accuracy of the DLPNO-DH approximation compared to conventional DHs on a large number of 7925 data points for thermochemistry and 239 data points for structural features, including main-group and transition-metal systems. It is shown that DLPNO-DH-DFT can be applied successfully to perform energy calculations and geometry optimizations for large molecules at a drastically reduced computational cost. Furthermore, PNO space extrapolation is shown to be applicable, similar to its DLPNO-CCSD(T) counterpart, to reduce the remaining error.

Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA

An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.

Metastable Charged Dimers in Organometallic Species: A Look into Hydrogen Bonding between Metallocene Derivatives

Multiply charged complexes bound by noncovalent interactions have been previously described in the literature, although they were mostly focused on organic and main group inorganic systems. In this work, we show that similar complexes can also be found for organometallic systems containing transition metals and deepen in the reasons behind the existence of these species. We have studied the structures, binding energies, and dissociation profiles in the gas phase of a series of charged hydrogen-bonded dimers of metallocene (Ru, Co, Rh, and Mn) derivatives isoelectronic with the ferrocene dimer. Our results indicate that the carboxylic acid-containing dimers are more strongly bonded and present larger barriers to dissociation than the amide ones and that the cationic complexes tend to be more stable than the anionic ones. Additionally, we describe for the first time the symmetric proton transfer that can occur while in the metastable phase. Finally, we use a density-based energy decomposition analysis to shine light on the nature of the interaction between the dimers.

Stereochemistry in All Its Shapes and Forms: The 56th Bürgenstock Conference

Acknowledging the crucial role of stereochemistry in fields as diverse as total synthesis, synthetic methodology, spectroscopy, and the study of the origin of life, the 56th SCS Conference on Stereochemistry, better known as the BÃ1/4rgenstock Conference, brought together a diverse range of chemistry expertise in Brunnen, Switzerland.

Theoretical investigation of anion perfluorocubane

CONTEXT

In this work, we did a theoretical exploration of C8F8 (Ib) and its anion radical analogue (IIb) in this work. By investigating the thermochemistry of electron capture, we find that the free energy associated with the conversion of C8H8 (Ia) into its anion radical analogue IIa is of the order of + 92.83 kcal.mol-1, while the conversion of Ib into IIb is - 6.42 kcal.mol-1. Therefore, species IIb is thermodynamically more stable than its neutral analogue. Natural bond orbitals (NBO) analyses revealed that compound Ib exhibits a relative electronic stability as a function of intramolecular delocalisations of the type [Formula: see text] of the order of 2.70 kcal.mol-1. Similar delocalizations for Ia are energetically lower (1.45 kcal.mol-1). Topological analyses of compounds Ib and IIb indicate that the addition of an electron to Ib enhances the covalency of the C-C bond, as can be seen by the reduction in the ellipticity of the C-C bond. The opposite is observed for Ia, whose addition of the electron (leading to IIa) reduces the covalency of the C-C bond. By comparing the free and packaged forms of the species, it is found that, in the crystalline form, the system will present greater relative stability due to the dispersive interactions involved, as evidenced by non-covalent interactions (NCI) analysis. Finally, it was possible to verify that the manifestation of the current density with a lower paratropic and less antiaromatic character in Ib and IIb point to C8F8 as a strong candidate for electron capture.

METHODS

Geometry optimization calculations were carried out, for all monomer structures using the hybrid functional B3LYP-D3 and the 6-31+G(d,p) basis set. To determine the formation thermochemistry of the ions, electronic energy corrections was performed using the DLPNO-CCSD(T)/aug-cc-pVTZ/C method. Starting from the optimised forms, shielding, nuclear magnetic resonance (NMR) spectra employing gauge-independent atomic orbital (GIAO), and NBO calculations were performed for these monomers, using the PBE0 functional and the pCSseg-2 atomic basis set. The magnetochemical analysis of ring currents was performed using the GIMIC formalism. For the topological analysis, it was applied the combination DLPNO-CCSD(T)/aug-cc-pVTZ/C, previously used for correcting the electronic energy.

Correlation Consistent Basis Sets for Explicitly Correlated Theory: The Transition Metals

We present correlation consistent basis sets for explicitly correlated (F12) calculations, denoted VnZ(-PP)-F12-wis (n = D,T), for the d-block elements. The cc-pVDZ-F12-wis basis set is contracted to [8s7p5d2f] for the 3d-block, while its ECP counterpart for the 4d and 5d-blocks, cc-pVDZ-PP-F12-wis, is contracted to [6s6p5d2f]. The corresponding contracted sizes for cc-pVTZ(-PP)-F12-wis are [9s8p6d3f2g] for the 3d-block elements and [7s7p6d3f2g] for the 4d and 5d-block elements. Our VnZ(-PP)-F12-wis basis sets are evaluated on challenging test sets for metal-organic barrier heights (MOBH35) and group-11 metal clusters (CUAGAU-2). In F12 calculations, they are found to be about as close to the complete basis set limit as the combination of standard cc-pVnZ-F12 on main-group elements with the standard aug-cc-pV(n+1)Z(-PP) basis sets on the transition metal(s). While our basis sets are somewhat more compact than aug-cc-pV(n+1)Z(-PP), the CPU time benefit is negligible for catalytic complexes that contain only one or two transition metals among dozens of main-group elements; however, it is somewhat more significant for metal clusters.

Spin-orbit charge transfer from guanine and 9-methylguanine radical cations to nitric oxide radicals and the induced triplet-to-singlet intersystem crossing

Nitric oxide (●NO) participates in many biological activities, including enhancing DNA radiosensitivity in ionizing radiation-based radiotherapy. To help understand the radiosensitization of ●NO, we report reaction dynamics between ●NO and the radical cations of guanine (a 9HG●+ conformer) and 9-methylguanine (9MG●+). On the basis of the formation of 9HG●+ and 9MG●+ in the gas phase and the collisions of the radical cations with ●NO in a guided-ion beam mass spectrometer, the charge transfer reactions of 9HG●+ and 9MG●+ with ●NO were examined. For both reactions, the kinetic energy-dependent product ion cross sections revealed a threshold energy that is 0.24 (or 0.37) eV above the 0 K product 9HG (or 9MG) + NO+ asymptote. To interrogate this abnormal threshold behavior, the reaction potential energy surface for [9MG + NO]+ was mapped out at closed-shell singlet, open-shell singlet, and triplet states using density functional and coupled cluster theories. The results showed that the charge transfer reaction requires the interaction of a triplet-state surface originating from a reactant-like precursor complex 3[9MG●+(↑)⋅(↑)●NO] with a closed-shell singlet-state surface evolving from a charge-transferred complex 1[9MG⋅NO+]. During the reaction, an electron is transferred from π∗(NO) to perpendicular π∗(9MG), which introduces a change in orbital angular momentum. The latter offsets the change in electron spin angular momentum and facilitates intersystem crossing. The reaction threshold in excess of the 0 K thermochemistry and the low charge-transfer efficiency are rationalized by the vibrational excitation in the product ion NO+ and the kinetic shift arising from a long-lived triplet intermediate.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Mechanistic Analysis of Alkyne Haloboration: A DFT, MP2, and DLPNO-CCSD(T) Study

Stereocontrol of the alkyne haloboration reaction has received attention in many experimental but few theoretical studies. Here we present a detailed quantum-chemical study of mechanisms leading to Z versus E isomers of haloboration products, considering acetylene and propyne combined with BCl3, BBr3, and BI3. Calculations using B3LYP-D3, MP2, and DLPNO-CCSD(T) methods are used to study polar reactions between the alkyne and BX3 in the absence and presence of an additional halide anion whose content in the reaction mixture can be controlled experimentally. The formation of anti-haloboration products via radical mechanisms is also explored, namely, by adding BX3 to (Z)-halovinyl radical. For the anti-haloboration of propyne, the radical route is prohibited by the regiochemistry of the initiating halopropenyl radical, while the polar route is unlikely due to a competitive allene generation. In contrast, energetically accessible routes exist for both syn- and anti-bromoboration of acetylene; hence, careful control of reaction conditions is necessary to steer the stereochemical outcome. Methodologically, MP2 results correspond better to the DLPNO-CCSD(T) energies than the B3LYP-D3 results in terms of both reaction barrier heights and relative ordering of energetically close stationary points.

On the Reactivity of Mes*P(PMe3 ) towards Aluminum(I) Compounds - Evidence for the Intermediate Formation of Phosphaalumenes

Phosphaalumenes are the heavier isoelectronic analogs of alkynes and have eluded facile synthesis until recently. We have reported that the combination of a phosphinidene transfer agent, Ar TerP(PMe3 ) (Ar Ter=2,6-Ar2 -C6 H3 ), with (Cp*Al)4 (Cp*=C5 (CH3 )5 ) afforded the phosphaalumenes Ar TerPAlCp* as isolable, violet, thermally stable compounds. In here we describe attempts to utilize Mes*P(PMe3 ) (Mes*=2,4,6-tBu3 -C6 H2 ) as a phosphinidene source in combination with different Al(I) precursors, namely Dip NacnacAl (Dip Nacnac=HC[C(Me)NDip]2 , Dip=2,6-iPr2 -C6 H3 ), (Cp*Al)4 and Cp3t Al (Cp3t =1,2,4-tBu3 -C5 H2 ). In all cases the formation of phosphaalumenes was not observed, however, their intermediate formation is indicated by formation of the dimer [Cp*Al(μ-PMes*)]2 (2) and C-H-bond activation products along the putative P=Al bond, giving unusual 1,2-P,Al-tetrahydronaphtalene derivatives 1 and 4, clearly underlining the role the sterically demanding group on phosphorus plays in these transformations. The reactivity studies are supported by theoretical studies, demonstrating a thermodynamic preference for the C-H activation products. Additionally, we show that there are potential pitfalls in the synthesis of Cp*2 AlH, the precursor to make (Cp*Al)4 and give recommendations how to circumvent these.

Barrier Heights for Diels-Alder Transition States Leading to Pentacyclic Adducts: A Benchmark Study of Crowded, Strained Transition States of Large Molecules

Theoretical characterization of reactions of complex molecules depends on providing consistent accuracy for the relative energies of intermediates and transition states. Here we employ the DLPNO-CCSD(T) method with core-valence correlation, large basis sets, and extrapolation to the CBS limit to provide benchmark values for Diels-Alder transition states leading to competitive strained pentacyclic adducts. We then used those benchmarks to test a diverse set of wave function and density functional methods for the absolute and relative barrier heights of these transition states. Our results show that only a few of the tested density functionals can predict the absolute barrier heights satisfactorily, although relative barrier heights are more accurate. The most accurate functionals tested are ωB97M-V, M11plus, ωB97X-V, PBE-D3(0), M11, and MN15 with MUDs from best estimates less than 3.0 kcal. These findings can guide selection of density functionals for future studies of crowded, strained transition states of large molecules.

An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics

We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm^-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.

The entropic penalty for associative reactions and their physical treatment during routine computations

A systematic study of the entropic penalty for associative reactions is presented. It is shown that computed solution-phase Gibbs free energies typically overestimate entropic contributions. This entropic penalty for associative reactions in solution, i.e., if the number of particles decreases along the reaction coordinate (sum of stoichiometric numbers ), originates from the insufficient treatment of entropic effects by implicit solvent models. We propose an additive correction scheme to Gibbs free energies that is suitable for routine applications by non-expert users. This correction is based on Garza's formalism for the solution-phase entropy [A. J. Garza, J. Chem. Theory Comput., 2019, 15, 3204.] that is physically sound and embedded into an efficient black-box type algorithm. To critically evaluate the entropic penalty and its proposed treatment, we compiled an experimental benchmark set of 31 Δ_rG and 22 in 15 different solvents. Using a representative best-practice computational protocol (at wave function theory (WFT) based DLPNO-CCSD(T) and density functional theory (DFT) based revDSD-PBEP86-D4 level with an implicit solvent model), we determined a sizeable entropic penalty ranging from 2-11 kcal mol^-1. Using the correction scheme presented herein, the entropic penalty is corrected to the chemical accuracy of ≤1 kcal mol^-1 (WFT and DFT). The same applies to at the WFT level. Barriers at the DFT level are overestimated by 2 kcal mol^-1 (classic) and underestimated by 2 kcal mol^-1 (corrected). This effect is attributed to the finding that barriers computed at the DFT level are systematically 2-3 kcal mol^-1 lower than barriers obtained with WFT.

Synthesis, Characterization, and Properties of High-Energy Fillers Derived from Nitroisobutylglycerol

Herein we report a comprehensive laboratory synthesis of a series of energetic azidonitrate derivatives (ANDP, SMX, AMDNNM, NIBTN, NPN, 2-nitro-1,3-dinitro-oxypropane) starting from the readily available nitroisobutylglycerol. This simple protocol allows obtaining the high-energy additives from the available precursor in yields higher than those reported using safe and simple operations not presented in previous works. A detailed characterization of the physical, chemical, and energetic properties including impact sensitivity and thermal behavior of these species was performed for the systematic evaluation and comparison of the corresponding class of energetic compounds.

On the Dynamic Behavior of Pacman Phosphanes─A Case of Cooperativity and Redox Isomerism

In solution, the Pacman chlorophosphane (2Cl) shows fast exchange of the endo/exo-orientation of the two P-Cl bonds in the molecule featuring cooperativity. Experimental and quantum mechanical investigations of the inversion on the phosphorus(III) centers reveal a crucial role of chloride ions in the dynamic process. To confirm the results, the homologous Pacman halogen-phosphanes 2X were prepared by halogen exchange reactions (X = F, Br, and I). Besides accelerated dynamic behavior for the heavier analogues, significant differences in the molecular structure are caused by the halogen exchange reactions, including the formation of an endo-endo substituted Pacman fluorophosphane as well as dicationic species by phosphorus halogen bond dissociation. The latter process can be regarded as redox isomerism since two P^III atoms in 2X become P^V centers in the dications.

Polyhedral Dicobaltadithiaboranes and Dicobaltdiselenaboranes as Examples of Bimetallic Nido Structures without Bridging Hydrogens

The geometries and energetics of the n-vertex polyhedral dicobaltadithiaboranes and dicobaltadiselenaboranes Cp2Co2E2Bn−4Hn−4 (E = S, Se; n = 8 to 12) have been investigated via the density functional theory. Most of the lowest-energy structures in these systems are generated from the (n + 1)-vertex most spherical closo deltahedra by removal of a single vertex, leading to a tetragonal, pentagonal, or hexagonal face depending on the degree of the vertex removed. In all of these low-energy structures, the chalcogen atoms are located at the vertices of the non-triangular face. Alternatively, the central polyhedron in most of the 12-vertex structures can be derived from a Co2E2B8 icosahedron with adjacent chalcogen (E) vertices by breaking the E–E edge and 1 or more E–B edges to create a hexagonal face. Examples of the arachno polyhedra with two tetragonal and/or pentagonal faces derived from the removal of two vertices from isocloso deltahedra were found among the set of lowest-energy Cp2Co2E2Bn−4Hn−4 (E = S, Se; n = 8 and 12) structures.

Electron Attachment to DNA: The Protective Role of Amino Acids

We have studied the effect of amino acids on the electron attachment properties of a DNA nucleobase, with cytosine as a model system. The equation of motion coupled cluster theory with an extended basis set has been used to simulate the electron-attached state of the DNA model system. Arginine, alanine, lysine, and glycine are the four amino acids considered to investigate their role in electron attachment to a DNA nucleobase. The electron attachment to cytosine in all the four cytosine-amino acid gas-phase dimer complexes follows a doorway mechanism, where the electron gets transferred from the initial dipole-bound doorway state to the final nucleobase-bound state through the mixing of electronic and nuclear degrees of freedom. When cytosine is bulk-solvated with glycine, the glycine-bound state acts as the doorway state, where the initial electron density is localized on the bulk amino acid and away from the nucleobase, thus leading to the physical shielding of the nucleobase from the incoming electron. At the same time, the presence of amino acids can increase the stability of the nucleobase-bound anionic state, which can suppress the sugar-phosphate bond rupture caused by dissociative electron attachment to DNA.

The perfluoroadamantoxy aluminate as an ideal weakly coordinating anion? - synthesis and first applications

Weakly coordinating anions (WCAs) facilitate the stabilization and isolation of highly reactive and almost "naked" cations. Alkoxyaluminate-based WCAs such as [Al(OC(CF₃)₃)₄]^- ([pf]^-) are widely used due to their synthetic accessibility and their high stability. However, small cations are still able to coordinate the oxygen atoms of the [pf]^- anion or even to abstract an alkoxy ligand. The novel WCA [Al(OC₁₀F₁₅)₄]^- ([pfAd]^-; OC₁₀F₁₅ = perfluoro-1-adamantoxy) is characterized by a very rigid core framework, thus indicating a higher stability towards fluoride-ion abstraction (DFT calculations) and providing hope to generate less disordered crystal structures. The [pfAd]^- anion was generated by the reaction of the highly acidic alcohol perfluoro-1-adamantanol C₁₀F₁₅OH with LiAlH₄ in o-DFB. Li[pfAd] could not be synthesized free of impurities (and still contains unreacted alcohol). Yet, starting from contaminated Li[pfAd], the very useful pure salts Ag[pfAd], [Ph₃C][pfAd] and [H(OEt₂)₂][pfAd] could be synthesized. The salts were characterized by NMR spectroscopy, single-crystal X-ray diffraction and IR spectroscopy. Additionally, [NO][pfAd] could be synthesized containing alcohol impurities but nonetheless enabled the synthesis of the salt P₉⁺[pfAd]^-. The synthesis of Tl[pfAd] in a mixture of H₂O/acetone/o-DFB demonstrated the water stability of the [pfAd]^- anion.

Aromaticity in P8 allotropes and (CH)8 analogues: significance of their 40 valence electrons?

The currently unknown phosphorus allotrope P₈ is of interest since its 40 total valence electrons is a "magic number" corresponding to a filled 1S²1P⁶1D¹⁰1S²1F¹⁴2P⁶ shell such as found in the relatively stable main group element clusters Al₁₃^- and Ge₉^4-. However, P₈ still remains as an elusive structure not realized experimentally. The lowest energy P₈ structure by a margin of ∼9 kcal mol^-1 is shown by density functional theory to be a cuneane analogue with no PP double bonds and two each of P₅, P₄, and P₃ rings. Higher energy P₈ structures are polycyclic systems having at most a single PP double bond. These P₈ systems are not "carbon copies" of the corresponding (CH)₈ hydrocarbons with exactly one hydrogen atom bonded to each carbon atom. Thus the lowest energy (CH)₈ structure is cyclooctatetraene with four CC bonds followed by benzocyclobutene with three CC bonds. The cuneane (CH)₈ structure is a relatively high energy isomer lying ∼36 kcal mol^-1 above cyclooctatetraene. The cubane P₈ and (CH)₈ structures are even higher energy structures, lying ∼37 and ∼74 kcal mol^-1 in energy above the corresponding global minima. Our results demonstrate differences in medium sized aggregates of elemental phosphorus and isolobal hydrocarbon species.

A new look at the structure of the neutral Au18 cluster: hollow versus filled golden cage

The geometry of the neutral Au₁₈ gold cluster was probed by a combination of quantum chemical calculations and far-infrared multiple photon dissociation (FIR-MPD) spectroscopy of a Kr messenger complex. Two low-lying isomers are identified to potentially contribute to the experimental IR spectrum, both being derived from a star-like Au₁₇ structure upon capping with one extra Au atom either inside (18_1) or outside (18_5) the star. In particular, the present detection of structure 18_1 by DFT computations where a golden cage encapsulates an endohedral Au atom, is intriguing as a stable core-shell isomer has, to our knowledge, never been found before for such small neutral gold clusters. DFT and local coupled-cluster (DLPNO and PNO-CCSD(T)) computations indicate that both Au₁₈ isomers are close to each other, within ∼3 kcal mol^-1, on the energy scale. Although the exact energy ordering is again method-dependent and remains, at present, inconclusive, the most striking spectral signatures of both isomers are related to vibrational modes localized at atoms capping the inner pentaprism sub-structure that result in prominent peaks centered at ∼80 cm^-1, close to the most prominent experimental feature found at 78 cm^-1. The calculated IR spectra of both core-shell and hollow isomers are very similar to each other and both agree comparably well with the experimental FIR-MPD spectra of the Au₁₈Kr_1,2 complexes.

Tensor Hypercontraction Form of the Perturbative Triples Energy in Coupled-Cluster Theory

We present the working equations for a reduced-scaling method of evaluating the perturbative triples (T) energy in coupled-cluster theory, through the tensor hypercontraction (THC) of the triples amplitudes (t_ijk^abc). Through our method, we can reduce the scaling of the (T) energy from the traditional O(N7) to a more modest O(N5). We also discuss implementation details to aid future research, development, and software realization of this method. Additionally, we show that this method yields submillihartree (mEh) differences from CCSD(T) when evaluating absolute energies and sub-0.1 kcal/mol energy differences when evaluating relative energies. Finally, we demonstrate that this method converges to the true CCSD(T) energy through the systematic increasing of the rank or eigenvalue tolerance of the orthogonal projector, as well as exhibiting sublinear to linear error growth with respect to system size.