DFT study on cooperativity in the interactions of hydrazoic acid clusters

Density functional theory at the B3LYP level with the 6‐311G** basis set is performed to calculate the systems consisting of up to four hydrazoic acid molecules. The dimers are found to exhibit cyclic and chain structures with N … H contacts at ca. 2.1–2.7 Å. However, there are only cyclic structures with N … H contacts at ca. 2.0–2.3 Å and 2.0–2.1 Å in the trimer and tetramer, respectively. Hydrogen bond distances in the trimer and tetramer are shorter than those in the cyclic dimer as a result of the stronger interaction between molecules. The contribution of cooperative effect to the interaction energy is significant. After the correction of the basis set superposition error and zero‐point energy, the binding energies are ?10.69, ?29.34, and ?54.26 kJ·mol^?1 for the most stable dimer, trimer, and tetramer, respectively. The calculated IR shifts for N? H stretching mode increase with the size of the cluster growths, reaching more than 200 cm^?1 in the tetramer. For the most stable clusters, the transition from the monomer to dimer, dimer to trimer, and trimer to tetramer involve changes of ?14.40, ?25.68, and ?31.88 kJ·mol^?1 for the enthalpies at 298.15 K and 1atm, respectively. We also perform Mulliken populations analysis and find the Mulliken populations on intermolecular N … H increasing in the sequence of the dimer, trimer, and tetramer. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 279–286, 2003     

An Ab Initio Study of Conformational Energies in Dioxo-, Trioxo-, and Tetroxo-Dithianes

Ab initio Hartree–Fock calculations at the HF/6-31G* level of theory for geometry optimization and the MP2/6-31G*//HF/6-31G* and B3LYP/6-311G(2df,p)//HF/6-31G* levels for a single point total energy calculation are reported for the important energy-minimum conformations of 1,1-dioxo-thiane (2), 1,1-dioxo-1,2-dithiane (3), 1,1-dioxo-1,3-dithiane (4), 1,1-dioxo-1,4-dithiane (5), 1,1,2-trioxo-1,2-dithiane (6), 1,1,3-trioxo-1,3-dithiane (7), 1,1,4-trioxo-1,4-dithiane (8), 1,1,2,2-tetroxo-1,2-dithiane (9), 1,1,3,3-tetroxo-1,3-dithiane (10), and 1,1,4,4-tetroxo-1,4-dithiane (11). According to the MP2/6-31G*//HF/6-31G* calculations, compound 5 is more stable than 3 and 4 by 7.8 and 8.9 kJ mol^?1, respectively. The axial geometries of 6 and 8 are more stable than the equatorial forms by 21.4 and 19.1 kJ mol^?1, respectively, but the equatorial form of 7 is 4.1 kJ mol^?1 more stable than the axial geometry. Compound 11 is more stable than 9 and 10 by 49.3 and 31.0 kJ mol^?1, respectively.     

MP2, DFT‐D,and PCM study of the HMB–TCNE complex: Thermodynamics,electric properties,and solvent effects

