Conformations of dimethoxymethane: matrix isolation infrared and ab initio studies

Conformations of dimethoxymethane (DMM) were studied using matrix isolation infrared spectroscopy. DMM was trapped in an argon matrix using an effusive source at 298, 388 and 430 K. Experiments were also done using a supersonic jet source to look for conformational cooling in the expansion process. As a result of these experiments, spectrally resolved infrared features of the ground and first higher energy conformer of DMM have been recorded, for the first time. The experimental studies were supported by ab initio computations performed at HF and B3LYP levels, using a 6-31++G** basis set. Computationally, four minima were identified corresponding to conformers with GG, TG, G+G- and TT structures. The computed frequencies at the B3LYP level were found to compare well with the experimental matrix isolation frequencies, leading to a definitive assignment of the infrared features of DMM, for the GG and TG conformers. At the B3LYP/6-31++G** level, the energy difference between the GG and TG conformers was computed to be 2.30 kcal mol(-1). The barrier for conformation interconversion, TG-->GG level was calculated to be 0.95 kcal mol(-1). This value is consistent with the experimental observation that the spectral features due to the TG conformer disappeared in the matrix on annealing.     

Temperature-induced dynamical conformational disorder in 4-vinyl benzoic acid molecular crystals: a molecular simulation study

Extensive molecular simulations are carried out as a function of temperature to understand and quantify the conformational disorder in molecular crystals of 4-vinyl benzoic acid. The conformational disorder is found to be dynamic and associated with a flip-flop motion of vinyl groups. The population of minor conformer is less than 3% up to 300 K and is 13.2% at 350 K and these results are consistent with the experimental observations. At still higher temperatures, the population of minor conformer increases up to 25%. The evolution of structure at both molecular and unit-cell level of the molecular crystal as a function of temperature has been characterized by various quantities such as radial distribution functions, average cell parameters, volume, and interaction energies. The van't Hoff plot shows a nonlinear behavior at lower temperatures as it has been reported recently by Ogawa and co-workers in the case of stilbene, azobenzene, and N-(4-methylbenzylidene)-4-methylaniline molecular crystals. A set of rigid body simulations were also carried out to quantify the effect of conformational disorder on structural quantities such as unit-cell volume and interaction energy. The anomalous shrinkage of vinyl C=C bond length as a function of temperature has been explained by combining the results of simulations and a set of constrained optimizations using ab initio electronic structure calculations for various molecular structures differing in torsional angle.     

Quantum simulation of a hydrated noradrenaline analog with the torsional path integral method

An extended version of the torsional path integral Monte Carlo (TPIMC) method is presented and shown to be useful for studying the conformation of flexible molecules in solvated clusters. The new technique is applied to the hydrated clusters of the 2-amino-1-phenyl-ethanol (APE) molecule. APE + nH2O clusters with n = 0-4 are studied at 100 and 300 K using both classical and quantum simulations. Only at the lower temperature is the hydration number n found to impact the conformational distribution of the APE molecule. This is shown to be a result of the temperature-dependent balance between the internal energy and entropy contributions to the relative conformer free energies. Furthermore, at 100 K, large quantum effects are observed in the calculated conformer populations. A particularly large quantum shift of 30% of the total population is calculated for the APE + 2H2O cluster, which is explained in terms of the relative zero point energy of the lowest-energy hydrated structures for this cluster. Finally, qualitative agreement is found between the reported calculations and recent spectroscopy experiments on the hydrated clusters of APE, including an entropically driven preference for the formation of AG-type hydrated structures and the formation of a water "droplet" in the APE + 4H2O cluster.     

Infrared photodissociation of a water molecule from a flexible molecule-H2O complex: rates and conformational product yields following XH stretch excitation

