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.

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.

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)2 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)2 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)2 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-d6), 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)2 and (pyrazine)(benzene) in a supersonic jet. For (pyrazine)2, 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)2, 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.

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.

Exploring the structure and properties of molecular clusters with accuracy using the ab initio methods is a resource intensive task due to the increasing cost of the ab initio methods and the number of distinct conformers as the size increases. The energy landscape of methanol clusters has been previously explored using computationally efficient empirical models to collect a database of structurally distinct minima, followed by re-optimization using ab initio methods. In this work, 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. Picking 684 local minima of (CH3OH)5 to (CH3OH)8 from the existing database, we can generate ∼51 000 data points of one-body, two-body, three-body and four-body molecular systems to train an ML model to reach a mean absolute error (MAE) of 3.19 kJ mol-1 (in energy) and 2.48 kJ mol-1 Å-1 (in forces) tested against ab initio calculations up to (CH3OH)14. This ML model is then used to create a database of low energy isomers of (CH3OH)n (n = 15-20). The proposed scheme can be applied to other hydrogen bonded molecular clusters with an accuracy of first-principles methods and computational speed of empirical force-fields.

Collision-induced dissociation of Na+-tagged ketohexoses: experimental and computational studies on fructose.

Collision-induced dissociation tandem mass spectrometry (CID-MSn) and computational investigation at the MP2/6-311+G(d,p) level of theory have been employed to study Na+-tagged fructose, an example of a ketohexose featuring four cyclic isomers: α-fructofuranose (αFruf), ß-fructofuranose (ßFruf), α-fructopyranose (αFrup), and ß-fructopyranose (ßFrup). The four isomers can be separated by high-performance liquid chromatography (HPLC) and they show different mass spectra, indicating that CID-MSn can distinguish the different fructose forms. Based on a simulation using a micro-kinetic model, we have obtained an overview of the mechanisms for the different dissociation pathways. It has been demonstrated that the preference for the C-C cleavage over the competing isomerization of linear fructose is the main reason for the previously reported differences between the CID-MS spectra of aldohexoses and ketohexoses. In addition, the kinetic modeling helped to confirm the assignment of the different measured mass spectra to the different fructose isomers. The previously reported assignment based on the peak intensities in the HPLC chromatogram had left some open questions as the preference for the dehydration channels did not always follow trends previously observed for aldohexoses. Setting up the kinetic model further enabled us to directly compare the computational and experimental results, which indicated that the model can reproduce most trends in the differences between the dissociation pathways of the four cyclic fructose isomers.

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.

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.

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

Collision-induced dissociation (CID) of α-xylose and ß-xylose were studied using mass spectrometry and quantum chemistry calculations. Three dissociation channels, namely loss of metal ions, dehydration, and cross-ring dissociation were found. The major dissociation channel of sodium adducts is the loss of sodium ions, and the minor dissociation channels are dehydration and cross-ring dissociation. By contrast, dehydration and cross-ring dissociation are the major dissociation channels of lithium adducts, and the corresponding dissociation mechanisms can be used to determine the anomericity and linkages of xylose in oligosaccharides. These mechanisms include (1) the dehydration branching ratio can be used to differentiate the anomericity of xylose and xylose in oligosaccharides because α-xylose has a larger branching ratio of dehydration than ß-xylose, (2) various cross-ring dissociation reactions can be used to identify linkage positions. The oligosaccharide with xylose at the reducing end is predicted to undergo 0,2X, 0,3X, and 0,2A cross-ring dissociation for the 1 → 2, 1 → 3, and 1 → 4 linkages, respectively. Application of these mechanisms to determine the anomericity and linkage positions of xylobiose was demonstrated.

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 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(H2O)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(H2O)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(H2O)n upon the VUV ionization, which was followed by the elimination of water molecules. Among the recorded IR spectra, only one spectrum monitored, (Pyd)2+ cation, showed a well-resolved structure. The spectrum was analyzed by comparing with the simulated ones of possible stable isomers of (Pyd)2(H2O)n, which were obtained with quantum-chemical calculations. Most of the calculated (Pyd)2(H2O)n clusters had the characteristic structure in which H2O or (H2O)2 forms a hydrogen-bonded bridge between two pyridines to form the π-stacked (Pyd)2, and an additional H2O molecule(s) extends the H-bonded network. The π-stacked (Pyd)2(H2O)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(H2O)n was found to be less pronounced in the case of (Pyd)m(NH3)n studied previously.

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.

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.

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

The aggregation of phenylacetylene in the gas phase was investigated by selectively recording the IR spectra of clusters consisting of up to six monomer units. Analysis of the IR spectra with the aid of B97-D3/aug-cc-pVDZ level calculations reveals the formation of an anti-parallel π-stacked structure of the dimer and a hitherto unknown assembly of clusters incorporating exclusively aromatic C-Hπ interactions between various units of the trimer and higher clusters. The aggregation behaviour of phenylacetylene in the gas phase is fundamentally different from benzene, phenol and aniline vis-à-vis their crystal structures. The structures of the three known polymorphic crystals can be reconciled by the formation of supramolecular synthons with acetylenic C-Hπ interactions, which is preferred over energetically favored aromatic C-Hπ interactions. Furthermore, the small (phenylacetylene)n [n = 3-6] clusters, the structures incorporating aromatic C-Hπ interactions, can be envisaged as liquid-like aggregates which under varied conditions lead to the formation of multiple polymorphs during in situ cryo-crystallization.

Unexpected Dissociation Mechanism of Sodiated N-Acetylglucosamine and N-Acetylgalactosamine.

