Dispersion-Corrected Mean-Field Electronic Structure Methods

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

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

Translational eigenstates of He@C60 from four-dimensional ab initio potential energy surfaces interpolated using Gaussian process regression

We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two ab initio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.

Effective Descriptor for Screening Single-Molecule Conductance Switches

The adsorption-type molecular switch exhibits bistable states with an equivalently long lifetime at the organic/inorganic interface, promising reliable switching behavior and superior assembly ability in the electronic circuits at the molecular scale. However, the number of reported adsorption-type molecular switches is currently less than 10, and exploring these molecular switches poses a formidable challenge due to the intricate interplay occurring at the interface. To address this challenge, we have developed a model enabling the identification of diverse molecular switches on metal surfaces based on easily accessible physical characteristics. These characteristics primarily include the metal valency electron concentration, the work function of metal surfaces, and the electronegativity difference of molecules. Using this model, we identified 56 new molecular switches. Employing the gradient descent algorithm and statistical linear discriminant analysis, we constructed an explicit descriptor that establishes a relationship between the interfacial structure and chemical environment and the stability of molecular switches. The model's accuracy was validated through density functional theory calculations, achieving a 90% accuracy for aromatic molecular switches. The conductive switching behaviors were further confirmed by nonequilibrium Green's function transport calculations.

Electronic Structure of Metalloporphenes, Antiaromatic Analogues of Graphene

Zinc porphene is a two-dimensional material made of fully fused zinc porphyrins in a tetragonal lattice. It has a fully conjugated π-system, making it similar to graphene. Zinc porphene has recently been synthesized, and a combination of rough conductivity measurements and infrared and Raman spectroscopies all suggested that it is a semiconductor (Magnera, T.F. et al. Porphene and Porphite as Porphyrin Analogs of Graphene and Graphite, Nat. Commun.2023, 14, 6308). This is in contrast with all previous predictions of its electronic structure, which indicated metallic conductivity. We show that the gap-opening in zinc porphene is caused by a Peierls distortion of its unit cell from square to rectangular, thus giving the first account of its electronic structure in agreement with the experiment. Accounting for this distortion requires proper treatment of electron delocalization, which can be done using hybrid functionals with a substantial amount of exact exchange. Such a functional, PBE38, is then applied to predict the properties of many first transition row metalloporphenes, some of which have already been prepared. We find that changing the metal strongly affects the electronic structure of metalloporphenes, resulting in a rich variety of both metallic conductors and semiconductors, which may be of great interest to molecular electronics and spintronics. Properties of these materials are mostly governed by the extent of the Peierls distortion and the number of electrons in their π-system, analogous to changes in aromaticity observed in cyclic conjugated molecules upon oxidation or reduction. These results give an account of how the concept of antiaromaticity can be extended to periodic systems.

Closed-shell d10-d10 mechanochromic [AuPh(CNPh)]n complex: quantum chemistry electronic and optical properties

The electronic structure, spectroscopic properties, and solid state chemistry of monomer and dimers of [AuPh(CNPh)] complex were studied at post-Hartree-Fock (MP2, SCS-MP2, and CC2) and density functional theory levels. The absorption spectra of these complexes were calculated using single excitation time-dependent (TD) methods at DFT, CC2, and SCS-CC2 levels. The influences of the bulk are accounted for at the PBE-D3 level, incorporating dispersion effects. The calculated values agree with the experimental range, where absorption and emission energies reproduce experimental trends with large Stokes shifts. The aurophilic interaction is identified as a key factor influencing the spectroscopic and structural properties of these complexes. The intermetallic interactions were found as the main factor responsible for MMCT electronic transitions in the models studied.

Mid-infrared spectroscopy of 1-cyanonaphthalene cation for astrochemical consideration

We present the low temperature gas-phase vibrational spectrum of ionised 1-cyanonaphthalene (1-CNN+) in the mid-infrared region. Experimentally, 1-CNN+ ions are cooled below 10 K in a cryogenic ion trapping apparatus, tagged with He atoms and probed with tuneable radiation. Quantum-chemical calculations are carried out at a density functional theory level. The spectrum is dominated by the CN-stretch at 4.516 μm, with weaker CH modes near 3.2 μm.

Excited state hydrogen or proton transfer pathways in microsolvated n-cyanoindole fluorescent probes

The sensitivity of the fluorescence properties of n-cyanoindole (n-CNI) fluorescent probes to the microenvironment makes them potential reporters of protein conformation and hydration. The fluorescence intensity of 5-CNI, 6-CNI, and 7-CNI is quenched when exposed to water solvent whereas substitution on position 4 of indoles dramatically increases it. A potential mechanism for this sensitivity to water may be similar to that found in indole. The fluorescence of indole is found to be quenched when interacting with water and ammonia solvent molecules via radiationless decay through an S1 (πσ*)/S0 conical intersection caused by excited state proton or hydrogen transfer to the solvent molecules. In this study we examine this fluorescence quenching mechanism along the N-H bond stretch of n-CNI probes using water cluster models and quantum mechanical calculations of the excited states. We find that n-CNI-(H2O)1-2 clusters form cyclic or non-cyclic structures via hydrogen bonds which lead to different photochemical reaction paths that can potentially quench the fluorescence by undergoing internal conversion from S1 to S0. However, the existence of a high energy barrier along the potential energy surface of the S1 state in most cases prevents this from occurring. We show that substitution on position 4 leads to the highest energy barrier that prevents the fluorophore from accessing these non-radiative channels, in agreement with its high intensity. We also find that the energy barrier in the S1 state of non-cyclic 5-CNI-(H2O)1-2 excited complexes decreases as the number of water molecules increases, which suggests great sensitivity of the fluorescence quenching on the aqueous environment.

Electronic Properties of Graphene Nano-Parallelograms: A Thermally Assisted Occupation DFT Computational Study

In this computational study, we investigate the electronic properties of zigzag graphene nano-parallelograms (GNPs), which are parallelogram-shaped graphene nanoribbons of various widths and lengths, using thermally assisted occupation density functional theory (TAO-DFT). Our calculations revealed a monotonic decrease in the singlet-triplet energy gap as the GNP length increased. The GNPs possessed singlet ground states for all the cases examined. With the increase of GNP length, the vertical ionization potential and fundamental gap decreased monotonically, while the vertical electron affinity increased monotonically. Some of the GNPs studied were found to possess fundamental gaps in the range of 1-3 eV, lying in the ideal region relevant to solar energy applications. Besides, as the GNP length increased, the symmetrized von Neumann entropy increased monotonically, denoting an increase in the degree of the multi-reference character associated with the ground state GNPs. The occupation numbers and real-space representation of active orbitals indicated that there was a transition from the nonradical nature of the shorter GNPs to the increasing polyradical nature of the longer GNPs. In addition, the edge/corner localization of the active orbitals was found for the wider and longer GNPs.

