Density Functional Theory and Hydrogen Bonds: Are We There Yet?

The Peculiar H-Bonding Network of 4-Methylcatechol: A Coupled Diffraction and In Silico Study

The crystal structure of 4-methylcatechol (4MEC) has, to date, never been solved, despite its very simple chemical formula C7O2H8 and the many possible applications envisaged for this molecule. In this work, this gap is filled and the structure of 4MEC is obtained by combining X-ray powder diffraction and first principle calculations to carefully locate hydrogen atoms. Two molecules are present in the asymmetric unit. Hirshfeld analysis confirmed the reliability of the solved structure, since the two molecules show rather different environments and H-bond interactions of different directionality and strength. The packing is characterised by a peculiar hydrogen bond network with hydroxyl nests formed by two adjacent octagonal frameworks. It is noteworthy that the observed short contacts suggest strong inter-molecular interactions, further confirmed by strong inter-crystalline aggregation observed by microscopic images, indicating the growth, in many crystallization attempts, of single aggregates taller than half a centimetre and, often, with spherical shapes. These peculiarities are induced by the presence of methyl group in 4MEC, since the parent compound catechol, despite its chemical similarity, shows a standard layered packing alternating hydrophobic and polar layers. Finally, the complexity and peculiarity of the packing and crystal growth features explain why a single crystal could not be obtained for a standard structural analysis.

Tetrel Bonding of the Carbenium Ion Forms a Pentacoordinate Carbon Atom

As a flat trigonal species, the CR3 + carbenium ion contains a pair of deep π-holes above and below its molecular plane. In the case of CH3 + a first base will form a covalent bond with the central C, making the combined species tetrahedral. Approach of a second base to the opposite side results in a longer but rather strong noncovalent tetrel bond (TB). While CMe3 + can also form a similar asymmetric complex with a pair of bases, it also has the capacity to form a pair of nearly equivalent TBs, such that the resulting symmetric trigonal bipyramid configuration is only slightly higher in energy. When the three substituents on the central C are phenyl rings, the symmetric configuration with two TBs predominates. These tetrel bonds are quite strong, reaching up to 20 kcal/mol. Adding OPH2 or OCH substituents to the phenyl rings permits the formation of intramolecular C⋅⋅O TBs to the central C, very similar in many respects to the case where these TBs are intermolecular.

Types of noncovalent bonds within complexes of thiazole with CF4 and SiF4

The five-membered heteroaromatic thiazole molecule contains a number of electron-rich regions that could attract an electrophile, namely the N and S lone pairs that lie in the molecular plane, and π-system areas above the plane. The possibility of each of these sites engaging in a tetrel bond (TB) with CF4 and SiF4, as well as geometries that encompass a CH⋯F H-bond, was explored via DFT calculations. There are a number of minima that occur in the pairing of thiazole with CF4 that are very close in energy, but these complexes are weakly bound by less than 2 kcal mol-1 and the presence of a true TB is questionable. The inclusion of zero-point vibrational energies alters the energetic ordering, which is further modified when entropic effects are added. The preferred geometry would thus be sensitive to the temperature of an experiment. Replacement of CF4 by SiF4 leaves intact most of the configurations, and their tight energetic clustering, the ordering of which is again altered as the temperature rises. But there is one exception in that by far the most tightly bound complex involves a strong Si⋯N TB between SiF4 and the lone pair of the thiazole N, with an interaction energy of 30 kcal mol-1. Even accounting for its high deformation energy and entropic considerations, this structure remains as clearly the most stable at any temperature.

Wolfium bonds in homodimers of MX4Y (M = Mo, W; X = F, Cl, Br; Y = O, S, Se)

The term "wolfium bond" has been recently introduced to describe the noncovalent attraction between an atom of group 6 and a nucleophile via a σ-hole binding site. Crystal structures commonly contain a motif wherein two MX4Y units are arranged in close proximity, where M represents either Mo or W, and X and Y refer to halogen and chalcogen atoms respectively. DFT calculations were thus applied to a wide range of homodimers of these molecules so as to assess their preferred arrangements, and to characterize the types of bonding that are present in each in a systematic manner. The most stable Dual-X configuration is symmetric and contains a pair of equivalent M⋯X bonds. The interaction energies range from -8 to -29 kcal mol-1, and are largest for X = F, Y = O, and M = W. The X electron donor is replaced by Y, and the two wolfium bonds are reduced to one, in the less stable Mono-Y structure, with interaction energies between -2 and -10 kcal mol-1. There is some question as to whether the weaker bonds of this type constitute true wolfium bonds.

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.

Influence of Internal Angular Arrangement on Pnicogen Bond Strength

The three Z-X covalent bonds of a ZX3 unit (Z = P, As, Sb, Bi) are normally arranged in a pyramidal structure. Quantum chemical calculations show that pnicogen bonds (ZBs) to the central Z are weakened if ZX3 is flattened, as in the opening of an umbrella. The partial closing of the umbrella has the opposite effect of substantially strengthening these ZBs, even amounting to a 2- or 3-fold magnification in certain cases. The strongest such bonds, wherein Sb and Bi are in a strained configuration within a ZO3CH model system, have interaction energies of 20 kcal/mol with an NH3 base. Most of these systems, whether flattened or more pyramidal, are capable of engaging in three ZBs simultaneously, despite a certain amount of negative cooperativity.