The mechanism for the collision-induced dissociation (CID) of two sodiated N-acetylhexosamines (HexNAc), N-acetylglucosamine (GlcNAc), and N-acetylgalactosamine (GalNAc), was studied using quantum-chemistry calculations and resonance excitation in a low-pressure linear ion trap. Experimental results show that the major dissociation channel of the isotope labeled [1-18O, D5]-HexNAc is the dehydration by eliminating HDO, where OD comes from the OD group at C3. Dissociation channels of minor importance include the 0,2A cross-ring dissociation. No difference has been observed between the CID spectra of the α- and ß-anomers of the same HexNAc. At variance, the CID spectra of GlcNAc and GalNAc showed some differences, which can be used to distinguish the two structures. It was observed in CID experiments involving disaccharides with a HexNAc at the nonreducing end that a ß-HexNAc shows a larger dissociation branching ratio for the glycosidic bond cleavage than the α-anomer. This finding can be exploited for the rapid identification of the anomeric configuration at the glycosidic bond of HexNAc-R' (R' = sugar) structures. The experimental observations indicating that the dissociation mechanisms of HexNAcs are significantly different from those of hexoses were explained by quantum-chemistry calculations. Calculations show that ring opening is the major channel for HexNAcs in a ring form. After ring opening, dehydration shows the lowest barrier. In contrast, the glycosidic bond cleavage becomes the major channel for HexNAcs at the nonreducing end of a disaccharide. This reaction has a lower barrier for ß-HexNAcs as compared with the barrier of the corresponding α-anomers, consistent with the higher branching ratio for ß-HexNAcs observed in experiment.

Toward Closing the Gap between Hexoses and N-Acetlyhexosamines: Experimental and Computational Studies on the Collision-Induced Dissociation of Hexosamines.

Motivated by the fundamental difference in the reactivity of hexoses and N-acetylhexosamines under collision-induced dissociation (CID) mass spectrometry conditions, we have investigated the CID of two hexosamines, glucosamine (GlcN) and galactosamine (GalN), experimentally and computationally. Both hexosamines undergo ring-opening and then dissociate via the 0,2A and the 0,3A (0,3X) cross-ring cleavage channels. The preference for the ring-opening is similar to the behavior of N-acetylhexosamines and explains why the two anomers of the same sugar give the same mass spectrum. While the spectrum for GlcN is dominated by the 0,2A signal, the signal intensities for both 0,2A and the 0,3A (0,3X) dissociation channels are comparable for GalN, which allows GlcN and GalN to be distinguished easily. Calculations at MP2 level of theory indicate that this is related to the differences in the relative barrier heights for the 0,2A and the 0,3A (0,3X) cross-ring cleavage channels. This, in return, reflects the circumstance that the 0,2A cross-ring cleavage barriers are different for the two sugars, while the barriers of all other dissociation channels are comparable. While the mechanisms of the cross-ring dissociation channels of hexoses are well described using the retro-aldol mechanism in the literature, this study proposes a new mechanism for the 0,3A (0,3X) cross-ring cleavage of hexosamines that involves the formation of an epoxy intermediate or a zwitterionic intermediate.

Infrared spectra of neutral dimethylamine clusters: An infrared-vacuum ultraviolet spectroscopic and anharmonic vibrational calculation study.

Infrared-vacuum ultraviolet (IR-VUV) spectra of neutral dimethylamine clusters, (DMA)n (n = 2-5), were measured in the spectral range of 2600-3700 cm-1. The experimental IR-VUV spectra show NH stretch modes gradually redshift to 3200-3250 cm-1 with the increase in the cluster size and complex Fermi Resonance (FR) pattern of the CH3 group in the 2800-3000 cm-1 region. Ab initio anharmonic vibrational calculations were performed on low-energy conformers of (DMA)2 and (DMA)3 to examine vibrational coupling among CH/NH and to understand the Fermi resonance pattern in the observed spectra features. We found that the redshift of NH stretching mode with the size of DMA cluster is moderate, and the overtone of NH bending modes is expected to overlap in frequency with the CH stretching fundamental modes. The FR in CH3 groups is originated from the strong coupling between CH stretching fundamental and bending overtone within a CH3 group. Well-resolved experimental spectra also enable us to compare the performance of ab initio anharmonic algorithms at different levels.

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.

Because of the hydrogen bond coordination properties of alcohols, their possible hydrogen bond network structures are categorized into only a few types. Therefore, gas phase clusters of alcohols can be a very simple model system to examine the properties of hydrogen bond networks, such as structural development with cluster size and temperature dependence. In this perspective, we focus on the structural study of protonated short-chain alcohol clusters, whose excess protons (charge) enable size-selective spectroscopy in combination with mass spectrometric techniques. Size-selective infrared spectroscopy and a theoretical multi-scale isomer search were applied to protonated clusters of methanol, which is a prototype of short-chain alcohols, and their hydrogen bond network development is elucidated in detail. Complete isomer population switching with increasing temperature was predicted by the quantum harmonic superposition approximation and this isomer switching was evidenced by the remarkable temperature (internal vibrational energy) dependence of the observed infrared spectra. The characteristics of the temperature dependence of protonated methanol were compared with those of water and neutral methanol. In addition, possible hydrogen bond networks of methanolated ions were discussed on the basis of the results for protonated methanol. Stepwise changes in the internal energy of clusters with inert gas tagging are demonstrated. Convergence of the hydrogen bond network to the bulk-like network in large clusters is also discussed. The hydrogen bond structures of the protonated clusters of longer normal alkyl chain alcohols (ethanol, 1-propanol, 1-butanol, and 1-pentanol) are determined by comparison of their infrared spectra with those of the protonated methanol clusters. It is demonstrated that the normal alkyl chain interferes only slightly with the most stable hydrogen bond structure, although a few exceptional cases were also found. These exception cases serve as good model systems for further theoretical and computational studies.