Dispersion-Corrected Mean-Field Electronic Structure Methods

The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations

The elements of the p-block of the periodic table are of high interest in various chemical and technical applications like frustrated Lewis-pairs (FLP) or opto-electronics. However, high-quality benchmark data to assess approximate density functional theory (DFT) for their theoretical description are sparse. In this work, we present a benchmark set of 604 dimerization energies of 302 "inorganic benzenes" composed of all non-carbon p-block elements of main groups III to VI up to polonium. This so-called IHD302 test set comprises two classes of structures formed by covalent bonding and by weaker donor-acceptor (WDA) interactions, respectively. Generating reliable reference data with ab initio methods is challenging due to large electron correlation contributions, core-valence correlation effects, and especially the slow basis set convergence. To compute reference values for these dimerization reactions, after thorough testing, we applied a computational protocol using state-of-the-art explicitly correlated local coupled cluster theory termed PNO-LCCSD(T)-F12/cc-VTZ-PP-F12(corr.). It includes a basis set correction at the PNO-LMP2-F12/aug-cc-pwCVTZ level. Based on these reference data, we assess 26 DFT methods in combination with three different dispersion corrections and the def2-QZVPP basis set, five composite DFT approaches, and five semi-empirical quantum mechanical methods. For the covalent dimerizations, the r2SCAN-D4 meta-GGA, the r2SCAN0-D4 and ωB97M-V hybrids, and the revDSD-PBEP86-D4 double-hybrid functional are found to be the best-performing methods among the evaluated functionals of the respective class. However, since def2 basis sets for the 4th period are not associated to relativistic pseudo-potentials, we obtained significant errors in the covalent dimerization energies (up to 6 kcal mol-1) for molecules containing p-block elements of the 4th period. Significant improvements were achieved for systems containing 4th row elements by using ECP10MDF pseudopotentials along with re-contracted aug-cc-pVQZ-PP-KS basis sets introduced in this work with the contraction coefficients taken from atomic DFT (PBE0) calculations. Overall, the IHD302 set represents a challenge to contemporary quantum chemical methods. This is due to a large number of spatially close p-element bonds which are underrepresented in other benchmark sets, and the partial covalent bonding character for the WDA interactions. The IHD302 set may be helpful to develop more robust and transferable approximate quantum chemical methods in the future.

Local Energy Decomposition Analysis of London Dispersion Effects: From Simple Model Dimers to Complex Biomolecular Assemblies

ConspectusLondon dispersion (LD) forces are ubiquitous in chemistry, playing a pivotal role in a wide range of chemical processes. For example, they influence the structure of molecular crystals, the selectivity of organocatalytic transformations, and the formation of biomolecular assemblies. Harnessing these forces for chemical applications requires consistent quantification of the LD energy across a broad and diverse spectrum of chemical scenarios. Despite the great progress made in recent years in the development of experimental strategies for LD quantification, quantum chemical methods remain one of the most useful tools in the hand of chemists for the study of these weak interactions. Unfortunately, the accurate quantification of LD effects in complex systems poses many challenges for electronic structure theories. One of the problems stems from the fact that LD forces originate from long-range electronic dynamic correlation, and hence, their rigorous description requires the use of complex, highly correlated wave function-based methods. These methods typically feature a steep scaling with the system size, limiting their applicability to small model systems. Another core challenge lies in disentangling short-range from long-range dynamic correlation, which from a rigorous quantum mechanical perspective is not possible.In this Account, we describe our research endeavors in the development of broadly applicable computational methods for LD quantification in molecular chemistry as well as challenging applications of these schemes in various domains of chemical research. Our strategy lies in the use of local correlation theories to reduce the computational cost associated with complex electronic structure methods while providing at the same time a simple means of decomposition of dynamic correlation into its long-range and short-range components. In particular, the local energy decomposition (LED) scheme at the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) level has emerged as a powerful tool in our research, offering a clear-cut quantitative definition of the LD energy that remains valid across a plethora of different chemical scenarios. Typical applications of this scheme are examined, encompassing protein-ligand interactions and reactivity studies involving many fragments and complex electronic structures. In addition, our research also involves the development of novel cost-effective methodologies, which exploit the LED definition of the LD energy, for LD energy quantification that are, in principle, applicable to systems with thousands of atoms. The Hartree-Fock plus London Dispersion (HFLD) scheme, correcting the HF interaction energy using an approximate CCSD(T)-based LD energy, is a useful, parameter-free electronic structure method for the study of LD effects in systems with hundreds of molecular fragments. However, the usefulness of the LED scheme reaches beyond providing an interpretation of the calculated DLPNO-CCSD(T) or DLPNO-MP2 interaction energies. For example, the dispersion energies obtained from the LED can be fruitfully used in order to parametrize semiempirical dispersion models. We will demonstrate this in the context of an emerging semiempirical method, namely, the Natural Orbital Tied Constructed Hamiltonian (NOTCH) method. NOTCH incorporates LED-derived LD energies and shows promising accuracy at a minimum amount of empiricism. Thus, it holds substantial promise for large and complex systems.

