Density fitting of intramonomer correlation effects in symmetry-adapted perturbation theory

Stacking Interactions of Druglike Heterocycles with Nucleobases

Stacking interactions contribute significantly to the interaction of small molecules with RNA, and harnessing the power of these interactions will likely prove important in the development of RNA-targeting inhibitors. To this end, we present a comprehensive computational analysis of stacking interactions between a set of 54 druglike heterocycles and the natural nucleobases. We first show that heterocycle choice can tune the strength of stacking interactions with nucleobases over a large range and that heterocycles favor stacked geometries that cluster around a discrete set of stacking loci characteristic of each nucleobase. Symmetry-adapted perturbation theory results indicate that the strengths of these interactions are modulated primarily by electrostatic and dispersion effects. Based on this, we present a multivariate predictive model of the maximum strength of stacking interactions between a given heterocycle and nucleobase that depends on molecular descriptors derived from the electrostatic potential. These descriptors can be readily computed using density functional theory or predicted directly from atom connectivity (e.g., SMILES). This model is used to predict the maximum possible stacking interactions of a set of 1854 druglike heterocycles with the natural nucleobases. Finally, we show that trivial modifications of standard (fixed-charge) molecular mechanics force fields reduce errors in predicted stacking interaction energies from around 2 kcal/mol to below 1 kcal/mol, providing a pragmatic means of predicting more reliable stacking interaction energies using existing computational workflows. We also analyze the stacking interactions between ribocil and a bacterial riboswitch, showing that two of the three aromatic heterocyclic components engage in near-optimal stacking interactions with binding site nucleobases.

Competition between Hydrogen and Chalcogen Bonding in Homodimers of Chalcogen Hydrides (H2X)2, X = O, S, Se, Te

The structural and chemical bonding motifs manifested by the competition between hydrogen and chalcogen bonding in the homodimers of chalcogen hydrides (H2X)2, where X = O, S, Se, Te have been characterized using high-level electronic structure calculations and analysis of the electron density based on Quasi-atomic orbital (QUAO) and the Symmetry-adapted perturbation theory (SAPT) methods. The QUAO analysis clearly identifies a three-center interaction responsible for either hydrogen or chalcogen bonds: in the former, the σ-bond connecting the donor and hydrogen atom participating in the hydrogen bond interacts with the lone pair on the nucleophile acceptor via the hydrogen atom, while in the latter this same σ-bond interacts with the nucleophile lone pair via the donor chalcogen. The number of minimum energy structures increase dramatically from one for (H2O)2, three for (H2S)2, four for (H2Se)2, and finally six for (H2Te)2. The emergence of the chalcogen-bonded arrangements appears for (H2S)2 with their subsequent energetic stabilization over the hydrogen-bonded minima manifesting in (H2Se)2 and (H2Te)2. In particular, one of the (H2S)2 , two of the (H2Se)2, and three of the (H2Te)2 dimers are chalcogen bonded. Induction plays a small but important role in stabilizing hydrogen over chalcogen-bonded structures, while dispersion is more important for chalcogen bonds.

What is the Exchange Repulsion Energy? Insight by Partitioning into Physically Meaningful Contributions

It is shown that the exchange repulsion energy, Exr, can be rationalized by partitioning the respective energy expression for two systems with Hartree-Fock orbitals into physically meaningful contributions. A division of Exr into a positive kinetic and a negative potential part is possible, but these contributions correlate only poorly with the actual exchange repulsion energy. A more meaningful partitioning is derived, where all kinetic energy contributions are collected in a term that vanishes for exact Hartree-Fock orbitals due to their stationarity conditions. The remaining terms can be distinguished into an exchange integral contribution as well as contributions to the repulsion energy with two, three and four orbital indices. The forms, relationships and absolute sizes of these terms suggest an intuitive partitioning of the exchange repulsion energy into Molecular Orbital Pair Contributions to the Exchange repulsion energy (MOPCE). Insight into the analytic form and quantitative size of these contributions is provided by considering the3 Σ u + ( 1 σ g 1 σ u ) ${^3 \Sigma _u^+ (1\sigma _g 1\sigma _u )}$ state of the H2 molecule, the water dimer, as well as an argon atom interacting with Cl2 and N2.

Methyl formate synthesis via S N Acyl esterification on interstellar ice mantles

CONTEXT

Methyl formate (MF) has been detected in several interstellar environments, but whether or not the formation of this molecule takes place in the gas phase or on the ices of interstellar grains is still unclear. In this study, we explore the synthesis of methyl formate through the nucleophilic acyl substitution (SN Acyl) reaction between methanol (CH3 OH) and formic acid (HCOOH) on amorphous solid water, which is the main component of interstellar ice mantles.

METHODS

Using density functional theory (DFT), we model MF formation by sampling HCOOH in different catalytic sites on the water clusters with CH3 OH, and vice versa, for initial reactant configurations. We select the initial binding modes from the binding energy distributions of both reactant species. We assess the energy and synchronicity of the reaction by analyzing the reaction mechanisms through intrinsic reaction coordinate (IRC) energy, reaction force, and reaction electronic flux profiles. Using Wiberg bond order derivatives, we identify reaction events linked to hidden transition states that are encountered along the reaction coordinate.

In Silico Screening of CO2-Dipeptide Interactions for Bioinspired Carbon Capture

Carbon capture, sequestration and utilization offers a viable solution for reducing the total amount of atmospheric CO2 concentrations. On an industrial scale, amine-based solvents are extensively employed for CO2 capture through chemisorption. Nevertheless, this method is marked by the high cost associated with solvent regeneration, high vapor pressure, and the corrosive and toxic attributes of by-products, such as nitrosamines. An alternative approach is the biomimicry of sustainable materials that have strong affinity and selectivity for CO2. Bioinspired approaches, such as those based on naturally occurring amino acids, have been proposed for direct air capture methodologies. In this study, we present a database consisting of 960 dipeptide molecular structures, composed of the 20 naturally occurring amino acids. Those structures were analyzed with a novel computational workflow presented in this work that considers certain interaction sites that determine CO2 affinity. Density functional theory (DFT) and symmetry-adapted perturbation theory (SAPT) computations were performed for the calculation of CO2 interaction energies, which allowed to limit our search space to 400 unique dipeptide structures. Using this computational workflow, we provide statistical insights into dipeptides and their affinity for CO2 binding, as well as design principles that can further enhance CO2 capture through cooperative binding.

Carbonless DNA

