Different Ways of Hydrogen Bonding in Water - Why Does Warm Water Freeze Faster than Cold Water?

Systematic Study of Different Types of Interactions in α-, β- and γ-Cyclodextrin: Quantum Chemical Investigation

In this work, comprehensive ab initio quantum chemical calculations using the DFT level of theory were performed to characterize the stabilization interactions (H-bonding and hyperconjugation effects) of two stable symmetrical conformations of α-, β-, and γ-cyclodextrins (CDs). For this purpose, we analyzed the electron density using "Atom in molecules" (AIM), "Natural Bond Orbital" (NBO), and energy decomposition method (CECA) in 3D and in Hilbert space. We also calculated the H-bond lengths and OH vibrational frequencies. In every investigated CD, the quantum chemical descriptors characterizing the strength of the interactions between the H-bonds of the primary OH (or hydroxymethyl) and secondary OH groups are examined by comparing the same quantity calculated for ethylene glycol, α-d-glucose (α-d-Glcp) and a water cluster as reference systems. By using these external standards, we can characterize more quantitatively the properties of these bonds (e.g., strength). We have demonstrated that bond critical points (BCP) of intra-unit H-bonds are absent in cyclodextrins, similar to α-d-Glcp and ethylene glycol. In contrast, the CECA analysis showed the existence of an exchange (bond-like) interaction between the interacting O…H atoms. Consequently, the exchange interaction refers to a chemical bond, namely the H-bond between two atoms, unlike BCP, which is not suitable for its detection.

Mechanistic Insights into S-Depalmitolyse Activity of Cln5 Protein Linked to Neurodegeneration and Batten Disease: A QM/MM Study

Ceroid lipofuscinosis neuronal protein 5 (Cln5) is encoded by the CLN5 gene. The genetic variants of this gene are associated with the CLN5 form of Batten disease. Recently, the first crystal structure of Cln5 was reported. Cln5 shows cysteine palmitoyl thioesterase S-depalmitoylation activity, which was explored via fluorescent emission spectroscopy utilizing the fluorescent probe DDP-5. In this work, the mechanism of the reaction between Cln5 and DDP-5 was studied computationally by applying a QM/MM methodology at the ωB97X-D/6-31G(d,p):AMBER level. The results of our study clearly demonstrate the critical role of the catalytic triad Cys280-His166-Glu183 in S-depalmitoylation activity. This is evidenced through a comparison of the pathways catalyzed by the Cys280-His166-Glu183 triad and those with only Cys280 involved. The computed reaction barriers are in agreement with the catalytic efficiency. The calculated Gibb's free-energy profile suggests that S-depalmitoylation is a rate-limiting step compared to the preceding S-palmitoylation, with barriers of 26.1 and 25.3 kcal/mol, respectively. The energetics were complemented by monitoring the fluctuations in the electron density distribution through NBO charges and bond strength alterations via local mode stretching force constants during the catalytic pathways. This comprehensive protocol led to a more holistic picture of the reaction mechanism at the atomic level. It forms the foundation for future studies on the effects of gene mutations on both the S-palmitoylation and S-depalmitoylation steps, providing valuable data for the further development of enzyme replacement therapy, which is currently the only FDA-approved therapy for childhood neurodegenerative diseases, including Batten disease.

Mpemba effect for a Brownian particle trapped in a single well potential

The Mpemba effect refers to the counterintuitive phenomenon of a hotter system equilibrating faster than a colder system when both are quenched to the same low temperature. For a Brownian particle trapped in a piecewise linear single well potential that is devoid of any other metastable minima, we show the existence of the Mpemba effect for a wide range of parameters through an exact solution. This result challenges the prevalent explanation of the Mpemba effect that requires the energy landscape to be rugged with multiple minima. We also demonstrate the existence of inverse and strong Mpemba effects.

Mpemba effect in a Langevin system: Population statistics, metastability, and other exact results