Balance between Physical Interpretability and Energetic Predictability in Widely Used Dispersion-Corrected Density Functionals

We assess the performance of different dispersion models for several popular density functionals across a diverse set of noncovalent systems, ranging from the benzene dimer to molecular crystals. By analyzing the interaction energies and their individual components, we demonstrate that there exists variability across different systems for empirical dispersion models, which are calibrated for reproducing the interaction energies of specific systems. Thus, parameter fitting may undermine the underlying physics, as dispersion models rely on error compensation among the different components of the interaction energy. Energy decomposition analyses reveal that, the accuracy of revPBE-D3 for some aqueous systems originates from significant compensation between dispersion and charge transfer energies. However, revPBE-D3 is less accurate in describing systems where error compensation is incomplete, such as the benzene dimer. Such cases highlight the propensity for unpredictable behavior in various dispersion-corrected density functionals across a wide range of molecular systems, akin to the behavior of force fields. On the other hand, we find that SCAN-rVV10, a targeted-dispersion approach, affords significant reductions in errors associated with the lattice energies of molecular crystals, while it has limited accuracy in reproducing structural properties. Given the ubiquitous nature of noncovalent interactions and the key role of density functional theory in computational sciences, the future development of dispersion models should prioritize the faithful description of the dispersion energy, a shift that promises greater accuracy in capturing the underlying physics across diverse molecular and extended systems.

Regulating the Content of Donor Unit in Donor-Acceptor Covalent Triazine Frameworks for Promoting Photocatalytic H2 Production

Using their own triazine groups as natural receptors, the introduction of various donor units to construct donor-receptor configuration in covalent triazine frameworks (CTFs) has been shown to be an effective strategy to improve photocatalytic activity. In this work, the effect of donor unit content (D-content) on the photoelectric properties and photocatalytic activity of CTFs was thoroughly investigated. Four analogous CTFs with different D-content have been rationally designed and synthesized, in which the bithiophene (Btp) as the donor unit and triazine as the acceptor unit. And CTF-Btp with the highest D-content showed the best photocatalytic activity. The experimental and theoretical results indicated this improvement is attributed to stronger visible light absorption capacity and higher photoinduced charge carrier separation efficiency. This study elucidates the relationship between the structural features of CTFs with varying D-content and their photocatalytic activity, offering a promising strategy for developing efficient photocatalysts.

Chemiluminescence of a Firefly Luciferin Analogue Reveals that Formation of the Key Intermediate Responsible for Excited State Generation Occurs on a Fully Concerted Step

The chemiluminescence (CL) reaction of eight different 2-(4-hydroxyphenyl)-4,5-dihydrothiazole-4-carboxylate esters with an organic superbase and oxygen was investigated through a kinetic and computational study. These esters are all analogues to the luciferin substrate involved in efficient firefly bioluminescence. The kinetic data obtained from CL emission and light absorption assays were used in the context of linear free energy relationships (LFER); we obtained the Hammett reaction constant ρ = +1.62 ± 0.09 and the Brønsted constant βlg = -0.39 ± 0.04. These observations from LFER, together with activation parameters obtained from Arrhenius plots, suggest that the formation of the high-energy intermediate (HEI) 1,2-dioxetanone occurs via a concerted mechanism during the rate-determining step of the reaction. Calculations performed using density functional theory support a late transition state for HEI formation within the reaction mechanism pathway, which was described considering geometric parameters, Wiberg bond indices from natural bond order analysis, and the atomic charges derived from the electrostatic potential.

Influence of Steric and Dispersion Interactions on the Thermochemistry of Crowded (Fluoro)alkyl Compounds

ConspectusAlkanes play a pivotal role in industrial, environmental, and biological processes. They are characterized by their carbon-carbon single-bond structure, remarkable stability, and conformational diversity. Fluorination of such compounds imparts unique physicochemical properties that often enhance pharmacokinetic profiles, metabolic stability, and receptor interactions while keeping beneficial properties. However, such per- and polyfluoroalkyl substances (PFAS) show a persistent presence in the environment and potential adverse health effects, which propelled them to the forefront of global environmental and health discussions. Alkyl compounds are also prototypical for stereoelectronic (SE) effects that are widely applied in chemistry. Substituents are typically described as electron-density-donating/withdrawing and/or responsible for sterically interacting with reagents or strategic groups in the molecule. That alkane branching can result in higher stability compared to less-branched isomers has been investigated in detail also by testing quantum chemical methods, in particular density functional theory (DFT). Alkane branching results in spatially compact structures with close intramolecular contacts so that at a specific size the detailed balance of attractive London dispersion and covalent versus repulsive Pauli exchange interactions shifts to new, chemically unfragile situations. This may lead to dissociation at room temperature and opens the central question: what is the smallest crowed alkane that cannot be made synthetically? In this Account, we try to shed light on the interplay among the various (free) energy components for crowded (fluoro)alkane dissociation. In this context, homolytic cleavage of the central C-C bond in a series of model alkanes of increasing size with tert-butyl (tBu), adamantyl (Ad), and [1.1.1]propellanyl (Prop) substituents is investigated. Reference energies are calculated at the PNO-LCCSD(T)-F12b level and used to benchmark the performance of contemporary DFT functionals. In line with previous conclusions, the application of dispersion corrections to density functionals is mandatory. For crowed structures, the accurate description of the midrange correlation effects, specifically repulsive van der Waals interactions, is crucial, and we observed that the density-dependent VV10 correction is superior to D4 in this context, although the asymptotic region is better described by the latter. The best available dispersion-inclusive functionals show systematic and reasonably small residual errors and can be safely applied to large systems (>100 atoms), for which coupled cluster methods with large basis sets are not computationally feasible anymore. For qualitatively correct predictions of synthetic accessibility under equilibrium conditions (free energy), the inclusion of thermostatistical (entropy) contributions is also essential. According to our results, tetra-tert-butylmethane (C17tBu) is the largest and most crowded system with a positive dissociation free energy and should be synthesizable. The difference between hydrogenated and perfluorinated systems originates from the increase in the steric repulsion of spatially close substituents, which is not compensated to the same extent by attractive orbital and dispersion interactions. A sometimes-assumed similar steric demand for fluorine and hydrogen atoms is not corroborated by our investigations on crowded systems. Perfluorination is found to substantially decrease thermal stability, rendering perfluorinated hexamethylethane (C8tBuF) the last potentially stable representative.