Carbonless DNA was designed by replacing all carbon atoms in the standard DNA building blocks with boron and nitrogen, ensuring isoelectronicity. Electronic structure quantum chemistry methods (DFT(ωB97XD)/aug-cc-pVDZ) were employed to study both the individual building blocks and the larger carbon-free DNA fragments. The reliability of the results was validated by comparing selected structures and binding energies using more accurate methods such as MP2, CCSD, and SAPT2+3(CCD)δMP2. Carbonless analogs of DNA components, including cytosine, thymine, guanine, adenine, and deoxyribose, were investigated, showing strong resemblance to the carbon-based versions in terms of spatial structure, polarity, and molecular interaction capabilities. Complementary base pairs of the carbonless analogs exhibited a similar number and length of hydrogen bonds as those found in their carbon-containing counterparts, with binding energies for A-T and G-C analogs remaining comparable. Carbonless DNA fragments containing two and six base pairs were studied, revealing double-helix structures analogous to natural DNA. Structural parameters such as fragment size, hydrogen bond lengths, and rise per base pair were consistent with those observed in unmodified DNA. Docking simulations with a 12 base pair fragment and netropsin as a ligand indicated a slight shift in binding preference for the carbonless DNA through the minor groove, with an approximate 25% increase in binding affinity compared to natural DNA.

Evidence of strong O-H⋯C interactions involving apical pyramidane carbon atoms as hydrogen atom acceptors

Using high-level quantum chemical calculations, we predicted a strong O-H⋯C interaction between the apical carbon atoms of pyramidane and its derivatives and water molecules. Analysis of calculated electrostatic potential maps showed that there are areas of strong negative potential above apical carbon atoms in all studied structures. The results of quantum chemical calculations showed that the O-H⋯C interaction between the hydrogen atom of water and the apical carbon atom of pyramidane derivatives with four -CH3 substituents is unexpectedly strong, ΔECCSD(T)/CBS = -7.43 kcal mol-1. The strong hydrogen bonds were also predicted in the case of unsubstituted pyramidane (ΔECCSD(T)/CBS = -6.41 kcal mol-1) and pyramidane with four -OH substituents (ΔECCSD(T)/CBS = -5.87 kcal mol-1). Although there are not many crystal structures of pyramidane-like molecules, we extracted examples of pyramidal-shaped molecules with apical carbon atoms from the Cambridge Structural Database and analyzed their hydrogen-bonding patterns. Analysis of crystal structures confirmed the existence of short non-covalent contacts between apical carbon atoms and neighboring hydrogen atoms.

Simulations of γ-Valerolactone Solvents and Electrolytes for Lithium Batteries Using Polarizable Molecular Dynamics

In this paper, we present a molecular dynamics study of the structural and dynamical properties of γ-valerolactone (GVL) both as a standalone solvent and in electrolyte formulations for electrochemistry applications. This study involves developing a new parameterization of a polarizable forcefield and applying it to simulate pure GVL and selected salt solutions. The forcefield was validated with experimental bulk data and quantum mechanical calculations, with excellent agreement obtained in both cases. Specifically, two 1M electrolyte solutions of lithium bis(fluorosulfonyl)imide and lithium bis(oxalate)borate in GVL were simulated, focusing on their ionic transport and highlighting ion solvation structure. Ion pairing in the electrolytes was also investigated through enhanced sampling molecular dynamics, obtaining a detailed picture of the ion dynamics in the GVL solution.

How Flexible Is the Hydrogen Sulfide Molecule Structure? Influence of Hydrogen Sulfide Molecule Geometry on Its Hydrogen Bonds

The geometry of hydrogen sulfide was studied by calculating potential energy surface (PES) with over 1800 configurations. The calculations were performed at very accurate CCSD(T)/aug-cc-pvz5 level. The most stable geometry on the PES has bond angle (H-S-H) of 92.40° and bond length (S-H) of 1.338 Å. The PES shows that hydrogen sulfide is a quite flexible molecule. Namely, it can change the bonding angle (H-S-H) in the range of 15.6° (from 84.6° to 100.2°) and the bond lengths (S-H) in the range of 0.082 Å (from 1.299 Å to 1.381 Å) with an energy increase of only 1.0 kcal/mol. An influence of hydrogen sulfide geometry on its hydrogen bonds was studied on several hydrogen sulfide/hydrogen sulfide and water/hydrogen sulfide dimers. It showed that the change of hydrogen sulfide geometry does not influence the strength of hydrogen bond. Fully optimized geometries in gas and water solution phases revealed structural differences of both monomers and dimers in gas phase and water phase. SAPT analysis of the optimized dimer geometries showed that in all the dimers electrostatic is the most dominant contribution, while, in the dimers with hydrogen sulfide, the influence of dispersion contribution becomes quite pronounced.

Canonical coupled cluster binding benchmark for nanoscale noncovalent complexes at the hundred-atom scale

In this study, we introduce two datasets for nanoscale noncovalent binding, featuring complexes at the hundred-atom scale, benchmarked using coupled cluster with single, double, and perturbative triple [CCSD(T)] excitations extrapolated to the complete basis set (CBS) limit. The first dataset, L14, comprises 14 complexes with canonical CCSD(T)/CBS benchmarks, extending the applicability of CCSD(T)/CBS binding benchmarks to systems as large as 113 atoms. The second dataset, vL11, consists of 11 even larger complexes, evaluated using the local CCSD(T)/CBS method with stringent thresholds, covering systems up to 174 atoms. We compare binding energies obtained from local CCSD(T) and fixed-node diffusion Monte Carlo (FN-DMC), which have previously shown discrepancies exceeding the chemical accuracy threshold of 1 kcal/mol in large complexes, with the new canonical CCSD(T)/CBS results. While local CCSD(T)/CBS agrees with canonical CCSD(T)/CBS within binding uncertainties, FN-DMC consistently underestimates binding energies in π-π complexes by over 1 kcal/mol. Potential sources of error in canonical CCSD(T)/CBS are discussed, and we argue that the observed discrepancies are unlikely to originate from CCSD(T) itself. Instead, the fixed-node approximation in FN-DMC warrants further investigation to elucidate these binding discrepancies. Using these datasets as reference, we evaluate the performance of various electronic structure methods, semi-empirical approaches, and machine learning potentials for nanoscale complexes. Based on computational accuracy and stability across system sizes, we recommend MP2+aiD(CCD), PBE0+D4, and ωB97X-3c as reliable methods for investigating noncovalent interactions in nanoscale complexes, maintaining their promising performance observed in smaller systems.

Design and Evaluation of Pyridinyl Sulfonyl Piperazine LpxH Inhibitors with Potent Antibiotic Activity Against Enterobacterales

