An efficient and near linear scaling pair natural orbital based local coupled cluster method

Computational investigations of stable multiple-cage-occupancy He clathrate-like hydrostructures

One of the several possibilities offered by the interesting clathrate hydrates is the opportunity to encapsulate several atoms or molecules, in such a way that more efficient storage materials could be explored or new molecules that otherwise do not exist could be created. These types of applications are receiving growing attention from technologists and chemists, given the future positive implications that they entail. In this context, we investigated the multiple cage occupancy of helium clathrate hydrates, to establish stable novel hydrate structures or ones similar to those predicted previously by experimental and theoretical studies. To this purpose, we analyzed the feasibility of including an increased number of He atoms inside the small (D) and large (H) cages of the sII structure through first-principles properly assessed density functional approaches. On the one hand, we have computed energetic and structural properties, in which we examined the guest-host and guest-guest interactions in both individual and two-adjacent clathrate-like sII cages by means of binding and evaporation energies. On the other hand, we have carried out a thermodynamical analysis on the stability of such He-containing hydrostructures in terms of changes in enthalpy, ΔH, Gibbs free energy, ΔG, and entropy, ΔS, during their formation process at various temperature and pressure values. In this way, we have been able to make a comparison with experiments, reaffirming the ability of computational DFT approaches to describe such weak guest-host interactions. In principle, the most stable structure involves the encapsulation of one and four He atoms inside the D and H sII cages, respectively; however, more He atoms could be entrapped under lower temperature and/or higher pressure thermodynamic conditions. We foresee such accurate computational quantum chemistry approaches contributing to the current emerging machine-learning model development.

Validation of Field-Dependent Ion-Solvent Cluster Modeling via Direct Measurement of Cluster Size Distributions

Ion mobility spectrometry is widely used in analytical chemistry, either as a stand-alone technique or coupled to mass spectrometry. Ions in the gas phase tend to form loosely bound clusters with surrounding solvent vapors, artificially increasing the collisional cross section and the mass of the ion. This, in turn, affects ion mobility and influences separation. Further, ion-solvent clusters play an important role in most ionization mechanisms occurring in the gas phase. Consequently, a deeper understanding of ion-solvent cluster association and dissociation processes is desirable to guide experimental design and interpretation. A few computational models exist, which aim to describe the amount of clustering as a function of the reduced electric field strength, bath gas pressure and temperature, and the chemical species probed. It is especially challenging to model ion mobility under high reduced electrical field strengths due to the nonthermal conditions created by the field. In this work, we aim to validate a recently proposed first-principles model by comparing its predictions with direct measurements of cluster size distributions over a range of 20-120 Td as observed using a High Kinetic Energy Ion Mobility Spectrometer coupled to a mass spectrometer (HiKE-IMS-MS). By studying H⁺(H₂O)_n, [MeOH + H + n(H₂O)]⁺, [ACE + H + n(H₂O)]⁺, and [PhNH₂ + H + n(H₂O)]⁺ as test systems, we find very good agreement between model and experiment, supporting the validity of the computational workflow. Further, the detailed information gained from the modeling yields important insights into the cluster dynamics within the HiKE-IMS, allowing for better interpretation of the measured ion mobility spectra.

Limited Role of Malonic Acid in Sulfuric Acid-Dimethylamine New Particle Formation

Aerosols play an important role in climate and air quality; however, the mechanisms behind aerosol particle formation in the atmosphere are poorly understood. Studies have identified sulfuric acid, water, oxidized organics, and ammonia/amines as key precursors for forming aerosol particles in the atmosphere. Theoretical and experimental investigations have indicated that other species, such as organic acids, may be involved in atmospheric nucleation and growth of freshly formed aerosol particles. Organic acids, such as dicarboxylic acids, which are abundant in the atmosphere, have been measured in ultrafine aerosol particles. These observations suggest that organic acids may contribute to new particle formation in the atmosphere but their role remains ambiguous. This study examines how malonic acid interacts with sulfuric acid and dimethylamine to form new particles at warm boundary layer conditions using experimental observations from a laminar flow reactor and quantum chemical calculations coupled with cluster dynamics simulations. Observations reveal that malonic acid does not contribute to the initial steps (formation of <1 nm diameter particle) of nucleation with sulfuric acid-dimethylamine. In addition, malonic acid was found to not participate in the subsequent growth of the freshly nucleated 1 nm particles from sulfuric acid-dimethylamine reactions to diameters of 2 nm.

First-Principles Modeling of Preferential Solvation in Mixed-Modifier Differential Mobility Spectrometry

Differential mobility spectrometry (DMS) separates ions based on mobility differences between high and low electric field conditions. To enhance resolution, solvents such as water and acetonitrile are often used to modify the collision environment and take advantage of differing dynamic clustering behavior between analytes that coelute in hard-sphere environments (e.g., N₂). When binary solvent mixtures are used to modify the DMS environment, one solvent can have a dominant influence over the other with respect to ion trajectories. For example, for quinoline derivatives, a 9:1 water:acetonitrile solvent mixture exhibited identical behavior to an environment containing only acetonitrile as a modifier. It was hypothesized that this effect arises due to the significantly different binding strengths of the two solvents. Here, we utilize a first-principles model of DMS to study analytes in single and binary solvent mixtures and explore the effects governing the dominance of one solvent over the other. Computed DMS dispersion curves of quinoline derivatives are in excellent agreement with those measured experimentally. For mixed-modifier environments, the predicted cluster populations show a clear preferential solvation of ions with the stronger binding solvent. The influence of ion-solvent binding energies, solvent concentration, and solvent molecule size is discussed in the context of the observed DMS behavior. This work can guide the usage of binary solvent mixtures for improving ion separations in cases where compounds coelute in pure N₂ and in single-solvent modifier environments. Moreover, our results indicate that binary solvent mixtures can be used to create a relative scale for solvent binding energies.

