Searching low-energy conformers of neutral and protonated di-, tri-, and tetra-glycine using first-principles accuracy assisted by the use of neural network potentials

In the last ten years, combinations of state-of-the-art gas-phase spectroscopies and quantum chemistry calculations have suggested several intuitive trends in the structure of small polypeptides that may not hold true. For example, the preference for the cis form of the peptide bond and multiple protonated sites was proposed by comparing experimental spectra with low-energy minima obtained from limited structural sampling using various density functional theory methods. For understanding the structures of polypeptides, extensive sampling of their configurational space with high-accuracy computational methods is required. In this work, we demonstrated the use of deep-learning neural network potential (DL-NNP) to assist in exploring the structure and energy landscape of di-, tri-, and tetra-glycine with the accuracy of high-level quantum chemistry methods, and low-energy conformers of small polypeptides can be efficiently located. We hope that the structures of these polypeptides we found and our preliminary analysis will stimulate further experimental investigations.

A first-principles exploration of the conformational space of sodiated di-saccharides assisted by semi-empirical methods and neural network potentials

Previous exploration of the conformational space of sodiated mono-saccharides using a random search algorithm leads to ∼103 structurally distinct conformers covering an energy range of ∼150 kJ mol-1. Thus, it is reasonable to expect that the number of distinct conformers for a given disaccharide would be on the order of 106. Efficient identification of distinct conformers at the first-principles level has been demonstrated with the assistance of neural network potential (NNP) with an accuracy of ∼1 kJ mol-1 compared to DFT. Leveraging a local minima database of neutral and sodiated glucose (Glc), we develop algorithms to systematically explore the conformation landscape of 19 Glc-based sodiated disaccharides. To accelerate the exploration, the NNP method is implemented. The NNP achieves an accuracy of ∼2.3 kJ mol-1 compared to DFT, offering a comparable quality to that of DFT. Through a multi-model approach integrating DFTB3, NNP and DFT, we can rapidly locate low-energy disaccharide conformers at the first-principles level. The methodology we show here can be used to efficiently explore the potential energy landscape of any di-saccharides when first-principles accuracy is required.

Structures of (Pyrazine)2 and (Pyrazine)(Benzene) Dimers Investigated with Infrared-Vacuum Ultraviolet Spectroscopy and Quantum-Chemical Calculations: Competition among π-π, CH···π, and CH···N Interactions

The structures of a pyrazine dimer (pyrazine)₂ and (pyrazine)(benzene) hetero-dimer cooled in a supersonic beam were investigated by the measurement of the infrared spectra in the C-H stretching region with infrared-vacuum ultraviolet (IR-VUV) spectroscopy and quantum-chemical calculations. The stabilization energy calculation at the CCSD(T)/aug-cc-pVTZ level of theory predicted three isomers for (pyrazine)₂ and three for (pyrazine)(benzene) with energy within 6 kJ/mol. Among them, the cross-displaced π-π stacked structure is the most stable in both dimers. In the observed IR spectra, both dimers exhibited two intense bands near 3065 cm^-1, with intervals of 8 cm^-1 in (pyrazine)₂ and 11 cm^-1 in (pyrazine)(benzene), while only one band appeared in the monomer. For (pyrazine)(benzene), we also measured the IR spectrum of (pyrazine)(benzene-d₆), where the interval of the two bands was unchanged. The analysis of the observed IR spectra with anharmonic calculations suggested the coexistence of three isomers of (pyrazine)₂ and (pyrazine)(benzene) in a supersonic jet. For (pyrazine)₂, the two isomers which were previously assigned to the H-bonded planar and the π-π stacked structures respectively were reassigned to the cross-displaced π-π stacked and T-shaped structures, respectively. In addition, the quantum chemical calculation and IR-VUV spectral measurement suggested the coexistence of the H-bonded planar isomer in the jet. For (pyrazine)(benzene), the IR spectrum of the (pyrazine) site showed a similar spectral pattern to that of (pyrazine)₂, especially the split at ∼3065 cm^-1. However, the anharmonic analysis suggested that they are assigned to the different vibrational motions of (pyrazine). The anharmonic vibrational analysis is essential to associate the observed IR spectra with the correct structures of the dimer.

A first-principles exploration of the conformational space of sodiated pyranose assisted by neural network potentials