Enterobacterales, a large order of Gram-negative bacteria, including Escherichia coli and Klebsiella pneumoniae, are major causes of urinary tract and gastrointestinal infections, pneumonia, and other diseases in healthcare settings and communities. ESBL-producing Enterobacterales and carbapenem-resistant Enterobacterales can break down commonly used antibiotics, with some strains being resistant to all available antibiotics. This public health threat necessitates the development of novel antibiotics, ideally targeting new pathways in these bacteria. Gram-negative bacteria possess an outer membrane enriched with lipid A, a saccharolipid that serves as the membrane anchor of lipopolysaccharides and the active component of the bacterial endotoxin, causing septic shock. The biosynthesis of lipid A is crucial for the viability of Gram-negative bacteria, and as an essential enzyme in this process, LpxH has emerged as a promising target for developing novel antibiotics against multidrug-resistant Gram-negative pathogens. Here, we report the development of pyridinyl sulfonyl piperazine LpxH inhibitors. Among them, ortho-substituted pyridinyl compounds significantly boost LpxH inhibition and antibiotic activity over the original phenyl series. Structural and QM/MM analyses reveal that these improved activities are primarily due to the enhanced interaction between F141 of the LpxH insertion lid and the pyridinyl group. Incorporation of the N-methyl-N-phenyl-methanesulfonamide moiety into the pyridinyl sulfonyl piperazine backbone results in JH-LPH-106 and JH-LPH-107, both of which exhibit potent antibiotic activity against wild-type Enterobacterales such as K. pneumoniae and E. coli. JH-LPH-107 exhibits a low rate of spontaneous resistance and a high safety window in vitro, rendering it an excellent lead for further clinical development.

Beryllium as a Base: Complexes of Be(CO)3 with HX (X=F, Cl, Br, CN, NC, CCH, OH)

Beryllium chemistry is typically governed by its electron deficient character, but in some compounds it can act as a base. In order to understand better the unusual basicity of Be, we have systematically explored the complexes of one such compound, Be(CO)3, towards several hydrogen bond donors HX (X=F, Cl, Br, CN, NC, CCH, OH). For all complexes we find three different minima, two hydrogen bonded minima (to the Be or O atoms), and one weak beryllium bonded minimum. Further characterization of the interactions using a topological analysis of the electron density and Symmetry Adapted Perturbation Theory (SAPT) provide insight into the nature of these interactions. Overall these results highlight the capability of certain beryllium compounds to act as either a weak Lewis acid or, unconventionally, a Lewis base whose basicity towards hydrogen bonding is comparable to that of π systems.

Anisole-Water and Anisole-Ammonia Complexes in Ground and Excited (S1) States: A Multiconfigurational Symmetry-Adapted Perturbation Theory (SAPT) Study

Binary complexes of anisole have long been considered paradigm systems for studying microsolvation in both the ground and electronically excited states. We report a symmetry-adapted perturbation theory (SAPT) analysis of intermolecular interactions in anisole-water and anisole-ammonia complexes within the framework of the multireference SAPT(CAS) method. Upon the S1 ← S0 electronic transition, the hydrogen bond in the anisole-water dimer is weakened, which SAPT(CAS) shows to be determined by changes in the electrostatic energy. As a result, the water complex becomes less stable in the relaxed S1 state despite decreased Pauli repulsion. Stronger binding of the anisole-ammonia complex following electronic excitation is more nuanced and results from counteracting shifts in the repulsive (exchange) and attractive (electrostatic, induction and dispersion) forces. In particular, we show that the formation of additional binding N-H···π contacts in the relaxed S1 geometry is possible due to reduced Pauli repulsion in the excited state. The SAPT(CAS) interaction energies have been validated against the coupled cluster (CC) results and experimentally determined shifts of the S1 ← S0 anisole band. While for the hydrogen-bonded anisole-water dimer SAPT(CAS) and CC shifts are in excellent agreement, for ammonia SAPT(CAS) is only qualitatively correct.

Hydrogen bonding patterns and cooperativity in polyproline II helical bundles

Hydrogen bond cooperativity (HBC) plays an important role in stabilizing protein assemblies built by α-helices and β-sheets, the most common secondary structures. However, whether HBC exists in other types of protein secondary structures such as polyproline II (PPII) helices remains unexplored. This is intriguing, since PPII systems as assembling blocks are continuously emerging across multiple fields. Here, using a combination of computational chemistry tools and molecular modeling corroborated by experimental observables, we characterize the distinct H-bonding patterns present in PPII helical bundles and establish that HBC stabilizes intermolecular PPII helices as seen in other protein assemblies such as amyloid fibrils. In addition to cooperative interactions in canonical CO···HN H-bonds, we show that analogous interactions in non-canonical CO···HαCα H-bonds are relevant in Gly-rich PPII bundles, thus compensating for the inability of glycine residues to create hydrophobic cores. Our results provide a mechanistic explanation for the assembly of these bundles.

Photoelectron Spectroscopy and Computational Study on Microsolvated [B10H10]2- Clusters and Comparisons to Their [B12H12]2- Analogues

Microhydrated closo-boranes have attracted great interest due to their superchaotropic activity related to the well-known Hofmeister effect and important applications in biomedical and battery fields. In this work, we report a combined negative ion photoelectron spectroscopy and quantum chemical investigation on hydrated closo-decaborate clusters [B10H10]2-·nH2O (n = 1-7) with a direct comparison to their analogues [B12H12]2-·nH2O and free water clusters. A single H2O molecule is found to be sufficient to stabilize the intrinsically unstable [B10H10]2- dianion. The first two water molecules strongly interact with the solute forming B-H···H-O dihydrogen bonds while additional water molecules show substantially reduced binding energies. Unlike [B12H12]2-·nH2O possessing a highly structured water network with the attached H2O molecules arranged in a unified pattern by maximizing B-H···H-O dihydrogen bonding, distinct structural arrangements of the water clusters within [B10H10]2-·nH2O are achieved with the water cluster networks from trimer to heptamer resembling free water clusters. Such a distinct difference arises from the variations in size, symmetry, and charge distributions between these two dianions. The present finding again confirms the structural diversity of hydrogen-bonding networks in microhydrated closo-boranes and enriches our understanding of aqueous borate chemistry.

An efficient method by combining different basis sets and SAPT levels