Infrared-ultraviolet hole-burning and hole-filling spectroscopies have been used to study IR-induced dissociation of the tryptamine.H2O and tryptamine.D2O complexes. Upon complexation of a single water molecule, the seven conformational isomers of tryptamine collapse to a single structure that retains the same ethylamine side chain conformation present in the most highly populated conformer of tryptamine monomer. Infrared excitation of the tryptamine.H2O complex was carried out using a series of infrared absorptions spanning the range of 2470-3715 cm-1. The authors have determined the conformational product yield over this range and the dissociation rate near threshold, where it is slow enough to be measured by our methods. The observed threshold for dissociation occurred at 2872 cm-1 in tryptamine.H2O and at 2869 cm-1 in tryptamine.D2O, with no dissociation occurring on the time scale of the experiment (approximately 2 micros) at 2745 cm-1. The dissociation time constants varied from approximately 200 ns for the 2869 cm-1 band of tryptamine.D2O to approximately 25 ns for the 2872 cm-1 band of tryptamine.H2O. This large isotope dependence is associated with a zero-point energy effect that increases the binding energy of the deuterated complex by approximately 190 cm-1, thereby reducing the excess energy available at the same excitation energy. At all higher energies, the dissociation lifetime was shorter than the pulse duration of our lasers (8 ns). At all wavelengths, the observed products in the presence of collisions are dominated by conformers A and B of tryptamine monomer, with small contributions from the other minor conformers. In addition, right at threshold (2869 cm-1), tryptamine.D2O dissociates exclusively to conformer A in the absence of collisions with helium, while both A and B conformational products are observed in the presence of collisions with helium. Using resolution-of-identity approximation to second-order Moller-Plesset binding energies extrapolated to the complete basis set limit and harmonic vibrational frequencies and transition states calculated at the density functional limit B3LYP/6-31+G* level of theory, Rice-Ramsperger-Kassel-Marcus (RRKM) predictions for the dissociation, isomerization, and water shuttling rates as a function of energy are made. At threshold, the experimental dissociation rate is almost 10(3) faster than RRKM predictions. Reasons for this apparent non-RRKM behavior will be discussed.     

An ab initio potential energy surface for the υ4 vibrational mode of hydrogen peroxide

In this paper, the levels and the torsional microstates of hydrogen peroxide are determined from fully optimized ab initio calculations using a nuclear model in one dimension. Calculations have been performed at the MP2 level with the 6-311 G(2df,2pd), 6-31 1+G(2df,2pd), cc-pVTZ and AUG-cc-pVTZ basis sets including polarization orbitals and diffuse functions. The most stable conformation, calculated with the MP2/AUG-cc-pVTZ approach, is a trans–gauche conformer lying at 67.5° from the trans structure. By using the same level of calculations, the heights of the trans and cis barriers have been determined to be 386.5 and 2643.8 cm⁻¹ in a good agreement with the experimental data. The variational torsional levels split into four components by the tunnelling effect of the barriers. The splitting of the fundamental level caused by the trans barrier has been found to be 11.8683 cm⁻¹, whereas the splitting caused by the cis barrier is insignificant under n=2. Current ab initio energies confirm the experimental assignments and verify the separability of the torsion from the rest of the vibrations. However, the experimental relation of dependence on the torsion of the rotational constants cannot be reproduced in one-dimension and depends on several additional vibrational effects.     

Matrix isolation infrared and ab initio study of the conformations of 2,2-dimethoxypropane

Conformations of 2,2-dimethoxypropane (DMP) were studied using matrix isolation infrared spectroscopy. An effusive source maintained at different temperatures (298, 388 and 430 K) was used to deposit DMP in a nitrogen matrix. As a result of these experiments, spectrally resolved infrared features of the ground and first higher energy conformer of DMP have been recorded, for the first time. The experimental studies were supported by ab initio computations performed at B3LYP/6-31++G** level. Computationally, four minima were identified corresponding to conformers with G+/-G-/+, TG+/-, G+/-G+/- and TT structures. The computed frequencies at the B3LYP level were found to compare well with the experimental matrix isolation frequencies, leading to a definitive assignment of the infrared features of DMP, for the G+/-G-/+ and TG+/- conformers. At the B3LYP/6-31++G** level, the energy difference between the G+/-G-/+ and TG+/- conformers was computed to be 3.25 kcal x mol(-1). The barrier for conformation interconversion, TG+/--->G+/-G-/+, was calculated to be 1.29 kcal x mol(-1). The magnitude of this barrier is consistent with the experimental observation that the spectral features due to the TG+/- decreased considerably in intensity when the matrix was annealed.     