Sampling the conformational space of monosaccharides using the first-principles methods is important and as a database of local minima provides a solid base for interpreting experimental measurements such as infrared photo-dissociation (IRPD) spectroscopy or collision-induced dissociation (CID). IRPD emphasizes low-energy conformers and CID can distinguish conformers with distinct reaction pathways. A typical computational approach is to engage empirical or semi-empirical methods to sample the conformational space first, and only selected minima are reoptimized at first-principles levels. In this work, we propose a computational scheme to explore the configurational space of 12 types of sodiated pyranoses with the assistance of a neural network potential (NNP). We demonstrated that it is possible to train an NNP based on the density functional calculations extracted from a previous study on sodiated glucose (Glc), galactose (Gal), and mannose (Man). This NNP yields a better description of the other five types of aldohexoses than the four types of ketohexoses. We further show that such a discrepancy in the accuracy of NNP can be resolved by an active learning scheme where the NNP model is engaged in generating the data and has itself updated. Through this iterative process, we can locate more than 17 000 distinct local minima at the B3LYP/6-311+G(d,p) level and an NNP with an accuracy of 1 kJ mol^-1 was created, which can be used for further studies.

A self-adapting first-principles exploration on the dissociation mechanism in sodiated aldohexose pyranoses assisted with neural network potentials

Understanding the mechanism of collision-induced dissociation (CID) in mono-saccharides with density functional theory (DFT) is challenging because of many possible reaction paths that originate from their high structural diversity. To search for the transition state (TS) from the huge number of conformers, we propose a three-step search scheme with the assistance of neural network potential (NNP). The search starts from a cross-checking of sugars, to a global search of all possible channels, and in the end, an exhaustive exploration around the low-lying channels. The cross-checking step quickly adapts the NNP from the studied molecules to the target ones. The other two steps utilize the adapted NNP to find the available pathways via random sampling of the structures. The study of the CID reactions in all eight types of aldohexose pyranoses was applied using the search scheme. The DFT calculations on AH-0 (Glc, Gal, and Man) in the previous study were utilized to construct an NNP and provide the TS structure database for searching AH-1 (All, Alt, Gul, Ido, and Tal). In total, we identified around 5200 TSs in AH-0 and AH-1, and the final NNP covers an energy range of more than 500 kJ mol^-1 with a mean absolute error of energy less than 4 kJ mol^-1. The search scheme is useful not only for saccharides but also for highly flexible bio-molecules.

Hydrogen bond networks of dimethylsulfoxide (DMSO) pentamer

Understanding of clusters of dimethylsulfoxide (DMSO) is important in several applications in Chemistry. Despite its importance, very few studies of DMSO clusters, (DMSO)_n, have been reported in comparison to systems such as water clusters or methanol clusters. In order to provide further understanding of DMSO clusters, we investigated the structures and non-covalent interactions of the (DMSO)_n, n=5. Therefore, the potential energy surface (PES) of the DMSO pentamer has been examined using classical molecular dynamics. The structures generated using classical molecular dynamics are further optimized at the PW6B95D3/aug-cc-pVDZ level of theory. To comprehend the non-covalent bondings in the DMSO pentamer, we carried out a quantum theory of atoms in molecule (QTAIM) analysis. In addition, the effects of temperature on the structural stability is investigated between 20 and 500K. It comes out that seven different kind of non-covalent bondings can be found in DMSO pentamers.

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.

Capturing the potential energy landscape of large size molecular clusters from atomic interactions up to a 4-body system using deep learning

We propose a new method that utilizes the database of stable conformers and borrow the fragmentation concept of many-body-expansion (MBE) methods in ab initio methods to train a deep-learning machine learning (ML) model using SchNet.

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.

Is Dissociation of HCl in DMSO Clusters Bistable?