In symmetry-adapted perturbation theory (SAPT), accurate calculations on non-covalent interaction (NCI) for large complexes with more than 50 atoms are time-consuming using large basis sets. More efficient ones with smaller basis sets usually result in poor prediction in terms of dispersion and overall energies. In this study, we propose two composite methods with baseline calculated at SAPT2/aug-cc-pVDZ and SAPT2/aug-cc-pVTZ with dispersion term corrected at SAPT2+ level using bond functions and smaller basis set with δ MP2 corrections respectively. Benchmark results on representative NCI data sets, such as S22, S66, and so forth, show significant improvements on the accuracy compared to the original SAPT Silver standard and comparable to SAPT Gold standard in some cases with much less computational cost.

A comprehensive study on three typical photoacid generators using photoelectron spectroscopy and ab initio calculations

Conducting a comprehensive molecular-level evaluation of a photoacid generator (PAG) and its subsequent impact on lithography performance can facilitate the rational design of a promising 193 nm photoresist tailored to specific requirements. In this study, we integrated spectroscopy and computational techniques to meticulously investigate the pivotal factors of three prototypical PAG anions, p-toluenesulfonate (pTS-), 2-(trifluoromethyl)benzene-1-sulfonate (TFMBS-), and triflate (TF-), in the lithography process. Our findings reveal a significant redshift in the absorption spectra caused by specific PAG anions, attributed to their involvement in electronic transition processes, thereby enhancing the transparency of the standard PAG cation, triphenylsulfonium (TPS+), particularly at ∼193 nm. Furthermore, the electronic stability of PAG anions can be enhanced by solvent effects with varying degrees of strength. We observed the lowest vertical detachment energy of 6.6 eV of pTS- in PGMEA solution based on the polarizable continuum model, which prevents anion loss at 193 nm lithography. In addition, our findings indicate gas-phase proton affinity values of 316.4 kcal/mol for pTS-, 308.1 kcal/mol for TFMBS-, and 303.2 kcal/mol for TF-, which suggest the increasing acidity strength, yet even the weakest acid pTS- is still stronger than strong acid HBr. The photolysis of TPS+-based PAG, TPS+·pTS-, generated an excited state leading to homolysis bond cleavage with the lowest reaction energy of 83 kcal/mol. Overall, the PAG anion pTS- displayed moderate acidity, possessed the lowest photolysis reaction energy, and demonstrated an appropriate redshift. These properties collectively render it a promising candidate for an effective acid producer.

A first-principles alternative to empirical solvent parameters

The use of solvents is ubiquitous in chemistry. Empirical parameters, such as the Kamlet-Taft parameters and Gutmann donor/acceptor numbers, have long been used to predict and quantify the effects solvents have on chemical phenomena. Collectively however, such parameters are unsatisfactory, since each describes ultimately the same non-covalent solute-solvent and solute-solute interactions in completely disparate ways. Here we hypothesise that empirical solvent parameters are essentially proxy measures of the electrostatic terms that dominate solvent-solute interactions. On the basis of this hypothesis, we develop a new fundamental descriptor of these interactions, , and show that it is a self-consistent, probe-free, first principles alternative to established empirical solvent parameters.

Quantum Mechanics Characterization of Non-Covalent Interaction in Nucleotide Fragments

Accurate calculation of non-covalent interaction energies in nucleotides is crucial for understanding the driving forces governing nucleic acid structure and function, as well as developing advanced molecular mechanics forcefields or machine learning potentials tailored to nucleic acids. Here, we dissect the nucleotides' structure into three main constituents: nucleobases (A, G, C, T, and U), sugar moieties (ribose and deoxyribose), and phosphate group. The interactions among these fragments and between fragments and water were analyzed. Different quantum mechanical methods were compared for their accuracy in capturing the interaction energy. The non-covalent interaction energy was decomposed into electrostatics, exchange-repulsion, dispersion, and induction using two ab initio methods: Symmetry-Adapted Perturbation Theory (SAPT) and Absolutely Localized Molecular Orbitals (ALMO). These calculations provide a benchmark for different QM methods, in addition to providing a valuable understanding of the roles of various intermolecular forces in hydrogen bonding and aromatic stacking. With SAPT, a higher theory level and/or larger basis set did not necessarily give more accuracy. It is hard to know which combination would be best for a given system. In contrast, ALMO EDA2 did not show dependence on theory level or basis set; additionally, it is faster.

In-depth exploration of catalytic sites on amorphous solid water: I. The astrosynthesis of aminomethanol

Chemical processes taking place on ice-grain mantles are pivotal to the complex chemistry of interstellar environments. In this study, we conducted a comprehensive analysis of the catalytic effects of an amorphous solid water (ASW) surface on the reaction between ammonia (NH3) and formaldehyde (H2CO) to form aminomethanol (NH2CH2OH) using density functional theory. We identified potential catalytic sites based on the binding energy distribution of NH3 and H2CO reactants, on a set-of-clusters surface model composed of 22 water molecules and found a total of 14 reaction paths. Our results indicate that the catalytic sites can be categorized into four groups, depending on the interactions of the carbonyl oxygen and the amino group with the ice surface in the reactant complex. A detailed analysis of the reaction mechanism using Intrinsic Reaction Coordinate and reaction force analysis, revealed three distinct chemical events for this reaction: formation of the C-N bond, breaking of the N-H bond, and formation of the O-H hydroxyl bond. Depending on the type of catalytic site, these events can occur within a single, concerted, albeit asynchronous, step, or can be isolated in a step-wise mechanism, with the lowest overall transition state energy observed at 1.3 kcal mol-1. A key requirement for the low-energy mechanism is the presence of a pair of dangling OH bonds on the surface, found at 5% of the potential catalytic sites on an ASW porous surface.

Simulation of Linear and Cyclic Alkanes with Second-Order Møller-Plesset Perturbation Theory through Adaptive Force Matching

Predicting ensemble properties, such as density and heat of vaporization, of small hydrocarbons is challenging due to the dispersion-dominated weak interactions between these molecules. With the adaptive force matching (AFM) method, the bonded and short-range nonbonded interactions are fitted to second-order Møller-Plesset perturbation theory (MP2) references computed with the def2-TZVP basis set. The dispersion is modeled using symmetry adapted perturbation theory (SAPT) at MP4 accuracy using the def2-TZVPD basis set. A new charge matrix decomposition technique is described to obtain partial charges in AFM. Although the models developed do not have any empirical parameters, several properties of the resulting models are compared with experiments as validations. The density, heat of vaporization, pressure dependence of density, diffusion constants, and surface tensions all show quantitative agreement with experiments. Although the density shows a very small systematic error, which could be due to missing three-body dispersion, the heat of vaporization agrees with experiments of within 0.5%. The paper shows that AFM can be used as a reliable tool to enable simulations at post-Hartree-Fock quality at the cost of molecular mechanics force fields.