Thermodynamic conformational analysis and structural stability of the nicotinic analgesic ABT-594

This work presents a theoretical study of the nicotinic analgesic ABT-594. We describe its neutral (precursor) and protonated (active) forms in vacuum and aqueous solution at the MP2/cc-pVDZ level. A conformational analysis is performed on the two torsional angles describing the orientation of the azetidinyl group and the azetidinylmethoxy moiety. To account for entropic effects, a thermostatistical study of conformational populations at physiological temperature is carried out. In the neutral form, conformer I is found as the most populated in vacuum and solution. Here, the nitrogen of the azetidinyl group is far from the electron pairs of the oxygen and the pyridinic nitrogen. In the protonated form, conformer VIII is the most stable in vacuum and solution. Now, the additional proton on the azetidinyl group is oriented toward the electron lone pairs of oxygen. The structural stability of conformers I and VIII is considered through the atoms in molecules theory. The conformer I, in the neutral forms, is stabilized by an intramolecular hydrogen bond. The preference of conformer VIII in the protonated forms is explained by the higher strength of its intramolecular hydrogen bond over the cation-pi interaction found in conformer I. The effect of the interaction energy with the receptor on the conformational preferences of protonated ABT-594 is simulated. The result is that the population of conformers associated to the rotation of the azetidinyl group increases. So, the molecule can easily adopt the optimal internitrogen separation for interaction with the receptor.     

Photolytic Decarbonylation of Oxalyl Diisothiocyanate in Solid Argon Matrices to syn-anti Carbonyl Diisothiocyanate and Its Isomerization

Oxalyl diisothiocyanate, ((CO)NCS)₂, has been studied in solid argon matrices at 4.2 K with the aid of infrared (IR) spectroscopy. The spectra show mainly signals attributed to the most stable anti-anti conformer, which is corroborated by comparison to computed anharmonic fundamental IR transitions. Upon irradiation with 254 nm UV light, oxalyl diisothiocyanate eliminates carbon monoxide under formation of carbonyl diisothiocyanate, CO(NCS)₂. This reaction is only slightly exothermic by 0.4 kcal mol⁻¹ at the DLPNO-CCSD(T)/def2-QZVPP//B3LYP-D3/def2-TZVPP level of theory. Remarkably, photolysis produces mostly the less stable syn-anti conformer of carbonyl diisothiocyanate. Subsequent annealing at 30 K for two minutes results in a structural relaxation to the 0.7 kcal mol⁻¹ more stable syn-syn conformer confirming a low torsional barrier height between the isomers.     

Ab initio determination of optical rotatory dispersion in the conformationally flexible molecule (R)-epichlorohydrin

Ab initio optical rotation data from linear-response coupled-cluster and density-functional methods are compared to both gas-phase and liquid-phase polarimetry data for the small, conformationally flexible molecule epichlorohydrin. Three energy minima exist along the C-C-C-Cl dihedral angle, each with strong, antagonistic specific rotations ranging from ca. -450 to +500 deg/[dm (g/mL)] at 355 nm. Density-functional theory (specifically the B3LYP functional) consistently overestimates the optical rotations of each conformer relative to coupled-cluster theory (in agreement with our earlier observations for conformationally rigid species), and we attribute this to density-functional theory's underestimation of the lowest-lying excitation energies of epichlorohydrin. Length- and velocity-gauge formulations of the coupled-cluster response function lead to slightly different specific rotations (ca. 7% at short wavelengths). We have determined well-converged Gibbs free energy differences among the conformers using complete-basis-set extrapolations of coupled-cluster energies including triple excitations to obtain Boltzmann-averaged specific rotations for comparison to the gas-phase results. The length-gauge coupled-cluster data agree remarkably well with experiment, with the velocity-gauge coupled-cluster and density-functional data bracketing the experimental results from below and above, respectively. Liquid-phase conformer populations reported earlier by Polavarapu and co-workers from combined infrared absorption and theoretical analyses differ markedly from the gas-phase populations, particularly for polar solvents. Nevertheless, Boltzmann-averaged specific rotations from both coupled-cluster and density-functional calculations agree well with the corresponding experimental intrinsic rotations, although the theoretical specific rotations for the individual conformers do not take solvent effects into account. PCM-based estimates of conformer populations lead to poor agreement with experiment.     

