A Tale of Two Carenes: Intrinsic Optical Activity ... - ACS Publications

Aug 13, 2012 - Lykourgos Bougas , Dimitris Sofikitis , Georgios E. Katsoprinakis ... K. Spiliotis , Paraskevas Tzallas , Benoit Loppinet , T. Peter Ra...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCA

A Tale of Two Carenes: Intrinsic Optical Activity and Large-Amplitude Nuclear Displacement Priyanka Lahiri, Kenneth B. Wiberg, and Patrick H. Vaccaro* Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States S Supporting Information *

ABSTRACT: The specific rotation for two isomeric members of the terpene family, (S)(+)-2-carene and (S)-(+)-3-carene, has been investigated under complementary solvated and isolated conditions, where the latter vapor-phase work has been performed at excitation wavelengths of 355 and 633 nm by means of ultrasensitive cavity ring-down polarimetry (CRDP). Linear-response computations of dispersive optical activity built upon analogous density-functional (B3LYP/aug-cc-pVTZ) and coupled-cluster (CCSD/aug-cc-pVDZ) levels of theory have been enlisted to unravel the structural and electronic origins of observed behavior. The six-membered portion of the bicyclic skeleton in the nominally rigid 3-carene system is predicted to be near-planar in nature, with calculated and measured rotatory powers for the isolated (gas-phase) species shown to be in excellent agreement. In contrast, the inherent flexibility of 2-carene gives rise to two quasidegenerate conformations that are interconnected by a large-amplitude ringpuckering motion and exhibit antagonistic chiroptical properties. Various approaches to simulate the intrinsic response evoked from a thermally equilibrated ensemble of gaseous (S)-(+)-2-carene molecules have been considered, including implicit averaging over independent conformers and explicit (albeit restricted) averaging over nuclear degrees of freedom. A polarizable continuum model for implicit solvation was found to describe solvent-dependent trends reasonably well in the case of (S)-(+)-2-carene, but failed to reproduce the specific-rotation patterns emerging from polarimetric studies of (S)-(+)-3-carene.

I. INTRODUCTION Recent years have witnessed a veritable renaissance in the application of chiroptical probes to problems of molecular structure and dynamics,1 with numerous experimental and theoretical efforts highlighting the ability to determine the absolute stereochemical configuration and to extract ancillary conformational information for diverse chemical species. To a large extent, this resurgence of interest has been stimulated by the advent of robust and efficient quantum-chemical paradigms for reliably predicting the spectral response (magnitude and sign) evoked from a chiral molecule; 2 however, the complications incurred by extrinsic environmental perturbations (e.g., solvation processes) and intrinsic nuclear couplings (e.g., vibrational motion as well as large-amplitude displacements) still present formidable challenges.3 The present work strives to elucidate the provenance of dispersive electronic optical activity (circular birefringence or CB)4,5 by examining the wavelength-dependent specific rotation (optical rotatory dispersion or ORD) for two isomeric members of the terpene family under complementary solution-phase and vapor-phase conditions, where the latter isolated-molecule studies have exploited the ultrasensitive techniques of cavity ring-down polarimetry (CRDP).6,7 Figure 1 depicts the two molecules of interest for the present work, (a) (1S,6R)-3,7,7-trimethylbicyclo[4.1.0]hept-2-ene and (b) (1S,6R)-3,7,7-trimethylbicyclo[4.1.0]hept-3-ene, which will be referred to as (S)-(+)-2-carene and (S)-(+)-3-carene in the ensuing discussion. These species clearly are structural isomers of one another that differ primarily by position of the olefinic © 2012 American Chemical Society

Figure 1. Structures of the targeted carenes. The minimum-energy structures of (a) (S)-(+)-2-carene and (b) (S)-(+)-3-carene, as predicted by the B3LYP/apVTZ level of theory, are depicted, with the numbering scheme used to distinguish the carbon atoms composing the bicyclic ring system being superimposed.

moiety within the six-membered portion of the bicyclic ring systema modification in chromophore topology that might be expected to affect the nature and location of excited electronic states only slightly. Nevertheless, previous studies performed in the condensed phase have documented the substantial changes in optical activity imbued by this modest alteration in skeletal bonding patterns, with the highest specificrotation values reported for sodium D-line (589.3 nm) Received: April 5, 2012 Revised: August 9, 2012 Published: August 13, 2012 9516

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

excitation of neat (S)-(+)-2-carene8 and neat (S)-(+)-3-carene9 being +97.7 deg dm−1 (g/mL)−1 and +17.7deg dm−1 (g/mL)−1, respectively. The molecular geometries in Figure 1 display subtle differences in puckering of the carbon skeleton accompanied by reorientation of the methyl groups to minimize steric hindrance. Such structural preferences have been the subject of numerous experimental and theoretical efforts that, oftentimes, have led to contradictory conclusions. An early investigation of Kerr constants in 3-carene suggested an equilibrium to exist between boat and inverted-boat conformations,10 with the former dominating under ambient conditions. This assertion was supported, in part, by the product distribution generated from the sensitized photo-oxygenation of 3-carene,11 which was believed to stem from a rapidly interconverting 1:1 mixture of the aforementioned conformers; however, the interference patterns observed in gas-phase measurements of electron diffraction could be simulated well by considering only the boat motif.12 While the temperature dependence of chemical shifts observed for the methyl protons of 3-carene was interpreted in terms of a conformational equilibrium in which the dominant boat form was stabilized by an estimated 6.2 kJ mol−1 relative to its inverted-boat counterpart,13 analogous proton-NMR studies performed at higher magnetic-field strengths argued that the six-membered portion of the bicyclic ring is essentially planar in nature.14 Conflicting assertions also have emerged from theoretical analyses of 3-carene. While early molecular-mechanics simulations by Favini and co-workers15 reported a stable planar arrangement for the six-fold carbon ring, geometry optimizations performed at minimal Hartree−Fock (HF/STO-3G) levels of theory suggested the inverted-boat conformer to reside 10.7 kJ mol−1 lower in energy than the boat form,16 a prediction at odds with prior experimental findings. A subsequent MMP2 molecular-mechanics study performed in conjunction with measurements of NMR proton−proton coupling constants lent further support to the existence of a planar ring motif.17 In the most recent computational investigations of 3-carene, Stephens et al.18 exploited an extensive Monte Carlo search algorithm based upon the MM94 force field to locate a single, globally stable conformation, the validity of which was confirmed by examining the potential surface as a function of the C 3C 4C 5 C 6 dihedral angle (cf., Figure 1). This coordinate was found to be essentially nil (0°) at the minimum-energy configuration, in keeping with the expectation for a planar six-membered carbon skeleton. The refined structure obtained from this work subsequently was employed for density-functional calculations of chiroptical properties by means of B3LYP/6-31G* and B3LYP/aug-cc-pVDZ linearresponse methods, leading to an estimated specific-rotation value of +31.9 deg dm−1 (g/mL)−1 for the sodium D-line excitation of isolated (solvent-free) (S)-(+)-3-carene molecules. In contrast to 3-carene, few experimental and theoretical studies of conformational dynamics have been reported for 2carene. Two stable forms adopting pseudochair and pseudoboat motifs have been proposed for this system,13 with chemicalshift measurements and chemical-reactivity arguments suggesting the former to be more stable than the latter. Molecularmechanics simulations15 and minimal (HF/STO-3G) geometry-optimization procedures19 supported this conjecture, predicting the attendant energy difference (between conformers) to be 8.8 and 4.9 kJ mol−1, respectively. The ensuing