Reducing the Runtime of Fault-Tolerant Quantum Simulations in Chemistry through Symmetry-Compressed Double Factorization

Quantum phase estimation based on qubitization is the state-of-the-art fault-tolerant quantum algorithm for computing ground-state energies in chemical applications. In this context, the 1-norm of the Hamiltonian plays a fundamental role in determining the total number of required iterations and also the overall computational cost. In this work, we introduce the symmetry-compressed double factorization (SCDF) approach, which combines a CDF of the Hamiltonian with the symmetry shift technique, significantly reducing the 1-norm value. The effectiveness of this approach is demonstrated numerically by considering various benchmark systems, including the FeMoco molecule, cytochrome P450, and hydrogen chains of different sizes. To compare the efficiency of SCDF to other methods in absolute terms, we estimate Toffoli gate requirements, which dominate the execution time on fault-tolerant quantum computers. For the systems considered here, SCDF leads to a sizable reduction of the Toffoli gate count in comparison to other variants of DF or even tensor hypercontraction, which is usually regarded as the most efficient approach for qubitization.

Hexakispyrazolylethane: New Strategy for Stabilization of Hexaarylethane

Hexakis(4-trimethylsilylpyrazol-1-yl)ethane was synthesized by the oxidative dimerization of tris(4-trimethylsilylpyrazol-1-yl)methane. Single-crystal X-ray structural analysis of hexakis(4-trimethylsilylpyrazol-1-yl)ethane showed that the ethane C-C bond (1.623(4) Å) is shorter than that in hexaphenylethane (1.67(3) Å). In solution, hexakis(4-trimethylsilylpyrazol-1-yl)ethane existed as a single species, contrastive that conventional hexaphenylethanes can keep the central C-C bond only by the aid of additional bridges between the two triarylmethyl units. Theoretical calculations indicated that the tris(pyrazol-1-yl)methyl radical, which is anticipated to be under equilibrium with hexakis(pyrazol-1-yl)ethane, is less stable than trityl radicals due to lack of delocalization of the radicals. Furthermore, introduction of pyrazole groups allowed additional bridging between the two triarylmethyl moieties through metal coordination to the adjacent N atoms: hexakis(4-trimethylsilylpyrazol-1-yl)ethane exhibited unique coordination to three Ag atoms affording a hexaarylethane analog bearing three N-Ag-N bridges.

Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Revisiting conventional noncovalent interactions towards a complete understanding: from tetrel to pnicogen, chalcogen, and halogen bond

Typical noncovalent interactions, including tetrel (TtB), pnicogen (PniB), chalcogen (ChalB), and halogen bonds (HalB), were systematically re-investigated by modeling the N⋯Z interactions (Z = Si, P, S, Cl) between NH3 - as a nucleophilic, and SiF4, PF3, SF2, and ClF - as electrophilic components, employing highly reliable ab initio methods. The characteristics of N⋯Z interactions when Z goes from Si to Cl, were examined through their changes in stability, vibrational spectroscopy, electron density, and natural orbital analyses. The binding energies of these complexes at CCSD(T)/CBS indicate that NH3 tends to hold tightly most with ClF (-34.7 kJ mol-1) and SiF4 (-23.7 kJ mol-1) to form N⋯Cl HalB and N⋯Si TtB, respectively. Remarkably, the interaction energies obtained from various approaches imply that the strength of these noncovalent interactions follows the order: N⋯Si TtB > N⋯Cl HalB > N⋯S ChalB > N⋯P PniB, that differs the order of their corresponding complex stability. The conventional N⋯Z noncovalent interactions are characterized by the local vibrational frequencies of 351, 126, 167, and 261 cm-1 for TtB, PniB, ChalB, and HalB, respectively. The SAPT2+(3)dMP2 calculations demonstrate that the primary force controlling their strength retains the electrostatic term. Accompanied by the stronger strength of N⋯Si TtB and N⋯Cl HalB, the AIM and NBO results state that they are partly covalent in nature with amounts of 18.57% and 27.53%, respectively. Among various analysis approaches, the force constant of the local N⋯Z stretching vibration is shown to be most accurate in describing the noncovalent interactions.

Locking water molecules via ternary O-H⋯O intramolecular hydrogen bonds in perhydroxylated closo-dodecaborate

A multitude of applications related to perhydroxylated closo-dodecaborate B12(OH)122- in the condensed phase are inseparable from the fundamental mechanisms underlying the high water orientation selectivity based on the base B12(OH)122-. Herein, we directly compare the structural evolution of water clusters, ranging from monomer to hexamer, oriented by functional groups in the bases B12H122-, B12H11OH2- and B12(OH)122- using multiple theoretical methods. A significant revelation is made regarding B12(OH)122-: each additional water molecule is locked into the intramolecular hydrogen bond B-O-H ternary ring in an embedded form. This new pattern of water cluster growth suggests that B-(H-O)⋯H-O interactions prevail over the competition from water-hydrogen bonds (O⋯H-O), distinguishing it from the behavior observed in B12H122- and B12H11OH2- bases, in which competition arises from a mixed competing model involving dihydrogen bonds (B-H⋯H-O), conventional hydrogen bonds (B-(H-O)⋯H-O) and water hydrogen bonds (O⋯H-O). Through aqueous solvation and ab initio molecular dynamics analysis, we further demonstrate the largest water clusters in the first hydrated shell with exceptional thermodynamic stability around B12(OH)122-. These findings provide a solid scientific foundation for the design of boron cluster chemistry incorporating hydroxyl-group-modified borate salts with potential implications for various applications.

Experimental Quantification of Halogen⋅⋅⋅Arene van der Waals Contacts

Crystallographic and computational studies suggest the occurrence of favourable interactions between polarizable arenes and halogen atoms. However, the systematic experimental quantification of halogen⋅⋅⋅arene interactions in solution has been hindered by the large variance in the steric demands of the halogens. Here we have synthesized molecular balances to quantify halogen⋅⋅⋅arene contacts in 17 solvents and solvent mixtures using 1 H NMR spectroscopy. Calculations indicate that favourable halogen⋅⋅⋅arene interactions arise from London dispersion in the gas phase. In contrast, comparison of our experimental measurements with partitioned SAPT0 energies indicate that dispersion is sufficiently attenuated by the solvent that the halogen⋅⋅⋅arene interaction trend was instead aligned with increasing exchange repulsion as the halogen increased in size (ΔGX ⋅⋅⋅Ph =0 to +1.5 kJ mol-1 ). Halogen⋅⋅⋅arene contacts were slightly less disfavoured in solvents with higher solvophobicities and lower polarizabilities, but strikingly, were always less favoured than CH3 ⋅⋅⋅arene contacts (ΔGMe ⋅⋅⋅Ph =0 to -1.4 kJ mol-1 ).

