Water Clusters in

Jan 6, 2017 - Infrared Spectroscopy of Methanol and Methanol/Water Clusters in Helium Nanodroplets: The OH Stretching Region. Media I. Sulaiman, Sheng...
2 downloads 14 Views 1MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Infrared Spectroscopy of Methanol and Methanol/Water Clusters in Helium Nanodroplets: The OH Stretching Region Media I. Sulaiman, Shengfu Yang, and Andrew Michael Ellis J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.6b11170 • Publication Date (Web): 06 Jan 2017 Downloaded from http://pubs.acs.org on January 10, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Infrared Spectroscopy of Methanol and Methanol/Water Clusters in Helium Nanodroplets: the OH Stretching Region

Media I. Sulaiman, Shengfu Yang* and Andrew M. Ellis*

Department of Chemistry, University of Leicester, University Road, Leicester, LE1 7RH, UK

Email: [email protected]; [email protected]

Manuscript submitted to Journal of Physical Chemistry A

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 28

Abstract Infrared (IR) spectra of methanol clusters in helium nanodroplets are reported in the OH stretching region for the first time. A simple series of intense bands are seen which almost perfectly match previous gas phase studies of these clusters and which are consistent with cyclic structures for the trimer and larger clusters. This finding differs from an earlier report of (CH3OH)n clusters in helium nanodroplets, which focused on the CO stretching region and concluded that while the trimer was cyclic, the tetramer and pentamer adopted branched structures based on a cyclic trimer core. We also present preliminary data for small (CH3OH)n(H2O) clusters, and in particular we report the first IR spectra for (CH3OH)2(H2O) and (CH3OH)3(H2O). Supporting ab initio calculations suggest that, like the pure methanol clusters, cyclic structures are adopted by these mixed solvent clusters in helium droplets.

2 ACS Paragon Plus Environment

Page 3 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction The molecules in liquid or solid methanol are held together by a mixture of hydrogen bonding and van der Waals interactions. A detailed exploration of the nature of these intermolecular forces and how they evolve as molecules are progressively added is possible through the study of clusters of methanol. There have been many spectroscopic studies of methanol clusters, particularly in the gas phase1-14 but also in inert matrices at low temperature.15-19 The experimental work, supported by theoretical predictions,20-26 shows that the smallest cluster, the dimer, exists as an open structure with a single hydrogen bond holding the two molecules together: one methanol molecule acts as a proton donor and the other as an acceptor. For the trimer and larger clusters the adoption of cyclic structures maximizes the number of hydrogen bonds and is energetically preferred over chain-like structures. IR and Raman spectra of these clusters in supersonically-cooled gas expansions and molecular beams have found clear signals for the cyclic trimer and tetramer in the OH stretching region, while the signals for the pentamer and larger clusters seem to merge into a single broad OH stretching band.2,3,5-14 No evidence has been found for non-cyclic methanol clusters larger than the dimer in the gas phase. Preferential production of cyclic isomers for the trimer and larger clusters has also been noted in experiments using solid nitrogen and argon matrices.15,19 Liquid helium nanodroplets provide an alternative matrix environment in which to study molecules and their clusters. Characteristics of helium nanodroplets as a matrix in which to embed molecules include a low steady-state temperature (0.4 K), rapid cooling of dopants by evaporative loss of helium atoms from the surface of the superfluid helium, and very small perturbations to the vibrational spectra of the dopants by the surrounding helium.27-29 In essence, the helium droplets provide a cold, gas-like, environment. Methanol and its clusters were among the earliest of molecules studied in helium nanodroplets. IR

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 28

spectra of small methanol clusters were recorded by Behrens et al. in the CO stretching region.30 The spectra were recorded using a depletion method, in which IR absorption was registered by a subsequent decline in ion signal produced by electron ionization of the helium droplets (see also below). This approach confers an element of mass selectivity and can be used to selectively remove contributions of neutral clusters to an IR spectrum smaller than a given size. A strong band was identified which is in good agreement with the absorption band known for the cyclic trimer, (CH3OH)3, in the gas phase.2 However, the bands assigned to the tetramer and pentamer were spread over a wider range than predicted in density functional theory (DFT) calculations for the cyclic isomers. It was concluded that the methanol trimer adopts a cyclic structure in superfluid helium, just like the gas phase, whereas the tetramer and pentamer have structures based upon a cyclic trimer core with one external methanol molecule (tetramer) or a dangling chain of two methanol molecules (pentamer). The possibility of different structures being found in gas phase and helium nanodroplet experiments should come as no surprise. The rapid cooling in helium droplets makes it possible to trap clusters in potential energy wells lying above the global potential energy minimum. A classic example of this is the case of water clusters in helium nanodroplets. For water the cyclic trimer, tetramer and pentamer represent the lowest energy structures.31 However, the situation changes for the hexamer. Here more three-dimensional structures, the prism and cage structures, have a lower energy than the cyclic isomer.31-33 Nevertheless, Nauta and Miller were able to show that the cyclic version of the water hexamer is the primary product in helium nanodroplets.34 The work by Behrens et al. seems to suggest an almost inverse situation for methanol clusters where, instead of forming fully cyclic isomers, branched structures are adopted in liquid helium which are metastable.30 Since liquid helium nanodroplets provide a useful environment in which to study hydrogenbonded clusters, and can sometimes allow access to parts of the intermolecular potential