analyses build upon such prior efforts to elucidate the structural preferences of (S)-(+)-2-carene and (S)-(+)-3-carene, which will be shown to manifest themselves prominently through distinct chiroptical properties.

II. EXPERIMENTAL PROCEDURES AND COMPUTATIONAL METHODS Requisite samples of (S)-(+)-2-carene and (S)-(+)-3-carene were obtained from a commercial source (Sigma-Aldrich) with specified chemical purities of 97% and 99%, respectively, and used without further treatment. Corresponding percentages of enantiomeric excess (% ee) were estimated to be >92% and >97% on the basis of solvated and neat specific-rotation measurements performed at the sodium D-line (589.3 nm). Materials employed for vapor-phase studies were sealed in glass vessels having grease-free (Teflon) vacuum stopcocks and subjected to at least three freeze−pump−thaw cycles so as to minimize potential contamination from entrained (atmospheric) gases. Measurements performed on isolated (vapor-phase) species exploited the basic CRDP instrumentation and analysis procedures reported elsewhere,6,7 with the most salient details being discussed briefly here. Requisite visible and ultraviolet light was generated by pumping a high-resolution pulsed dye laser (Lambda Physik FL3002E; 0.035 cm−1 bandwidth) with the second harmonic of a Nd:YAG system (Spectra-Physics GCR-4-20; 20 pps repetition rate; ∼10 ns pulse duration). While the fundamental dye-laser output was used for 633 nm work, studies performed at 355 nm necessitated the additional step of frequency doubling (Inrad Autotracker III; BBO crystal) followed by isolation of the emerging ultraviolet radiation. The resulting linearly polarized 633/355 nm beam propagated through a variable attenuator before entering a Keplerian telescope where it was spatially filtered and refocused to match the mode structure of a ring-down cavity having a total length of L = 170.7284(35) cm, where the value in parentheses denotes the one standard-deviation uncertainty for the final two significant digits. The mode-matched excitation light passed though a circular polarizer (consisting of a tandem calcite prism and λ/4retarder) prior to entering the CRDP linear resonator assembly through the planar rear surface of one cavity mirror. The concave mirrors employed for 355 and 633 nm studies (Los Gatos Research; 6 m radius of curvature) had specified intensity reflectivities of ≥99.95% and ≥99.99%, respectively. Two intracavity quarter-wave plates of composite zero-order construction (Alpine Research Laboratories; optically contacted with antireflection coatings exhibiting [α]Tλ,gas, an assertion diametrically 9526

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

extracted for 355 and 633 nm excitation are plotted against the Onsager dielectric-constant function:55,56 f (ε) =

ε−1 2ε + 1

(8)

which spans the range 0 ≤ f(ε) ≤ 0.5 as the surrounding medium varies from being completely nonpolar (ε = 1) to being infinitely polar (ε ≫ 1). In particular, this quantity stems from the progenitor of all implicit solvation models known as Onsager reaction-field theory,55 where it mediates the solvent− solute stabilization energy engendered by the lowest-order multipole moment (i.e., the permanent electric dipole moment) of the neutral solute molecule (vide infra). Excluding the anomalous behavior of chloroform that has been attributed to dominant quadrupolar interactions,32,57 a strong linear relationship is found between [α]Tλ and f(ε) for the remaining solvents, with extrapolation to the origin (ε = 1) suggesting the attendant vapor-phase specific rotations to be 77.5 and 14.1 deg dm−1 (g/mL)−1 at 355 and 633 nm, respectively. These predictions differ markedly from the present CRDP measurements of intrinsic optical activity, reinforcing prior assertions that a serious incongruity exists between dispersive chiroptical studies performed under solvated and isolated conditions.7 Table 9 contains a subset of the optical-activity data acquired for (S)-(+)-2-carene in five different solvents, with Figure 8 depicting associated ORD profiles. In contrast to 3-carene, solvation now causes a modest increase in the magnitude of rotatory power. When compared to their vapor-phase counterparts (cf., Table 5 or 6), solution-phase [α]Tλ parameters extrapolated for 355 nm (633 nm) excitation rise by +26.3% (+16.8%) and +10.6% (+4.2%) in the presence of cyclohexane and acetonitrile media, thus making the latter, once again, the best mimic for isolated-molecule measurements. Also of note is the overall tendency for specific rotation to decrease in proportion to solvent polarity; however, both ethanol and chloroform appear to deviate from this pattern. The inherent nonrigidity of the bicyclic framework in 2-carene can be expected to complicate interpretation of these results, leading to patterns that reflect the modification of chiroptical response by the surrounding environment as well as the differential stabilization experienced by individual conformers.37 Some insight can be gleaned by considering the nature of solvation processes, which often are dominated by electrostatic effects built upon the lowest-order multipole moments of the entrained solute. The most rudimentary implementation of Onsager reaction-field theory predicts a solute−solvent interaction energy that scales as Esol ∝ −f(ε)μ2,55,56 where μ denotes the magnitude of the permanent electric dipole moment for the solute. Consequently, species possessing larger values of μ will be relaxed to a greater extent than their lesspolar analogues, with the overall degree of stabilization increasing in proportion to the polarity of the solvent as encoded in the dielectric-constant function f(ε). While the magnitude of the permanent electric dipole moment estimated for isolated 3-carene molecules is rather small (0.169 D from B3LYP/apVTZ), this quantity grows to 0.364 and 0.358 D for the (A) and (B) forms of 2-carene, respectively, and rises further to 0.374 D for the transition-state configuration. The differences in μ among these 2-carene structures are slight; however, the preferential stabilization of polar species in high-dielectric surroundings might be expected to increase the energy separation between conformers and (potentially) decrease the accompanying barrier to isomer-

Figure 8. ORD of solvated (S)-(+)-2-carene. The specific-rotation values measured at discrete excitation wavelengths (closed symbols and solid curves) for (S)-(+)-2-carene are shown in the top panel for a variety of solvents. The accompanying vapor-phase data at 355 and 633 nm (isolated symbols) highlight the role of environmental perturbations. The bottom panel depicts ORD profiles predicted by a thermal conformer-averaging analyses (T = 300 K) based on structures optimized at the B3LYP/apVTZ level of theory in the presence of implicit (PCM) acetonitrile solvation. These results follow from B3LYP/apVTZ linear-response calculations, with the energy (ΔEη) and free-energy (ΔGη) estimates for individual conformers being obtained from three computational methods: B3LYP/apVTZ, G3B3, and G3.

ization. Such assertions are mostly in accord with the computational results compiled in Table 10, which highlight the solvent-dependent energies of stationary points optimized at the B3LYP/apVTZ level of theory in the presence of a PCM treatment for solvation. In particular, the B3LYP/apVTZ, G3B3, and G3 approaches concur in predicting the relative energy of 2-carene(B) to climb upon moving from a weakly polar cyclohexane environment to a strongly polar acetonitrile medium. Table 11 contains wavelength-resolved values of specific rotation computed for the three stationary points of groundstate (S)-(+)-2-carene molecules entrained in the five solvents targeted by the present study. These results follow from application of the B3LYP/apVTZ linear-response formalism to geometries optimized at the same level of density-functional theory, with both analyses being performed in the presence of a PCM treatment for solvation. As such, the tabulated quantities embody changes effected in both structural and chiroptical properties. The putative role of the latter was accessed by 9527

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

Table 10. Predicted Relative Energy for Stationary Points of Solvated (S)-(+)-2-Carenea relative energy (cm−1)

population fraction

source

solvent

conformer (B)

transition state

conformer (A)

conformer (B)

ΔEη [B3LYP/apVTZ]

acetonitrile ethanol chloroform dibutyl ether cyclohexane

57.79 56.74 44.80 38.88 32.18

449.00 449.02 449.20 449.24 449.22

0.569 0.568 0.554 0.546 0.539

0.431 0.432 0.446 0.454 0.461

ΔGη [G3B3]

acetonitrile ethanol chloroform dibutyl ether cyclohexane

80.33 79.01 63.87 56.62 48.50

731.51 731.07 725.36 723.39 721.63

0.595 0.594 0.576 0.567 0.558

0.405 0.406 0.424 0.433 0.442

ΔGη [G3]

acetonitrile ethanol chloroform dibutyl ether cyclohexane

109.74 109.08 100.96 96.13 90.86

··· ··· ··· ··· ···

0.629 0.628 0.619 0.613 0.607

0.371 0.372 0.381 0.387 0.393

a Relative energies for the (A) and (B) conformers of solvated 2-carene molecules, as well as for the associated transition-state (TS) configuration, are tabulated for five different solvents. The predicted energy (ΔEη) and free-energy (ΔGη) parameters follow from B3LYP/apVTZ, G3B3, and G3 calculations performed on geometries optimized at the B3LYP/apVTZ level of theory, with all analyses incorporating implicit PCM solvation. The listed fractional populations of 2-carene(A) and 2-carene(B) are evaluated for thermal equilibrium at T = 300 K.

Table 11. Predicted Specific Rotation for Stationary Points of Solvated (S)-(+)-2-Carenea specific optical rotation [deg dm−1 (g/mL)−1] species

wavelength (nm)

acetonitrile

ethanol

chloroform

dibutyl ether

cyclohexane

(S)-2-carene A

355.00 365.02 435.83 546.07 589.30 633.00

−231.29 −208.39 −117.80 −64.86 −53.94 −45.60

−234.53 −211.30 −119.40 −65.72 −54.65 −46.20

−262.08 −235.89 −132.56 −72.59 −60.29 −50.91

−268.30 −241.30 −135.08 −73.77 −61.23 −51.68

−281.76 −253.28 −141.39 −77.02 −63.89 −53.90

(S)-2-carene TS

355.00 365.02 435.83 546.07 589.30 633.00

739.23 663.85 362.54 190.65 156.51 130.88

738.49 663.17 362.16 190.45 156.33 130.73

746.34 669.84 364.81 191.45 157.08 131.31

731.54 656.69 358.02 188.03 154.30 129.00

731.04 656.11 357.37 187.56 153.90 128.64

(S)-2-carene B

355.00 365.02 435.83 546.07 589.30 633.00

1223.26 1102.40 613.89 328.66 271.03 227.50

1229.02 1107.49 616.48 329.94 272.07 228.36

1254.54 1129.93 627.55 335.31 276.39 231.92

1238.59 1115.68 620.01 331.45 273.25 229.30

1247.06 1123.07 623.49 333.09 274.56 230.37

a The wavelength-resolved specific rotation computed for the solvated (A) and (B) conformers of (S)-(+)-2-carene, as well as for the attendant transition state (TS), are presented for five solvents. These results follow from application of B3LYP/apVTZ linear-response methods to structures refined at the same level of density-functional theory, with all analyses incorporating a PCM treatment of implicit solvation.

repeating calculations of dispersive optical activity without the PCM ansatz while still retaining solvent-optimized configurations of the nuclear framework. This procedure led to [α]Tλ parameters (for the unrelaxed “gas-phase” species) that increased in magnitude by ≤10%, again contradicting the behavior anticipated from the Lorentz local-field correction of eq 7. At a given excitation wavelength, the magnitude of specific rotation computed for each stationary point of (S)-(+)-2-carene

tends to decrease with increasing dielectric constant of the surroundings; however, clear deviations from this pattern are evident for chloroform and ethanol. The most favorable comparison between B3LYP/apVTZ solution-phase and vapor-phase predictions of dispersive optical activity are found in the case of acetonitrile, where solvated [α]λT parameters for 2-carene(A) differ from their isolated-molecule counterparts (cf., Table 3) by near-negligible amounts. In contrast, analogous acetonitrile results for 2-carene(B) and 29528

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

activity performed at the B3LYP/apVTZ and CCSD/apVDZ (MVG) levels of theory were found to be in excellent agreement with each other and to reproduce isolated-molecule measurements of specific rotation in a quantitative fashion. As such, the intrinsic chiroptical properties of (S)-(+)-3-carene would appear to be described well by purely electronic effects, with cumulative contributions arising from nuclear degrees of freedom (e.g., small-amplitude vibrational motion) ostensibly being of negligible consequence. Two nearly degenerate conformers, which interconvert along a large-amplitude ring-puckering coordinate that presents a potential barrier of finite (400−800 cm−1) height, have been identified for (S)-(+)-2-carene. These species are predicted to have antagonistic chiroptical properties, with the wavelengthresolved (positive) rotatory power for 2-carene(B) being four times larger in magnitude than the corresponding (negative) quantity for 2-carene(A) while the transition-state configuration (which supports a quasiplanar six-membered ring akin to that in 3-carene) displays an intermediate (albeit positive) behavior. Despite the uniformity attained among theoretical calculations in (S)-(+)-3-carene, coupled-cluster (MVG) estimates of dispersive optical activity for all stationary points in (S)(+)-2-carene are almost a factor of 2 smaller than their densityfunctional counterparts at each targeted wavelength. Various methods to account for the inherent structural flexibility of 2-carene have been considered, including Boltzmann-weighted averaging over properties of independent conformers and restricted vibrational averaging over torsional displacement of the bicyclic ring system. A conformer-averaging procedure, whereby relative free-energies estimated for each equilibrium configuration from composite G3B3 analyses (ΔGη) are combined with corresponding B3LYP/apVTZ eq ), affords the best calculations of specific rotation ([α]λ,η simulation of observed isolated-molecule behavior, with use of analogous CCSD/apVDZ estimates for [α]eq λ,η significantly reducing the quality of theoretical results. The overall agreement achieved between the conformer-averaging and vibration-averaging treatments confirms the reliability of the former scheme (despite the rudimentary nature of its underlying assumptions)40 and reinforces the important role played by the large-amplitude, ring-puckering degree of freedom. The specific rotation of (S)-(+)-3-carene drops markedly upon solvation, with acetonitrile offering the best mimic for isolated-molecule behavior. While dispersive optical-activity measurements in the vapor phase are reproduced well by B3LYP/apVTZ and CCSD/apVDZ calculations, attempts to describe solution-phase results by incorporating a PCM ansatz into the linear-response framework of density functional theory proved to be less successful. Such analyses were found to shift solution-phase ORD curves uniformly in the wrong direction (toward larger magnitudes) relative to their vapor-phase counterpart; however, the overall growth in observed specific rotation with decreasing solvent polarity was reproduced successfully, including the deviation of a chloroform medium from this basic trend. In contrast to (S)-(+)-3-carene, (S)-(+)-2-carene undergoes a modest increase in rotatory power when introduced into the solution phase, with acetonitrile, once again, found to best duplicate the intrinsic (vapor-phase) response. The observed specific rotation tends to increase with decreasing solvent polarity; however, both ethanol and chloroform deviate from this general pattern. The PCM treatment (used to calculate

carene(TS) deviate by 15% or more. Given the discussion above, such observations can be attributed primarily to the structural relaxation brought about by solute−solvent interactions. The condensed-phase specific rotation exhibited by a thermally equilibrated ensemble of (S)-(+)-2-carene molecules has been simulated under the conformer-averaging ansatz of eq 3, with requisite fractional populations fη(T) following from B3LYP/apVTZ, G3B3, and G3 calculations performed in the presence of implicit (PCM) solvation (cf., Table 10). In the case of acetonitrile, excellent agreement between experiment and theory is obtained by using B3LYP/apVTZ estimates of ΔEη, which give a deviation of +2.2% for sodium D-line excitation. Analogous treatments based on more robust G3B3 and G3 predictions of ΔGη increase the magnitude of discrepancies at 589.3 nm to −7.9% and −20.8%, respectively. Similar levels of accord are found for the other solvents and the tendency for [α]Tλ at a fixed wavelength to decrease in proportion to the polarity of the surroundings is reproduced uniformly. The solvent dependence of dispersive optical activity in (S)(+)-2-carene is summarized by the top panels of Figure 8, which plot [α]Tλ parameters extrapolated for 355 and 633 nm excitation against the Onsager dielectric-constant function, f(ε). In contrast to the kindred analysis for (S)-(+)-3-carene, the graphs now display an inverse (negative) slope with the points for chloroform and ethanol departing markedly from those of other media. The latter behavior reflects the putative action of nondipolar effects32,57 as well as the influence of other processes (e.g., specific solute−solvent couplings) that mediate conformational and chiroptical properties. A linear relationship is obtained among the acetonitrile, di-n-butyl ether, and cyclohexane results, with the projected specific rotations for the nominally isolated solute molecule (ε = 1) being 468.8 and 84.8 deg dm−1(g/mL)−1 at 355 and 633 nm, respectively. Once again, CRDP measurements of intrinsic rotatory power depart appreciably from these estimates (by +36.8% at 355 nm and +26.2% at 633 nm), displaying characteristics more reminiscent of those obtained in strongly polar environments.

IV. SUMMARY AND CONCLUSIONS A synergistic experimental and computational investigation of dispersive optical activity has been performed for two isomeric members of the terpene family, (S)-(+)-2-carene and (S)(+)-3-carene, under ambient (room-temperature) conditions. The specific rotation of isolated molecules in the vapor phase was measured at 355 and 633 nm by exploiting the ultrasensitive methods of cavity ring-down polarimetry.6,7 Complementary condensed-phase work was conducted at discrete excitation wavelengths spanning the visible and ultraviolet regions of the spectrum by using solvent media selected to display a wide range of chemical and physical properties. Linear-response calculations based on the densityfunctional (B3LYP/apVTZ) and coupled-cluster (CCSD/ apVDZ) frameworks were enlisted to unravel the provenance of observed chiroptical behavior, with a polarizable continuum model (PCM) for implicit solvation being used to explore the putative roles of solute−solvent interactions. In the case of (S)-(+)-3-carene, quantum-chemical calculations confirmed the presence of a single, low-lying conformer where the six-membered portion of the bicyclic ring is held in a near-planar configuration by a nominally rigid molecular framework. Linear-response predictions of dispersive optical 9529

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

ΔGη and [α]eq λ,η for the conformer-averaging procedure) appears to afford a reasonable description for solute−solvent interactions in the case of (S)-(+)-2-carene, reflecting, in part, differential stabilization of conformers that exhibit antagonistic chiroptical properties. Nevertheless, analyses of continuumdielectric scaling based on a rudimentary implementation of Onsager reaction-field theory reinforce previous assertions that a serious incongruity exists between specific-rotation values extracted from solvated and isolated environments,7 with the latter routinely displaying characteristics reminiscent of those encountered in nonintuitive solvents of high polarity. The transposition of single/double carbon−carbon bonds between the isomers of carene does alter electronic structure slightly, with TDDFT(B3LYP)/apVTZ and EOM-CCSD/ apVDZ calculations showing the lowest-lying excited singlet state in 3-carene to reside at 5.30 and 6.00 eV, respectively, while the analogous feature in 2-carene(A) is predicted to be roughly 6% lower in energy (cf., Tables S6 and S7 in the Supporting Information). Nevertheless, for the nonresonant excitation wavelengths examined during the present study, the dispersive chiroptical response is not dominated by a particular (electronic) eigenstate, but instead reflects a subtle superposition of contributions arising from the entire manifold of excited states. As such, the marked differences observed in the rotatory powers of the targeted carene species might be attributable to factors other than those of purely electronic origin. Since the excellent agreement found between theory and experiment for the specific rotation of isolated (S)-(+)-3-carene molecules would tend to discount effects arising from smallamplitude vibrations, the pronounced ring-puckering motion that mediates interconversion between the quasidegenerate (A) and (B) conformers of (S)-(+)-2-carene becomes a prime suspect. Indeed, inspection of geometries predicted for the two equilibrium configurations (cf., Figure 3) reveals the sixmembered portion of the bicyclic-ring system to be twisted helically in opposite directions, thereby affording a viable explanation for the antagonistic chiroptical properties exhibited by these isomeric forms. Further computational and experimental efforts designed to unravel the dependence of electronic optical activity on nuclear degrees of freedom are underway, with special emphasis being directed toward molecular frameworks that are subject to large-amplitude displacements.

internal coordinates, with the corresponding determinant of g being denoted by g  ∥g∥. Straightforward expansion and simplification of eq A1 yields59 T̂ =T̂′ + V̂ ′ ⎧⎛ jk ⎞⎛ ⎞ ⎪ ∂g ⎟⎜ ∂ ln g ⎟ ⎨⎜ ⎜ ⎜ ⎟ ∂q ∂qk ⎟⎠ j,k j,k ⎪ ⎩⎝ j ⎠⎝ ⎤⎫ ⎡⎛ 2 ⎞ ⎞ ⎛ 1 ⎜ ∂ ln g ⎟⎛⎜ ∂ ln g ⎞⎟⎥⎪ jk ⎢⎜ ∂ ln g ⎟ ⎬ +g ⎜ + ⎜ ⎢ ∂q ∂q ⎟ 4 ⎝ ∂qj ⎟⎠⎜⎝ ∂qk ⎟⎠⎥⎦⎪ ⎣⎝ j k ⎠ ⎭

=−

ℏ2 −1/4 g 2

∑ j,k

∂ 1/2 jk ∂ −1/4 g g g ∂qj ∂qk



∂ jk ∂ ℏ2 g + 8 ∂qj ∂qk

3N − 3



where the symmetric nature of the metric tensor, gjk=gkj, has been exploited to arrive at the final expression. The second term in eq A2 represents an effective potential (or pseudopotential), V̂ ′, that arises from the kinetic energy operator. Although this quantity can be incorporated into the actual potential energy, V̂  V(q), its magnitude, which reflects the variation of elements gjk and determinant g with respect to internal coordinates, often is negligibly small in comparison to that of V̂ .60 Under such circumstances, the kinetic energy operator becomes T̂ ≈ T̂ ′ = −

ℏ2 2

3N − 3

∑ j,k

⎧ ⎫ ⎛ jk ⎞ ⎪ jk ∂ 2 ∂g ⎟ ∂ ⎪ ⎨g ⎬ + ⎜⎜ ⎟ ⎪ ∂qj∂qk ⎝ ∂qj ⎠ ∂qk ⎪ ⎩ ⎭

(A3)

where the second term in the summation can be discounted if gjk displays only a weak dependence on qj. Likewise, the mixed partial derivatives of the first term would vanish if the transformation between Cartesian and curvilinear coordinates leads to a diagonal metric tensor, gjk = gjj δjk (where δjk denotes the canonical Kronecker delta symbol), thus allowing (gjj)−1 to be associated with an effective mass factor for displacement along qj. By construction, the kinetic energy operator of eq A3 embodies both rotational and vibrational degrees of freedom, with g capable of being partitioned into two blocks (representing pure-rotational and pure-vibrational contributions) that are coupled by off-diagonal terms.46,60 The functional form of gjk for pure-rotational motion and for rotation−vibration interaction depends on how the body-fixed axis system is attached to the molecular framework; however, this is not the case for pure-vibrational elements and for the metric-tensor determinant.46,59 Consequently, any convenient choice of axis embedding can be used to evaluate the latter quantities.

APPENDIX A For evaluation of the kinetic energy associated with the torsional isomerization of 2-carene, it is assumed that the center-of-mass motion for a free molecule comprised of N atoms has been separated, leaving 3N − 3 internal degrees of freedom to partition among rotations and vibrations. As first elaborated by Podolsky,58 the position representation for the quantum-mechanical operator of kinetic energy, T̂ , can be recast in terms of curvilinear internal coordinates, {qj; j = 1, 2, ..., 3N − 3}, by means of T̂ = −

3N − 3

(A2)



3N − 3

ℏ2 2



ASSOCIATED CONTENT

S Supporting Information *

Structure for (S)-(+)-3-carene (Table S1); structure for (S)(+)-2-carene(A) (Table S2); structure for (S)-(+)-2-carene(B) (Table S3); structure for (S)-(+)-2-carene(TS) (Table S4); origin dependence for coupled-cluster predictions of specific rotation (Table S5); density-functional predictions for excited electronic states (Table S6); coupled-cluster predictions for excited electronic states (Table S7); and temperature dependence of chiroptical response in (S)-(+)-2-carene (Figure S1). This material is available free of charge via the Internet at http://pubs.acs.org.

(A1)

which is appropriate for an accompanying rovibrational wave function, ψ(q), normalized with respect to the volume element dq = dq1dq2···dq3N−3. The quantity gjk signifies an element of the contravariant mass-weighted metric tensor, g, that reflects the conversion between rectangular-Cartesian and curvilinear9530

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A



Article

(11) Schenck, G. O.; Schroeter, S.; Ohloff, G. Photosensitized Oxidation of (+)-3-Carene. Chem. Ind. (London, U.K.) 1962, No. 10, 459. Gollnick, K.; Schroeter, S.; Ohloff, G.; Schade, G.; Schenck, G. O. Zur Photosensibilisierten O2-Ü bertragung auf (+)-Caren-(3). Liebigs Ann. Chem. 1965, 687 (1), 14−25. Dassler, H.-G. Ü ber Δ3-CarenHydro-Peroxyd. Naturwissenschaften 1960, 47 (7), 157. (12) Naumov, V. A.; Bezzubov, V. M. Electron Diffraction Study of the Structure of the Δ3-Carene Molecule. Dokl. Akad. Nauk SSSR 1966, 171 (3), 634−636. (13) Acharya, S. P. Conformations of 3-carene and 2-carene. Their Conformational Preference and the Reactivity of the Double Bond. Tetrahedron Lett. 1966, 7 (34), 4117−4122. (14) Abraham, R. J.; Cooper, M. A.; Whittaker, D. The NMR Spectra and Conformations of Cyclic Compounds. VIII. Car-3-ene. Org. Magn. Reson. 1973, 5 (11), 515−516. (15) Favini, G.; Sottocornola, M.; Todeshchini, R. Conformation of Bicyclo[n.1.0] Derivatives. J. Mol. Struct.: THEOCHEM 1982, 90 (1− 2), 165−176. (16) Tylli, H. A Conformational Study of (+)-Car-3-ene: Fully Optimized STO-3G Structures. J. Mol. Struct.: THEOCHEM 1987, 151, 71−79. (17) Tkachev, A. V.; Denisov, A. Y. Is the Six-Membered Ring of the 3-Carene Molecule Planar? Mendeleev Commun. 1991, 1 (3), 98−99. (18) McCann, D. M.; Stephens, P. J. Determination of Absolute Configuration Using Density Functional Theory Calculations of Optical Rotation and Electronic Circular Dichroism: Chiral Alkenes. J. Org. Chem. 2006, 71 (16), 6074−6098. (19) Tylli, H. A Molecular Orbital Study of (+)-Car-2-ene: Fully Optimized STO-3G Structures. J. Mol. Struct: THEOCHEM 1988, 181 (1−2), 11−18. (20) Polavarapu, P. L.; Petrovic, A.; Wang, F. Intrinsic Rotation and Molecular Structure. Chirality 2003, 15 (S1), S143−S149. (21) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al.; Gaussian 09, Revision C.01; Gaussian, Inc.: Wallingford, CT, 2009. (22) Baboul, A. G.; Curtiss, L. A.; Redfern, P. C.; Raghavachari, K. Gaussian-3 Theory Using Density Functional Geometries and ZeroPoint Energies. J. Chem. Phys. 1999, 110 (16), 7650−7657. (23) Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Rassolov, V.; Pople, J. A. Gaussian-3 (G3) Theory for Molecules Containing First and Second-Row Atoms. J. Chem. Phys. 1998, 109 (18), 7764−7776. (24) Ochterski, J. W.; Petersson, G. A. A Complete Basis Set Model Chemistry. V. Extensions to Six or More Heavy Atoms. J. Chem. Phys. 1996, 104 (7), 2598−2619. (25) Curtiss, L. A.; Redfern, P. C.; Raghavachari, K. Gaussian-4 Theory. J. Chem. Phys. 2007, 126 (8), 084108/1−12. (26) Tomasi, J.; Mennucci, B.; Cancès, E. The IEF Version of the PCM Solvation Method: An Overview of a New Method Addressed to Study Molecular Solutes at the QM Ab Initio Level. J. Mol. Struct: THEOCHEM 1999, 464 (1−3), 211−226. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105 (8), 2999−3093. (27) Crawford, T. D.; Sherrill, C. D.; Valeev, E. F.; Fermann, J. T.; King, R. A.; Leininger, M. L.; Brown, S. T.; Janssen, C. L.; Seidl, E. T.; Kenny, J. P.; Allen, W. D. PSI3: An Open-Source Ab Initio Electronic Structure Package. J. Comput. Chem. 2007, 28 (9), 1610−1616. (28) Stanton, J. F.; Gauss, J.; Harding, M. E.; Szalay, P. G.; CFOUR, Coupled-Cluster Techniques for Computational Chemistry (June 2012 release; http://www.cfour.de); Harding, M. E.; Metzroth, T.; Gauss, J. Parallel Calculation of CCSD and CCSD(T) Analytic First and Second Derivatives. J. Chem. Theory Comput. 2008, 4 (1), 64−74. (29) London, F. Théorie Quantique des Courants Interatomiques dans les Combinaisons Aromatiques. J. Phys. Radium 1937, 8 (10), 397−409. Ditchfield, R. Self-Consistent Perturbation Theory of Diamagnetism. Mol. Phys. 1974, 27 (4), 789−807. Wolinski, K.; Hinton, J. F.; Pulay, P. Efficient Implementation of the GaugeIndependent Atomic Orbital Method for NMR Chemical-Shift Calculations. J. Am. Chem. Soc. 1990, 112 (23), 8251−8260.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Telephone: (203) 4323975. Fax: (203) 432-6144. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This study was performed under the auspices of grants CHE0809856 and CHE-1112239 awarded by the Chemical Structures, Dynamics, and Mechanisms Program in the Directorate for Mathematical and Physical Sciences of the United States National Science Foundation. The authors wish to thank Prof. T. Daniel Crawford (Virginia Tech) for technical assistance with the Psi3.4 quantum-chemistry package. Computational resources utilized in this work were supported in part by the Yale University Faculty of Arts and Sciences High Performance Computing Center and by the National Science Foundation under grant CNS-0821132, which partially funded the acquisition of requisite computer facilities.



REFERENCES

(1) Polavarapu, P. L. Renaissance in Chiroptical Spectroscopic Methods for Molecular Structure Determination. Chem. Rec. 2007, 7 (2), 125−136. (2) Pecul, M.; Ruud, K. The Ab Initio Calculation of Optical Rotation and Electronic Circular Dichroism. Adv. Quantum Chem. 2005, 50, 185−212. Crawford, T. D. Ab Initio Calculation of Molecular Chiroptical Properties. Theor. Chem. Acc. 2006, 115 (4), 227−245. Crawford, T. D.; Tam, M. C.; Abrams, M. L. The Current State of Ab Initio Calculations of Optical Rotation and Electronic Circular Dichroism Spectroscopy. J. Phys. Chem. A 2007, 111 (48), 12057− 12068. Autschbach, J. Computing Chiroptical Properties with FirstPrinciples Theoretical Methods: Background and Illustrative Examples. Chirality 2009, 21 (1E), E116−E152. (3) Vaccaro, P. H. Chapter 11: Optical Rotation and Intrinsic Optical Activity. In Comprehensive Chiroptical Spectroscopy (Vol. I: Instrumentation, Methodologies, and Theroretical Simulations), Berova, N., Polavarapu, P. L., Nakanishi, K., Woody, R. W., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, 2012; Vol. 1, pp 275−323. (4) Polavarapu, P. L. Optical Rotation: Recent Advances in Determining the Absolute Configuration. Chirality 2002, 14 (10), 768−781. Polavarapu, P. L. Erratum: Optical Rotation: Recent Advances in Determining the Absolute Configuration. Chirality 2002, 15(3), 284−285. (5) Barron, L. D. Molecular Light Scattering and Optical Activity, 2nd ed.; Cambridge University Press: Cambridge, UK, 2004. (6) Müller, T.; Wiberg, K. B.; Vaccaro, P. H. Cavity Ring-Down Polarimetry (CRDP): A New Scheme for Probing Circular Birefringence and Circular Dichroism in the Gas Phase. J. Phys. Chem. A 2000, 104 (25), 5959−5968. Müller, T.; Wiberg, K. B.; Vaccaro, P. H.; Cheeseman, J. R.; Frisch, M. J. Cavity Ring-Down Polarimetry (CRDP): Theoretical and Experimental Characterization. J. Opt. Soc. Am. B 2002, 19 (1), 125−141. (7) Wilson, S. M.; Wiberg, K. B.; Cheeseman, J. R.; Frisch, M. J.; Vaccaro, P. H. Nonresonant Optical Activity of Isolated Organic Molecules. J. Phys. Chem. A 2005, 109 (51), 11752−11764. (8) Ohloff, G.; Schulte-Elte, K. H.; Giersch, W. Die Ü berführung von (+)-Δ3 Caren in (+)-Δ4-Caren. Helv. Chim. Acta 1965, 48 (7), 1665− 1668. (9) Jadhav, P. K.; Prasad, J. V. N. V.; Brown, H. C. A Convenient Method for Upgrading the Enantiomeric Purities of (+)-Longifolene and (+)-3-Carene to Materials Approaching 100% ee. J. Org. Chem. 1985, 50 (17), 3203−3206. ́ (10) Vereshchagin, A. N.; Vulfson, S. G. Conformation of Δ(su3) Carene. J. Struct. Chem. 1968, 9 (6), 978−980. 9531

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

Helgaker, T.; Jørgensen, P. An Electronic Hamiltonian for Origin Independent Calculations of Magnetic Properties. J. Chem. Phys. 1991, 95 (4), 2595−2601. Cheeseman, J. R.; Frisch, M. J.; Devlin, F. J.; Stephens, P. J. Hartree-Fock and Density Functional Theory Ab Initio Calculation of Optical Rotation Using GIAOs: Basis Set Dependence. J. Phys. Chem. A 2000, 104 (5), 1039−1046. (30) Olsen, J.; Jørgensen, P. Linear and Nonlinear Response Functions for an Exact State and for an MCSCF State. J. Chem. Phys. 1985, 82 (7), 3235−3264. Grimme, S.; Furche, F.; Ahlrichs, R. An Improved Method for Density-Functional Calculations of the Frequency-Dependent Optical Rotation. Chem. Phys. Lett. 2002, 361 (3−4), 321−328. (31) Pedersen, T. B.; Koch, H.; Boman, L.; de Meras, A. M. J. S. Origin Invariant Calculation of Optical Rotation Without Recourse to London Orbitals. Chem. Phys. Lett. 2004, 393 (4−6), 319−326. (32) Wiberg, K. B.; Wang, Y.; Wilson, S. M.; Vaccaro, P. H.; Jorgensen, W. L.; Crawford, T. D.; Abrams, M. L.; Cheeseman, J. R.; Luderer, M. Optical Rotatory Dispersion of 2,3-Hexadiene and 2,3Pentadiene. J. Phys. Chem. A 2008, 112 (11), 2415−2422. (33) Ruud, K.; Taylor, P. R.; Åstrand, P.-O. Zero-Point Vibrational Effects on Optical Rotation. Chem. Phys. Lett. 2001, 337 (1−3), 217− 223. Ruud, K.; Zanasi, R. The Importance of Molecular Vibrations: The Sign Change of the Optical Rotation of Methyloxirane. Angew. Chem., Int. Ed. 2005, 44 (23), 3594−3596. (34) Mort, B. C.; Autschbach, J. Magnitude of Zero-Point Vibrational Corrections to Optical Rotation in Rigid Organic Molecules: A TimeDependent Density Functional Study. J. Phys. Chem. A 2005, 109 (38), 8617−8623. Mort, B. C.; Autschbach, J. Temperature Dependence of the Optical Rotation of Fenchone Calcuated by Vibrational Averaging. J. Phys. Chem. A 2006, 110 (40), 11381−11383. Mort, B. C.; Autschbach, J. Temperature Dependence of the Optical Rotation in Six Bicyclic Organic Molecules Calculated by Vibrational Averaging. Chem. Phys. Chem. 2007, 8 (4), 605−616. (35) Kongsted, J.; Pedersen, T. B.; Jensen, L.; Hansen, A. E.; Mikkelsen, K. V. Coupled Cluster and Density Functional Theory Studies of the Vibrational Contribution to the Optical Rotation of (S)Propylene Oxide. J. Am. Chem. Soc. 2006, 128 (3), 976−982. Pedersen, T. B.; Kongsted, J.; Crawford, T. D. Gas Phase Optical Rotation Calculated from Coupled Cluster Theory with Zero-Point Vibrational Corrections from Density Functional Theory. Chirality 2009, 21 (1E), E68−E75. Pedersen, T. B.; Kongsted, J.; Crawford, T. D.; Ruud, K. On the Importance of Vibrational Contributions to Small-Angle Optical Rotation: Fluoro-Oxirane in Gas Phase and Solution. J. Chem. Phys. 2009, 130 (3), 034310 /1−7. (36) Wiberg, K. B.; Vaccaro, P. H.; Cheeseman, J. R. Conformational Effects on Optical Rotation. 3-Substituted 1-Butenes. J. Am. Chem. Soc. 2003, 125 (7), 1888−1896. Wiberg, K. B.; Wang, Y. G.; Vaccaro, P. H.; Cheeseman, J. R.; Luderer, M. R. Conformational Effects on Optical Rotation. 2-Substituted Butanes. J. Phys. Chem. A 2005, 109 (15), 3405−3410. (37) Wilson, S. M.; Wiberg, K. B.; Murphy, M. J.; Vaccaro, P. H. The Effects of Conformation and Solvation on Optical Rotation: Substituted Epoxides. Chirality 2008, 20 (3−4), 357−369. (38) Wang, F.; Polavarapu, P. L. Conformational Stability of (+)-Epichlorohydrin. J. Phys. Chem. A 2000, 104 (26), 6189−6196. He, J.; Petrovic, A.; Polavarapu, P. L. Quantitative Determination of Conformer Populations: Assessment of Specific Rotation, Vibrational Absorption, and Vibrational Circular Dichroism in Substituted Butynes. J. Phys. Chem. A 2004, 108 (10), 1671−1680. He, J.; Petrovic, A. G.; Polavarapu, P. L. Determining the Conformer Populations of (R)-(+)-3-Methylcyclopentanone Using Vibrational Absorption, Vibrational Circular Dichroism, and Specific Rotation. J. Phys. Chem. A 2004, 108 (52), 20451−20457. Polavarapu, P. L.; Jeirath, N.; Walia, S. Conformational Sensitivity of Chiroptical Spectroscopic Methods: 6,6′-Dibromo-1,1′-Bi-2-Naphthol. J. Phys. Chem. A 2009, 113 (18), 5423−5431. (39) Tam, M. C.; Crawford, T. D. Ab Initio Determination of Optical Rotatory Dispersion in the Conformationally Flexible Molecule (R)Epichlorohydrin. J. Phys. Chem. A 2006, 110 (6), 2290−2298.

(40) Crawford, T. D.; Allen, W. D. Optical Activity in Conformationally Flexible Molecules: A Theoretical Study of Large-Amplitude Vibrational Averaging in (R)-3-Chloro-1-Butene. Mol. Phys. 2009, 107 (8−12), 1041−1057. (41) Grimme, S.; Bahlmann, A.; Haufe, G. Ab Initio Calculations for the Optical Rotations of Conformationally Flexible Molecules: A Case Study on Six-, Seven-, and Eight-Membered Fluorinated Cycloalkanol Esters. Chirality 2002, 14 (10), 793−797. Pecul, M.; Ruud, K.; Rizzo, A.; Helgaker, T. Conformational Effects on the Optical Rotation of Alanine and Proline. J. Phys. Chem. A 2004, 108 (19), 4269−4276. Marchesan, D.; Coriani, S.; Forzato, C.; Nitti, P.; Pitacco, G.; Ruud, K. Optical Rotation Calculations of a Highly Flexible Molecule: The Case of Paraconic Acid. J. Phys. Chem. A 2005, 109 (7), 1449−1453. (42) Wiberg, K. B.; Wang, Y.; Wilson, S. M.; Vaccaro, P. H.; Cheeseman, J. R. Sum-Over-States Calculation of the Specific Rotations of Some Substituted Oxiranes, Chloropropionitrile, Ethane, and Norbornenone. J. Phys. Chem. A 2006, 110 (51), 13995−14002. (43) Although rigorous in terms of contributions arising from the ground electronic potential energy surface, the averaging procedure of eq 2 does not account for effects attributable to the manifolds of vibrational eigenstates supported by electronically excited states. (44) Detailed analyses have shown that the C2C3C4C5 dihedral angle, ϕ ≡ τC2C3−C4−C5, maps isomorphically onto the intrinsic reaction coordinate (IRC) for the 2-carene isomerization process, while still retaining the ability to project the molecular framework pragmatically beyond the stable equilibrium points for conformers (A) and (B) [as required for quantization of the aperiodic potential well, V(ϕ)]. (45) Allen, W. D.; Bodi, A.; Szalay, V.; Császár, A. G. Adiabatic Approximations to Internal Rotation. J. Chem. Phys. 2006, 124 (22), 224310/1−9. (46) Castro, M. E.; Niño, A.; Muñoz-Caro, C. Gmat. A Software Tool for the Computation of the Rovibrational G Matrix. Comput. Phys. Commun. 2009, 180 (7), 1183−1187. Castro, M. E.; Niño, A.; Muñoz-Caro, C. Heuristic Computation of the Rovibrational G Matrix in Optimized Molecule-Fixed Axes. Gmat 2.1. Comput. Phys. Commun. 2010, 181 (8), 1471−1475. (47) Berova, N.; Nakanishi, K.; Woody, R. W. Circular Dichroism: Principles and Applications; John Wiley and Sons, Inc.: New York, 2000; Eliel, E. L.; Wilen, S. H.; Doyle, M. P. Chapter 12: Chiroptical Properties. In Basic Organic Stereochemistry; John Wiley and Sons, Inc.: New York, 2001; pp 534−609. (48) Kumata, Y.; Furukawa, J.; Fueno, T. Effect of Solvent on the Optical Rotation of Propylene Oxide. Bull. Chem. Soc. Jpn. 1970, 43 (12), 3920−3921. Fischer, A. T.; Compton, R. N.; Pagni, R. M. Solvent Effects on the Optical Rotation of (S)-(−)-α-Methylbenzylamine. J. Phys. Chem. A 2006, 110 (22), 7067−7071. Lambert, J.; Compton, R. N.; Crawford, T. D. The Optical Activity of Carvone: A Theoretical and Experimental Investigation. J. Chem. Phys. 2012, 136 (11), 114512/1−12. (49) Wiberg, K. B.; Wang, Y. G.; Vaccaro, P. H.; Cheeseman, J. R.; Trucks, G.; Frisch, M. J. Optical Activity of 1-Butene, Butane, and Related Hydrocarbons. J. Phys. Chem. A 2004, 108 (1), 32−38. (50) Cramer, C. J.; Truhlar, D. G. Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics. Chem. Rev. 1999, 99 (8), 2161−2200. Cramer, C. J.; Truhlar, D. G. A Universal Approach to Solvation Modeling. Acc. Chem. Res. 2008, 41 (6), 760−768. Mennucci, B.; Cappelli, C.; Cammi, R.; Tomasi, J. Modeling Solvent Effects on Chiroptical Properties. Chirality 2011, 23 (9), 717−729. (51) Mukhopadhyay, P.; Zuber, G.; Goldsmith, M.-R.; Wipf, P.; Beratan, D. N. Solvent Effect on Optical Rotation: A Case Study of Methyloxirane in Water. Chem. Phys. Chem 2006, 7 (12), 2483−2486. Mukhopadhyay, P.; Zuber, G.; Wipf, P.; Beratan, D. N. Contribution of a Solute’s Chiral Solvent Imprint to Optical Rotation. Angew. Chem., Int. Ed. 2007, 46 (34), 6450−6452. (52) Wiberg, K. B.; Wang, Y. G.; Wilson, S. M.; Vaccaro, P. H.; Cheeseman, J. R. Chiroptical Properties of 2-Chloropropionitrile. J. Phys. Chem. A 2005, 109 (15), 3448−3453. 9532

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533

The Journal of Physical Chemistry A

Article

(53) Stephens, P. J.; Devlin, F. J.; Cheeseman, J. R.; Frisch, M. J. Calculation of Optical Rotation Using Density Functional Theory. J. Phys. Chem. A 2001, 105 (22), 5356−5371. Mennucci, B.; Tomasi, J.; Cammi, R.; Cheeseman, J. R.; Frisch, M. J.; Devlin, F. J.; Gabriel, S.; Stephens, P. J. Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules. J. Phys. Chem. A 2002, 106 (25), 6102−6113. Stephens, P. J.; McCann, D. M.; Cheeseman, J. R.; Frisch, M. J. Determination of Absolute Configurations of Chiral Molecules Using Ab Initio Time-Dependent Density Functional Theory Calculations of Optical Rotation: How Reliable Are Absolute Configurations Obtained for Molecules with Small Rotations? Chirality 2005, 17 (S1), S52−S64. (54) Condon, E. U. Theories of Optical Rotary Power. Rev. Mod. Phys. 1937, 9 (4), 432−457. Kauzmann, W.; Walter, J.; Eyring, H. Theories of Optical Rotatory Power. Chem. Rev. 1940, 26 (3), 339− 407. (55) Onsager, L. Electric Moments of Molecules in Liquids. J. Am. Chem. Soc. 1936, 58 (8), 1486−1493. (56) Wong, M. W.; Frisch, M. J.; Wiberg, K. B. Solvent Effects. 1. The Mediation of Electrostatic Effects by Solvents. J. Am. Chem. Soc. 1991, 113 (13), 4776−4782. (57) Foresmann, J. B.; Keith, T. A.; Wiberg, K. B.; Snoonian, J.; Frisch, M. J. Solvent Effects. 5. Influence of Cavity Shape, Truncation of Electrostatics, and Electron Correlation on Ab Initio Reaction Field Calculations. J. Phys. Chem. 1996, 100 (40), 16098−16104. (58) Podolsky, B. Quantum-Mechanically Correct Form of Hamiltonian Function for Conservative Systems. Phys. Rev. 1928, 32 (5), 812−816. (59) Pickett, H. M. Vibration-Rotation Interactions and Choice of Rotating Axes for Polyatomic Molecules. J. Chem. Phys. 1972, 56 (4), 1715−1723. (60) Niño, A.; Muñoz-Caro, C. Computation of Kinetic Energy Constants for Large Range Internal Motions in Molecules. Comput. Chem. 1994, 18 (1), 27−32.

9533

dx.doi.org/10.1021/jp303270d | J. Phys. Chem. A 2012, 116, 9516−9533