Geometry, thermodynamic, and electric properties of the π‐EDA complex between hexamethylbenzene (HMB) and tetracyanoethylene (TCNE) are investigated at the MP2/6‐31G* and, partly, DFT‐D/6‐31G* levels. Solvent effects on the properties are evaluated using the PCM model. Fully optimized HMB–TCNE geometry in gas phase is a stacking complex with an interplanar distance 2.87 × 10^?10 m and the corresponding BSSE corrected interaction energy is ?51.3 kJ mol^?1. As expected, the interplanar distance is much shorter in comparison with HF and DFT results. However the crystal structures of both (HMB)₂–TCNE and HMB–TCNE complexes have interplanar distances somewhat larger (3.18 and 3.28 × 10^?10 m, respectively) than our MP2 gas phase value. Our estimate of the distance in CCl₄ on the basis of PCM solvent effect study is also larger (3.06–3.16 × 10^?10 m). The calculated enthalpy, entropy, Gibbs energy, and equilibrium constant of HMB–TCNE complex formation in gas phase are: ΔH⁰ = ?61.59 kJ mol^?1, ΔS = ?143 J mol^?1 K^?1, ΔG⁰ = ?18.97 kJ mol^?1, and K = 2,100 dm³ mol^?1. Experimental data, however, measured in CCl₄ are significantly lower: ΔH⁰ = ?34 kJ mol^?1, ΔS = ?70.4 J mol^?1 K^?1, ΔG⁰ = ?13.01 kJ mol^?1, and K = 190 dm³ mol^?1. The differences are caused by solvation effects which stabilize more the isolated components than the complex. The total solvent destabilization of Gibbs energy of the complex relatively to that of components is equal to 5.9 kJ mol^?1 which is very close to our PCM value 6.5 kJ mol^?1. MP2/6‐31G* dipole moment and polarizabilities are in reasonable agreement with experiment (3.56 D versus 2.8 D for dipole moment). The difference here is due to solvent effect which enlarges interplanar distance and thus decreases dipole moment value. The MP2/6‐31G* study supplemented by DFT‐D parameterization for enthalpy calculation, and by the PCM approach to include solvent effect seems to be proper tools to elucidate the properties of π‐EDA complexes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Conformational Energies in Organic Six-Membered Cyclic Sulfoxides

Ab initio Hartree–Fock calculations at the HF/6?31 G* level of theory for geometry optimization and the MP2/6?31 G*//HF/6?31 G* and B3LYP/6-311G(2df,p)//HF/6?31 G* levels for a single point total energy calculation are reported for the important energy-minimum conformations of 1-oxo-thiane (1), 1-oxo-1,2-dithiane (2), 1-oxo-1,3-dithiane (3), 1-oxo-1,4-dithiane (4), 1,2-dioxo-1,2-dithiane (5), 1,3-dioxo-1,3-dithiane (6), and 1,4-dioxo-1,4-dithiane (7). According to the MP2/6-31G*//HF/6-31G* calculations, while the axial conformations of compounds 1, 2, and 4 are more stable than the equatorial forms by 6.0, 20.0, and 9.9 kJ mol^?1, respectively, the equatorial geometry of 3 is 3.0 kJ mol^?1 more stable than the axial form. The diaxial conformations of 5 and 7 are calculated to have similar energies, but the diaxial form of 6 is about 43 kJ mol^?1 less stable than that of 5 or 7.     

C?X⋅⋅⋅O Halogen Bonding: Interactions of Trifluoromethyl Halides with Dimethyl Ether

The formation of weakly bound molecular complexes between dimethyl ether (DME) and the trifluoromethyl halides CF₃Cl, CF₃Br and CF₃I dissolved in liquid argon and in liquid krypton is investigated, using Raman and FTIR spectroscopy. For all halides evidence is found for the formation of C? X???O halogen‐bonded 1:1 complexes. At higher concentrations of CF₃Br, a weak absorption due to a 1:2 complex is also observed. Using spectra recorded at temperatures between 87 and 125 K, the complexation enthalpies for the complexes are determined to be ?6.8(3) kJ mol^?1 (DME?CF₃Cl), ?10.2(1) kJ mol^?1 (DME?CF₃Br), ?15.5(1) kJ mol^?1 (DME?CF₃I), and ?17.8(5) kJ mol^?1 [DME(?CF₃Br)₂]. Structural and spectral information on the complexes is obtained from ab initio calculations at the MP2/ 6‐311++G(d,p) and MP2/6‐311++G(d,p)+LanL2DZ* levels. By applying Monte Carlo free energy perturbation calculations to account for the solvent influences, and statistical thermodynamics to estimate the zero‐point vibrational and thermal influences, the ab initio complexation energies are converted into complexation enthalpies for the solutions in liquid argon. The resulting values are compared with the experimental data deduced from the cryosolutions.     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Pi-pi interaction in pyridine

