Conformational stability of peroxynitric acid molecule and the barriers to internal rotation about the O−O and N−O bonds

The structure of conformers and potential curves of the internal rotation (PCR) about the O?O and N?O bonds in peroxynitric acid (PNA) were calculated by the unrestricted Hartree-Fock-Roothaan method. The standard valence-split 6–31G and 6–31G^* basis sets were used. The presence of two maxima on every curve has been shown. To refine the values of barriers to the internal rotation in the regions of minima and maxima of PCR, calculations taking into account the electron correlation energy have been carried out at the second- and fourth-order Møller-Plesset level of perturbation theory (MP2 and MP4, respectively). At the MP4/6-61G^* level of approximation, the barriers to the rotation about the O?O bond are equal to 8.6 kJ mol^?1 and 14.7 kJ mol^?1, and both barriers to the rotation about the N?O bond are equal to 33.5 kJ mol^?1. The results are compared with those published for PCR in hydrogen peroxide and peroxynitric acid.     

Effects of different GIAO and CSGT models and basis sets on 2-aryl-1,3,4-oxadiazole derivatives

The direct molecular structure implementations of the gage-including atomic orbital (GIAO), individual gages for atoms in molecules (IGAIM) and continuous set of gage transformations (CSGT) methods for calculating nuclear magnetic shielding tensors at both the Hartree-Fock (HF) and density functional (B3LYP) levels of theory with 6-31G(d), 6-311G(d), 6-31++G(d,p), 6-311++G(d,p), and 6-311++G(df,pd) basis sets are presented. Dependence on the ¹H and ¹³C NMR chemical shifts on the choice of method and basis set have been investigated. Also, these chemical shifts of 2-aryl-1,3,4-oxadiazoles 5a–g have been performed related to dihedral angles (C4–C3–C2–O) of two conformers. The optimized molecular geometries and ¹H and ¹³C chemical shift values of 2-aryl-1,3,4-oxadiazoles 5a–g in the ground state have been obtained. The linear correlation coefficients of ¹³C NMR chemical shifts for these molecules were given. The new nuclear magnetic shielding tensors of tetramethylsilane (TMS) were calculated. The data of 2-aryl-1,3,4-oxadiazole derivatives display significant molecular structure and NMR analysis. Also, these provide the basis for future design of efficient materials having the 1,3,4-oxadiazole core.     

Ab initio and infrared spectral study of 4′-substituted phenylthiolbenzoates

The geometry of the 4′-cyano-(4′-CNPTB) and the 4′-methoxy-(4′-MePTB) phenylthiolbenzoates were obtained by ab initio calculations employing 6–31G basis set at Hartee-Fock level of theory. The results predict an extended form of the molecules and torsional angle between the phenyl rings at 90.85(6)⁰ and 90.87(3)⁰, respectively. On the basis of vibrational analysis the frequency assignment was carried out. The calculated frequencies were compared with the experimental IR spectral data in carbon tetrachloride, carbon disulfide solutions and in solid state.     

Calculation of the hydrogen binding energies of complexes between formamide and adenine-thymine pair

The hydrogen bonding of complexes formed between formamide and adenine-thymine base pair has been completely investigated in the present study. In order to gain deeper insight into structure, charge distribution, and energies of complexes, we have investigated them using density functional theory at 6–311++G(d, p), 6–31G, 6–31+G(d), and 6–31++G(d, p) level. Seven stable cyclic structures (ATF1–ATF7) being involved in the interaction has been found on the potential energy surface. In addition, for further correction of interaction energies between adenine—thymine and formamide, the basis set superposition error associated with the hydrogen bond energy has been computed via the counterpoise method using the individual bases as fragments. The infrared spectrum frequencies, IR intensities and the vibrational frequency shifts are reported.     

Experimental and computational studies on (<Emphasis Type="Italic">Z</Emphasis>)-1-((4-phenylamino)phenylamino)-methylene)naphthalen-2(1<Emphasis Type="Italic">H</Emphasis>)-one

The Schiff base compound (Z)-1-((4-phenylamino)phenylamino)methylene)naphthalen-2(1H)-one has been synthesized and characterized by IR, UV–Vis, and X-ray single-crystal determination. Molecular geometry from X-ray experiment of the title compound in the ground state have been compared using the Hartree–Fock (HF) and density functional method (B3LYP) with 6−31G(d,p) basis set. Calculated results show that density functional theory DFT and HF can well reproduce the structure of the title compound. Using the time-dependent density functional theory (TD-DFT) and Hartree–Fock (TD-HF) methods, electronic absorption spectra of the title compound have been predicted and a good agreement with the TD-DFT method and experimental ones is determined. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6−31G(d,p) basis set by applying the polarizable continuum model (PCM). The total energy of the title compound decreases with increasing polarity of the solvent. In addition, DFT calculations of the title compound, molecular electrostatic potential (MEP), natural bond orbital analysis (NBO), and non-linear optical (NLO) properties were performed at B3LYP/6−31G(d,p) level of theory.     

