Folding of the Hydrogen Bond Network of H+(CH3OH)7 with Rare Gas Tagging

Hydrogen bond network structures of protonated 2,2,2-trifluoroethanol/ethanol mixed clusters probed by infrared spectroscopy combined with a deep-learning structure sampling approach: the origin of the linear type network preference in protonated fluoroalcohol clusters

While preferential hydrogen bond network structures of cold protonated alcohol clusters H+(ROH)n are generally switched from a linear type to a cyclic one at n = 4-5, those of protonated 2,2,2-trifluoroethanol (TFE) clusters maintain linear type structures at least in the size range of n = 3-7. To explore the origin of the strong linear type network preference of H+(TFE)n, infrared spectra of protonated mixed clusters H+(TFE)m(ethanol)n (m + n = 5) were measured. An efficient structure sampling technique using parallelized basin-hopping algorithms and deep-learning neural network potentials is developed to search for essential isomers of the mixed clusters. Vibrational simulations based on the harmonic superposition approximation were compared with the observed spectra to identify the major isomer component at each mixing ratio. It was found that the formation of the cyclic structure occurs only in n ≥ 3 of the mixed clusters, in which the proton solvating sites and the double acceptor site are occupied by ethanol. The crucial role of the stability of the double acceptor site in the cyclic structure formation is discussed.

Proton quantal delocalization and H/D translocations in (MeOH)nH+ (n = 2, 3)

In this study, we present results from path integral molecular dynamics simulations that describe the characteristics of the quantum spatial delocalizations of protons participating in OH bonds in (MeOH)2H+ and in (MeOH)3H+. The characterization was carried out by examining the overall structures of the corresponding isomorphic polymers. To introduce full flexibility in the force treatment, we have adopted a neural network fitting procedure based on second-order Møller-Plesset perturbation theory predictions. For the dimer case, we found that the spatial extent of the shared connective proton can be portrayed in terms of a prolate-like structure with typical dimensions of ∼0.1 Å. On the other hand, the dangling polymers lie confined within a thin spherical layer, spread over length scales of the order of ∼0.25 Å. In contrast, connective protons in (MeOH)3H+ exhibit larger delocalizations along the O-H bond and more localized ones along perpendicular directions, compared to their dangling counterparts. We also examined the characteristics of the relative propensities of H and D isotopes to be localized in dangling and connective positions. Physical interpretations of the different thermodynamic trends are provided in terms of the local geometrical characteristics and of the strengths of the corresponding intermolecular connectivities.

Hydrogen bond network structures of protonated dimethylamine clusters H+(DMA)n (n = 3-7)

Infrared spectroscopy of protonated dimethylamine clusters, H+(DMA)n, (n = 3-7), and their Ar-tagged clusters was performed in the NH and CH stretching vibrational region to explore their hydrogen bond network structures. A stable isomer search and vibrational spectral simulations of the clusters were also carried out to support the interpretations of the observed spectra. Weakly hydrogen-bonded NH stretching vibrational bands, which are characteristic of cyclic structures of small-sized protonated clusters, are observed in the spectra of the Ar-tagged clusters of n ≥ 5, while only linear chain type structures are suggested for the Ar-tagged clusters of n = 3-4 and the bare clusters of all the sizes. These results demonstrate that the size and temperature dependence of the hydrogen bond network structures of the protonated dimethylamine clusters is analogous to that of protonated monohydric alcohol clusters.

Network Models of BACE-1 Inhibitors: Exploring Structural and Biochemical Relationships

This study investigates the clustering patterns of human β-secretase 1 (BACE-1) inhibitors using complex network methodologies based on various distance functions, including Euclidean, Tanimoto, Hamming, and Levenshtein distances. Molecular descriptor vectors such as molecular mass, Merck Molecular Force Field (MMFF) energy, Crippen partition coefficient (ClogP), Crippen molar refractivity (MR), eccentricity, Kappa indices, Synthetic Accessibility Score, Topological Polar Surface Area (TPSA), and 2D/3D autocorrelation entropies are employed to capture the diverse properties of these inhibitors. The Euclidean distance network demonstrates the most reliable clustering results, with strong agreement metrics and minimal information loss, indicating its robustness in capturing essential structural and physicochemical properties. Tanimoto and Hamming distance networks yield valuable clustering outcomes, albeit with moderate performance, while the Levenshtein distance network shows significant discrepancies. The analysis of eigenvector centrality across different networks identifies key inhibitors acting as hubs, which are likely critical in biochemical pathways. Community detection results highlight distinct clustering patterns, with well-defined communities providing insights into the functional and structural groupings of BACE-1 inhibitors. The study also conducts non-parametric tests, revealing significant differences in molecular descriptors, validating the clustering methodology. Despite its limitations, including reliance on specific descriptors and computational complexity, this study offers a comprehensive framework for understanding molecular interactions and guiding therapeutic interventions. Future research could integrate additional descriptors, advanced machine learning techniques, and dynamic network analysis to enhance clustering accuracy and applicability.