pi-pi Interaction in pyridine dimer and trimer has been investigated in different geometries and orientations at the ab initio (HF, MP2) and DFT (B3LYP) levels of theory using various basis sets (6-31G, 6-31G, 6-311++G) and corrected for basis set superposition error (BSSE). While the HF and DFT calculations show the pyridine dimer and the trimer to be unstable with respect to the monomer, the MP2 calculations show them to be clearly stable, thus emphasizing the need to include electron correlation while determining stacking interaction in such systems. The calculated MP2/6-311++G binding energy (100% BSSE corrected) of the parallel-sandwich, antiparallel-sandwich, parallel-displaced, antiparallel-displaced, T-up and T-down geometries for pyridine dimer are 1.53, 3.05, 2.39, 3.97, 1.91, 1.47 kcal/mol, respectively. The results show the antiparallel-displaced geometry to be the most stable. The binding energies for the trimer in parallel-sandwich, antiparallel-sandwich, and antiparallel-displaced geometry are found to be 3.18, 6.14, and 8.04 kcal/mol, respectively.     

An investigation on the reaction mechanism of the F2+Cl2→2ClF using the B3LYP method

The reaction mechanism of F₂+Cl₂→2ClF has been investigated with the density functional theory at the B3LYP/6‐311G* level. Six transition states have been found for the three possible reaction paths and verified by the normal mode vibrational and IRC analyses. Ab initio MP2/6‐311G* geometry optimizations and CCSD(T)/6‐311G(2df)//MP2/6‐311G* single‐point energy calculations have been performed for comparison. It is found that when the F₂ (or Cl₂) molecule decomposes into atoms first and then the F (or Cl) atom reacts with the molecule Cl₂ (or F₂) nearly along the molecular axis, the energy barrier is very low. The calculated energy barrier of F attacking Cl₂ is zero and that of Cl attacking F₂ is only 15.57 kJ?mol^?1 at the B3LYP level. However, the calculated dissociation energies of F₂ and Cl₂ are as high as 145.40 and 192.48 kJ?mol^?1, respectively. When the reaction proceeds through a bimolecular reaction mechanism, two four‐center transition states are obtained and the lower energy barrier is 218.69 kJ?mol^?1. Therefore, the title reaction F₂+Cl₂→2ClF is most probably initiated from the atomization of the F₂ molecule and terminated by the reaction of F attacking Cl₂ nearly along the Cl? Cl bond. MP2 calculations lead to the same conclusion, but the geometry of TS and the energy barrier are somewhat different. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2002     

Cooperative effects in water tetramers. Comparison between the empirical many-body model TCPE and ab initio calculations

Ab initio calculations at the MP2/6-311+G(2d,2p) and the empirical water many-body model TCPE have been applied to the study of four water tetramers corresponding to various molecular arrangements. For the cyclic tetramer (where each molecule is simultaneously donor and acceptor of hydrogen bonds, HBs), cooperative effects have been shown from ab initio computations to be stabilizing and to represent a contribution in the binding energy of 9 kcal mol⁻¹, while for the tetramer where only two molecules are simultaneously donor and acceptor of HB, such effects are stabilizing by only 1.5 kcal mol⁻¹. At last, for the tetramer where no molecule is simultaneously donor and acceptor of HBs, cooperative effects are smoothly destabilizing. TCPE predictions have been shown to be in good agreement with these ab initio estimates, both in terms of binding energy and cooperative effect contribution, which exhibits the accuracy of this potential.     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Dimethylnitroamine Clusters

Ab initio SCF and Mφller-Plesset correlation correction methods in combination with counterpose procedure for BSSE correction have been applied to the theroetical studying of dimethylnitroamine and its dimers and trimers.Three optimized stable dimers and two trimers have been obtained.The corrected binding energies of the most stable dimer and trimer were predicted to be -24.68kJ/mol and -47.27kJ/mol,respectively at the MP2/6-31G^*//HF/6-31G^* level.The proportion of correlated interation energies to their total interaction energies for all clusters was at least 29.3 percent,and the BSSE of ΔE(MP2) was at least 10.0kJ/mol.Dispersion and/or electrostatic force were dominant in all clusters.There exist cooperative effects in both the chain and the cyclic trimers.The vibrational frequencies associated with N-O stretches or wags exhibit slight red shifts,but the modes associated with the motion of hydrogen atoms of the methyl group show somewhat blue shifts with respect to those of monomer.Thermodynamic properties of dimethylnitroamine and its clusters at different temperatures have been calculated on the basis of vibrational analyses.The changes of the Gibbs free energies for the aggregation from monomer to the most stable dimer and trimer were predicted to be 14.37kJ/mol and 30.40kJ/mol,respectively,at 1 atm and 298.15K.     