Hydroxy-Functionalized Ionic Liquids under Pressure: The Influence on Hydrogen Bonding between Ions of Opposite and Like Charges

Hydroxy functionalization of cations in ionic liquids (ILs) can lead to formation of hydrogen bonds between their OH groups, resulting in so-called (c-c) H-bonds. Thereby, the (c-c) H-bonds compete with regular H-bonds (c-a) between the OH groups and the anions. Polarizable cations, weakly interacting anions, and long alkyl chains at the cation support the propensity for the formation of (c-c) H-bonds. At low temperatures, the equilibrium between (c-c) and (c-a) H-bonds is strongly shifted in favor of the cation-cation interaction. Herein, we clarify the pressure dependence on (c-c) and (c-a) H-bond distributions in the IL 1-(2-hydroxyethyl)-3-methylimidazolium hexafluorophosphate [HOC2C1Im][PF6], in mixtures of [HOC2C1Im][PF6] with the nonhydroxy-functionalized IL 1-propyl-3-methylimidazolium hexafluorophosphate [C3C1Im][PF6] and in [HOC2C1Im][PF6] including trace amounts of water. The infrared (IR) spectra provide clear evidence that the (c-c) H-bonds diminish with increasing pressure in favor of the (c-a) H-bonds. Adding trace amounts of water results in enhanced (c-c) clustering due to cooperative effects. At ambient pressure, the water molecules are involved in the (c-c) H-bond motifs. Increasing pressure leads to squeezing them out of H-bond clusters, finally resulting in demixing of water and the IL at the microscopic level.

A Strain-Promoted Divergent Chemical Steroidation Unveils Potent Anti-Inflammatory Pseudo-Steroidal Glycosides

The development of novel agents with immunoregulatory effects is a keen way to combat the growing threat of inflammatory storms to global health. To synthesize pseudo-steroidal glycosides tethered by ether bonds with promising immunomodulatory potential, we develop herein a highly effective deoxygenative functionalization of a novel steroidal donor (steroidation) facilitated by strain-release, leveraging cost-effective and readily available Sc(OTf)3 catalysis. This transformation produces a transient steroid-3-yl carbocation which readily reacts with O-, C-, N-, S-, and P-nucleophiles to generate structurally diverse steroid derivatives. DFT calculations were performed to shed light on the mechanistic details of the regioselectivity, underlying an acceptor-dependent steroidation mode. This approach can be readily extended to the etherification of sugar alcohols to enable the achievement of a diversity-oriented, pipeline-like synthesis of pseudo-steroidal glycosides in good to excellent yields with complete stereo- and regiospecific control for anti-inflammatory agent discovery. Immunological studies have demonstrated that a meticulously designed cholesteryl disaccharide can significantly suppress interleukin-6 secretion in macrophages, exhibiting up to 99% inhibition rates compared to the negative control. These findings affirm the potential of pseudo-steroidal glycosides as a prospective category of lead agents for the development of novel anti-inflammatory drugs.

Understanding the effect of density functional choice and van der Waals treatment on predicting the binding configuration, loading, and stability of amine-grafted metal organic frameworks