Matrix-isolation study and ab initio calculations of the structure and spectra of hydroxyacetone

The structure of hydroxyacetone (HA) isolated in an argon matrix (at 12 K) and in a neat solid phase (at 12-175 K) was characterized by using infrared (IR) spectroscopy. The interpretation of the experimental results was supported by high-level quantum chemical calculations, undertaken by using both ab initio (MP2) and density functional theory methods. A potential-energy surface scan, carried out at the MP2/6-311++G(d,p) level of theory, predicted four nonequivalent minima, Cc, Tt, Tg, and Ct, all of them doubly degenerate by symmetry. The energy barriers for conversion between most of the symmetrically related structures and also between some of the nonequivalent minima (e.g., Tg --> Tt and Ct --> Tt) are very small and stay below the zero-point vibrational level associated with the isomerization coordinate in the higher-energy form in each pair. Therefore, only Cc and Tt conformers have physical significance, with populations of 99 and 1%, respectively, in gas phase at room temperature. For the matrix-isolated compound, only the most stable Cc conformer was observed. On the other hand, the polarizable continuum model calculations indicated that in water solution, the population of Tt and Ct conformers might be high enough (ca. 6 and 11%, respectively) to enable their experimental detection, thus supporting the conclusions of a previous IR spectroscopy study [ Spectrochim. Acta A 2005, 61, 477] in which the presence of more than one HA conformer in aqueous solution was postulated. The signatures of these minor conformers, however, do not appear in the spectra of the neat HA crystal, and the crystal structure was rationalized in terms of centrosymmetric hydrogen-bonded dimers consisting of two Cc-like units. Finally, we calculated (1)H, (13)C, and (17)O NMR chemical shifts at different levels of theory and found them to agree with available experimental data.     

Impact of molecular conformation on barriers to internal methyl rotation: the rotational spectrum of m-methylbenzaldehyde

The ground state spectrum of m-methylbenzaldehyde (m-MBA) was measured with a chirped-pulse Fourier transform microwave (CP-FTMW) spectrometer. The methyl rotor on m-MBA introduces an internal rotation barrier, which leads to splitting of the torsional energy level degeneracy into A and E states. Ab initio calculations predict a low torsional barrier for both the O-cis and O-trans conformers, resulting in a large doublet splitting up to several gigahertz in the frequency spectrum. The rotational constants, distortion terms, and V(3) values for both species have been determined from the ground state rotational spectrum using the BELGI-C(s) fitting program. There are significant differences in the torsional potential for the O-cis and O-trans m-MBA conformers. Molecular orbitals and resonance structures for each conformer are analyzed to understand the difference in torsional barrier height as well as the irregular shape of the O-trans torsional potential.     

Evidence of conformational equilibrium of 1-ethyl-3-methylimidazolium in its ionic liquid salts: Raman spectroscopic study and quantum chemical calculations

Raman spectra of liquid 1-ethyl-3-methylimidazolium (EMI+) salts, EMI(+)BF4-, EMI(+)PF6-, EMI(+)CF3SO3-, and EMI(+)N(CF3SO2)2-, were measured over the frequency range 200-1600 cm(-1). In the range 200-500 cm(-1), we found five bands originating from the EMI+ ion at 241, 297, 387, 430, and 448 cm(-1). However, the 448 cm(-1) band could hardly be reproduced by theoretical calculations in terms of a given EMI+ conformer, implying that the band originates from another conformer. This is expected because the EMI+ involves an ethyl group bound to the N atom of the imidazolium ring, and the ethyl group can rotate along the C-N bond to yield conformers. The torsion energy for the rotation was then theoretically calculated. Two local minima with an energy difference of ca. 2 kJ mol(-1) were found, suggesting that two conformers are present in equilibrium. Full geometry optimizations followed by normal frequency analyses indicate that the two conformers are those with planar and nonplanar ethyl groups against the imidazolium ring plane, and the nonplanar conformer is favorable. It elucidates that bands at 241, 297, 387, and 430 cm(-1) mainly originate from the nonplanar conformer, whereas the 448 cm(-1) band does originate from the planar conformer. Indeed, the enthalpy for conformational change from nonplanar to planar EMI+ experimentally obtained by analyzing band intensities of the conformers at varying temperatures is practically the same as that evaluated by theoretical calculations. We thus conclude that the EMI+ ion exists as either a nonplanar or planar conformer in equilibrium in its liquid salts.     

