Conformational Stability from Variable Temperature Infrared Spectra of

Mar 2, 2011 - Department of Chemistry, University of Missouri-Kansas City, Kansas City, ... Department of Chemistry and Biochemistry, College of Charl...
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 Ab Initio Calculations of Cyclopropylisocyanate James R. Durig,*,‡ Sarah Xiaohua Zhou,†,‡ Gamil A. Guirgis,§ and Charles J. Wurrey‡ ‡

Department of Chemistry, University of Missouri-Kansas City, Kansas City, Missouri 64110, United States Department of Chemistry and Biochemistry, College of Charleston, Charleston, South Carolina 29424, United States

§

bS Supporting Information ABSTRACT: Infrared spectra (4000 to 400 cm-1) of the gas and variable temperature xenon solutions, and the Raman spectrum of the liquid have been recorded for cyclopropylisocyanate. The enthalpy difference has been determined to be 77 ( 8 cm-1 (0.92 ( 0.10 kJ/mol) with the trans form more stable than the cis conformer with 59 ( 2% present at ambient temperature. By utilizing three rotational constants for each conformer, combined with structural parameters predicted from MP2(full)/6-311þG(d, p) calculations, the adjusted r0 parameters have been obtained. Heavy atom structural parameters for the trans [cis] conformers are the following: distances (Å) (C-C2,3) = 1.509(3) [1.509(3)], (C2-C3) = 1.523(3) [1.521(3)], (C-N) = 1.412(3) [1.411(3)], (NdC) =1.214(3) [1.212(3)], (CdO) = 1.163(3) [1.164(3)]; angles () — CCN = 116.7(5) [120.1(5)], — CNC = 136.3(5) [137.6(5)]. The centrifugal distortion constants have been predicted from ab initio and DFT calculations and are compared to the experimentally determined values.

’ INTRODUCTION Organoisocyanates and organoisothiocyanates provide interesting challenges to structural scientists for the experimental determination of their conformational stabilities and structural parameters because of the large CNC(X) angle, which gives rise to many of the molecules having low-frequency, large amplitude bending vibrations that may be very anharmonic. Additionally, the barrier to internal rotation about the bond to the CNdC(dX) group may be relatively small, which can result in nearly free internal rotation of the NCX moiety. Also, because of the very low frequency of the bending mode there may be a very large number of excited vibrational states populated even at dry ice temperature, which is the temperature usually used to record the normal microwave rotational spectra of such molecules. It is also difficult to obtain good structural parameters for many of these molecules from electron diffraction studies. The nearly free rotation of the NCX moiety for more than 50% of the molecules at ambient temperature, along with a large number of anharmonic bending vibrations of very low frequency in many excited states, result in a large number of undetermined parameters, making it difficult to determine conformational stabilities. Further confusion about the structure of several of these molecules has arisen from ab initio predictions of the structures and conformational stabilities of organo-NCX compounds from calculations with relatively small basis sets which often do not agree with the most stable conformer or the stable conformers determined experimentally. Also, in some cases the potential r 2011 American Chemical Society

wells are so shallow that they cannot accommodate a vibrational state, hence, the microwave determined rotational constants will be different from the conformation which has the re minimum. Therefore, we have been reinvestigating the conformational structural stabilities of a number of these molecules, particularly by recording variable temperature vibrational spectra of the gas or of liquid rare gas solutions. Additionally, we have obtained the heavy atom structural parameters as well as ab initio predictions of the conformational stabilities by utilizing much larger basis sets than those that were previously used and which predicted the most stable conformer incorrectly. We initiated these recent studies by investigations of methylisothiocyanate (CH3NCS)1 and methylisocyanate (CH3NCO),2 where free internal rotation was observed for the methyl rotors for both molecules and the r0 structural parameters were determined. These studies have been followed with spectroscopic and theoretical studies of ethylisothiocyanate (CH3CH2NCS)3 and ethylisocyanate (CH3CH2NCO).4 For the thio- compound, it was concluded from the microwave spectra5 that the most stable conformer was the trans form, whereas from ab initio calculations6 at the MP2 level with various basis sets up to TZVP, it was concluded that the gauche form was the most stable form. However, from a more recent spectroscopic and theoretical study7 using larger basis sets it was conclusively shown that the most stable form was the cis conformer. A Received: September 17, 2010 Revised: January 13, 2011 Published: March 02, 2011 2297

dx.doi.org/10.1021/jp108879u | J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A similar study of ethylisocyanate with large basis sets predicted the depth of the cis well to range from 11 cm-1 to a high value of 31 cm-1 and with the largest basis set (425 for the MP2(full)ccPVQZ calculation) predicting a value of 19 cm-1. This is a strong indication that the gauche wells could not contain a bound vibrational energy level. Also, the spectroscopic data was only consistent with the most stable conformer being the cis form. As a continuation of these studies of pseudohalogen organic compounds we have again investigated the conformational stability and determined the structural parameters of cyclopropylisocyanate. In an earlier microwave, infrared and Raman spectral study8 it was concluded that there were two conformers present at ambient temperature, which were the most stable trans form and the cis form (Figure 1). However, from a later theoretical study,9 the authors claimed that the trans and cis forms previously reported are actually two gauche forms with dihedral [CNC(1)C(2)] angles of 144.9 and 324.3 based on MP2/631G(d) calculations. There was little information on how the calculations were performed and a figure in the article implied that this dihedral angle should be 180 for the cis form. There was also a later microwave study.10 Since in the earlier study,8 the value of the A rotational constants for both conformers had significant uncertainties, a much larger number of transitions were assigned and measured for both conformers. For the cis conformer, it was not possible to fit the five usually determined centrifugal distortion constants by the standard equations, which lead the authors to question the conformation for the form assigned as the cis conformer. Thus, there is still a question concerning the stable conformers for cyclopropylisocyanate. Additionally, in the earlier study8 the enthalpy difference was determined for the liquid and only estimates were made for the r(C-N), — CNC, and — CCN parameters. Therefore, we have carried out a variable temperature infrared study of a xenon solution of cyclopropylisocyanate which should give an enthalpy value near that of the gas. Also, with three well-determined rotational constants it should be possible to obtain reliable structural parameters, and with ab initio calculations using much larger basis sets the predicted conformational stability is expected to be more reliable. The results of this spectroscopic and theoretical study are reported herein.

’ EXPERIMENT AND THEORETICAL CALCULATIONS The sample of cyclopropylisocyanate was prepared via a twostep procedure according to the method of Kricheldorf and Regel.11,12 In the first step, an intermediate product of diphenyl diazido silane was formed from dichloro diphenyl silane and sodium azide. The final product is achieved through the Curtius rearrangement of the intermediate formed from reaction with

ARTICLE

cyclopropyl carbonyl chloride. The sample was purified by trapto-trap distillation and the final purification was obtained by using a low-temperature, low-pressure sublimation column. The sample identity was verified by NMR and infrared spectroscopic data. The mid-infrared spectra from 3500 to 400 cm-1 of the gas shown in Figure 2A was recorded on a Perkin-Elmer model 2000 Fourier transform spectrometer equipped with a nichrome wire source, Ge/CsI beamsplitter and DTGS detector. The spectrum of the gas was obtained with the samples contained in 12 cm cells equipped with CsI windows. Atmospheric water vapor was removed from the spectrometer chamber by purging with dry nitrogen. Interferograms obtained after 128 scans for the gas sample and the background reference were transformed by using a boxcar apodization function with theoretical resolutions of 0.5 cm-1 for the gaseous sample. The mid-infrared spectra of cyclopropylisocyanate dissolved in liquefied xenon (Figure 2B) were recorded on a Bruker model IFS-66 Fourier interferometer equipped with a Globar source, Ge/KBr beamsplitter, and DTGS detector. The interferograms were recorded at variable temperatures ranging from -55 to 100 C with 100 scans and transformed by a Blackman-Harris apodization function with a theoretical resolution of 1.0 cm-1. The temperature studies in liquefied xenon were carried out in a specially designed cryostat cell, which is composed of a copper cell with a 4 cm path length and wedged silicon windows sealed to the cell with indium gaskets. The temperature was monitored by two platinum thermoresistors and the cell was cooled by the vapors from boiling liquid nitrogen. All of the observed fundamental

Figure 2. Infrared spectra of cyclopropylisocyanate: (A) gas; (B) xenon solution at -100 C. The asterisk indicates the spectral region at ∼600 cm-1 has nearly zero energy due to the absorption from the silicon windows.

Figure 1. Possible conformations of cyclopropylisocyanate, δ indicates the dihedral angle CdN-C-H4. 2298

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

modes for the trans and cis conformers in the infrared and Raman spectra are listed in Tables 1 and 2, respectively. The Raman spectra of the liquid (Figure 3) was recorded on a Spex model 1403 spectrophotometer equipped with a SpectraPhysics model 2017 argon ion laser operating on the 514.5 nm line. The laser power used was 0.5 W with a spectral bandpass of 3 cm-1. The spectrum of the liquid was recorded with the sample sealed in a Pyrex glass capillary. The Raman spectrum of the gas (Figure 4A) at ambient vapor pressure was recorded by using the standard Cary multipass cell with 1-W laser power at the sample. The LCAO-MO-SCF calculations were performed with the Gaussian 03 program13 by using Gaussian-type basis functions. The energy minima with respect to nuclear coordinates were obtained by simultaneous relaxation of all geometric parameters consistent with symmetry restrictions using the gradient method of Pulay.14 A number of basis sets starting from 6-31G(d), and increasing to 6-311þG(3df,3pd), were employed at the level of Møller-Plesset perturbation theory15 to the second order (MP2), as well as hybrid density functional theory by the B3LYP method, to obtain energy differences (Table 3) among the three most likely conformers of cyclopropylisocyanate. At all levels of calculation carried out, with and without diffuse functions, predicted conformational stabilities varied extensively with the size of the basis set.

To aid in making the vibrational assignment (Tables 1 and 2), we have carried out a normal coordinate analysis by utilizing the force fields obtained from the Gaussian 03 program at the MP2(full)/6-31G(d) level. The internal coordinates used to calculate the G and B matrices for cyclopropylisocyanate are listed along with the structural parameters in Table 5. By using the B matrix, the force field in Cartesian coordinates was converted to a force field in internal coordinates.16 Subsequently, scaling factors of 0.88 for the CH stretches and 0.90 for all other modes were used, along with the geometric average of scaling factors for interaction force constants, to obtain the fixed scaled force fields and the resultant wavenumbers (Table 1). A set of symmetry coordinates (Table 1S) was used to determine the corresponding potential energy distributions (P.E.D.s). The observed and calculated wavenumbers of cyclopropylisocyanate along with the calculated infrared intensities, Raman activities, depolarization ratios, and P.E.D.s are given in Tables 1 and 2. To identify the fundamental vibrations for possible conformers of cyclopropylisocyanate, we have simulated the Raman spectrum (Figures 3 and 4) from the scaled ab initio MP2(full)/ 6-31G(d) results. The evaluation of Raman activity by using the analytical gradient method has been developed.17,18 The activity Sj can be expressed as: Sj = gj (45R2j þ 7β2j ), where gj is the degeneracy of the vibrational mode j, Rj is the derivative of the

Table 1. Calculated and Observed Vibrational Frequencies (cm-1) of Cyclopropylisocyanate for the trans Conformer infrared vib. No. 0

A

A00

approx. description

MP2/6-31G(d) MP2 scaleda IR int.b Raman act.c dp ratio gas

Xe

Raman

contour

gasd liquid

PEDe

A

C

ν1

CH2 antisym stretch

3320

3114

5.2

34.9

0.73

3109 3097 3107 3095

99S1

43

57

ν2

CH stretch

3229

3029

11.0

125.6

0.12

3036 3034 3034 3025

89S2, 10S3

35

65

ν3

CH2 symmetric stretch

3218

3019

0.9

98.3

0.06

88S3, 11S2

59

41

ν4

NCO antisym stretch

2410

2286

938.4

3.7

0.16

2281 2277 2289 2281

98S4

99

0

ν5

CH2 deformation

1568

1488

32.8

29.2

0.34

1480 1476 1482 1472

62S5, 18S8

100 0

ν6

NCO symmetric stretch

1515

1437

18.1

56.1

0.23

1446 1448 1449 1445

39S6, 26S5, 25S13

80

20

ν7

CH bend

1402

1330

53.0

8.3

0.74

1344 1342

99

1

ν8 ν9

ring breathing CH2 twist

1272 1176

1207 1115

0.3 4.2

16.9 4.1

0.15 0.62

1205 1201 1201 1200 1116 1115

65S8 38S9, 21S12, 16S7

96 30

4 69

ν10

CH2 wag

1111

1054

11.6

0.3

0.75

1032 1026 1034

92S10

72

28

ν11

ring deformation

979

929

7.4

9.1

0.60

918

67S11, 13S13

10

90

ν12

CH2 rock

826

784

1.4

3.0

0.28

∼800 ∼800 797

67S12, 13S13

0

100

ν13

ringC-N stretch

761

724

7.9

11.7

0.42

733

734 730

23S13, 34S9, 16S11

2

98

ν14

NCO bend

616

614

23.1

2.0

0.65

630

635 632

82S14

52

48

ν15

CC-N bend

371

371

12.4

3.5

0.31

369

367 376

77S15

93

7

ν16 ν17

CNC bend CH2 antisym stretch

130 3309

130 3104

6.7 0.4

2.6 76.5

0.72 0.75

126d

100 0

3095

88S16 99S17

ν18

CH2 symmetric stretch

3215

3016

6.6

27.6

0.75

3021

100S18

ν19

CH2 deformation

1515

1437

3.2

10.0

0.75

1445

100S19

ν20

CH2 twist

1236

1173

0.2

7.8

0.75

1167

ν21

CH bend

1158

1099

1.0

1.9

0.75

1091

ν22

CH2 wag

1110

1053

3.6

0.2

0.75

1050 1045

ν23

ring deformation

990

939

11.4

11.7

0.75

936

∼930

932

64S23, 14S26, 12S24

ν24 ν25

CH2 rock NCO bend

854 543

810 543

8.7 15.6

5.2 0.7

0.75 0.75

811 566

810 566

819 570

32S24, 30S23, 20S21, 16S20 100S25

ν26

R-NCO bend

413

413

2.9

0.4

0.75

416

418

ν27

ringC-N bend

12

12

0.5

3.0

0.75

3018

915 735

∼1350 64S7, 20S6

920 913

1172

52S20, 42S24 69S21, 22S20 94S22

87S26, 10S24 100S27

a Force constant scaling factors: 0.88 for CH stretches; 0.9 for all other modes. b Infrared intensities in km/mol. c Raman activities in Å4/u. d Ref 8. e PEDs from MP2 scaled calculation: values less than 10% are omitted.

2299

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

Table 2. Calculated and Observed Vibrational Frequencies (cm-1) of Cyclopropylisocyanate for the cis Conformer infrared

0

A

A00

vib. No.

approx. description

MP2/6-31G(d) MP2 scaleda IR int.b Raman act.c dp ratio gas

Raman

contour

Xe gasd liquid

PEDe

A

B

ν1

CH2 antisymmetric stretch

3313

3108

7.3

32.1

0.50

3102 3097 3107 3095 99S1

89

11

ν2

CH stretch

3259

3057

9.7

100.2

0.26

3053 3047 3054 3053 98S2

59

41

ν3

CH2 symmetric stretch

3214

3015

0.9

135.4

0.06

99S3

57

43

ν4

NCO antisymmetric stretch

2406

2282

797.6

2.2

0.09

98S4

88

12

ν5

CH2 deformation

1564

1484

9.2

17.3

0.52

1465 1472 1465 1470 77S5, 15S8

23

77

ν6

NCO symmetric stretch

1515

1437

15.7

48.5

0.17

1446 1448

1445 53S6, 33S13

93

7

ν7

CH bend

1439

1366

9.8

7.4

0.53

79S7

97

3

ν8 ν9

ring breathing CH2 twist

1264 1154

1199 1095

3.4 27.3

16.0 9.9

0.19 0.75

1197 1198 1106 1101

42 97

58 3

ν10

CH2 wag

1107

1050

19.5

0.3

0.73

1048 1045

ν11

ring deformation

971

920

0.9

11.9

0.51

ν12

CH2 rock

822

780

8.7

1.6

0.24

ν13

ringC-N stretch

763

725

14.7

9.7

ν14

NCO bend

624

620

24.1

ν15

CC-N bend

424

424

ν16 ν17

CNC bend CH2 antisymmetric stretch

107 3302

ν18

CH2 symmetric stretch

ν19

CH2 deformation

ν20

CH2 twist

1235

1172

0.5

7.9

0.75

1175 1167

1172 45S20, 44S24, 10S21

ν21

CH bend

1169

1109

0.7

1.2

0.75

1112 1112

68S21, 27S20

ν22

CH2 wag

1108

1051

3.7

0.1

0.75

1031 1026

94S22

ν23

ring deformation

991

940

10.7

8.2

0.75

930 928

63S23, 14S26, 12S24

ν24 ν25

CH2 rock NCO bend

856 551

812 551

8.6 18.2

5.9 0.5

0.75 0.75

805 805 571 566

33S24, 31S23, 17S20, 16S21 100S25

ν26

R-NCO bend

415

315

2.4

0.6

0.75

408 405

88S26

ν27

ringC-N bend

46

46

1.1

2.4

0.75

∼40d

100S27

1200 63S8, 11S9 1105 36S9, 22S12, 14S7, 14S6 94S10

9

91

65S11, 17S13, 11S6

0

100

788 787 792 792

66S12, 12S13, 10S9

81

19

0.32

725 726

26S13, 26S9, 23S11

33

67

1.0

0.74

615 613 623

78S14, 11S16

31

69

15.6

1.4

0.49

72S15

64

36

107 3098

3.6 0.5

3.2 73.2

0.58 0.75

88S16, 11S15 3095 100S17

92

8

3210

3011

6.8

21.0

0.75

1520

1442

3.5

8.2

0.75

920 726

107

100S18 1448

1445 100S19

a Force constant scaling factors: 0.88 for CH stretches; 0.9 for all other modes. b Infrared intensities in km/mol. c Raman activities in Å4/u. d Ref 8. e PEDs from MP2 scaled calculation: values less than 10% are omitted.

isotropic polarizability, and βj is that of the anisotropic polarizability. The Raman scattering cross sections, ∂σj/∂Ω, which are proportional to Raman activities, can be calculated from the scattering activities as well as the predicted wavenumbers for each normal mode, by using the relationship:19,20 ∂σj/∂Ω = [(2π)4/45] [(ν0 - νj)4/(1 - exp(-hcνj/kT))] [h/(8π2cνj)] Sj, where ν0 is the excitation wavenumber, νj is the vibrational wavenumber of the jth normal mode, and Sj is the corresponding Raman scattering activity. To obtain the polarized Raman scattering cross sections, the polarizabilities are incorporated into Sj by multiplying Sj by (1 - Fj)/(1 þ Fj), where Fj is the depolarization ratio of the jth normal mode. The Raman scattering cross sections and calculated wavenumbers obtained from the scaled ab initio force fields were used together with a Lorentzian function to obtain the simulated Raman spectra. To further support the vibrational assignments, the infrared spectrum (Figure 5) was predicted by using scaled wavenumbers from MP2(full)/6-31G(d) results. Infrared intensities were calculated based on the dipole moment derivatives with respect to Cartesian coordinates. The derivatives were transformed into normal coordinate derivatives by (∂μu/∂Q i) = Σ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/∂Q i)2 þ (∂μy/∂Q i)2 þ (∂μz/∂Q i)2]. Vibrational Assignment. To determine the enthalpy difference between the two conformers it is necessary to provide vibrational assignments for the fundamentals of each conformer in the region where the conformer pairs will be chosen for the enthalpy determination. Because the lower frequency region will have the fewest number of possible overtones or combination bands, the assignments for the fundamentals in this region will be very important to make before selecting pairs of conformer bands. However, because of the large amplitude motion of the NCO and R-NCO bends many of the bands below 700 cm-1 are very broad in the xenon solution, which makes it difficult to confidently assign them to the correct conformer. However, in the region between 700 to 1400 cm-1 most of the fundamentals can be indentified for each conformer; but in several cases the similar vibrations for the two conformers are observed within a very few wavenumbers from each other as predicted from the ab initio calculations. To make the assignments, considerable reliance was placed on the predicted fundamental frequencies and intensities from the ab initio calculations as well as the infrared band contours in the gas. From these experimental data most of the fundamentals in the 700 to 1400 cm-1 could be assigned with confidence. 2300

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

Table 3. Calculated Electronic Energiesa (hartree) and Energy Differencesb (cm-1) for the trans, cis, and the skew Conformers of Cyclopropylisocyanate number of basis functions

trans

cis

skewc

MP2(full)/6-31G(d)

100

0.531149

52

105

MP2(full)/cc-PVDZ

109

0.569719

84

135

MP2(full)/6-31þG(d)

124

0.548210

95

173

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

115 139

0.572075 0.588575

57 97

112 184

MP2(full)/TZVP

144

0.778199

14

110

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

138

0.767042

36

90

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

162

0.776779

7

124

MP2(full)/aug-cc-PVDZ

183

0.623931

19

84

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

183

0.839084

75

122

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

207

0.847659

54

121

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

250 274

0.943928 0.951234

39 16

97 96

method/basis set

Figure 3. Observed (room temperature) and predicted (MP2(full)/ 6-31G(d) at 25 C) Raman spectra of cyclopropylisocyanate: (A) liquid; (B) predicted spectrum of the mixture of cis and trans conformers with ΔH = 77 cm-1; (C) predicted spectrum of the pure cis conformer; (D) predicted spectrum of the pure trans conformer; (E) solid. Asterisk indicates that the intensity has been reduced by half.

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

295

0.968773

72

129

MP2(full)/aug-cc-PVTZ

391

1.960505

53

133

B3LYP/6-31G(d)

100

1.372666

203

237

B3LYP/6-31þG(d)

124

1.384722

250

284

B3LYP/6-311G(d,p)

138

1.453174

257

279

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

162

1.458582

232

264

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

183 207

1.461373 1.466843

271 226

285 250

B3LYP/6-311G(2df,2pd)

250

1.472328

243

268

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

274

1.477358

209

241

B3LYP/6-311G(3df,3pd)

295

1.476357

256

285

B3LYP/6-311þG(3df,3pd)

319

1.480813

198

233

a Energy of trans conformer is given as -(E þ 284) hartree. b All energy differences are relative to the energy of the trans conformer. c Transition state.

Figure 4. Observed (room temperature) and predicted (MP2(full)/ 6-31G(d) at 25 C) Raman spectra of cyclopropylisocyanate: (A) gas; (B) predicted spectrum of the mixture of cis and trans conformers with ΔH = 77 cm-1; (C) predicted spectrum of the pure cis conformer; and (D) predicted spectrum of the pure trans conformer.

The assignments for the fundamentals for the trans form are shown in Table 1 and for the cis conformer in Table 2. For some of the fundamentals with B-type band contours it was very difficult to determine the band center, but for some cases the infrared spectrum of the xenon solution could aid significantly in the band center determination. Nevertheless, there is a broad

nondescript band centered at 1106 cm-1 where two fundamentals for each conformer are predicted. The assignments of the fundamentals are somewhat arbitrary because there is little band contour or intensity information to aid in the identification of the band centers. A similar problem arises in the region 1130 to 1220 cm-1, which makes it difficult to confidently assign the four fundamentals predicted in this range at 1207 (ν8, trans), 1173 (ν20, trans), 1199 (ν8, cis), and 1172 (ν20, cis). All except the 1199 cm-1 band are predicted to have very low intensities and the band contours in this region are very ill defined. Therefore, the infrared spectrum of the xenon solution and Raman spectrum of the liquid were needed to determine the band centers. Only two bands (1198 and 1167 cm-1) can be readily observed in the xenon solution while only one very strong band (1200 cm-1) is observed in the Raman spectrum of the liquid. The band in the xenon solution at 1198 cm-1 is assigned as the ν8 (cis) mode and the 1167 cm-1 bands is assigned as the ν20 (cis) mode. With these assignments, corresponding band centers were identified in the infrared spectrum of the gas at 1197 and 1175 cm-1, respectively. The corresponding Raman band at 1200 cm-1 in the spectrum of the liquid supports these assignments and the weaker Raman line at 1172 cm-1 is consistent with the lower frequency assignment. Enthalpy Difference. The initial ΔH determination8 was made from a single conformer pair from Raman data of the liquid. Unfortunately, the band for the cis conformer was due to stray light so its intensity did not decrease when the temperature 2301

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

Figure 5. Observed and predicted (MP2(full)/6-31G(d)) infrared spectra of cyclopropylisocyanate: (A) xenon solution at -100 C; (B) predicted spectrum of the mixture of cis and trans conformers with ΔH = 77 cm-1; (C) predicted spectrum of the pure cis conformer; and (D) predicted spectrum of the pure trans conformer.

was lowered. Since it did not decrease, the ΔH was determined from essentially the increase of the fundamental due to the trans form when the temperature was lowered so the reported value is expected to be slightly smaller than the actual value. Therefore we have determined the enthalpy difference by a variable temperature study of the xenon solution. For the best results the lowest frequency fundamentals should be used for the conformer pairs since they have the lowest probability of having intensity contributions from overtones or combination bands. For the cis conformer the bands at 726 and 787 cm-1 have been confidently assigned to this form and for the trans form the fundamentals at 810, 915, 928, and 1026 cm-1 were utilized (Figure 2B). These bands were selected because they are well resolved and well separated, with relatively flat baselines, and they are without predicted underlying bands from the other conformer. From these fundamentals for each conformer, eight conformer pairs were selected for the conformational stability determination. The intensities of the infrared bands were measured as a function of temperature and their ratios were determined. By application of the van’t Hoff equation, -ln K = ΔH/(RT) - ΔS/R, ΔH was determined from a plot of -ln K versus 1/T (Figure 1S), where ΔH/R is the slope of the line and K is substituted with the appropriate intensity ratios, that is, Itrans/Icis. It is assumed that ΔH and ΔS are not functions of temperature in the temperature range studied. The resulting values with statistical uncertainties of one sigma for each set are listed in Table 4. The average of the eight values is 77 ( 3 cm-1 (0.92 ( 0.04 kJ/mol)) where the error limit is derived from the statistical standard deviation of one sigma of the measured intensity data, where the data are taken as a single set. Although the statistical uncertainties are relatively small, they do not take into account possible contribution from combination or overtone bands from other conformers contributing to the

ARTICLE

measured fundamental band intensities. The variations of ΔH values are undoubtedly due to these types of interferences but by taking eight conformer pairs it is expected that these effects are nearly canceled. However, the nature of the technique utilized for enthalpy determination is not expected to provide data better than 10%. Therefore, a more realistic ΔH value is 77 ( 8 cm-1 (0.92 ( 0.10 kJ/mol)). From this ΔH value, a realistic estimation of the abundance of the less stable cis conformer present at ambient temperature is 41 ( 2%. Structural Parameters. In the initial microwave study8 only the r(C1-N), — C2,3C1N, and — C-NdC were determined for both conformers from the three rotational constants for each form with the other parameters estimated from related molecules. The larger C-N distance and the smaller CCN angle of the trans conformer compared to the corresponding parameters of the cis form are consistent with the expected difference of these parameters. However, the two A rotational constants had rather large uncertainties which could have affected the values of the three determined parameters. Additionally, the values estimated from those obtained from other molecules may not be comparable to the values of similar parameters for cyclopropylisocyanate. Furthermore, it is suggested that by combining ab initio predicted values as constraints more reliable structural parameters could probably be obtained. We have found that good structural parameters for hydrocarbons and many substituted ones can be determined by adjusting the structural parameters obtained from the ab initio MP2/6-311þG(d,p) calculations to fit the rotational constants obtained from microwave experimental data by using a computer program “A&M” (ab initio and microwave) developed21 in our laboratory. 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 four heavy atom distances for cyclopropylisocyanate, bond angles and torsional angles in the same set keep their difference in degrees. This assumption is based on the fact that the errors from ab initio calculations are systematic. Additionally, we have also shown that the difference in predicted distances and angles from the ab initio calculations for different conformers of the same molecule can usually be used as one parameter, namely, the ab initio predicted difference except for some dihedral angles. However, for the trans and cis conformers, these dihedral angles are zero and 180 so they are not variable parameters. Thus, there are six heavy atom parameters to be adjusted, so it should be possible to obtain “adjusted r0” structural parameters for cyclopropylisocyanate by utilizing the six determined microwave rotational constants from the two conformers since the sets are the same for the two forms. We22 have also shown that ab initio MP2/6-311þG(d,p) calculations predict the r0 structural parameters for more than 50 carbon-hydrogen distances to better than 0.002 Å compared to the experimentally determined values from isolated CH stretching frequencies, which were compared to previously determined23 values from earlier microwave studies. Therefore, all of the carbon-hydrogen distances were taken from the MP2/ 6-311þG(d,p) predicted values for cyclopropylisocyanate. Therefore, a complete structure has been determined for each of the two conformers. The determined adjusted r0 parameters are listed in Table 5 and the final fit of the rotational constants is shown in Table 6. The differences are quite small with all six rotational constants fit to better than 1.0 MHz. The bond distances that changed the most during the fitting are r(C2-C3), r(C1-N), and r(CdO), 2302

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

Table 4. Temperature and Intensity Ratios of the Conformer Bands of Cyclopropylisocyanate from the Infrared Spectra of the Liquid Xenon Solution T (C)

1/T (  10-3 K-1)

I810/I726

I810/I787

I928/I726

I928/I787

I1026/I726

I1026/I787

I1342/I726

I1342/I787

liquid

-60.0

4.6915

2.6141

12.7273

1.9835

4.9618

8.0745

36.1111

5.3114

23.3764

xenon

-65.0 -70.0

4.8042 4.9225

2.6281 2.6721

13.1950 13.4683

2.0410 2.0661

4.7980 4.9505

8.0116 8.0674

37.5000 38.5366

5.3696 5.3331

23.8948 23.9598

-75.0

5.0467

2.6770

13.1596

2.0891

5.0292

8.2652

39.0986

5.4850

24.0773

-80.0

5.1773

2.7893

13.6409

2.0722

5.1135

8.5417

38.8182

5.5966

24.2691

-85.0

5.3149

2.8239

13.9234

2.0954

5.2184

8.7174

39.8649

5.6052

25.1982

-90.0

5.4600

2.8493

14.0000

2.1005

5.2857

8.6462

40.7333

5.6632

25.6448

-95.0

5.6132

2.9066

14.5783

2.1923

5.3826

8.8091

41.7391

5.7683

26.3534

-100.0

5.7753

2.9042

14.6962

2.2556

5.5291

8.8769

41.7722

5.8333

26.6806

77 ( 7

87 ( 9

66 ( 10

82 ( 9

73 ( 9

88 ( 9

63 ( 5

87 ( 7

ΔHa

a

Average value obtained by utilizing all data as a single set gives ΔH = 77 ( 3 cm-1 (0.92 ( 0.03 kJ/mol) with the trans conformer more stable.

Table 5. Structural Parameters (Å and degree), Rotational Constants (MHz), and Dipole Moments (debye) for the trans and cis Conformers, Respectively, of Cyclopropylisocyanate from Ab Initio Calculations (6-311þG(d,p) Basis Set) and Experimental Data trans parametera

int. coord.

MP2

B3LYP

r(C1-C2,3)

R1, R2

1.504

1.507

r(C2-C3)

R3

1.512

1.510

r(C1-N)

R4

1.425

r(NdC)

R5

r(CdO) r(C1-H4) r(C2,3-H8,10)

MWc

cis gauche-1b

gauche-2b

trans

adjusted r0d trans

cis

1.520

1.509(3)

1.509(3)

1.515

1.523(3)

1.521(3)

1.417(10)

1.407(14)

1.412(3)

1.411(3)

1.216

1.210

1.210

1.214(3)

1.212(3)

1.175

1.191

1.170

1.170

1.163(3)

1.164(3)

1.082 1.083

1.079 1.079

1.079 1.086

1.079 1.086

1.085(2) 1.083(2)

1.082(2) 1.084(2)

1.083

1.083

1.079

1.082

1.082

60.4

60.2

60.2

MP2

B3LYP

1.499

1.505

1.508

1.501

1.520

1.504

1.516

1.513

1.506

1.515

1.425

1.421

1.423

1.421

1.419

1.216

1.200

1.215

1.215

1.198

R6

1.179

1.174

1.190

1.180

r5 r1, r3

1.085 1.083

1.084 1.082

1.081 1.078/1.079

1.082 1.084

r(C2,3-H9,11)

r2, r4

1.083

1.083

1.079

— C2C1C3

π1

60.4

60.1

60.2

cis

1.083(2)

1.083(2)

60.6(5)

60.5(5)

— C2,3C1N

ε1, ε2

117.2

118.1

117.4

120.2

121.1

120.6

117.9(11)

120.3(24)

116.7(5)

120.1(5)

— C-NdC

ψ

134.3

139.0

137.4

135.3

142.1

137.9

136.9(38)

138.6(13)

136.3(5)

137.6(5)

172.7

172.6

172.6

172.2(5)

173.0(5)

115.4(5)

112.0(5)

117.9(5) 117.2(5)

117.7(5) 117.3(5)

— NdC = O

θ

172.9

174.1

172.7

172.9

174.2

— NC1H4

ζ

115.4

114.5

115.2

112.0

111.5

— C2,3C1H — C2,3C3,2H10,8

σ1, σ2 β3, β1

117.9 117.2

117.7 117.6

117.9 118.0/117.5

117.6 117.2

117.1 117.5

— C2,3C3,2H11,9

β4, β2

118.4

118.9

118.3/118.9

118.5

119.0

118.5(5)

118.4(5)

— C1C2,3H8,10

R1, R3

116.2

116.7

116.4/116.3

116.6

117.2

116.4/116.7

116.2(5)

116.7(5)

— C1C2,3H9,11

R2, R4

117.6

117.6

117.8/117.8

117.9

117.8

118.0/118.0

117.6(5)

117.8(5)

CdN-C1-C2

145.6

145.4

144.9

324.4

324.1

324.3

145.6(5)

324.4(5)

A

17016

17190

17218

9997

10706

10346

16691(14)

10161(13)

16941.7

10214.9

B

1769

1758

1753

2208

2085

2139

1784.303(3)

2186.797(10)

1783.4

2186.6

C |μa|

1706 2.997

1693 2.966

1691

2114 3.288

2029 3.217

2063

1716.133(3) 2.56(2)

2106.242(10) 2.720(4)

1717.0

2106.5

0.079

0.076

|μc|

1.061

0.867

|μtot|

3.179

3.090

|μb|

117.5

116.1

116.1

0.17(0.01) 0.71(3)

2.91

3.289

3.218

2.99

2.65(2)

2.726(1)

a

cis and trans refers to the orientation of the hydrogen relative H4. b Ref 9. Optimized geometry determined from MP2/6-31G(d). c Structural parameter from ref 8. d Rotational constants utilized for the fitting are from ref 10.

these adjusted r0 parameters differ by at least 0.010 Å from their predicted values for both of the conformers. The largest adjustment for the bond angle is associated with the C-NdC angle, where the adjusted r0 is at least 2 larger than the MP2(full)/ 6-311þG(d, p) predicted value. The predicted structural parameters from the MP2(full)6-311þG(d,p) calculations for both the trans and cis conformers

are listed in Table 5 along with those from B3LYP/6-311þG(d,p) calculations. These are compared to those previously reported for the two conformers indicated as gauche-1 and gauche-2, which were obtained9 from ab initio predictions from MP2/6-31G(d,p) calculations. By using the dihedral CdN-C-C2 angle as an indication of the conformation, the gauche-1 form reported differs by only 0.7 from this angle for the trans form and only 2303

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

Table 6. Comparison of Rotational Constants (MHz) Obtained from Modified Ab Initio MP2(full)/6-311þG(d,p) Predictions, Microwave Spectra, and the Adjusted Structural Parameters for Cyclopropylisocyanate conformer

rotational constant

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

microwavea

adjusted r0

|Δ|

trans

A

17018.0

16941.889(8)

16941.4

0.5

B C

1769.3 1706.2

1784.3121(6) 1716.1278(6)

1783.4 1717.0

0.9 0.9

A

10006.4

10215(6)

10214.9

0.1

B

2206.2

2186.859(6)

2186.6

0.3

C

2113.2

2106.172(6)

2106.5

0.3

cis

a

Rotational constants obtained from a fourth order centrifugal distortion fit, ref 10.

0.1 difference for the assigned gauche-2 form for this angle for the cis conformer. Further clarification concerning the incorrect dihedral angle that was used for determining the gauche forms versus cis or trans will be provided in the discussion

’ DISCUSSION The vibrational assignments given earlier8 were primarily based on infrared band contours, Raman depolarization data, and group frequencies. For those presented in the current study there is significantly more information, which includes the predicted frequencies for the fundamentals, the infrared intensities, and Raman activities, with probably the most important data being the infrared spectra of variable temperature xenon solution from which the bands for the two different conformers could be identified. Therefore, a few of the previously assigned bands as fundamentals, such as the 1005 cm-1 band due to stray light, as well as the bands at 823 and 351 cm-1, have been reassigned. There are some differences because the band centers have different frequencies in some cases. A major problem arises due to the significant breadth of several of the fundamentals of the trans conformer, particularly those arising mainly from motions of the NCO group. Most of the fundamentals for the cis conformer are very broad, which is probably due to only two energy levels in the potential well for this conformer since transitions for only one excited state were observed in the microwave spectrum. For the trans conformer four excited states of the low frequency NCO bend (torsion) were obtained. Therefore many of the fundamentals for the cis conformer have bands with nondescript contours. The ab initio predicted frequencies for the fundamentals of the trans form with two scaling factors of 0.88 for the CH stretches, 0.90 for the CH bends, and heavy atom stretches and no scaling for the other modes agreed with the modes of A0 symmetry with an average frequency difference of 8.0 cm-1, which is 0.6% error. For the A00 modes, the average difference is 9.1 cm-1 (0.8%), which indicates the value of using the relatively small basis set for aiding the assignment. The descriptions for the vibrational motions usually have only one or two symmetry coordinate contributions with four modes having contributions from three symmetry coordinates and only one (ν24) with four contributors with the largest one contributing only 32%. This latter mode has nearly equal contributions from the CH2 rock and the ring deformation. Also, it should be noted that ν13 is described as the C-N stretch but only 23% of S23 contributes to the 730 cm-1 band assigned to this motion, which also has a significant contribution from the NCO symmetric stretch, and lesser amounts from ν11 and ν12. Nevertheless, for most of the vibrations the approximate description gives the major contribution to that vibration.

The enthalpy value obtained from the xenon solution is expected to be close to the value for the gas since only small interactions are expected to occur between the dissolved molecules and the surrounding noble gas atoms.24-27 Further, the “pseudo gas phase” spectrum shows only small frequency shifts compared to the spectrum of the gas. A significant advantage of this type of cryogenic spectroscopic study is that the conformer bands are better resolved in comparison with those in the spectrum of the gas. This is particularly important since most of the conformer bands for this molecule are expected to be observed within a few wavenumbers of each other. The different conformer bands which can be clearly identified in the spectral region between 1100 and 600 cm-1 of the xenon solution are shown in Figure 1 and the ones indicated are due to the individual conformers. The extremely small uncertainty in ΔH of (3 cm-1 (0.04 kJ/mol) is probably smaller than the technique justifies but it is the statistical value. Nevertheless, the determined value is certainly more accurate than the value that could have been obtained for the gas. There is support for the very small enthalpy value from the ab initio predicted values where the MP2(full)/aug-cc-PVTZ calculation with the largest basis set of 391 wave functions predicts an energy difference of only 20 cm-1. It is expected that this estimated value is the most reliable one because of the large basis set.28 Although large basis sets are necessary for the determination of the energy difference between conformers, much smaller ones can be utilized for the determination of quadrupole coupling constants and usually for the centrifugal distortion constant as well. Satisfactory quadrupole coupling constant values have been obtained for both the trans and cis conformer (Table 7) and centrifugal distortion constant values for the trans form. The structural parameters obtained from the earlier microwave study8 were obtained from the three rotational constants for each conformer. The parameters for the hydrogen atoms were fixed at the value reported from these parameters for cyclopropyl chloride29 and the structure for the NCO group from the reported parameters of vinyl isocyanate30,31 and isocyanic acid.32 The r(C1-C2) value was estimated from the planar moment33 and the C2-C3 distances of 1.515 Å for this molecule were taken from those of cyclopropylisothiocyanate (electron diffraction34). With these parameters fixed, the r(C-N), — CNC, and — CCN were obtained for both conformers by diagnostic least-squares adjustment.35 The previously assumed and determined parameters are listed in Table 5 and there are some significant differences from the parameters obtained in this study. For example, the assumed r(C1-C2,3) and r(CdO) seem entirely too long, whereas r(C1-H) is too short. Also, the determined difference of 0.010 Å between the two values of r(C1-N) is too 2304

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

Table 7. Rotational (MHz), Centrifugal Distortion (kHz), and Quadrupole Coupling Constants for the trans and cis Conformers of Cyclopropylisocyanate trans

cis

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

a

a

MW

MWa

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

A

17221

17016

17191

16941.889(8)

10359

9997

10705

10215(6)

B

1753

1769

1758

1784.3121(6)

2137

2208

2085

2186.859(6)

C DJ

1691 0.115

1706 0.114

1693 0.160

1716.1278(6) 0.153(2)

2062 1.504

2114 1.448

2029 1.550

2106.172(6) 2.01(1)

DJK

23.670

47.713

2.319

52.28(3)

3.140

-5.393

7.270

29.5(3)

DK

-2.843

-22.405

19.965

-6(3)

57.49

51.83

72.69

n.d

D1

-0.002

-0.001

-0.002

-0.027(1)

-0.22

-0.23

-0.22

-0.31(2)

D2

-0.002

-0.011

0.017

0.07

0.03

0.07

χaa

2.9179

2.9715

2.9504

2.6306(26)

2.6501

2.8092

2.8517

2.5647(49)

χbb

-1.7909

-1.8545

-1.5362

-1.3839(31)

-0.7920

-0.8843

-1.2248

-1.0624(75)

χcc

-1.1270

-1.1170

-1.4142

-1.2467(31)

-1.8581

-1.9249

-1.6269

-1.5022(75)

Ref 10; n.d.: not determined.

large because the predicted ab initio value is only 0.001 Å! Although the ab initio predicted value for the distance for one of the conformers may differ significantly from the experimental value, approximately the same difference is expected for the other conformer, so the predicted difference is usually quite accurately predicted. However, the listed uncertainty for the r(C1-N) distance could be interpreted as being essentially the same distance. The other parameter where there is a significant difference between the estimated value and the determined distance in this study is r(CdO), which is significantly shorter from this study. However, the previously determined — C2,3C1N and — C-NdC values are in good agreement with the value obtained for the corresponding parameters in the current study. In the second microwave study,10 the scientists were mainly interested in the value of the centrifugal distortion constants. The experimentally determined values of the quadratic distortion constants for the trans form are listed in Table 7 along with the value obtained from the force constants obtained from the ab initio and the density functional theory calculations. Often the prediction from B3LYP calculations are in better agreement than those from the MP2 prediction but in this case MP2/6-311þ G(d,p) values were quite superior, particularly for DJK and DK. However, for the cis form, reasonable centrifugal distortion values could not be obtained even by using the fourth order fit. In fact, these investigators10 initially suggested in the introduction that the conformational behavior of cyclopropylisocyanate raised some doubts about the existence of the cis form. Also, the authors10 attempted to explain the failure of the centrifugal distortion analysis for this conformer by using the basic assumptions of the conventional distortion theory for a rigid rotor Hamiltonian. For this Hamiltonian, the centrifugal distortion effects are described by a harmonic potential function and the theory allows only the small displacements of the internal coordinates from the equilibrium values to be considered. They explain that “These conditions are obviously not fulfilled for “cis” cyclopropylisocyanate. Because of the rigidity of the cyclopropyl frame and of the NCO group, especially the torsion around the Cframe-N bond, allows for flexibility of the molecule. If we assume a very flat and anharmonic potential near the “cis” conformation with a low barrier to the trans conformation, the failure of the centrifugal distortion analysis is plausible.”10 Nevertheless, the problem was not the question of cis being the correct

conformer, it is due to the need for a potential function that has been described as like an “inverted champagne bottle.”36 This type of potential simply arises from the large amplitude bending motion of the C-NCO linkage near the top of the potential well. It would be of interest to use such a potential to assign high J-level values, particularly values larger than J = 2. It was noted10 earlier that for K-1 = 3 the transitions showed a positive deviation, whereas those with K- > 3 were characterized by negative deviation. The authors of the theoretical paper9 seemingly did not take into consideration the information from the microwave studies.8,10 The dipole moment components determined8 from the Stark effects clearly indicated that the trans form had a plane of symmetry with |μb| = 0 and similarly for the cis conformer with |μc| = 0. The optimization of the two possible gauche equilibria9 from the MP2/6-31G(d,p) calculation gives the trans and cis conformations unless the gauche conformers are fixed as the structures. Apparently, from the earlier theoretical studies9 the investigators really did obtain the trans and cis conformers. In fact, the authors9 even acknowledged the similarity of the two “gauche equilibrium structures” with the trans and cis forms of the previous study.8 For the actual gauche and skew forms, the dipole moments obtained would have all three components unequal to zero (Table 2S). A more appropriate choice of the dihedral angle used would be CdN-C-H4 instead of CdN-C-C2. When CdN-C-H4 is used as the dihedral angle, the equilibrium structure for the trans and cis conformers would have a plane of symmetry and it would also clearly show how the NCO group is oriented relative to the ring if the molecule is in the two “gauche equilibrium structures.” We have carried out the MP2(full)/6-31G(d,p) calculation utilized in the earlier theoretical study.9 In Figure 6, the potential functions of MP2(full)/6-31G(d,p) and MP2(full)/6-311G(2d,2p) are shown, and the smaller basis set gives a much broader and shallower potential well for the trans conformer than the MP2(full)/6-311G(2d,2p) calculation. The depth of the well for the cis form is approximately the same from both calculations. Due to the breadth of this potential function from the smaller basis set, it is possible that the method utilized would make it difficult to determine the minimum because only two points were used between 90 and 180. Nevertheless, the MP2(full)6-31G(d, p) predicted the minimum at the trans and cis position. 2305

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

Figure 6. Potential energy curve of cyclopropylisocyanate as a function of the dihedral angle CdN-C1-H4 obtained from MP2(full)/ 6-31G(d,p) calculations (dashed line) and MP2(full)/6-311G(2d,2p) calculations (solid line).

If a more appropriate dihedral angle was used for indicating the conformers, the MP2(full)6-31G(d,p) from the theoretical study9 would have given the same result as this current study.

’ ASSOCIATED CONTENT

bS

Supporting Information. Table S1, symmetry coordinates for cyclopropylisocyanate; Table 2S, a comparison of cyclopropylisocyanate structural parameters (Å and degree), rotational constants (MHz), and dipole moments (debye) for the trans, cis, gauche, skew, gauche-1, and gauche-2 conformers from MP2/6-31G(d,p) calculation; and Figure 1S, van’t Hoff plot of the eight infrared conformer bands utilized in the determination of enthalpy difference value for cyclopropylisocyanate dissolved in liquid xenon solution. 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. Author Contributions †

Taken in part from the dissertation of S. X. Zhou, which will be submitted in partial fulfillment for the Ph.D. degree.

’ ACKNOWLEDGMENT J.R.D. acknowledges the University of Kansas City Trustees for a Faculty award and S.X.Z. acknowledges the Graduate Assistance Fund of the UMKC Women’s Council (Linda Hood Talbott Award, Eleanor Brantley Schwartz Award, Samuel L. and Peni P. Colville Award) for partial financial support of this research. ’ REFERENCES (1) Zheng, C.; Guirgis, G. A.; Deeb, H.; Durig, J. R. On the Structural Parameters and Vibrational Spectra of CH3NCS, SiH3NCS and GeH3NCS. J. Mol. Struct. 2007, 829, 88–110. (2) Zhou, S. X.; Durig, J. R. The r0 Structural Parameters, Vibrational Spectra, Ab Initio Calculations and Barriers to Internal

ARTICLE

Rotation and Linearity of Methylisocyanate. J. Mol. Struct. 2009, 924, 111–119. (3) Durig, J. R.; Zheng, C. Conformational Stability, Structural Parameters and Vibrational Assignment from Variable Temperature Infrared Spectra of Krypton Solutions and Ab Initio Calculations of Ethylisothiocyanate. Spectrochim. Acta 2007, 68, 783–795. (4) Durig, J. R.; Zhou, S. X.; Zheng, C.; Durig, D. T. Conformational Stability, Structural Parameters and Vibrational Assignment from Variable Temperature Infrared Spectra of Xenon and Krypton Solutions and Ab Initio Calculations of Ethylisocyanate. J. Mol. Struct. 2010, 971, 23–32. (5) Sakaizumi, T.; Ohashi, O.; Yamaguchi, I. Microwave Spectrum of Ethyl Isothiocyanate. Bull. Chem. Soc. Jpn. 1976, 49, 948–953. (6) Palmer, M. H.; Nelson, A. D. The Structures of the Azido-, Isocyanato- and Isothiocyanato-Derivatives of Methane and Silane and their Derivatives. A Comparison of Ab Initio with Experimental Results. J. Mol. Struct. 2004, 689, 161–173. (7) Durig, J. R.; Zheng, C. Conformational Stability, Structural Parameters and Vibrational Assignment from Variable Temperature Infrared Spectra of Krypton Solutions and Ab Initio calculations of Ethylisothiocyanate. Spectrochim. Acta 2007, 68, 783–795. (8) Durig, J. R.; Berry, R. J.; Wurrey, C. J. Microwave, Infrared and Raman Spectra, Conformational Stability, Structure, Dipole Moment and Vibrational Assignment for Cyclopropylisocyanate. J. Am. Chem. Soc. 1988, 110, 718–726. (9) Nguyen, M. T.; Hajnal, M. R.; Vanquickenborne, L. Q.; Ha, T.-K.; Stohner, J. Structure and Conformation of Chlorosulfonylisocyanate and Cyclopropylisocyanate. J. Chem. Soc., Faraday Trans. 1993, 89, 2381–2384. (10) Heldmann, C.; Dreizler, H. Investigation of the Microwave Spectrum of Cyclopropyl Isocyanate. An Example for the Failure of Centrifugal Distortion Theory. Z. Naturforsch. 1990, 45, 1175–1184. :: (11) Kricheldorf, H. R.; Regel, W. Reaktionen mit Silylaziden, 4. Uber Cyclopropylisocyanate. Chem. Ber. 1973, 106, 3753–3764. (12) Kricheldorf, H. R.; Regel, W. Reaktionen mit Silylaziden, 5. Notiz €uber die Acylierung und Alkylierung von Trimethylsilylazid. Chem. Ber. 1973, 106, 3765–3768. (13) 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.; 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. (14) Pulay, P. Ab Initio Calculation of Force Constants and Equilibrium Geometries. I. Theory. Mol. Phys. 1969, 17, 197–204. (15) Møller, C.; Plesset, M. S. Note on an Approximation Treatment for Many-Electron Systems. Phys. Rev. 1934, 46, 618–622. (16) Guirgis, G. A; Zhu, X.; Yu, Z.; Durig, J. R. Raman and Infrared Spectra, Conformational Stability, Normal Coordinate Analysis, Vibrational Assignment, and Ab Initio Calculations of 3,3-Difluorobutene. J. Phys. Chem. A 2000, 104, 4383–4393. (17) Frisch, M. J.; Yamaguchi, Y.; Gaw, J. F.; Schaefer, H. F., III; Binkley, J. S. Analytic Raman Intensities from Molecular Electronic Wave Functions. J. Chem. Phys. 1986, 84, 531–532. (18) Amos, R. D. Calculation of Polarizability Derivatives Using Analytic Gradient Methods. Chem. Phys. Lett. 1986, 124, 376–381. 2306

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307

The Journal of Physical Chemistry A

ARTICLE

(19) Polavarapu, P. L. Ab Initio Vibrational Raman and Raman Optical Activity Spectra. J. Phys. Chem. 1990, 94, 8106–8112. (20) Chantry, G. W. In The Raman Effect; Anderson, A., Ed.; Marcel Dekker, Inc.: New York, NY, 1971; Vol. 1, Chap. 2. (21) Van der Veken, B. J.; Herrebout, W. A.; Durig, D. T.; Zhao, W.; Durig, J. R. Conformational Stability of 3-Fluoropropene in Rare Gas Solutions from Temperature-Dependent FT-IR Spectra and Ab Initio Calculations. J. Phys. Chem. A 1999, 103, 1976–1985. (22) Durig, J. R.; Ng, K. W.; Zheng, C.; Shen, S. Comparison of Ab Initio MP2/6-311þG(d,p) Predicted Carbon-Hydrogen Bond Distances with Experimentally Determined r0 (C-H) Distances. Struct. Chem. 2004, 15, 149–157. (23) McKean, D. C. New Light on the Stretching Vibrations, Lengths and Strengths of CH, SiH and GeH Bonds. J. Mol. Struct. 1984, 113, 251–266. (24) Bulanin, M. O. Infrared Spectroscopy in Liquefied Gases. J. Mol. Struct. 1973, 19, 59–79. (25) Van der Veken, B. J.; DeMunck, F. R. An Infrared Study of Monomeric and Oligomeric (n = 2, 3, and 4) Hydrogen Chloride in Liquified Noble Gases. J. Chem. Phys. 1992, 97, 3060–3072. (26) Herrebout, W. A.; van der Veken, B. J.; Wang, A.; Durig, J. R. Enthalpy Difference between Conformers of n-Butane and the Potential Function Governing Conformational Interchange. J. Phys. Chem. 1995, 99, 578–585. (27) Bulanin, M. O. Spectroscopy of Molecules in Liquid Noble Gases. J. Mol. Struct. 1995, 347, 73–82. (28) Galabov, B; Kim, S.; Xie, Y.; Schaefer, H. F., III; Leininger, M. L.; Durig, J. R. Conformations of Allyl Amine: Theory vs Experiment. J. Phys. Chem. A 2008, 112, 2120–2124. (29) Schwendeman, R. H.; Jacobs, G. D.; Krigas., T. M. Molecular Structure of Cyclpropyl Chloride. J. Chem. Phys. 1964, 40, 1022–1028. (30) Kirby, C.; Kroto, H. W. Microwave and Photoelectron Study of cis- and trans-Isocyanatoethene, CH2:CHNCO (Vinyl Isocyanate). J. Mol. Spectrosc. 1978, 70, 216–228. (31) Bouchy, A.; Roussy, G. Geometrical and Electronic Structure of Two Conformers of Vinyl Isocyanate by Microwave Spectroscopy. J. Mol. Spectrosc. 1977, 68, 156–65. (32) Yamada, K. Molecular Structure and Centrifugal Distortion Constants of Isocyanic Acid from the Microwave, Millimeterwave, and Far Infrared Spectra. J. Mol. Spectrosc. 1980, 79, 323–346. (33) Penn, R. E.; Bogg, J. E. Substituent-Induced Asymmetry of Cyclopropane Ring. J. Chem. Soc., Chem. Commun. 1972, 666–667. (34) Durig, J. R.; Sullivan, J. F.; Berry, R. J.; Cradock, S. An Electron Diffraction and Ab Initio Investigation of Cyclopropylisothiocyanate. J. Chem. Phys. 1987, 86, 4313–4318. (35) Groner, P.; Church, J. S.; Li, Y. S.; Durig, J. R. Spectra and Structure of Organophosphorus Compounds. XXIII. Microwave Spectra, Electric Dipole Moment and Molecular Structure of Two Conformers of Ethyldifluorophosphine. J. Chem. Phys. 1985, 82, 3894–3902. (36) Winnewisser, M.; Winnewisser, B. P.; Medvedev, I. R.; de Lucia, F. C.; Ross, S. C.; Bates, L. M. The Hidden Kernel of Molecular QuasiLinearity: Quantum Monodromy. J. Mol. Struct. 2006, 798, 1–26.

2307

dx.doi.org/10.1021/jp108879u |J. Phys. Chem. A 2011, 115, 2297–2307