In-depth unveiling the interfacial adsorption mechanism of triazine derivatives as corrosion inhibitors for carbon steel in carbon dioxide saturated oilfield produced water

In this work, two triazine derivatives (BTT-1 and BTT-2) were synthesized by the simple one-step condensation of three components and used as high-efficient corrosion inhibitors to deal with the corrosion issue of carbon steel (CS) in petroleum industry. Electrochemical tests indicate that both BTT-1 and BTT-2 present superior inhibition performance with the inhibition efficiency of 97.9 % and 98.4 % at a low concentration of 0.18 mM, respectively. Quantum chemical calculations reveal that compared to BTT-1 molecule with a butyl chain, the introduction of benzyl group endows BTT-2 molecule with more adsorption sites, which favors the adsorption of BTT-2 molecule on CS surface. Furthermore, the GFN-xTB calculations demonstrate that BTT-1 and BTT-2 could adsorb on CS surface through the formation of Fe-N and Fe-S bonds. Compared to BTT-1, BTT-2 exhibits stronger adsorption on CS surface by forming more and shorter bonds with a more negative adsorption energy, which accounts for the better inhibitive performance of BTT-2.

Current and future machine learning approaches for modeling atmospheric cluster formation

The formation of strongly bound atmospheric molecular clusters is the first step towards forming new aerosol particles. Recent advances in the application of machine learning models open an enormous opportunity for complementing expensive quantum chemical calculations with efficient machine learning predictions. In this Perspective, we present how data-driven approaches can be applied to accelerate cluster configurational sampling, thereby greatly increasing the number of chemically relevant systems that can be covered.

Quantifying the Intrinsic Strength of C-H⋯O Intermolecular Interactions

It has been recognized that the C-H⋯O structural motif can be present in destabilizing as well as highly stabilizing intermolecular environments. Thus, it should be of interest to describe the strength of the C-H⋯O hydrogen bond for constant structural factors so that this intrinsic strength can be quantified and compared to other types of interactions. This description is provided here for C_2h-symmetric dimers of acrylic acid by means of the calculations that employ the coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] together with an extrapolation to the complete basis set (CBS) limit. Dimers featuring the C-H⋯O and O-H⋯O hydrogens bonds are carefully investigated in a wide range of intermolecular separations by the CCSD(T)/CBS approach, and also by the symmetry-adapted perturbation theory (SAPT) method, which is based on the density-functional theory (DFT) treatment of monomers. While the nature of these two types of hydrogen bonding is very similar according to the SAPT-DFT/CBS calculations and on the basis of a comparison of the intermolecular potential curves, the intrinsic strength of the C-H⋯O interaction is found to be about a quarter of its O-H⋯O counterpart that is less than one might anticipate.

Measuring Electron Correlation: The Impact of Symmetry and Orbital Transformations

In this perspective, the various measures of electron correlation used in wave function theory, density functional theory and quantum information theory are briefly reviewed. We then focus on a more traditional metric based on dominant weights in the full configuration solution and discuss its behavior with respect to the choice of the N-electron and the one-electron basis. The impact of symmetry is discussed, and we emphasize that the distinction among determinants, configuration state functions and configurations as reference functions is useful because the latter incorporate spin-coupling into the reference and should thus reduce the complexity of the wave function expansion. The corresponding notions of single determinant, single spin-coupling and single configuration wave functions are discussed and the effect of orbital rotations on the multireference character is reviewed by analyzing a simple model system. In molecular systems, the extent of correlation effects should be limited by finite system size and in most cases the appropriate choices of one-electron and N-electron bases should be able to incorporate these into a low-complexity reference function, often a single configurational one.

No Transition Metals Required - Oxygen Promoted Synthesis of Imines from Primary Alcohols and Amines under Ambient Conditions

The synthesis of imines denotes a cornerstone in organic chemistry. The use of alcohols as renewable substituents for carbonyl-functionality represents an attractive opportunity. Consequently, carbonyl moieties can be in situ generated from alcohols upon transition-metal catalysis under inert atmosphere. Alternatively, bases can be utilized under aerobic conditions. In this context, we report the synthesis of imines from benzyl alcohols and anilines, promoted by KOt Bu under aerobic conditions at room temperature, in the absence of any transition-metal catalyst. A detailed investigation of the radical mechanism of the underlying reaction is presented. This reveals a complex reaction network fully supporting the experimental findings.

Algebraic Diagrammatic Construction Theory for Simulating Charged Excited States and Photoelectron Spectra

Charged excitations are electronic transitions that involve a change in the total charge of a molecule or material. Understanding the properties and reactivity of charged species requires insights from theoretical calculations that can accurately describe orbital relaxation and electron correlation effects in open-shell electronic states. In this Review, we describe the current state of algebraic diagrammatic construction (ADC) theory for simulating charged excitations and its recent developments. We start with a short overview of ADC formalism for the one-particle Green's function, including its single- and multireference formulations and extension to periodic systems. Next, we focus on the capabilities of ADC methods and discuss recent findings about their accuracy for calculating a wide range of excited-state properties. We conclude our Review by outlining possible directions for future developments of this theoretical approach.