Direct measurement of the energy thresholds to conformational isomerization in tryptamine: experiment and theory

The methods of stimulated emission pumping-hole filling spectroscopy (SEP-HFS) and stimulated emission pumping population transfer spectroscopy (SEP-PTS) were applied to the conformation-specific study of conformational isomerization in tryptamine [TRA, 3-(2-aminoethyl)indole]. These experimental methods employ stimulated emission pumping to selectively excite a fraction of the population of a single conformation of TRA to well-defined ground-state vibrational levels. This produces single conformations with well-defined internal energy, tunable over a range of energies from near the zero-point level to well above the lowest barriers to conformational isomerization. When the SEP step overcomes a barrier to isomerization, a fraction of the excited population isomerizes to form that product. By carrying out SEP excitation early in a supersonic expansion, these product molecules are subsequently cooled to their zero-point vibrational levels, where they can be detected downstream with a third tunable laser that probes the ground-state population of a particular product conformer via a unique ultraviolet transition using laser-induced fluorescence. The population transfer spectra (recorded by tuning the SEP dump laser while holding the pump and probe lasers fixed) exhibit sharp onsets that directly determine the energy thresholds for conformational isomerization in a given reactant-product conformer pair. In the absence of tunneling effects, the first observed transition in a given X-Y PTS constitutes an upper bound to the energy barrier to conformational isomerization, while the last transition not observed constitutes a lower bound. The bounds for isomerizing conformer A of tryptamine to B(688-748 cm(-1)), C(1)(860-1000 cm(-1)), C(2)(1219-1316 cm(-1)), D(1219-1282 cm(-1)), E(1219-1316 cm(-1)), and F(688-748 cm(-1)) are determined. In addition, thresholds for isomerizing from B to A(<1562 cm(-1)), B to F(562-688 cm(-1)), and out of C(2) to B(<747 cm(-1)) are also determined. The A-->B and B-->A transitions are used to place bounds on the relative energies of minima B relative to A, with B lying at least 126 cm(-1) above A. The corresponding barriers have been computed using both density functional and second-order many-body perturbation theory methods in order to establish the level of theory needed to reproduce experimental results. While most of the computed barriers match experiment well, the barriers for the A-F and B-F transitions are too high by almost a factor of 2. Possible reasons for this discrepancy are discussed.     

Molecular structure and conformation of chloronitromethane as determined by gas-phase electron diffraction and theoretical calculations

The molecular structure of chloronitromethane was studied in the gas phase at a nozzle-tip temperature of 373 K. The experimental data were interpreted using a dynamic model where the molecules are undergoing torsional motion governed by a potential function: V = V2/2x(1 - cos 2tau) + V4/2x(1 - cos 4tau) with V2 = 0.81(30) and V4 = 0.12(40) kcal/mol (tau is the dihedral angle between the C-Cl and N-O bond). The conformer with a zero degree dihedral angle is the most stable conformer. Comparison with results from HF/MP2/B3LYP 6-311G(d,p) calculations were made. The important geometrical parameter values (for the eclipsed form) obtained from least-squares refinements are the following: r(C-H) = 1.061(18)A, r(C-N) = 1.509 (5)A, r(N-O) = 1.223(1)A, r(C-Cl) = 1.742(2)A, angleClCN = 115.2(7) degrees, angleO4NC = 118.9(10) degrees, angleO5NC = 114.9(16) degrees, and angleClCH 115(4) degrees.     

Vibrational and theoretical studies of N-(fluorosulfonyl)imidosulfurous difluoride, FSO2N=SF2

