Conformational Stability from Variable-Temperature Infrared Spectra

May 27, 2011 - James R. Durig,* Ahmed M. El-Defrawy, Arindam Ganguly, Savitha S. Panikar, ... Department of Chemistry, University of Missouri-Kansas C...
0 downloads 0 Views 1MB Size
ARTICLE pubs.acs.org/JPCA

Conformational Stability from Variable-Temperature Infrared Spectra of Xenon Solutions, r0 Structural Parameters, and Vibrational Assignment of Pyrrolidine† James R. Durig,* Ahmed M. El-Defrawy, Arindam Ganguly, Savitha S. Panikar, and Mamdouh S. Soliman‡ Department of Chemistry, University of Missouri-Kansas City, Kansas City, Missouri 64110, United States

bS Supporting Information ABSTRACT: The infrared spectra of gaseous and variable-temperature liquid xenon solutions of pyrrolidine have been recorded. The enthalpy difference has been determined to be 109 ( 11 cm1 (1.30 ( 0.13 kJ mol1) with the envelope-equatorial conformer more stable than the twist form with 37 ( 3% present at ambient temperature. Ab initio calculations utilizing various basis sets up to MP2(full)/aug-ccpVTZ have been used to predict the conformational stabilities, energy at the equatorialaxial saddle point, and barriers to planarity. From previously reported microwave rotational constants along with MP2(full)/6-311þG(d,p) predicted structural values, adjusted r0 parameters have been obtained for both conformers. Heavy atom distances (Å) of equatorial[twist] conformer are as follows: N1C2 = 1.469(3)[1.476(3)], N1C3 = 1.469(3)[1.479(3)], C2C4 = 1.541(3)[1.556(3)], C3C5 = 1.541(3)[1.544(3)], C4C5 = 1.556(3)[1.543(3)]; and angles (deg) — N1C2C4 = 102.5(5)[107.6(5)], — N1C3C5 = 102.5(5)[105.4(5)], — C2C4C5 = 104.3(5)[104.6(5)], — C3C5C4 = 104.3(5)[103.7(5)], — C2N1C3 = 104.1(5)[103.9(5)], τC2C4C5C3 = 0.0(5)[13.5(5)]. A complete vibrational assignment is proposed for both conformers.

’ INTRODUCTION The saturated five-membered ring molecules are structurally interesting, since they can be in one or more of the three stable conformers, i.e., twisted, envelope, and planar forms. Furthermore, the energy difference between these forms can be very small with the barrier for the molecule to go from one form to another sufficiently small, where instead of going through the planar form, it may go through a series of twisted forms where one atom moves out perpendicular to a hypothetical plane of the rest of the atoms. This concept was called pseudorotation14 and was first experimentally observed for cyclopentane.4 As a continuation of our earlier investigations of five-membered rings, we more recently initiated variable-temperature infrared studies of rare gas solutions of several molecules some of which we had investigated earlier.5,6 We have determined the enthalpy difference between the two stable, envelope-axial and envelopeequatorial, forms of cyclopentyl chloride and cyclopentyl bromide by this method and the experimental values are 145 ( 15 cm1 (1.73 ( 0.18 kJ mol1) and 233 ( 23 cm1 (2.79 ( 0.28 kJ mol1), respectively, with the axial conformer as the more stable form for both molecules. Also, from studies of the variable-temperature infrared spectra of rare gas solutions of silacyclopentane7 and germacyclopentane,8 we found these molecules have only one stable twist (C2) conformer present in all physical states. From the initial vibrational study9,10 of cyclopentyl fluoride it was concluded that two conformers were present in the fluid states and they were the envelope-axial and envelope-equatorial conformers with the equatorial form the r 2011 American Chemical Society

more stable conformer. However, from a later vibrational investigation11 of the Raman spectrum of the liquid, it was concluded that there was only one form present and it was the equatorial form. This conclusion was supported by a theoretical prediction by a CNDO/2 calculation12 which predicted only the equatorial conformer as the stable form. However, we recently carried out13 a study of the variable temperature infrared spectrum of xenon solutions of cyclopentyl fluoride and determined that there was only one conformer present in the fluid states but it was the twisted (C1) form. This conclusion was supported by ab initio predictions as well as by the determined microwave rotational constants. Therefore, it is expected that the conformational conclusions of the structural stabilities of some of the other substituted five-membered rings may be in error. Another saturated five-membered ring molecule of considerable conformational interest is pyrrolidine, c-C4H8NH, which has received much scientific interest due to its presence in biologically important molecules such as peptides, proteins, and amino acids. It has been the subject of several experimental and theoretical conformational studies. From the first conventional microwave study, Caminati et al.14 reported that the pyrrolidine molecule exists in an envelope-axial form but there was no spectroscopic evidence for the envelope-equatorial form. In a subsequent electron diffraction study,15 which was supported Received: January 22, 2011 Revised: May 25, 2011 Published: May 27, 2011 7473

dx.doi.org/10.1021/jp200692t | J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A by HF/4-21N* ab initio calculations, again the envelope-axial form was reported to be the more stable conformer with the envelope-equatorial conformer higher in energy by 339 cm1 (4.05 kJ mol1) and the barrier to pseudorotation was predicted to be 580 cm1 (6.94 kJ mol1). However, from a more recent microwave study, Caminati et al.16 discovered transitions of the equatorial conformer during the investigation of rotational free jet spectrum of pyrrolidinewater adduct. In this microwave investigation, it was observed that the axial conformer relaxes to the equatorial conformer when the cooling conditions were increased, which clearly shows that the equatorial form is the more stable conformer. To estimate the enthalpy difference from the relative intensity of the microwave lines the dipole moments were predicted by utilizing the vectorial compositions of the bond moments. An arbitrary uncertainty value of the dipole moments was taken to be 0.4 D. These predicted dipole moment values for both conformers were utilized (only the |μc|) which gave an estimated enthalpy difference to be ΔE0,0 (=Eax  Eeq) = 80 ( 300 cm1. However, by combining this value with the cooling effects of the jet which shows the equatorial conformer to be more stable, the predicted enthalpy difference of 0 < ΔE0,0 < 200 cm1 was obtained.16 In order to obtain an experimentally determined enthalpy difference, we initiated an investigation of the vibrational spectrum of pyrrolidine with a study of the infrared spectra of the gas and xenon solutions at variable temperatures. To obtain reliable enthalpy determinations from variable-temperature infrared spectrum in xenon solutions, it is essential to have a confident vibrational assignment for both forms particularly in the spectral region where the conformer pair or pairs are selected for the enthalpy determination. Therefore, assignments need to be made for the fundamentals for each of the conformers which is facilitated by utilizing ab initio calculations at the MP2 level with full electron correlation by the perturbation method.17 In the last three vibrational investigations,1820 no assignments were made for any of the fundamentals of the less stable conformer. Additionally, the structure of the less stable conformer has not been definitively determined since it was proposed to be the axial form as the most stable conformer from the initial microwave14 and electron diffraction studies15 and it seems to be accepted without experimental verifications. Also, we were interested in obtaining the r0 structural parameters for both the stable conformers. The results for these spectroscopic and theoretical studies are reported herein.

’ EXPERIMENTAL AND THEORETICAL METHODS The sample of pyrrolidine was purchased from Aldrich Chemical Co., with stated purity of 98%. The sample was further purified by a low-temperature, low-pressure fractionation column, and the purity of the sample was verified by comparing the infrared spectrum with previously reported ones.1820 The mid-infrared spectrum of the gas was obtained from 3500 to 300 cm1 on a Perkin-Elmer model 2000 Fourier transform spectrometer equipped with a Ge/CsI beamsplitter and a DTGS detector. Atmospheric water vapor was removed from the spectrometer housing by purging with dry nitrogen. The theoretical resolution used to obtain the spectrum of the gas was 0.5 cm1. Sixty-four interferograms were added and transformed with a boxcar truncation function. The mid-infrared spectrum of the solid was obtained with the sample deposited on a silicon substrate at 77 K, and multiple annealings were performed in order to obtain a good polycrystalline sample.

ARTICLE

The far-infrared spectrum of the sample was recorded on the previously described Perkin-Elmer 2000 spectrometer. A grid beamsplitter was used to record the spectrum of the gas with the sample contained in a 12 cm cell equipped with polyethylene windows. The spectra were recorded at a spectral resolution of 0.5 cm1 and, typically, 256 scans were used for both the sample and the reference data to give a satisfactory signal-to-noise ratio. The interferograms were averaged and then transformed with a boxcar truncation function. All of the observed bands of the fundamentals for both conformers are listed in Tables 1 and 2. The mid-infrared spectra (3500400 cm1) (Figure 1A) of the sample dissolved in liquefied xenon as a function of temperature, ranging from 55 to 100 C, were recorded on a Bruker model IFS-66 Fourier transform spectrometer equipped with a globar source, a Ge/KBr beamsplitter, and a DTGS detector. For all spectra, 100 interferograms were collected at 1.0 cm1 resolution, averaged, and transformed with a boxcar truncation function. For these studies, a specially designed cryostat cell was used. It consists of a copper cell with a path length of 4 cm with wedged silicon windows sealed to the cell with indium gaskets. The copper cell was enclosed in an evacuated chamber fitted with KBr windows. The temperature was maintained with boiling liquid nitrogen and monitored by two Pt thermoresistors. The LCAO-MO-SCF restricted HartreeFock calculations were performed with the Gaussian 03 program21 by using Gaussian-type basis functions. The energy minima with respect to nuclear coordinates were obtained by the simultaneous relaxation of all geometric parameters by using the gradient method of Pulay.22 A variety of basis sets with and without diffuse functions were employed with the MøllerPlesset perturbation method17 to the second-order (MP2(full)) as well as with the density functional theory by the B3LYP method. The predicted conformational energy differences are listed in Table 3. The infrared spectra were predicted from the MP2(full)/ 6-31G(d) calculations. The predicted scaled frequencies were used together with a Lorentzian function to obtain the calculated spectra. Infrared intensities determined from MP2(full)/ 6-31G(d) calculations were obtained based on the dipole moment derivatives with respect to Cartesian coordinates. The derivatives were transformed with respect to normal coordinates by (∂μu/∂Qi) = ∑j(∂μu/∂Xj)Lij, where Qi is the ith normal coordinate, Xj is the jth Cartesian displacement coordinate, and Lij is the transformation matrix between the Cartesian displacement coordinates and the normal coordinates. The infrared intensities were then calculated by (Nπ)/(3c2)[(∂μx/ ∂Qi)2 þ (∂μy/∂Qi)2 þ (∂μz/∂Qi)2]. In Figure 1, a comparison of experimental and simulated infrared spectra of pyrrolidine is shown. The predicted spectrum is in good agreement with the experimental spectrum which shows the utility of the scaled predicted frequencies and predicted intensities for supporting the vibrational assignment. In order to obtain a complete description of the molecular motions involved in the fundamental modes of pyrrolidine, a normal-coordinate analysis was carried out. The force field in Cartesian coordinates was obtained with the Gaussian 03 program at the MP2(full) level with the 6-31G(d) basis set. The internal coordinates used to calculate the G and B matrices are given in Table 4 with the atomic numbering shown in Figure 2. By using the B matrix,23 the force field in Cartesian coordinates was converted to a force field in internal coordinates. Subsequently, scaling factors of 0.88 for CH2 stretches and CH2 7474

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

ARTICLE

Table 1. Calculateda and Observed Frequencies (cm1) for Envelope-Equatorial (Cs) Form of Pyrrolidine infrared sym block A

0

A00

vib no.

approx description

ab initio

fixed scaledb

IR intensity

band contour

gas

liq Xe

PEDc

A

B

C

ν1

NH str

3535

3353

0.4

3351

3348

100S1

2



98

ν2

β-CH2 antisym str

3200

3002

45.2

2986

2985

98S2

4



96

ν3

R-CH2 antisym str

3158

2962

9.6

2971

2965

60S3, 33S5

2



98

ν4

β-CH2 sym str

3138

2944

34.0

2944

2943

94S4

100





ν5

R-CH2 sym str

3015

2828

115.2

2828

2826

66S5, 34S2

3



97

ν6

R-CH2 def

1591

1492

0.8

1476

1475

77S6, 21S7

3



97

ν7

β-CH2 def

1568

1471

4.4

1463

1458

78S7, 20S6

71



29

ν8 ν9

R-CH2 wag β-CH2 wag

1439 1364

1365 1294

1.8 1.9

1351 1288

1348 1286

79S8 58S9, 26S10

81 96

 

19 4

ν10

R-CH2 twist

1284

1218

0.8

1222

1221

52S10, 17S9, 12S13

12



88

ν11

β-CH2 twist

1268

1202

8.2

1203

1197

55S11, 16S13

1



99

ν12

ring def

1111

1054

0.7

1054

1041

18S12, 32S14, 14S9,

1



99

ν13

R-CH2 rock

1041

987

6.3

984

977

30S13, 21S12, 16S16, 17S11

3



97

ν14

ring def/breathing

982

931

2.5

933

929

27S14, 37S15, 12S11

44



56

ν15

ring breathing/def

941

893

6.5

896

898

58S15, 22S14

7



93

ν16 ν17

NH in-plane bend/ring def β-CH2 rock

921 798

874 757

69.5 22.9

792 739

791 733

36S16, 47S12 70S17, 10S16

30 45

 

70 55

ν18

ring def

596

565

43.3

568

568

76S18, 12S17, 10S13

41



59

ν19

ring puckering

318

302

4.9

302

288

98S19

100





ν20

β-CH2 antisym str

3182

2985

0.9

2981

2980

95S20



100



ν21

R-CH2 antisym str

3155

2959

52.9

2967

2966

62S21, 32S23



100



ν22

β-CH2 sym str

3130

2937

16.0

2937

2936

98S22



100



ν23

R-CH2 sym str

3013

2826

35.1

2826

2824

67S23, 33S21



100



ν24 ν25

R-CH2 def β-CH2 def

1573 1547

1476 1451

1.2 0.9

1472 1450

1464 1449

93S24 93S25

 

100 100

 

ν26

NH out-of-plane bend

1469

1394

5.6

1400

1399

56S26, 26S28, 10S33



100



ν27

β-CH2 wag

1366

1296

9.3

1283

1282

48S27, 22S29, 13S26



100



ν28

R-CH2 wag

1341

1272

15.0

1281

1281

42S28, 33S27, 13S29



100



ν29

β-CH2 twist

1283

1217

7.7

1208

1209

57S29, 13S34, 11S28



100



ν30

R-CH2 twist

1233

1170

0.8

1174

1172

63S30, 10S31



100



ν31

ring def

1169

1109

10.0

1114

1112

42S31, 25S32



100



ν32 ν33

R-CH2/β-CH2 rocks ring def

1140 971

1081 921

0.5 0.3

1080 916

1076 910

30S32, 29S34, 18S31, 68S33, 12S31

 

100 100

 

ν34

β-CH2/R-CH2 rocks

900

854

2.7

855

853

35S34, 31S32, 12S35



100



ν35

ring def

639

606

0.6

600

598

73S35, 10S32



100



ν36

ring twist

67

64

0.2

65



95S36



100



MP2(full)/6-31G(d) ab initio calculations, scaled frequencies, infrared intensities (km mol1), and potential energy distributions (PEDs). str = stretch; def = deformation. b Scaled frequencies with scaling factors of 0.88 for CH2 stretches and CH2 deformations and 0.90 for all other modes. c Symmetry coordinates with PED contribution less than 10% are omitted. a

deformations and 0.90 for all other coordinates were applied, along with the geometric average of scaling factors for interaction force constants, to obtain the fixed scaled force field and resultant wavenumbers. Comparisons between the observed and calculated wavenumbers, along with the calculated infrared intensities and potential energy distributions for equatorial and twist forms, are listed in Tables 1 and 2, respectively.

’ VIBRATIONAL ASSIGNMENT To obtain reliable enthalpy determinations from variabletemperature infrared spectra, it is essential to have confident assignments for both forms particularly in the spectral region where the

conformer pairs are being selected for the enthalpy determination. With significantly more information than what was utilized in making the earlier vibrational assignment such as ab initio predicted intensities, fundamental frequencies for the second conformer, and gas-phase contours, we have utilized all of these data as well as infrared spectra from xenon solutions where the sharp bands make it possible to identify closely spaced fundamentals. We expected to be able to assign all of the fundamentals for the more stable conformer and many of those for the second conformer particularly in the “fingerprint” spectral region (Figure 3). The assignments in the carbonhydrogen stretching region are nearly the same as those given earlier18 for the gas of the more stable conformer but it is possible to assign several of the 7475

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

ARTICLE

Table 2. Calculateda and Observed Frequencies (cm1) for Twist (C1) Form of Pyrrolidine infrared vib no.

approx description

ab initio

fixed scaledb IR intensity gas

band contour PEDc

liq. Xe

ν1

NH str

3500

3320

0.4

3324

ν2

β-CH2 antisym str

3190

2993

64.7

(2986)

ν3

R-CH2 antisym str

3181

2984

17.1

2981

ν4

R-CH2 antisym str

3174

2977

0.8

2971

ν5

β-CH2 antisym str

3164

2968

12.3

2961

ν6

β-CH2 sym str

3128

2934

28.2

ν7

β-CH2 sym str

3114

2922

ν8 ν9

R-CH2 sym str R-CH2 sym str

3113 3108

2920 2915

ν10

β-CH2 def

1580

ν11

R-CH2 def

1562

3322

100S1

A

B

C

15

2

83

3

1

96

72S3, 18S5

7

81

12

57S4, 31S2

10

67

23

76S5, 16S3

12

85

3

2932

78S6, 16S7

70

30



41.4

2922

67S7, 18S6

20

58

22

20.3 34.6

2925 2909

53S8, 24S9, 16S7 63S9, 27S8

39 6

14 67

47 27

1482

1.1

(1467)

63S10, 35S11

91

6

3

1465

3.6

1457

54S11, 36S10

22

3

75

65S12, 24S13, 10S11

2

90

8

70S13, 27S12

4



96

65S2, 30S4

1456

ν12

R-CH2 def

1558

1462

1.8

(1457)

ν13

β-CH2 def

1549

1453

0.5

1463

ν14

NH out-of-plane bend

1500

1423

7.3

1419

1419

78S14, 11S21

1

99



ν15

R-CH2 wag

1410

1338

2.7

1341

1336

84S15

88



12

ν16 ν17

R-CH2 twist R-CH2 twist

1386 1371

1315 1301

0.6 1.7

1315 1302

1315 1297

24S16, 24S18,15S17,12S20, 11S19 48S17, 12S16,10S18, 10S19

88 90

5 10

7 

ν18

β-CH2 wag

1340

1271

1.0

1271

1269

ν19

β-CH2 wag

1329

1260

1.0

1261

1457

55S18, 36S16



99

1

34S19,15S27, 10S16

18

81

1

ν20

β-CH2 twist

1292

1226

2.0

1226

53S20, 16S22, 10S27

30

70



ν21

R-CH2 wag

1255

1190

3.4

1188

1182

36S21, 13S19, 10S22

85

14

1

ν22

β-CH2 twist/R-CH2 wag

1247

1183

4.7

1186

1182

17S22, 24S21, 13S31, 10S20

94

3

3

ν23

ring def

1148

1089

20.3

1089

1083

64S23, 22S28



100



ν24 ν25

ring def β-CH2/R-CH2 rocks

1083 1080

1028 1025

0.7 2.4

1020 1013

39S24, 11S29, 11S19 30S25, 25S31, 17S21, 12S14

35 4

41 95

24 1

ν26

R-CH2 rock/β-CH2 twist

1028

975

25.0

975

962

27S26, 27S22, 20S24

94

1

5

ν27

ring def/breathing

973

923

10.5

923

919

25S27, 52S29, 18S30

94



6

ν28

ring def

958

908

0.5

909

59S28, 19S23

81

17

2

ν29

ring breathing/NH bend

930

883

8.1

885

883

24S29, 51S30, 12S24

98





ν30

NH in-plane bend/R-CH2

894

848

84.0

848

842

15S30, 31S26, 12S34, 11S24

90



10

784

15S31, 23S25, 18S26, 12S32 61S32, 13S34, 10S25

84 3

1 23

15 74 23

rock ν31 ν32

R-CH2/β-CH2 rocks β-CH2 rock

884 817

838 775

33.5 1.8

838 787

ν33

ring def

664

630

1.7

630

39S33, 32S34

28

49

ν34

ring def

616

584

0.9

591

28S34, 30S33, 13S32, 10S27

21

64

15

ν35

ring puckering

308

292

8.0

288

91S35

21



79

ν36

ring twist

38

36

0.3

39

91S36

5

28

67

MP2(full)/6-31G(d) ab initio calculations, scaled frequencies, infrared intensities (km mol1), and potential energy distributions (PEDs). str = stretch; def = deformation. b Scaled frequencies with scaling factors of 0.88 for CH2 stretches and CH2 deformations and 0.90 for all other modes. c Symmetry coordinates with PED contribution less than 10% are omitted. a

fundamentals to the second conformer in this spectral region. Beginning with the CH2 deformations and the remaining bands in the “fingerprint” region, the assignments for many of the fundamentals of the more stable conformer are different than those given earlier since many of the corresponding vibrations of the second conformer are quite pronounced. Most of the CH2 bending vibrations have expected frequencies based on group frequencies but there is significant mixing among the CH2 rocks and heavy atom ring vibrations. There is also considerable mixing with the NH bends which makes the description sometimes arbitrary. Nevertheless, three-fourths of the vibrations have major contributions of 50% or greater so that the descriptions

are useful for comparison of corresponding vibrations in similar molecules. On this basis, the assignments of several of the ring modes are of interest with the ring breathing vibrations at 896 and 885 cm1 with the latter one for the second conformer but extensively mixed with the NH bend. This latter band was assigned earlier18 to this vibration but the other four bands previously assigned as ring modes in the 8001100 cm1 region do not agree with the proposed assignment presented in this study. Also, only one of the lower frequency assigned ring modes agrees with those previously proposed18 but it should be noted that the earliest assigned ones made based on ab initio prediction were for the axial conformer as the most stable form. After the 7476

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

ARTICLE

second microwave investigation was reported16 with the equatorial conformer as the more stable form, ab initio predictions were used to estimate the frequencies for the equatorial form and the bands from the earlier reported spectra were reassigned to the equatorial form but none of them were assigned as arising from the less stable conformer. The other assignment reported18 for

Figure 1. Comparison of experimental and calculated infrared spectra of pyrrolidine: (A) observed infrared spectrum in xenon; (B) simulated infrared spectrum of a mixture of envelope-equatorial and twist conformers at 100 C with ΔH = 109 cm1; (C) simulated infrared spectrum of twist (C1) form; (D) simulated infrared spectrum of equatorial (Cs) form.

the equatorial form was made from the spectrum of the liquid where significant association is expected due to hydrogen bonding of the amine. Therefore, bands at 923, 909, and 885 cm1 are now assigned as ring modes for the second conformer, whereas, all three bands had been assigned18 as fundamentals for the equatorial form with two of them ring modes. Thus, only one of the five ring modes above 800 cm1 assigned in the current study agrees with those assigned earlier18 but two of the four ring modes below 700 cm1 agree with those assigned18 earlier. It should be noted that for the equatorial form only one of the vibrations has contributions from four symmetry coordinates but 10 of them have significant contributions from three symmetry coordinates. This mixing makes it impossible to assign the NH in-plane bend, CH2 rocks, and three of the ring deformations to mainly a single vibration which has resulted in several different vibrations contributing significantly to the NH bend. Therefore, some of the descriptions of these modes are more for bookkeeping than to convey a description of the major atom motions involved. Although the descriptions of several of the band contours in the fingerprint spectral region are ambiguous, definitive assignments for most of them to one or the other of the conformers can be confidently made. This is possible from the spectral data of the xenon solution where the bands are relatively sharp and the band centers for the broad B-type bands for the gas can be clearly identified. For example, the broad band in the 10001050 cm1 region of the xenon solution has clear bands at 1020 and 1013 cm1 which can be assigned to the second conformer whereas the band at 1041 cm1 must be assigned to the equatorial form. This example shows that the relatively weak broad bands can be assigned to the individual conformers. With the vibrations which are assigned to the individual conformers and not badly overlapped, it should be possible to obtain the enthalpy difference between the conformers from variable-temperature studies of the infrared spectra of the xenon solutions. Relatively small interactions are expected to occur between xenon and the sample and, thus, small frequency shifts are anticipated when passing from gas to the liquefied xenon solutions. A significant advantage of this

Table 3. Calculated Electronic Energies (Hartrees) and Energy Differencesa (cm1) for the Envelope-Equatorial (Cs), Twist (C1), Envelope-Axial (Transition State) (Cs), and Planar (Transition State) (Cs) Conformers of Pyrrolidine method/basis set MP2(full)/6-31G(d) MP2(full)/6-31þG(d)

no. of basis functions

envelopeequatorial

twist

ΔEEeqT

ΔEEeqEax

planar

ΔEEeqP

93

211.850047

211.850122

16

211.850115

15

211.840341

2130

113

211.863522

211.862925

131

211.862888

139

211.853918

2108

MP2(full)/6-311G(d,p)

144

212.070241

212.070741

110

212.070548

67

212.060696

2095

MP2(full)/6-311þG(d,p)

164

212.077306

212.077220

19

212.077050

56

212.068326

1971

MP2(full)/6-311G(2d,2p)

196

212.132264

212.132212

11

212.132159

23

212.122841

2068

MP2(full)/6-311þG(2d,2p) MP2(full)/6-311G(2df,2pd)

216 276

212.137810 212.217237

212.137472 212.216998

74 52

212.137446 212.216869

80 81

212.128537 212.207550

2035 2126

MP2(full)/6-311þG(2df,2pd)

296

212.222209

212.221736

104

212.221664

119

212.212711

2085

MP2(full)/6-311þþG(2df,2pd)

305

212.222513

212.222009

111

212.221944

125

212.212996

2089

MP2(full)/aug-cc-pVTZ

565

212.232642

212.232319

71

212.232291

77

212.223483

82 ( 24

averageb

96 ( 24

2010 2069 ( 46

93

212.581750

212.582016

58

212.582023

60

212.574451

1602

B3LYP/6-31þG(d)

113

212.591062

212.590634

94

212.590641

92

212.583712

1613

B3LYP/6-311G(d,p) B3LYP/6-311þG(d,p)

144 164

212.642045 212.646368

212.642152 212.646034

23 73

212.642160 212.646043

25 71

212.634865 212.639520

1576 1503

B3LYP/6-31G(d)

a

envelopeaxial

B3LYP/6-311G(2d,2p)

196

212.651253

212.651322

15

212.651331

17

212.644739

1430

B3LYP/6-311þG(2d,2p)

216

212.655218

212.654989

50

212.655000

48

212.648727

1425

Energy differences of the twist, axial, and planar forms are relative to the equatorial form. b Average values are for the five largest basis sets. 7477

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

ARTICLE

Table 4. Structural Parameters (Å and deg) and Rotational Constants (MHz) for Equatorial and Twist Forms of Pyrrolidine adjusted r0

MP2(full)/6-311þG(d,p) struct parameters

a

a

internal coord

equatorial

twist

ED axial

equatorial

twist

r(NC2)

R1

1.463

1.473

1.469(10)

1.469(3)

1.476(3)

r(NC3)

R2

1.463

1.466

1.469(10)

1.469(3)

1.479(3)

r(C2C4) r(C3C5)

R3 R4

1.536 1.536

1.557 1.535

1.543(8) 1.543(8)

1.541(3) 1.541(3)

1.556(3) 1.544(3)

r(C4C5)

R5

1.551

1.542

1.543(8)

1.556(3)

1.543(3)

r(NH6)

r1

1.016

1.017

1.020

1.016(2)

1.017(2)

r(C2H7)

r2

1.104

1.093

1.090(4)

1.104(2)

1.093(2)

r(C2H9)

r3

1.093

1.094

1.090(4)

1.093(2)

1.094(2)

r(C3H10)

r4

1.104

1.093

1.090(4)

1.104(2)

1.093(2)

r(C3H8)

r5

1.093

1.096

1.090(4)

1.093(2)

1.097(2)

r(C5H11) r(C5H13)

r6 r7

1.092 1.092

1.095 1.093

1.090(4) 1.090(4)

1.092(2) 1.092(2)

1.095(2) 1.093(2)

r(C4H14)

r8

1.092

1.093

1.090(4)

1.092(2)

1.093(2)

r(C4H12)

r9

1.092

1.094

1.090(4)

1.092(2)

1.093(2)

— N1C2C4

ϕ1

102.5

107.7

104.6

102.5(5)

107.6(5)

— N1C3C5

ϕ2

102.5

105.4

104.6

102.5(5)

105.4(5)

— C2C4C5

ϕ3

104.2

104.2

104.9

104.3(5)

104.6(5)

— C3C5C4

ϕ4

104.2

102.8

104.9

104.3(5)

103.7(5)

— C2N1C3 — C2N1H6

ψ σ1

103.8 111.8

104.1 108.4

105.2(35) 107.0

104.1(5) 111.7(5)

103.9(5) 108.4(5)

— C3N1H6

σ2

111.8

107.9

111.7(5)

108.0(5)

— H7C2N1

R1

112.0

110.0

112.0(5)

110.0(5)

— H9C2N1

R3

111.1

108.3

111.1(5)

108.3(5)

— H10C3N1

R2

112.0

110.7

112.0(5)

110.7(5)

— H8C3N1

R4

111.1

108.2

111.1(5)

108.2(5)

— H7C2C4

R5

109.6

112.9

109.6(5)

113.0(5)

— H9C2C4 — H10C3C5

R6 R7

113.1 109.6

110.1 113.9

113.1(5) 109.6(5)

110.1(5) 113.9(5)

— H8C3C5

R8

113.1

110.2

113.1(5)

110.2(5)

— H7C2H9

γ1

108.5

107.7

108.5(5)

107.7(5)

— H10C3H8

γ2

108.5

108.2

108.5(5)

108.2(5)

— H11C5C3

β1

111.6

108.9

111.6(5)

109.0(5)

— H13C5C3

β3

110.1

113.6

110.1(5)

113.6(5)

— H14C4C2

β2

111.6

112.0

111.6(5)

112.0(5)

— H12C4C2 — H11C5C4

β4 β5

110.1 112.6

110.6 110.4

110.1(5) 112.5(5)

110.6(5) 110.0(5)

— H13C5C4

β6

110.4

112.9

110.4(5)

112.4(5)

— H14C4C5

β7

112.6

111.8

112.5(5)

111.6(5)

— H12C4C5

β8

110.4

110.8

110.4(5)

110.6(5)

— H11C5H13

ε1

107.8

108.1

107.9(5)

108.1(5)

— H14C4H12

ε2

107.8

107.5

107.9(5)

107.5(5)

τN1C2C4C5

τ1

27.5

4.25

27.0(5)

10.0(5)

τN1C3C5C4 τC2C4C5C3

τ2 τ3

27.5 0.0

36.9 19.2

27.0(5) 0.0

32.8(5) 13.5(5)

τC2N1C3C5

τ4

46.5

40.1

45.8(5)

39.4(5)

τC3N1C2C4

τ5

46.5

27.3

45.8(5)

30.6(5)

A

6914.6

6884.4

6865.2

6833.9

B

6828.3

6715.9

6792.6

6678.6

C

3934.6

3924.9

3902.1

3888.9

39.0(14)

Values taken from ref 15. 7478

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

ARTICLE

Figure 2. Planar (Cs) conformer of pyrrolidine showing atom numbering.

Figure 4. Temperature (55 to 100 C) dependent mid-infrared spectrum in the 525600 cm1 region of pyrrolidine dissolved in liquid xenon.

Figure 3. Mid-infrared spectra (5001500 cm1) of pyrrolidine: (A) gas in transmittance; (B) liquid xenon solution in absorbance.

study is that the conformer bands are better resolved in comparison with those in the infrared spectrum of the gas. From ab initio calculations, the dipole moments of the two conformers are predicted to have similar values (Table 3) and the molecular sizes of the rotamers are nearly the same, so that the ΔH value obtained from the temperature-dependent FT-IR study is expected to be comparable to the value for the gas.2428

’ CONFORMATIONAL STABILITY Once the vibrational assignments were made for most of the observed bands in the spectral region from 400 to 1400 cm1, the bands to be used to obtain the enthalpy difference need to be selected. Lowest frequency bands are best selected since the number of overtone or combination bands will increase significantly for the higher wavenumber bands. Thus, the first fundamental considered for the equatorial conformer was ν18 which is mainly due to a ring deformation. However, the second conformer has a fundamental on the high-frequency side as well as ν35 for the equatorial conformer. Also, there is a distinct band on the low-frequency side with a pronounced Q branch which is tentatively assigned as the overtone of ν35 of the less stable conformer. There is also a second weaker band slightly below the overtone band and the breadth of ν18 makes it difficult to obtain the relative intensity change with the decrease in temperature. Thus, the value of the uncertainty when this band is used is larger than normally obtained. Nevertheless, this fundamental is

mainly (76%) due to a ring mode which is desirable since bands associated with an NH mode can be significantly affected by association in the solution. Another fundamental for the equatorial conformer is the 1041 cm1 band which has a highfrequency shoulder but it is sufficiently separated so its intensity can be reasonably well measured. The third equatorial fundamental chosen was the 1112 cm1 band which is relatively sharp in the xenon solution and has significant contribution from S31(42%), a ring mode. The final two bands chosen were the 791 (47S12, ring deformation) and 1282 cm1 where both have sufficient intensity to be used with the most intense band being used for the second conformer. The bands chosen for the second conformer were the very pronounced 1083 cm1 which is a ring deformation and the weak doublet at 1013 and 1020 cm1 with the latter one also a ring deformation. The five equatorial bands and the three from the second conformer made it possible to obtain six individual enthalpy differences. The intensities of the chosen infrared bands were measured as a function of temperature (Figure 4) and their ratios were determined. By application of the van’t Hoff equation ln K = ΔH/(RT)  ΔS/R, the enthalpy differences were determined from a plot of ln K versus 1/T, where ΔH/R is the slope of the line and K is substituted with the appropriate intensity ratios, i.e., IformA/IformB. It was assumed that ΔH is not a function of temperature in the range studied. By combining the intensity of the bands of the equatorial conformer with those of the second form, six individual ΔH values were obtained where different values ranged from a high value of 134 cm1 to the lowest one of 71 cm1 but with different statistical uncertainties ranging from 41 to 6 cm1 with an arithmetic average of 24 cm1. However, by combining these six pairs into a single set the determined ΔH value 7479

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

ARTICLE

is 109 ( 11 cm1 (1.30 ( 0.13 kJ mol1) (Table 5). This error limit was derived from the statistical standard deviation of 1σ of the measured intensity data where the data from the six pairs were taken as a single set. These error limits do not take into account small associations with the liquid xenon or the interference of overtones and combination bands in near coincidence with the measured fundamentals. The variations are undoubtedly due to these types of interferences but by taking six pairs it is hoped that the effect of such interferences would be minimized. From these results the abundance of the less stable conformer at ambient temperature is estimated to be 37 ( 3%. The experimentally determined enthalpy difference is in satisfactory agreement with the energy difference predicted from the ab initio calculations for a second conformer which could be either the envelope-axial or the twist form. However, the ab initio calculations have a negative frequency for the axial form which indicates that it is a first-order saddle point and not a stable conformer. However, the energy differences are relatively small from the calculations with the predictions from the MP2(full)/ aug-cc-pVTZ expected to be the best estimate.

changes, 11 of the 15 rotational constants are fit to 0.4 MHz or better with the other four differing by 0.6, 0.7, and 1.3 MHz (Table 6). The estimated uncertainties for the individual parameters are listed in Table 4. For the second conformer which is predicted from the ab initio calculations to be the twisted form with no symmetry, there are Table 6. Comparison of Rotational Constants (MHz) Obtained from ab Initio MP2(full)/6-311þG(d,p) Predictions, Microwave Spectra, and Adjusted r0 Structural Parameters for Pyrrolidine conformer

isotopomer

equatorial

C4H8NH

rotational const

obsda

calcd

|Δ|

A

6864.7

6865.2

0.6

B

6791.9

6792.6

0.7

C

3902.3

3902.1

0.2

A

6865.6

6865.2

0.4

B C

6405.8 3771.1

6405.6 3771.3

0.2 0.2

A

6803.2

6803.6

0.4

B

6709.7

6709.6

0.1

C

3858.5

3857.2

1.3

A

6845.1

6845.4

0.3

B

6661.0

6661.0

0.0

C

3851.6

3852.3

0.7

A B

6856.1 6665.0

6855.9 6665.3

0.2 0.3

C

3862.3

3862.7

0.4

A

6834.5

6833.9

0.6

B

6677.8

6678.5

0.7

C

3888.1

3888.9

0.8

A

6707.4

6707.6

0.2

B

6403.8

6404.5

0.7

C A

3834.6 6735.6

3834.3 6735.7

0.3 0.1

B

6635.8

6634.9

0.9

C

3843.9

3843.2

0.7

A

6803.6

6804.3

0.7

B

6560.0

6559.0

1.0

C

3838.4

3838.6

0.2

C4H8ND

2-13C

’ STRUCTURAL PARAMETERS We29 have shown that ab initio MP2(full)/6-311þG(d,p) calculations predict the r0 structural parameters for 50 substituted hydrocarbons for the carbonhydrogen distances better than 0.002 Å compared to the experimentally determined values from isolated CH stretching frequencies30 which agree with previously determined values from earlier microwave studies. Therefore, all of the carbonhydrogen distances can be taken from the MP2(full)/6-311þG(d,p) predicted values for pyrrolidine. We have also shown31 that we can obtain good structural parameters by adjusting the structural parameters obtained from the ab initio calculations to fit the rotational constants obtained from microwave experimental data. For the pyrrolidine equatorial conformer there are three heavy atom distances, three heavy atom angles, and two dihedral angles to be determined from the 12 microwave determined (excluding the N-D isotopomer) rotational32 constants by adjusting the predicted parameters from the MP2(full)/6-311þG(d,p) calculation. The adjusted r0 parameters are given in Table 4, and the three heavy atom distances are increased by 0.006 and 0.005 Å and — CCC increased by 0.1, and the — CNC is increased by 0.3 and one of the dihedral angles changes by 0.7. With these very small

3-13C

C4H815NH

C4H8NH

twist

C4H8ND

2-13C

3-13C

a

Values for the rotational constants taken from ref 32.

Table 5. Temperature and Intensity Ratios of the Equatorial and Twist Bands of Pyrrolidine T (C) liquid xenon

ΔHa (cm1)

1/T (103 K1)

I568/I1083 

I791/I1083

I1282/I1083

I1112/I1020

I1041/I1020

I1041/I1013

55.0

4.584





0.360

2.409

2.154

60.0

4.692

1.333

0.899

0.309

0.357

2.207

1.919

65.0

4.804

1.280

0.842

0.314

0.345

2.442

1.921

70.0 75.0

4.923 5.047

1.250 1.231

0.823 0.881

0.283 0.332

0.345 0.337

2.154 2.355

1.824 2.010

80.0

5.177

1.250

0.809

0.274

0.332

2.189

1.931

85.0

5.315

1.232

0.845

0.303

0.327

2.276

1.846

90.0

4.460

1.121

0.802

0.271

0.326

1.979

1.742

95.0

5.613

1.118

0.766

0.262

0.322

2.025

1.778

100.0

5.775

1.068

0.765

0.263

0.320

1.813

1.615

130 ( 16

88 ( 21

119 ( 41

71 ( 6

134 ( 32

116 ( 27

a Average value ΔH = 109 ( 11 cm1 (1.30 ( 0.13 kJ mol1) with the equatorial conformer the more stable form and the statistical uncertainty (1σ) obtained by utilizing all of the data as a single set.

7480

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A now five heavy atom distances, five angles, and five dihedral angles which must be determined from the 12 microwave reported rotational constants. In order to reduce the number of independent variables, the structural parameters are separated into sets according to their types. Bond lengths in the same set keep their relative ratio which results in only three heavy atom distances for the twisted form. Also, the bond angles in the same set keep their differences in degrees which reduces them to three. This assumption is based on the fact that the errors from ab initio calculations are systematic. By this process there are 11 parameters to be adjusted so it should be possible to obtain “adjusted r0” structural parameters for the twisted conformer by utilizing the 12 determined microwave rotational constants from this conformer. The parameters obtained for the twist conformer are listed in Table 4 along with the estimated uncertainties of the parameters. From the microwave studies14,16 of pyrrolidine the centrifugal distortion constants were reported for both conformers. From the predicted force constants by utilizing two different basis sets from MP2(full) ab initio calculations and with the same basis sets from density functional calculations, the centrifugal distortion constants have been calculated for the envelope-equatorial and twist conformers (Table S1 in the Supporting Information). In the case of the envelope-equatorial conformer, only the ΔJ and δJ values were obtained from the microwave study16 and the remaining quartic centrifugal distortion parameters were fixed to zero because they could not be determined from the experimental data. The three distortion constants which were set to zero have been predicted to have significant values (Table S1) which makes the values reported for the other two distortion constants with statistical uncertainties unrealistic. The experimental distortion constants for the twist conformer were originally attributed incorrectly to the envelope-axial conformer.14 This misassignment has been rectified in Table S1 but the agreement of the predicted values from the MP2(full)/6-311þG(d,p) calculations is only satisfactory for the ΔJ and ΔK constants but the other three constants have very large differences (2 orders of magnitude).

’ DISCUSSION The frequency predictions from the ab initio calculations for the A0 modes of the envelope-equatorial conformer are on the average within 5.6 cm1 which represents an error of 0.4%. For the A00 modes for this conformer the predictions are within 4.9 cm1 which is an error of 0.3%. These predictions were obtained with only two scaling factors with the ab initio calculation from the relatively small 6-31G(d) basis set by second-order perturbation with full electron correlation. Similar predictions are expected for the second conformer, so the present assignment for this less stable form is expected to be very good. For the envelope-equatorial conformer which has a plane of symmetry, the mixing of the normal modes is relatively minor except for the two highest frequency ring modes and the NH bend in the A0 symmetry block. Except for these four modes the approximate descriptions given represent a reasonable description of the major atomic motions for the vibrations. However, for the twisted form where there are no symmetry elements, several of the modes are given approximate descriptions more for bookkeeping than to convey the major atomic displacements. For example, ν16 has contributions from five different symmetry coordinates with contributions greater than 10% and six of them with contributions from four symmetry coordinates. Therefore, the

ARTICLE

Figure 5. ν30 and ν31 modes of the twist form predicted and observed at 848 and 838 cm1of the gas.

descriptions for the modes for the twisted form are not as meaningful as those for the equatorial conformer. We investigated the difference in the predicted frequencies for the axial and twist forms and found most of the values were only 2 or 3 cm1 apart except for ν9, ν17, ν18, ν29, ν35, and ν36 which differ by 5 or 6 cm1 (Table S2 in the Supporting Information). However, the ν35 and ν18 which are ring modes are predicted to be 20 cm1 lower and 16 cm1 higher, respectively, for the axial form than the corresponding modes for the twist form. It should be noted that these modes are predicted at 630 and 584 cm1 and observed at 630 and 591 cm1 for the twist conformer whereas they are predicted at 610 and 600 cm1, respectively, for the axial form. These ab initio predictions support the twist conformer as the second conformer. A more interesting observation is the predicted and observed intensities of the ν30 and ν31 modes of the twist form where they are predicted at 848 and 838 cm1, and observed at exactly these frequencies. Their predicted intensities are 84.0 and 33.5 km mol1 and both with A/C contours with a strong Q-branch (Figure 5). However, for these two modes for the axial conformer they are predicted 848 and 834 cm1 with intensities of 116.5 and 0.9 km mol1, respectively, with A/C contour for the higher frequency band but a B-type contour for this conformer which has a plane of symmetry. Therefore, these vibrational data provide convincing experimental evidence that the less stable conformer is the twist form. The data obtained from the variable low temperature study in xenon has provided some excellent information for determining the correct assignment for the normal modes of both of the conformers. It is clear that the previous proposed assignments are in error in several cases. The lack of agreement between those proposed earlier and the ones obtained from the current studies arises from the much larger amount of data available in this study such as the band contours, the relative predicted intensities from the ab initio calculation, and the intensity changes from the variable temperature study. Additionally, much of the earlier data 7481

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A

Figure 6. Mid-infrared spectra (3001500 cm1) of pyrrolidine: (A) gas; (B) liquid xenon solution; (C) solid.

were obtained from condensed-phase samples rather than for the gas. For the liquid there must be significant association which can significantly shift the band frequencies as well as change the intensities. Association can also be present in the solution but with liquid xenon the 4 cm cell makes it possible to utilize a very dilute solution of approximately 104 M which significantly reduces the association. For the solid the lower frequency region is drastically changed as can be seen in Figure 6 where the pronounced fundamentals in the gas [xenon] at 568 [568], 739 [733], and 792 [791] cm1 either disappear or shift to ∼900 cm1. These spectra dramatically show that an assignment made from spectral data from the condensed phases cannot be compared to one for the gas for the monomer. Therefore, we believe that the current proposed vibrational assignments for the two conformers are based on good scientific evidence for the monomeric species. An attempt was made to obtain the conformational stability by utilizing two bands observed in the far-infrared spectrum of a krypton solution which was initially thought to be due to the two different conformers. However, on determining the enthalpy difference for these two bands a value of 550 cm1 was obtained and, therefore, it was concluded that this enthalpy value was not between the two forms of pyrrolidine but was due to the dimerization arising from the NH association from one molecule to the electron pair on other molecules based on the much lower temperature of the krypton solution. This enthalpy value was certainly too large based on the amount of the second conformer present at ambient temperature as well as the predicted energy difference from the ab initio and density function theory calculations. The enthalpy difference of 109 ( 11 cm1 (1.30 ( 0.13 kJ mol1) has the expected 10% uncertainty which is usually found by using the variable-temperature technique of rare gas solutions. The spread in values is not unusual particularly when the higher frequencies in the fingerprint region have to be used, so there are many opportunities for some interferences for each conformer. Also, even with the very low concentration of the sample one does expect some association as the temperature is lowered which could lead to slight lowering of the experimentally

ARTICLE

determined enthalpy value. However, the predicted value from the ab initio calculations can vary extensively from the different calculations. For example, it has been shown that only theoretical methods that incorporate correlation effect in conjunction with large basis sets provided satisfactory predictions for allyl amine.33 The structural parameters (Table 4) determined by utilizing the previously reported microwave rotational constants along with the predicted parameters from MP2(full)/6-311þG(d,p) calculations have provided excellent distances and angles which we believe to have small uncertainties such as 0.003 Å for the heavy atom distances, 0.002 Å for the CH and NH distances, and 0.5 for the angles. These parameters are probably as good as can be obtained experimentally by any technique for the gas phase. It is interesting to compare the CN distance of pyrrolidine with other secondary amines. For example, the CN distance for dimethylamine has the reported34 rs value of 1.462 ( 0.005 Å and the CNC angle of 112.2 ( 0.2. The r0 parameters for ethylmethylamine35 are 1.462(5) and 1.461(5) Å for the two CN distances and 112.1(5) for the CNC angle which are essentially the same as those for the dimethyl compound. Disregarding the uncertainties, these CN distances are slightly smaller than the value of 1.469(3) Å for pyrrolidine. Of course, the CNC angle is much smaller at 104.1(5), as expected, since the angle in the ring must be much smaller than those in straight-chain molecules. There are no structural parameters for N-methylpyrrolidine obtained from spectroscopic data but they were reported from an electron diffraction study36 where the ra values are 1.455(3) Å for the CN, 1.542(4) Å for the two CC, and 1.555(4) Å for the unique CC distance and the CNC angle is 107.4(17) for the envelope-equatorial form. These electron diffraction parameters provided rotational constants which satisfactorily predicted the rotational constants obtained later from the microwave study.37 Also, these parameters are in good agreement with those obtained in the current study for this conformer of pyrrolidine. It is rather surprising that the second most stable conformer for pyrrolidine is the twisted form rather than the envelope-axial form which has been found for most of the previously studied fivemembered rings which have two stable conformers. Therefore, it would be of interest to investigate whether there are other five membered rings with two conformers where the equatorial form is the more stable conformer and the second most stable form the twisted conformer. For example, N-methylpyrrolidine should have a second conformer but from the several spectroscopic studies there has not been a second conformer reported. Therefore, studies of some additional N-substituted pyrrolidines would be of considerable interest. If there are additional ones with the second conformer the twisted form it would be of interest to learn what the major factors are which contribute to the difference.

’ ASSOCIATED CONTENT

bS

Supporting Information. Table S1, quadratic centrifugal distortion constants (kHz) for isotopomers of twist and equatorial pyrrolidine.Table S2, calculated frequencies (cm1) for the envelopeaxial (Cs) and twist (C1) forms of pyrrolidine. This material is available free of charge via the Internet at http://pubs.acs.org.

’ AUTHOR INFORMATION Corresponding Author

*Phone: 01-816-235-6038. Fax: 01 816-235-2290. E-mail: durigj@ umkc.edu. 7482

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483

The Journal of Physical Chemistry A Present Addresses ‡

Department of Chemistry, Mansoura University, Mansoura, Egypt.

Notes †

Taken in part from the theses of A.M.E. and A.G., which were submitted in partial fulfillments of their Ph.D. degrees, and in part from the thesis of S.S.P., which will be submitted in partial fulfillment of her Ph.D. degree.

’ ACKNOWLEDGMENT J.R.D. acknowledges the University of Missouri-Kansas City for a Faculty Research Grant for partial financial support of this research. ’ REFERENCES (1) Pitzer, K. S.; Donath, W. E. J. Am. Chem. Soc. 1959, 81, 3213– 3218. (2) Laane, J. In Vibrational Spectra and Structure; Durig, J. R., Ed., Marcel Dekker: New York, 1972; Vol. 1, Chapter 2. (3) Harris, D. O.; Engerholm, G. G.; Tolman, C. A.; Luntz, A. C.; Keller, R. A.; Kim, H.; Gwinn, W. D. J. Chem. Phys. 1969, 50, 2438–2445. (4) Durig, J. R.; Wertz, D. W. J. Chem. Phys. 1968, 49, 2118–2121. (5) Badawi, H. M.; Herrebout, W. A.; Zheng, C.; Mohamed, T. A.; van der Veken, B. J.; Durig, J. R. Struct. Chem. 2003, 14, 617–635. (6) Badawi, H. M.; Herrebout, W. A.; Mohamed, T. A.; van der Veken, B. J.; Sullivan, J. F.; Durig, D. T.; Zheng, C.; Kalasinsky, K. S.; Durig, J. R. J. Mol. Struct. 2003, 645, 89–107. (7) Guirgis, G. A.; El Defrawy, A. M.; Gounev, T. K.; Soliman, M. S.; Durig, J. R. J. Mol. Struct. 2007, 832, 73–83. (8) Guirgis, G. A.; El Defrawy, A. M.; Gounev, T. K.; Soliman, M. S.; Durig, J. R. J. Mol. Struct. 2007, 834, 17–29. (9) Ekejiuba, I. O. C.; Hallam, H. E. Spectrochim. Acta 1970, 26A, 59–66. (10) Ekejiuba, I. O. C.; Hallam, H. E. Spectrochim. Acta 1970, 26A, 67–75. (11) Wertz, D. W.; Shasky, W. E. J. Chem. Phys. 1971, 55, 2422. (12) Carreira, L. A. J. Chem. Phys. 1971, 55, 181. (13) Durig, J. R.; El Defrawy, A. M.; Ganguly, A.; Gounev, T. K.; Guirgis, G. A. J. Phys. Chem. A 2009, 113, 9675–9683. (14) Caminati, W.; Oberhammer, H.; Pfafferott, G.; Filgueira, R. R.; Gomez, C. H. J. Mol. Spectrosc. 1984, 106, 217–226. (15) Pfafferott, G.; Oberhammer, H.; Boggs, J. E.; Caminati, W. J. Am. Chem. Soc. 1985, 107, 2305–2309. (16) Caminati, W.; Dell’Erba, A.; Maccaferri, G.; Favero, P. G. J. Mol. Spectrosc. 1998, 191, 45–48. (17) Møller, C.; Plesset, M. S. Phys. Rev. 1934, 46, 618–622. (18) Billes, F.; Geidel, E. Spectrochim. Acta A 1997, 53, 2537–2551. (19) El Gogary, T. M.; Soliman, M. S. Spectrochim. Acta A 2001, 57, 2647–2657. (20) Carballeira, L.; Perez-Juste, I.; Alsenoy, C. V. J. Phys. Chem. A 2002, 106, 3873–3884. (21) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, Jr., J. A.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.;

ARTICLE

Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian 03 (Revision E.01); Gaussian, Inc.: Wallingford, CT, 2004. (22) Pulay, P. Mol. Phys. 1969, 17, 197–204. (23) Guirgis, G. A; Zhu, X.; Yu, Z.; Durig, J. R. J. Phys. Chem. A 2000, 104, 4383–4393. (24) Bulanin, M. O. J. Mol. Struct. 1973, 19, 59–79. (25) van der Veken, B. J.; DeMunck, F. R. J. Chem. Phys. 1992, 97, 3060–3072. (26) Herrebout, W. A.; van der Veken, B. J.; Wang, A.; Durig, J. R. J. Phys. Chem. 1995, 99, 578–585. (27) Bulanin, M. O. J. Mol. Struct. 1995, 347, 73–82. (28) Herrebout, W. A.; van der Veken, B. J. J. Phys. Chem. 1996, 100, 9671–9677. (29) Durig, J. R.; Ng, K. W.; Zheng, C.; Shen, S. Struct. Chem. 2004, 15, 149–157. (30) McKean, D. C. J. Mol. Struct. 1984, 113, 251–266. (31) van der Veken, B. J.; Herrebout, W. A.; Durig, D. T.; Zhao, W.; Durig, J. R. J. Phys. Chem. A 1999, 103, 1976–1985. (32) Velino, B.; Millemaggi, A.; Dell’Erba, A.; Caminati, W. J. Mol. Struct. 2001, 599, 89–93. (33) Galabov, B.; Kim, S.; Xie, Y.; Schaefer, H. F., III; Leininger, M. L.; Durig, J. R. J. Phys. Chem. A 2008, 112, 2120–2124. (34) Wollrab, J. E.; Laurie, V. W. J. Chem. Phys. 1968, 48, 5058–5066. (35) Durig, J. R.; Zheng, C.; Herrebout, W. A.; van der Veken, B. J. J. Mol. Struct. 2002, 641, 207–224. (36) Pfafferott, G.; Oberhammer, H.; Boggs, J. E. J. Am. Chem. Soc. 1985, 107, 2309–2313. (37) Caminati, W.; Scappini, F. J. Mol. Spectrosc. 1986, 117, 184–194.

7483

dx.doi.org/10.1021/jp200692t |J. Phys. Chem. A 2011, 115, 7473–7483