Dissociation of HCl embedded in dimethyl sulfoxide (DMSO) clusters was investigated by projecting the solvent electric field along the HCl bond using B3LYP-D3/6-31+G(d) and MP2/6-31+G(d,p) levels of theory. A large number of distinct structures (about 1500) consisting of up to five DMSO molecules were considered in the present work for statistical reliability. The B3LYP-D3 calculations reveal that the dissociation of HCl embedded in DMSO clusters requires a critical electric field of 138 MV cm^-1 along the H-Cl bond. However, a large number of exceptions wherein the electric field values much higher than the critical electric field of 138 MV cm^-1 did not result in dissociation of HCl were observed, in addition to several cases wherein the HCl dissociates with an electric field less than the critical electric field. On the other hand, the MP2 level calculations reveal that the critical electric field for HCl dissociation is about 181 MV cm^-1 with almost no exceptions. A comparison of calculations carried out using the MP2 and the B3LYP-D3 levels suggests that the dissociation of HCl embedded in DMSO clusters is bistable at the B3LYP-D3 level, which is an artifact, suggesting that care must be exercised in interpreting the processes of proton transfer. The answer to the question raised as the title of this paper is NO.

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.

Structures of Pyridine-Water Clusters Studied with Infrared-Vacuum Ultraviolet Spectroscopy

The infrared (IR) spectra of the O-H stretching vibrations of pyridine-water clusters (Pyd)_m(H₂O)_n, with m, n = 1-4, have been investigated with infrared-vacuum ultraviolet (VUV) spectroscopy under a jet-cooled condition. The time-of-flight mass spectrum of (Pyd)_m(H₂O)_n⁺ by VUV ionization at ∼9 eV showed an unusual intensity pattern with very weak ion signals for m = 1 and 2 and stronger signals for m ≥ 3. This unusual mass pattern was explained by a drastic structural change of (Pyd)_m(H₂O)_n upon the VUV ionization, which was followed by the elimination of water molecules. Among the recorded IR spectra, only one spectrum monitored, (Pyd)₂⁺ cation, showed a well-resolved structure. The spectrum was analyzed by comparing with the simulated ones of possible stable isomers of (Pyd)₂(H₂O)_n, which were obtained with quantum-chemical calculations. Most of the calculated (Pyd)₂(H₂O)_n clusters had the characteristic structure in which H₂O or (H₂O)₂ forms a hydrogen-bonded bridge between two pyridines to form the π-stacked (Pyd)₂, and an additional H₂O molecule(s) extends the H-bonded network. The π-stacked (Pyd)₂(H₂O)_n moiety is very stable and is thought to exist as a local structure in a pyridine/water mixed solution. The Fermi resonance between the O-H stretch fundamentals and the overtones of the O-H bending vibrations in (Pyd)_m(H₂O)_n was found to be less pronounced in the case of (Pyd)_m(NH₃)_n studied previously.

Identification of Anomericity and Linkage of Arabinose and Ribose through Collision-Induced Dissociation

Arabinose and ribose are two common pentoses that exist in both furanose and pyranose forms in plant and bacteria oligosaccharides. In this study, each pentose isomer, namely α-furanose, β-furanose, α-pyranose, and β-pyranose, was first separated through high-performance liquid chromatography followed by an investigation of collision-induced dissociation in an ion trap mass spectrometer. The major dissociation channels, dehydration and cross-ring dissociation, were analyzed by using high-level quantum chemistry calculations and transition state theory. The branching ratio of major dissociation channels was governed by two geometrical features: one being the cis or trans configuration of O1 and O2 atoms determining dehydration preferability and the other being the number of hydroxyl groups on the same side of the ring as the O1 atom determining the preferability of cross-ring dissociation. The relative branching ratios of the major channels were used to identify anomericity and the linkages of arabinose and ribose. Arabinose in the β-configuration and ribose in the α-configuration are predicted to have larger relative dehydration branching ratios than arabinose in the α-configuration and ribose in the β-configuration, respectively. Arabinose and ribose at the reducing end of oligosaccharides with 1 → 2 (pyranose and furanose), 1 → 3 (pyranose and furanose), 1 → 4 (pyranose only), and 1 → 5 (furanose only) linkages are predicted to undergo ^0,2X, ^0,3X, ^0,2A, and ^0,2A/^0,3A cross-ring dissociation, respectively. Application of the dissociation mechanism to the disaccharide linkage determination is demonstrated.

Disentangling the Complex Vibrational Spectra of Hydrogen-Bonded Clusters of 2-Pyridone with Ab Initio Structural Search and Anharmonic Analysis