A theoretical study of the interactions between N,N-dimethylformamide and xylenes

The ab initio and density functional (DFT) methods were performed on binary systems of N,N-dimethylformamide (DMF) with xylenes (o-, or m-, or p-xylene), and seven stable configurations were obtained with no imaginary frequencies. To obtain the interaction energies of these complexes, single-point energy calculations with basis set superposition error (BSSE) correction were carried out at B3LYP/6-31G* and MP2/6-31G* levels. The structures, Chelpg (charges from electrostatic potentials using a grid-based method) charge distribution and bond characteristics of the mentioned complexes were calculated. The results indicated the presence of double C–H···O hydrogen bonds between DMF and xylenes in these complexes and the interaction energies of hydrogen bonding between DMF and xylene systems decreased in the following sequence: DMF–o-xylene: a1 > DMF–m-xylene: b1 > DMF–p-xylene: c1.     

Structure of N-chloromethyllactams and features of the interaction of atoms in them as a result of ab initio calculations

Quantum-chemical calculations have been carried out by the RHF/6-31G(d) and MP2/6-31+G(d) methods of molecules of N-chloromethylpyrrolidone, N-chloromethylcaprolactam, N-chloromethyl-succinimide, and N-chloromethylphthalimide with full optimization of their geometry, and also N-chloromethylpyrrolidone molecule by the RHF/6-31G(d) method at various angles of rotation of the CH₂Cl group around the C―N bond. It was shown that the lower frequencies of the ³⁵Cl NQR of the first two molecules in comparison with the later are mainly determined by the high populations of the p _σ -orbitals of their Cl atoms. The population of the orbitals of the unshared electron pair of the N atom is practically unchanged on rotating the CH₂Cl group, but the N atom polarizes the C―Cl bond in the indicated molecule. This does not confirm the supposed p,σ*-conjugation in the Cl―C―N grouping. Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 10, pp. 1537–1544, October, 2008.     

Ab initio study of molecular structure and internal rotation in methyldicyanophosphine and methyldiisocyanophosphine. The Raman spectrum of methyldicyanophosphine and its interpretation with the use of scaling ofab initio force fields

The equilibrium geometric parameters and structures of transition states of internal rotation for the molecules of methyldicyanophospine MeP(CN)₂ and its isocyano analog MeP(NC)₂ were calculated by the RHF and MP2 methods with the 6–31G^* and 6–31G^** basis sets. At the MP2 level, the total energy of cyanide is −35 kcal mol⁻¹ lower than that of isocyanide and the barriers to internal rotation of methyl group for MeP(CN)₂ and MeP(NC)₂ are 2.2 and 2.7 kcal mol⁻¹, respectively. For both molecules, the one-dimensionalab initio potential functions of internal rotation approximated by a truncated Fourier series were used to determine the frequencies of torsional transitions by solving direct vibrational problems for a non-rigid model. The Raman spectrum of crystalline MeP(CN)₂ was recorded in the range 3500–50 cm⁻¹. The vibrational spectra of this compound were interpreted by scalingab initio force fields calculated by the RHF and MP2 method. The vibrational spectrum of methyldiisocyanophosphine was predicted with the use of the obtained scale factors. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1703–1714, September, 1998.     

Mechanisms of inversion of bond configuration at the tetrahedral boron atom in five-membered chelate cycles

Mechanisms of inversion of the bond configuration at the tetrahedral boron center in five-membered chelate cycles of the 1,3,2-oxazaborolidine and 1,3,2-oxazaborolidene molecules were studied by theab initio MP2(full)/6-31G^** method. It was shown that enantiotopomerization occurs by a dissociative mechanism with the cleavage of the B←N bond and the formation of acyclic intermediates with tricoordinate planar boron atom. The calculated energy barriers to inversion of tetrahedral bond configurations at boron centers in the two chelate complexes are equal to 13.1 and 15.4 kcal mol⁻¹, respectively. In contrast to 1,3,2-oxazaborolidine, internal rotation about the B−O bond in its unsaturated analog makes an appreciable contribution to the reaction coordinate. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 2, pp. 250–255 February, 1999.     