Experimental confirmation of the Badger-Bauer rule in the protonated methanol clusters: weak hydrogen bond formation as a measure of terminal OH acidity in hydrogen bond networks

We report a linear correlation between the OH stretch frequency shift of the protonated methanol cluster, H+(MeOH)n, upon the π-hydrogen bond formation with benzene and the enthalpy change in clustering of H+(MeOH)n to H+(MeOH)n+1. This result suggests a new method to explore hydrogen bond strength in hydrogen bond networks.

Infrared spectroscopy and theoretical structure analyses of protonated fluoroalcohol clusters: the impact of fluorination on the hydrogen bond networks

To explore the impact of fluorination on the hydrogen bond networks of protonated alkylalcohols, infrared spectroscopy and theoretical computations of protonated 2,2,2-trifluoroethanol clusters, H⁺(TFE)_n, (n = 4-7), were performed. It has been demonstrated that the development of the hydrogen bond networks from a linear type to cyclic types occurs in this size region for the protonated alkylalcohol clusters. In contrast, infrared spectroscopy of H⁺(TFE)_n in the OH/CH stretch region clearly indicated that the linear type structures are held in the whole size range, irrespective of temperature of the clusters. The extensive stable isomer structure search of H⁺(TFE)_n based on our latest sampling approach supported the strong preference of the linear type hydrogen bond networks. Detailed analyses of the free OH stretching vibrational bands evidenced the intra- and intermolecular OH⋯FC interactions in the clusters. In addition, infrared spectra of protonated clusters of 2,2-difluoroethanol, 2,2-difluoropropanol, and 3,3,3-trifluoropropanol were measured for n = 4 and 5, and their spectra also indicated the effective inhibition of the cyclic hydrogen bond network formation by the fluorination.

How many methanol molecules effectively solvate an excess proton in the gas phase? Infrared spectroscopy of H+(methanol)n-benzene clusters

An excess proton in a hydrogen-bonded system enhances the strength of hydrogen bonds of the surrounding molecules. The extent of this influence can be a measure of the number of molecules effectively solvating the excess proton. Such extent in methanol has been discussed by the observation of the π-hydrogen-bonded OH stretch bands of the terminal sites of protonated methanol clusters, H⁺(methanol)_n, in benzene solutions, and it has been concluded that ∼8 molecules effectively solvate the excess proton (Stoyanov et al., Chem. Eur. J. 2008, 14, 3596-3604). In the present study, we performed infrared spectroscopy of H⁺(methanol)_n-benzene clusters in the gas phase. The cluster size and hydrogen-bonded network structure are identified by the tandem mass spectrometric technique and the comparison of the observed infrared spectra with density functional theory calculations. Though changes of the preferred hydrogen bond network type occur with the increase of cluster size in the gas phase clusters, the observed size dependence of the π-hydrogen bonded OH frequency agrees well with that in the benzene solutions. This means that the observations in both the gas and condensed phases catch the same physical essence of the excess proton solvation by methanol.

Structures, binding energies and non-covalent interactions of furan clusters

Understanding of the furan solvent is subjected to the knowledge of the structures of the furan clusters and interactions taking place therein. Although, furan clusters can be very important to determine the dynamics and the properties of the furan solvent, there has been only a few investigations reported on furan dimer. In this work, we have explored the potential energy surfaces (PESs) of the furan clusters using two incremental levels of theory. Structures have been initially generated using classical molecular dynamics followed by full optimization at the MP2/aug-cc-pVDZ level of theory. The results show that the most stable structure of the furan dimer has a stacking configuration while that of the trimer has a cyclic configuration. We have noted that the structures of the furan tetramer have no definite configurations. In addition, we have performed a quantum theory of atoms in molecule (QTAIM) analysis to identify all possible non-covalent interactions of the furan clusters. The results show that six different types of non-covalent interactions can be identified in furan clusters. We have noted that the CH⋯C and CH⋯O hydrogen bondings are the strongest non-covalent interactions while the H⋯H bonding interaction is found to be the weakest. Furthermore, we have assessed the performance of ten DFT functionals in calculating the binding energies of the furan clusters. The ten DFT functionals (M05, M05-2X, M06, M06-2X, M08HX, PBE0, ωB97XD, PW6B95D3, APFD and MN15) have been benchmarked to DLPNO-CCSD(T)/CBS. The functionals M05-2X and M06 are recommended for further affordable investigations of the furan clusters.

