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

Actinide-Actinide Bonding: Electron Delocalisation and σ-Aromaticity in the Tri-Thorium Cluster [{Th(η8 -C8 H8 )(μ-Cl)2 }3 K2 ]

The tri-thorium cluster [{Th(η8 -C8 H8 )(μ3 -Cl)2 }3 {K(THF)2 }2 ]∞ (Nature 2021, 598, 72-75) was reported to feature intriguing σ-aromatic bonding between the thorium atoms, a mode of metal-metal bonding unique in the actinide series. However, the presence of this bonding motif has since been challenged by others. Here, we computationally explore electron delocalisation in a molecular cluster fragment of [{Th(η8 -C8 H8 )(μ3 -Cl)2 }3 {K(THF)2 }2 ]∞ and examine its responses to an applied magnetic field using a variety of methods. We also discuss the importance of the choice of basis set for the Th atoms and issues regarding locating QTAIM bond critical points. When taken together, the computed data consistently suggest the presence of delocalised Th-Th bonding and Th3 σ-aromaticity.

Insights into the catalytic activity of boron-doped thiazoles in the Diels-Alder reaction

The role of boron-doped thiazoles as a Lewis acid catalyst in [4+2] cycloaddition reaction between 1,3-butadiene and acrolein has been addressed. Three different organic heterocycles were designed to study their catalytic activity. It has been observed that these heterocycles efficiently work as catalysts than the well-known Lewis acid BF3. All the reactions follow the normal electron demand process and are exothermic. Different conceptual DFT-based reactivity descriptors and electronic structure principles such as maximum hardness and minimum electrophilicity lend additional support to the feasibility of the reaction mechanism. The reaction force (RF), reaction electronic flux (REF), and its different components exhibit a detailed electronic activity throughout the reaction.

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.

Production of Dihydrogen Using Ammonia Borane as Reagent and Pyrazole as Catalyst