A quest for ideal electric field-driven MX@C70 endohedral fullerene memristors: which MX fits the best?

Endohedral fullerenes with a dipolar molecule enclosed in the fullerene cage have great potential in molecular electronics, such as diodes, switches, or molecular memristors. Here, we study a series of model systems based on MX@D_5h(1)-C₇₀ (M = a metal or hydrogen, X = a halogen or a chalcogen) endohedral fullerenes to identify potential molecular memristors and to derive a general formula for rapid identification of potential memristors among analogous MX@C_n systems. To obtain sufficiently accurate results for switching barriers and encapsulation energies, we perform a benchmark of ten DFT functionals against ab initio SCS-MP2 and DLPNO-CCSD(T) methods at the complete basis set limit. The whole series is then investigated using the PBE0 functional which was found to be the most efficient vs. the ab initio methods. Nine of the 34 MX@C₇₀ molecules studied are predicted to have suitable switching barriers to be considered as potential candidates for molecular switches and memristors. We have identified several structure-property relationships for the switching barrier and response of the systems to the electric field, in particular the dependence of the switching barrier on the available space for M-X switching and faster response of the system to the electric field with a larger dipole moment of MX and MX@C₇₀.

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

Microsolvated clusters of multiply charged anions play a crucial role in atmospheric chemistry and some of them were previously registered experimentally. At the same time, there are no experimental observations of [CO₃·(H₂O)_n]^2-. The reasons for this may be related to the thermodynamical or kinetical instability of microsolvated CO₃^2- toward autoionization or autoprotonation processes. In this study we theoretically investigate the potential stability of the [CO₃·(H₂O)_n]^2- microsolvated clusters from both perspectives - thermodynamic and kinetic - and we claim they are stable toward autoionization and kinetically semi-stable toward autoprotonation. In addition, the behaviour of CO₃^2- anions in bulk water solvent was analysed to highlight important precautions for synthetic purposes.

Effective Reconstruction of Expectation Values from Ab Initio Quantum Embedding

Quantum embedding is an appealing route to fragment a large interacting quantum system into several smaller auxiliary "cluster" problems to exploit the locality of the correlated physics. In this work, we critically review approaches to recombine these fragmented solutions in order to compute nonlocal expectation values, including the total energy. Starting from the democratic partitioning of expectation values used in density matrix embedding theory, we motivate and develop a number of alternative approaches, numerically demonstrating their efficiency and improved accuracy as a function of increasing cluster size for both energetics and nonlocal two-body observables in molecular and solid state systems. These approaches consider the N-representability of the resulting expectation values via an implicit global wave function across the clusters, as well as the importance of including contributions to expectation values spanning multiple fragments simultaneously, thereby alleviating the fundamental locality approximation of the embedding. We clearly demonstrate the value of these introduced functionals for reliable extraction of observables and robust and systematic convergence as the cluster size increases, allowing for significantly smaller clusters to be used for a desired accuracy compared to traditional approaches in ab initio wave function quantum embedding.

Modeling coarse-grained van der Waals interactions using dipole-coupled anisotropic quantum Drude oscillators

The Quantum Drude Oscillator (QDO) model is a promising candidate for accurately calculating the van der Waals (vdW) interaction. Anisotropic QDO models have recently been used to represent quantum fluctuations of molecular fragments rather than that of single atoms. While this model promises accurate calculation of vdW energy, there is significant room for improvements, such as incorporating a proper fragmentation method, higher-order dispersion corrections, and so forth. The present work attempts to gauge dipole-dipole interactions' ability without fragmentation. A suitable anisotropic damping function is also introduced to work with anisotropic QDO. This revised model accurately predicts the binding energies of vdW complexes for most of the systems considered. This work indicates the limit of dipole approximation for an anisotropic QDO-based model.

New insights into the mechanism of synergetic photoredox/copper(i)-catalyzed carbocyanation of 1,3-dienes: a DFT study

This work presents a DFT-based computational study to understand the mechanism, and regio- and enantioselectivities in the synergetic photoredox/copper(i)-catalyzed carbocyanation of 1,3-dienes with alkyl redox-active esters. The calculated results show an unprecedented copper catalytic mechanism, where the reaction follows a catalytic cycle involving CuI-only catalysis, instead of a Cu(i)/Cu(ii)/Cu(iii)/Cu(i) cycle as proposed in the experimental study. Moreover, it is found that the critical step involves the reaction of the cyanocopper(i) species with an allyl cation rather than the cyanocopper(ii) species reacting with an allyl radical as proposed in the experiment, and that the photocatalyst is regenerated via single electron transfer from the allyl radical to the oxidized photocatalyst. In the newly proposed photoredox/copper(i) catalysis, the reaction consists of four stages: (i) generation of the copper(i) active catalyst, (ii) formation of an allyl radical with oxidative quenching of the photoexcited species, (iii) generation of an allylcopper complex accompanied by the regeneration of the photocatalyst, and (iv) formation of the allyl cyanide product with the regeneration of the copper(i) active catalyst. The cyanation of the allyl cation is calculated to be the regio- and enantioselectivity-determining step. The excellent regio- and stereoselectivities are attributed to the favorable CH-π interaction between the substrate and catalyst as well as the small distortion of the substrate.

Kinetic, Thermodynamic, and Dynamic Control in Normal vs. Cross [2 + 2] Cycloadditions of Ene-Keteniminium Ions: Computational Understanding, Prediction, and Experimental Verification