4 ACS Paragon Plus Environment

Page 5 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

energy surface that are not sampled in gas phase experiments, it is important to understand how this environment can affect the structures of molecular clusters. Although IR spectra of methanol clusters in liquid helium nanodroplets have been reported in the CO stretching region,30 the only previous study in the OH stretching region focused on the methanol monomer.35 Since the OH stretching region is highly sensitive to the nature of the hydrogen bonding, and therefore the structures adopted, we present here IR spectra of methanol clusters in this spectral region. Contrary to the findings by Behrens et al., our data are entirely consistent with gas phase studies and show cyclic structures for the tetramer and pentamer. Furthermore, we present preliminary IR data for clusters of methanol with a single added water molecule, and in particular make the first observation of IR bands of (CH3OH)2H2O and (CH3OH)3H2O. These mixed clusters are also found to adopt cyclic structures.

Experimental Full details of the apparatus can be found in previous publications.36,37 In brief, helium nanodroplets with a mean size close to 5000 helium atoms were formed by expanding precooled gaseous helium into a vacuum through a 5 µm pinhole. The droplets then passed through a pick-up cell, where dopants were added. Experiments were carried out both on pure methanol clusters and for mixed methanol/water clusters. For the latter, water and methanol were pre-mixed in the liquid state and the resulting vapor was then added to the pick-up cell. The flow of vapor into the pick-up cell was controlled by a needle valve. Infrared spectra were recorded via a depletion method which exploits the evaporative loss of helium atoms as the absorbed energy becomes dispersed into the droplet.27-29 This effect causes a decrease in the droplet size and therefore its geometric cross section, which can be registered by a drop in the signal identified via electron ionization mass spectrometry

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

downstream of the IR excitation zone. In the current experiments mass-selective ion detection was used in order to provide a degree of mass selectivity: this aids spectral assignment, as mentioned earlier and also detailed later.

Computational details Ab initio calculations were performed in support of the experiments. The calculations used MP2 methodology combined with aug-cc-pVDZ basis sets centred on the various atoms. The only exception was for the (3)+2 structure of the methanol pentamer (see later), which would only converge to a stable structure using the DFT/B3LYP method (also with an aug-cc-pVDZ basis set). The principal aim here was to predict the IR spectra for small methanol clusters and mixed methanol/water clusters in order to confirm or eliminate specific IR band assignments. Plausible equilibrium structures were first calculated using geometry optimization and then IR absorption spectra were calculated using the double harmonic approximation, which means that the vibrational frequencies need to be scaled for any meaningful comparison with experiment. Since our focus later is on the stretching bands from OH groups, and especially those involved in hydrogen bonds (bonded OH stretches), we have used a comparison of the well-known bonded OH band of methanol dimer with the calculated equivalent to determine a suitable scaling factor (0.972). This scaling factor was then applied to the calculated harmonic vibrational frequencies of all other clusters, including the mixed methanol-water clusters. The corresponding scaling factor for the DFT calculation on the (3)+2 structure of the methanol pentamer was 0.976. All calculations were performed using Gaussian 0938 on the Slater cluster at the UK National Service for Computational Chemistry Software (NSCCS). Images of the equilibrium structures of all of the clusters calculated in this work are provided in the Supporting Information.

6 ACS Paragon Plus Environment

Page 7 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Results and Discussion

1.

Methanol clusters

Figure 1 shows several IR spectra recorded for pure methanol clusters in the OH stretching region. These spectra differ in the mass channel used for detection, the size of the helium droplets used, and the quantity of added methanol. The focus in this paper is on the OH stretching region above 3000 cm-1, although the spectra in Figure 1(a) and 1(b) extend to lower wavenumbers and so also show CH stretching structure. The lowest mass channel used for detection was m/z 33 and corresponds to detection of (CH3OH)H+ ions. The IR spectrum recorded in this mass channel is shown in trace (a) and consists of three strong OH stretching bands at 3295, 3477 and 3571 cm-1, along with a much weaker peak at 3686 cm-1. The assignment of these bands is straightforward based on a comparison with previous IR work and is summarized in Figure 1 and in Table 1. The two highest frequency and sharpest peaks derive from methanol dimer. This cluster is held together by a single hydrogen bond and so there is a dangling OH, which is the source of the weak band at 3686 cm-1. The other OH group is involved in the hydrogen bond, which weakens the OH bond and induces a significant red-shift of the corresponding IR band. Hydrogen bond formation also strongly increases the transition dipole moment. The two other, broader, peaks seen in Figure 1 near 3480 and 3300 cm-1 match the expected positions for the cyclic trimer and tetramer, respectively, of methanol (see Table 1). In these cyclic structures all of the OH groups are involved in hydrogen bonds and so there is no free OH stretch at higher frequency.

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 28