Metal organic frameworks (MOFs) are crystalline, three-dimensional structures with high surface areas and tunable porosities. Made from metal nodes connected by organic linkers, the exact properties of a given MOF are determined by node and linker choice. MOFs hold promise for numerous applications, including gas capture and storage. M2(4,4'-dioxidobiphenyl-3,3'-dicarboxylate)-henceforth simply M2(dobpdc), with M = Mg, Mn, Fe, Co, Ni, Cu, or Zn-is regarded as one of the most promising structures for CO2 capture applications. Further modification of the MOF with diamines or tetramines can significantly boost gas species selectivity, a necessity for the ultra-dilute CO2 concentrations in the direct-air capture of CO2. There are countless potential diamines and tetramines, paving the way for a vast number of potential sorbents to be probed for CO2 adsorption properties. The number of amines and their configuration in the MOF pore are key drivers of CO2 adsorption capacity and kinetics, and so a validation of computational prediction of these quantities is required to suitably use computational methods in the discovery and screening of amine-functionalized sorbents. In this work, we study the predictive accuracy of density functional theory and related calculations on amine loading and configuration for one diamine and two tetramines. In particular, we explore the Perdew-Burke-Ernzerhof (PBE) functional and its formulation for solids (PBEsol) with and without the Grimme-D2 and Grimme-D3 pairwise corrections (PBE+D2/3 and PBEsol+D2/3), two revised PBE functionals with the Grimme-D2 and Grimme-D3 pairwise corrections (RPBE+D2/3 and revPBE+D2/3), and the nonlocal van der Waals correlation (vdW-DF2) functional. We also investigate a universal graph deep learning interatomic potential's (M3GNet) predictive accuracy for loading and configuration. These results allow us to identify a useful screening procedure for configuration prediction that has a coarse component for quick evaluation and a higher accuracy component for detailed analysis. Our general observation is that the neural network-based potential can be used as a high-level and rapid screening tool, whereas PBEsol+D3 gives a completely qualitatively predictive picture across all systems studied, and can thus be used for high accuracy motif predictions. We close by briefly exploring the predictions of relative thermal stability for the different functionals and dispersion corrections.

pyCOFBuilder: A Python Package for Automated Creation of Covalent Organic Framework Models Based on the Reticular Approach

Covalent organic frameworks (COFs) have gained significant popularity in recent years due to their unique ability to provide a high surface area and customizable pore geometry and chemistry, making them an ideal choice for a wide range of applications. However, exploring COFs experimentally can be arduous and time-consuming due to their immense number of potential structures. As a result, computational high-throughput studies have become an attractive option. Nevertheless, generating COF structures can also be a challenging and time-consuming task. To address this challenge, here, we introduce the pyCOFBuilder, an open-source Python package designed to facilitate the generation of COF structures for computational studies. The pyCOFBuilder software provides an easy-to-use set of functionalities to generate COF structures following the reticular approach. In this paper, we describe the implementation, main features, and capabilities of the pyCOFBuilder, demonstrating its utility for generating COF structures with varying topologies and chemical properties. pyCOFBuilder is freely available on GitHub at https://github.com/lipelopesoliveira/pyCOFBuilder.

Collaboration on Machine-Learned Potentials with IPSuite: A Modular Framework for Learning-on-the-Fly

The field of machine learning potentials has experienced a rapid surge in progress, thanks to advances in machine learning theory, algorithms, and hardware capabilities. While the underlying methods are continuously evolving, the infrastructure for their deployment has lagged. The community, due to these rapid developments, frequently finds itself split into groups built around different implementations of machine-learned potentials. In this work, we introduce IPSuite, a Python-driven software package designed to connect different methods and algorithms from the comprehensive field of machine-learned potentials into a single platform while also providing a collaborative infrastructure, helping ensure reproducibility. Furthermore, the data management infrastructure of the IPSuite code enables simple model sharing and deployment in simulations. Currently, IPSuite supports six state-of-the-art machine learning approaches for the fitting of interatomic potentials as well as a variety of methods for the selection of training data, running of ab initio calculations, learning-on-the-fly strategies, model evaluation, and simulation deployment.

Predicting valence tautomerism in diverse cobalt-dioxolene complexes: elucidation of the role of ligands and solvent

