Steric and Electronic Structure of Selenoanisole: A Quantum-Chemical Study

The potential functions of internal rotation around the C_sp2-Se bond in selenoanisole were ob- tained by quantum-chemical calculations in the approximations HF/3-21G(d), HF/6-31G(d), MP2(f)/6-31G(d), and B3LYP/6-31G(d). The calculations were performed in the range of variation of the torsion angle (between the planes of the benzene ring and C_sp2-Se-C_sp3 bonds) from 0° to 90° with 15° step. The energy minimum is in the region of the orthogonal conformation ( 90°), and the energy maximum, in the region of the planar form ( 0°). The rotation barriers (kJ mol^{- 1}) are as follows: HF/3-21G(d), 9.20; HF/6-31G(d), 13.13; MP2(f)/6-31G(d), 10.25; and B3LYP/6-31G(d), 6.41. The geometric parameters, Koopmans ionization potentials, and dipole moments are given. The energies, degrees of hybridization, populations of the lone electron pairs of Se, energies of their interaction with the antibonding * orbitals of the benzene ring, and electron density distributions were determined in terms of the natural bond orbital approach.     

Quantum-Chemical Study of Anisole Molecule

The potential functions of internal rotation around the C -O bond in the C₆H₅OCH₃ molecule were obtained by HF/6-31G(d), MP2(f)/6-31G(d), and B3LYP/6-31(d) calculations. Hartree-Fock calculations reveal a fourfold barrier to internal rotation around the C -O bond. The MP2 and B3LYP calculations reveal a twofold barrier with a height of 7.78 and 10.70 kJ mol^{- 1}, respectively (corrected for the zero vibration energy). The molecular geometries, first Koopmans ionization potentials, and dipole moments are reported. Calculations for liquid anisole in the self-consistent reactive field (SCRF) continual model give the results that only slightly differ from the results obtained for the isolated molecule in a vacuum. Within the framework of the Natural Bond Orbitals formalism, the following parameters were determined: energy, degree of hybridization, and population of oxygen lone electron pairs and energy of their interaction with antibonding ^* orbitals of the aromatic ring.     

An Ab Initio Molecular Orbital Theory and Density Functional Theory (DFT) Study of Conformers and Rotamers of 4-Substituted (Methyl, Hydroxymethyl, Methanoyl, Ethanoyl, Cyano, Fluoro, Chloro, Bromo, Acetoxy) Tetrahydro-2H-thiopyran-1,1-dioxides

The structures and energies of axial and equatorial conformers and rotamers of 4-substituted tetrahydro-2H-thiopyran-1,1-dioxides (tetrahydrothiopyran-1,1-dioxides, thiacyclohexane-1,1-dioxides, thiane-1,1-dioxides, and 1,1-dioxothianes; CH₃, CH₂OH, CHO, COCH₃, CN, F, Cl, Br, and OCOCH₃) were calculated using the hybrid density functionals B3LYP, B3P86, and B3PW91, as well as MP2 and the 6-31G(d), 6-31G(2d), 6-31G(3d), 6-31G(d,p), and 6-31+G(d) basis sets. MP2/6-31+G(d)/ /HF/6-31+G(d) [–G° = 1.73 kcal/mol], B3P86/6-31G(d) [–G° = 1.75 kcal/mol], and B3PW91/6-31G(d) [–G° = 1.85 kcal/mol] gave conformational free energy (G°) values at 180 K for 4-methyltetrahydro-2H-thiopyran-1,1-dioxide which were similar to the reported experimental values for methylcyclohexane (–G° = 1.80 kcal/mol), 4-methyltetrahydro-2H-thiopyran (–G° = 1.80 kcal/mol), and other 4-methyl-substituted heterocycles. All levels of theory showed that the conformational preferences of the 4-methanoyl (4-formyl), 4-ethanoyl (4-acetyl), and 4-cyano substituents were small. The HF calculations gave conformational free energy (G°) values for 4-chlorotetrahydro-2H-thiopyran-1,1dioxide which were closer to the experimental value than the MP2 and density functional methods. The best agreement with available experimental data for 4-bromotetrahydro-2H-thiopyran-1,1-dioxide was obtained from the HF/6-31G(2d), HF/6-31G(3d), and B3LYP/6-31G(2d) calculations, and, for 4-acetoxytetrahydro-2H-thiopyran-1,1-dioxide, from the HF/6–31G(3d) calculations. The conformational free energies (G°) and relative energies (E) of the conformers and rotamers have been compared with the correspondingly substituted cyclohexanes and tetrahydro-2H-thiopyrans and are discussed in terms of dipole–dipole (electrostatic) interactions and repulsive nonbonded interactions (steric) in the most stable axial and equatorial conformers. The axial S=O bond lengths are shorter than the equatorial S=O bond lengths and the C2–C3 bond lengths in the substituents with carbon-bonded to the ring are shorter than the C3–C4 and C4–C-5 bond lengths. In contrast, the C2–C3 bond lengths in the 4-halogen and 4-acetoxy substituents are longer than the C3–C4 and C4–C-5 bond lengths.     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Potential Functions of Internal Rotation around the Csp2-X Bonds and Intermolecular Interactions in the Compounds CH2 = CHXCH3 (X = O,S, Se)