Fig. 1. Infrared spectra of methanol clusters recorded by detection of (a) (CH3OH)H+ at m/z 33, (b) (CH3OH)3H+ at m/z 97 and (c) (CH3OH)3H+ at m/z 97 but at a higher methanol pick-up pressure than in (b).

Behrens et al. have previously reported that the fragmentation of methanol cluster ions in helium nanodroplets following electron ionization is dominated by the following processes,30

(CH3OH)n+



(CH3OH)n-1H+ + CH3O

(1a)



(CH3OH)n-2H+ + CH3O + CH3OH

(1b)

8 ACS Paragon Plus Environment

Page 9 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

According to this scheme, if (CH3OH)H+ is detected any depletion induced by IR absorption should derive from either neutral (CH3OH)2 or (CH3OH)3. In Figure 1(a) we see bands from both of these clusters, but also a substantial contribution from (CH3OH)4. There is also some weak absorption at the expected position for the pentamer (see below), and so the fragmentation process is more extensive than implied in the scheme above. Nevertheless, a marked change is seen when detection switches to (CH3OH)3H+ at m/z 97. In this IR spectrum, shown in blue in Figure 1(b), the IR bands attributed to methanol dimer and trimer are now absent, as would be expected since the smallest neutral cluster that can contribute to the ion signal at m/z 97 is the tetramer. At this higher mass the tetramer absorption is clearly seen along with another equally intense band at a lower frequency, which can be assigned to the pentamer. At higher methanol pick-up pressures, as shown in Figure 1(c), the peak assigned to the pentamer is significantly broadened. This suggests underlying and unresolved contributions from larger clusters, such as the hexamer.

Positions (cm-1) and assignments of OH stretching bands of methanol clusters.

Table 1.

This work

CRDS7

FTIR8

3686

3683.8

3571

3574.4

3575

(CH3OH)2

3477

3473.2

3474

(CH3OH)3

3295

3278.5a)

3293

(CH3OH)4

(CH3OH)2

3239 a)

Assigned carrier

(CH3OH)5

Thought to be a typographical error in the original reference.7

There can be no doubt that we are seeing the cyclic versions of the methanol trimer, tetramer and pentamer. The band positions are in excellent agreement with previous work in the gas phase, where the assignment to cyclic clusters has been firmly established.7,8

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 28

However, in their helium droplet experiment Behrens et al. interpreted their IR data in the CO stretching region rather differently.30 For the trimer the conclusion was the same, namely a cyclic structure is formed. However, absorption features assigned to the tetramer were attributed to a (3)+1 structure, which consists of a cyclic trimer core with one additional methanol molecule hydrogen bonded to the ring. Likewise, the pentamer signal was assigned to a (3)+2 structure, where a chain of two methanol molecules is attached to the trimer ring. We can eliminate these isomers under our experimental conditions by using ab initio calculations, whose findings are shown in Figure 2. The predicted IR spectra for the fully cyclic clusters are simple and are in good agreement with the observed IR spectra. On the other hand, the (3)+1 and (3)+2 branched structures are predicted to have intense bands at positions where no absorption is seen in the IR spectra. The predicted spectra of the branched structures are more complicated than for the fully cyclic structures, as expected since there are several quite distinct OH environments which will deliver substantially different stretching frequencies. This is not consistent with our spectroscopic observations and so we conclude that homodromic cyclic methanol clusters dominate in our experiments.

10 ACS Paragon Plus Environment

Page 11 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Fig. 2. Comparison of the infrared spectrum of methanol clusters in the OH stretching region recorded by detection of (CH3OH)H+ at m/z 33 (upper trace) with ab initio predictions (lower trace) for the fully cyclic trimer and tetramer along with the (3)+1 tetramer and (3)+2 pentamer. The predicted band positions for fully cyclic trimer and tetramer are in excellent agreement with experimental band positions, unlike the predictions for the (3)+1 and (3)+2 isomers.

