Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Theoretical Study on the Tautomerization Reaction of Formamide to Formamidic Acid Encapsulated in a Nano-sized Molecular Container

The hybrid calculations with ONIOM(B3LYP/6-31G*:AM1) method were carried out on the tautomerization reaction of formamide to formamidic acid in the microcontainer-encapsulated state. The free-state tautomerization process was also investigated with B3LYP/6-31G**//B3LYP/6-31G* method for the purpose of comparison. Bare tautomerization, H2O-assisted(single-H2O or multiple-H2O) and self-assisted mechanisms were all taken into consideration for the encapsulated state. The results show that only bare tautomerization and single-H2O catalysis mechanisms are possible to the encapsulated for-mamide tautomerization owing to the container's size confinement effect. Geometrical changes in the complexed container and guest molecules are discussed to deeply understand the complex's structural properties. The bare tautomerization barrier in the encapsulated state increases by 23.826 kJ/mol, ac-counting for 12% of the corresponding total energy barrier in the free state, and the increased values for the single-H2O catalysis process are 12.958 kJ/mol, accounting for 16% of the corresponding total energy barrier, respectively. This finding suggests that the encapsulation can make the tautomerization process slightly difficult.     

Conformational study of the structure of free 18-crown-6

A conformational search was performed for 18-crown-6 using the CONLEX method at the MM3 level. To have a more accurate energy order of the predicted conformations, the predicted conformations were geometry optimized at the HF/STO-3G level and the 198 lowest energy conformations, according to the HF/STO-3G energy order, were geometry optimized at the HF/6-31+G level. In addition, the 47 nonredundant lowest energy conformations, according to the MP2/6-31+G energy order at the HF/6-31+G optimized geometry, hereafter the MP2/6-31+G//HF/6-31+G energy order, were geometry optimized at the B3LYP/6-31+G level. According to the MP2/6-31+G//B3LYP/6-31+G energy order, three conformations had energies lower than the experimentally known Ci conformation of 18c6. At the MP2/6-31+G//B3LYP/6-31+G level, the S6 lowest energy conformation is more stable by 1.96 kcal/mol than this Ci conformation. This was confirmed by results at the MP2/6-31+G level with an energy difference of 1.84 kcal/mol. Comparison between the structure of the S6 conformation of 18c6 and the S4 lowest energy conformation of 12-crown-4, as well as other important conformations of both molecules, is made. It is concluded that the correlation energy is necessary to have an accurate energy order of the predicted conformations. A rationalization of the conformational energy order in terms of the hydrogen bonding and conformational dihedral angles is given. It is also suggested that to have a better energy order of the predicted conformations at the MM3 level, better empirical force fields corresponding to the hydrogen bond interactions are needed.     

An Ab Initio Molecular Orbital Study of Isophosphinines

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single-point total energy calculation are reported for phosphinine and 13 isophosphinines 7-19 . Isomers 7-11 with an allenic system are calculated to be 8-18 kcal mol m 1 more stable than structures 12-17 with an acetylenic moiety. The calculated energy difference (66.19 kcal mol m 1 ) between phosphinine and the most stable isophosphinine (1-phospha-1,2,4-cyclohexatriene, 10 ) is smaller than the difference (78.96 kcal mol m 1 ) between benzene and the most stable isobenzene (cyclohexa-1,2,4-triene, 2 ). The isophosphinines 18 and 19 , with a butatriene moiety, are calculated to be the least stable isomers.     

Dimers of and tautomerism between 2-pyrimidinethiol and 2(1H)-pyrimidinethione: a density functional theory (DFT) study