Novel conjugated microporous polymers for efficient tetracycline adsorption: insights from theoretical investigations

This paper presents a detailed theoretical understanding of the noncovalent interactions between antibiotics tetracycline and conjugated microporous polymer (CMP), which is important to understand the recent experimental finding of efficient removal of antibiotics by CMP materials. We show that the co-work of π-π and H-π interactions determines the final equilibrium structures, when a tetracycline molecule spontaneously adsorbs to the surface or within the pores of the CMP network at physisorption distances. The binding energies for tetracycline/CMP systems are calculated to be -0.31 ∼ -1.15 eV, demonstrating the reliability of the adsorption. The electronic structures of CMP nanostructures remain basically undamaged upon the tetracycline adsorption. The replacement of benzothiadiazole unit with S and N heteroatoms to the phenyl moiety in the linker effectively enhanced the molecular polarity of CMP molecule and increases the interaction area between tetracycline and CMP network, consequently enhancing the average binding energies notably. Our calculations provide useful theoretical guidance for design of novel carbon-based porous adsorbents with good adsorption performance to remove residual tetracycline and other antibiotics in water.

Theoretical investigation of asphaltene molecules in crude oil viscoelasticity enhancement

Understanding the mechanisms of viscosity enhancement in crude oil phases is crucial for optimizing extraction and transportation processes. The enhanced viscosity mechanism of crude oil phase can be attributed to the intricate intermolecular interactions between asphaltene molecules. However, the molecular mechanism of the viscosification of asphaltene molecules in crude oil is not yet to be fully understood. In this work, molecular dynamics simulations were employed to investigate the dynamic behavior and viscosification mechanism of asphaltene molecules in complex oil phases. Research suggests that the neutral surface of asphaltenes features abundant positive and negative electrostatic potential regions, facilitating complementary pairing between these areas. This significantly augments electrostatic interactions among asphaltene molecules. Besides, the expansive nonpolar expanse on the normal asphaltene surface facilitates interactions between asphaltenes and crude oil molecules. This leads the crude oil viscosity of the system containing normal asphaltene is higher than that of the system containing acidic asphaltene under the same mass fraction (382 μ Pa·s for AAsp and 416 μ Pa·s for NAsp). This work provides insight into the viscosity enhancement mechanisms in crude oil phases and is helpful in improving the efficiency of crude oil extraction and transportation.

Design and Exploration of Benzene Like Azobis Triazoles for Long-range Push-Pull Photo-Switching Attributes

This research paper presents a comprehensive study on the design and photovoltaic parameters of azobenzene type 24 photo switches (PSs) of triazole by density functional theory (DFT). The focus was on investigating how to create a long-range push-pull effect of different substituents on the PS properties for their application in photovoltaics by further substituent decoration. Their range of values for the maximum wavelength (λmax) ranged 315-556 nm while their HOMO-LUMO energies (Egaps) were 0.57-6.35eV. The stability of the PS was evaluated by measuring hardness (η) and softness (σ) values. Additionally, photovoltaic parameters such as open-circuit voltage (Voc), short-circuit current density (Jsc), fill factor (FF), and maximum power (Pmax) were calculated to assess the performance of the PS as photovoltaic materials. The results revealed that PSs 6 exhibited promising photovoltaic parameters to include Voc values ranging from 0.4-1.63eV, FF values ranging from 0.5438-0.929, Jsc values ranging from 19.27-50.75 mA/cm2, and Pmax values ranging from 14.72-75.91W. This indicates its potential as an efficient light-harvesting material for photovoltaic applications. Moreover, this study presents a pioneering investigation on the correlation between rotational velocity (R) and Mayer bond index (MBI) for the first time. The findings revealed a significant correlation between R and MBI, providing valuable insights into the structural dynamics of the PS. This novel finding opens up new avenues for understanding the structural dynamics of PS and their potential applications in various fields, including photovoltaics. The study provides valuable insights into the structure-property relationships of azobenzene-based PS and their suitability for photovoltaic devices. Further investigations are warranted to optimize the design of the PS, enhance their photovoltaic performance, and explore the underlying mechanisms of the correlation between R and MBIs.

Contributions of London Dispersion Forces to Solution-Phase Association Processes

ConspectusDespite their ubiquity and early discovery, London dispersion forces are often overlooked. This is due, in part, to the difficulty in assessing their contributions to molecular and polymeric structure, stability, properties, and reactivities. However, recent advances in modeling have revealed that dispersion interactions play an important role in many important chemical and biological processes. Experimental confirmation of their impact in solution has been challenging, leading to controversies about their relative importance.In the course of studying noncovalent interactions using molecular devices, our understanding and appreciation for the importance of dispersion interactions have evolved. This Account follows this intellectual journey by using examples from the literature. The goals are twofold: to describe recent advances in understanding the interaction and to provide guidance to researchers studying weak noncovalent interactions. However, first, the experimental methods for measuring the effects of dispersion interactions and the strategies for isolating their influence are described. These include the design of molecular devices to measure these weak noncovalent interactions and the strategies to disentangle the solvation, solvophobic, and dispersion components of the resulting equilibria.The literature examples are organized around five fundamental questions. (1) Do dispersion interactions have a measurable effect on solution equilibria? (2) To what extent do solvents attenuate or compensate for dispersion interactions? (3) To what extent do the solvation and solvophobic terms influence the dispersion equilibria? (4) Can we predict whether a system will form attractive dispersion or repulsive steric interactions? (5) Can the dispersion term be isolated and interrogated? We were often surprised by the answers to these questions. In each case, we describe how the systems were designed to address these questions and discuss possible interpretations of the results.While dispersion interactions in solution were weak (usually <1 kcal/mol), their influence on complexation and conformational equilibria can be observed and measured. This underscores the significance of these interactions in molecular recognition, coordination chemistry, reaction design, and catalysis. The solvent components of the dispersion equilibria can also be significant. Therefore, the isolation of the dispersion contributions from the solvation and solvophobic effects represents an ongoing challenge. The experimental studies also provide important benchmarks and offer valuable insights to help refine the next generation of computational solvent models.

Second quantization of many-body dispersion interactions for chemical and biological systems