Anharmonic Coupling Revealed by the Vibrational Spectra of Solvated Protonated Methanol: Fermi Resonance, Combination Bands, and Isotope Effect

Intriguing vibrational features of solvated protonated methanol between 2400-3800 cm^-1 are recorded by infrared predissociation spectroscopy. Positions of absorption bands corresponding to OH stretching modes are sensitive to changes in solvation environments, thus leading to changes in these vibrational features. Two anharmonic coupling mechanisms, Fermi resonance (FR) contributed by bending overtones and combination band (CB) associated with intermolecular stretching modes, are known to lead to band splitting of OH stretching fundamentals in solvated hydronium and ammonium. Theoretical analyses based on the ab initio anharmonic algorithm not only well reproduce the experimentally observed features but also elucidate the magnitudes of such couplings and the resulting interplay between these two mechanisms, which provide convincing assignments of the spectral patterns. Moreover, while the hydroxyl group plays the leading role in all the above-mentioned features, the role of the methyl group is also analyzed. Through the H/D isotope substitution, we identify overtones of the methyl-hydroxyl rocking modes and their participation in FR.

Effects of mixing between short-chain and branched-chain alcohols in protonated clusters

The previous analysis of the neat protonated branched-chain alcohol clusters revealed the impact of steric repulsion and dispersion of the bulky alkyl group on the hydrogen-bonded (H-bonded) structures and their temperature-dependence. To further understand the influence of the alkyl groups in H-bonded clusters, we studied the mixing of the two extremes of alcohols, methanol (MeOH) and tert-butyl alcohol (t-BuOH), with an excess proton. Infrared spectroscopy and a structural search of first principles calculations on the size-selected clusters H⁺(MeOH)_m(t-BuOH)_t (m + t = 4 and 5) were conducted. Temperature-dependence of the dominant H-bonded structures was explored by the Ar-tagging technique and quantum harmonic superposition approach. By introducing the dispersion-corrected density functional theory methods, it was shown that the effects of dispersion due to the bulky alkyl groups in the mixed clusters cannot be ignored for t≥ 2. The computational results qualitatively depicted the characteristics of the observed IR spectra, but overestimation of the temperature-dependence with dispersion correction was clearly seen due to the unbalanced correction between linear H-bonded structures and compact cyclic ones. These results demonstrate the importance of extensive investigation and benchmarks on different levels of theory, and that a properly sampled structure database is crucial to evaluate theoretical models.

Competition between hydrogen bonds and van der Waals forces in intermolecular structure formation of protonated branched-chain alcohol clusters

To investigate the influence of bulky alkyl groups on hydrogen-bonded (H-bonded) network structures of alcohols, infrared (IR) spectra of protonated clusters of 2-propanol (2-PrOH) and tert-butyl alcohol (t-BuOH) were observed in the OH and CH stretch regions. In addition, by varying the tag species, the temperature dependence profile of the isomer population of H+(t-BuOH)n was revealed. An extensive search for stable isomers was performed using dispersion-corrected density functional theory methods, and temperature-dependent IR spectral simulations were done on the basis of the harmonic superposition approximation. The computational results qualitatively agreed with the observed size and temperature dependence of the H-bonded network structures of these protonated bulky alcohol clusters. However, the difficulty in the quantitative evaluation of dispersion was also demonstrated. It was shown that H+(2-PrOH)n (n = 4-7) have essentially the same network structures as the protonated normal alcohol clusters studied so far. On the other hand, H+(t-BuOH)n (n = 4-8) showed a clear preference for the smaller-membered ring structures, that is very different from the preference of the protonated normal alcohol clusters. The origin of the different structure preferences was discussed in terms of the steric effect and dispersion.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

Influence of argon and D2 tagging on the hydrogen bond network in Cs+(H2O)3; kinetic trapping below 40 K

Tuning cluster ion conformations between 12 and 21 K.

Quantitative Temperature Dependence of the Microscopic Hydration Structures Investigated by Ultraviolet Photodissociation Spectroscopy of Hydrated Phenol Cations