Systematic Study of Heteroarene Stacking Using a Congeneric Set of Molecular Glues for Procaspase-6

Heteroaromatic stacking interactions are important in drug binding, supramolecular chemistry, and materials science, making protein-ligand model systems of these interactions of considerable interest. Here we studied 30 congeneric ligands that each present a distinct heteroarene for stacking between tyrosine residues at the dimer interface of procaspase-6. Complex X-ray crystal structures of 10 analogs showed that stacking geometries were well conserved, while high-accuracy computations showed that heteroarene stacking energy was well correlated with predicted overall ligand binding energies. Empirically determined KD values in this system thus provide a useful measure of heteroarene stacking with tyrosine. Stacking energies are discussed in the context of torsional strain, the number and positioning of heteroatoms, tautomeric state, and coaxial orientation of heteroarene in the stack. Overall, this study provides an extensive data set of empirical and high-level computed binding energies in a versatile new protein-ligand system amenable to studies of other intermolecular interactions.

Noncovalent Interactions in Hydrated Nitrosonium Ion Clusters Mediated by Hydrogen-Bonded Water Networks

As important species in the D region of the ionosphere, hydrated nitrosonium ion clusters [NO⁺(H₂O)_n] are also archetypal and concise models to illustrate effects of different solvent shells. We have investigated noncovalent interactions in NO⁺(H₂O)₃ and NO⁺(H₂O)₄ isomers with high levels of ab initio and symmetry-adapted perturbation theory (SAPT) methods. On the basis of our computations, the exchange energies become much more repulsive, whereas the induction energies are significantly more attractive for the noncovalent interactions of NO⁺ with hydrogen-bonded water chains. Combined with analyses of the electron densities for the NO⁺(H₂O)₃ and NO⁺(H₂O)₄ isomers, we propose that the counteracting effect of the exchange and induction energies could be deemed as an index for the tendency to form the HO-NO covalent bond. Moreover, we have also found that the third-order induction terms are very important to evaluate reasonable charge transfer energies with the SAPT computations.

Transition Metal Sensing with Nitrogenated Holey Graphene: A First-Principles Investigation

The toxicity of transition metals, including copper(II), manganese(II), iron(II), zinc(II), hexavalent chromium, and cobalt(II), at elevated concentrations presents a significant threat to living organisms. Thus, the development of efficient sensors capable of detecting these metals is of utmost importance. This study explores the utilization of two-dimensional nitrogenated holey graphene (C₂N) nanosheet as a sensor for toxic transition metals. The C₂N nanosheet's periodic shape and standard pore size render it well suited for adsorbing transition metals. The interaction energies between transition metals and C₂N nanosheets were calculated in both gas and solvent phases and were found to primarily result from physisorption, except for manganese and iron which exhibited chemisorption. To assess the interactions, we employed NCI, SAPT0, and QTAIM analyses, as well as FMO and NBO analysis, to examine the electronic properties of the TM@C₂N system. Our results indicated that the adsorption of copper and chromium significantly reduced the HOMO-LUMO energy gap of C₂N and significantly increased its electrical conductivity, confirming the high sensitivity of C₂N towards copper and chromium. The sensitivity test further confirmed the superior sensitivity and selectivity of C₂N towards copper. These findings offer valuable insight into the design and development of sensors for the detection of toxic transition metals.

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

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

Static and Dynamical Quantum Studies of CX3-AlX2 and CSiX3-BX2 (X = F, Cl, Br) Complexes with Hydrocyanic Acid: Unusual Behavior of Strong π-Hole at Triel Center

The set of TX₃-TrX₂ (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX₃-AlX₂) and least (SiX₃-BX₂) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a "residual π-hole" site that was able to drive second ligand attachment, giving rise to the two "unusual trimers" examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the "unusual trimers". Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.

Isomer Stability Dependence of Hydrogen-Bonded Benzoxazole Clusters on Solvent Molecules

The structures of hydrogen-bonded benzoxazole clusters with methanol and ammonia, BO-(CH₃OH)_n (n = 1-3) and BO-(NH₃)_n (n = 1, 2), in a supersonic jet have been investigated by measuring the S₁-S₀ electronic spectra and isomer-selected vibrational spectra with the aid of quantum chemical calculations. Similar to BO-(H₂O)₁, two isomers of BO-(NH₃)₁ were observed, which form two types of hydrogen bond networks starting from the CH bond at the 2-/7-position to the nitrogen atom of BO (C2HN/C7HN). The relative stability of these isomers strongly depends on solvent molecules. Natural bond orbital analysis reveals that the OH···N hydrogen bond is dominant in BO-(H₂O)₁ and that intermolecular interaction between the CH group and the nitrogen atom of ammonia, especially C2H···N, is significantly enhanced, resulting in a more stable C2HN isomer. Symmetry-adapted perturbation theory calculations indicate that the dispersion interaction between the methyl group of methanol and π electron cloud on the BO ring is responsible for the extreme stability of the C7HN of BO-(CH₃OH)₁. Furthermore, using time-dependent density functional theory calculations, the isomer tendency of the electronic transition shifts from the monomer origin is reproduced and it is proposed that the significant blue shift in C2HN is due to the shortened C2H bond length upon electronic excitation.

Importance of Noncovalent Interactions Involving Sulfur Atoms in Thiopeptide Antibiotics─Glycothiohexide α and Nocathiacin I

Noncovalent interactions involving sulfur atoms play essential roles in protein structure and function by significantly contributing to protein stability, folding, and biological activity. Sulfur is a highly polarizable atom that can participate in many types of noncovalent interactions including hydrogen bonding, sulfur-π interactions, and S-lone pair interactions, but the impact of these sulfur-based interactions on molecular recognition and drug design is still often underappreciated. Here, we examine, using quantum chemical calculations, the roles of sulfur-based noncovalent interactions in complex naturally occurring molecules representative of thiopeptide antibiotics: glycothiohexide α and its close structural analogue nocathiacin I. While donor-acceptor orbital interactions make only very small contributions, electrostatic and dispersion contributions are predicted to be significant in many cases. In pursuit of understanding the magnitudes and nature of these noncovalent interactions, we made potential structural modifications that could significantly expand the chemical space of effective thiopeptide antibiotics.