The potential functions of internal rotation around the Csp ²-X bonds in the CH₂ = CHXCH₃ molecules (X = O, S, Se) were determined by MP2/6-31G* and MP2/6-31G** calculations. The stationary points were identified by solution of vibrational problems. The rotation barriers were evaluated taking into account the zero-point vibration energy. The intramolecular interactions were considered in terms of the method of natural bond orbitals. The degrees of hybridization, energies, and populations of the orbitals of the lone electron pairs of the O, S, and Se atoms, the energies of their donor-acceptor interaction with the antibonding * and * orbitals of the double bond, and the natural atomic charges in various conformations were determined.     

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     

Why betaine crystallizes in high local Cs symmetry. An ab initio MO and DFT study of anhydrous betaine and betaine monohydrate

A theoretical study of the structure, charge distribution, rotational barrier and fundamental vibrations of anhydrous betaine (CH₃)₃NCH₂COO (trimethylglycine) was carried out and compared with available experimental data. Calculations were carried out at HF, MP2 and B3LYP levels using a 6-31+G(d,p) basis set. The calculated rotational barrier of the betaine carboxylic group is 40.5 kJ/mol at the MP4(SDQ)/6-311G(d,p)//HF/6-31+G(d,p) level of theory. The rotation of the carboxylic group changes the molecule from a highly symmetric (C _s ) conformation into a twisted conformation resulting in shortening of the molecule by about 50 pm. Natural population analysis (NPA) indicates intramolecular interaction between the carboxylic oxygen and the nearest methyl hydrogens resulting in internal hydrogen bonding. MP4(SDQ)/6-311G(d,p) single-point NPA calculations on a betaine monohydrate model taken from the X-ray geometry show an expected weakening in the internal hydrogen bond. Calculations explain why betaine preferentially crystallizes in high local C _s symmetry. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 7 December 1998     

Torsion potentials and electronic structure of trifluoromethoxy- and trifluoromethylthiobenzene: an ab initio study