Cooperative effects in extended hydrogen bonded systems involving O?H groups. Ab initio studies of the cyclic S4 water tetramer

The additional energy stabilization due to cooperative effects was calculated in extended hydrogen bonded systems O? H ?O? H ?O? H with unidirectional (homodromic) orientation of the O? H groups. Ab initio restricted Hartree Fock, MP2 and MP3 calculations with geometry optimization and BSSE correction have been performed using the GAUSSIAN 83 program package for the ground states of the linear water dimer with C_s symmetry and the cyclic water tetramer with S₄ symmetry. The latter represents the smallest possible, experimentally observed cooperative structure. A new definition for a cooperativity parameter is proposed. The definition is based on the two-body, non-neighbour interaction energy, plus three- and four-body contributions, including one-body deformation terms in relation to the total interaction energy of the water tetramer. The advantage of this definition is its independence of the reference system, which is necessary in complicated molecular systems with an undefined number of hydrogen bonds, such as disordered or flip-flop systems. According to this definition the energy gain based on cooperativity in the S₄ water tetramer is 29% with the MP3/6-31G** approximation, (30% with HF/4-31G* and 46% with HF/3-21G). The largest contribution of 18% is due to the three-body term on the MP3/6-31G** level, followed by the two-body, non-neighbour term with 11%. The four-body term and the deformation term are in the order of 1% and cancel each other because they have opposite sign.     

Proton and sodium ion affinities of glycine and its sodium salt in the gas phase. Ab initio calculations

The energies of protonation and Na⁺ cationization of glycine (GLY) and its (GLY ? H + Na) salt in the gas phase were calculated using ab inltio calculations. The proton affinity of GLY, valued at the MP2/6–31G*//3-21G level, is 937 kJ mol^?1. The amino function is confirmed to be the most favourable site of protonation: ‘proton affinities’ of the carbonyl and hydroxyl functions are calculated to be 75 and 180 kJ mol^?1, respectively, lower than that of NH₂ at the MP2/6-31G*//3–21G level. Calculations performed up to the MP2/6–31G*//3–21G level give the Na⁺ affinity of GLY as 189 kJ mol^?1 and the H⁺ and Na⁺ affinities of (GLY – H + Na) as 1079 and 298 kJ mol^?1, respectively. The geometries of all neutral and protonated species optimized with the 3–21G basis set are described. Both H* and Na⁺ cations complex preferably between the nitrogen atom and the carbonyl oxygen atom, leading to pseudo-five-membered ring structures in which Na? O and Na? N bonds lengths are greater than 2 Å.     

Density Functional Theory Study of Hydrogen Bonds of Bipyridine with 1,3,5-Benzenetricarboxylic Acid

The hydrogen‐bonded dimer and trimer formed between 1,3,5‐benzenetricarboxylic acid and bipyridine have been investigated using a density functional theory (DFT) method and 6‐31++G** basis set. The interaction energies are ?45.783 and ?89.998 kJ·mol^?1 for the most stable dimer and trimer, respectively, after the basis set superposition error and zero‐point corrections. The formation of O–H...N hydrogen bonds makes O–H symmetric stretching modes in the dimer and trimer red‐shifted relative to those of the 1,3,5‐benzenetricarboxylic acid monomer. The natural bond orbit analysis shows that the inter‐molecular charge transfers are 0.60475e and 1.20225e for the dimer and trimer, respectively. Thermodynamic analysis indicates that the formation of trimer is an exothermic and spontaneous process at low and room temperature. A supramolecule can be constructed through the strong N···H–O intermolecular hydrogen bonds between bipyridine and 1,3,5‐benzenetricarboxylic acid, which is in good agreement with the experimental results.     

Ab Initio Study of Ring Flipping of the Overcrowded Peri-Substituted Naphthalenes