Transition between the Noncovalency and Covalency of σ-Hole Bonds

The properties of the bond between a N-ligand and a Lewis acid containing a σ-hole are studied by quantum chemical methods. Interactions considered include pnicogen bonds involving SbX5, PX5, and PX3, where X represents any of the halogen atoms F, Cl, Br, or I. Also studied are the tetrel bonds of PbX4 and SiX4, as well as the chalcogen bond involving TeOX4. Both NH3 and NCH are applied as two possible bases of differing potency. Some of the bonds are very strong with interaction energies easily exceeding 25 kcal/mol and with AIM bond critical point densities much higher than 0.04 au, suggesting their classification as coordinate covalent bonds. The pentavalent SbX5 and PX5 fall into this category when combined with NH3, as does TeOX4. Although the tetrel bonds involving PbX4 are only slightly weaker, they are probably better viewed as a strong noncovalent bond on the cusp of covalency. Changing the internal bonding of hypervalent SbX5 to the more conventional SbX3 weakens the interaction to a classical noncovalent pnicogen bond. Reducing the base nucleophilicity from NH3 to NCH weakens the bonds so that they are clearly noncovalent.

Concerted versus stepwise proton transfer reactions in the [2, 2'-bipyridyl]-3-3'-diol molecule: A static and dynamic ab-initio investigation

The double proton transfer (PT) reaction has been investigated in the [2,2'-bipyridyl]-3-3'-diol, a complex molecule where the proton movements is coupled to significant rearrangement of the electronic structure. Moreover, the reaction could be concerted, that is the two protons are exchanged simultaneously, or stepwise, where the two protons are transferred sequentially. To this end, a static exploration of the potential energy surface (PES) was carried together with the analysis of the free-energy surface (FES), both surfaces being evaluated at density functional theory level and different exchange-correlation functionals. While the concerted mechanism has been clearly discharged, the characteristics of the stepwise PT significantly depends on the chosen functionals, some suggesting a clear stepwise mechanism characterized by a stable reaction intermediates and two transitions states, whereas other approaches propend for a asynchronous PT, with a single TS. These features appear on both PES and FES, albeit some differences appears due to their different nature.

Relation between Halogen Bond Strength and IR and NMR Spectroscopic Markers

The relationship between the strength of a halogen bond (XB) and various IR and NMR spectroscopic quantities is assessed through DFT calculations. Three different Lewis acids place a Br or I atom on a phenyl ring; each is paired with a collection of N and O bases of varying electron donor power. The weakest of the XBs display a C-X bond contraction coupled with a blue shift in the associated frequency, whereas the reverse trends occur for the stronger bonds. The best correlations with the XB interaction energy are observed with the NMR shielding of the C atom directly bonded to X and the coupling constants involving the C-X bond and the C-H/F bond that lies ortho to the X substituent, but these correlations are not accurate enough for the quantitative assessment of energy. These correlations tend to improve as the Lewis acid becomes more potent, which makes for a wider range of XB strengths.

DFT Model of Elemental Sulfur in Sulfolane Solutions

The structure and the thermodynamic and optical (UV) properties of elemental sulfur solution in sulfolane (Sl) have been studied using density functional theory methods. The cyclic molecular form of sulfur (S8 "crown") was found using PBE1PBE/6-311+G(d,p) approximation in combination with a polarizable continuum model (the integral equation formalism variant) to exist in sulfolane medium as a Sl-S8-Sl solvate. It has been theoretically established that sulfur can form stable (S8)n clusters in concentrated solutions. An increase in the extent of association (n) of the sulfur cluster leads to a decrease in the extinction coefficient [TD-DFT(TPSSTPSS/6-311+G(d,p))] of the most intense absorption maximum lying at about 50,000 cm-1 while maintaining the shape of the remaining part of the spectrum. The observed pattern qualitatively expresses the spectral regularities of solutions with different concentrations of sulfur in sulfolane. It has been proposed that a model of the absorption spectrum of elemental sulfur suggests a minor contribution of the S12 molecular form (G298°((S12)2) - G298°((S8)3) ≈ -15.5 kJ mol-1). The findings of the study will provide deeper insights into the transformation of molecular forms of sulfur and more precisely analyze processes involving sulfur as an acting species using electronic spectroscopy.

Solvent-induced local environment effect in plasmonic catalysis

Solvents are known to affect the local surface plasmon resonance of metal nanoparticles; however, how solvents can be used to manipulate the interfacial charge and energy transfer in plasmonic catalysis remains to be explored. Here, using NH3 decomposition on a Ru-doped Cu surface as an example, we report density functional theory (DFT) and delta self-consistent field (SCF) calculations, through which we investigate the effect of different protic solvent molecules on interfacial charge transfer by calculating excitation energy of an electronic transition between the metal and the molecular reactant. We find that the H-bonds between water and NH3 can alter the direct interfacial charge transfer due to the shift of the molecular frontier orbitals with respect to the metal Fermi level. These effects are also observed when the H-bonds are formed between methanol (or phenol) and ammonia. We show that the solvent possessing stronger basicity induces a more pronounced effect on the excitation energy. This work thus provides valuable insights for tuning the excitation energy and controlling different routes to channel the photon energy into plasmonic catalysis.

Heavy pnicogen atoms as electron donors in sigma-hole bonds

DFT calculations evaluate the strength of σ-hole bonds formed by ZH3 and ZMe3 (Z = N, P, As, Sb) acting as electron donor. Bond types considered include H-bond, halogen, chalcogen, pnicogen, and tetrel bond to perfluorinated Lewis acids FH, FBr, F2Se F3As, F4Ge, respectively, as well as their monofluorinated analogues. All of the Z atoms can engage in bonds of at least moderate strength, varying from 3 to more than 40 kcal mol-1. In most cases, N forms the strongest bonds, but the falloff from P to Sb is quite mild. However, this pattern is not characteristic of all cases, as for example in the halogen bonds, where the heavier Z atoms are comparable to, or even stronger than N. Most of the bonds are strengthened by replacing the three H atoms of ZH3 by methyl groups, better simulating the situation that would be generally encountered. Structural and NMR shielding data ought to facilitate the identification of these bonds within crystals or in solution.

Competition between Binding to Various Sites of Substituted Imidazoliums

The imidazolium cation has a number of different sites that can interact with a nucleophile. Adding a halogen atom (X) or a chalcogen (YH) group introduces the possibility of an NX···nuc halogen or NY···nuc chalcogen bond, which competes against the various H-bonds (NH and CH donors) as well as the lone pair···π interaction wherein the nucleophile lies above the plane of the cation. Substituted imidazoliums are paired with the NH₃ base, and the various different complexes are evaluated by density functional theory (DFT) calculations. The strength of XB and YB increases quickly along with the size and polarizability of the X/Y atom, and this sort of bond is the strongest for the heavier Br, I, Se, and Te atoms, followed by the NH···N H-bond, but this order reverses for Cl and S. The various CH···N H-bonds are comparable to one another and to the lone pair···π bond, all with interaction energies of 10-13 kcal/mol, values which show very little dependence upon the substituent placed on the imidazolium.

A Paradigm Shift from 2D to 3D: Surface Supramolecular Assemblies and Their Electronic Properties Explored by Scanning Tunneling Microscopy and Spectroscopy

Exploring supramolecular architectures at surfaces plays an increasingly important role in contemporary science, especially for molecular electronics. A paradigm of research interest in this context is shifting from 2D to 3D that is expanding from monolayer, bilayers, to multilayers. Taking advantage of its high-resolution insight into monolayers and a few layers, scanning tunneling microscopy/spectroscopy (STM/STS) turns out a powerful tool for analyzing such thin films on a solid surface. This review summarizes the representative efforts of STM/STS studies of layered supramolecular assemblies and their unique electronic properties, especially at the liquid-solid interface. The superiority of the 3D molecular networks at surfaces is elucidated and an outlook on the challenges that still lie ahead is provided. This review not only highlights the profound progress in 3D supramolecular assemblies but also provides researchers with unusual concepts to design surface supramolecular structures with increasing complexity and desired functionality.

Enhancement of Halogen Bond Strength by Intramolecular H-Bonds

Quantum calculations study the potential of an intramolecular H-bond between the halogen atom (X) of a halobenzene and a substituent placed ortho to it, to amplify the ability of X to engage in a halogen bond (XB) with a Lewis base. H-bonding substituents NH₂, CH₂CH₂OH, CH₂OH, OH, and COOH were added to halobenzenes (X = Cl, Br, I). The amino group had little effect, but those containing OH increased the CX···N XB energy to a NH₃ nucleophile by about 0.5 kcal/mol; the increment associated with COOH is larger, nearly 2 kcal/mol. These energy increments were approximately doubled if two such H-bonding substituents are present. Combining a pair of ortho COOH groups with an electron-withdrawing NO₂ group in the para position has a particularly large effect, raising the XB energy by about 4 kcal/mol, which can amount to as much as a 4-fold magnification.

Fluorescent sensing of non-steroidal anti-inflammatory drugs naproxen and ketoprofen by dansylated squaramide-based receptors

Bis-squaramide receptors L1-L4 bearing a dansyl moiety were synthesised and their potential applications as fluorescent probes towards non steroidal anti-inflammatory drugs naproxen and ketoprofen was investigated. A detailed photophysical characterization in CH₃CN/DMSO solution (9 : 1 v/v) was conducted and demonstrated that the two macrocyclic receptors L1 and L2 show good sensitivity towards ketoprofen with an ON-OFF fluorescent response, while the two open chain receptors L3 and L4 behave similarly with the three guests considered. DFT theoretical calculations carried out on L2 and L4 as model receptors allowed to propose a possible coordination mode towards the guests. Finally, ¹H-NMR spectroscopy in DMSO-d₆/0.5% water solution demonstrated that the four receptors interact with the considered guests via H-bonds.

Does a halogen bond require positive potential on the acid and negative potential on the base?

It is usually expected that formation of a halogen bond (XB) requires that a region of positive electrostatic potential associated with a σ or π-hole on the Lewis acid will interact with the negative potential of the base, either a lone pair or π-bond region. Quantum calculations of model systems suggest this not to be necessary. The placement of electron-withdrawing substituents on the base can reverse the sign of the potential in its lone pair or π-bond region to positive, and this base can nonetheless engage in a XB with the positive σ-hole of a Lewis acid. The reverse scenario is also possible in certain circumstances, as a negatively charged σ-hole can form a XB with the negative lone pair region of a base. Despite these classical Coulombic repulsions, the overall electrostatic interaction is attractive in these XBs, albeit only weakly so. The strengths of these bonds are surprisingly insensitive to changes in the partner molecule. For example, even a wide range in the depth of the σ-hole of the approaching acid yields only a minimal change in the strength of the XB to a base with a positive potential.

Effect of a Substituent on the Properties of Salicylaldehyde Hydrazone Derivatives

The present study investigates the effect of the substitution of salicylaldehyde hydrazones at two selected positions, i.e., the para-position with regard to the proton-donating and proton-accepting centers forming the hydrogen bridge. A detailed analysis of structural data obtained by theoretical approaches and X-ray experiments, together with original resonance Hammett's constants, indicates that the strength of the intramolecular hydrogen bonding present in salicylaldehyde hydrazones can be selectively modulated by substitution of the parent molecular system with the chemical group of known π-electron-donating or -accepting properties. Our findings provide an insight into planning synthesis pathways for salicylaldehyde hydrazone species and predicting their result with regard to their H-bonding and related physical and chemical properties.

Eutectic Mixture Formation and Relaxation Dynamics of Coamorphous Mixtures of Two Benzodiazepine Drugs

The formation of coamorphous mixtures of pharmaceuticals is an interesting strategy to improve the solubility and bioavailability of drugs, while at the same time enhancing the kinetic stability of the resulting binary glass and allowing the simultaneous administration of two active principles. In this contribution, we describe kinetically stable amorphous binary mixtures of two commercial active pharmaceutical ingredients, diazepam and nordazepam, of which the latter, besides being administered as a drug on its own, is also the main active metabolite of the other in the human body. We report the eutectic equilibrium-phase diagram of the binary mixture, which is found to be characterized by an experimental eutectic composition of 0.18 molar fraction of nordazepam, with a eutectic melting point of T_e = 395.4 ± 1.2 K. The two compounds are barely miscible in the crystalline phase. The mechanically obtained mixtures were melted and supercooled to study the glass-transition and molecular-relaxation dynamics of amorphous mixtures at the corresponding concentration. The glass-transition temperature was always higher than room temperature and varied linearly with composition. The T_e was lower than the onset of thermal decomposition of either compound (pure nordazepam decomposes upon melting and pure diazepam well above its melting point), thus implying that the eutectic liquid and glass can be obtained without any degradation of the drugs. The eutectic glass was kinetically stable against crystallization for at least a few months. The relaxation processes of the amorphous mixtures were studied by dielectric spectroscopy, which provided evidence for a single structural (α) relaxation, a single Johari-Goldstein (β) relaxation, and a ring-inversion conformational relaxation of the diazepinic ring, occurring on the same timescale in both drugs. We further characterized both the binary mixtures and pure compounds by FTIR spectroscopy and first-principles density functional theory (DFT) simulations to analyze intermolecular interactions. The DFT calculations confirm the presence of strong attractive forces within the heteromolecular dimer, leading to large dimer interaction energies of the order of -0.1 eV.

Adjusting the balance between hydrogen and chalcogen bonds

A complex is assembled which pairs a carboxyl group of X₁COOH with a 1,2,5-chalcogenadiazole ring containing substituents on its C atoms. The OH of the carboxyl group donates a proton to a N atom of the ring to form a OH⋯N H-bond (HB), while its carbonyl O engages in a Y⋯O chalcogen bond (ChB) with the ring in which Y = S, Se, Te. The ChB is strengthened by enlarging the size of the Y atom from S to Se to Te. Placement of an electron-withdrawing group (EWG) X₁ on the acid strengthens the HB while weakening the ChB; the reverse occurs when EWGs are placed on the ring. By selection of the proper substituents on the two units, it is possible to achieve a near perfect balance between the strengths of these two bonds. These bond strengths are also reflected in the NMR spectroscopic properties of the chemical shielding of the various atoms and the coupling between the nuclei directly involved in each bond.

Testing of Exchange-Correlation Functionals of DFT for a Reliable Description of the Electron Density Distribution in Organic Molecules

Researchers carrying out calculations using the DFT method face the problem of the correct choice of the exchange-correlation functional to describe the quantities they are interested in. This article deals with benchmark calculations aimed at testing various exchange-correlation functionals in terms of a reliable description of the electron density distribution in molecules. For this purpose, 30 functionals representing all rungs of Jacob's Ladder are selected and then the values of some QTAIM-based parameters are compared with their reference equivalents obtained at the CCSD/aug-cc-pVTZ level of theory. The presented results show that the DFT method undoubtedly has the greatest problems with a reliable description of the electron density distribution in multiple strongly polar bonds, such as C=O, and bonds associated with large electron charge delocalization. The performance of the tested functionals turned out to be unsystematic. Nevertheless, in terms of a reliable general description of QTAIM-based parameters, the M11, SVWN, BHHLYP, M06-HF, and, to a slightly lesser extent, also BLYP, B3LYP, and X3LYP functionals turned out to be the worst. It is alarming to find the most popular B3LYP functional in this group. On the other hand, in the case of the electron density at the bond critical point, being the most important QTAIM-based parameter, the M06-HF functional is especially discouraged due to the very poor description of the C=O bond. On the contrary, the VSXC, M06-L, SOGGA11-X, M06-2X, MN12-SX, and, to a slightly lesser extent, also TPSS, TPSSh, and B1B95 perform well in this respect. Particularly noteworthy is the overwhelming performance of double hybrids in terms of reliable values of bond delocalization indices. The results show that there is no clear improvement in the reliability of describing the electron density distribution with climbing Jacob's Ladder, as top-ranked double hybrids are also, in some cases, able to produce poor values compared to CCSD.

Anisotropy, segmental dynamics and polymorphism of crystalline biogenic carboxylic acids

Carboxylic acids of the Krebs cycle possess invaluable biochemical significance. Still, there are severe gaps in the availability of thermodynamic and crystallographic data, as well as ambiguities prevailing in the literature on the thermodynamic characterization and polymorph ranking. Providing an unambiguous description of the structure, thermodynamics and polymorphism of their neat crystalline phases requires a complex multidisciplinary approach. This work presents results of an extensive investigation of the structural anisotropy of the thermal expansion and local dynamics within these crystals, obtained from a beneficial cooperation of NMR crystallography and ab initio calculations of non-covalent interactions. The observed structural anisotropy and spin-lattice relaxation times are traced to large spatial variations in the strength of molecular interactions in the crystal lattice, especially in the orientation of the hydrogen bonds. A completely resolved crystal structure for oxaloacetic acid is reported for the first time. Thanks to multi-instrumental calorimetric effort, this work clarifies phase behavior, determines third-law entropies of the crystals, and states definitive polymorph ranking for succinic and fumaric acids. These thermodynamic observations are then interpreted in terms of first-principles quasi-harmonic calculations of cohesive properties. A sophisticated model capturing electronic, thermal, and configurational-entropic effects on the crystal structure approaches captures the subtle Gibbs energy differences governing polymorph ranking for succinic and fumaric acids, representing another success story of computational chemistry.

Properties and Stabilities of Cyclic and Open Chains of Halogen Bonds

Open and cyclic chains from two to eight units of ICl and IF are constructed and examined by density functional theory (DFT) calculations. These chains contain either I···I or I···X halogen bonds (XBs) where X refers to Cl or F. The closed rings are more stable than the open chains due to the presence of an additional XB and enhanced cooperativity. This pattern is true even for most trimers where there is sizable geometric distortion in the rings. I···F rings are generally more stable than the corresponding I···I cycles as the I···F bond is stronger than I···I even in the simple dimer. However, I···I rings are comparable in energy to I···Cl. It is possible to construct I···I rings of at least as large as eight units, which are held together exclusively by XBs. On the other hand, the maximum possible size of I···X rings is 6. Red shifts are observed in the I-X stretching frequency bands, which magnify as the chain, both cyclic and open, grows longer. The NMR chemical shielding of the I atoms increases for I···I chains but diminishes when I···Cl bonds are present.

Uniting Nonempirical and Empirical Density Functional Approximation Strategies Using Constraint-Based Regularization

In this work, we present a general framework that unites the two primary strategies for constructing density functional approximations (DFAs): nonempirical (NE) constraint satisfaction and empirical (E) data-driven optimization. The proposed method employs B-splines, bell-shaped spline functions with compact support, to construct each inhomogeneity correction factor (ICF). This choice offers several distinct advantages over traditional polynomial expansions by enabling explicit enforcement of linear and nonlinear constraints as well as ICF smoothness using Tikhonov and penalized B-splines (P-splines) regularization. As proof-of-concept, we use the so-called CASE (constrained and smoothed empirical) framework to construct a constraint-satisfying and data-driven global hybrid that exhibits enhanced performance across a diverse set of chemical properties. We argue that the CASE approach can be used to generate DFAs that maintain the physical rigor and transferability of NE-DFAs while leveraging high-quality quantum-mechanical data to remove the arbitrariness of ansatz selection and improve performance.

Proton or Carbene Transfer? On the Dark and Light Reaction of Diazoalkanes with Alcohols

The formal alkylation reaction of OH groups with diazoalkanes under catalyst-free reaction conditions finds broad application in organic synthesis. However, even today, this reaction is mainly limited to the use of diazomethane as reaction partner. In this combined experimental and theoretical study, we aim at a fundamental understanding of the reaction of diazoalkanes with alcohols to make this transformation amenable to a generalized approach towards formal alkylation reactions of alcohols with diazoalkanes. Experimental and theoretical studies suggest a direct proton transfer only in exceptional cases. In a more general setting, such O-H functionalization proceed both under dark and photochemical conditions via a key hydrogen-bonded singlet carbene intermediate that undergoes a protonation-addition mechanism. We conclude with applications of this approach in O-H functionalization reactions of alcohols, including simple fluorinated, halogenated and aliphatic alcohols and showcase functional-group tolerance of this method in the reaction of biologically active and pharmaceutically relevant alcohols.

The radicals of quercetin-derived antioxidants in Triton X-100 micelles

We have employed photoionization with a pulsed laser (5 ns, 355 nm) as a direct access to the radicals of quercetin, five of its monoethers and three of its diethers...

In Search of the Most Stable Molecular Configuration of Heptakis(2,6-O-dimethyl)-β-cyclodextrin and Its Complex with Mianserin: A Comparison of the B3LYP-GD2 and M062X-GD3 Results

Cyclodextrins are well known for their ability to form stable, highly soluble complexes with various substances, which makes them widely used as excipients in food, cosmetics, and pharmaceuticals. In this work, properties of heptakis(2,6-O-dimethyl)-β-cyclodextrin (DM-β-CD) in vacuo and in water, as well as its ability to bind the antidepressant drug mianserin (MIA) in aqueous solution, are investigated computationally. The results are shown to depend strongly on the density functional theory (DFT) applied. The most stable conformers of DM-β-CD found with the B3LYP-GD2 method differ from these indicated by M062X-GD3 and other functionals. According to the latter, two crystal structures, ZULQAY and BOYFOK03, optimized in vacuo and in water, respectively, have the lowest energy. Both the B3LYP-GD2 and M062X-GD3 results show that all tested inclusion and noninclusion complexes of MIA:DM-β-CD in stoichiometry 1:1 are stable in water. However, the structures and their energetic properties obtained with each method differ: in the most stable configurations, different aromatic rings of MIA are embedded inside DM-β-CD, and the corresponding complexation energies (calculated with the 6-31++G(d,p) basis set and corrected for the basis set superposition error) are -29.6 (B3LYP-GD2) and -23.9 (M062X-GD3) kcal/mol. The NMR spectra of DM-β-CD and MIA:DM-β-CD are also compared.

Predicting pKa Values of Quinols and Related Aromatic Compounds with Multiple OH Groups

Quinonoid compounds play central roles as redox-active agents in photosynthesis and respiration and are also promising replacements for inorganic materials currently used in batteries. To design new quinonoid compounds and predict their state of protonation and redox behavior under various conditions, their pK_a values must be known. Methods that can predict the pK_a values of simple phenols cannot reliably handle complex analogues in which multiple OH groups are present and may form intramolecular hydrogen bonds. We have therefore developed a straightforward method based on a linear relationship between experimental pK_a values and calculated differences in energy between quinols and their deprotonated forms. Simple adjustments allow reliable predictions of pK_a values when intramolecular hydrogen bonds are present. Our approach has been validated by showing that predicted and experimental values for over 100 quinols and related compounds differ by an average of only 0.3 units. This accuracy makes it possible to select proper pK_a values when experimental data vary, predict the acidity of quinols and related compounds before they are made, and determine the sites and orders of deprotonation in complex structures with multiple OH groups.

Dependence of the substituent energy on the level of theory

Most often, the substituent effects are described using rather troublesome Hammett constants. Quite recently, it has been proposed to use the so-called substituent energy, which is based on total energies of the X-substituted polycyclic aromatic hydrocarbon and phenyl. This article concerns the influence of the applied level of theory (i.e., both the basis set and the method) on the determined values of the substituent energies. For this purpose, the energies of the NH₂ and NO₂ groups in 16 unique positions of naphthalene, anthracene, tetracene, phenanthrene, and pyrene were calculated using 10 different basis sets and 23 various exchange-correlation functionals representing the entire Jacob's Ladder, from local, through gradient- and meta-gradient-corrected, to hybrid and double-hybrid ones. Additionally, using even larger 6-311++G(2df,2p) basis set, the energies of NH₂ , NO₂ , CN, and Cl were also computed. Both the basis set and the method used have little effect on the substituent energy if the substituent is in the benzene-like position. On the contrary, the effect of the level of theory is pronounced especially in the case of the most spatially crowded 4-substituted phenanthrene. Substituent energies have been shown to be very useful theoretical parameters describing the proximity effect in the substituted derivatives of polycyclic aromatic hydrocarbons.

Relative energetics of CH3CH2O, CH3CHOH, and CH2CH2OH radical products from ethanol dehydrogenation

This study has examined the relative energetics of nine stationary points associated with the three different radical isomers generated by removing a H atom from ethanol at the O atom (ethoxy, CH₃CH₂O), the α C atom (CH₃CHOH), and the β C atom (CH₂CH₂OH). For the first time, CCSD(T) geometry optimizations and harmonic vibrational frequency computations with the cc-pVTZ and aug-cc-pVTZ basis sets have been carried out to characterize two unique minima for each isomer along with three transition state structures with C_s symmetry. Explicitly correlated CCSD(T) computations were also performed to estimate the relative energetics of these nine stationary points near the complete basis set limit. These benchmark results were used to assess the performance of various density functional theory (DFT) and wave function theory methods, and they will help guide method selection for future studies of alcohols and their radicals. The structures generated by abstracting H from the α C atom have significantly lower electronic energies (by at least 7 kcal mol^-1) than the CH₃CH₂O and CH₂CH₂OH radicals. Although previously reported as a minimum on the ground-state surface, the ²A″ C_s structure of the ethoxy radical was found to be a transition state in this study with MP2, CCSD(T), and a number of DFT methods. An implicit solvation model used in conjunction with DFT and MP2 methods did not qualitatively change the relative energies of the isomers, but the results suggest that the local minima for the CH₃CHOH and CH₂CH₂OH radicals could become more energetically competitive in condensed phase environments, such as liquid water and ethanol.

Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals

This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N⁶) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.

BEP-Like Correction of Nonequilibrium Thermodynamic Parameters of the Solvent-Assisted Reactions: The DFT and Ab Initio Study of Hydration, Peroxidation, and Enolization of Acetone and 1,1,1-Trifluoroacetone in Aqueous Solutions

The mechanisms of enolization and reactions of nucleophilic addition to carbonyl compounds were analyzed by density functional theory (DFT) (PBE1PBE) and ab initio (DLPNO-CCSD(T)) level of theory using the interaction of water and hydrogen peroxide with acetone and 1,1,1-trifluoroacetone (TFA) as the reference reactions. The transition states of the studied reactions were localized within the integrated approach that includes both the dielectric continuum theory (polarizable continuum model (PCM)) and the cyclic or two-cluster explicit solvation models. The considered models provide proton transfer in the enolization, hydration, and peroxidation reactions by the Grotthuss mechanism. It is shown that the calculated activation parameters at a sufficiently high level of theory and a sufficiently flexible solvation model can be additionally refined using the Bell-Evans-Polanyi (BEP)-like correction (in a form of the Bell-Evans-Polanyi equation), which is linear scaling of the model potential energy surface according to the equilibrium parameters of the reference reaction (experiment or high-level calculation). Quite good correspondence of the corrected and reference activation parameters and the lower sensitivity of the calculation results to the choice of the solvation model indicate the high reliability of the proposed BEP-like correction technique.

Properties of the Pt(111)/electrolyte electrochemical interface studied with a hybrid DFT-solvation approach

Self-consistent modeling of the interface between solid metal electrode and liquid electrolyte is a crucial challenge in computational electrochemistry. In this contribution, we adopt the effective screening medium reference interaction site method (ESM-RISM) to study the charged interface between a Pt(111) surface that is partially covered with chemisorbed oxygen and an aqueous acidic electrolyte. This method proves to be well suited to describe the chemisorption and charging state of the interface at controlled electrode potential. We present an in-depth assessment of the ESM-RISM parameterization and of the importance of computing near-surface water molecules explicitly at the quantum mechanical level. We found that ESM-RISM is able to reproduce some key interface properties, including the peculiar, non-monotonic charging relation of the Pt(111)/electrolyte interface. The comparison with independent theoretical models and explicit simulations of the interface reveals strengths and limitations of ESM-RISM for modeling electrochemical interfaces.

Structural Water Stabilizes Protein Motifs in Liquid Protein Phase: The Folding Mechanism of Short β-Sheets Coupled to Phase Transition

Macromolecular associates, such as membraneless organelles or lipid-protein assemblies, provide a hydrophobic environment, i.e., a liquid protein phase (LP), where folding preferences can be drastically altered. LP as well as the associated phase change from water (W) is an intriguing phenomenon related to numerous biological processes and also possesses potential in nanotechnological applications. However, the energetic effects of a hydrophobic yet water-containing environment on protein folding are poorly understood. Here, we focus on small β-sheets, the key motifs of proteins, undergoing structural changes in liquid–liquid phase separation (LLPS) and also model the mechanism of energy-coupled unfolding, e.g., in proteases, during W → LP transition. Due to the importance of the accurate description for hydrogen bonding patterns, the employed models were studied by using quantum mechanical calculations. The results demonstrate that unfolding is energetically less favored in LP by ~0.3–0.5 kcal·mol⁻¹ per residue in which the difference further increased by the presence of explicit structural water molecules, where the folded state was preferred by ~1.2–2.3 kcal·mol⁻¹ per residue relative to that in W. Energetics at the LP/W interfaces was also addressed by theoretical isodesmic reactions. While the models predict folded state preference in LP, the unfolding from LP to W renders the process highly favorable since the unfolded end state has >1 kcal·mol⁻¹ per residue excess stabilization.

Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins

Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.

Delocalization-Assisted Transport through Nucleic Acids in Molecular Junctions

The flow of charge through molecules is central to the function of supramolecular machines, and charge transport in nucleic acids is implicated in molecular signaling and DNA repair. We examine the transport of electrons through nucleic acids to understand the interplay of resonant and nonresonant charge carrier transport mechanisms. This study reports STM break junction measurements of peptide nucleic acids (PNAs) with a G-block structure and contrasts the findings with previous results for DNA duplexes. The conductance of G-block PNA duplexes is much higher than that of the corresponding DNA duplexes of the same sequence; however, they do not display the strong even-odd dependence conductance oscillations found in G-block DNA. Theoretical analysis finds that the conductance oscillation magnitude in PNA is suppressed because of the increased level of electronic coupling interaction between G-blocks in PNA and the stronger PNA-electrode interaction compared to that in DNA duplexes. The strong interactions in the G-block PNA duplexes produce molecular conductances as high as 3% G₀, where G₀ is the quantum of conductance, for 5 nm duplexes.

Insights into co-amorphous systems in therapeutic drug delivery

Co-amorphous (CAM) systems are promising drug-delivery systems in the arena of therapeutic drug delivery, addressing the poor aqueous solubility of drugs by enhancing solubility and thereby improving the oral bioavailability and therapeutic effect of the drug. A CAM system is a single-phase homogeneous blend of two or more low molecular weight molecules that can be drug–drug or drug–co-former, stabilized via intermolecular interactions, adding the benefit of thermodynamic stability. This review covers the fundamentals of CAM systems and recent advances in formulation development. In particular, we strive to address the theoretical, molecular, technical and biopharmaceutical aspects, advantages over polymeric amorphous solid dispersions, mechanisms of stabilization of amorphous forms, insights into unexplored in silico tools in excipient selection and regulatory viewpoints.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Tuning of ORR activity through the stabilization of the adsorbates by hydrogen bonding with substituent groups

Metallocorroles and metalloporphyrins (M-N-C) are some of the best alternative molecular catalysts for the replacement of the expensive platinum-group metals (PGM) in oxygen reduction reaction (ORR) catalysis in polymer electrolyte membrane (PEM) fuel cells. To date, Co-based corroles have shown the best performance, but still suffer from the poor stability and the toxicity of the Co metal. Mn-N-C are more stable than Co-N-C, and are also less reactive towards peroxide formation. In this work, using first-principles density functional theory calculations, we study the improvement of the Mn-based corrole ORR activity by exploiting hydrogen bonding with substituent groups to modify the adsorption energies of the ORR intermediates and obtain higher onset potential (Vonset) values. We found that by using phenyl acetic acid as a substituent, Vonset increased from 0.54 V for the unsubstituted corrole to ∼0.9 V which is competitive with the Vonset of the Co-based corroles. Our results suggest that hydrogen bonding with substituent groups should be considered in the analysis and design of the reactivity of active sites in non-PGM ORR catalysts.

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo-7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm ) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