The reason for the different conclusions reached in this study versus that by Behrens et al.30 is unclear. We broadly agree with their theoretical predictions of the IR absorption behaviour in the CO stretching region and their conclusions about the cluster structures on that basis alone do not seem unreasonable. We note that Behrens et al.30 claimed to be using helium droplets with a mean size of 2700 helium atoms, which is somewhat smaller than the mean droplet size used here (roughly 5000 helium atoms). Nevertheless, we would be surprised if this led to different cluster structures. Another possibility is the neglect of 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 28

anharmonicity in the calculated CO stretching structure by Behrens et al. Their harmonic oscillator calculations predict a number of CO stretching bands within a relatively narrow 10 cm-1 range, and our own calculations confirm this. However, a slightly wider spread of these features over 20 cm-1 range would be sufficient to account for the band structure seen by Behrens et al. and anharmonic coupling might deliver this. This is clearly a matter worthy of further investigation.

2.

Methanol + water

We have also performed some limited experiments with a mixture of methanol and water. The aim here, rather than being an exhaustive investigation of mixed methanol/water clusters, was to add a small quantity of water to form methanol clusters with one water molecule (primarily) and to see if the resulting mixed clusters remain cyclic. To do, this, an 80:20 (by volume) liquid methanol/water mixture maintained at 2 °C was used for doping the pick-up cell. Given that methanol has a much higher vapor pressure than water, this gave a methanol/water ratio of ∼ 30:1 in the vapor and therefore ensured that the pick-up of two or more water molecules had negligible probability. Figure 3 shows IR depletion spectra recorded at two different masses, m/z 33 and m/z 51. The spectrum at m/z 33 contains contributions from several different pure methanol clusters, as illustrated in the figure, with the dimer and trimer being particularly prominent. However, it also contains two additional peaks marked by arrows, a stronger one at 3525 cm-1 and a weaker one at 3380 cm-1. Both of these peaks become more intense when recorded at m/z 51. Since the ion detected at this mass is (CH3OH)(H2O)H+, pure methanol clusters are eliminated from the m/z 51 spectrum, as seen by the disappearance of the methanol dimer peak. The remaining signals must come from one or more neutral clusters containing at least one water molecule. Given the deliberately low partial pressure of added water, we expect

12 ACS Paragon Plus Environment

Page 13 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that the dominant mixed clusters will contain only one water molecule. Note also that the signal/noise ratio is compromised for the IR spectrum recorded at m/z 51 because only a small subset of the doped helium droplets will contain a single water molecule.

Fig. 3. Infrared spectra obtained for a methanol/water vapour mixture (∼30:1) with detection at m/z 33 ((CH3OH)H+) and m/z 51 ((CH3OH)(H2O)H+). The arrows in the m/z 33 spectrum refer to peaks discussed in the main text. In the lowest trace theoretical predictions of the IR band structure for cyclic (CH3OH)2(H2O) (red lines) and (CH3OH)3(H2O) (blue lines) are shown.

We can assign the absorption bands in the m/z 51 spectrum in Figure 3 with the aid of ab initio calculations and we find that all of the observed features can be explained rather well assuming cyclic structures. The two most intense bands, at ca. 3300 and 3480 cm-1, coincide with bands seen in the pure methanol spectrum for the cyclic tetramer and trimer, respectively. According to the ab initio predictions we also expect strong bands at roughly 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

these positions for cyclic (CH3OH)3H2O and (CH3OH)2H2O clusters. However, theory also predicts additional bands near 3400 cm-1 (for both (CH3OH)2H2O and (CH3OH)3H2O) and 3530 cm-1 ((CH3OH)2H2O only). These predictions fit with the observation of new bands in the IR spectrum at 3380 and 3525 cm-1. Furthermore, the relative band intensities predicted by the ab initio calculations are in good agreement with the experimental relative band intensities. We therefore conclude that we are observing IR bands of the cyclic mixed methanol/water clusters (CH3OH)2H2O and (CH3OH)3H2O. Although mixed methanol/water clusters have been studied previously, most of the focus has been on the heterodimer, (CH3OH)H2O,39-47 including its observation in helium nanodroplets during a study of the OH(CH3OH) radical-molecule complex.48 In contrast, experimental and theoretical studies of larger (CH3OH)nH2O clusters are relatively few.43,44,49-53 The only spectroscopic work was a combined IR and Raman study by Nedić et al. in which a single Raman band was assigned to (CH3OH)2H2O in the OH stretching region and tentative assignments were made of Raman bands to (CH3OH)3H2O and (CH3OH)4H2O.51 However, no IR bands were assigned to these clusters and so the present study represents, as far as we are aware, the first spectroscopic observation of IR bands of (CH3OH)2H2O and (CH3OH)3H2O. The Raman band reported by Nedić et al. for (CH3OH)2H2O was observed at 3425 cm-1. This band almost certainly corresponds to the lowest frequency OH stretching band in the simulation for (CH3OH)2H2O in Figure 3 (the lowest frequency line in red), which the ab initio calculations predict to be at 3410 cm-1. The vibration responsible is essentially a symmetrical ring breathing vibration and is the weakest of the four IR OH stretching bands but is the most intense (by a factor of more than three) of the Raman active bands of (CH3OH)2H2O, according to the ab initio calculations. This presumably explains why Nedić et al. were able to see the 3425 cm-1 band in their Raman spectrum but we see nothing above the noise level at this position in our IR spectrum. For the