Theoretical study of the structural character of weakly bonding silicon carbonyl complexes

The structures, properties and the bonding character for sub-carbonyl Si, SiCO and Si(CO)₂, in singlet and triplet states have been investigated using complete-active-space self-consistent field (CASSCF), density functional theory and second-order M?ller–Plesset methods with a 6-311+G* basis set. The results indicate that the SiCO species possesses a ³∑⁻ ground state, and the singlet ¹Δ excited state is higher in energy than the ³∑⁻ state by 17.3 kcalmol⁻¹ at the CASSCF–MP2/6-311+G* level and by 16.4 kcalmol⁻¹ at the CCSD(T)/6-311+G* level. The SiCO ground state may be classified as silene (carbonylsilene), and its CO^δ− moiety possesses CO⁻ property. The formation of SiCO causes the weakening of CO bonds. The Si–C bond consists of a weak σ bond and two weak π bonds. Although the Si–C bond length is similar to that of typical Si–C bonds, the bond strength is weaker than the Si–C bonds in Si-containing alkanes; the calculated dissociation energy is 26.2 kcalmol⁻¹ at the CCSD(T)/6-311+G* level. The corresponding bending potential-energy surface is flat; therefore, the SiCO molecule is facile. For the bicarbonyl Si systems, Si(CO)₂, there exist two V-type structures for both states. The stablest state is the singlet state (¹A₁), and may be referred to the ground state. The triplet state (³B₁) is energetically higher in energy than the ¹A₁ state by about 40 kcalmol⁻¹ at the CCSD(T)/6-311 + G* level. The bond lengths in the ¹A₁ state are very close to those of the SiCO species, but the SiCO moieties are bent by about 10°, and the CSiC angles are only about 78°. The corresponding ³B₁ state has a CSiC angle of about 54° and a SiCO angle of about 165°, but its Si–C and C–O bonds are longer than those in the ¹A₁ state by about 0.07 and 0.03 ?, respectively. This Si(CO)₂ (¹A₁) has essentially silene character and should be referred to as a bicarbonyl silene. Comparison of the CO dissociation energies of SiCO and Si(CO)₂ in their ground states indicates that the first CO dissociation energy of Si(CO)₂ is smaller by about 7 kcalmol⁻¹ than that of SiCO; the average one over both CO groups is also smaller than that of SiCO. A detailed bonding analysis shows that the possibility is small for the existence of polycarbonyl Si with more than three CO. This prediction may also be true for similar carbonyl complexes containing other nonmetal and non-transition-metal atoms or clusters. Received: 17 April 2002 / Accepted: 11 August 2002 / Published online: 4 November 2002 Acknowledgements. This work was supported by the National Natural Science Foundation of China (29973022) and the Foundation for Key Teachers in University of the State Ministry of Education of China. Correspondence to: Y. Bu e-mail: byx@sdu.edu.ch     

AB Initio calculations and vibrational spectra of CH3SC≡CH and CH3SC≡CSCH3

Methylthio- (MTE) and bis-methylthioethyne (BMTE) molecules are calculated by the SCF MO method (geometry optimization, basis set 6–31G*/MP2). The calculated internal rotation barriers of methyl groups are 7.12 kJ/lmole for MTE and 12.86 kJ/mole for BMTE (both groups are simultaneously rotated). The s-gosh-orientation of the thiomethyl fragments corresponds to a stable conformation of BMTE. The estimated values of the s-cis- and s-trans-barriers of mutual rotation of SCH₃ groups about the axis of the C≡C bond are 13.61 and 12.54 kJ/mole, respectively. Conformationally sensitive MOs and vibration frequencies are established. An analysis of the experimental IR absorption and Raman spectra and the calculated vibrational spectrum makes it possible to conclude that in the liquid phase the BMTE molecules also have an s-gosh-conformation. Translated fromZhumal Strukturnoi Khimii, Vol. 39, No. 4, pp. 602–609, July–August, 1998.     

Ab initio studies of tetrafluoroethylene and tetratrifluoromethylethylene molecules in the ground and “twisted” states