The many-body dispersion (MBD) framework is a successful approach for modeling the long-range electronic correlation energy and optical response of systems with thousands of atoms. Inspired by field theory, here we develop a second-quantized MBD formalism (SQ-MBD) that recasts a system of atomic quantum Drude oscillators in a Fock-space representation. SQ-MBD provides: (i) tools for projecting observables (interaction energy, transition multipoles, polarizability tensors) on coarse-grained representations of the atomistic system ranging from single atoms to large structural motifs, (ii) a quantum-information framework to analyze correlations and (non)separability among fragments in a given molecular complex, and (iii) a path toward the applicability of the MBD framework to molecular complexes with even larger number of atoms. The SQ-MBD approach offers conceptual insights into quantum fluctuations in molecular systems and enables direct coupling of collective plasmon-like MBD degrees of freedom with arbitrary environments, providing a tractable computational framework to treat dispersion interactions and polarization response in intricate systems.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Combined Spectroscopic, Thermodynamic, and Theoretical Approach for Detecting and Quantifying Hydrogen Bonding and Dispersion Interaction in Ionic Liquids

ConspectusIonic liquids (ILs) are attracting increasing interest in science and engineering due to their unique properties that can be tailored for specific applications. Clearly, a better understanding of their behavior on the microscopic scale will help to elucidate macroscopic fluid phenomena and thereby promote potential applications. The advantageous properties of these innovative fluids arise from the delicate balance of Coulomb interactions, hydrogen bonding, and dispersion forces. The development of these properties requires a fundamental understanding of the strength, location, and direction of the different types of interactions and their contribution to the overall phase behavior. Contrary to expectations, hydrogen bonding and dispersion interactions have a significant influence on the structure, dynamics, and phase behavior of ILs.The synergy between experimental and theoretical methods has now advanced to a stage where hydrogen bonds and dispersion effects as well as the competition between the two can be studied in detail. In this account, we demonstrate that a suitable combination of spectroscopic, thermodynamic, and theoretical methods enables the detection, dissection, and quantification of noncovalent interactions, even in complex systems such as ionic liquids. This approach encompasses far-infrared vibrational spectroscopy (FIR), various thermodynamic methods for determining enthalpies of vaporization, and quantum chemical techniques that allow us to switch dispersion contributions on or off when calculating the energies and spectroscopic properties of clusters.We briefly discuss these experimental and theoretical methods, before providing various examples illustrating how the mélange of Coulomb interaction, hydrogen bonds, and dispersion forces can be analyzed, and their individual contributions quantified. First, we demonstrated that both hydrogen bonding and dispersion interactions are manifested in the FIR spectra and can be quantified by observed shifts of characteristic spectral signatures. Through the selection of suitable protic ionic liquids (PILs) featuring anions with varying interaction strengths and alkyl chain lengths, we were able to demonstrate that dispersion interactions can compete with hydrogen bonding. The resultant transition enthalpy serves as a measure of the dispersion interaction. Contrary to expectations, PILs possess lower enthalpies of vaporization compared with aprotic ILs (AILs). The reason for this is simple: In protic ILs, ion pairs carry both the hydrogen bond and attractive dispersion between the cation and anion into the gas phase. By utilizing a well-curated set of protic ILs and molecular analogues, we successfully disentangled Coulomb interaction, hydrogen bonding, and dispersion interaction through purely thermodynamic methods.

Improved Description of Intra- and Intermolecular Interactions through Dispersion-Corrected Second-Order Møller-Plesset Perturbation Theory

ConspectusThe quantum chemical modeling of organic crystals and other molecular condensed-phase problems requires computationally affordable electronic structure methods which can simultaneously describe intramolecular conformational energies and intermolecular interactions accurately. To achieve this, we have developed a spin-component-scaled, dispersion-corrected second-order Møller-Plesset perturbation theory (SCS-MP2D) model. SCS-MP2D augments canonical MP2 with a dispersion correction which removes the uncoupled Hartree-Fock dispersion energy present in canonical MP2 and replaces it with a more reliable coupled Kohn-Sham treatment, all evaluated within the framework of Grimme's D3 dispersion model. The spin-component scaling is then used to improve the description of the residual (nondispersion) portion of the correlation energy.The SCS-MP2D model improves upon earlier corrected MP2 models in a few ways. Compared to the highly successful dispersion-corrected MP2C model, which is based solely on intermolecular perturbation theory, the SCS-MP2D dispersion correction improves the description of both inter- and intramolecular interactions. The dispersion correction can also be evaluated with trivial computational cost, and nuclear analytic gradients are computed readily to enable geometry optimizations. In contrast to earlier spin-component scaling MP2 models, the optimal spin-component scaling coefficients are only mildly sensitive to the choice of training data, and a single global parametrization of the model can describe both thermochemistry and noncovalent interactions.The resulting dispersion-corrected, spin-component-scaled MP2 (SCS-MP2D) model predicts conformational energies and intermolecular interactions with accuracy comparable to or better than those of many range-separated and double-hybrid density functionals, as is demonstrated on a variety of benchmark tests. Among the functionals considered here, only the revDSD-PBEP86-D3(BJ) functional gives consistently smaller errors in benchmark tests. The results presented also hint that further improvements of SCS-MP2D may be possible through a more robust fitting procedure for the seven empirical parameters.To demonstrate the performance of SCS-MP2D further, several applications to molecular crystal problems are presented. The three chosen examples all represent cases where density-driven delocalization error causes GGA or hybrid density functionals to artificially stabilize crystals exhibiting more extended π-conjugation. Our pragmatic strategy addresses the delocalization error by combining a periodic density functional theory (DFT) treatment of the infinite lattice with intramolecular/conformational energy corrections computed with SCS-MP2D. For the anticancer drug axitinib, applying the SCS-MP2D conformational energy correction produces crystal polymorph stabilities that are consistent with experiment, in contrast to earlier studies. For the crystal structure prediction of the ROY molecule, so named for its colorful red, orange, and yellow crystals, this approach leads to the first plausible crystal energy landscape, and it reveals that the lowest-energy polymorphs have already been found experimentally. Finally, in the context of photomechanical crystals, which transform light into mechanical work, these techniques are used to predict the structural transformations and extract design principles for maximizing the work performed.

The pseudotorsional space of RNA

The characterization of the conformational landscape of the RNA backbone is rather complex due to the ability of RNA to assume a large variety of conformations. These backbone conformations can be depicted by pseudotorsional angles linking RNA backbone atoms, from which Ramachandran-like plots can be built. We explore here different definitions of these pseudotorsional angles, finding that the most accurate ones are the traditional η (eta) and θ (theta) angles, which represent the relative position of RNA backbone atoms P and C4'. We explore the distribution of η - θ in known experimental structures, comparing the pseudotorsional space generated with structures determined exclusively by one experimental technique. We found that the complete picture only appears when combining data from different sources. The maps provide a quite comprehensive representation of the RNA accessible space, which can be used in RNA-structural predictions. Finally, our results highlight that protein interactions lead to significant changes in the population of the η - θ space, pointing toward the role of induced-fit mechanisms in protein-RNA recognition.