Vibrational spectra of 1:1 clusters of 2-pyridone (2PY) with water, ammonia, and other hydrogen bond-forming molecules have been measured by several experimental groups over the past two decades. Complex vibrational signatures associated with the N-H stretching fundamental at 3 μm are often observed. Several anharmonic coupling schemes have been proposed; however, the origin of these commonly seen complex features remains unclear. In this work, we present our theoretical analysis on the structure and vibrational spectra of these clusters using ab initio random search and ab initio anharmonic algorithms, respectively. Low-energy conformers were found to be hydrogen-bonded clusters and their vibrational spectra at 3 μm were simulated with ab initio anharmonic algorithms. We demonstrate that simple anharmonic mechanisms of Fermi resonance (FR), coupling between NH stretching modes, and overtone/combinations of skeleton modes of 2PY can lead to the complex vibrational signatures observed experimentally. Since this vibrational coupling scheme is inherent to the cis-amides with adjacent N-H and C═O groups when a hydrogen bond is formed with 2PY as the donor and acceptor, we believe that such a phenomenon is general to other hydrogen-bonding systems with the same functional group.

Dipole moment enhanced π-π stacking in fluorophenylacetylenes is carried over from gas-phase dimers to crystal structures propagated through liquid like clusters

The aggregates of monofluorinated phenylacetylenes in the gas-phase, investigated using the IR-UV double resonance spectroscopic method in combination with extensive structural search and electronic structure calculations, reveal the formation of liquid-like clusters with a π-stacked dimeric core. The structural assignment based on the IR spectra in the acetylenic and aromatic C-H stretching regions suggests that, unlike the parent non-fluorinated phenylacetylene, the substitution of a F atom on the phenyl ring increases the dipole moment, leading to robustness in the formation of a ππ stacked dimer, which propagates incorporating C-Hπ_{Ar/Ac} and C-HF interactions involving both acetylenic and aromatic C-H groups. The structural evolution of fluorophenylacetylene aggregates in the gas phase shows marginal effects due to fluorine atom position on the phenyl ring, with substitution in the para-position tending towards phenylacetylene. The present study signifies that the ππ stacked dimers act as a nucleus for the growth of higher clusters to which other molecular units are added predominantly via the {Ar}_C-Hπ_{Ar} type of interaction and the dominant interactions present in the crystal structures gradually emerge with increasing cluster size. Based on these features, gas-phase clusters of fluorophenylacetylene are hypothesized as "liquid-like clusters" acting as intermediates in the generation of various polymorphic forms starting from a ππ stacked dimer as the core molecular unit.

Size of the hydrogen bond network in liquid methanol: a quantum cluster equilibrium model with extensive structure search

Studies have debated what is a favorable cluster size in liquid methanol. Applications of the quantum cluster equilibrium (QCE) model on a limited set of cluster structures have demonstrated the dominance of cyclic hexamers in liquid methanol. In this study, we examined the aforementioned question by integrating our implementation of QCE with a molecular-dynamics-based structural searching scheme. QCE simulations were performed using a database comprising extensively searched stable conformers of (MeOH)n for n = 2-14, which were optimized by B3LYP/6-31+G(d,p) with and without the dispersion correction. Our analysis indicated that an octamer structure can contribute significantly to cluster probability. By reoptimizing selected conformers with high probability at the MP2 level, we found that the aforementioned octamer became the dominant species due to favorable vibrational free energy, which was attributed to modes of intermolecular vibration.

Collision-induced dissociation of xylose and its applications in linkage and anomericity identification

Different dehydration barrier heights result in different branching ratio, a simple and fast anomeric configuration identification for xylose.

IR-VUV spectroscopy of pyridine dimers, trimers and pyridine-ammonia complexes in a supersonic jet