The Mpemba effect is a fingerprint of the anomalous relaxation phenomenon wherein an initially hotter system equilibrates faster than an initially colder system when both are quenched to the same low temperature. Experiments on a single colloidal particle trapped in a carefully shaped double well potential have demonstrated this effect recently [A. Kumar and J. Bechhoefer, Nature 584, 64 (2020)]. In a similar vein, here, we consider a piece-wise linear double well potential that allows us to demonstrate the Mpemba effect using an exact analysis based on the spectral decomposition of the corresponding Fokker-Planck equation. We elucidate the role of the metastable states in the energy landscape as well as the initial population statistics of the particles in showcasing the Mpemba effect. Crucially, our findings indicate that neither the metastability nor the asymmetry in the potential is a necessary or a sufficient condition for the Mpemba effect to be observed.

Relaxation Shortcuts through Boundary Coupling

When a hot system cools down faster than an equivalent cold one, it exhibits the Mpemba effect (ME). This counterintuitive phenomenon was observed in several systems including water, magnetic alloys, and polymers. In most experiments the system is coupled to the bath through its boundaries, but all theories so far assumed bulk coupling. Here we build a general framework to characterize anomalous relaxations through boundary coupling, and present two emblematic setups: a diffusing particle and an Ising antiferromagnet. In the latter, we show that the ME can survive even arbitrarily weak couplings.

Hyperacceleration of quantum thermalization dynamics by bypassing long-lived coherences: An analytical treatment

We develop a perturbative technique for solving Markovian quantum dissipative dynamics, with the perturbation parameter being a small gap in the eigenspectrum. As an example, we apply the technique and straightforwardly obtain analytically the dynamics of a three-level system with quasidegenerate excited states, where quantum coherences persist for very long times, proportional to the inverse of the energy splitting squared. We then show how to bypass this long-lived coherent dynamics and accelerate the relaxation to thermal equilibration in a hyper-exponential manner, a Markovian quantum-assisted Mpemba-like effect. This hyperacceleration of the equilibration process manifests if the initial state is carefully prepared, such that its coherences precisely store the amount of population relaxing from the initial condition to the equilibrium state. Our analytical method for solving quantum dissipative dynamics readily provides equilibration timescales, and as such it reveals how coherent and incoherent effects interlace in the dynamics. It further advises on how to accelerate relaxation processes, which is desirable when long-lived quantum coherences stagnate dynamics.

Radial distribution and hydrogen bonded network graphs of alcohol-aniline binary mixture

CONTEXT

Hydrogen bonds play a vital role in the stability and functioning of biomolecules. Suitable binary liquids are often used as prototypes for the study of biologically significant hydrogen bond studies and their intricate networks. Often, such systems show deviations in their physico-chemical properties from ideal conditions. As a continuation of our research interest in biologically important hydrogen-bonded systems, this paper reports the classical molecular dynamic studies on mixtures of aniline with 8 primary alcohols (C_RH_2R+1-OH, R = 1 to 8) for the complete concentration range. The energetics results indicate the predominance of OH--O interactions over other hydrogen bonds. Structures in the network are analyzed using radial distribution function (RDF), hydrogen bond statistics, and graph theoretical analysis (GTA). Coordination numbers, hydrogen bond statistics, and GTA show a bunching of alcohol-alcohol hydrogen bonds for lower aniline concentrations, while the aniline-aniline interactions are not affected by changes in the concentration.

METHODS

Interaction energies are calculated using B3LYP/6-311G++(d, p) density functional theory using Gaussian-09. The molecular dynamics simulations are carried out using GROMACS (V 2020.6) with the OPLS/AA force field and the simulation box is visualized using VMD. The NetworkX Python package is used for GTA calculation.

Dihydrogen Bonding-Seen through the Eyes of Vibrational Spectroscopy

In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms-in-Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode-mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H3+(H2)n systems, we conclude that the bond strength in these crystal-like structures makes them interesting for potential hydrogen storage applications.

Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study

Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.