Potential functions of internal rotation about the C_sp²X bonds in molecules C₆H₅XCF₃ (X=O, S) were calculated at the second-order Møller-Plesset perturbation level of theory with 6-31G(d) basis set. The profile of the potential function and the rotation barrier (ΔE^#=3.0 kJ/mol) found for C₆H₅OCF₃ suggest that, depending on experimental conditions, there can be either free rotation about the C_sp²O bond or the conformational equilibrium is shifted to the side of the orthogonal form. The rotational barrier for C₆H₅SCF₃ is 14.7 kJ/mol and the molecule exists in the stable orthogonal conformation. The nature of hybridization, energy and population of lone electron pairs (LPs) on the oxygen and sulfur atoms were considered by using the Natural Bond Orbital (NBO) method. The energy of interactions of the LPs with antibonding π^∗-orbitals of the aromatic moiety were estimated for different conformations. The distribution of electron density in the molecules was discussed. The results were compared with analogous calculations on the molecules C₆H₅XCH₃.     

An ab initio study of the structure,dissociation energy,and heat of formation of Na2S

The equilibrium geometries and fundamental frequencies of Na₂S are calculated at HF, MP2(FC, FU), and MP3 with the 6–31G(d) basis set and at HF and MP2(FC, FU) with the 6–31G(d) basis set, respectively. The total energy at MP2(FU)/6–31G(d)-optimized geometry is computed at MP4 with 6–311G(d, p), 6–311 + G(d, p), and 6–311G(2df, p), at QCISD(T)/6–311G(d, p), and at MP2/6–311G(3df, 2p) levels, respectively. The dissociation energy, the atomization energy, and the heat of formation for Na₂S are evaluated using the G1 and G2 models. The calculated results indicated that Na₂S in its ground state was a bent structure (C_2v). Electron correlation corrections on the bending angle are very significant. The equilibrium geometrical parameters are R_e(Na-S) = 2.45 Å and ∠Na-S-Na = 111.13° at the MP2(FU)/6–31G(d) level. The theoretically estimated dissociation energy, total atomization energy, and heat of formation are 67.07, 117.55, and 0.35 kcal mol⁻¹, respectively, at 298.15 K. © 1997 John Wiley & Sons, Inc.     

Quantum-Chemical Study of Electronic and Steric Structure of Vinyltetrazoles: I. 2-Substituted 5-Vinyltetrazoles

The total Mulliken charges on the carbon atoms of the vinyl group, populations of S-trans-(N¹)conformers, and internal rotation energies were calculated ab initio (HF/6-31G**, MP2/6-31G**, and MP2/6-31G**//AM1) for a series of 2R-5-vinyltetrazoles (R = CH₃, C₂H₅, i-C₃H₇, t-C₄H₉, C₆H₅). The calculation results were compared to the available experimental data.     

On the additivity of basis set effects in some simple fluorine containing systems

The reactions F + H₂ → HF + H, HF → H + F, F → F⁺ + e^? and F + e^? → F^? were used as simple test cases to assess the additivity of basis set effects on reaction energetics computed at the MP4 level. The 6-31G and 6-311G basis sets were augmented with 1, 2, and 3 sets of polarization functions, higher angular momentum polarization functions, and diffuse functions (27 basis sets from 6-31Gd, p) to 6-31 ++ G(3df, 3pd) and likewise for the 6-311G series). For both series substantial nonadditivity was found between diffuse functions on the heavy atom and multiple polarization functions (e.g., 6-31 + G(3d, 3p) vs. 6-31 + G(d, p) and 6-31G(3d, 3p)). For the 6-311G series there is an extra nonadditivity between d functions on hydrogen and multiple polarization functions. Provided that these interactions are taken into account, the remaining basis set effects are additive to within ±0.5 kcal/mol for the reactions considered. Large basis set MP4 calculations can also be estimated to within ±0.5 kcal/mol using MP2 calculations, est. E_MP4(6-31 ++ G(3df, 3pd)) ≈ E_MP4(6-31G(d, p)) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31G(d, p)) or E_MP4(6-31 + G(d, p) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31 + G(d, p)) and likewise for the 6-311G series.     

An ab initio study of the structures and enthalpies of the hydrogen-bonded complexes of the acids H2O,H2S,HCN, and HCl with the anions OH−, SH−, CN−, and Cl−

Ab initio calculations at second-order Møller-Plesset perturbation theory with the 6-31 + G(d,p) basis set have been performed to determine the equilibrium structures and energies of a series of negative-ion hydrogen-bonded complexes with H₂O, H₂S, HCN, and HCl as proton donors and OH^–, SH^–, CN^–, and Cl^– as proton acceptors. The computed stabilization enthalpies of these complexes are in agreement to within the experimental error of 1 kcal mol^–1 with the gas-phase hydrogen bond enthalpies, except for HOHOH^–, in which case the difference is 1.8 kcal mol^–1. The structures of these complexes exhibit linear hydrogen bonds and directed lone pairs of electrons except for complexes with H₂O as the proton donor, in which cases the hydrogen bonds deviate slightly from linearity. All of the complexes have equilibrium structures in which the hydrogen-bonded proton is nonsymmetrically bound, although the symmetric structures of HOHOH^– and ClHCl^– are only slightly less bound than the equilibrium structures. MP2/6-31 + G(d,p) hydrogen bond energies calculated at optimized MP2/B-31 + G(d,p) and at optimized HF/6-31G(d) geometries are similar. Using HF/6-31G(d) frequencies to evaluate zero-point and thermal vibrational energies does not introduce significant error into the computed hydrogen bond enthalpies of these complexes provided that the hydrogen-bonded proton is definitely nonsymmetrically bound at both Hartree-Fock and MP2.     