The ability of molecular switches to reversibly interconvert between different forms promises potential applications at the scale of single molecules up to bulk materials. One type of molecular switch comprises cobalt-dioxolene compounds that exhibit thermally-induced valence tautomerism (VT) interconversions between low spin Co(iii)-catecholate (LS-CoIII-cat) and high spin Co(ii)-semiquinonate (HS-CoII-sq) forms. Two families of these compounds have been investigated for decades but have generally been considered separately: neutral [Co(diox)(sq)(N2L)] and cationic [Co(diox)(N4L)]+ complexes (diox = generic dioxolene, N2L/N4L = bidentate/tetradentate N-donor ancillary ligand). Computational identification of promising new candidate compounds prior to experimental exploration is beneficial for environmental and cost considerations but requires a thorough understanding of the underlying thermochemical parameters that influence the switching. Herein, we report a robust approach for the analysis of both cobalt-dioxolene families, which involved a quantitative density functional theory-based study benchmarked with reliable quasi-experimental references. The best-performing M06L-D4/def2-TZVPP level of theory has subsequently been verified by the synthesis and experimental investigation of three new complexes, two of which exhibit thermally-induced VT, while the third remains in the LS-CoIII-cat form across all temperatures, in agreement with prediction. Valence tautomerism in solution is markedly solvent-dependent, but the origin of this has not been definitively established. We have extended our computational approach to elucidate the correlation of VT transition temperature with solvent stabilisation energy and change in dipole moment. This new understanding may inform the development of VT compounds for applications in soft materials including films, gels, and polymers.

Initial decomposition pathways of 1,1-diamino-2,2-dinitroethylene (α-FOX-7) in the condensed phase

The initial decomposition pathways of α-FOX-7 in the condensed phase (crystal) were investigated via density functional theory. Calculations were carried out using three FOX-7 systems with increasing complexity from 1-layer (sheet) via 2-layer (surface) to 3-layer (bulk). The encapsulated environment of the central α-FOX-7 molecule, where decomposition takes place, is reconstructed by neighbouring molecules following a crystal structure. A minimal number of neighbouring molecules that have an impact on the energetics of decomposition are identified among all surrounding molecules. The results show that the presence of intermolecular hydrogen bonds due to the encapsulated environment in the condensed phase decreases the sensitivity of α-FOX-7, i.e. it increases the barrier of decomposition, but it does not alter the initial decomposition pathways of the reaction compared to the gas phase. Moreover, increasing the complexity of the system from a single gas phase molecule via sheet and surface to bulk increases the decomposition barriers. The calculations reveal a remarkable agreement with experimental data [A. M. Turner, Y. Luo, J. H. Marks, R. Sun, J. T. Lechner, T. M. Klapötke and R. I. Kaiser, Exploring the Photochemistry of Solid 1, 1-Diamino-2, 2-Dinitroethylene (FOX-7) Spanning Simple Bon Ruptures, Nitro-to-Nitrite Isomerization, and Nonadiabatic Dynamics, J. Phys. Chem. A, 2022, 126, 29, 4747-4761] and suggest that the initial decomposition of α-FOX-7 likely takes place at the surface of the crystal.

Review and Theoretical Analysis of Fluorinated Radicals in Direct CAr-H Functionalization of (Hetero)arenes

We highlight key contributions in the field of direct radical CAr- H (hetero)aromatic functionalization involving fluorinated radicals. A compilation of Functional Group Transfer Reagents and their diverse activation mechanisms leading to the release of radicals are discussed. The substrate scope for each radical is analyzed and classified into three categories according to the electronic properties of the substrates. Density functional theory computational analysis provides insights into the chemical reactivity of several fluorinated radicals through their electrophilicity and nucleophilicity parameters. Theoretical analysis of their reduction potentials also highlights the remarkable correlation between electrophilicity and oxidizing ability. It is also established that highly fluorinated radicals (e.g. ⋅OCF3) are capable of engaging in single-electron transfer (SET) processes rather than radical addition, which is in good agreement with experimental literature data. A reactivity scale, based on activation barrier of addition of these radicals to benzene is also elaborated using the high accuracy DLPNO-(U)CCSD(T) method.

Realizing the effect of s-block metals on a charge transfer crystal of indol-2-one for enhanced NLO responses with efficient energetic offsets

CONTEXT

Due to their unique photophysical properties, organic charge transfer crystals are becoming promising materials for next-generation optoelectronic devices. This research paper explores the impact of s-block metals on a charge transfer crystal of indol-2-one for enhanced nonlinear optical (NLO) responses with efficient energetic offsets. The study reveals that alkali metals can enhance NLO performance due to their free electrons.