Almost all reported intramolecular [2 + 2] reactions of ene-keteniminium ions gave normal [2 + 2] products with a fused bicycle framework, but not cross [2 + 2] products with a bicyclo[3.1.1]heptane skeleton, a highly pursued bioisostere in pharmaceutical chemistry. How to rationalize this and design new cross [2 + 2] reactions? Theoretical studies using density functional theory, high-level ab initio single-point energy calculations, and molecular dynamics showed that this [2 + 2] reaction has all three patterns of regiochemical control: the reaction is controlled either kinetically, thermodynamically, or dynamically. A carbocation model of forming endo and exo carbocations has been proposed to rationalize the reaction outcomes, revealing that the tethers (between alkenes and keteniminium ions), substituents (on the alkenes), and alkene configurations in ene-keteniminium ions play critical roles. These understandings were further used to predict that introducing a substituent in the terminal position of alkene with a trans configuration in ene-keteniminium ions can realize the cross [2 + 2] reaction, which is dynamically controlled for alkyl substituents or kinetically controlled for aryl substituents. These and more other predictions were realized experimentally, and many cross [2 + 2] products with a bicyclo[3.1.1]heptane skeleton can be achieved. Both molecular dynamics and new experiments have also been applied to correct a key but misassigned [2 + 2] product reported in the literature, further supporting the insightful mechanisms reported here.

Additivity of Diene Substituent Gibbs Free Energy Contributions for Diels-Alder Reactions between Me2C=CMe2 and Substituted Cyclopentadienes

Systematic computational studies of pericyclic Diels-Alder reactions between (H₃C)₂C=C(CH₃)₂, 1, and all permutations of substituted cyclopentadienes c-C₅R¹R²R³R⁴R^5aR^5b (R = H, CH₃, CF₃, F) allowed isolation of substitutional effects on Gibbs free energy barrier heights and reaction Gibbs free energies. "Average Substitution Gibbs Free Energy Correction" ΔG _ASC# ^‡/ΔG _ASC# values for each substituent in each position appeared to be additive. Substituent effects on barriers showed interesting contrasts. Methyl substitution at positions 5a and 5b increased barriers significantly, while substitution at all other positions had essentially no impact. In contrast, fluoro substitution at positions 5a and 5b lowered barriers more than substitution at other positions. Trifluoromethyl substitution mixed these effects, in that substitution at positions 5a and 5b increased barriers, but substitution at other positions lowered them. Despite the variances, ΔG _ASC# ^‡/ΔG _ASC# values allowed reliable prediction of barriers and exergonicities for reactions between 1 and highly substituted cyclopentadienes, and between 1 and cyclopentadienes with random mixtures of CH₃/CF₃/F substituents. ΔG _ASC# ^‡/ΔG _ASC# values were correlated with steric considerations and quantum theory of atoms in molecules (QTAIM) calculations. Overall, the ASC values provide a resource for predicting which Diels-Alder reactions of this type should occur at rapid rates and/or give stable bicyclic products.

Spin-Forbidden Addition of Molecular Oxygen to Stable Enol Intermediates-Decarboxylation of 2-Methyl-1-tetralone-2-carboxylic Acid

The deprotonation of an organic substrate is a common preactivation step for the enzymatic cofactorless addition of O2 to this substrate, as it promotes charge-transfer between the two partners, inducing intersystem crossing between the triplet and singlet states involved in the process. Nevertheless, the spin-forbidden addition of O2 to uncharged ligands has also been observed in the laboratory, and the detailed mechanism of how the system circumvents the spin-forbiddenness of the reaction is still unknown. One of these examples is the cofactorless peroxidation of 2-methyl-3,4-dihydro-1-naphthol, which will be studied computationally using single and multi-reference electronic structure calculations. Our results show that the preferred mechanism is that in which O2 picks a proton from the substrate in the triplet state, and subsequently hops to the singlet state in which the product is stable. For this reaction, the formation of the radical pair is associated with a higher barrier than that associated with the intersystem crossing, even though the absence of the negative charge leads to relatively small values of the spin-orbit coupling.

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.

Investigating the Thermodynamics and Kinetics of Catechin Pyrolysis for Environmentally Friendly Binders

The thermodynamics and kinetics of the pyrolysis of (+)-catechin, a building block of the condensed tannins found in recipes for sustainable binders, are evaluated at the DLPNO-CCSD(T) level and compared to other methods from quantum chemistry. Using the climbing image nudged elastic band method coupled with transition state optimization, minimum energy paths and highest-energy transition states are identified for the first two pyrolysis steps, a catechol split-off with subsequent dehydrogenation. While the catechol split-off path was very smooth, the dehydrogenation featured an additional transition state in the form of an OH group rotation. The combined reaction was judged endothermic in the range of 0 to 1250 K and exergonic at 1000 K and above. It is shown that the catechol split-off is the rate-determining step of the pyrolysis of catechin, which is equivalent to kinetic inhibition at all investigated temperatures.

Quasi-periodic scattering of topological edge states induced by the vacancies in chloridized gallium bismuthide nanoribbons

The chloridized gallium bismuthide was predicted to be a two-dimensional topological insulator with large topological band gap. It may be beneficial for achieving the quantum spin Hall effect and its related applications at high temperatures. To better understand the quantum transport in topological nanoribbons, we investigated the effect of vacancy on the quantum transport of topological edge states in the armchair chloridized gallium bismuthide nanoribbons by combining density functional theory and nonequilibrium Green's function. The results suggest the vacancies at center are more likely to cause the scattering of topological edge states. The average scattering is insensitive to the enlargement of vacancy along the transport direction. More interestingly, the obvious scattering of topological edge states can only be found at some special energies, and these special energies are distributed quasi-periodically. The quasi-periodic scattering may be used as a kind of fingerprint of vacancies. Our studies may be helpful for the application of topological nanoribbons.