14 ACS Paragon Plus Environment

Page 15 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

other bands there is good agreement between theory and experiment, which shows that in liquid helium droplets the clusters resulting from a homogeneous mixture of methanol and water vapor derive from self-assembly into fully cyclic clusters. Calculations suggest that these

cyclic

structures

are

the

lowest

energy

isomers

of

(CH3OH)2H2O

and

(CH3OH)3H2O.43,44,49 Finally, we note that there has been much interest over many decades in trying to understand why liquid methanol-water mixtures show a negative entropy of mixing. The prevailing view for several decades was the so-called “iceberg model”, in which water molecules were assumed to cluster together on the hydrophobic methyl groups.54 However, more recent experimental and theoretical work tends to favour an explanation in terms of incomplete mixing of the water and the methanol,55-59 although the debate continues right up to the present.60 The study in this paper cannot shed any light on this issue, except to say that there is no evidence of segregation of the methanol and water in the very limited experiments performed here, i.e. no structure based on a water molecule dangling from a cyclic methanol ring. It appears that the water molecule simply slots into the cyclic structure adopted by pure methanol clusters. However, it would be interesting to determine if segregation occurs in mixed methanol/water clusters with more equal quantities of both solvents, such as (CH3OH)3(H2O)3. Several possible structures are conceivable for this particular mixed cluster and experiments on such clusters are viable with the approach described in this paper.

Conclusions Infrared spectra in the OH stretching region have been recorded for pure methanol clusters trapped in liquid helium nanodroplets. The data are in excellent agreement with previous gas phase studies and show that the trimer, tetramer and pentamer adopt fully cyclic structures in helium nanodroplets. This conclusion is in marked contrast to an earlier study of these

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 28

clusters in helium nanodroplets, which focused on the CO stretching region and came to the conclusion that while the trimer is cyclic, the tetramer and pentamer adopt structures based upon a cyclic trimer, in which additional methanol molecules dangle from this ring. We have also reported the first IR spectra of (CH3OH)2H2O and (CH3OH)3H2O. In this preliminary work there is good agreement between the observed spectra and the predicted spectra from ab initio calculations when fully cyclic structures are assumed. This work could be readily extended to other combinations of methanol and water and therefore to more fully explore the structural landscape and hydrogen bonding in these mixed solvent clusters.

Supporting information Images of the equilibrium structures of the various clusters discussed in this paper, which were derived from ab initio calculations, are contained within the Supporting Information along with the Cartesian coordinates of the atoms.

Acknowledgements The authors wish to thank the UK Engineering and Physical Sciences Research Council (EPSRC, grant EP/J021342/1) and the Leverhulme Trust (grant numbers RPG-2012-552 and RPG-2012-740) for grants in support of this work. The authors are also grateful for the scholarship provided to MIS by the Human Capacity Development Program from the Kurdistan Government.

References

16 ACS Paragon Plus Environment

Page 17 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1.

Odutola, J. A.; Viswanathan, R.; Dyke, T. R. Molecular beam electric deflection behavior and polarity of hydrogen-bonded complexes of ROH, RSH, and RNH. J. Am. Chem. Soc. 1979, 101, 4787–4792.

2.

Buck, U.; Gu, X. J.; Hobein, M.; Lauenstein, C.; Rudolph, A. Infrared photodissociation spectra of size-selected (CH3OH)n clusters from n = 2 to n = 8. J. Phys. Chem. 1988, 92, 5561- 5562.

3.

Huisken, F.; Kulcke, A.; Lausch, C.; Lisy, J. M. Dissociation of small methanol clusters after excitation of the O–H stretch vibration at 2.7 µ. J. Chem. Phys. 1991, 95, 3924-3929.

4.

Lovas, F. J.; Belov, S. P.; Tretyakov, M. Y.; Stahl, W.; Suenram, R. D. The a-type K = 0 microwave spectrum of the methanol dimer. J. Molec. Spectrosc. 1995, 170, 478492.

5.

Huisken, F.; Kaloudis, M.; Koch, M.; Werhahn, O. Experimental study of the O-H ring vibrations of the methanol trimer. J. Chem. Phys. 1996, 105, 8965-8968.