Intramolecular Hypervalent X←X Interaction in Mono- and Bicyclic Heteroelement Compounds Containing Se and Te Atoms

Chalcogen-containing heterapentalene and quasimonocyclic compounds having SeÄSeÄSe and TeÄTeÄTe triads or SeÄSe and TeÄTe diads were studied by the ab initio [MP2(full)/6-31G**, MP2(fc)/6-31+G**, and MP2(fc)/LANL2DZ] and DFT methods (B3LYP/6-31G^**, B3LYP/6-31+G^**, and B3LYP/LANL2DZ). Heterapentalene compounds were found to be stable as planar bicyclic structures having a C _2v symmetry. The stability of quasimonocycic -chalcogenovinyl aldehydes increases with increase in the electron-acceptor power of the substituent at the X atom.     

Chemical Reactivity of Oxo‐ and Aza‐γ‐lactam Rings

Different mechanisms for the alkaline hydrolysis of oxo and aza‐γ‐lactam rings have been studied by ab initio calculations at the MP2/6‐31+G*//MP2/6‐31+G* and B3LYP/6‐31+G*//B3LYP/6‐31+G* levels. The tetrahedral intermediate can undergo two different reactions, the cleavage of the C₂−N₂ bond (the classical mechanism) and the cleavage of the C₂−X₆ bond (X=O, N). Both compounds present similar energy barriers for the classical fragmentation, and show considerably lower barriers for the alternative mechanism. Because of this reactivity, the compounds studied are expected to be β‐lactamase inhibitors.     

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.     

Structures and stability of N9, N9 − and N9 + clusters

 Ab initio molecular orbital calculations for N₉, N⁻ ₉ and N⁺ ₉ isomers were carried out at the HF/ 6-31G*, B3PW91/6-31G*, B3LYP/6-31G* and MP2/ 6-31G* levels of theory. Stable equilibrium geometric structures were determined by harmonic vibrational frequency analyses at the HF/6-31G*, B3PW91/6-31G* and B3LYP/6-31G* levels of theory. The most stable free-radical N₉ cluster is structure 1 with C _{2 v} symmetry and that of anion N⁻ ₉ is structure 3 with C _s symmetry. Only one stable structure of the N⁺ ₉ cation with C _{2 v} symmetry was predicted. Their potential application as high-energy-density materials has been examined. Received: 15 June 1999 / Accepted: 11 October 1999 / Published online: 14 March 2000     

Rotational barriers of disilane,hexafluorodisilane, and hexamethyldisilane: Ab initio,density functional,and molecular mechanics (MM3) studies

We investigated structures, vibrational frequencies, and rotational barriers of disilane (Si₂H₆), hexafluorodisilane (Si₂F₆), and hexamethyldisilane (Si₂Me₆) by using ab initio molecular orbital and density functional theories. We employed four different levels of theories (i.e., HF/6–31G*, MP2/6–31G*, BLYP/6–31G*, and B3LYP/6–31G*) to optimize the structures and to calculate the vibrational frequencies (except for Si₂Me₆ at MP2/6–31G*). MP2/6–31G* calculations reproduce experimental bond lengths well, while BLYP/6–31G* calculations largely overestimate some bond lengths. Vibrational frequencies from density functional theories (BLYP/6–31G* and B3LYP/6–31G*) were in reasonably good agreement with experimental values without employing additional correction factors. We calculated the ΔG^‡(298 K) values of the internal rotation by correcting zero-point vibration energies, thermal vibration energies, and entropies. We performed CISD/6–31G*//MP2/6–31G* calculations and found the ΔG^‡(298 K) values for the internal rotation of Si₂H₆, Si₂F₆, and Si₂Me₆ to be 1.36, 2.06, and 2.69 kcal/mol, respectively. The performance of this level was verified by using G2 and G2(MP2) methods in Si₂H₆. According to our theoretical results, the ΔG^‡(298 K) values were marginally greater than the ΔE^‡(0 K) values in Si₂F₆ and Si₂Me₆ due to the contribution of the entropy. In Si₂H₆ the ΔE^‡(0 K) and ΔG^‡(298 K) values were coincidently similar due to a cancellation of two opposing contributions between zero-point and thermal vibrational energies, and entropies. Our calculated ΔG^‡(298 K) values were in good agreement with experimental values published recently. In addition, we also performed MM3 calculations on Si₂H₆ and Si₂Me₆. MM3 calculated rotational barriers and thermodynamic properties were compared with high level ab initio results. Based on this comparison, MM3 calculations reproduced high level ab initio results in rotational barriers and thermodynamic properties of Si₂H₆ derivatives including vibrational energies and entropies, although large errors exist in some vibrational frequencies. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1523–1533, 1997     