Ab inintio molecular orbital and density functional theory method were used to investigate the structural and dynamic behavior of 1,8-di-tert-butyl naphthalene (1), 1,8-bis(trimethylsilyl)naphthalene (2), 1,8-bis(trimethylgermyl)naphthalene (3), and 1,8-bis(trimethylstannyl)naphthalene (4). HF/3-21G//HF/3-21G results revealed that the ring flipping barrier height of compound 1–4 is 92.59, 32.13, 26.76, and 15.46 kJ mol^?1 respectively. The obtained results show that the transition state structure for ring flipping of the bulky-groups is in a planar form with naphthalene ring. Contrary to compound 1, the ring flipping of compounds 2–4 occurred easily at room temperature. Also, MP2/3-21G//HF/3-21G energy calculation, show that the enantiomerization energy of compounds 1–4 are 97.99, 33.24, 26.80, and 15.38 kJ·mol^?1 respectively. The required energy for ring inversion of compounds 1–4 are 85.09, 27.26, 21.54, and 10.21 kJ mol^?1 respectively, as calculated by B3LYP/3-21G//HF/3-21G method. It can be concluded that the lower energy barrier of the ring flipping of compounds 2–4 is related to the increasing of the bond lengths of Si—C, Ge—C, and Sn—C, in contrast to C—C bond.     

Hydrogen Bonding in Phosphine Oxide/Phosphate–Phenol Complexes

To develop a new solvent‐impregnated resin (SIR) system for the removal of phenols and thiophenols from water, complex formation by hydrogen bonding of phosphine oxides and phosphates is studied using isothermal titration calorimetry (ITC) and quantum chemical modeling. Six different computational methods are used: B3LYP, M06‐2X, MP2, spin component‐scaled (SCS) MP2 [all four with 6‐311+G(d,p) basis set], a complete basis set extrapolation at the MP2 level (MP2/CBS), and the composite CBS‐Q model. This reveals a range of binding enthalpies (ΔH) for phenol–phosphine oxide and phenol–phosphate complexes and their thio analogues. Both structural (bond lengths/angles) and electronic elements (charges, bond orders) are studied. Furthermore, solvent effects are investigated theoretically by the PCM solvent model and experimentally via ITC. From our calculations, a trialkylphosphine oxide is found to be the most promising extractant for phenol in SIRs, yielding ΔH=?14.5 and ?9.8 kcal mol^?1 with phenol and thiophenol, respectively (MP2/CBS), without dimer formation that would hamper the phenol complexation. In ITC measurements, the ΔH of this complex was most negative in the noncoordinating solvent cyclohexane, and slightly less so in π–π interacting solvents such as benzene. The strongest binding is found for the dimethyl phosphate–phenol complex [?15.1 kcal mol^?1 (MP2/CBS)], due to the formation of two H‐bonds (P?O???H‐O‐ and P‐O‐H???O‐H); however, dimer formation of these phosphates competes with complexation of phenol, and would thus hamper their use in industrial extractions. CBS‐Q calculations display erroneous trends for sulfur compounds, and are found to be unsuitable. Computationally relatively cheap SCS‐MP2 and M06‐2X calculations did accurately agree with the much more elaborate MP2/CBS method, with an average deviation of less than 1 kcal mol^?1.     

A theoretical study on the intermolecular interaction of energetic system-Nitromethane dimer

Three optimized geometries of nitromethane dimer have been obtained at the HF/6-31G level.Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero point energy.Computed results indicate that the cyclic structure of (CH3NO2)2 is the most stable of three optimized geometries,whose corrected binding energyis 17.29 kJ mol-1 at the MP4SDTQ/6-31G//HF/6-31G level.In the optimized structures of nitromethane dimer,the inter-molecular hydrogen bond has not been found; and the charge-transfer interaction between CH3NO2 subsystems is weak; and the correlation interaction energy makes a little contribution to the intermolecular interaction energy of the dimer.     

CH3F分子间相互作用的ab initio研究

在MP2/6-311+ +G(3d,3p)电子相关校正水平上,对CH₃F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证,发现3种势能面上有极小点的构型.进一步在高级电子相关校正的MP4S-DTQ、CCSD(T)/6-311+ +G(3df,3pd)方法水平上,对其中总能量最小的构型进行精确计算,得到二聚体的结合能为-9.707kJ/mol.研究结果支持了由光谱实验结果推测的构型,解释了CH₃F二聚体的谐振频率的多样性.     