Density functional theory (BLYP, B3LYP, B3P86, B3PW91) with the 6-31+G(d,p), 6-311+G(d,p), and cc-pVTZ basis sets has been used to calculate structural parameters, relative energies, and vibrational spectra of 2-pyrimidinethiol (1) and 2(1H)-pyrimidinethione (2) and their hydrogen-bonded homodimers (C(2) 3, C(2h) [4](double dagger), C(2h) 5), monohydrates, and dihydrates and a heterodimer (6). Several transition state structures proposed for the tautomerization process have also been examined. At the B3PW91/6-311+G(d,p)//B3PW91/6-31+G(d,p) level of theory 2-pyrimidinethiol (1) is predicted to be 3.41 kcal/mol more stable (E(rel)) than 2(1H)-pyrimidinethione (2) in the gas phase and 2 is predicted to be 6.47 kcal/mol more stable than 1 in aqueous medium. An unfavorable planar intramolecular strained four center transition state (TS1) for the tautomerization of 1 and 2 in the gas-phase lies 29.07 kcal/mol higher in energy than 2-pyrimidinethiol (1). The C(2) 2-pyrimidinethiol dimer (3) is 6.84 kcal/mol lower in energy than the C(2) homodimer transition state structure ([11](double dagger)) that connects dimers 3 and 4. Transition state [11](double dagger) provides a facile pathway for tautomerization between 1 and 2 in the gas phase (monomer-dimer promoted tautomerization). The hydrogen bonded 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O structures are predicted to be 1.27 and 1.55 kcal/mol, respectively, higher in energy than 2(1H)-pyrimidinethione- - -H(2)O and 2(1H)-pyrimidinethione- - -2H(2)O. Water promoted tautomerization via cyclic transition states involving one water molecule (TS- - -H(2)O, [12](double dagger)) and two water molecules (TS- - -2H(2)O, [13](double dagger)) lie 11.42 and 11.44 kcal/mol, respectively, higher in energy than 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O. Thus, the hydrated transition states [12](double dagger) and [13](double dagger) are involved in the tautomerism between 1 and 2 in aqueous medium.     

Computational studies of the structure and vibrational spectra of hexafluorodisilane

Ab initio calculations predict that D_3d symmetry of Si₂F₆ is more stable than D_3h symmetry. The calculated potential barrier to internal rotation was 0.77, 0.73 and 0.78 kcal/mol using HF/6-31G*, B3LYP/6-31G* and MP2/6-31G* methods respectively, which was in good agreement with the experimental value between 0.51±0.10 and 0.73±0.14 kcal/mol. The optimized geometries, harmonic force fields, infrared intensities, Raman activities, and vibrational frequencies are reported for D_3d symmetry of Si₂F₆ from HF/6-31G* and B3LYP/6-31G*. A normal coordinate analysis was carried out. The average error between the scaled DFT frequencies obtained from the B3LYP/6-31G* calculation and observed frequencies was 4.2 cm⁻¹ and the average error between the scaled HF and observed frequencies was 2.2 cm⁻¹.     

Heat of hydrogenation of 1,5-dehydroquadricyclane. A computational and experimental study of a highly pyramidalized alkene

The radical anion of the highly pyramidalized alkene 1,5-dehydroquadricyclane (1) was generated in the gas phase from the Squires reaction of 1,5-bis(trimethylsilyl)quadricyclane with F-/F2. The electron binding energy and proton affinity of 1*- were determined by bracketing experiments to be 0.6 +/- 0.1 eV and 386 +/- 5 kcal/mol, respectively. These values are in good agreement with values predicted by density functional theory (B3LYP/6-31+G*) and ab initio (CASPT2/6-31+G*) calculations. The experimental heat of hydrogenation of 1, obtained from a thermochemical cycle, was found to be 91 +/- 9 kcal/mol. This value of deltaH(H2) leads to values of 67 +/- 9 kcal/mol for the olefin strain energy (OSE) of 1, 172 +/- 9 kcal/mol for its heat of formation, and 23 +/- 9 kcal/mol for its pi bond dissociation enthalpy. Since the retro-Diels-Alder reaction of neutral 1 is computed to be highly exothermic, the finding that 1*- apparently does not undergo a retro-Diels-Alder reaction is of particular interest. The B3LYP/6-31+G* optimized geometry of 1 suggests that the bonding in this alkene is partially delocalized, presumably because the highly pyramidalized double bond in 1 interacts with the distal cyclopropane bonds in a manner that eventually leads to a retro-Diels-Alder reaction. The good agreement of the B3LYP and (2/2)CASPT2 values for the heat of hydrogenation and OSE of 1 with the experimentally derived values provides indirect evidence for the correctness of the B3LYP prediction that the equilibrium geometry of 1 lies part way along the reaction coordinate to the transition structure for the retro-Diels-Alder reaction.     