The Ultrashort Spike-Ring Interaction in Substituted Iron Maiden Molecules

The in forms of molecular iron maidens are known for their unique ultrashort interaction between the apical hydrogen atom or its small substituent and the surface of the benzene ring. It is generally believed that this forced ultrashort X⋯π contact is associated with high steric hindrance, which is responsible for specific properties of iron maiden molecules. The main aim of this article is to investigate the influence of significant charge enrichment or depletion of the benzene ring on the characteristics of the ultrashort C-X⋯π contact in iron maiden molecules. For this purpose, three strongly electron-donating (-NH2) or strongly electron-withdrawing (-CN) groups were inserted into the benzene ring of in-[34,10][7]metacyclophane and its halogenated (X = F, Cl, Br) derivatives. It is shown that, despite such extremely electron-donating or electron-accepting properties, the considered iron maiden molecules surprisingly reveal quite high resistance to changes in electronic properties.

Highly Structured Water Networks in Microhydrated Dodecaborate Clusters

We report a combined photoelectron spectroscopy and theoretical investigation of a series of size-selected hydrated closo-dodecaborate clusters B₁₂X₁₂^2-·nH₂O (X = H, F, or I; n = 1-6). Distinct structural arrangements of water clusters from monomer to hexamer can be achieved by using different B₁₂X₁₂^2- bases, illustrating the evident solute specificity. Because B-H···H-O dihydrogen bonds are stronger than O···H-O hydrogen bonds in water, the added water molecules are arranged in a unified binding mode by forming highly structured water networks manipulated by B₁₂H₁₂^2-. As a comparison, the hydrated B₁₂F₁₂^2- clusters display similar water evolution for n values of 1 and 2 but different binding modes for larger clusters, while water networks in B₁₂I₁₂^2- share similarities with the free water clusters. This finding provides a consistent picture of the structural diversity of hydrogen bonding networks in microhydrated dodecaborates and a molecular-level understanding of microsolvation dynamics in aqueous borate chemistry.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r^-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

A Theoretical Study of the Halogen Bond between Heteronuclear Halogen and Benzene

Halogen bonds play an important role in many fields, such as biological systems, drug design and crystal engineering. In this work, the structural characteristics of the halogen bond between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene were studied using density functional theory. The structures of the complexes between heteronuclear halogen and benzene have Cs symmetry. The interaction energies of the complexes between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene range from -27.80 to -37.18 kJ/mol, increasing with the increases in the polarity between the atoms of X and D, and are proportional to the angles of a between the Z axis and the covalent bond of heteronuclear halogen. The electron density (ρ) and corresponding Laplacian (∇²ρ) values indicate that the interaction of the heteronuclear halogen and benzene is a typical long-range weak interaction similar to a hydrogen bond. Independent gradient model analysis suggests that the van der Waals is the main interaction between the complexes of heteronuclear halogen and benzene. Symmetry-adapted perturbation theory analysis suggests that the electrostatic interaction is the dominant part in the complexes of C₆H₆⋯ClF, C₆H₆⋯ICl, C₆H₆⋯BrF and C₆H₆⋯IF, and the dispersion interaction is the main part in the complexes of C₆H₆⋯BrCl, C₆H₆⋯IBr.

Structures of the (Imidazole)nH+ ... Ar (n=1,2,3) complexes determined from IR spectroscopy and quantum chemical calculations

Here, we present new cryogenic infrared spectra of the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + (n=1,2,3) ions. The data was obtained using helium tagging infrared predissociation spectroscopy. The new results were compared with the data obtained by Gerardi et al. (Chem. Phys. Lett. 501:172–178, 2011) using the same technique but with argon as a tag. Comparison of the two experiments, assisted by theoretical calculations, allowed us to evaluate the preferable attachment positions of argon to the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + frame. Argon attaches to nitrogen-bonded hydrogen in the case of the (Imidazole)H$$^+$$ + ion, while in (Imidazole)$$_{2}\mathrm{H}^{+}$$ 2 H + and (Imidazole)$$_{3}\mathrm{H}^{+}$$ 3 H + the preferred docking sites for the argon are in the center of the complex. This conclusion is supported by analyzing the spectral features attributed to the N–H stretching vibrations. Symmetry adapted perturbation theory (SAPT) analysis of the non-covalent forces between argon and the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + (n=1,2,3) frame revealed that this switch of docking preference with increasing complex size is caused by an interplay between induction and dispersion interactions.

Mechanism and Origin of Remote Stereocontrol in the Organocatalytic Enantioselective Formal C(sp2)–H Alkylation Using Nitroalkanes as Alkylating Agents

Experimental ¹³C kinetic isotope effects (KIEs) and density functional theory (DFT) calculations are used to evaluate the mechanism and origin of enantioselectivity in the formal C(sp²)-H alkylative desymmetrization of cyclopentene-1,3-diones using nitroalkanes as the alkylating agent. An unusual combination of an inverse (∼0.980) and a normal (∼1.033) KIE is observed on the bond-forming carbon atoms of the cyclopentene-1,3-dione and nitroalkane, respectively. These data provide strong support for a mechanism involving reversible carbon-carbon bond formation followed by rate- and enantioselectivity-determining nitro group elimination. The theoretical free-energy profile and the predicted KIEs indicate that this elimination event occurs via an E1cB pathway. The origin of remote stereocontrol is evaluated by distortion-interaction and SAPT0 analyses of the E1cB transition states leading to both enantiomers.

Water: new aspect of hydrogen bonding in the solid state

All water-water contacts in the crystal structures from the Cambridge Structural Database with d OO ≤ 4.0 Å have been found. These contacts were analysed on the basis of their geometries and interaction energies from CCSD(T)/CBS calculations. The results show 6729 attractive water-water contacts, of which 4717 are classical hydrogen bonds (d OH ≤ 3.0 Å and α ≥ 120°) with most being stronger than -3.3 kcal mol-1. Beyond the region of these hydrogen bonds, there is a large number of attractive interactions (2062). The majority are antiparallel dipolar interactions, where the O-H bonds of two water molecules lying in parallel planes are oriented antiparallel to each other. Developing geometric criteria for these antiparallel dipoles (β1, β2 ≥ 160°, 80 ≤ α ≤ 140° and T HOHO > 40°) yielded 1282 attractive contacts. The interaction energies of these antiparallel oriented water molecules are up to -4.7 kcal mol-1, while most of the contacts have interaction energies in the range -0.9 to -2.1 kcal mol-1. This study suggests that the geometric criteria for defining attractive water-water interactions should be broader than the classical hydrogen-bonding criteria, a change that may reveal undiscovered and unappreciated interactions controlling molecular structure and chemistry.

Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C60 Dimer, DNA-Ellipticine, and HIV-Indinavir