Water Charge Transfer Accelerates Criegee Intermediate Reaction with H2O- Radical Anion at the Aqueous Interface

Criegee intermediates (CIs) are important carbonyl oxides that may react with atmospheric trace chemicals and impact the global climate. The CI reaction with water has been widely studied and is a main channel for trapping CIs in the troposphere. Previous experimental and computational reports have largely focused on reaction kinetic processes in various CI-water reactions. The molecular-level origin of CI's interfacial reactivity at the water microdroplet surface (e.g., as found in aerosols and clouds) is unclear. In this study, by employing the quantum mechanical/molecular mechanical (QM/MM) Born-Oppenheimer molecular dynamics with the local second-order Møller-Plesset perturbation theory, our computational results reveal a substantial water charge transfer up to ∼20% per water, which creates the surface H₂O⁺/H₂O^- radical pairs to enhance the CH₂OO and anti-CH₃CHOO reactivity with water: the resulting strong CI-H₂O^- electrostatic attraction at the microdroplet surface facilitates the nucleophilic attack to the CI carbonyl by water, which may counteract the apolar hindrance of the substituent to accelerate the CI-water reaction. Our statistical analysis of the molecular dynamics trajectories further resolves a relatively long-lived bound CI(H₂O^-) intermediate state at the air/water interface, which has not been observed in gaseous CI reactions. This work provides insights into what may alter the oxidizing power of the troposphere by the next larger CIs than simple CH₂OO and implicates a new perspective on the role of interfacial water charge transfer in accelerating molecular reactions at aqueous interfaces.

Secondary Electron Attachment-Induced Radiation Damage to Genetic Materials

Reactions of radiation-produced secondary electrons (SEs) with biomacromolecules (e.g., DNA) are considered one of the primary causes of radiation-induced cell death. In this Review, we summarize the latest developments in the modeling of SE attachment-induced radiation damage. The initial attachment of electrons to genetic materials has traditionally been attributed to the temporary bound or resonance states. Recent studies have, however, indicated an alternative possibility with two steps. First, the dipole-bound states act as a doorway for electron capture. Subsequently, the electron gets transferred to the valence-bound state, in which the electron is localized on the nucleobase. The transfer from the dipole-bound to valence-bound state happens through a mixing of electronic and nuclear degrees of freedom. In the presence of aqueous media, the water-bound states act as the doorway state, which is similar to that of the presolvated electron. Electron transfer from the initial doorway state to the nucleobase-bound state in the presence of bulk aqueous media happens on an ultrafast time scale, and it can account for the decrease in DNA strand breaks in aqueous environments. Analyses of the theoretically obtained results along with experimental data have also been discussed.

Identification of Ge≡O Triple Bond in Ge6O- Cluster: Anion Photoelectron Spectroscopy and Theoretical Calculations

Unlike C≡O, which is common in coordination chemistry and organometallic chemistry, little is known about Si≡O or Ge≡O compounds. Here we report a Ge₆O^- cluster featuring a Ge≡O triple bond. The structural and chemical bonding properties of Ge₆O^-/0 are investigated using anion photoelectron spectroscopy and theoretical calculations. Two nearly degenerate isomers have been found for Ge₆O^-. The lowest-energy structure (6A) can be viewed as an O atom bonding with a tetragonal bipyramidal Ge₆. The second one (6B) can be considered as an O atom interacting with a capped trigonal bipyramidal Ge₆. Chemical bonding analyses reveal that Ge₆O^- (6A) can be viewed as a Ge≡O unit interacting with a σ antiaromatic C_2v symmetric tetragonal pyramidal Ge₅^3- moiety. Comparisons of the chemical bonding in Ge₆O^- (6A) with that in Ge₅CO^- and Ge₅MnO^- indicate the similar behavior of Ge≡O to C≡O and Mn≡O in its bonding to the Ge₅^3- and Ge₅^4- moieties.

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.

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.

Performance of Density Functionals and Semiempirical 3c Methods for Small Gold-Thiolate Clusters

In light of the recent surge in computational studies of gold thiolate clusters, we present a comparison of popular density functionals (DFAs) and three-part corrected methods (3c-methods) on their performance by taking a data set named as AuSR18 consisting of 18 isomers of Au_n(SCH₃)_m (m ≤ n = 1-3). We have compared the efficiency and accuracy of the DFAs and 3c-methods in geometry optimization with RI-SCS-MP2 as the reference method. Similarly, the performance for accurate and efficient energy evaluation was compared with DLPNO-CCSD(T) as the reference method. The lowest energy structure among the isomers of the largest stoichiometry from our data set, AuSR18, i.e., Au₃(SCH₃)₃, is considered to evaluate the computational time for SCF and gradient evaluations. Alongside this, the numbers of optimization steps to locate the most stable minima of Au₃(SCH₃)₃ are compared to assess the efficiency of the methods. A comparison of relevant bond lengths with the reference geometries was made to estimate the accuracy in geometry optimization. Some methods, such as LC-BLYP, ωB97M-D3BJ, M06-2X, and PBEh-3c, could not locate many of the minima found by most of the other methods; thus, the versatility in locating various minima is also an important criterion in choosing a method for the given project. To determine the accuracy of the methods, we compared the relative energies of the isomers in each stoichiometry and the interaction energy of the gold core with the ligands. The dependence of basis set size and relativistic effects on energies are also compared. The following are some of the highlights. TPSS has shown accuracy, while mPWPW shows comparable speed and accuracy. For the relative energies of the clusters, the hybrid range-separated DFAs are the best option. CAM-B3LYP excels, whereas B3LYP performs poorly. Overall, LC-BLYP is a balanced performer considering both the geometry and relative stability of the structures, but it lacks diversity. The 3c-methods, although fast, are less impressive in relative stability.