6.

Buck, U.; Ettischer, I. Vibrational predissociation spectra of size selected methanol clusters: new experimental results. J. Chem. Phys. 1998, 108, 23-38.

7.

Provencal, R. A.; Paul, J. B.; Roth, K.; Chapo, C.; Casaes, R. N.; Saykally, R. J.; Tschumper, G. S.; Schaefer III, H. F. Infrared cavity ringdown spectroscopy of methanol clusters: single donor hydrogen bonding. J. Chem. Phys. 1999, 110, 42584267.

8.

Häber, T.; Schmitt, U.; Suhm, M. A. FTIR-spectroscopy of molecular clusters in pulsed supersonic slit-jet expansions. Phys. Chem. Chem. Phys. 1999, 1, 5573-5582.

9.

Buck, U.; Huisken, F. Infrared spectroscopy of size-selected water and methanol clusters. Chem. Rev. 2000, 100, 3863-3890.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

10.

Page 18 of 28

Steinbach, C.; Fárnik, M.; Ettischer, I.; Siebers, J.; Buck, U. Isomeric transitions in size-selected methanol hexamers probed by OH-stretch spectroscopy. Phys. Chem. Chem. Phys. 2006, 8, 2752-2758.

11.

Fu, H. B.; Hu, Y. J.; Bernstein, E. R. IR + vacuum ultraviolet (118 nm) nonresonant ionization spectroscopy of methanol monomers and clusters: neutral cluster distribution and size-specific detection of the OH stretch vibrations. J. Chem. Phys.

2006, 124, 024302. 12.

Hu, Y. J.; Fu, H. B.; Bernstein, E. R. Infrared plus vacuum ultraviolet spectroscopy of neutral and ionic methanol monomers and clusters: New experimental results. J. Chem. Phys. 2006, 125, 154306.

13.

Larsen, R. W.; Zielke, P.; Suhm, M. A. Hydrogen-bonded OH stretching modes of methanol clusters: a combined IR and Raman isotopomer study. J. Chem. Phys. 2007, 126, 194307.

14.

Han, H. –L.; Camacho, C.; Witek, H. A.; Lee, Y. –P. Infrared absorption of methanol clusters (CH3OH)n with n = 2-6 recorded with a time-of-flight mass spectrometer using infrared depletion and vacuum-ultraviolet ionization. J. Chem. Phys. 2011, 134, 144309.

15.

Coussan, S.; Loutellier, A.; Perchard, J. P.; Racine, S.; Peremans, A.; Tadjeddine, A.; Zheng, W. Q. Infrared laser induced isomerization of methanol polymers trapped in nitrogen matrix. I. Trimers. J. Chem. Phys. 1997, 107, 6526-6540.

16.

Coussan, S.; Bouteiller, Y.; Loutellier, A.; Perchard, J. P.; Racine, S.; Peremans, A.; Zheng, W. Q.; Tadjeddine, A. Infrared photoisomerization of the methanol dimer trapped in argon matrix; monochromatic irradiation experiments and DFT calculations. Chem. Phys. 1997, 219, 221-234.

18 ACS Paragon Plus Environment

Page 19 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

17.

Coussan, S.; Loutellier, A.; Perchard, J. P.; Racine, S.; Peremans, A.; Tadjeddine, A.; Zheng, W. Q. IR-induced interconversions between five conformers of methanol dimer trapped in nitrogen matrix. Chem. Phys. 1997, 223, 279-292.

18.

Perchard, J. P.; Mielke, Z. Anharmonicity and hydrogen bonding I. A near-infrared study of methanol trapped in nitrogen and argon matrices, Chem. Phys. 2001, 264, 221-234.

19.

Doroshenko, I. Y. Matrix isolation study of the formation of methanol cluster structures in the spectral region of C–O and O–H stretch vibrations. Low Temp. Phys.

2011, 37, 604-608. 20.

Hagemeister, F. C.; Gruenloh, C. J.; Zwier, T. S. Density functional theory Calculations of the structures, binding energies, and infrared spectra of methanol clusters. J. Phys. Chem. A 1998, 102, 82-94.

21.

David, J.; Guerra, D.; Restrepo, A. Structural characterization of the (methanol)4 potential energy surface. J. Phys. Chem. A 2009, 113, 10167-10173.

22.

Do, H.; Besley, N. A. Structural optimization of molecular clusters with density functional theory combined with basin hopping. J. Chem. Phys. 2012, 137, 134106.

23.

Kazachenko, S.; Bulusu, S.; Thakkar, A. J. Methanol clusters (CH3OH)n: putative global minimum-energy structures from model potentials and dispersion-corrected density functional theory. J. Chem. Phys. 2013, 138, 224303.