In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C₆₀ dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C₆₀. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r²SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.

Origin of the Microsolvation Effect on the Central Barriers of SN2 Reactions

We have quantitatively analyzed the microsolvation effect on the central barriers of microsolvated bimolecular nucleophilic substitution (S_N2) reactions by means of a two-step energy decomposition procedure. According to the first energy decompositions, an obvious increase in the central barrier for a microsolvated S_N2 reaction against its unsolvated counterpart can be mainly ascribed to the fact that the interaction between the solute and the conjunct solvent becomes less attractive from the reactant complex to the transition state. On the basis of the second energy decompositions with symmetry-adapted perturbation theory, this less attractive interaction in the transition state is primarily due to the interplay of the changes in the electrostatic, exchange, and induction components. However, the contribution of the change for the dispersion component is relatively small. A distinct linear correlation has also been observed between the changes of the total interaction energies and those of the corresponding electrostatic components for the microsolvated S_N2 reactions studied in this work. Moreover, the two-step energy decomposition procedure employed in this work is expected to be extensively applied to the gas phase reactions mediated by molecules or clusters.

Adsorption of Industrial Gases (CH4, CO2, and CO) on Olympicene: A DFT and CCSD(T) Investigation

Olympicene C₁₉H₁₂, an organic semiconductor, is investigated as an adsorption material for toxic industrial gas molecules such as CH₄, CO₂, and CO. A deep insight of complexation of CH₄, CO₂, and CO with olympicene (analyte@OLY) was obtained by interaction energy, symmetry-adopted perturbation theory (SAPT2+), quantum theory of atoms in molecules (QTAIM), density of states (DOS), noncovalent interaction (NCI), and frontier molecular orbital and natural bond orbital analysis. Domain-based local pair natural orbital coupled cluster theory single-point energy calculations were performed using the cc-pVTZ basis set in combination with corresponding auxiliary cc-pVTZ/JK and cc-pVTZ/C basis sets. For all property calculations of doped olympicene complexes, the ωB97M-V functional was employed. The stability trend for interaction energies is CO₂@OLY > CH₄@OLY > CO@OLY. QTAIM and NCI analysis confirmed the presence of NCIs, where the dispersion factor (in CH₄@OLY) has the highest contribution, as revealed from SAPT2+. The chemical sensitivity of the system was evidenced by the origination of new energy states in DOS spectra. The recovery time for the analyte@OLY complex was calculated at 300 K, and an excellent recovery response was observed. All results evidently indicated weak interactions of the olympicene surface with CH₄, CO₂, and CO.

External Electric Field Effect on the Strength of σ-Hole Interactions: A Theoretical Perspective in Like⋯Like Carbon-Containing Complexes

For the first time, σ-hole interactions within like⋯like carbon-containing complexes were investigated, in both the absence and presence of the external electric field (EEF). The effects of the directionality and strength of the utilized EEF were thoroughly unveiled in the (F-C-F₃)₂, (F-C-H₃)₂, and (H-C-F₃)₂ complexes. In the absence of the EEF, favorable interaction energies, with negative values, are denoted for the (F-C-F₃)₂ and (H-C-F₃)₂ complexes, whereas the (F-C-H₃)₂ complex exhibits unfavorable interactions. Remarkably, the strength of the applied EEF exhibits a prominent role in turning the repulsive forces within the latter complex into attractive ones. The symmetrical nature of the considered like⋯like carbon-containing complexes eradicated the effect of directionality of the EEF. The quantum theory of atoms in molecules (QTAIM), and the noncovalent interaction (NCI) index, ensured the occurrence of the attractive forces, and also outlined the substantial contributions of the three coplanar atoms to the total strength of the studied complexes. Symmetry-adapted perturbation theory (SAPT) results show the dispersion-driven nature of the interactions.

Radial Potential Energy Functions of Linear Halogen-Bonded Complexes YX···ClF (YX = FB, OC, SC, N2) and the Effects of Substituting X by Second-Row Analogues: Mulliken Inner and Outer Complexes

Energies of linear, halogen-bonded complexes in the isoelectronic series YX···ClF (YX = FB, OC, or N₂) are calculated at several levels of theory as a function of the intermolecular distance r(X···Cl) to yield radial potential energy functions. When YX = OC, a secondary minimum is observed corresponding to lengthened and shortened distances r(ClF) and r(CCl), respectively, relative to the primary minimum, suggesting a significant contribution from the Mulliken inner complex structure [O=C–Cl]⁺···F^–. A conventional weak, halogen-bond complex OC···ClF occurs at the primary minimum. For YX = FB, the primary minimum corresponds to the inner complex [F=B–Cl]⁺···F^–, while the outer complex FB···ClF is at the secondary minimum. The effects on the potential energy function of systematic substitution of Y and X by second-row congeners and of reversing the order of X and Y are also investigated. Symmetry-adapted perturbation theory and natural population analyses are applied to further understand the nature of the various halogen-bond interactions.

σ-Hole and LP-Hole Interactions of Pnicogen···Pnicogen Homodimers under the External Electric Field Effect: A Quantum Mechanical Study

σ-Hole and lone-pair (lp)-hole interactions within σ-hole···σ-hole, σ-hole···lp-hole, and lp-hole···lp-hole configurations were comparatively investigated on the pnicogen···pnicogen homodimers (PCl3)2, for the first time, under field-free conditions and the influence of the external electric field (EEF). The electrostatic potential calculations emphasized the impressive versatility of the examined PCl3 monomers to participate in σ-hole and lp-hole pnicogen interactions. Crucially, the sizes of σ-hole and lp-hole were enlarged under the influence of the positively directed EEF and decreased in the case of reverse direction. Interestingly, the energetic quantities unveiled more favorability of the σ-hole···lp-hole configuration of the pnicogen···pnicogen homodimers, with significant negative interaction energies, than σ-hole···σ-hole and lp-hole···lp-hole configurations. Quantum theory of atoms in molecules and noncovalent interaction index analyses were adopted to elucidate the nature and origin of the considered interactions, ensuring their closed shell nature and the occurrence of attractive forces within the studied homodimers. Symmetry-adapted perturbation theory-based energy decomposition analysis alluded to the dispersion force as the main physical component beyond the occurrence of the examined interactions. The obtained findings would be considered as a fundamental underpinning for forthcoming studies pertinent to chemistry, materials science, and crystal engineering.

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car-Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion-the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.