Ab initio calculations with full optimization of geometry have been carried out with the 6–31 G^* basis set on tetrafluoroethylene (with the unrestricted Hartree-Fock method—UHF and the second-order Moller-Plesset perturbation theory—MP2) and tetratrifluoromethylethylene (with UHF) molecules in the singlet ground and triplet biradical states. The symmetry of the tetrafluoroethylene molecule in the triplet biradical state was demonstrated to differ from that of ethylene (D _2d ) due to the deviation of fluorine atoms from CCFF plane. The MP2 optimized geometries of ethylene and tetrafluoroethylene were used for higher level calculations (MP3, MP4, CCSD). The energy of the ground state singlet-biradical triplet splitting decreases in the series: ethylene>tetrafluoroethylene> tetratrifluoromethylethylene. These data on energy splitting explain the increase in reactivity toward the [2+2]-cycloaddition on going from ethylene to tetrafluoroethylene. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 605–607, April, 1998.     

Linear and non-linear optical properties of some donor–acceptor oxadiazoles by ab initio Hartree-Fock calculations

The molecular hyperpolarizability of some donor–acceptor oxadiazoles was investigated using ab initio methods. Ab initio optimizations were performed at the Hartree-Fock level using different basis sets, starting with the minimal basis set, and then split valence sets. The first hyperpolarizabilities were calculated at the Hartree-Fock level employing the corresponding basis sets using Gaussian 98W. In general, the first hyperpolarizability is dependent on the choice of method and basis set. In order to understand this phenomenon in the context of molecular orbital picture, we examined the molecular HOMOs and molecular LUMOs generated via HF/6-31G level. It has also been calculated the polarizability, anisotropy of polarizability and ground state dipole moment of all the molecules. Several of these oxadiazoles display significant second-order molecular nonlinearity, β(8.57–195.05 × 10⁻³⁰ esu) and provide the basis for future design of efficient nonlinear optical materials having the oxadiazole core.     

Study of <Emphasis Type="Italic">cis</Emphasis>–<Emphasis Type="Italic">trans</Emphasis> isomerization mechanism of 3,3′-azobenzene disulphonate in the lowest singlet and triplet electronic states by density functional theory

Abstract  

The cis–trans isomerization pathways of 3,3′-azobenzene disulphonate in the S₀ and T₁ states are studied by DFT method at the B3LYP/6-31G(d,p) level. In the S₀ state, the cis–trans isomerization concerns the complex pathway that is characterized by the inversion of one NNC angle combined with rotation around the NC bond, and the three sequential transition states are also found on the potential energy profile. Therefore, the cis–trans isomerization of 3,3′-azobenzene disulphonate can be understood in terms of a pathway involving successive rotation, inversion, and rotation processes. The energy barrier of the S₀ state is 22.79 kcal mol⁻¹. In the T₁ state, the isomerization mainly concerns the rotational pathway around the NN double bond, and the two isomers are connected through only one transition state. The isomerization of the T₁ state is related to a lower energy barrier, 5.02 kcal mol⁻¹, but requires a change in spin-multiplicity.     

Molecular structures of acetylene derivatives of tin

The geometry and force fields of the bis(trimethylstannyl)acetylene molecule (a conformer withD _3d symmetry corresponding to a minimum of the total energy of the molecule) were calculated by the RHF and MP2(fc) methods. The effective core potential in SBK form with the optimized 31G^* valence basis set was employed in the case of Sn atoms. The 6–31G^** and 6–311G^** basis sets were used for carbon and hydrogen atoms. Vibrational spectra of the light and perdeuterated isotopomers of bis(trimethylstannyl)acetylene were interpreted using the procedure of scaling the quantum-chemical force fields. For Part 5, see Ref. 1. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 616–626, April, 2000.     

Structure and DFT calculations of 2-{[3-methyl-3-phenyl-cyclobutyl)-thiazol-2-yl]-hydrazonomethyl}-phenol

The title compound 2-{[3-Methyl-3-phenyl-cyclobutyl)-thiazol-2-yl]-hydrazonomethyl}-phenol (C₂₁H₂₁N₃S₁O₁) crystallizes in the P-1 triclinic space group with a = 5.8880(4) ?, b = 9.5618(5) ?, c = 17.0484(10) ?, α = 80.214(5)°, β = 80.532(5)°, γ = 80.116(5)°. In addition to molecular geometry and packing from X-ray experiment, we have also calculated the molecular geometry and vibrational frequencies of the title compound in the ground state using density functional theory DFT (B3LYP) with the 6–31G(d,p) basis set. Apart from this, the structure of the title compound is characterized by ¹H NMR, ¹³C NMR, IR and UV-vis. Spectra, and the experimental emission energies are compared with the HOMO-LUMO energy gaps calculated by the DFT method.     

Ab initio calculations of atomic interaction in the ClCH2OCH3 molecule