Catalytic artificial nitroalkane oxidases - a way towards organocatalytic umpolung

Nitroalkane oxidases (NAOs) are flavoenzymes that catalyse the oxidation of nitroalkanes to their corresponding carbonyl compounds while producing nitrite anions. Herein, we present an artificial catalytic system using flavins or ethylene-bridged flavinium salts that works via an NAO-like process. Under conditions optimised in terms of solvent, base, temperature and oxygen pressure, primary nitroalkanes were transformed to aldehydes. In our system, aldehydes immediately reacted with other nitroalkane molecules to form β-nitroalcohols. The reduced flavin catalyst was re-oxidised by oxygen. An alternative mechanism towards β-nitroalcohols via 5-(2-nitrobutyl)-1,5-dihydroflavin was suggested through quantum chemical calculations and by trapping and characterising this dihydroflavin intermediate. Interestingly, 5-(2-nitrobutyl)-1,5-dihydroflavin is an analogue of the flavin adenine dinucleotide adduct previously observed in an NAO X-ray structure. In both mechanistic pathways, flavin-5-iminium species is formed by nitroalkanide addition to flavin. This process represents flavin-based umpolung of an original donor to an acceptor.

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.

Exploring the Accuracy Limits of PNO-Based Local Coupled-Cluster Calculations for Transition-Metal Complexes

While the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples (DLPNO-CCSD(T)) has proven instrumental for computing energies and properties of large and complex systems accurately, calculations on first-row transition metals with a complex electronic structure remain challenging. In this work, we identify and address the two main error sources that influence the DLPNO-CCSD(T) accuracy in this context, namely, (i) correlation effects from the 3s and 3p semicore orbitals and (ii) dynamic correlation-induced orbital relaxation (DCIOR) effects that are not described by the local MP2 guess. We present a computational strategy that allows us to completely eliminate the DLPNO error associated with semicore correlation effects, while increasing, at the same time, the efficiency of the method. As regards the DCIOR effects, we introduce a diagnostic for estimating the deviation between DLPNO-CCSD(T) and canonical CCSD(T) for systems with significant orbital relaxation.

Clusteromics V: Organic Enhanced Atmospheric Cluster Formation

Formic acid (FA) is a prominent candidate for organic enhanced nucleation due to its high abundance and stabilizing effect on smaller clusters. Its role in new particle formation is studied through the use of state-of-the-art quantum chemical methods on the cluster systems (acid)_1-2(FA)₁(base)_1-2 with the acids being sulfuric acid (SA)/methanesulfonic acid (MSA) and the bases consisting of ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). A funneling approach is used to determine the cluster structures with initial configurations generated through the ABCluster program, followed by semiempirical PM7 and ωB97X-D/6-31++G(d,p) calculations. The final binding free energy is calculated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory using the quasi-harmonic approximation. Cluster dynamics simulations show that FA has a minuscule or negligible effect on the MSA-FA-base systems as well as most of the SA-FA-base systems. The SA-FA-DMA cluster system shows the highest influence from FA with an enhancement of 21%, compared to its non-FA counterpart.

Formation of 3-Oxa- and 3-Thiacyclohexyne from Ring Expansion of Heterocyclic Alkylidene Carbenes: A Mechanistic Study

The rearrangement pathways of two alkylidene carbenes appended to an oxa or thiacyclopentane into the corresponding heterocyclohexynes were elucidated using 13C-labeling experiments. Both carbenes exhibited a preference for migration of the allylic carbon bound to the heteroatom. Anomeric interactions involving a heteroatom lone pair and antibonding orbital of the migrating bond and inductive destabilization of the minor migratory pathway are discussed as plausible reasons for the observed trends.

Atmospheric Sulfuric Acid-Multi-Base New Particle Formation Revealed through Quantum Chemistry Enhanced by Machine Learning

The formation of molecular clusters and secondary aerosols in the atmosphere has a significant impact on the climate. Studies typically focus on the new particle formation (NPF) of sulfuric acid (SA) with a single base molecule (e.g., dimethylamine or ammonia). In this work, we examine the combinations and synergy of several bases. Specifically, we used computational quantum chemistry to perform configurational sampling (CS) of (SA)_0-4(base)_0-4 clusters with five different types of bases: ammonia (AM), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). Overall, we studied 316 different clusters. We used a traditional multilevel funnelling sampling approach augmented by a machine-learning (ML) step. The ML made the CS of these clusters possible by significantly enhancing the speed and quality of the search for the lowest free energy configurations. Subsequently, the cluster thermodynamics properties were evaluated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory. The calculated binding free energies were used to evaluate the cluster stabilities for population dynamics simulations. The resultant SA-driven NPF rates and synergies of the studied bases are presented to show that DMA and EDA act as nucleators (although EDA becomes weak in large clusters), TMA acts as a catalyzer, and AM/MA is often overshadowed by strong bases.

The role of hydration in atmospheric salt particle formation