Intramolecular Hydrogen Bonds in Tip-Functionalized Single-Walled Carbon Nanotubes as pH-Sensitive Gates

Since their discovery, carbon nanotubes and other related nanomaterials are in the spotlight due to their unique molecular structures and properties, having a wide range of applications. The cage-like structure of carbon nanotubes is especially appealing as a route to confine molecules, isolating them from the solvent medium. This study aims to explore and characterize, through density functional theory (DFT) calculations, covalent tip-functionalization of single-walled carbon nanotubes (SWCNTS) with carboxymethyl moieties that establish pH sensitive molecular gates. The response of the molecular gate to pH fluctuations arises from variations in the noncovalent interactions between functionalized groups, which depend on the extent of protonation, leading to conformational changes. Overall, the hydrogen bonds present in the molecular models under study, as evaluated through topological analysis and pK_a calculations, suggest that functionalized SWCNTs may be suitable for the design of drug delivery systems to enhance the efficiency of some pharmacological treatments, or even in the area of catalysis and separation processes, through their incorporation in nanocomposites.

Density Functional Theory in the Prediction of Mutagenicity: A Perspective

As a field, computational toxicology is concerned with using in silico models to predict and understand the origins of toxicity. It is fast, relatively inexpensive, and avoids the ethical conundrum of using animals in scientific experimentation. In this perspective, we discuss the importance of computational models in toxicology, with a specific focus on the different model types that can be used in predictive toxicological approaches toward mutagenicity (SARs and QSARs). We then focus on how quantum chemical methods, such as density functional theory (DFT), have previously been used in the prediction of mutagenicity. It is then discussed how DFT allows for the development of new chemical descriptors that focus on capturing the steric and energetic effects that influence toxicological reactions. We hope to demonstrate the role that DFT plays in understanding the fundamental, intrinsic chemistry of toxicological reactions in predictive toxicology.

NMR crystallography of molecular organics

Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.

Microhydration of a Carbonyl Group: How does the Molecular Electrostatic Potential (MESP) Impact the Formation of (H2O)n:(R2C═O)Complexes?

The presence of a carbonyl group in a molecule usually leads to the identification of a π-hole on the molecular electrostatic potential (MESP) of the species. How does this electrophilic site influence the formation of microhydrated complexes? To address this point, a panel of R₂CO solutes with various MESPs was selected, and we identified the structures and properties of several complexes containing one, two, three and six water molecules. The following solutes were considered in the present study: H₂CO, F₂CO, Cl₂CO,(NC)₂CO and H₂C═CO. Geometry optimizations and frequency calculations were carried out at the LC-ωPBE/6-311++G(d,p) level, with the GD3BJ empirical correction for dispersion. For a number of n water molecules around the R₂CO solute, the structure and the features of the most stable (H₂O)_n:(R₂CO) complexes are highly dependent on the MESP of the isolated R₂CO solute. The formation of pi-hole bondings appears to play a decisive role in the initiation of a three-dimensional organization of water molecules around the solute.