Theoretical chemistry (DLPNO-CCSD(T)/def2-TZVP//M06-2x/aug-cc-pVDZ) was used to design a system based on ammonia boranes catalyzed by pyrazoles with the aim of producing dihydrogen, nowadays of high interest as clean fuel. The reactivity of ammonia borane and cyclotriborazane were investigated, including catalytic activation through 1H-pyrazole, 4-methoxy-1H-pyrazole, and 4-nitro-1H-pyrazole. The results point toward a catalytic cycle by which, at the same time, ammonia borane can initially store and then, through catalysis, produce dihydrogen and amino borane. Subsequently, amino borane can trimerize to form cyclotriborazane that, in presence of the same catalyst, can also produce dihydrogen. This study proposes therefore a consistent progress in using environmentally sustainable (metal free) catalysts to efficiently extract dihydrogen from small B-N bonded molecules.

Phosphorescent Properties of Heteroleptic Ir(III) Complexes: Uncovering Their Emissive Species

In this contribution, we assess the computational machinery to calculate the phosphorescence properties of a large pool of heteroleptic [Ir(C^N)2(N^N)]+ complexes (where N^N is an ancillary ligand and C^N is a cyclometalating ligand) including their phosphorescent rates and their emission spectra. Efficient computational protocols are next proposed. Specifically, different flavors of DFT functionals were benchmarked against DLPNO-CCSD(T) for the phosphorescence energies. The transition density matrix and decomposition analysis of the emitting triplet excited state enable us to categorize the studied complexes into different cases, from predominant triplet ligand-centered (3LC) character to predominant charge-transfer (3CT) character, either of metal-to-ligand charge transfer (3MLCT), ligand-to-ligand charge transfer (3LLCT), or a combination of the two. We have also calculated the vibronically resolved phosphorescent spectra and rates. Ir(III) complexes with predominant 3CT character are characterized by less vibronically resolved bands as compared to those with predominant 3LC character. Furthermore, some of the complexes are characterized by close-lying triplet excited states so that the calculation of their phosphorescence properties poses additional challenges. In these scenarios, it is necessary to perform geometry optimizations of higher-lying triplet excited states (i.e., Tn). We demonstrate that in the latter scenarios all of the close-lying triplet species must be considered to recover the shape of the experimental emission spectra. The global analysis of computed emission energies, shape of the computed emission spectra, computed rates, etc. enable us to unambiguously pinpoint for the first time the triplet states involved in the emission process and to provide a general classification of Ir(III) complexes with regard to their phosphorescence properties.

On the Intermolecular Interactions in Thiophene-Cored Single-Stacking Junctions

There have been attempts, both experimental and based on density-functional theory (DFT) modeling, at understanding the factors that govern the electronic conductance behavior of single-stacking junctions formed by pi-conjugated materials in nanogaps. Here, a reliable description of relevant stacked configurations of some thiophene-cored systems is provided by means of high-level quantum chemical approaches. The minimal structures of these configurations, which are found using the dispersion-corrected DFT approach, are employed in calculations that apply the coupled cluster method with singles, doubles and perturbative triples [CCSD(T)] and extrapolations to the complete basis set (CBS) limit in order to reliably quantify the strength of intermolecular binding, while their physical origin is investigated using the DFT-based symmetry-adapted perturbation theory (SAPT) of intermolecular interactions. In particular, for symmetrized S-Tn dimers (where "S" and "T" denote a thiomethyl-containing anchor group and a thiophene segment comprising "n" units, respectively), the CCSD(T)/CBS interaction energies are found to increase linearly with n ≤ 6, and significant conformational differences between the flanking 2-thiophene group in S-T1 and S-T2 are described by the CCSD(T)/CBS and SAPT/CBS computations. These results are put into the context of previous work on charge transport properties of S-Tn and other types of supramolecular junctions.

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

Experimental and Theoretical Exploration of the Kinetics and Thermodynamics of the Nucleophile-Induced Fragmentation of Ylidenenorbornadiene Carboxylates

Ylidenenorbornadienes (YNDs), prepared by [4 + 2] cycloadditions between fulvenes and acetylene carboxylates, react with thiol nucleophiles to yield mixtures of four to eight diastereomers depending on the symmetry of the YND substrate. The mixtures of diastereomers fragment via a retro-[4 + 2] cycloaddition with a large variation in rate, with half-lives ranging from 16 to 11,000 min at 80 °C. The diastereomer-enriched samples of propane thiol adducts [YND-propanethiol (PTs)] were isolated and identified by nuclear Overhauser effect spectroscopy (NOESY) correlations. Simulated kinetics were used to extrapolate the rate constants of individual diastereomers from the observed rate data, and it correlated well with rate constants measured directly and from isolated diastereomer-enriched samples. The individual diastereomers of a model system fragment at differing rates with half-lives ranging from 5 to 44 min in CDCl3. Density functional theory calculations were performed to investigate the mechanism of fragmentation and support an asynchronous retro-[4 + 2] cycloaddition transition state. The computations generally correlated well with the observed free energies of activation for four diastereomers of the model system as a whole, within 2.6 kcal/mol. However, the observed order of the fragmentation rates across the set of diastereomers deviated from the computational results. YNDs display wide variability in the rate of fragmentation, dependent on the stereoelectronics of the ylidene substituents. A Hammett study showed that the electron-rich aromatic rings attached to the ylidene bridge increase the fragmentation rate, while electron-deficient systems slow fragmentation rates.

Accurate Evaluation of Dispersion Energies at Coupled Cluster Level to Understand the Substituent Effects in Am(III) and Eu(III) Complexes

The effect of cyclic and aromatic substituents on the complexation behavior of phosphine oxide ligands with Am(III) and Eu(III) was investigated at density functional theory (DFT) and domain-based local pair natural orbital coupled-cluster (DLPNO-CC) levels. Combining DFT with accurate coupled cluster methods, we have evaluated the dispersion energy contributions to the complexation energies for trivalent Am and Eu complexes for the first time. Irrespective of the nature of substituents on the P atom, the electronic structure of the P═O group remains identical in all of the ligands. The study reveals the importance of dispersion interactions during complexation and is estimated to be more significant for Am(III) than for Eu(III) complexes.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers

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

Rh-Catalyzed [4 + 1] Reaction of Cyclopropyl-Capped Dienes (but not Common Dienes) and Carbon Monoxide: Reaction Development and Mechanistic Study

Transition-metal-catalyzed [4 + 1] reaction of dienes and carbon monoxide (CO) is the most straightforward and easily envisioned cyclization for the synthesis of five-membered carbocycles, which are ubiquitously found in natural products and functional molecules. Unfortunately, no test of this reaction was reported, and consequently, chemists do not know whether such kind of reaction works or not. Herein, we report that the [4 + 1] reaction of common dienes and CO cannot work, at least under the catalysis of [Rh(cod)Cl]2. However, using cyclopropyl-capped dienes (also named allylidenecyclopropanes) as substrates, the corresponding [4 + 1] reaction with CO proceeds smoothly in the presence of [Rh(cod)Cl]2. This [4 + 1] reaction, with a broad scope, provides efficient access to five-membered carbocyclic compounds of spiro[2.4]hept-6-en-4-ones. The [4 + 1] cycloadducts can be further transformed into other molecules by using the unique chemistry of cyclopropyl groups present in these molecules. The mechanism of this [4 + 1] reaction has been investigated by quantum chemical calculations, uncovering that cyclopropyl-capped dienes are strained dienes and the oxidative cyclization step in the [4 + 1] catalytic cycle can release this (angular) strain both kinetically and thermodynamically. The strain release in this step then propagates to all followed CO coordination/CO insertion/reductive elimination steps in the [4 + 1] catalytic cycle, helping the realization of this cycloaddition reaction. In contrast, common dienes (including cyclobutyl-capped dienes) do not have such advantages and their [4 + 1] reaction suffers from energy penalty in all steps involved in the [4 + 1] catalytic cycle. The reactivity of ene-allenes for the [4 + 1] reaction with CO is also discussed.

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.

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger-Weber potentials. The model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. The model enables simulation of early stage Al(OH)₃ nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

Metal-free catalytic hydroboration of imine with pinacolborane using a pincer-type phosphorus compound: mechanistic insight and improvement of the reaction

A mechanistic study of the catalytic hydroboration of imine using a pincer-type phosphorus compound 1NP was performed through the combination of DFT and DLPNO-CCSD(T) calculations. The reaction proceeds through a phosphorus-ligand cooperative catalytic cycle, where the phosphorus center and triamide ligand work in a synergistic manner. First, the pinB-H bond activation by 1NP occurs through the cooperative functions of the phosphorus center and the triamide ligand, leading to a phosphorus-hydride intermediate 2NP. This is the rate-determining step, with the Gibbs energy barrier and Gibbs reaction energy of 25.3 and -17.0 kcal mol^-1, respectively. Subsequently, the hydroboration of phenylmethanimine takes place through a concerted transition state through the cooperative function of the phosphorus center and the triamide ligand. It leads to the final hydroborated product 4 with the regeneration of 1NP. Our computational results reveal that the experimentally isolated intermediate 3NP is a resting state of the reaction. It is formed through the B-N bond activation of 4 by 1NP, rather than via the insertion of the CN double bond of phenylmethanimine into the P-H bond of 2NP. However, this side reaction can be suppressed by utilizing a planar phosphorus compound AcrDipp-1NP as the catalyst, which features steric-demanding substituents on the chelated N atom of the ligand.

Low-energy electron interaction with 2-(trifluoromethyl)acrylic acid, a potential component for EUVL resist material

Motivated by the current introduction of extreme ultraviolet lithography (EUVL) into chip manufacturing processes, and thus the transition to electron-induced chemistry within the respective resist materials, we have studied low energy electron-induced fragmentation of 2-(trifluoromethyl)acrylic acid (TFMAA). This compound is chosen as a potential resist component, whereby fluorination enhances the EUV adsorption and may at the same time promote electron-induced dissociation. Dissociative ionization and dissociative electron attachment are studied, and to aid the interpretation of the observed fragmentation channels, the respective threshold values are calculated at the DFT and coupled cluster level of theory. Not surprisingly, we find significantly more extensive fragmentation in DI than in DEA and in fact, the only significant DEA fragmentation channel is the cleavage of HF from the parent molecule upon electron attachment. Rearrangement and new bond formation are substantial in DI and are, in fact, similar to DEA, mainly associated with HF formation. The observed fragmentation reactions are discussed in relation to the underlying reactions and potential implications for the suitability of TFMAA as a component of EUVL resist materials.

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations

The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method has been proven to provide accurate single-point energies at a fraction of the cost of canonical CCSD(T) calculations. However, the desired "chemical accuracy" can only be obtained with a large PNO space and extended basis set. We present a simple yet accurate and efficient correction scheme based on a perturbative approach. Here, in addition to DLPNO-CCSD(T) energy, one calculates DLPNO-MP2 correlation energy with the same settings as in the preceding coupled-cluster calculation. In the next step, the canonical MP2 correlation energy is obtained in the same orbital basis. This can be efficiently performed for essentially all molecule sizes accessible with the DLPNO-CCSD(T) method. By taking the difference between the canonical MP2 and DLPNO-MP2 energies, we obtain a correction term that can be added to the DLPNO-CCSD(T) correlation energy. This way, one can obtain the total correlation energy close to the limit of the complete PNO space (cPNO). The presented approach allows us to significantly increase the accuracy of the DLPNO-CCSD(T) method for both closed- and open-shell systems. The latter are known to be especially challenging for locally correlated methods. Unlike the previously developed PNO extrapolation procedure by Altun, Neese, and Bistoni ( J. Chem. Theory Comput. 2020, 16, 6142-6149), this strategy allows us to get the DLPNO-CCSD(T) correlation energy at the cPNO limit in a cost-efficient way, resulting in a minimal overall increase in calculation time as compared to the uncorrected method.

Coupled cluster theory on modern heterogeneous supercomputers

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

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.