The Raman spectrum of liquid N-(fluorosulfonyl)imidosulfurous difluoride FSO2N = SF2 and the IR spectrum of its vapour phase were recorded. The observed features in combination with the theoretical studies indicate the existence at room temperature of only one conformer. In accord with previous results obtained by electron-diffraction analysis, this single conformer possesses C1 symmetry in which the SF2 group is oriented syn with respect to the N-S single bond. Theoretical vibrational spectra were also determined using ab initio and density functional theory (DFT) calculations at different levels of approximation. For all except one of the torsional modes, experimental wavenumbers were obtained. A subsequent normal coordinate analysis was performed using a torsional wavenumber calculated by theoretical methods.     

Conformations of trimethoxymethylsilane: matrix isolation infrared and ab initio studies

Conformations of trimethoxymethylsilane were studied using matrix isolation infrared spectroscopy and ab initio computations. Trimethoxymethylsilane was trapped in both argon and nitrogen matrixes using heated nozzle effusive sources and a supersonic jet source, in an effort to alter the conformational population in the matrix. Ab initio calculations were carried out at the HF and B3LYP level using 6-31++G basis set to support our experimental observations. The frequencies computed at the B3LYP level was found to fit well with our experimental data. A conformer with a C1(g(+/-)g(+/-)t) structure was predicted by our computations to be the ground state conformer.     

Structural selection by microsolvation: conformational locking of tryptamine

The conformational space of tryptamine has been thoroughly investigated using rotationally resolved laser-induced fluorescence spectroscopy. Six conformers could be identified on the basis of the inertial parameters of several deuterated isotopomers. Upon attaching a single water molecule, the conformational space collapses into a single conformer. For the hydrogen-bonded water cluster, this conformer is identified unambiguously as tryptamine A. In the complex, the water molecule acts as proton donor with respect to the amino group. An additional interaction with one of the aromatic C-H bonds selectively stabilizes the observed conformer more than all other conformers. Ab initio calculations confirm much larger energy differences between the conformers of the water complex than between those of the monomers.     

Hammett correlations in the chemistry of 3-phenylpropyl radicals

The energetics and kinetics of the reaction of variously substituted benzyl radicals with a model alkene were calculated at the G3(MP2)-RAD//B3-LYP/6-31G(d) level of theory to determine whether such reactions are amenable to Hammett analysis. The reactions were studied both in the gas phase and in toluene solution in the temperature range 298-353 K; calculations include 1D-hindered rotor corrections for low frequency torsional modes, and the solvation energies were calculated using COSMO-RS at the BP/TZP level of theory. The addition reaction was found to be dominated by radical stabilization effects, but under circumstances where olefin substituent effects were decoupled from aryl substituent effects, a modest polar effect comes into play, which is enhanced by solvation. Reasonable correlations with empirical substituent parameters such as Hammett σ and σ(?) are observed for the enthalpy of activation, but additional entropic factors act to decrease the degree of correlation with respect to free energies and rate coefficients, confirming hypotheses from earlier experimental work. Substituent effects on the reverse β-fragmentation reaction, and potential cyclization of the 3-phenylpropyl radicals formed by addition are also discussed.     

High resolution millimeter-wave spectroscopy of vinyltellurol

The millimeter-wave rotational spectrum of vinyltellurol has been recorded and assigned for the first time. To support the spectrum assignment, high level ab initio calculations have been carried out. Geometries, total electronic energies, and harmonic vibrational frequencies have been determined at the MP2 level. A small-core relativistic pseudopotential basis set (cc-pVTZ-PP) was employed to describe the tellurium atom. Two stable conformers, synperiplanar (sp) and anticlinal (ac), have been identified. The sp conformer is planar with a small negative inertia defect of -0.025 u ?(2). The ac conformer was found to be nonplanar with a C-C-Te-H dihedral angle of about 140° from sp. This conformer exhibits a large amplitude motion associated with the torsion about the C-Te bond. The barrier to internal rotation is about 1 kJ/mol, according to the theoretical calculations. For the ac conformation, a torsional potential function consisting of quartic and quadratic terms of the torsional angle has been partially determined from the observed rotational constants.