24.

Nishimura, Y.; Lee, Y. –P.; Irle, S.; Witek, H. A. Critical interpretation of CH– and OH– stretching regions for infrared spectra of methanol clusters (CH3OH)n (n = 2–5) using self-consistent-charge density functional tight-binding molecular dynamics simulations. J. Chem. Phys. 2014, 141, 094303.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

25.

Page 20 of 28

Heger, M.; Suhm, M. A.; Mata, R. A. Towards the binding energy and vibrational red shift of the simplest organic hydrogen bond: Harmonic constraints for methanol dimer. J. Chem. Phys. 2014, 141, 101105.

26.

Albrecht, L.; Boyd, R. J. Atomic energy analysis of cooperativity, anti-cooperativity, and non-cooperativity in small clusters of methanol, water, and formaldehyde. Comp. Theor. Chem. 2015, 1053, 328-336.

27.

Toennies, J. P.; Vilesov, A. F. Angew. Chem. Int. Ed. 2004, 43, 2622-2648.

28.

Yang, S.; Ellis, A. M. Helium droplets: a chemistry perspective. Chem. Soc. Rev.

2013, 42, 472-484. 29.

Callegari, C.; Ernst, W. E. Helium Droplets as Nanocryostats for Molecular Spectroscopy—from the Vacuum Ultraviolet to the Microwave Regime. Handbook of High-resolution Spectroscopy, Wiley, 2011.

30.

Behrens, M.; Fröchtenicht, R.; Hartmann, M.; Siebers, J. –G.; Buck, U.; Hagemeister, F. C. Vibrational spectroscopy of methanol and acetonitrile clusters in cold helium droplets. J. Chem. Phys. 1999, 111, 2436-2443.

31.

Bates, D. M.; Tschumper, G. S. CCSD(T) Complete basis set limit relative energies for low-lying water hexamer structures. J. Phys. Chem. A 2009, 113, 3555-3559.

32.

Temelso, B.; Archer, K. A.; Shields, G. C. Benchmark structures and binding energies of small water clusters with anharmonicity corrections, J. Phys. Chem. A 2011, 115, 12034-12046.

33.

Pérez, C.; Muckle, M. T.; Zaleski, D. P.; Seifert, N. A.; Temelso, B.; Shields, G. C.; Kisiel, Z.; Pate, B. H. Structures of cage, prism, and book isomers of water hexamer from broadband rotational spectroscopy. Science 2012, 336, 897-901.

34.

Nauta, K.; Miller, R. E. Formation of cyclic water hexamer in liquid helium: The smallest piece of ice. Science 2000, 287, 293-295.

20 ACS Paragon Plus Environment

Page 21 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

35.

Raston. P. L.; Douberly, G. E.; Jäger, W. Single and double resonance spectroscopy of methanol embedded in superfluid helium nanodroplets. J. Chem. Phys. 2014, 141, 044301.

36.

Shepperson, B.; Tandy, J.; Boatwright, A.; Feng, C.; Spence, D.; Shirley, A.; Yang, S.; Ellis, A. M. Electronic spectroscopy of toluene in helium nanodroplets: evidence for a long-lived excited state. J. Phys. Chem. A 2013, 117, 13591-13595.

37.

Tandy, J.; Feng, C.; Boatwright, A.; Sarma, G.; Sadoon, A. M.; Shirley, A.; Das Neves Rodrigues, N.; Cunningham, E. M.; Yang, S.; Ellis, A. M. Infrared spectroscopy of salt-water complexes. J. Chem. Phys. 2016, 144, 121103.

38.

M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, et al., Gaussian 09, Revision E.01, Gaussian, Inc., Wallingford CT, 2009.

39.

Huisken, F.; Stemmler, M. On the structure of the methanol-water dimer. Chem. Phys. Lett. 1991, 180, 332-338.

40.

Bakkas, N.; Bouteiller, Y.; Loutellier, A.; Perchard, J. P.; Racine, S. The watermethanol complexes. I. A matrix isolation study and an ab initio calculation on the 11 species. J. Chem. Phys. 1993, 99, 3335-3342.

41.

Bakkas, N.; Bouteiller, Y.; Loutellier, A.; Perchard, J. P.; Racine, S. The watermethanol complexes. Matrix induced structural conversion of the 1-1 species. Chem. Phys. Lett. 1995, 232, 90-98.

42.

Stockman, P. A.; Blake, G. A.; Lovas, F. J.; Suenram, R. D. Microwave rotationtunneling spectroscopy of the water-methanol dimer: Direct structural proof for the strongest bound conformation. J. Chem. Phys. 1997, 107, 3782-3790.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

43.

Page 22 of 28