To discuss the temperature effect on microscopic hydration structures in clusters, relative populations of the isomers having different hydration structures at well-defined temperatures are quite important. In the present study, we measured ultraviolet photodissociation spectra of the temperature-controlled hydrated phenol cation [PhOH(H₂O)₅]⁺ trapped in the 22-pole ion trap. Two isomers having a distinct hydration motif with each other are identified in the spectra, and a clear change in the relative populations is observed in the temperature range from 30 to 150 K. This behavior is quantitatively interpreted by statistical mechanical estimation based on density functional theory calculations. A ring with tail-type hydration motif is dominant in cold conditions, whereas a chain-like motif is dominant in hot conditions. The present study provides very quantitative information about the temperature effect on the microscopic hydration structures.

Structures of protonated hydrogen sulfide clusters, H+(H2S)n, highlighting the nature of sulfur-centered intermolecular interactions

Though H₂S has the same hydrogen bond coordination property as H₂O, intermolecular structures of H⁺(H₂S)_n are very different from those of H⁺(H₂O)_n, indicating the competition among hydrogen bond and other intermolecular interactions.

Exploration of hydrogen bond networks and potential energy surfaces of methanol clusters using a two-stage clustering algorithm

A two-stage algorithm based both on the similarity in shape and hydrogen bond network is developed to explore the potential energy surface of methanol clusters.

Structure prediction of the solid forms of methanol: an ab initio random structure searching approach

Liquid methanol and methanol clusters have been comprehensively studied to reveal their local structure and hydrogen bond networks.

Structures and spectroscopy of protonated ammonia clusters at different temperatures

Protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules.

Competing Insertion and External Binding Motifs in Hydrated Neurotransmitters: Infrared Spectra of Protonated Phenylethylamine Monohydrate

Hydration has a drastic impact on the structure and function of flexible biomolecules, such as aromatic ethylamino neurotransmitters. The structure of monohydrated protonated phenylethylamine (H(+) PEA-H2 O) is investigated by infrared photodissociation (IRPD) spectroscopy of cold cluster ions by using rare-gas (Rg=Ne and Ar) tagging and dispersion-corrected density functional theory calculations at the B3LYP-D3/aug-cc-pVTZ level. Monohydration of this prototypical neurotransmitter gives an insight into the first step of the formation of its solvation shell, especially regarding the competition between intra- and intermolecular interactions. The spectra of Rg-tagged H(+) PEA-H2 O reveal the presence of a stable insertion structure in which the water molecule is located between the positively charged ammonium group and the phenyl ring of H(+) PEA, acting both as a hydrogen bond acceptor (NH(+) ⋅⋅⋅O) and donor (OH⋅⋅⋅π). Two other nearly equivalent isomers, in which water is externally H bonded to one of the free NH groups, are also identified. The balance between insertion and external hydration strongly depends on temperature.

Theoretical vibrational spectra of OH(-)(H2O)2: the effect of quantum distribution and vibrational coupling

We performed ab initio path integral molecular dynamics simulations for the hydroxide-water cluster, OH(-)(H2O)2, at 50 K, 100 K, and 150 K to investigate its flexible structure. From our simulations, we found that nuclear quantum effects enhance hydroxide hydrogen atom inversion and the conformational change between isomers occurs by simultaneous rotation of the free hydrogen atom. We propose the importance of including the transition state conformer with C2 symmetry, for the description of this system at temperatures realized in predissociation experiments. Temperature dependence of relative populations of each conformer along with multidimensional vibrational calculations were used to simulate the vibrational spectra and compare with the experimental spectra of Johnson and coworkers. We assign the doublet peaks seen in the experiment at 2500 to 3000 cm(-1), as the mixture of the ionic hydrogen bonded OH stretching overtone, ionic hydrogen bonded OH bending overtone, and the combination band of the ionic hydrogen bonded OH stretch and bend, which are modulated by the van der Waals OO vibrations. We concluded that for OH(-)(H2O)2, the vibrational couplings between the ionic hydrogen bonded motion and floppy modes contribute to the broadening of peaks observed in the 2500 to 3000 cm(-1) region.

An infrared spectroscopic and theoretical study on (CH3)3N–H+–(H2O)n, n = 1–22: highly polarized hydrogen bond networks of hydrated clusters

Highly polarized water networks are found in the micro hydaration of protonated trimethylamine.

Hydrogen-bonded ring closing and opening of protonated methanol clusters H+(CH3OH)n (n = 4–8) with the inert gas tagging

Temperature dependence of hydrogen bond network structures of protonated methanol clusters is explored by IR spectroscopy and DFT simulations.

Infrared spectroscopy of large-sized neutral and protonated ammonia clusters

Size selective IR spectroscopy shows the nature of hydrogen bond networks in neutral and protonated ammonia clusters.