Ab initio calculations of both ClCH₂OCH₃ and ClCH₂CH₃ molecules and various ClCH₂OCH₃ structures with fixed angles of rotation of the methoxy group about the C−O bond were performed by the restricted Hartree-Fock method in the valence-split 6–31 G^* basis set with full optimization of the geometry. The populations of the valent p-orbitals of the chlorine atoms in these molecules have been analyzed. The³⁵Cl NQR frequencies and the asymmetry parameters of the electric field gradient (EFG) at the³⁵Cl nuclei have been calculated. Good agreement with experimental NQR frequencies was obtained for the calculations where only the populations of the less diffuse 3p-components of these orbitals were used. The³⁵Cl NQR frequency in ClCH₂OCH₃ is lower than that in ClCH₂CH₃ due to the higher population of the less diffuse component of the p_σ-orbital of the Cl atom in the former molecule. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 431–434, March, 1997.     

Force field parameters for sulfates and sulfamates based on ab initio calculations: Extensions of AMBER and CHARMm fields

Ab initio self-consistent field (SCF) Hartree-Fock calculations of sulfates R? O? SO₃(?1) (R = Me, Et, i-Pr) and sulfamates R? NHSO₃(?1) (R = H, Me, Et, i-Pr) were performed at the 4-31G(*S*N) //3-21G(*S*N) basis set levels, where asterisks indicate d functions on sulfur and nitrogen atoms. These standard levels were determined by comparing calculation results with several basis sets up to MP2/6-31G*//6-31G*. Several conformations per compound were studied to obtain molecular geometries, rotational barriers, and potential derived point charges. In methyl sulfate, the rotational barrier around the C? O bond is 1.6 kcal/mol at the MP2 level and 1.4 kcal/mol at the standard level. Its ground state has one of three HCOS torsion angles trans and one of three COSO torsion angles trans. Rotation over 60° around the single O? S bond in the sulfate group costs 2.5 kcal/mol at the MP2 and 2.1 kcal/mol at the standard level. For ethyl sulfate, the calculated rotational barrier in going from the ground state, which has its CCOS torsion angle trans, to the syn-periplanar conformation (CCOS torsion angle cis) is 4.8 kcal/mol. However, a much lower barrier of 0.7 kcal/mol leads to a secondary gauchelike conformation about 0.4 kcal/mol above the ground state, with the CCOS torsion angle at 87.6°. Again, one of the COSO torsion angles is trans in the ground state, and the rotational barrier for a 60° rotation of the sulfate group amounts to 1.8 kcal/mol. For methyl sulfamate, the rotational barriers are 2.5 kcal/mol around the C? N bond and 3.3 kcal/mol around the N? S bond. This is noteworthy because sulfamate itself has a calculated rotational barrier around the N? S bond of only 1.7 kcal/mol. These and other data were used to parameterize the well-known empirical force fields AMBER and CHARMm. When the new fields were tested by means of vibrational frequency calculations at the 6-31G*//6-31G* level for methyl sulfate, sulfamate, and methyl sulfamate ground states, the frequencies compared favorably with the AMBER and CHARMm calculated frequencies. The transferability of the force parameters to β-D -glucose-6-sulfate and isopropyl sulfate appears to be better than to isopropyl sulfamate. © 1995 by John Wiley & Sons, Inc.     

Quantum chemical and experimental studies on the structure and vibrational spectra of substituted 2-pyranones

A systematic study on the structural characteristics of the 2-pyranone ring containing molecules with bromine, nitrile, and amide substituents at the C-3 position in the ring is conducted in the electronic ground (S ₀) state by DFT calculations using the B3LYP/6-311++G** method. The geometrical structure of the bromine substituted compound, which shows potent hepatoprotective activity, is studied both in the ground (S ₀) and first excited singlet (S ₁) states using RHF/6-311++G** and CIS/6-311++G** methods respectively. The molecules are found to exist in two isomeric forms gauche and trans that have the enthalpy difference of less than 3.32 kcal/mol; the latter is the preferred orientation in the gaseous phase. The S ₁ state is a ¹(π,π*) state that arises π-electron transfer from the region of a double bond in the pyranone ring to the region of the internuclear bond connecting the 2-pyranone and benzene rings. A complete vibrational analysis is conducted for the 3-bromo-6-(4-chlorophenyl)-4-thiomethyl-2H-pyran-2-one molecule based on the experimental infrared spectra in the 50–4000 cm⁻¹ region and DFT/6- 311++G** computations of vibrational frequencies for the gauche and trans isomeric forms. Spectral assignments based on the potential energy distribution along the internal coordinates confirm the nonplanar structure of the molecule.