González, L.; Mó, O.; Yáñez, Y. High level ab initio and density functional theory studies on methanol–water dimers and cyclic methanol(water)2 trimer. J. Chem. Phys.

1998, 109, 139-150. 44.

Iosue, J. L.; Benoit, D. M.; Clary, D. C. Diffusion Monte Carlo simulations of methanol-water clusters. Chem. Phys. Lett. 1999, 301, 275-280.

45.

Moskowitz, J. W.; Bačić, Z.; Sarsa, A.; Schmidt, K. E. Relative stabilities of the two isomers of the methanol-water dimer: The effects of the internal rotations of the hydroxyl and methyl groups of methanol. J. Chem. Phys. 2001, 114, 10294-10299.

46.

Yaramoto, N.; Nishino, Y.; Miyoshi, E. Classical trajectory calculations of intramolecular vibrational energy redistribution. I. Methanol-water complex.

J.

Chem. Phys. 2004, 121, 2058-2066. 47.

Fileti, E. E.; Canuto, S. Calculated infrared spectra of hydrogen-bonded methanolwater, water-methanol, and methanol-methanol complexes. Int. J. Quantum Chem.

2005, 104, 808-815. 48.

Hernandez, F. J.; Brice, J. T.; Leavitt, C. M.; Pino, G. A.; Douberly, G. E. Infrared spectroscopy of OH··CH3OH: Hydrogen-bonded intermediate along the hydrogen abstraction reaction path J. Phys. Chem. A 2015, 119, 8125-8132.

49.

Mejía, S. M.; Espinal, J. F.; Mondragón, F. Cooperative effects on the structure and stability of (ethanol)3-water, (methanol)3-water heterotetramers and (ethanol)4, (methanol)4 tetramers. J. Mol. Struct: THEOCHEM 2009, 901, 186-193.

50.

Mandal, A.; Prakash, M.; Kumar, R. M.; Parthasarathi, R.; Subramanian, V. Ab initio and DFT studies on methanol-water clusters. J. Phys. Chem. A 2010, 114, 2250-2258.

51.

Nedić, M.; Wassermann, T.; Larsen, R. W.; Suhm, M. A. A combined Raman and infrared jet study of mixed methanol-water and ethanol-water clusters. Phys. Chem. Chem. Phys. 2011, 13, 14050-14063.

22 ACS Paragon Plus Environment

Page 23 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

52.

Mejía, S. M.; Florez, E.; Mondragon, F. An orbital and electron density analysis of weak interactions in ethanol-water, methanol-water, ethanol and methanol small clusters. J. Chem. Phys. 2012, 136, 144306.

53.

Mejía, S. M.; Espinal, J. F.; Mills, M. J. L.; Mondragón, F. The role of OH…O and CH…O hydrogen bonds and H…H interactions in ethanol/methanol–water heterohexamers. J. Mol. Model. 2016, 22, 181-190.

54.

Frank, H. S.; Evans, M. W. Free volume and entropy in condensed systems. III. Entropy in binary liquid mixtures; partial molar entropy in dilute solutions; structure and thermodynamics in aqueous electrolytes. J. Chem. Phys. 1945, 13, 507-532.

55.

Dixit, S.; Crain, J.; Poon, W. C. K.; Finney, J. L.; Soper, A. K. Molecular segregation observed in a concentrated alcohol-water solution. Nature 2002, 416, 829-832.

56.

Guo, J. –H.; Luo, Y.; Augustsson, A.; Kashtanov, S.; Rubensson, J. -E.; Shuh, D. K.; Ågren, H.; Nordgren, J. Molecular structure of alcohol-water mixtures. Phys. Rev. Lett. 2003, 91, 157401.

57.

Pascal, T. A.; Goddard III, W. A. Hydrophobic segregation, phase transitions and the anomalous thermodynamics of water/methanol mixtures. J. Phys. Chem. B 2012, 116, 13905-13912.

58.

Benson, S. P.; Pleiss, J. Incomplete mixing versus clathrate-like structures: a molecular view on hydrophobicity in methanol-water mixtures, J. Mol. Model. 2013, 19, 3427-3436.

59.

Matisz, G.; Kelterer, A. –M.; Fabian, W. M. F.; Kunsági-Máté, S. Structural properties of methanol–water binary mixtures within the quantum cluster equilibrium model. Phys. Chem. Chem. Phys. 2015, 17, 8467-8479.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

60.

Page 24 of 28

Lam, R. K.; Smith, J. W.; Saykally, R. J. Hydrogen bonding interactions in wateralcohol mixtures from X-ray absorption spectroscopy. J. Chem. Phys. 2016, 144, 191103.

24 ACS Paragon Plus Environment

Page 25 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ToC Graphic

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1 272x208mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2 289x202mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3 289x202mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 28