Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)
Article
An Electron Momentum Spectroscopy Investigation on Molecular Conformations of Ethanol Considering Vibrational Effects YaGuo Tang, Xu Shan, Shanshan Niu, Zhaohui Liu, Enliang Wang, Noboru Watanabe, Masakazu Yamazaki, Masahiko Takahashi, and Xiangjun Chen J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.6b10009 • Publication Date (Web): 20 Dec 2016 Downloaded from http://pubs.acs.org on December 23, 2016
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 44
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
An Electron Momentum Spectroscopy Investigation on Molecular Conformations of Ethanol Considering Vibrational Effects Yaguo Tang,† Xu Shan,† Shanshan Niu,† Zhaohui Liu,† Enliang Wang,† Noboru Watanabe,‡ Masakazu Yamazaki,‡ Masahiko Takahashi,‡ and Xiangjun Chen*,†
†
Hefei National Laboratory for Physical Sciences at the Microscale and Department
of Modern Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China ‡
Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, Sendai 980-8577, Japan
ABSTRACT The interpretation of experimental electron momentum distributions (EMDs) of ethanol, one of the simplest molecules having conformers, confuses researchers for years. High-level calculations of Dyson orbital EMDs by thermally averaging the gauche and trans conformers as well as molecular dynamical simulations failed to quantitatively reproduce the experiments for some of the outer valence orbitals. In this work, the valence shell electron binding energy spectrum and EMDs of ethanol are revisited by the high-sensitivity electron momentum spectrometer employing symmetric non-coplanar geometry at an incident energy of 1200 eV plus binding energy, together with a detailed analysis of the influence of vibrational motions on the EMDs for the two conformers employing a harmonic analytical quantum mechanical
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
(HAQM) approach by taking into account all the vibrational modes. The significant discrepancies between theories and experiments in the previous works have now been interpreted quantitatively, indicating that the vibrational effect plays a significant role in reproducing the experimental results, not only through the low-frequency OH and CH3 torsion modes but also through other high-frequency ones. The rational explanation of experimental momentum profiles provides solid evidence that the trans conformer is slightly more stable than the gauche conformer, in accordance with the thermodynamical predictions and other experiments. The case of ethanol demonstrates the significance of considering vibrational effects when performing conformational study on flexible molecules using electron momentum spectroscopy.
1. INTRODUCTION Conformational phenomena are ubiquitous for polyatomic molecules in nature. The conformational study, including the confirmation of the existence of certain conformers and the determination of the relative stability, has attracted much interest for researchers due to the rather discrepant chemical reactivity of different conformers. Based on the fact that conformers with different geometries generally exhibit quite different electronic wavefunctions for corresponding orbitals, electron momentum spectroscopy (EMS),1-5 a robust technique for exploring the electronic structures of atoms and molecules, has been applied to this field successfully6-24 for years owing to its unique ability of imaging electron momentum distributions (EMDs) for individual molecular orbitals (MOs). The EMS technique is based on the high energy electron
2 ACS Paragon Plus Environment
Page 2 of 44
Page 3 of 44
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
impact (e, 2e) reaction near the Bethe ridge. Within a series of approximations including Born-Oppenheimer, binary encounter, weak-coupling and plane wave impulse approximations, the triple differential cross section (TDCS) of the (e, 2e) reaction is proportional to the modulus square of the wavefunction of the ionized orbital in momentum space, i.e., electron momentum profile. Many meaningful EMS experiments have been carried out to investigate molecular conformations,6-24 and even to determine the relative abundances of conformers with the aid of theoretical simulations.16 A vital point when employing EMS on conformational study is the precise prediction of momentum profiles for each conformer in theory because conclusions are generally made through comparison between experimental and theoretical ones. Previous studies have proved that several significant factors, including electron correlation,25-27 distorted wave effects,28,29 and nuclear dynamics in the initial and final states,30-42 have to be taken into account when interpreting EMS results due to the possible invalidity of the approximations used in EMS and the complexity of polyatomic molecules. Therefore, careful treatment should be carried out to produce theoretical momentum profiles with high quality in conformational investigation through EMS. Ethanol, one of the most common and important organic molecules, has two conformers, gauche conformer of C1 symmetry and trans conformer of Cs symmetry, corresponding to the two minima of the potential energy surface about the C-C-O-H internal torsional angle.43,44 Confusion still exists for the understanding of the
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
experimental momentum profiles of outer valence orbitals of ethanol in the previous EMS works.14,15,19 In the EMS study of Ning et al.,15 momentum profiles of the five bands in the outer-valence region were extracted. The theoretical calculation with 60.7/39.3 gauche/trans abundance ratio reproduced the experimental EMDs well for the second, fourth and fifth bands but failed to reproduce those for the first and third bands to a large extent. The first point one may naturally consider is whether the relative abundances adopted are correct since the EMDs for the first and third bands seem to be well interpreted using abundance ratio of 80/20.18,19 However, various experiments including microwave spectroscopy45,46, far infrared spectroscopy47 and X-ray photoelectron spectroscopy48 as well as theoretical calculations with high accuracy14 exhibit that the trans conformer is slightly more stable than the gauche one and that the relative abundances of the gauche and trans conformers are nearly 61% and 39%, respectively, considering the inherent twofold degeneracy of the gauche conformation. The 80/20 ratio is not appropriate in that it destroys the agreement for other bands dramatically and not in accordance with the thermodynamic prediction and other experimental results. The discrepancies may also be attributed to the distorted wave effect which sometimes results in turn-up in EMD at low momentum. However, the consistent EMS results at incident energies of 1200 eV and 2400 eV proved that the distorted wave effect has little influence on the momentum profiles of ethanol,15 and can be excluded from the possible reasons for the deviation. The next possible reason for the discrepancies can be the inadequate consideration of electron correlation. High-level calculation using one particle Green’s function theory in
4 ACS Paragon Plus Environment
Page 4 of 44
Page 5 of 44
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
conjunction with the so called third-order Algebraic diagrammatic construction scheme (ADC(3))15 showed that electron correlation can be well described by the B3LYP Kohn-Sham orbitals with suitably large basis sets, and only a very limited improvement for the highest occupied molecular orbital (HOMO) was achieved when using Dyson orbitals from the ADC(3) calculations. Finally, nuclear dynamics in the initial ground and final ionized states may also be responsible for the discrepancies. Nuclear dynamics in the ground state in terms of the OH and CH3 torsions were simulated,19 and gave a feeble improvement for HOMO. Hajgato et al. proposed19 that the structural relaxation induced by the ionization process may be one of the tempting explanations for HOMO. To sum up, by taking all these factors into account, the discrepancies between experiments and theoretical calculations still exist. In this work, we performed a new EMS experiment along with theoretical calculations involving vibrational effects, to check the possible reasons of the disagreement between experimental and theoretical EMDs on both experimental and theoretical aspects. On the experimental side, the inferred momentum profiles from (e, 2e) angular distributions of TDCSs should be checked employing another independent instrument to avoid systematic errors. During the data analysis, the complication of binding energies due to molecular conformation gives rise to the difficulty of extracting accurate momentum profiles. The symmetries of individual orbitals should be checked to make sure whether using one Gaussian peak to represent two or more experimentally unresolved ionization bands in the deconvolution process is suitable. On the theoretical side, nuclear dynamics in the initial ground state, or vibrational
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
effects, considering all degrees of freedom have not been examined until now. Recent investigations30-40 have already demonstrated that vibrational motions do make a noticeable influence and play important roles in explaining the substantial unexpected turn-up in the low momentum region of experimental momentum profiles for specific orbitals of polyatomic molecules. Watanabe et al.’s harmonic analytical quantum mechanical (HAQM) approach method33,34,39 provided a simple and practical way to estimate nuclear dynamics in the ground state including all degrees of freedom. It turns out that vibrational modes with extremely high frequency can also influence EMDs greatly, e.g. the effects of the v5 mode (~3000 cm-1) on 1b3g orbital for ethylene33. Our recent work on oxetane40 indicated that both the low-frequency and high-frequency vibrational modes may have a non-negligible influence on EMDs. In this paper, all vibrational modes are taken into account to perform a comprehensive analysis of vibrational effects on the two conformers of ethanol, not just the OH and CH3 torsion modes.
2. EXPERIMENTAL AND THEORETICAL BACKGROUND EMS is a coincidence experiment based on the kinematically complete (e, 2e) reaction, in which a fast incident electron is scattered by a target atom or molecule and meanwhile an electron is knocked out from the target. The conservation of energy and momentum gives the binding energy εf and the momentum p of the target orbital electron as1-5
ε f = E0 − E1 − E 2 ,
(1)
6 ACS Paragon Plus Environment
Page 6 of 44
Page 7 of 44
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
p = p1 + p2 − p0 ,
(2)
where E0, E1, and E2 are energies, and p0, p1, and p2 are momenta of the incident and two outgoing electrons, respectively. At symmetric non-coplanar geometry, the magnitude of momentum can be expressed as 2
p=
φ 2 ( p0 − 2 p1 cos θ ) + 2 p1 sin θ sin , 2
(3)
where θ is the polar angle and ϕ is the relative azimuthal angle between the two outgoing electrons. The high-sensitivity EMS apparatus used in this experiment has been described in detail elsewhere49 and only a brief introduction is presented here. The typical energy of the incident electron is 1200 eV plus binding energy. The symmetric non-coplanar kinematics is employed, in which the ejected and scattered electrons with equal polar angles ( θ1 = θ2 = θ = 45o ) and energies ( E1 = E2 ) are analyzed by a 90º sector, 2π spherical electrostatic analyzer. The two electrons are then detected in coincidence by a position sensitive detector placed at the exit plane of the analyzer. With arriving positions of the two outgoing electrons the binding energy and magnitude of the target electron momentum can be deduced. Several modifications have been implemented. The double-half wedge-strip anode of 55 mm in diameter was replaced by a delay line anode (HEX 120 from RoentDek Handels GmbH) to fully utilize the potential of the analyzer. Moreover, the passing energy of the analyzer was decreased from 600 eV to 200 eV using retarding lenses to achieve a better energy resolution. The energy and momentum resolution of the present EMS apparatus are determined to be ~1.3 eV (full width at half maximum (FWHM)) and 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 44
~0.2 a.u. (∆θ = 0.8º, ∆ϕ = 2.0º), respectively, by measuring Ar 3p orbital. On the theoretical aspect, within the binary encounter and plane wave impulse approximations for randomly oriented atoms or molecules the TDCS for EMS can be expressed as5
σ EMS =
d 3σ pp pp = (2π )4 1 2 f ee Σav pf i = (2π ) 4 1 2 f ee M ( p) , (4) d Ω1d Ω2 dE1 p0 p0
where fee is the e-e collision factor and Σav represents an average over initial i a sum over finial
f
and
degenerate states. M(p) is the relevant structure factor
corresponding to the momentum profile. According to Born-Oppenheimer approximation, the total wavefunction of a molecular system can be separated as a product of electronic, vibrational, and rotational wavefunctions. Ignoring rotational motions and employing the closure relation of the final vibrational eigenstates considering the limited energy resolution of current EMS apparatuses, we can simplify the relevant structure factor as33,34 M ( p ) = ∑ v pv (T ) ∫ χ iv (Q ) ρ f ( p; Q ) dQ , (5) 2
where pv (T ) is the population of the eigenstate v in the ground state at temperature T and χ iv (Q ) is the vibrational wavefunction. ρ f ( p; Q ) is the spherically averaged EMD at a given geometry Q, which reads
ρ f ( p; Q )=
1 2 S f (Q ) ∫ ϕ f ( p; Q ) d Ω p , (6) 4π
where Sf (Q ) is the pole strength or spectroscopic factor and ϕ f ( p; Q ) is the normalized Dyson orbital. If we ignore the change of ρ f ( p; Q ) due to vibrational displacements, in other words, vibrational effects are neglected, the relevant structure factor is no other than the spherically averaged EMD at equilibrium geometry, 8 ACS Paragon Plus Environment
Page 9 of 44
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
indicating the capability of EMS probing MOs in momentum space. Watanabe et al.’s HAQM method33,34 makes it available to estimate vibrational effects from every vibrational mode on momentum profiles and to understand which kind of and how vibrational motions influence EMDs. In this method, the rigid-rotator harmonic oscillator (RRHO) approximation is adopted, which can be expressed as
χiv (Q) = ∏ L ξv (QL ) ,
(7)
L
where ξ vL (QL ) is the harmonic function describing molecular vibrations. Then vibrational effects on momentum profiles can be written as a sum of the contribution of each vibrational mode if the cross-terms between vibrational modes are neglected:33,34 M ( p) = ρ f ( p; Q0 ) + ∑ L ∆M L ( p) = ρ f ( p; Q0 ) + ∑ L
(∑
2
vL
pvL (T ) ∫ ξ vL ( QL ) ρf ( p; Q0 + QL qˆL )dQL − ρf ( p; Q0 )
) , (8)
where ∆ M L ( p ) is the contribution to vibrational effects from the L-th normal mode. QL is the displacement from the equilibrium geometry Q0, and qˆ L is the unit vector along the normal coordinate. On the basis of this equation and high-level ab initio calculations, vibrational effects can be evaluated involving all vibrational degrees of freedom. Technically, the structural parameters of the two conformers in the ground state are optimized by density functional theory (DFT) along with the B3LYP functional50,51 and Dunning’s aug-cc-pVTZ basis set52,53 using Gaussian 03 program.54 The B3LYP geometries are good enough compared with experimental data as displayed by Morini et al.14 Then the vibrational frequencies (Table 1 and 2) and
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
normal coordinates of each conformer are computed at the B3LYP/aug-cc-pVTZ level of theory, where the computed frequencies are in good agreement with experimental results44 if the scaling factor of 0.968 from ref 55 is employed. In particular, the torsional normal modes from B3LYP/aug-cc-pVTZ appear as the combinations of the two torsional internal coordinates, labeled as v20 and v21 for both conformers. The OH torsion is dominant for the v20 mode, while the v21 mode is mainly the CH3 torsion. The next step is to fit ρ f ( p; Q0 + QL qˆ L ) as a function of QL employing fourth-order polynomial function to compute the contributions of each mode. In order to do this,
ρ f ( p; Q0 + QL qˆ L ) is calculated by B3LYP/aug-cc-pVTZ where the displacement QL 2
is chosen so that ξv (QL ) / ξv (0) L L
2
= 3/4, 1/2, 1/4, 1/8, and 1/16, respectively. Here
the target Kohn-Sham approximation56 is adopted and the pole strengths are assumed to be constant free of vibrational variations for all valence orbitals and to be unity for outer valence orbitals, which is approximately valid in general. In fact, the ADC(3) calculations14 exhibited that Kohn-Sham orbitals are good approximations of Dyson orbitals and the pole strengths of outer-valence orbitals of ethanol are nearly the same (0.90~0.92) for different orbitals and conformers. Finally the theoretical momentum profiles considering vibrational effects for each conformer are calculated according to eq 8 using the Boltzmann distributions of vibrational eigenstates at room temperature (298 K). Actually, the temperature of the sample in the reaction center is hard to be determined since it may cool down through expansion. However, on the one hand, the relative abundance would not change by more than ~2% if the gas was cooled to even 200 K as the two conformers are nearly isoenergetic. On the other hand, the results of
10 ACS Paragon Plus Environment
Page 10 of 44
Page 11 of 44
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
vibrational effects vary quite a little when a lower temperature, like 270 K, is employed. Noting that the temperature will not decrease too much for a normal diffusive gas beam, the temperature lowering effect is ignored. The relative abundances of 61% and 39% for the gauche and trans conformers are adopted when computing thermally averaged momentum profiles.
3. RESULTS AND DISCUSSION A. Binding energy spectra The B3LYP/aug-cc-pVTZ theory gives the electronic configurations of ethanol as
(core)6 (4a)2 (5a)2 (6a)2 (7a)2 (8a)2 (9a)2 (10a)2 (11a)2 (12a)2 (13a)2 for the gauche conformer and
(core)6 (4a ')2 (5a ')2 (6a ')2 (7a ')2 (1a")2 (8a ')2 (9a ')2 (2a")2 (10a)2 (3a")2 for the trans one, respectively. The two-dimensional electron density map (2D map) of binding energy and relative azimuthal angle of ethanol in our experiment is displayed in Figure 1(a), from which the symmetries of orbitals can be clearly recognized. The total BES obtained by integrating the 2D map over all relative azimuthal angles is shown in Figure 1(b), which is consistent with the previous work15 but with better statistics. The available binding energy region in our measurement mainly covers the whole outer-valence orbitals and two inner-valence ones. Careful analysis must be taken to obtain meaningful momentum profiles as a result of the inherent complication of ionization transitions in the case of the
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
conformationally versatile ethanol molecule. The limited energy resolution of 1.3 eV makes it difficult to resolve ionization bands of all orbitals considering also the overlap due to the two conformers. Efforts are made to resolve separate orbital sets from the angle analysis of experimental ionization intensities. Five resolvable bands (A-E) with remarkable different orbital symmetries in the outer valence region and two complete and independent inner valence bands can be clearly recognized from the BES and 2D map. To give a quantitative analysis of those bands, a deconvolution process employing 9 Gaussian peaks is carried out to infer momentum profiles for particular orbital sets with the help of ionization potentials from high resolution photoelectron spectra (PES)57,58 and high-level theoretical calculations15 shown in Table 3. The peak positions and widths of Gaussian functions mainly come from the PES data, with small adjustments for compensating the asymmetric Franck-Condon envelops and considering molecular conformations and the instrumental energy resolution. In order to have an intuitive sight of how MOs of the two conformers are grouped to contribute to the different ionization bands, a diagram is drawn in Figure 2 to show the ionization potentials of different MOs from the ADC(3) calculation.15 In the outer valence region, band A (peak 1) at ~ 10.7 eV corresponds to HOMO {13a, 3a"}, whose momentum distribution can be easily recognized as s-p type character in the 2D map. Band B (peak 2) at ~ 12.1 eV stands for the next orbitals of the two conformers {12a, 10a'}. Since the momentum profiles of 11a+2a" and 10a+9a' have similar shapes with very high intensity at low momentum while that of 8a' orbital
12 ACS Paragon Plus Environment
Page 12 of 44
Page 13 of 44
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
exhibits typical p type character, we adopt two Gaussian peaks (peak 3 and 4) at 13.5 and 14.4 eV to extract the momentum profile of band C, corresponding to the orbital set {11a, 2a", 10a, 9a', 8a'}, where peak 4 corresponds to the 8a' orbital alone. Analogically, another two Gaussian peaks (peak 5 and 6) at 15.5 and 16.4 eV are used to obtain {9a, 8a, 1a"} momentum profile, referring to band D. Peak 6 is necessary since the shape of momentum profile of the 1a" orbital, whose ionization energy is ~ 0.5 eV larger than that of the 8a orbital, differs from those of the 9a and 8a orbitals. Band E (peak 7), representing {7a, 7a'}, contributes to the shoulder at ~17.5 eV of BES. Band F and G, corresponding to the completely resolved peaks at 20.7 eV and 24.2 eV, characterize {6a, 6a'} and {5a, 5a'}, respectively. Through plotting of the fitted area of the 9 Gaussian peaks in the BES at different azimuthal angles as a function of target electron momentum p (i.e. azimuthal angle), the momentum profiles of the resolvable 7 bands in the valence region can be achieved. Note that momentum distributions for peak 3, 4, 5, or 6 alone is unreliable due to the severe peak overlaps.
B. Momentum profiles Experimental momentum profiles for bands A~G obtained from the deconvolution process described above are presented in Figure 3~7, in which the error bars represents the overall error of deconvolution and statistical uncertainties. The consistency with the previous EMS work15 confirms the reliability of the present experimental outcome. To compare with observed results, the momentum resolution of ∆θ = 0.8º, ∆ϕ = 2.0º is folded in theoretical momentum profiles employing the
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
Gaussian weighted planar grid method.59 Experimental data are normalized to theory, in which the global normalization factor comes from the comparison of the sum of momentum profiles of all outer valence orbitals of ethanol assuming that the one-particle picture of ionization is valid and pole strengths are all identical to 1 which is approximately true according to the ADC(3) calculation. The summed momentum profile of all the outer valence orbitals is presented in Figure 3, together with theoretical simulations. A distinct underestimation at low momentum is shown for the associated equilibrium geometry calculation (blue dashed line), which can be accounted for by vibrational effects unambiguously. On one hand, the gauche and trans conformers have similar momentum distributions either considering vibrational effects or using the equilibrium molecular structure. On the other hand, theoretical calculations at equilibrium geometry predict lower intensity at small momentum noticeably for both conformers. Therefore, the fairly good agreement between theory (red solid line) and experiment (solid circle) for total outer valence orbitals has been achieved when vibrational effects are comprehensively taken into account, exhibiting the significance of nuclear dynamics in the ground state. In the following, a detailed analysis on the relevant ionization bands is performed to reveal how vibrational motions affect EMDs and which modes are responsible for the enhancement of the intensity at low momentum. For the HOMO (13a of gauche and 3a" of trans), the experimental momentum profile of band A exhibits much higher intensity than the equilibrium-geometry B3LYP prediction at small momentum taking account of the 61%/39% thermal
14 ACS Paragon Plus Environment
Page 14 of 44
Page 15 of 44
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
abundances of the gauche and trans conformers, as shown in Figure 4(a), which is in line with the previous EMS works. The large discrepancy at small momentum is obviously decreased with a comprehensive consideration of vibrational effects. It mainly comes from the improvement of intensity at small momentum of the p-type 3a" orbital of trans ethanol. The 3a" orbital, which consists of oxygen 2p lone pair and CH2 bonds, is expected to have zero intensity at zero momentum at equilibrium geometry due to its anti-symmetry character about the CCO plane, keeping in mind that the value of the momentum-space wavefunction at p = 0 reads
ϕ f ( p; Q) p=0 = ( 2π )
−3/2
∫ϕ
f
(r; Q)dr for its Dirac-Fourier transform relation of the
corresponding position space wavefunction. The vibrational modes involving the OH and CH3 torsions (v20, v21), CH2 antisymmetric stretch (v15), and CH2 rock (v18) break the symmetry (Figure 4(d)) and give a positive contribution to the 3a" momentum profile at zero momentum (Figure 4(c)). For the sp-type 13a orbital of gauche ethanol, the contributions of vibrational modes are complicated and nearly smeared out, giving a slight positive contribution on balance (Figure 4(b)). As a result, the deformation of molecular structure by vibrational motions makes a positive modification on the distribution curve with equilibrium geometry at low momentum, which accounts for the discrepancy between experiment and theory principally. The imperfect description of electron correlation for the B3LYP method may be the origin of the remaining deviation since a higher intensity at near zero momentum than that of B3LYP was achieved by the ADC(3) calculation15 which is expected to produce more accurate Dyson orbitals. If we treat the electron correlation and vibrational effects as a first
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
order approximation, the total effects can explain the large turn-up of HOMO rationally. Another possible reason is the ultrafast reorganization of the molecular geometry with a sudden removal of an electron as proposed by Hajgato et al.19 It is worthwhile to compare our result with Hajgato et al.’s thermally averaged momentum profile of HOMO over exceedingly large sets of nuclear configurations in terms of different OH and CH3 torsional angles.19 The HAQM result suggests that the torsional vibrations decrease the intensity of momentum profile at zero momentum by 0.0036 for the 13a orbital and increase that by 0.0113 for the 3a" orbital. While the increase of intensity of the 3a" orbital at p ~ 0 through the torsional modes can be easily understood as analyzed above, readers may ask why there is a decrease of intensity for the 13a orbital. Ning et al.15 demonstrated the asymmetric change of intensity at zero momentum along the OH torsional angle for the 13a orbital which reaches a maximum intensity at 90º. It means that the intensity at p ~ 0 will decrease rapidly when the torsional angle moves towards 0 from the equilibrium geometry (~ 62º), while first rises then falls towards 120º. The OH torsional motion will definitely reduce the electron density at p ~ 0 considering the large-amplitude variations. Consequently, the intensity of HOMO at zero momentum will be improved by 0.0022 if the torsional modes are only involved in the HAQM method taking the 61/39 abundance ratio into account. This agrees with Hajgato et al.’s simulations quite well, proving the effectiveness of the HAQM method for the torsional modes which may suffer from anharmonicity due to its large vibrational amplitude when estimating vibrational effects on EMDs.
16 ACS Paragon Plus Environment
Page 16 of 44
Page 17 of 44
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
Likewise, vibrational motions also have a non-negligible effect on the next orbital of both conformers as displayed in Figure 5(a). Contrary to HOMO, the 12a orbital of the gauche conformer is calculated to show a p-type momentum profile with a maximum at p ~ 0.8 a.u., while the 10a' exhibits an obviously different s-p type momentum distribution. The thermally averaged theoretical momentum distribution at equilibrium geometry underestimates the intensity of band B at p = 0.4~1.0 a.u. noticeably. The discrepancy can be well interpreted by vibrational effects. The intensity in this region is enhanced for both conformers when vibrational motions of all modes are considered. It is interesting that nuclear variations including the OH and CH3 torsions increase the intensity of the momentum profile of the 12a orbital at p < 0.4 a.u. a lot (Figure 5(b)), while decrease that of 10a' of the trans conformer much (Figure 5(c)), and lead to only a small reduction of intensity at low momentum. For bands C and D, the conformation effect on binding energies makes it a little complicated to analyze. Band C is assigned to the orbital sets {11a, 10a, 2a", 9a', 8a'} while band D is assigned to orbital sets {9a, 8a, 1a"}, which are supported by the momentum profiles as shown in Figure 6(a-c). Experimental momentum profile of band C exhibits a much higher intensity at low momentum than the equilibrium geometry calculation. Theoretical calculations considering vibrational effects generally reproduce the experiment well, with a small underestimation and overestimation of intensities at low momentum for bands C and D, respectively. We attribute this small discrepancy to the deconvolution uncertainties due to the overlap of the two bands based on the fact that the HAQM prediction for the sum of bands C
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
and D agrees with the measurement perfectly (Figure 6(c)). This indicates that the EMDs of these orbitals are greatly influenced by the vibrational effects to have a higher intensity at low momentum, which accounts for the severe underestimation of the equilibrium geometry simulation on the observed distribution of band C. In Figure 6(d1-i1), we compare the theoretical momentum profiles for individual orbitals corresponding to bands C and D with and without vibrational effects. The contributions of each vibrational mode are also plotted. It should be noted that the 11a and 10a orbitals of gauche ethanol, which have very close ionization potentials, are coupled seriously by the v2, v4, v9, v14 modes. The mixture of momentum profiles makes the fourth order polynomial fitting of ρ f ( p; Q0 + QL qˆ L ) as a function of QL not suitable for separate orbitals, and only a sum of the two orbitals is given here. It is the same case for the 9a and 8a orbitals by the v12 and v20 modes. Generally, calculations with vibrational effects predict higher intensities near zero momentum for all, especially for the 2a" and 1a" orbitals of trans ethanol. The 2a" orbital, consisted of O 2p lone pair and CH bonds and showing an antisymmetric character about the CCO plane at equilibrium geometry (Figure 2), is expected to exhibit typical p type character. The torsional modes (v20, v21) understandably give a positive contribution as they destroy the symmetry seriously. The situation is the same for the 1a" orbital which is composed of CH bonds. The contributions of vibrational modes for other orbitals are complex and cannot be understood intuitively. Nevertheless, their influences are also significant for well reproducing experimental results. For band E (Figure 7(a)), the inferred momentum profile is well interpreted by
18 ACS Paragon Plus Environment
Page 18 of 44
Page 19 of 44
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
the sum of 7a and 7a' orbitals. The momentum profile of either conformer considering vibrational motions, including the torsional modes and other modes of higher frequency, exhibits higher intensity in the region of less than 1.0 a.u. but a similar shape with that of equilibrium geometry. It should be further noted that significant differences exist in the momentum profiles for every outer valence orbital for the gauche and trans conformers. The wonderful agreement between theory and experiment for the momentum profiles of bands B, C+D, and E, not only indicates the necessity of considering vibrational effects, but also confirms the correction of the 61/39 relative abundance ratio. That is to say, the conclusion that the trans conformer is more stable than the gauche stable is also drawn by the EMS results. In the inner valence region, momentum profiles of the two conformers are predicted to show the same shape and vibrational effects are ignorable for all the relevant orbitals. Theoretical calculations show a good agreement with experiment when pole strengths of 0.87 and 0.76 for the {6a, 6a'} and {5a, 5a'} orbitals respectively are taking into account (Figure 8(a-b)). The pole strength of 0.87 for {6a, 6a'} is in line with the value of 0.90 reported by Ning et al..15 The pole strength of 0.76 for {5a, 5a'} is a higher than Ning et al.’s value of 0.64, while agrees with that of ADC(3) predictions,15 0.70 / 0.91 = 0.77 and 0.67 / 0.91 = 0.74, for the 5a and 5a' orbitals, respectively, considering 0.91 pole strength of HOMO.
4. CONCLUSIONS
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
Page 20 of 44
The electron momentum profiles of valence orbitals of ethanol have been revisited by experiment using a newly conducted EMS as well as the theoretical simulations
involving
conformational
and
vibrational
effects.
The
B3LYP/aug-cc-pVTZ calculations taking into account all the vibrational modes through HAQM approach have quantitatively reproduced the experimental results well, except for the HOMO where observable underestimation still remains, although the agreement between the theory and the experiment has sufficiently been improved. The results unambiguously show that the substantial deviations between the theoretical calculations and measurements in the previous works for the ionization bands of the outer valence orbitals, especially HOMO, are attributed to the vibrational effects. The contributions of individual vibrational modes of the gauche and trans conformers are also analyzed comprehensively using the HAQM approach, which shows that besides the torsional modes, various vibrational modes of high or low frequency influence all outer valence orbitals more or less. In particular, the orbitals of the trans conformer with a" symmetry are affected by the torsional and other modes and exhibit a reasonable turn-up on the p-type momentum distributions at low momentum. Meanwhile, the rational explanation of experimental momentum profiles provides solid evidence that the trans conformer is slightly more stable than the gauche conformer, in accordance with thermodynamical predictions and other experiments. The case of ethanol demonstrates the significance of considering vibrational effects when performing conformational study on flexible molecules using EMS. Finally, it is worthwhile to note that by adopting the supersonic molecular beam
20 ACS Paragon Plus Environment
Page 21 of 44
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
in the future it is possible to perform the experiments at temperature near 0 K. At this temperature, the most stable conformer of ethanol will be dominant. Therefore, one can further confirm the relative stability of the two conformers and can extract the electron momentum profiles for each conformer. The vibrational effects on the individual conformers can then be demonstrated more clearly.
AUTHOR INFORMATION Corresponding Author *E-mail:
[email protected]. *Telephone: +86-551-63601170 Notes The authors declare no competing financial interest.
ACKNOWLEDGEMENT This work was supported by the National Natural Science Foundation of China (Grant Nos. 11534011, 11327404). The authors also gratefully acknowledge Professor C. E. Brion from the University of British Columbia (UBC) in Canada for giving us the HEMS programs.
REFERENCE (1) Brion, C. E. Looking at Orbitals in the Laboratory: The Experimental Investigation of Molecular Wavefunctions and Binding Energies by Electron
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
Momentum Spectroscopy. Int. J. Quantum Chem. 1986, 29, 1397-1428. (2) Coplan, M. A.; Moore, J. H.; Doering, J. P. (e,2e) Spectroscopy. Rev. Mod. Phys. 1994, 66, 985-1014. (3) Mccarthy, I. E.; Weigold, E. Electron Momentum Spectroscopy of Atoms and Molecules. Rep. Prog. Phys. 1991, 54, 789-879. (4) Takahashi, M. Looking at Molecular Orbitals in Three-Dimensional Form: From Dream to Reality. Bull. Chem. Soc. Jpn. 2009, 82, 751-777. (5) Weigold, E.; McCarthy, I. E. Electron Momentum Spectroscopy. Kluwer Academic / Plenum: New York, 1999. (6) Zheng, Y.; Neville, J. J.; Brion, C. E. Imaging the Electron Density in the Highest Occupied Molecular Orbital of Glycine. Science 1995, 270, 786-788. (7) Deleuze, M. S.; Pang, W. N.; Salam, A.; Shang, R. C. Probing Molecular Conformations with Electron Momentum Spectroscopy: The Case of n-Butane. J. Am. Chem. Soc. 2001, 123, 4049-4061. (8) Wang, F. Assessment of Quantum Mechanical Models Based on Resolved Orbital Momentum Distributions of n-Butane in the Outer Valence Shell. J. Phys. Chem. A 2003, 107, 10199-10207. (9) Wang, F.; Downton, M. Inner Valence Shell Bonding Mechanism of n-Butane Studied Using Orbital Momentum Distributions of Its Conformational Isomers. J. Phys. B: At., Mol. Opt. Phys. 2004, 37, 557-569. (10) Saha, S.; Wang, F.; Falzon, C. T.; Brunger, M. J. Coexistence of 1,3-Butadiene Conformers in Ionization Energies and Dyson Orbitals. J. Chem. Phys. 2005, 123,
22 ACS Paragon Plus Environment
Page 22 of 44
Page 23 of 44
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
124315. (11) Deleuze, M. S.; Knippenberg, S. Study of the Molecular Structure, Ionization Spectrum, and Electronic Wave Function of 1,3-Butadiene Using Electron Momentum Spectroscopy and Benchmark Dyson Orbital Theories. J. Chem. Phys. 2006, 125, 104309. (12) Huang, Y. R.; Knippenberg, S.; Hajgato, B.; Francois, J. P.; Deng, J. K.; Deleuze, M. S. Imaging momentum orbital densities of conformationally versatile molecules: A benchmark theoretical study of the molecular and electronic structures of dimethoxymethane. J. Phys. Chem. A 2007, 111, 5879-5897. (13) Yang, T. C.; Su, G. L.; Ning, C. G.; Deng, J. K.; Wang, F.; Zhang, S. F.; Ren, X. G.; Huang, Y. R. New Diagnostic of the Most Populated Conformer of Tetrahydrofuran in the Gas Phase. J. Phys. Chem. A 2007, 111, 4927-4933. (14) Morini, F.; Hajgato, B.; Deleuze, M. S.; Ning, C. G.; Deng, J. K. Benchmark Dyson Orbital Study of the Ionization Spectrum and Electron Momentum Distributions of Ethanol in Conformational Equilibrium. J. Phys. Chem. A 2008, 112, 9083-9096. (15) Ning, C. G.; Luo, Z. H.; Huang, Y. R.; Hajgato, B.; Morini, F.; Liu, K.; Zhang, S. F.; Deng, J. K.; Deleuze, M. S. Investigation of the Molecular Conformations of Ethanol Using Electron Momentum Spectroscopy. J. Phys. B: At., Mol. Opt. Phys. 2008, 41, 175103. (16) Wu, F.; Chen, X. J.; Shan, X.; Tian, S. X.; Li, Z. J.; Xu, K. Z. Conformational Stability of 1-Butene: An Electron Momentum Spectroscopy Investigation. J. Phys.
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
Chem. A 2008, 112, 4360-4366. (17) Xue, X. X.; Yan, M.; Wu, F.; Shan, X.; Xu, K. Z.; Chen, X. J. Electron Momentum Spectroscopy of Ethanethiol Complete Valence Shell. Chin. J. Chem. Phys. 2008, 21, 515-520. (18) Chen, X. J.; Wua, F.; Yan, M.; Li, H. B.; Tian, S. X.; Shan, X.; Wang, K. D.; Li, Z. J.; Xu, K. Z. Hyperconjugative Effect on the Electronic Wavefunctions of Ethanol. Chem. Phys. Lett. 2009, 472, 19-22. (19) Hajgato, B.; Deleuze, M. S.; Morini, F. Probing Nuclear Dynamics in Momentum Space: A New Interpretation of (e, 2e) Electron Impact Ionization Experiments on Ethanol. J. Phys. Chem. A 2009, 113, 7138-7154. (20) Luo, Z. H.; Ning, C. G.; Liu, K.; Huang, Y. R.; Deng, J. K. Electron Momentum Spectroscopy Study of A Conformationally Versatile Molecule: n-Propanol. J. Phys. B: At., Mol. Opt. Phys. 2009, 42, 165205. (21) Yan, M.; Shan, X.; Wu, F.; Xia, X.; Wang, K. D.; Xu, K. Z.; Chen, X. J. Electron Momentum Spectroscopy Study on Valence Electronic Structures of Ethylamine. J. Phys. Chem. A 2009, 113, 507-512. (22) Shojaei, S. H. R.; Filippo, M.; Bálazs, H.; Michael, S. D. Photoelectron and Electron Momentum Spectroscopy of 1-Butene at Benchmark Theoretical Levels. J. Phys. B: At., Mol. Opt. Phys. 2011, 44, 235101. (23) Shi, Y. F.; Shan, X.; Wang, E. L.; Yang, H. J.; Zhang, W.; Chen, X. J. Experimental and Theoretical Investigation on the Outer Valence Electronic Structure of Cyclopropylamine by (e, 2e) Electron Momentum Spectroscopy. J. Phys. Chem. A
24 ACS Paragon Plus Environment
Page 24 of 44
Page 25 of 44
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
2014, 118, 4484-4493. (24) Shi, Y. F.; Shan, X.; Wang, E. L.; Yang, H. J.; Zhang, W.; Chen, X. J. Electron Momentum Spectroscopy of Outer Valence Orbitals of 2-Fluoroethanol. Chin. J. Chem. Phys. 2015, 28, 35-42. (25) Huang, C. W.; Shan, X.; Zhang, Z.; Wang, E. L.; Li, Z. J.; Chen, X. J. High-Resolution Electron Momentum Spectroscopy of Valence Satellites of Carbon Disulfide. J. Chem. Phys. 2010, 133, 124303. (26) Miao, Y. R.; Ning, C. G.; Liu, K.; Deng, J. K. Dyson orbitals of N2O: Electron Momentum Spectroscopy and Symmetry Adapted Cluster-Configuration Interaction Calculations. J. Chem. Phys. 2011, 134, 204304. (27) Tian, Q.; Yang, J.; Shi, Y.; Shan, X.; Chen, X. Outer- and Inner-Valence Satellites of Carbon Dioxide: Electron Momentum Spectroscopy Compared with Symmetry-Adapted-Cluster Configuration Interaction General-R Calculations. J. Chem. Phys. 2012, 136, 094306. (28) Ren, X. G.; Ning, C. G.; Deng, J. K.; Zhang, S. F.; Su, G. L.; Huang, F.; Li, G. Q. Direct Observation of Distorted Wave Effects in Ethylene Using the (e,2e) Reaction. Phys. Rev. Lett. 2005, 94, 163201. (29) Ning, C. G.; Ren, X. G.; Deng, J. K.; Su, G. L.; Zhang, S. F.; Li, G. Q. Turn-up Effects at Low Momentum for the Highest Occupied Molecular Orbital of Oxygen at Various Impact Energies by Electron Momentum Spectroscopy. Phys. Rev. A 2006, 73, 022704. (30) Liu, K.; Ning, C. G.; Luo, Z. H.; Shi, L. L.; Deng, J. K. An Experimental and
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
Theoretical Study of the HOMO of W(CO)6: Vibrational Effects on the Electron Momentum Density Distribution. Chem. Phys. Lett. 2010, 497, 229-233. (31) Hajgató, B.; Morini, F.; Deleuze, M. S. Electron Momentum Spectroscopy of Metal Carbonyls: A Reinvestigation of the Role of Nuclear Dynamics. Theor. Chem. Acc. 2012, 131, 1-15. (32) Miao, Y. R.; Deng, J. K.; Ning, C. G. Vibrational Effects on the Electron Momentum Distributions of Valence Orbitals of Formamide. J. Chem. Phys. 2012, 136, 124302. (33) Watanabe, N.; Yamazaki, M.; Takahashi, M. Vibrational Effects on Valence Electron Momentum Distributions of Ethylene. J. Chem. Phys. 2012, 137, 114301. (34) Watanabe, N.; Yamazaki, M.; Takahashi, M. Vibrational Effects on Valence Electron Momentum Distributions of CH2F2. J. Chem. Phys. 2014, 141, 244314. (35) Shojaei, S. H. R.; Vandenbussche, J.; Deleuze, M. S.; Bultinck, P. Electron Momentum Spectroscopy of 1-Butene: A Theoretical Analysis Using Molecular Dynamics and Molecular Quantum Similarity. J. Phys. Chem. A 2013, 117, 8388-8398. (36) Yang, J.; Shan, X.; Zhang, Z.; Tang, Y.; Zhao, M.; Chen, X. Ring-Puckering Effects on Electron Momentum Distributions of Valence Orbitals of Oxetane. J. Phys. Chem. A 2014, 118, 11780-11786. (37) Morini, F.; Deleuze, M. S.; Watanabe, N.; Takahashi, M. Theoretical Study of Molecular Vibrations in Electron Momentum Spectroscopy Experiments on Furan: An Analytical Versus A Molecular Dynamical Approach. J. Chem. Phys. 2015, 142,
26 ACS Paragon Plus Environment
Page 26 of 44
Page 27 of 44
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
094308. (38) Morini, F.; Watanabe, N.; Kojima, M.; Deleuze, M. S.; Takahashi, M. Electron Momentum Spectroscopy of Dimethyl Ether Taking Account of Nuclear Dynamics in the Electronic Ground State. J. Chem. Phys. 2015, 143, 134309. (39) Farasat, M.; Shojaei, S. H. R.; Morini, F.; Golzan, M. M.; Deleuze, M. S. Electron Momentum Spectroscopy of Aniline Taking Account of Nuclear Dynamics in the Initial Electronic Ground State. J. Phys. B: At., Mol. Opt. Phys. 2016, 49, 075102. (40) Tang, Y.; Shan, X.; Yang, J.; Niu, S.; Zhang, Z.; Watanabe, N.; Yamazaki, M.; Takahashi, M.; Chen, X. Vibrational Effects on Electron Momentum Distributions of Outer-Valence Orbitals of Oxetane. J. Phys. Chem. A 2016, 120, 6855-6863. (41) Li, Z. J.; Chen, X. J.; Shan, X.; Liu, T.; Xu, K. Z. Electron Momentum Spectroscopy Study of Jahn-Teller Effect in Cyclopropane. J. Chem. Phys. 2009, 130, 054302. (42) Zhu, J. S.; Miao, Y. R.; Deng, J. K.; Ning, C. G. The Jahn-Teller Effect in the Electron Momentum Spectroscopy of Ammonia. J. Chem. Phys. 2012, 137, 174305. (43) Senent, M. L.; Smeyers, Y. G.; Domı́nguez-Gómez, R.; Villa, M. Ab Initio Determination of the Far Infrared Spectra of Some Isotopic Varieties of Ethanol. J. Chem. Phys. 2000, 112, 5809-5819. (44) Durig, J. R.; Deeb, H.; Darkhalil, I. D.; Klaassen, J. J.; Gounev, T. K.; Ganguly, A. The r0 Structural Parameters, Conformational Stability, Barriers to Internal Rotation, and Vibrational Assignments for Trans and Gauche Ethanol. J. Mol. Struct.
27 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
2011, 985, 202-210. (45) Kakar, R. K.; Quade, C. R. Microwave Rotational Spectrum and Internal Rotation in Gauche Ethyl Alcohol. J. Chem. Phys. 1980, 72, 4300-4307. (46) Pearson, J. C.; Sastry, K. V. L. N.; Herbst, E.; De Lucia, F. C. The Millimeterand Submillimeter-Wave Spectrum of Gauche-Ethyl Alcohol. J. Mol. Spectrosc. 1996, 175, 246-261. (47) Durig, J. R.; Larsen, R. A. Torsional Vibrations and Barriers to Internal Rotation for Ethanol and 2,2,2-Trifluoroethanol. J. Mol. Struct. 1990, 238, 195-222. (48) Abu-samha, M.; Børve, K. J. Franck-Condon Transitions in A System with Large-Amplitude Anharmonic Vibrations Coupled to A Harmonic-Oscillator Bath: Application to the C1s Photoelectron Spectrum of ethanol. Phys. Rev. A 2006, 74, 042508. (49) Tian, Q.; Wang, K.; Shan, X.; Chen, X. A high-Sensitivity Angle and Energy Dipersive Multichannel Electron Momentum Spectrometer with 2π Angle Range. Rev. Sci. Instrum. 2011, 82, 033110. (50) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula Into A Functional of The Electron Density. Phys. Rev. B 1988, 37, 785-789. (51) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648-5652. (52) Dunning, T. H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989,
28 ACS Paragon Plus Environment
Page 28 of 44
Page 29 of 44
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
90, 1007-1023. (53) Kendall, R. A.; Dunning, T. H.; Harrison, R. J. Electron Affinities of the First-Row Atoms Revisited. Systematic Basis Sets and Wave Functions. J. Chem. Phys. 1992, 96, 6796-6806. (54) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Rob, M. A.; Cheeseman, J. R.; Montgomery, Jr., J. A..; Vreven, T.; Kudin, K. N.; Burant, J. C., et al. Gaussian 03. Gaussian, Inc.: Wallingford, CT, 2003. (55) Fábri, C.; Szidarovszky, T.; Magyarfalvi, G.; Tarczay, G. Gas-Phase and Ar-Matrix SQM Scaling Factors for Various DFT Functionals with Basis Sets Including Polarization and Diffuse Functions. J. Phys. Chem. A 2011, 115, 4640−4649. (56) Duffy, P.; Chong, D. P.; Casida, M. E.; Salahub, D. R. Assessment of Kohn-Sham Density-Functional Orbitals as Approximate Dyson Orbitals for the Calculation of Electron-Momentum-Spectroscopy Scattering Cross Sections. Phys. Rev. A 1994, 50, 4707-4728. (57) von Niessen, W.; Bieri, G.; Åsbrink, L. 30.4-nm He (II) Photoelectron Spectra of Organic Molecules. J. Electron Spectrosc. Relat. Phenom. 1980, 21, 175-191. (58) Kimura, K.; Katsumata, S.; Achiba, Y.; Yamazaki, T.; Iwata, S. Handbook of HeI Photoelectron Spectra of Fundamental Organic Molecules. Halsted Press: New York, 1981. (59) Duffy, P.; Casida, M. E.; Brion, C. E.; Chong, D. P. Assessment of Gaussian-Weighted
Angular
Resolution
Functions
in
29 ACS Paragon Plus Environment
the
Comparison
of
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
Quantum-Mechanically
Calculated
Electron
Momentum
Experiment. Chem. Phys. 1992, 159, 347-363.
30 ACS Paragon Plus Environment
Page 30 of 44
Distributions
with
Page 31 of 44
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
Tables -1 Table 1. Vibrational frequencies (cm ) for gauche ethanol
a
mode
Main description a
Cal.
Exp.a
1
OH stretch
3812
3662
2
CH3 antisymmetric stretch
3098
2994
3
CH3 antisymmetric stretch
3084
2987
4
CH2 antisymmetric stretch
3063
2972
5
CH3 symmetric stretch
3018
2936
6
CH2 symmetric stretch
2988
2912
7
CH2 symmetric deformation
1517
1493
8
CH3 antisymmetric deformation
1493
1465
9
CH3 antisymmetric deformation
1488
1460
10
CH3 symmetric deformation
1415
1394
11
CH2 wag
1403
1373
12
COH bend
1368
1342
13
CH2 twist
1279
1249
14
CH2 rock/CH3 rock
1132
1117
15
CO stretch
1067
1066
16
CH3 rock/COH bend
1057
1058
17
CC stretch
882
879
18
CH3 rock
805
803
19
CCO bend
420
420
20
OH torsion
269
243
21
CH3 torsion
257
274
Reference 44.
31 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
Table 2. Vibrational frequencies (cm-1) for trans ethanol
a
Mode
Main description a
Cal.
Exp.a
1
OH stretch (A')
3828
3676
2
CH3 antisymmetric stretch (A')
3098
2992
3
CH3 symmetric stretch (A')
3032
2922
4
CH2 symmetric stretch (A')
2980
2888
5
CH2 symmetric deformation (A')
1526
1500
6
CH3 antisymmetric deformation (A')
1500
1480
7
CH2 wag (A')
1447
1450
8
CH3 symmetric deformation (A')
1405
1367
9
COH bend (A')
1264
1241
10
CO stretch (A')
1096
1090
11
CC stretch (A')
1028
1027
12
CH3 rock (A')
895
892
13
CCO bend (A')
418
418
14
CH3 antisymmetric stretch (A")
3101
2987
15
CH2 antisymmetric stretch (A")
3004
2901
16
CH3 antisymmetric deformation (A")
1483
1455
17
CH2 twist (A")
1299
1275
18
CH2 rock (A")
1179
1166
19
CH3 rock (A")
822
801
20
OH torsion (A")
276
203
21
CH3 torsion (A")
233
244
Reference 44.
32 ACS Paragon Plus Environment
Page 32 of 44
Page 33 of 44
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
Table 3. Ionization potentials of ethanol. The values in brackets are pole strengths, which are assumed to be 1 for the outer valence orbitals in EMS MO
He II PES a
He I PES b
EMS c
EMS d
A
13a/3a"
10.7
10.64
10.6(1.0)
B
12a/10a'
12.1
12.18
12.1(1.0)
11a/2a"
13.3
13.21
13.4(1.0)
10a/9a'
13.9
13.86
8a'
(14.5)
14.5
14.4(1.0)
16.0
15.85
15.9(1.0)
Band
C
9a D
8a 1a"
ADC(3)/cc-pVDZ++ c gauche
trans
10.7(1.0)
10.86(0.91)
10.89(0.91)
12.1(1.0)
12.28(0.92)
12.39(0.91)
13.51(0.92)
13.37(0.92)
13.87(0.91)
13.77(0.92)
13.5(1.0) 14.4(1.0) 15.5(1.0)
14.43(0.91) 15.34(0.91) 15.97(0.91)
16.4(1.0)
E
7a/7a'
17.4
F
6a/6a'
20.7
17.35
17.5(1.0)
17.4(1.0)
20.8(0.90)
20.7(0.87)
16.53(0.90) 17.51(0.90)
17.90(0.90)
21.20(0.85)
20.80(0.85)
23.86(0.03) 23.33(0.05) G
5a/5a'
24.2
24.2(0.64)
24.2(0.76)
23.35(0.08)
24.39(0.25)
24.48(0.14)
24.55(0.27)
24.60(0.02)
24.64(0.18)
24.63(0.26) 24.77(0.09) a
Reference 56. Reference 57. Reference 15. Present work. b
c
d
33 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 Captions Figure1. (a) Two-dimensional energy and relative azimuthal angle spectrum of ethanol. (b) Binding energy spectrum of ethanol. The dotted lines represent the nine Gaussian peaks corresponding to individual molecular orbital sets and the solid line is the sum. Figure 2. Diagram of molecular orbitals of the two conformers of ethanol. The orbital maps are plotted at the contour of 0.05, with different colors representing opposite phases. Ionization potentials are from the ADC(3) calculation.15 The curves above the ionization potential lines are theoretical momentum profiles in the momentum region of 0 ~ 3.5 a.u.by B3LYP/aug-cc-pVTZ. Figure 3. The sum of electron momentum profiles of outer valence orbitals of ethanol. Equilib. stands for the calculation at equilibrium geometry, and vib. means the computation involving vibrational effects. Figure 4. Electron momentum profiles of Band A. (a) Experimental momentum profiles with theoretical ones considering conformational and vibrational effects. (b) Contributions from each vibrational mode to the momentum profile of the 13a orbital of gauche ethanol. (c) Contributions from each vibrational mode to the momentum profile of the 3a" orbital of trans ethanol with a position-space orbital map at equilibrium geometry inside. (d) Nuclear displacement vectors of selected vibrational modes in which the red sphere represents the oxygen atom, and orbital maps of the 3a" orbital of trans ethanol distorted along those modes 2
2
ξv (QL ) / ξv (0) = 1/ 4 . The orbital maps are plotted at a contour of 0.05, with different L
L
colors representing opposite phases. Figure 5. Electron momentum profiles of Band B. (a) Experimental momentum profiles with
34 ACS Paragon Plus Environment
Page 34 of 44
Page 35 of 44
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
theoretical ones considering conformational and vibrational effects. (b) Contributions from each vibrational mode to the momentum profile of the 12a orbital of gauche ethanol. (c) Contributions from each vibrational mode to the momentum profile of the 10a' orbital of trans ethanol. Figure 6. Experimental momentum profiles with theoretical ones considering conformational and vibrational effects for band (a) C, (b) D and (c) C+D, together with theoretical momentum profiles at equilibrium geometry and considering vibrational effects for (d1) 11a+10a, (e1) 2a", (f1) 9a', (g1) 8a', (h1) 9a+8a and (i1) 1a" orbitals as well as contributions from each vibrational mode to the momentum profiles of these orbital (d2-i2). Note that the momentum resolution is only folded in theoretical momentum profiles in (a-c). Figure 7. Electron momentum profiles of Band E. (a) Experimental momentum profiles with theoretical ones considering conformational and vibrational effects. (b) Contributions from each vibrational mode to the momentum profile of the 7a orbital of gauche ethanol. (c) Contributions from each vibrational mode to the momentum profile of the 7a' orbital of trans ethanol. Figure 8. Experimental momentum profiles with theoretical ones considering conformational and vibrational effects for (a) band F and (b) band G in the inner valence region.
35 ACS Paragon Plus Environment
Page 36 of 44
1.000
40
(a)
2.078 4.317
20
8.970 18.64
0
38.72 80.45
-20
167.2 347.3
-40
721.6 1175
10
15
1 20.0k
Relative Intensity
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
Relative Azimuthal Angle (degree)
The Journal of Physical Chemistry
3 4 5 6 7
25
8
30
9
Ethanol
C
G
D
15.0k 10.0k
2
(b)
20
B E
A
F
5.0k 0.0 10
15
20
25
Binding Energy (eV)
ACS Paragon Plus Environment
30
Page 37 of 44
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
Gauche
10
Hf(eV)
Trans
13a
3a"
12a 11a 10a 9a 8a
10a' 2a" 9a' 8a'
15
1a"
7a
7a' 20
6a'
6a
5a
25
5a'
ACS Paragon Plus Environment
The Journal of Physical Chemistry
0.6
Band A+B+C+D+E 0.61×gauche+0.39×trans vib. 0.61×gauche vib. 0.39×trans vib. 0.61×gauche+0.39×trans equilib. 0.61×gauche equilib. 0.39×trans equilib.
0.5
Relative Intensity
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
Page 38 of 44
0.4 0.3 0.2 0.1 0.0 0.0
0.5
1.0
1.5
2.0
2.5
Momentum (a.u.)
ACS Paragon Plus Environment
3.0
3.5
Page 39 of 44
0.07
0.008
(a)
Band A 0.6113a+0.393a" vib. 0.6113a vib. 0.393a" vib. 0.6113a+0.393a" equilib. 0.6113a equilib. 0.393a" equilib.
Relative Intensity
0.05 0.04 0.03
Contributions of vibrational modes v6 v10
0.004
v12
v14
v19
v20
0.002
v21 sum of other modes total
0.000
0.02
-0.002
0.01 0.00 0.0
13a gauche ethanol
(b)
0.006
Intensity
0.06
0.5
1.0
1.5
2.0
2.5
3.0
3.5
-0.004 0.0
0.5
1.0
Momentum (a.u.) 0.020
(c)
1.5
2.0
2.5
3.0
Momentum (a.u.)
(d)
3a" trans ethanol Contributions of vibrational modes v15
0.015
v18
Intensity
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
The Journal of Physical Chemistry
v20
0.010
v21
v15
sum of other modes total
0.005
v18
0.000 0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Momentum (a.u.)
ACS Paragon Plus Environment
v20
v21
3.5
The Journal of Physical Chemistry
0.10
(a)
Band B 0.6112a+0.3910a' vib. 0.6112a vib. 0.3910a' vib. 0.6112a+0.3910a' equilib. 0.6112a equilib. 0.3910a' equilib.
Relative Intensity
0.08
0.06
0.04
0.02
0.00 0.0 0.030
0.5
1.0
(b)
0.025
1.5
2.0
3.0
3.5
Contributions of vibrational modes v4 v18 v20
Intensity
2.5
12a gauche ethanol
0.020
v21
sum of other modes total
0.015 0.010 0.005 0.000 0.0 0.02
0.5
1.0
1.5
2.0
2.5
3.0
3.5
(c)
0.00
Intensity
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 40 of 44
10a' trans ethanol
-0.02
Contributions of vibrational modes v14 v20
-0.04
v21 sum of other modes total
-0.06
-0.08 0.0
0.5
1.0
1.5
2.0
2.5
3.0
Momentum (a.u.)
ACS Paragon Plus Environment
3.5
Page 41 of 44
0.40 Band C 0.61(11a+10a) +0.39(2a"+9a'+8a') vib. 0.61(11a+10a) vib. 0.39(2a"+9a'+8a') vib. 0.61(11a+10a) +0.39(2a"+9a'+8a') equilib. 0.61(11a+10a) equilib. 0.39(2a"+9a'+8a') equilib.
Relative Intensity
0.16
0.12
0.08
0.16
(b)
Band D 0.61(9a+8a)+0.391a" vib. 0.61(9a+8a) vib. 0.391a" vib. 0.61(9a+8a)+0.391a" equilib. 0.61(9a+8a) equilib. 0.391a" equilib.
0.12
0.08
0.35
0.04
0.04
(c)
Band C+D 0.61×gauche+0.39×trans vib. 0.61×gauche vib. 0.39×trans vib. 0.61×gauche+0.39×trans equilib. 0.61×gauche vib. 0.39×trans vib.
0.30
Relative Intensity
(a)
Relative Intensity
0.20
0.25 0.20 0.15 0.10 0.05
0.00 0.0
0.5
1.0
1.5
2.0
2.5
3.0
0.00 0.0
3.5
0.5
1.0
Momentum (a.u.) 0.20
0.09
gauche ethanol 11a+10a vib. 11a+10a equilib. 11a equilib. 10a equilib.
(e1)
2.5
3.0
0.00 0.0
3.5
0.5
1.0
0.10
0.07
1.5
2.0
2.5
3.0
3.5
4.0
Momentum (a.u.)
2a" trans ethanol equilib. vib.
0.08
(f1)
9a' trans ethanol equilib. vib.
0.08
0.06
Intensity
Intensity
0.12
2.0
0.08
Intensity
(d1)
1.5
Momentum (a.u.)
0.10
0.16
0.05 0.04
0.06
0.04
0.03
0.04
0.02
0.02
0.01
(d2)
2.0
2.5
3.0
0.00 0.0
3.5
0.04
Contributions of vibrational modes v7 v13 v15
0.010
1.0
1.5
2.0
2.5
3.0
0.00 0.0
3.5
sum of other modes total
0.5
1.0
1.5
2.0
2.5
3.0
3.5
0.01
sum of other modes total
0.24
0.5
1.0
1.5
(h1)
2.0
2.5
3.0
-0.03 0.0
3.5
v11
v16
0.12
0.5
1.0
1.5
2.0
2.5
3.0
3.5
3.0
3.5
Momentum (a.u.)
0.20
(i1)
gauche ethanol 9a+8a vib. 8a+8a equilib. 9a equilib. 8a equilib.
0.16
0.04
3.5
sum of other modes total
0.20
Intensity
0.06
3.0
Contributions of vibrational modes v2 v6
-0.01
v9
Momentum (a.u.)
8a' trans ethanol equilib. vib.
2.5
9a' trans ethanol
-0.02
Momentum (a.u.)
(g1)
2.0
0.00
v21
0.00
0.08
1.5
Momentum (a.u.)
0.01
v20
0.02
-0.01 0.0
1.0
(f2)
Contributions of vibrational modes v14 v16 v18
0.5
0.02
2a" trans ethanol
v20
0.000
Intensity
(e2)
0.03
v21
0.005
-0.005 0.0
0.5
Momentum (a.u.)
11a+10a gauche ethanol
0.015
Intensity
1.5
Momentum (a.u.)
Intensity
1.0
trans ethanol 1a" vib. 1a" equilib.
0.16
0.12
Intensity
0.020
0.5
Intensity
0.00 0.0
0.08
0.08 0.02 0.04
0.04
0.5
1.0
1.5
2.0
2.5
3.0
0.00 0.0
3.5
0.5
1.0
Momentum (a.u.)
(g2)
8a' trans ethanol
0.02
0.02
0.01
2.0
2.5
3.0
0.00 0.0
3.5
0.5
1.0
v9
v6
v8
v10
sum of other modes total
0.01
1.5
2.0
2.5
Momentum (a.u.)
9a+8a gauche ethanol Contributions of vibrational modes v7 v14
(h2)
Contributions of vibrational modes v2 v7 Intensity
0.03
1.5
Momentum (a.u.)
v16
v17
v20
v21
0.02
Intensity
0.00 0.0
Intensity
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
sum of other modes total 0.00
1a" trans ethanol
(i2)
Contributions of vibrational modes v15 v20
0.01
v21 sum of other modes total
0.00
0.00
-0.01 0.0
0.5
1.0
1.5
2.0
Momentum (a.u.)
2.5
3.0
3.5
-0.01 0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Momentum (a.u.)
ACS Paragon Plus Environment
-0.01 0.0
0.5
1.0
1.5
2.0
Momentum (a.u.)
2.5
3.0
3.5
The Journal of Physical Chemistry
0.12
(a)
Band E 0.617a+0.397a' vib. 0.617a vib. 0.397a' vib. 0.617a+0.397a' equilib. 0.617a equilib. 0.397a' equilib.
Relative Intensity
0.10 0.08 0.06 0.04 0.02 0.00 0.0 0.012 0.010
0.5
1.0
(b)
1.5
v15
0.006
Intensity
2.0
2.5
3.0
3.5
7a gauche confomer Contributions of vibrational modes v6 v7
0.008
v20
v21
0.004
sum of other modes total
0.002 0.000 -0.002 0.0 0.012
0.5
1.0
(c)
1.5
2.0
2.5
3.0
3.5
7a' trans confomer Contributions of vibrational modes v11 v13
0.008
v15 Intensity
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 42 of 44
v20
v21
0.004
sum of other modes total 0.000
-0.004 0.0
0.5
1.0
1.5
2.0
2.5
3.0
Momentum (a.u.)
ACS Paragon Plus Environment
3.5
Page 43 of 44
0.14 0.12
(a)
Band F 0.616a+0.396a' vib. 0.87 0.616a vib. 0.87 0.396a' vib. 0.87 0.616a+0.396a' equilib. 0.87 0.616a equilib. 0.87 0.396a' equilib. 0.87
Relative Intensity
0.10 0.08 0.06 0.04 0.02 0.00 0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Momentum (a.u.) 0.35
(b)
Band G 0.615a+0.395a' vib. 0.76 0.615a vib. 0.76 0.395a' vib. 0.76 0.615a+0.395a' equilib. 0.76 0.615a equilib. 0.76 0.395a' equilib. 0.76
0.30
Relative Intensity
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
The Journal of Physical Chemistry
0.25 0.20 0.15 0.10 0.05 0.00 0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Momentum (a.u.)
ACS Paragon Plus Environment
The Journal of Physical Chemistry
0.07
HOMO Equilibrium geometry Vibrational effects
0.06
Relative Intensity
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
0.05 0.04 0.03
Gauche 13a 61%
0.02
Trans 3a" 39%
0.01 0.00 0.0
0.5
1.0
1.5
2.0
2.5
3.0
Momentum (a.u.)
ACS Paragon Plus Environment
3.5
Page 44 of 44