Is tetrahedral H42+ a minimum? Anomalous behavior of popular basis sets with the standard p exponents on hydrogen

The nature of the tetrahedral H₄²⁺ stationary point (minimum or triply degenerate saddle) depends remarkably upon the theoretical level employed. Harmonic vibrational analyses with, e.g., the 6-31G** (and 6-31 + +G**) and Dunning's [4s2p1d;2s1p] [D95(d,p)] basis sets using the standard p exponent suggest (erroneously) that the T_d geometry is a minimum at both the HF and MP2 levels. This is not the case at definitive higher levels. The C₃H₄²⁺ structure with an apical H is another example of the failure of the calculations with the 6-31G**, 6-311G**, and D95(d,p) basis sets. Even at MP2/6-31G** and MP2/ cc-pVDZ levels, the C_3v structure has no negative eigenvalues of the Hessian. Actually, this form is a second-order saddle point as shown by the MP2/6-31G** calculation with the optimized exponent. The D_4h methane dication structure is also an example of the misleading performance of the 6-31G** basis set. In all these cases, energy-optimized hydrogen p exponents give the correct results, i.e., those found with more extended treatments. Optimized values of the hydrogen polarization function exponents eliminate these defects in 6-31G** calculations. Species with higher coordinate hydrogens may also be calculated reliably by using more than one set of p functions on hydrogen [e.g., the 6-31G(d,2p) basis set]. Not all cases are critical. A survey of examples, also including some boron compounds, provides calibration. © 1993 John Wiley & Sons, Inc.     

Ab initio study of substituent effect on the addition of hydrogen fluoride to fluoroethylenes

An ab initio 3-21G study of the direct addition of HF to C₂H_nF_(4–n), with n = 0 to 4, has been performed to investigate the effect of the substituent on the reaction. Geometry optimization of all charge-transfer complexes and transition states has been done. Standard analysis of activation energies of addition reactions, vibrational and thermodynamical analysis, as well as Morokuma energy decomposition, BSSE correction, PMO analysis, and Pauling bond orders were used to explain the results. A subset of the reactions, including that of C₂H₄ as reference one and the two most favorable cases, was also studied at the MP2/6–31G(d,p)//HF/6–31G(d,p) level. The barriers so obtained are in agreement with the indirectly found from experimental data. It was found that the effect of the substituent is not monotonic for the additions. Decomposition of the interaction energy is shown to be adequate to explain this nonmonotonic behavior. The implications for laser chemistry of the addition of hydrogen halides to fluorosubstituted olefins is briefly discussed.     

Molecular structure and vibrational assignment of melaminium phthalate by density functional theory (DFT) and ab initio Hartree-Fock (HF) calculations

The molecular geometry, the normal mode frequencies and corresponding vibrational assignment of melaminium phthalate (C₃H₇N₆⁺·C₈H₅O₄⁻) in the ground state were performed by HF and B3LYP levels of theory using the 6-31G(d) basis set. The optimized bond length numbers with bond angles are in good agreement with the X-ray data. The vibrational spectra of melaminium phthalate which is calculated by HF and B3LYP methods, reproduces vibrational wave numbers with an accuracy which allows reliable vibrational assignments. The title compound has been studied in the 4000–100 cm⁻¹ region where the theoretical evaluation and assignment of all observed bands were made.