The infrared spectra of the C-H stretching vibrations of (pyridine)m, m = 1-3, and the N-H stretching vibrations of (pyridine)m-(NH3)n, m = 1, 2; n = 1-4, complexes were investigated by infrared (IR)-vacuum ultraviolet (VUV) spectroscopy under jet-cooled conditions. The ionization potential (IP0) of the pyridine monomer was determined to be 74 546 cm-1 (9.242 eV), while its complexes showed only smooth curves of the ionization thresholds at ∼9 eV, indicating large structural changes in the ionic form. The pyridine monomer exhibits five main features with several satellite bands in the C-H stretching region at 3000-3200 cm-1. Anharmonic calculations including Fermi-resonance were carried out to analyze the candidates of the overtone and combination bands which can couple to the C-H stretching fundamentals. For (pyridine)2 and (pyridine)3, most C-H bands are blue-shifted by 3-5 cm-1 from those of the monomer. The structures revealed by random searching algorithms with density functional methods indicate that the π-stacked structure is most stable for (pyridine)2, while (pyridine)3 prefers the structures stabilized by dipole-dipole and C-Hπ interactions. For the (pyridine)m-(NH3)n complexes, the mass spectrum exhibited a wide range distribution of the complexes. The observed IR spectra in the N-H stretching vibrations of the complexes showed four main bands in the 3200-3450 cm-1 region. These features are very similar to those of (NH3)n complexes, and the bands are assigned to the anti-symmetric N-H stretching band (ν3), the symmetric N-H stretching (ν1) band, and the first overtone bands of the N-H bending vibrations (2ν4). The anharmonic calculations including the Fermi-resonance between ν1 and 2ν4 well reproduced the observed spectra.

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.

Exploration of the potential energy surfaces of small ethanol clusters

The potential energy surfaces of small ethanol clusters, from dimer to pentamer, have been thoroughly explored using two different levels of theory. There is a clear relative energy gap between cyclic, linear and branched cyclic structures.

Anharmonic coupling behind vibrational spectra of solvated ammonium: lighting up overtone states by Fermi resonance through tuning solvation environments

Fascinating Fermi resonance bands emerge from anharmonic couplings between NH stretching fundamentals and bending overtones in ammonium-centered clusters.

The interplay and strength of the π⋯HF, C⋯HF, F⋯HF and F⋯HC hydrogen bonds upon the formation of multimolecular complexes based on C2H2⋯HF and C2H4⋯HF small dimers

The conception of this theoretical research was idealized aiming to unveil the intermolecular structures of complexes formed by acetylene or ethylene and hydrofluoric acid. At light of computational calculations by using the B3LYP/6-311++G(d,p) method, the geometries of the C₂H₂⋯(HF), C₂H₂⋯2(HF), C₂H₂⋯4(HF), C₂H₄⋯(HF), C₂H₄⋯2(HF) and C₂H₄⋯4(HF) hydrogen-bonded complexes were fully optimized. Moreover, the Post-Hartree-Fock calculations MP2/6-311++G(d,p), MP2/aug-cc-pVTZ, MP4(SDQ)/6-311++G(d,p) and CCSD/6-311++G(d,p) also were also used. The infrared spectra were analyzed in order to identify the new vibrational modes and frequencies of the proton donors shifted to red region. Through the modeling of charge-fluxes on the basis of the Quantum Theory of Atoms In Molecules (QTAIM) and, by contradicting the expectation of the hydrofluorination mechanisms of acetylene or ethylene, C⋯HF was recognized as a new type of hydrogen bond instead of the already well known π⋯H. The calculations of the Natural Bonding Orbital (NBO) and Charges derived from the Electrostatic Potential Grid-based (ChElPG) were also applied to interpret the shifting frequencies as well as measuring of the punctual charge-transfer after the formation of the complexes. Finally, the determination of the stabilization energy was carried out through the arguments of the Fock matrix in NBO basis and through the supermolecule approach. Also it is worthwhile to notice that some algebraic formulations were used for determining the electronic cooperative effect (CE).

A liquid crucible model for aggregation of phenylacetylene in the gas phase

Structural transformation from a π-stacked dimer to an aromatic C–H⋯π trimer and a tetramer.

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.

Vibrational spectra of small methylamine clusters accessed by an ab initio anharmonic approach

Anharmonic vibrational calculations on Methylamine (MMA) clusters suggest that the origin of the complexity between 2800 and 3000 cm^–1 is caused by Fermi resonance (FR) between the stretching and bending overtones of the CH₃ group. In trimer and tetramer, FR also causes complex spectra pattern in the NH₂ group.

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.

Collision-induced dissociation of sodiated glucose, galactose, and mannose, and the identification of anomeric configurations

Different dehydration barrier heights of cis and trans configurations between O1 and O2 provide a simple and fast anomeric configuration identification.

Solvation energies of the proton in methanol revisited and temperature effects

Various functionals assessing solvation free energies and enthalpies of the proton in methanol.