Do Optimally Tuned Range-Separated Hybrid Functionals Require a Reparametrization of the Dispersion Correction? It Depends

For ground- and excited-state studies of large molecules, it is the state of the art to combine (time-dependent) DFT with dispersion-corrected range-separated hybrid functionals (RSHs), which ensures an asymptotically correct description of exchange effects and London dispersion. Specifically for studying excited states, it is common practice to tune the range-separation parameter ω (optimal tuning), which can further improve the accuracy. However, since optimal tuning essentially changes the functional, it is unclear if and how much the parameters used for the dispersion correction depend on the chosen ω value. To answer this question, we explore this interdependency by refitting the DFT-D4 dispersion model for six established RSHs over a wide range of ω values (0.05-0.45 a0-1) using a set of noncovalently bound molecular complexes. The results reveal some surprising differences among the investigated functionals: While PBE-based RSHs and ωB97M-D4 generally exhibit a weak interdependency and robust performance over a wide range of ω values, B88-based RSHs, specifically LC-BLYP, are strongly affected. For these, even a minor reduction of ω from the default value manifests in strong systematic overbinding and poor performance in the typical range of optimally tuned ω values. Finally, we discuss strategies to mitigate these issues and reflect the results in the context of the employed D4 parameter optimization algorithm and fit set, outlining strategies for future improvements.

Interactions between Metals and Eudistomins of Ascidian Origin: A Computational Study

Ascidians are marine animals that adopt unusual techniques to deter predation. The three main methods are sequestration of unusual metals, high concentrations of sulfuric acid/sulfate ions in tunicate cells, and the presence of eudistomins. In this study, we hypothesize that ascidians sequester metals in their sulfate form, and the complexation of eudistomins with the metals could liberate the sulfate ion. Three representative metal aqua ions were chosen, viz., vanadyl, uranyl, and thorium ions, as well as four simple eudistomins which act as bidentate ligands, viz., eudistomin-W, debromoeudistomin-K, eudistomidin-C, and eudistomidin-B. By designing 7 model reactions, we tested our hypothesis using density functional theory (DFT) methods PBE-D3, BLYP, and B3LYP. The ΔG values of the model reactions provide strong support for our hypothesis. To verify the hypothesis further, we calculated the metal-eudistomin interactions with Be, Zn, and Pb. Based on our results, we suggest that ascidians may not prefer any particular metal. In addition, despite using different DFT functionals, we have observed similar ΔG values for each case. With our work, we have successfully used computational tools in our attempt to understand the unique behavior of ascidians.

Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments

The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

What do we know about bifunctional cage-like T8 silsesquioxanes? Theory versus lab routine

In this article, we explore theoretical validations of experimental findings pertaining to the classical corner-capping reactions of a commercially available heptaisobutyltrisilanol cage to mono-substituted phenylhepta(isobutyl)-POSS cages. Additionally, the process of opening a fully condensed cage is tracked to assess the possibility of isolating and separating the resulting isomers. The corner-capping reactions of potential silanotriols, both as monomers and dimers, and the impact of these structural motifs on their closing to bifunctional POSS cages are also investigated. Our studies highlight that analyzing experimental results alone, without incorporating complex theoretical investigations, does not offer a clear understanding of the reactions involving multiple simultaneously reacting substrates, which may also undergo further transformations, potentially complicating the conventional pathways of classic corner-opening/capping reactions.

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.

How the Metal Ion Affects the 1H NMR Chemical Shift Values of Schiff Base Metal Complexes: Rationalization by DFT Calculations

The chemical shift (CS) values obtained by 1H NMR spectroscopy for the hydrogen atoms of a tetradentate N2O2-substituted Salphen ligand (H2L1) are differently shifted in its complexes of nickel(II), palladium(II), platinum(II), and zinc(II), all bearing the same charge on the metal ions. To rationalize the observed trends, DFT calculations have been performed in the implicit d6-DMSO solvent in terms of the electronic effects induced by the metal ion and of the nature and strength of the metal-N and metal-O bonds. Overall, the results obtained point out that, in the complexes involving group 10 elements, the CS values show the greater shift when considering the two hydrogen atoms at a shorter distance from the coordinated metal center and follow the decreasing metal charge in the order Ni > Pd > Pt. This trend suggests a more covalent character of the ligand-metal bonds with the increase of the metal atomic number. Furthermore, a slightly poorer agreement between experimental and calculated data is observed in the presence of the nickel(II) ion. Such discrepancy is explained by the formation of stacked oligomers, aimed at minimizing the repulsive interactions with the polar DMSO solvent.

A Holistic View of the Interactions between Electron-Deficient Systems: Clustering of Beryllium and Magnesium Hydrides and Halides

In the search for common bonding patterns in pure and mixed clusters of beryllium and magnesium derivatives, the most stable dimers and trimers involving BeX2 and MgX2 (X = H, F, Cl) have been studied in the gas phase using B3LYP and M06-2X DFT methods and the G4 ab initio composite procedure. To obtain some insight into their structure, stability, and bonding characteristics, we have used two different energy decomposition formalisms, namely MBIE and LMO-EDA, in parallel with the analysis of the electron density with the help of QTAIM, ELF, NCIPLOT, and AdNDP approaches. Some interesting differences are already observed in the dimers, where the stability sequence observed for the hydrides differs entirely from that of the fluorides and chlorides. Trimers also show some peculiarities associated with the presence of compact trigonal cyclic structures that compete in stability with the more conventional hexagonal and linear forms. As observed for dimers, the stability of the trimers changes significantly from hydrides to fluorides or chlorides. Although some of these clusters were previously explored in the literature, the novelty of this work is to provide a holistic approach to the entire series of compounds by using chemical bonding tools, allowing us to understand the stability trends in detail and providing insights for a significant number of new, unexplored structures.

Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost

We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.

libMBD: A general-purpose package for scalable quantum many-body dispersion calculations

Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.

Assignment of the absolute configuration of molecules that are chiral by virtue of deuterium substitution using chiral tag molecular rotational resonance spectroscopy