乙烯、乙炔和苯与氨基锂间的π型锂键相互作用

分别在DFT-B3LYP和MP2/6-311++G**水平上求得乙烯、乙炔和苯与氨基锂锂键复合物势能面上的3个几何构型. 频率分析表明,3个构型均为稳定构型. 计算结果表明,形成锂键复合物后,质子供体N(2)- Li(4)的键长明显增大,且其伸缩振动的频率发生了不同程度的红移. 分别在乙烯…氨基锂、乙炔…氨基锂和苯…氨基锂三种复合物中,经MP2/6-311++G**水平计算的同时含基组重叠误差(BSSE)和零点振动能校正的单体间锂键相互作用能分别为-26.04、-24.86 和 -30.02 kJ·mol-1. 自然键轨道理论（NBO）分析表明单体间的弱相互作用属于π-s型锂键.     

Dihydrogen trioxide clusters, (HOOOH)n (n = 2-4), and the hydrogen-bonded complexes of HOOOH with acetone and dimethyl ether: implications for the decomposition of HOOOH

Hydrogen-bonded gas-phase molecular clusters of dihydrogen trioxide (HOOOH) have been investigated using DFT (B3LYP/6-311++G(3df,3pd)) and MP2/6-311++G(3df,3pd) methods. The binding energies, vibrational frequencies, and dipole moments for the various dimer, trimer, and tetramer structures, in which HOOOH acts as a proton donor as well as an acceptor, are reported. The stronger binding interaction in the HOOOH dimer, as compared to that in the analogous cyclic structure of the HOOH dimer, indicates that dihydrogen trioxide is a stronger acid than hydrogen peroxide. A new decomposition pathway for HOOOH was explored. Decomposition occurs via an eight-membered ring transition state for the intermolecular (slightly asynchronous) transfer of two protons between the HOOOH molecules, which form a cyclic dimer, to produce water and singlet oxygen (Delta (1)O 2). This autocatalytic decomposition appears to explain a relatively fast decomposition (Delta H a(298K) = 19.9 kcal/mol, B3LYP/6-311+G(d,p)) of HOOOH in nonpolar (inert) solvents, which might even compete with the water-assisted decomposition of this simplest of polyoxides (Delta H a(298K) = 18.8 kcal/mol for (H 2O) 2-assisted decomposition) in more polar solvents. The formation of relatively strongly hydrogen-bonded complexes between HOOOH and organic oxygen bases, HOOOH-B (B = acetone and dimethyl ether), strongly retards the decomposition in these bases as solvents, most likely by preventing such a proton transfer.     

AB INITIO MOLECULAR ORBITAL STUDY OF 1,2-DITHIANE AND 1,2,4,5-TETRATHIANE

Ab initio calculations at HF/6-31+G? level of theory for geometry optimization, and MP2/6-31+G?//HF/6-31+G? and B3LYP/6-31+G?//HF/6-31+G? levels for a single-point total energy calculation, are reported for the chair and twist conformations of 1,2-dithiane (1), 3,3,6,6-tetramethyl-1,2-dithiane (2), 1,2,4,5-tetrathiane (3), and 3,3,6,6-tetramethyl-1,2,4,5-tetrathiane (4). The C₂ symmetric chair conformations of 1 and 2 are calculated to be 21.9 and 8.6 kJ mol^?1 more stable than the corresponding twist forms. The calculated energy barriers for chair-to-twist processes in 1 and 2 are 56.3 and 72.8 kJ mol^?1, respectively. The C_2h symmetric chair conformation of 3 is 10.7 kJ mol^?1 more stable than the twist form. Interconversion of these forms takes place via a C₂ symmetric transition state, which is 67.5 kJ mol^?1 less stable than 3-Chair. The D₂ symmetric twist-boat conformation of 4 is calculated to be 4.0 kJ mol^?1 more stable than the C_2h symmetric chair form. The calculated strain energy for twist to chair process is 61.1 kJ mol^?1.