METHOD

The Perdew-Burke-Ernzerhof functional of DFT with dispersion correction (D3) was used, and the λmax values ranged between 596 and 669 nm, with the highest value for dichloromethane (DCM). Leveraging the unique properties of metals allowed for the development of nonlinear optical materials with improved performance and versatility. Softness (σ) values provide insight into electron density changes, with higher values indicating a greater tendency for changes and lower values indicating the opposite. The NLO results for the chromophores MMI1-MMI6 show varying linear polarizability (< α0 >) along with their first (β0) and second (γ0) hyperpolarizabilities. Chromophore MMI4 stands out with the highest NLO performance, having two potassium (K) atoms. Its < α0 > , β0, and γ0 values of 4.19, 7.09, and 17.43 (× 10-24 e.s.u), respectively, indicate a significant enhancement in NLO response compared to the other chromophores. The transitions involving (O20)LP → (C3-N5)π* and (O19)LP → (N12-C13)π* exhibit the highest level of stabilization, followed by (O23)π → (C10-C11)π*, while (C6-N12)π → (C6-C7)π* shows the lowest level of stabilization for chromophore MMI4. The present research work is facile in its nature, and it can be helpful for synthetic scientist to design the new materials for uniting crystal properties with metal doping for efficient NLO devices.

Novel thiazolidin-4-one benzenesulfonamide hybrids as PPARγ agonists: Design, synthesis and in vivo anti-diabetic evaluation

In the current study, two series of novel thiazolidin-4-one benzenesulfonamide arylidene hybrids 9a-l and 10a-f were designed, synthesized and tested in vitro for their PPARɣ agonistic activity. The phenethyl thiazolidin-4-one sulphonamide 9l showed the highest PPARɣ activation % by 41.7%. Whereas, the 3-methoxy- and 4-methyl-4-benzyloxy thiazolidin-4-one sulphonamides 9i, and 9k revealed moderate PPARɣ activation % of 31.7, and 32.8%, respectively, in addition, the 3-methoxy-3-benzyloxy thiazolidin-4-one sulphonamide 10d showed PPARɣ activation % of 33.7% compared to pioglitazone. Compounds 9b, 9i, 9k, 9l, and 10d revealed higher selectivity to PPARɣ over the PPARδ, and PPARα isoforms. An immunohistochemical study was performed in HepG-2 cells to confirm the PPARɣ protein expression for the most active compounds. Compounds 9i, 9k, and 10d showed higher PPARɣ expression than that of pioglitazone. Pharmacological studies were also performed to determine the anti-diabetic activity in rats at a dose of 36 mg/kg, and it was revealed that compounds 9i and 10d improved insulin secretion as well as anti-diabetic effects. The 3-methoxy-4-benzyloxy thiazolidin-4-one sulphonamide 9i showed a better anti-diabetic activity than pioglitazone. Moreover, it showed a rise in blood insulin by 4-folds and C-peptide levels by 48.8%, as well as improved insulin sensitivity. Moreover, compound 9i improved diabetic complications as evidenced by decreasing liver serum enzymes, restoration of total protein and kidney functions. Besides, it combated oxidative stress status and exerted anti-hyperlipidemic effect. Compound 9i showed a superior activity by normalizing some parameters and amelioration of pancreatic, hepatic, and renal histopathological alterations caused by STZ-induction of diabetes. Molecular docking studies, molecular dynamic simulations, and protein ligand interaction analysis were also performed for the newly synthesized compounds to investigate their predicted binding pattern and energies in PPARɣ binding site.

Distinction and Quantification of Noncovalent Dispersive and Hydrophobic Effects