Chiral tag molecular rotational resonance (MRR) spectroscopy is used to assign the absolute configuration of molecules that are chiral by virtue of deuterium substitution. Interest in the improved performance of deuterated active pharmaceutical ingredients has led to the development of precision deuteration reactions. These reactions often generate enantioisotopomer reaction products that pose challenges for chiral analysis. Chiral tag rotational spectroscopy uses noncovalent derivatization of the enantioisotopomer to create the diastereomers of the 1:1 molecular complexes of the analyte and a small, chiral molecule. Assignment of the absolute configuration requires high-confidence determinations of the structures of these weakly bound complexes. A general search method, CREST, is used to identify candidate geometries. Subsequent geometry optimization using dispersion corrected density functional theory gives equilibrium geometries with sufficient accuracy to identify the isomers of the chiral tag complexes produced in the pulsed jet expansion used to introduce the sample into the MRR spectrometer. Rotational constant scaling based on the fact that the diastereomers have the same equilibrium geometry gives accurate predictions allowing identification of the homochiral and heterochiral tag complexes and, therefore, assignment of absolute configuration. The method is successfully applied to three oxygenated substrates from enantioselective Cu-catalyzed alkene transfer hydrodeuteration reaction chemistry.

(3+2)-Cycloadditions of Levoglucosenone (LGO) with Fluorinated Nitrile Imines Derived from Trifluoroacetonitrile: An Experimental and Computational Study

The in situ-generated N-aryl nitrile imines derived from trifluoroacetonitrile smoothly undergo (3+2)-cycloadditions onto the enone fragment of the levoglucosenone molecule, yielding the corresponding, five-membered cycloadducts. In contrast to the 'classic' C(Ph),N(Ph) nitrile imine, reactions with fluorinated C(CF3),N(Ar) analogues lead to stable pyrazolines in a chemo- and stereoselective manner. Based on the result of X-ray single crystal diffraction analysis, their structures were established as exo-cycloadducts with the location of the N-Ar terminus of the 1,3-dipole at the α-position of the enone moiety. The DFT computation demonstrated that the observed reaction pathway results from the strong dominance of kinetic control over thermodynamic control.

Real-Time Extension of TAO-DFT

Thermally assisted occupation density functional theory (TAO-DFT) has been an efficient electronic structure method for studying the ground-state properties of large electronic systems with multi-reference character over the past few years. To explore the time-dependent (TD) properties of electronic systems (e.g., subject to an intense laser pulse), in this work, we propose a real-time (RT) extension of TAO-DFT, denoted as RT-TAO-DFT. Moreover, we employ RT-TAO-DFT to study the high-order harmonic generation (HHG) spectra and related TD properties of molecular hydrogen H2 at the equilibrium and stretched geometries, aligned along the polarization of an intense linearly polarized laser pulse. The TD properties obtained with RT-TAO-DFT are compared with those obtained with the widely used time-dependent Kohn-Sham (TDKS) method. In addition, issues related to the possible spin-symmetry breaking effects in the TD properties are discussed.

AmberTools

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Smooth particle mesh Ewald-integrated stochastic Lanczos many-body dispersion algorithm

We derive and implement an alternative formulation of the Stochastic Lanczos algorithm to be employed in connection with the Many-Body Dispersion model (MBD). Indeed, this formulation, which is only possible due to the Stochastic Lanczos' reliance on matrix-vector products, introduces generalized dipoles and fields. These key quantities allow for a state-of-the-art treatment of periodic boundary conditions via the O(Nlog(N)) Smooth Particle Mesh Ewald (SPME) approach which uses efficient fast Fourier transforms. This SPME-Lanczos algorithm drastically outperforms the standard replica method which is affected by a slow and conditionally convergence rate that limits an efficient and reliable inclusion of long-range periodic boundary conditions interactions in many-body dispersion modelling. The proposed algorithm inherits the embarrassingly parallelism of the original Stochastic Lanczos scheme, thus opening up for a fully converged and efficient periodic boundary conditions treatment of MBD approaches.

Diastereoselective synthesis and structure-affinity relationships of σ1 receptor ligands with spirocyclic scaffold

Spirocyclic scaffolds play an increasing role in drug discovery as they define a rigid three-dimensional space to increase specific interactions with protein binding sites. Herein, a spirocyclic center was introduced into the lead compound 1 to rigidify its flexible benzylaminoethyl side chain. The key step of the synthesis was the reaction of different α,β-unsaturated amides 6 and 13-16 with methyl acrylate in the presence of TBDMSOTf. DFT calculations explain the mechanism of this transformation as concerted Diels-Alder reaction (functionals B3LYP and TPSS) or double (aza)-Michael addition (functionals PBE and wB97X-D). After separation of the diastereomeric spirocyclic products 8 and 17-20, LiAlH4 reduction provided the spirocyclic hydroxymethyl piperidines 21a,b-25a,b showing low nanomolar σ1 affinity (Ki < 100 nM). trans-Configured ligands (a-series) showed higher or equal σ1 affinity and higher selectivity over σ2 receptors and GluN2B-NMDA receptors than their cis-configured analogs (b-series). The additional hydroxymethyl moiety brings the log D7.4 value in a promising range. The high σ1 affinity (Ki = 3.6 nM) and the low lipophilicity result in the highest lipophilic ligand efficiency for the dispiro compound 23a (LLE = 6.0). The spirocyclic compounds reported herein and in particular the dispiro compound 23a demonstrate that ligands containing a large number of sp3 C-atoms possess favorable pharmacological (σ1 receptor affinity, receptor selectivity) and physicochemical properties (log D7.4 value) resulting in promising LLE.

Electron density-based protocol to recover the interacting quantum atoms components of intermolecular binding energy

A new approach for obtaining interacting quantum atoms-defined components of binding energy of intermolecular interactions, which bypasses the use of standard six-dimensional integrals and two-particle reduced density matrix (2-RDM) reconstruction, is proposed. To examine this approach, three datasets calculated within the density functional theory framework using the def2-TZVP basis have been explored. The first two, containing 53 weakly bound bimolecular associates and 13 molecular clusters taken from the crystal, were used in protocol refinement, and the third one containing other 20 bimolecular and three cluster systems served as a validation reference. In addition, to verify the performance of the proposed approach on an exact 2-RDM, calculations within the coupled cluster formalism were performed for part of the first set systems using the cc-pVTZ basis set. The process of optimization of the proposed parametric model is considered, and the role of various energy contributions in the formation of non-covalent interactions is discussed with regard to the obtained trends.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields cosκos-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the N4 formal scaling, cosκos-SPL2 is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs. The accuracy of cosκos-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5.