Mechanistic Details of the Sharpless Epoxidation of Allylic Alcohols—A Combined URVA and Local Mode Study

In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationary points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.

Bonding in Nitrile Photo-dissociating Ruthenium Drug Candidates --A Local Vibrational Mode Study

In this work, we investigated bonding features 15 ruthenium complexes of the type [Ru(tpy)(L)-(CH3CN)]n+, containing the tridentate tpy ligand (tpy = 2,2':6',2'--terpyridine) and various bidentate ancillary ligands, 12 compounds originally synthesized by Loftus et al. (J. Phys. Chem. C 123, 10291-10299 (2019)) complemented with three additional complexes. The main focus of our work was to relate these local features to the experimental data of Loftus et al. which assess the efficiency of nitrile release in an indirect way via observed quantum yields for ruthenium water association after nitrile release. As a tool to quantitatively assess Ru-NC and Ru-L bonding we utilized the local vibrational mode analysis complemented by the topological analysis of the electron density and the natural bond orbital analysis. Interestingly, the stronger Ru-NC bonds have the greater observed quantum yields, leading to the conclusion that the observed quantum yields are a result of a complex interplay of several processes excluding a direct relationship between QY and Ru-NC or Ru-L bond strengths. We identified the ST splitting as one of the key players and not the Ru-NC bond strength, as one may have thought. In summary, this work has presented a modern computational tool set for the investigation of bonding features applied to nitrile photo-dissociating ruthenium drug candidates forming a valuable basis for future design and fine tuning of nitrile releasing ruthenium compounds, as well as for the understanding of how local properties affect overall experimental outcomes.

Mpemba effect - the effect of time

By analyzing the relation between time and speed, the relation between time and gravitational field, the gravitational redshift of photon and the black-body radiation theorem, the conclusion that time on an object is proportional to the fourth power of the absolute temperature of the object is obtained. Applying the above conclusion about the nature of time, the author analyzes the Mpemba effect and the inverse Mpemba effect, and reaches the following conclusion: the Mpemba effect is the time effect produced when heat flows from objects into space, and the “inverse” Mpemba effect is the time effect produced when heat flows from space into objects.

Uranium: The Nuclear Fuel Cycle and Beyond

This review summarizes the recent developments regarding the use of uranium as nuclear fuel, including recycling and health aspects, elucidated from a chemical point of view, i.e., emphasizing the rich uranium coordination chemistry, which has also raised interest in using uranium compounds in synthesis and catalysis. A number of novel uranium coordination features are addressed, such the emerging number of U(II) complexes and uranium nitride complexes as a promising class of materials for more efficient and safer nuclear fuels. The current discussion about uranium triple bonds is addressed by quantum chemical investigations using local vibrational mode force constants as quantitative bond strength descriptors based on vibrational spectroscopy. The local mode analysis of selected uranium nitrides, N≡U≡N, U≡N, N≡U=NH and N≡U=O, could confirm and quantify, for the first time, that these molecules exhibit a UN triple bond as hypothesized in the literature. We hope that this review will inspire the community interested in uranium chemistry and will serve as an incubator for fruitful collaborations between theory and experimentation in exploring the wealth of uranium chemistry.

BF3–Catalyzed Diels–Alder Reaction between Butadiene and Methyl Acrylate in Aqueous Solution—An URVA and Local Vibrational Mode Study

In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by the CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.

Anomalous heating in a colloidal system

A cold system normally takes longer to warm up than a cool system. Yet recent theoretical studies have suggested that the reverse may sometimes be possible. Here, using a colloidal particle in a heat bath, we present experimental evidence for this inverse Mpemba effect. By carefully choosing the energy landscape, we can make the cold system heat up exponentially faster than the heating rate of cool systems. While similar behavior has been seen in systems that cool down—the Mpemba effect—we find that entropic effects generally make anomalous heating harder to observe than anomalous cooling.

We report anomalous heating in a colloidal system, an experimental observation of the inverse Mpemba effect, where for two initial temperatures lower than the temperature of the thermal bath, the colder of the two systems heats up faster when coupled to the same thermal bath. For an overdamped, Brownian colloidal particle moving in a tilted double-well potential, we find a nonmonotonic dependence of the heating times on the initial temperature of the system. Entropic effects make the inverse Mpemba effect generically weaker—harder to observe—than the usual Mpemba effect (anomalous cooling). We also observe a strong version of anomalous heating, where a cold system heats up exponentially faster than systems prepared under slightly different conditions.

Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture

Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.

Assessing the Intrinsic Strengths of Ion–Solvent and Solvent–Solvent Interactions for Hydrated Mg2+ Clusters

Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Halogen Bonding Involving I2 and d8 Transition-Metal Pincer Complexes

We systematically investigated iodine–metal and iodine–iodine bonding in van Koten’s pincer complex and 19 modifications changing substituents and/or the transition metal with a PBE0–D3(BJ)/aug–cc–pVTZ/PP(M,I) model chemistry. As a novel tool for the quantitative assessment of the iodine–metal and iodine–iodine bond strength in these complexes we used the local mode analysis, originally introduced by Konkoli and Cremer, complemented with NBO and Bader’s QTAIM analyses. Our study reveals the major electronic effects in the catalytic activity of the M–I–I non-classical three-center bond of the pincer complex, which is involved in the oxidative addition of molecular iodine I2 to the metal center. According to our investigations the charge transfer from the metal to the σ* antibonding orbital of the I–I bond changes the 3c–4e character of the M–I–I three-center bond, which leads to weakening of the iodine I–I bond and strengthening of the metal–iodine M–I bond, facilitating in this way the oxidative addition of I2 to the metal. The charge transfer can be systematically modified by substitution at different places of the pincer complex and by different transition metals, changing the strength of both the M–I and the I2 bonds. We also modeled for the original pincer complex how solvents with different polarity influence the 3c–4e character of the M–I–I bond. Our results provide new guidelines for the design of pincer complexes with specific iodine–metal bond strengths and introduce the local vibrational mode analysis as an efficient tool to assess the bond strength in complexes.

Complexes of Formaldehyde and α-Dicarbonyls with Hydroxylamine: FTIR Matrix Isolation and Theoretical Study

The interactions of formaldehyde (FA), glyoxal (Gly) and methylglyoxal (MGly) with hydroxylamine (HA) isolated in solid argon and nitrogen were studied using FTIR spectroscopy and ab initio methods. The spectra analysis indicates the formation of two types of hydrogen-bonded complexes between carbonyl and hydroxylamine in the studied matrices. The cyclic planar complexes are stabilized by O–H⋯O(C), and C–H⋯N interactions and the nonplanar complexes are stabilized by O–H⋯O(C) bond. Formaldehyde was found to form with hydroxylamine, the cyclic planar complex and methylglyoxal, the nonplanar one in both argon and nitrogen matrices. In turn, glyoxal forms with hydroxylamine the most stable nonplanar complex in solid argon, whereas in solid nitrogen, both types of the complex are formed.

Comment on "Exploring nature and predicting strength of hydrogen bonds: A correlation analysis between atoms-in-molecules descriptors, binding energies, and energy components of symmetry-adapted perturbation theory"

We evaluate the correlation between binding energy (BE) and electron density ρ(r) at the bond critical point for 28 neutral hydrogen bonds, recently reported by Emamian and co-workers (J. Comput. Chem., 2019, 40, 2868). As an efficient tool, we use local stretching force constant k HB a derived from the local vibrational mode theory of Konkoli and Cremer. We compare the physical nature of BE versus k HB a , and provide an important explanation for cases with significant deviation in the BE- k HB a relation as well as in the BE-ρ(r) correlation. We also show that care has to be taken when different hydrogen bond strength measures are compared. The BE is a cumulative hydrogen bond strength measure while k HB a is a local measure of hydrogen bond strength covering different aspects of bonding. A simplified and unified description of hydrogen bonding is not always possible and needs an in-depth understanding of the systems involved.

Local vibrational mode analysis of ion-solvent and solvent-solvent interactions for hydrated Ca2+ clusters

Hydrated calcium ion clusters have received considerable attention due to their essential role in biological processes such as bone development, hormone regulation, blood coagulation, and neuronal signaling. To better understand the biological role of the cation, the interactions between the Ca²⁺ ions and water molecules have been frequently investigated. However, a quantitative measure for the intrinsic Ca-O (ion-solvent) and intermolecular hydrogen bond (solvent-solvent) interactions has been missing so far. Here, we report a topological electron density analysis and a natural population analysis to analyze the nature of these interactions for a set of 14 hydrated calcium clusters via local mode stretching force constants obtained at the ωB97X-D/6-311++G(d,p) level of theory. The results revealed that the strength of inner Ca-O interactions for Ca(H₂O)_n ²⁺ (n = 1-8) clusters correlates with the electron density. The application of a second hydration shell to Ca(H₂O)_n ²⁺ (n = 6-8) clusters resulted in stronger Ca-O interactions where a larger electron charge transfer between lp(O) of the first hydration shell and the lower valence of Ca prevailed. The strength of the intermolecular hydrogen bonds, formed between the first and second hydration shells, became stronger when the charge transfers between hydrogen bond (HB) donors and HB acceptors were enhanced. From the local mode stretching force constants of implicitly and explicitly solvated Ca²⁺, we found the six-coordinated cluster to possess the strongest stabilizations, and these results prove that the intrinsic bond strength measures for Ca-O and hydrogen bond interactions form new effective tools to predict the coordination number for the hydrated calcium ion clusters.

Modeling Hydrogen Release from Water with Borane and Alane Catalysts: A Unified Reaction Valley Approach

The unified reaction valley approach combined with the local vibrational mode and ring puckering analysis is applied to investigate the hydrogen evolution from water in the presence of small hydrides such as BH₃, metal hydrides as AlH₃, and their derivatives. We studied a series of reactions involving BH₃, AlH₃, B₂H₆, Al₂H₆, and AlH₃BH₃ with one- and two-water molecules, considering multiple reaction paths. In addition, the influence of the aqueous medium was examined. A general reaction mechanism was identified for most of the reactions. Those that deviate could be associated with unusually high reaction barriers with no hydrogen release. The charge transfer along the reaction path suggests that a viable hydrogen release is achieved when the catalyst adopts the role of a charge donor during the chemical processes. The puckering analysis showed that twistboat and boat forms are the predominant configurations in the case of an intermediate six-membered ring formation, which influences the activation barrier. The local mode analysis was used as a tool to detect the H-H bond formation as well as to probe catalyst regenerability. Based on the correlation between the activation energy and the change in the charge separation for cleaving O-H and B(Al)-H bonds, two promising subsets of reactions could be identified along with prescriptions for lowering the reaction barrier individually with electron-donating/withdrawing substituents.

Computational analysis of vibrational frequencies and rovibrational spectroscopic constants of hydrogen sulfide dimer using MP2 and CCSD(T)

Previous studies have shown that the weakly bonded H₂S dimer demands high level quantum chemical calculations to reproduce experimental values. We investigated the hydrogen bonding of H₂S dimer using MP2 and CCSD(T) levels of theory in combination with aug-cc-pV(D,T,Q)Z basis sets. More precisely, the binding energies, potential energy curves, rovibrational spectroscopic constants, decomposition lifetime, and normal vibrational frequencies were calculated. In addition, we introduced the local mode analysis of Konkoli-Cremer to quantify the hydrogen bonding in the H₂S dimer as well as providing for the first time the comprehensive decomposition of normal vibrational modes into local modes contributions, and a decomposition lifetime based on rate constant. The local mode force constant of the H₂S dimer hydrogen bond is smaller than that of the water dimer, in accordance with the weaker hydrogen bonding in the H₂S dimer.

PyVibMS: a PyMOL plugin for visualizing vibrations in molecules and solids

Visualizing vibrational motions calculated with different ab initio packages requires dedicated post-processing tools. Here, we present a PyMOL plugin called PyVibMS for visualizing the vibrational motions for both molecular and solid systems calculated by mainstream quantum chemical computer programs including Gaussian, Q-Chem, VASP, and CRYSTAL. Benefiting from the continuing development of the PyMOL platform, PyVibMS provides powerful functionalities and user-friendly interface. PyVibMS was written in Python and its open-source nature makes it flexible and sustainable. As an example, the motions of the Konkoli-Cremer local vibrational modes are shown in this work for the first time. PyVibMS is freely available at https://github.com/smutao/PyVibMS . Graphical abstract In this work, a PyMOL plugin named PyVibMS is developed to visualize molecular and lattice vibrations.

Observing the Mpemba effect with minimal bias and the value of the Mpemba effect to scientific outreach and engagement

The Mpemba effect is the assertion that it is quicker to cool or freeze water when the initial temperature is high. We define the Mpemba effect to have been observed when two samples of water, one initially warmer than the other, are cooled and the initially hotter sample freezes or cools to a prescribed temperature in less time; this assumes that the samples are of the same mass, are at least as pure as drinking water, the cooling is delivered to both samples in the same manner and the same level of insulation is applied to both. Under this definition, we enable repeatable observations of the Mpemba effect by systematically introducing increased nucleation sites (in our case by roughening the container walls with sandpaper) within the cooling environment of the initially warmer sample. We remain able to observe the Mpemba effect when the enthalpy of cooling and freezing the initially warmer sample exceeds that of the initially cooler sample by over 50%, corresponding to a difference in the initial temperatures of around 50 °C. The context of this study, initiated by a high-school student who then carried out all of the experiments, highlights the value of the Mpemba effect as a tool for scientific learning and engagement.

Mpemba effect in driven granular Maxwell gases

A Mpemba effect refers to the counterintuitive result that, when quenched to a low temperature, a system at higher temperature may equilibrate faster than one at intermediate temperatures. This effect has recently been demonstrated in driven granular gases, both for smooth as well as rough hard-sphere systems based on a perturbative analysis. In this paper, we consider the inelastic driven Maxwell gas, a simplified model for a granular gas, where the rate of collision is assumed to be independent of the relative velocity. Through an exact analysis, we determine the conditions under which the Mpemba effect is present in this model. For monodispersed gases, we show that the Mpemba effect is present only when the initial states are allowed to be nonstationary, while for bidispersed gases, it is present for some steady-state initial states. We also demonstrate the existence of the strong Mpemba effect for bidispersed Maxwell gas, wherein the system at higher temperature relaxes to a final steady state at an exponentially faster rate leading to smaller equilibration time.

Local Vibrational Mode Analysis of π–Hole Interactions between Aryl Donors and Small Molecule Acceptors

11 aryl–lone pair and three aryl–anion π –hole interactions are investigated, along with the argon–benzene dimer and water dimer as reference compounds, utilizing the local vibrational mode theory, originally introduced by Konkoli and Cremer, to quantify the strength of the π –hole interaction in terms of a new local vibrational mode stretching force constant between the two engaged monomers, which can be conveniently used to compare different π –hole systems. Several factors have emerged which influence strength of the π –hole interactions, including aryl substituent effects, the chemical nature of atoms composing the aryl rings/ π –hole acceptors, and secondary bonding interactions between donors/acceptors. Substituent effects indirectly affect the π –hole interaction strength, where electronegative aryl-substituents moderately increase π –hole interaction strength. N-aryl members significantly increase π –hole interaction strength, and anion acceptors bind more strongly with the π –hole compared to charge neutral acceptors (lone–pair donors). Secondary bonding interactions between the acceptor and the atoms in the aryl ring can increase π –hole interaction strength, while hydrogen bonding between the π –hole acceptor/donor can significantly increase or decrease strength of the π –hole interaction depending on the directionality of hydrogen bond donation. Work is in progress expanding this research on aryl π –hole interactions to a large number of systems, including halides, CO, and OCH3− as acceptors, in order to derive a general design protocol for new members of this interesting class of compounds.

Exploring the Mechanism of Catalysis with the Unified Reaction Valley Approach (URVA)—A Review

The unified reaction valley approach (URVA) differs from mainstream mechanistic studies, as it describes a chemical reaction via the reaction path and the surrounding reaction valley on the potential energy surface from the van der Waals region to the transition state and far out into the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting molecules is registered by a change in their normal vibrational modes and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction with curvature minima reflecting minimal change and curvature maxima, the location of important chemical events such as bond breaking/forming, charge polarization and transfer, rehybridization, etc. A unique decomposition of the path curvature into internal coordinate components provides comprehensive insights into the origins of the chemical changes taking place. After presenting the theoretical background of URVA, we discuss its application to four diverse catalytic processes: (i) the Rh catalyzed methanol carbonylation—the Monsanto process; (ii) the Sharpless epoxidation of allylic alcohols—transition to heterogenous catalysis; (iii) Au(I) assisted [3,3]-sigmatropic rearrangement of allyl acetate; and (iv) the Bacillus subtilis chorismate mutase catalyzed Claisen rearrangement—and show how URVA leads to a new protocol for fine-tuning of existing catalysts and the design of new efficient and eco-friendly catalysts. At the end of this article the pURVA software is introduced. The overall goal of this article is to introduce to the chemical community a new protocol for fine-tuning existing catalytic reactions while aiding in the design of modern and environmentally friendly catalysts.

In Situ Assessment of Intrinsic Strength of X-I⋯OA-Type Halogen Bonds in Molecular Crystals with Periodic Local Vibrational Mode Theory

Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1–0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger’s rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the σ-hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the n→σ∗(X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials.

Metal-Halogen Bonding Seen through the Eyes of Vibrational Spectroscopy

Incorporation of a metal center into halogen-bonded materials can efficiently fine-tune the strength of the halogen bonds and introduce new electronic functionalities. The metal atom can adopt two possible roles: serving as halogen acceptor or polarizing the halogen donor and acceptor groups. We investigated both scenarios for 23 metal–halogen dimers trans-M(Y₂)(NC₅H₄X-3)₂ with M = Pd(II), Pt(II); Y = F, Cl, Br; X = Cl, Br, I; and NC₅H₄X-3 = 3-halopyridine. As a new tool for the quantitative assessment of metal–halogen bonding, we introduced our local vibrational mode analysis, complemented by energy and electron density analyses and electrostatic potential studies at the density functional theory (DFT) and coupled-cluster single, double, and perturbative triple excitations (CCSD(T)) levels of theory. We could for the first time quantify the various attractive contacts and their contribution to the dimer stability and clarify the special role of halogen bonding in these systems. The largest contribution to the stability of the dimers is either due to halogen bonding or nonspecific interactions. Hydrogen bonding plays only a secondary role. The metal can only act as halogen acceptor when the monomer adopts a (quasi-)planar geometry. The best strategy to accomplish this is to substitute the halo-pyridine ring with a halo-diazole ring, which considerably strengthens halogen bonding. Our findings based on the local mode analysis provide a solid platform for fine-tuning of existing and for design of new metal–halogen-bonded materials.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH₂OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.

Multicentered hydrogen bonding in 1-[(1-de-oxy-β-d-fructo-pyranos-1-yl)aza-nium-yl]cyclo-pentane-carboxyl-ate ('d-fructose-cyclo-leucine')

The title compound, C12H21NO7, (I), is conformationally unstable; the predominant form present in its solution is the β-pyran-ose form (74.3%), followed by the β- and α-furan-oses (12.1 and 10.2%, respectively), α-pyran-ose (3.4%), and traces of the acyclic carbohydrate tautomer. In the crystalline state, the carbohydrate part of (I) adopts the 2 C 5 β-pyran-ose conformation, and the amino acid portion exists as a zwitterion, with the side chain cyclo-pentane ring assuming the E 9 envelope conformation. All heteroatoms are involved in hydrogen bonding that forms a system of anti-parallel infinite chains of fused R 3 3(6) and R 3 3(8) rings. The mol-ecule features extensive intra-molecular hydrogen bonding, which is uniquely multicentered and involves the carboxyl-ate, ammonium and carbohydrate hy-droxy groups. In contrast, the contribution of inter-molecular O⋯H/H⋯O contacts to the Hirshfeld surface is relatively low (38.4%), as compared to structures of other d-fructose-amino acids. The 1H NMR data suggest a slow rotation around the C1-C2 bond in (I), indicating that the intra-molecular heteroatom contacts survive in aqueous solution of the mol-ecule as well.

A Continuum from Halogen Bonds to Covalent Bonds: Where Do λ3 Iodanes Fit?

The intrinsic bonding nature of λ 3 -iodanes was investigated to determine where its hypervalent bonds fit along the spectrum between halogen bonding and covalent bonding. Density functional theory with an augmented Dunning valence triple zeta basis set ( ω B97X-D/aug-cc-pVTZ) coupled with vibrational spectroscopy was utilized to study a diverse set of 34 hypervalent iodine compounds. This level of theory was rationalized by comparing computational and experimental data for a small set of closely-related and well-studied iodine molecules and by a comparison with CCSD(T)/aug-cc-pVTZ results for a subset of the investigated iodine compounds. Axial bonds in λ 3 -iodanes fit between the three-center four-electron bond, as observed for the trihalide species IF 2 − and the covalent FI molecule. The equatorial bonds in λ 3 -iodanes are of a covalent nature. We explored how the equatorial ligand and axial substituents affect the chemical properties of λ 3 -iodanes by analyzing natural bond orbital charges, local vibrational modes, the covalent/electrostatic character, and the three-center four-electron bonding character. In summary, our results show for the first time that there is a smooth transition between halogen bonding → 3c–4e bonding in trihalides → 3c–4e bonding in hypervalent iodine compounds → covalent bonding, opening a manifold of new avenues for the design of hypervalent iodine compounds with specific properties.

Quantitative Assessment of Tetrel Bonding Utilizing Vibrational Spectroscopy

A set of 35 representative neutral and charged tetrel complexes was investigated with the objective of finding the factors that influence the strength of tetrel bonding involving single bonded C, Si, and Ge donors and double bonded C or Si donors. For the first time, we introduced an intrinsic bond strength measure for tetrel bonding, derived from calculated vibrational spectroscopy data obtained at the CCSD(T)/aug-cc-pVTZ level of theory and used this measure to rationalize and order the tetrel bonds. Our study revealed that the strength of tetrel bonds is affected by several factors, such as the magnitude of the σ-hole in the tetrel atom, the negative electrostatic potential at the lone pair of the tetrel-acceptor, the positive charge at the peripheral hydrogen of the tetrel-donor, the exchange-repulsion between the lone pair orbitals of the peripheral atoms of the tetrel-donor and the heteroatom of the tetrel-acceptor, and the stabilization brought about by electron delocalization. Thus, focusing on just one or two of these factors, in particular, the σ-hole description can only lead to an incomplete picture. Tetrel bonding covers a range of -1.4 to -26 kcal/mol, which can be strengthened by substituting the peripheral ligands with electron-withdrawing substituents and by positively charged tetrel-donors or negatively charged tetrel-acceptors.