New-particle formation from condensable acid and base molecules is a ubiquitous phenomenon in the atmosphere. The role of water in salt particle formation is not fully understood as it can stabilize or destabilize cluster structures, which leads to non-linear effects on cluster formation dynamics. In the studied systems, increased relative humidity can enhance the particle formation for up to four orders of magnitude in the case of nitric acid, but it can also slightly reduce the particle formation in the cases of sulfuric acid and methanesulfonic acid. As the effect of relative humidity in salt particle formation varies many orders of magnitude depending on the acid and base molecules, neglecting hydration or using the same value for different systems may introduce remarkable inaccuracies in large-scale models.

Strain-Release-Controlled [4 + 2 + 1] Reaction of Cyclopropyl-Capped Diene-ynes/Diene-enes and Carbon Monoxide Catalyzed by Rhodium

Achieving transition-metal-catalyzed reactions of diene-ynes/diene-enes and carbon monoxide (CO) to deliver [4 + 2 + 1] cycloadducts, rather than the kinetically favored [2 + 2 + 1] products, is challenging. Here, we report that this can be solved by adding a cyclopropyl (CP) cap to the diene moiety of the original substrates. The resulting CP-capped diene-ynes/diene-enes can react with CO under Rh catalysis to give [4 + 2 + 1] cycloadducts exclusively without forming [2 + 2 + 1] products. This reaction has a broad scope and can be used to synthesize useful 5/7 bicycles with a CP moiety. Of the same importance, the CP moiety in the [4 + 2 + 1] cycloadducts can act as an intermediate group for further transformations so that other challenging bicyclic 5/7 and tricyclic 5/7/5, 5/7/6, and 5/7/7 skeletons, some of which are widely found in natural products, can be accessed. The mechanism of this [4 + 2 + 1] reaction has been investigated by quantum chemical calculations, and the role of the CP group in avoiding the possible side [2 + 2 + 1] reaction has been identified, showing that the [4 + 2 + 1] is controlled by releasing the ring strain in the methylenecyclopropyl (MCP) group (about 7 kcal/mol) in the CP-capped dienes.

Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties?

Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for β and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.

On the Reactivity of Phosphaalumenes towards C-C Multiple Bonds

Heterocycles containing group 13 and 15 elements such as borazines are an integral part of organic, biomedical and materials chemistry. Surprisingly, heterocycles containing P and Al are rare. We have now utilized phosphaalumenes in reactions with alkynes, alkenes and conjugated double bond systems. With sterically demanding alkynes 1,2-phosphaalumetes were afforded, whereas the reaction with HCCH or HCCSiMe₃ gave 1,4-phosphaaluminabarrelenes. Using styrene saturated 1,2-phosphaalumates were formed, which reacted further with additional styrene to give different regio-isomers of 1,4-aluminaphosphorinanes. Using ethylene, a 1,4-aluminaphosphorinane is obtained, while with 1,3-butadiene a bicyclic system containing an aluminacyclopentane and a phosphirane unit was synthesized. The experimental work is supported by theoretical studies to shed light on the mechanism governing the formation of these heterocycles.

Inverse Electron-Demanding Diels-Alder Reactions in the Chemical Synthesis of Prenylated Indole Alkaloids Containing a Bicycle[2.2.2]diazaoctane Moiety: A Theoretical Study

The Diels-Alder reaction is believed to be a key step in the biosynthesis of prenylated indole alkaloids containing a bicycle[2.2.2]diazaoctane moiety. Many chemical syntheses of bicyclic structures by Diels-Alder reactions have been reported, but the reaction mechanism remains underexplored. We have carried out DFT calculations on both acid- and base-promoted Diels-Alder reactions in these syntheses and reveal that the reactions occur through an inverse-electron demand mechanism. We hope that the new mechanism is helpful for the mechanistic understanding of the biosynthesis of this class of important natural products.

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted-core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller-Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol^-1 relative to the experiment. The two-point QZ/5Z and TZ/QZ CBS 1 / l max 3 extrapolation schemes are inferior to their inverse power counterpart ( 1 / l max + 1 / 2 4 ) by 1.3 and 4.3 kJ mol^-1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol^-1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double-ζ/triple-ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

Sparsity of the electron repulsion integral tensor using different localized virtual orbital representations in local second-order Møller-Plesset theory

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree-Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller-Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the "ragged list") and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the "fixed sparsity pattern") for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising.

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

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

A study of the thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) with atomic chlorine (Cl) in the absence and presence of water, using electronic structure methods

The thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) + atomic chlorine (Cl) were investigated in the absence and presence of a single water molecule, using electronic structure methods. Stationary points on each reaction surface were located using density functional theory (DFT) with the M06-2X functional with aug-cc-pVDZ (aVDZ) and aug-cc-pVTZ (aVTZ) basis sets. Then fixed point calculations were carried out using the UM06-2X/aVTZ optimised stationary point geometries, with aug-cc-pVnZ basis sets (n = T and Q), using the coupled cluster method [CCSD(T)], as well as the domain-based local pair natural orbitals coupled cluster [DLPNO-UCCSD(T)] approach. Four reaction channels are possible, formation of (A) CH₃SCH₂ + HCl, (B) CH₃S + CH₃Cl, (C) CH₃SCl + CH₃, and (C') CH₃S(Cl)CH₃. The results show that, in the absence of water, channels A and C' are the dominant channels. In the presence of water, the calculations show that the reaction mechanisms for A and C formation change significantly. Channel A occurs via submerged TSs and is expected to be rapid. Channel B occurs via TSs which present significant energy barriers indicating that this channel is not significant in the presence of water relative to CH₃SCH₂ + HCl and DMS·Cl adduct formation, as is the case in the absence of water. Channel C was not considered as it is endothermic in the absence of water. In the presence of water, pathways which proceed via (a) DMS·H₂O + Cl, (b) Cl·H₂O + DMS and (c) DMS·Cl + H₂O were considered. It was found that under tropospheric conditions, reactions via pathway (b) are of minor importance relative to those that proceed via pathways (a) and (c). This study has shown that water changes the mechanisms of the DMS + Cl reactions significantly but the presence of water is not expected to affect the overall reaction rate coefficient under atmospheric conditions as the DMS + Cl reaction has a rate coefficient at room temperature close to the collisional limit.

XPS and quantum chemical analysis of 4Me-BODIPY derivatives

The results of a X-ray photoelectron spectroscopy (XPS) and steady-state absorption spectroscopy study of the electronic structure, and cationic and excited states of a series of 1,3,5,7-tetramethyl-substituted BODIPYs (4Me,2R-BODIPYs) are presented. The experimental data were interpreted using high-level ab initio quantum chemical computations, including the algebraic diagrammatic construction method for the polarization propagator of the second order (ADC(2)), the outer-valence Green's function (OVGF) method, the density functional (DFT) approach, and the time-dependent DFT (TD-DFT) approach. Substitution effects on the XPS and absorption spectra were determined for 2,6-positions of 4Me,2R-BODIPY pyrrole nuclei (R = H, Br, Bu, benzyl). A very satisfactory performance of the DFT Koopmans theorem analogue was demonstrated with respect to the energy intervals between the electronic levels of 4Me,2R-BODIPY above 13 eV (BHHLYP functional) and the values of the HOMO-LUMO energy gap (ωB97X functional).

Bis(amidophenolato)phosphonium: Si-H Hydride Abstraction and Phosphorus-Ligand Cooperative Activation of C-C Multiple Bonds

The first bis(amidophenolato)phosphonium salts are prepared and fully characterized. The perfluorinated derivative represents the strongest monocationic phosphorus Lewis acid on the fluoride and hydride ion affinity scale isolable to date. This affinity enables new reactions, such as hydride abstraction from Et₃ SiH, the first phosphaalkoxylation of an alkyne or a phosphorus catalyzed intramolecular hydroarylation. All properties and reactions are scrutinized by theory and experiment. Substantial σ- and π-acidity provides the required affinity for substrate activation, while phosphorus-ligand cooperativity substantially enriches the reactivity portfolio of phosphonium ions.

Sub-system self-consistency in coupled cluster theory

In this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

Investigation of the Structures and Chemical Bonding of Mn2Ge6- and Mn2Ge7- Clusters via Anion Photoelectron Spectroscopy and Theoretical Calculations

We present joint anion photoelectron spectroscopy and theoretical studies for Mn₂Ge₆^- and Mn₂Ge₇^-. Experimental results show that Mn₂Ge₆^- and Mn₂Ge₇^- have vertical electron detachment energies of 2.58 ± 0.08 and 2.88 ± 0.08 eV, respectively. Both Mn₂Ge₆^- and Mn₂Ge₆ have C_s symmetric structures with a Mn atom attached to a pentagonal bipyramid MnGe₆. Both Mn₂Ge₇^- and Mn₂Ge₇ have C_2v symmetric structures, which can be considered as two Mn₂Ge₄ tetragonal bipyramids sharing a MnGeMn face. According to chemical bonding analyses, Mn₂Ge₆ could be considered as a (Mn^II)(MnGe₆^2-) complex. Theoretical calculations predict that the extension of Mn₂Ge₇ to Mn_2nGe_4n+3 may be able to produce a new type of Mn-doped germanium nanostructure.

Computational Study of Mechanisms and Tether Length Effects of Rh-Catalyzed [3+2] and [3+2+1] Reactions of Ene/Yne-Vinylcyclopropanes

DFT calculations have been applied to study the mechanisms of [3+2] and [3+2+1] reactions of ene/yne-vinylcyclopropanes (shorted as ene/yne-VCPs). The [3+2] reactions of ene/yne-VCPs start from C-C cleavage of cyclopropane (CP cleavage) to form six-membered rhodacycle, followed by alkene/alkyne insertion and reductive elimination. The [3+2+1] reactions have two competing pathways, one is the [3+2+1] pathway (CP cleavage, ene/yne insertion, CO insertion and reductive elimination) and the other is the [3+1+2] pathway (CP cleavage, CO insertion, ene/yne insertion and reductive elimination). The length of tether in substrates affects the ene/yne insertion steps in these cycloadditions, making some reactions fail or changing the reaction pathways. The reasons for these tether length effects are discussed.

Thermochemistry of per- and polyfluoroalkyl substances

The determination of gas phase thermochemical properties of per- and polyfluoroalkyl substances (PFAS) is central to understanding the long-range transport behavior of PFAS in the atmosphere. Prior gas-phase studies have reported the properties of perfluorinated sulfonic acid (PFOS) and perfluorinated octanoic acid (PFOA). Here, this study reports the gas phase enthalpies of formation of short- and long-chain PFAS and their precursor molecules determined using density functional theory (DFT) and ab initio approaches. Two density functionals, two ab initio methods and an empirical method were used to compute enthalpies of formation with the total atomization approach and an isogyric reaction. The performance of the computational methods employed in this work were validated against the experimental enthalpies of linear alkanoic acids and perfluoroalkanes. The gas-phase determinations will be useful for future studies of PFAS in the atmosphere, and the methodological choices will be helpful in the study of other PFAS.