Mechanism and Selectivity of Copper-Catalyzed Bromination of Distal C(sp3)-H Bonds

Unactivated C(sp3)-H bonds are the most challenging substrate class for transition metal-catalyzed C-H halogenation. Recently, the Yu group [Liu, T.; Myers, M. C.; Yu, J. Q. Angew. Chem., Int. Ed.2017, 56 (1), 306-309] has demonstrated that a CuII/phenanthroline catalyst and BrN3, generated in situ from NBS and TMSN3 precursors, can achieve selective C-H bromination distal to a directing group. The current understanding of the mechanism of this reaction has left numerous questions unanswered. Here, we investigated the mechanism of Cu-catalyzed C(sp3)-H bromination with distal site selectivity using density functional theory calculations. We found that this reaction starts with the Br-atom transfer from BrN3 to the Cu center that occurs via a small energy barrier at the singlet-triplet state seam of crossing. In the course of this reaction, the presence of the N-H bond in the substrate is critically important and acts as a directing group for enhancing the stability of the catalyst-substrate interaction and for the recruitment of the substrate to the catalyst. The required C-centered radical substrate formation occurs via direct C-H dehydrogenation by the Cu-coordinated N3 radical, rather than via the previously proposed N-H bond dehydrogenation and then the 1,5-H transfer from the γ-(C-H) bond to the N-radical center pathway. The C-H bond activation by the azide radical is a regioselectivity-controlling step. The following bromination of the C-centered radical by the Cu-coordinated bromine completes the product formation. This reaction step is the rate-limiting step, occurs at the singlet-to-triplet state seam of the crossing point, and is exergonic.

Is there a Covalent Au(I)-Au(I) Bond in the trans-(AuX)2 (X = F, Cl, Br, I) Structure? A Post-Hartree-Fock and Density Functional Theory Study

We present an exhaustive exploration of the driving forces dominating the interaction between gold atoms in the trans-(AuX)2, where X is a halogen ligand. This work provides insights into the nature of the gold-gold contact in the trans-(AuX)2. The geometries and energies were calculated at the MP2, CCSD(T), and DFT-D3(BJ) (B3LYP, PBE, and TPSS) levels of theory. The results show a short Au-Au distance, typical of a covalent bond, but with a weak interaction energy associated with noncovalent interactions. It is established that the physical contributions from polarization and the electronic correlation forces are the most relevant at the post-Hartree-Fock level of theory. Also, the electrostatic term is attractive but with low contribution. Finally, the Wiberg indices and NBO analysis exposed a small covalent character between the gold atoms, revealing that this contribution is insufficient to explain the stability of the dimers. It is concluded that a sum of contributions makes it possible to establish an attraction between the gold atoms in the dimers studied beyond a classical aurophilic interaction.

Salt Confined in MIL-101(Cr)-Tailoring the Composite Sorbents for Efficient Atmospheric Water Harvesting

The adsorption method for atmospheric water harvesting (AWH) is considered as a promising heat-driven technology for potable water supply in arid regions. This research is focused on novel composite sorbents based on hygroscopic salts loaded in the pores of MIL-101(Cr) developed for AWH. The composites based on LiCl, LiBr, CaCl2 , and Ca(NO3 )2 were synthesized and comprehensively studied by SEM, XRD, N2 adsorption, and thermogravimetric methods. We evidence that the CaCl2 /MIL-101(Cr) composite demonstrates a high net water uptake of 0.52-0.59 g_(H2 O)/g_(composite) per cycle under conditions of Saudi Arabia and the Sahara desert as the reference regions with extra-dry climate, which exceeds the appropriate values for other adsorbents. It is shown that water adsorption on the composite cannot be presented as a combination of the adsorption on the components, thus indicating a synergistic effect. A detailed characterization of water coordination, mobility, and hydrogen bonding within the confined CaCl2 hydrates and salt solution using solid-state 2 H NMR spectroscopy has been performed. It is established that pore confinement promotes a prolonged transition to a dynamically melted state of the hydrated salt and a notable decrease of the melting temperature, which facilitates the molecular transport of water and causes the alteration of sorption properties of CaCl2 inside MIL-101 pores. Finally, the performance of AWH employing CaCl2 /MIL-101(Cr) was evaluated in terms of the fractions of water extracted and collected, and the specific energy consumption, demonstrating its high potential for AWH.

Sepiolite as a novel polysulfide trapper for energy applications: an electrochemical, X-ray spectroscopic and DFT study

Capacity retention is a critical property to enhance in electrochemical storage systems applied to renewable energy. In lithium-sulfur (Li-S) batteries, the capacity fade resulting from the shuttle effect of polysulfides is a major obstacle to their practical application. Sepiolite, an eco-friendly earth-abundant clay with suitable surface chemistry for anchoring and retaining various molecules and structures, was studied as a cathode additive to mitigate the shuttle effect using experimental and theoretical approaches. Electrochemical measurements, spectroscopy, and ab initio calculations were performed to describe the mechanism and interfaces involved in polysulfide retention using 2 wt% of sepiolite as an additive in Li-S batteries. The results showed that the addition of sepiolite significantly improved the capacity retention during battery cycling. Spectroscopic analysis revealed that the effective sepiolite-polysulfide interface was governed by oxidized sulfur species. Additionally, ab initio studies showed a highly exothermic adsorption both inside and outside the sepiolite pore. This study demonstrates the potential use of eco-friendly, low-cost, non-toxic, natural, and abundant materials as additives to increase capacity retention.

In Situ Cocrystallization via Spray Drying with Polymer as a Strategy to Prevent Cocrystal Dissociation