2-Pyridinethiol/2-pyridinethione tautomeric equilibrium. A comparative experimental and computational study

The gas phase and solvent dependent preference of the tautomerization between 2-pyridinethiol (2SH) and 2-pyridinethione (2S) has been assessed using variable temperature Fourier transform infrared (FTIR) experiments, as well as ab initio and density functional theory computations. No spectroscopic evidence (nu(S)(-)(H) stretch) for 2SH was observed in toluene, C(6)D(6), heptane, or methylene chloride solutions. Although, C(s)() 2SH is 2.61 kcal/mol more stable than C(s)() 2S (CCSD(T)/cc-pVTZ//B3LYP/6-311+G(3df,2p)+ZPE), cyclohexane solvent-field relative energies (IPCM-MP2/6-311+G(3df,2p)) favor 2S by 1.96 kcal/mol. This is in accord with the FTIR observations and in quantitative agreement with the -2.6 kcal/mol solution (toluene or C(6)D(6)) calorimetric enthalpy for the 2S/2SH tautomerization favoring the thione. As the intramolecular transition state for the 2S, 2SH tautomerization (2TS) lies 25 (CBS-Q) to 30 kcal/mol (CCSD/cc-pVTZ) higher in energy than either tautomer, tautomerization probably occurs in the hydrogen bonded dimer. The B3LYP/6-311+G(3df,2p) optimized C(2) 2SH dimer is 10.23 kcal/mol + ZPE higher in energy than the C(2)(h)() 2S dimer and is only 2.95 kcal/mol + ZPE lower in energy than the C(2) 2TS dimer transition state. Dimerization equilibrium measurements (FTIR, C(6)D(6)) over the temperature range 22-63 degrees C agree: K(eq)(298) = 165 +/- 40 M(-)(1), DeltaH = -7.0 +/- 0.7 kcal/mol, and DeltaS = -13.4 +/- 3.0 cal/(mol deg). The difference between experimental and B3LYP/6-311+G(3df,2p) [-34.62 cal/(mol deg)] entropy changes is due to solvent effects. The B3LYP/6-311+G(3df,2p) nucleus independent chemical shifts (NICS) are -8.8 and -3.5 ppm 1 A above the 2SH and 2S ring centers, respectively, and the thiol is aromatic. Although the thione is not aromatic, it is stabilized by the thioamide resonance. In solvent, the large 2S dipole, 2-3 times greater than 2SH, favors the thione tautomer and, in conclusion, 2S is thermodynamically more stable than 2SH in solution.     

An investigation of keto–enol tautomerism in N-vinyl-2-, and 3-pyrrolidinones using ab initio and density functional theory calculations

We investigated the protomeric tautomerism of N-vinyl-2-pyrrolidinone (NV2P) and N-vinyl-3-pyrrolidinone (NV3P) molecules, and considered the three neutral tautomers for NV2P (1a, 1b, 1c) and the five for NV3P (2a, 2b, 2c, 2d, 2e). Full geometry optimizations were carried out at the HF/6-31G**, HF/6-31+G** and B3LYP/6-31+G** levels in the gas phase and in water. Additionally, single-point MP2/6-31+G** calculations were performed on the HF/6-31+G** optimized geometries. The tautomerization processes in water (=78.54) were studied using the self-consistent reaction field theory. According to our estimations, the tautomerization of NV2P and NV3P molecules are not strongly influenced by the polar solvents, the improvement of the basis set quality and the electron correlation effects in the gas phase and in water. The calculated relative free energies (ΔG) predict that 1a and 2a are the energetically preferred tautomers in the gas phase and in water.     

Tautomeric and conformational equilibrium of 2-nitrosophenol and 9,10-phenanthrenequinonemonooxime: ab initio and NMR study

The tautomeric and conformational equilibrium of 2-nitrosophenol and 9,10-phenanthrenequinonemonooxime was studied by ab initio methods. The geometry optimizations of the structures investigated were done without any geometrical restrictions at HF/6-31G** and MP2/6-31G** levels of theory. The transition structures for tautomeric and rotameric conversions were located. To correct for electron correlation, single-point calculations were carried out up to MP4/6-311G*//MP2/6-31G* level of theory.

Ab initio calculations for 2-nitrosophenol in agreement with the available experimental data define the nitroso form as more stable. It was found that the influence of the correlation energy on the relative stabilities is smaller for the rotamers of the nitroso tautomer but substantially (4–6 kcal/mol) for the oxime forms. It was found that the barrier height of tautomerization reaction is 10.24 kcal/mol.

The structure of the 9,10-phenanthrenequinonemonooxime was studied by solid and liquid state NMR spectroscopy. Ab initio calculations in agreement with our experimental data predict that the compound exists as oxime tautomer and the syn-oxime is most stable. It was found that the solvent influence on the relative stabilities of both isomers: syn- and anti-oxime. While in chloroform solution the syn-oxime is preferred but in DMSO anti-oxime is more stable in energy.

At the MP4/6-311G*//MP2/6-31G**+ZPE level of theory the barrier of tautomerization was predicted to be 10.96 kcal/mol and the rotational barrier around the single C–O bond in the syn-oxime was found to be 7.57 kcal/mol. The rotation is facile and this explains the absence of nitroso tautomers in solution.     

A theoretical study on the origin of pi-facial stereoselectivity in the alkylation of enolates derived from 4-substituted gamma-butyrolactones

The alkylation of 4-methoxymethyl-gamma-butyrolactone enolate with methyl chloride was studied at the B3LYP/6-31+G* level. Conformer search of the free enolate gave 15 unique conformers within 5.39 kcal/mol. The transition structures for both anti- and syn-attacks of methyl chloride on these 15 conformers were located. In all cases, the anti-transition structures are more stable than the corresponding syn-ones. The alkylation of gamma-valerolactone was studied at the MP2, B3LYP, and HF levels of theory with the 6-31+G* basis set in the presence of Li+ and dimethyl ether molecules. Basis set effects were also examined by the comparison of the MP2 results with the 6-31+G*, 6-31+G**, and 6-311+G** basis sets in one case. This study shows that the main source of the anti-selectivity of 4-substituted gamma-butyrolactones is eclipsing strain in the syn-transition structures.     

Conformational analysis and the binding sites of nitrilotriacetamide: A computational study

The conformational analysis of nitrilotriacetamide has been carried out computationally, at both the semi‐empirical AM1 and density functional theory (DFT) (B3LYP/6‐31+G*) levels of theory. The lowest‐energy conformation predicted with the Monte Carlo search method, using the AM1 model, has two amide functionalities aligned on the same side; however, the DFT calculations at B3LYP/6‐31+G* predicted the global minimum with all three acetamide functionalities on the same side in the gas phase. In the aqueous phase, the DFT results predicted the orientations of amides similar to that of the reported crystal structure. The rotation barriers to transition to different low‐energy conformers of nitrilotriacetamide are lower in energy (5.0 kcal/mol) in water. The molecular electrostatic isopotentials (MESP) generated for the selected conformers at DFT level show that the nitrilotriacetamide could interact more effectively with the sodium chloride surface than that of its monomeric unit nitrilomonoacetamide. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Theoretical study of the structure and vibrational spectra of the Ci(2) conformation of 18-crown-6

Optimized geometries of 18-crown-6 (18ce6) were calculated at the HF/6-31G* and B3LYP/6-31G* levels of theory for the D3d, Ci(1) and Ci(2) conformations. At the B3LYP level, the Ci(2) optimized geometry was higher in energy by 23.3 and 18.8 kcal mol(-1) than the Ci(1) and D3d optimized geometries, respectively. Harmonic force field, vibrational frequencies and IR absorption intensities were calculated at the corresponding optimized geometry at the B3LYP level for the Ci(2) conformation. Scaled Ci(2) frequencies were compared with the experimental frequencies of free 18ce6, Ci(1) conformation, and 18ce6-urea complex, Ci(2) conformation. This comparison showed possible misassignments in the fundamental vibrational frequencies of 18ce6.     

Configurational and Conformational Properties of 1,3,7,9-Tetraphospha-Cyclododeca-1,2,7,8-tetraene: An Ab Initio Study and NBO Analysis

An investigation employing the ab initio molecular orbital (MO) and density functional theory (DFT) methods to calculate structural optimization and conformational interconversion pathways for the two diastereoisomeric forms, (±) and meso configurations of 1,3,7,9-tetraphospha-cyclododeca-1,2,7,8-tetraene (1) was undertaken. Two axial symmetrical conformations are found for (±)-1 configuration. (±)-1-TB axial symmetrical form is found to be about 0.35 and 0.99 kcal mol^?1 more stable than (±)-1-Crown axial symmetrical conformation, as calculated by HF/6-31G*//HF/6-31G* and B3LYP/6-31G*//HF/6-31G* levels of theory, respectively. The unsymmetrical meso-1-TBCC form is found to be the most stable geometry, among the various conformations of meso-1 configuration. HF/6-31G*//HF/6-31G* and B3LYP/6-31G*//HF/6-31G* results showed that between the two most stable conformations of (±) and meso configurations, (±)-1-TB is more stable than meso-1-TBCC by about 3.35 and 2.43 kcal mol^?1, respectively. In addition, MP2/6-31G* and B3LYP/6-311+G** results showed that the (±)-1-TB form is about 1.10 and 2.36 kcal mol^?1 more stable than the meso-1-TBCC form. Further, NBO results revealed that in the most stable form of meso configuration (meso-1-TBCC), the sum of the π* allenic antibonding orbital occupancies (Σ π *_occupancy) is greater than dl configuration ((±)-1-TB). Also, NBO results indicated that in the (±)-1-TB conformer, the sum of σ and π allenic moieties bonding orbital deviations (Σ σ _dev+Σ π _dev) from their normal values, is lower than in the meso-1-TBCC form.     

Understanding selenocysteine through conformational analysis,proton affinities,acidities and bond dissociation energies

Density functional methods have been employed to characterize the gas phase conformations of selenocysteine. The 33 stable conformers of selenocysteine have been located on the potential energy surface using density functional B3LYP/6‐31+G* method. The conformers are analyzed in terms of intramolecular hydrogen bonding interactions. The proton affinity, gas phase acidities, and bond dissociation energies have also been evaluated for different reactive sites of selenocysteine for the five lowest energy conformers at B3LYP/6‐311++G*//B3LYP/6‐31+G* level. Evaluation of these intrinsic properties reflects the antioxidant activity of selenium in selenocysteine. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Bowl-shaped hydrocarbons related to C60

Ab initio studies at the HF/6-31G* and B3LYP/6-31G* levels are reported for two bowl-shaped hydrocarbons related to C₆₀: C₃₀H₁₂ and C₃₆H₁₂, of C₃ and C_3v symmetry, respectively. The former has an approximate heat of formation of 211 kcal/mol. Bowl-to-bowl interconversion may occur through a planar (C_3h) form of ca. 64 kcal/mol greater energy having one imaginary vibrational frequency. The larger C₃₆H₁₂ bowl has a calculated ΔH°_f of 265 kcal/mol. Its HF/6-31G*, B3LYP/6-31G*, and MM3 bond lengths are in good agreement with a recent X-ray structure. Chemical shifts for both compounds calculated by the GIAO method are in good agreement with the measured NMR spectra. The observed ¹³C chemical shifts increase with the extent of pyramidalization. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 189–194, 1998