The possibilities of comparing computational results of noncovalent interactions with experimental data are discussed, first with respect to intramolecular interactions. For these a variety of experimental data such as heats of formation, crystal sublimation heats, comparison with energy minimized structures, and spectroscopic data are available, but until now largely have not found widespread application. Early force field and QM/MP2 calculations have already shown that the sublimation heats of hydrocarbons can be predicted with an accuracy of ±1%. Intermolecular interactions in solution or the gas phase are always accompanied by difficult to compute entropic contributions, like all associations between molecules. Experimentally observed T∆S values contribute 10% to 80% of the total ∆G, depending on interaction mechanisms within the complexes, such as, e.g., hydrogen bonding and ion pairing. Free energies ∆G derived from equilibrium measurements in solution allow us to define binding increments ∆∆G, which are additive and transferable to a variety of supramolecular complexes. Data from more than 90 equilibrium measurements of porphyrin receptors in water indicate that small alkanes do not bind to the hydrophobic flat surfaces within a measuring limit of ∆G = ±0.5 kJ/mol, and that 20 functions bearing heteroatoms show associations by dispersive interactions with up to ∆G = 8 kJ/mol, roughly as a function of their polarizability. Aromatic systems display size-dependent affinities ∆G as a linear function of the number of π-electrons.

Speciation of Technetium Dibutylphosphate in the Third Phase Formed in the TBP/HNO3 Solvent Extraction System

The speciation of technetium in the nitric acid/dibutylphosphoric acid (HDBP)-n-dodecane system was studied by extended X-ray absorption fine structure (EXAFS) spectroscopy and theoretical methods. Tetravalent technetium, produced by the hydrazine reduction of TcO4- in 3 M HNO3, was extracted by HDBP in n-dodecane (30% by volume). During extraction, the splitting of the organic phase into a heavy phase and a light phase was observed. EXAFS analysis is consistent with the presence of Tc(NO3)3(DBP)(HDBP)2 in the light phase and Tc(NO3)2(DBP)2(HDBP)2 in the heavy phase. Density functional theory calculations at the B3LYP/6-31G* level confirm the stability of the proposed species and indicate that stereoisomers -mer- and fac-Tc(NO3)3(DBP)(HDBP)2 for the light phase and cis- and trans-Tc(NO3)2(DBP)2(HDBP)2 for the heavy phase] could coexist in the system (in the n-dodecane solution). Mechanisms of formation of Tc(NO3)3(DBP)(HDBP)2 and Tc(NO3)2(DBP)2(HDBP)2 are proposed.

Synthesis and structural properties of para-diselenopyrazines

Recently, dithienopyrazines have emerged as promising building blocks in the field of materials science, showcasing their potential as hole-transport materials in organic electronic devices. Herein, we report the synthesis of its heavier analogues, the diselenopyrazines, along with an analysis of their optoelectronic and structural properties. In the acquired crystal structures, interesting molecular packing motifs suitable for potential device fabrication were observed. Additionally, short contact interactions are present in one of the investigated molecules.

Flexibility of Binding Site is Essential to the Ca2+ Selectivity in EF-Hand Calcium-Binding Proteins

High binding affinity and selectivity of metal ions are essential to the function of metalloproteins. Thus, understanding the factors that determine these binding characteristics is of major interest for both fundamental mechanistic investigations and guiding of the design of novel metalloproteins. In this work, we perform QM cluster model calculations and quantum mechanics/molecular mechanics (QM/MM) free energy simulations to understand the binding selectivity of Ca2+ and Mg2+ in the wild-type carp parvalbumin and its mutant. While a nonpolarizable MM model (CHARMM36) does not lead to the correct experimental trend, treatment of the metal binding site with the DFTB3 model in a QM/MM framework leads to relative binding free energies (ΔΔGbind) comparable with experimental data. For the wild-type (WT) protein, the calculated ΔΔGbind is ∼6.6 kcal/mol in comparison with the experimental value of 5.6 kcal/mol. The good agreement highlights the value of a QM description of the metal binding site and supports the role of electronic polarization and charge transfer to metal binding selectivity. For the D51A/E101D/F102W mutant, different binding site models lead to considerable variations in computed binding affinities. With a coordination number of seven for Ca2+, which is shown by QM/MM metadynamics simulations to be the dominant coordination number for the mutant, the calculated relative binding affinity is ∼4.8 kcal/mol, in fair agreement with the experimental value of 1.6 kcal/mol. The WT protein is observed to feature a flexible binding site that accommodates a range of coordination numbers for Ca2+, which is essential to the high binding selectivity for Ca2+ over Mg2+. In the mutant, the E101D mutation reduces the flexibility of the binding site and limits the dominant coordination number of Ca2+ to be seven, thereby leading to reduced binding selectivity against Mg2+. Our results highlight that the binding selectivity of metal ions depends on both the structural and dynamical properties of the protein binding site.

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.

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.