The aim of the present study was to investigate how different polymers affect the dissociation of cocrystals prepared by co-spray-drying active pharmaceutical ingredient (API), coformer, and polymer. Diclofenac acid-l-proline cocrystal (DPCC) was selected in this study as a model cocrystal due to its previously reported poor physical stability in a high-humidity environment. Polymers investigated include polyvinylpyrrolidone (PVP), poly(1-vinylpyrrolidone-co-vinyl acetate) (PVPVA), hydroxypropyl methyl cellulose, hydroxypropylmethylcellulose acetate succinate, ethyl cellulose, and Eudragit L-100. Terahertz Raman spectroscopy (THz Raman) and powder X-ray diffraction (PXRD) were used to monitor the cocrystal dissociation rate in a high-humidity environment. A Raman probe was used in situ to monitor the extent of the dissociation of DPCC and DPCC in crystalline solid dispersions (CSDs) with polymer when exposed to pH 6.8 phosphate buffer and water. The solubility of DPCC and solid dispersions of DPCC in pH 6.8 phosphate buffer and water was also measured. The dissociation of DPCC was water-mediated, and more than 60% of DPCC dissociated in 18 h at 40 °C and 95% RH. Interestingly, the physical stability of the cocrystal was effectively improved by producing CSDs with polymers. The inclusion of just 1 wt % polymer in a CSD with DPCC protected the cocrystal from dissociation over 18 h under the same conditions. Furthermore, the CSD with PVPVA was still partially stable, and the CSD with PVP was stable (undissociated) after 7 days. The superior stability of DPCC in CSDs with PVP and PVPVA was also demonstrated when systems were exposed to water or pH 6.8 phosphate buffer and resulted in higher dynamic solubility of the CSDs compared to DPCC alone. The improvement in physical stability of the cocrystal in CSDs was thought to be due to an efficient mixing between polymer and cocrystal at the molecular level provided by spray drying and in situ gelling of polymer. It is hypothesized that polymer chains could undergo gelling in situ and form a physical barrier, preventing cocrystal interaction with water, which contributes to slowing down the water-mediated dissociation.

Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations

Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.

Analyzing Fluoride Binding by Group 15 Lewis Acids: Pnictogen Bonding in the Pentavalent State

We report the results of a computational investigation into fluoride binding by a series of pentavalent pnictogen Lewis acids: pnictogen pentahalides (PnX5), tetraphenyl pnictogeniums (PnPh4+), and triphenyl pnictogen tetrachlorocatecholates (PnPh3Cat). Activation strain and energy decomposition analyses of the Lewis adducts not only clearly delineate the electrostatic and orbital contributions to these acid-base interactions but also highlight the importance of Pauli repulsion and molecular flexibility in determining relative Lewis acidity among the pnictogens.

Stability of Single Gold Atoms on Defective and Doped Diamond Surfaces

Polycrystalline boron-doped diamond (BDD) is widely used as a working electrode material in electrochemistry, and its properties, such as its stability, make it an appealing support material for nanostructures in electrocatalytic applications. Recent experiments have shown that electrodeposition can lead to the creation of stable small nanoclusters and even single gold adatoms on the BDD surfaces. We investigate the adsorption energy and kinetic stability of single gold atoms adsorbed onto an atomistic model of BDD surfaces by using density functional theory. The surface model is constructed using hybrid quantum mechanics/molecular mechanics embedding techniques and is based on an oxygen-terminated diamond (110) surface. We use the hybrid quantum mechanics/molecular mechanics method to assess the ability of different density functional approximations to predict the adsorption structure, energy, and barrier for diffusion on pristine and defective surfaces. We find that surface defects (vacancies and surface dopants) strongly anchor adatoms on vacancy sites. We further investigated the thermal stability of gold adatoms, which reveals high barriers associated with lateral diffusion away from the vacancy site. The result provides an explanation for the high stability of experimentally imaged single gold adatoms on BDD and a starting point to investigate the early stages of nucleation during metal surface deposition.

Evolutionary Algorithm-Based Crystal Structure Prediction of CuxZnyOz Ternary Oxides

Binary zinc(II) oxide (ZnO) and copper(II) oxide (CuO) are used in a number of applications, including optoelectronic and semiconductor applications. However, no crystal structures have been reported for ternary Cu-Zn-O oxides. In that context, we investigated the structural characteristics and thermodynamics of CuxZnyOz ternary oxides to map their experimental feasibility. We combined evolutionary crystal structure prediction and quantum chemical methods to investigate potential CuxZnyOz ternary oxides. The USPEX algorithm and density functional theory were used to screen over 4000 crystal structures with different stoichiometries. When comparing compositions with non-magnetic CuI ions, magnetic CuII ions, and mixed CuI-CuII compositions, the magnetic Cu2Zn2O4 system is thermodynamically the most favorable. At ambient pressures, the thermodynamically most favorable ternary crystal structure is still 2.8 kJ/mol per atom higher in Gibbs free energy compared to experimentally known binary phases. The results suggest that thermodynamics of the hypothetical CuxZnyOz ternary oxides should also be evaluated at high pressures. The predicted ternary materials are indirect band gap semiconductors.

Structures and Spectroscopic Properties of Polysulfide Radical Anions: A Theoretical Perspective

The potential involvement of polysulfide radical anions Sn•- is a recurring theme in discussions of the basic and applied chemistry of elemental sulfur. However, while the spectroscopic features for n = 2 and 3 are well-established, information on the structures and optical characteristics of the larger congeners (n = 4-8) is sparse. To aid identification of these ephemeral species we have performed PCM-corrected DFT calculations to establish the preferred geometries for Sn•- (n = 4-8) in the polar media in which they are typically generated. TD-DFT calculations were then used to determine the number, nature and energies of the electronic excitations possible for these species. Numerical reliability of the approach was tested by comparison of the predicted and experimental excitation energies found for S2•- and S3•-. The low-energy (near-IR) transitions found for the two acyclic isomers of S4•- (C2h and C2v symmetry) and for S5•- (Cs symmetry) can be understood by extension of the simple HMO π-only chain model that serves for S2•- and S3•-. By contrast, the excitations predicted for the quasi-cyclic structures Sn•- (n = 6-8) are better described in terms of σ → σ* processes within a localized 2c-3e manifold.

Low-Frequency Raman Spectroscopy of Pure and Cocrystallized Mycophenolic Acid

The aqueous solubility of solid-state pharmaceuticals can often be enhanced by cocrystallization with a coformer to create a binary cocrystal with preferred physical properties. Greater understanding of the internal and external forces that dictate molecular structure and intermolecular packing arrangements enables more efficient design of new cocrystals. Low-frequency (sub-200 cm^-1) Raman spectroscopy experiments and solid-state density functional theory simulations have been utilized together to investigate the crystal lattice vibrations of mycophenolic acid, an immunosuppressive drug, in its pure form and as a cocrystal with 2,2'-dipyridylamine. The lattice vibrations primarily consist of large-amplitude translations and rotations of the crystal components, thereby providing insights into the critical intermolecular forces governing cohesion of the molecular solids. The simulations reveal that despite mycophenolic acid having a significantly unfavorable conformation in the cocrystal as compared to the pure solid, the cocrystal exhibits greater thermodynamic stability over a wide temperature range. The energetic penalty due to the conformational strain is more than compensated for by the strong intermolecular forces between the drug and 2,2'-dipyridylamine. Quantifying the balance of internal and external energy factors in cocrystal formation indicates a path forward in the development of future mycophenolic acid cocrystals.