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.     

Theoretical studies on electron delocalisation in selenourea

Ab initio and density functional calculations have been performed on the different possible structures of selenourea(su), urea(u) and thiourea(tu) to understand the extent of delocalisation in selenourea in comparison to urea and thiourea. Selenourea(su-1) withC ₂ symmetry has the minima on the potential energy surface at MP2(fu)/6-31+G* level. The C-N rotational barrier in selenourea is 8.69 kcal/mol, which is 0.29 and 0.11 kcal/mol more than that of urea and thiourea respectively at MP2(fu)/6-31+G* level. N-inversion barrier is 0.55 kcal/mol at MP2(fu)6-31+G* level. NBO analysis has been carried out to understand the nature of different interactions responsible for the electron delocalisation.     

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.     

Ab initio and molecular mechanics calculations on imine derivatives: A study of the rotational barriers and the development of MM2 parameters

The application of molecular mechanics methods for the study of structures has become a standard approach to conformational analysis. The MM2 force field has been extended to include imines. In general, a diverse group of aliphatic and aromatic imine structures can now be treated. The rotational energy profiles and barriers to N-substituted imines about the C_sp2? C_sp2 single bonds adjacent to C?N functional group were calculated via ab initio MO theory. Information obtained from the quantum mechanical calculations at the 3-21G, 6-31G*, and MP4/6-31G* //6-31G* levels was used both to study the phenomena involved and to parameterize MM2. The syn-anti isomerization was also studied, and the mechanistic pathways have been evaluated. In cases where the comparison with experimental data can be made, the agreements are good.     

Theoretical Investigations on Fluorine-Substituted Ethylene Dications C2HnF4-n 2+(n = 0–4)

The optimized geometries and energies of fluorine-substituted ethylene dications C₂H_nF_4-n ²⁺ (n = 0–4) have been investigated by means of ab initio methods. At the MP3/6-31G**//6-31G* + zero-point energy level of theory, the results predict that C₂F₄²⁺ and C₂HF₃²⁺ are planar, while C₂H₄²⁺, C₂H₃F²⁺ and 1,1—C₂H₂F₂²⁺ prefer a perpendicular geometry. For 1,2—C₂H₂F₂²⁺ an energy difference of only 0.3 kcal/mol is found between the (trans) planar and perpendicular structure. The stabilizations attributed to hyperconjugation, fluorine lone-pair donation, and (C? F) double-bond conjugation are discussed. A comparison is made for the C? C and C? F stretching frequencies determined at 6-31G*//6-31G* between the neutral and dicationic species. The theoretically determined ionization energies for the vertical process N⁺ → N²⁺ at the MP3/6-31G*//3-21G level are compared with experimental Q_min values.     

The structure, Raman spectra, and force fields of methylbromosilanes CH3SiClnBr3−n (n = 0, 1, 2)

The laser-excited Raman spectra of liquid CH₃SiCl_nBr_3?n (n = 0, 1, 2) were studied. Quantumchemical calculations of these substances with geometry optimization were performed to determine their harmonic force fields and vibrational frequencies. The calculations were made using the HF/6-31G* and HF/6-311++G** approximations and density functional theory at the B3LYP/6-31G* and B3LYP/6-311++G** levels. An interpretation of the spectra was suggested and the calculated force fields were discussed in comparison with the data on related compounds.     

Predictive Abilities of Scaled Quantum Mechanical Molecular Force Fields: Application to 1,3-Butadiene

The positions of some IR bands of the s-trans-1,3-butadiene-h ₆ and -1,1,2-d ₃ isotopomers in the gas phase have been measured using a Brucker IFS 120 HR spectrometer with a resolution of 2 cm^–1. The structural parameters of the s-trans- and s-gauche-1,3-butadiene conformers were optimized completely at the MP2/6-31G* theoretical level and their MP2/6-31G*//MP2/6-31G* quantum mechanical force fields (QMFFs) were calculated. Using only the experimental vibrational frequencies of s-trans-1,3-butadiene-h ₆ the QMFF of the s-trans conformer was corrected by Pulay's scaling method (eight scale factors were involved). The scaled QMFF was used to calculate the mean vibrational amplitudes and the Coriolis coupling constants of s-trans-1,3-butadiene-h ₆ and the vibrational frequencies of 12 of its deuterated isotopomers. The set of scale factors obtained for correction of the s-trans QMFF was transferred to the QMFF of the s-gauche conformer. Its theoretical vibrational spectrum and those of some deuterated and C¹³ isotopomers were calculated. The ability of this scaling approach (transferring of scale factors) to predict the vibrational frequencies of rotational conformers and their isotopomers, as well as other molecular characteristics, and to permit detection of perturbations of the experimental bands are discussed.     

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     

Molecular orbital studies of polyethylene deformation

Young's modulus E for polyethylene in the chain direction is calculated with molecular orbital theory applied to n-alkanes C₃H₈ through n-C₁₃H₂₈ and analyzed with the cluster-difference method. Semiempirical CNDO, MNDO, and AM1 models and ab initio HF/STO-3G, HF/6-31G, HF/6-31G*, and MP2/6-31G* models are used. Cluster-difference results, when extrapolated to infinite chain length, give E in good agreement with moduli evaluated with molecular cluster or crystal orbital methods, provided minimal basis sets are employed. E decreases from 495 GPa (CNDO) to 336 GPa (MP2/6-31G*) as the level of theory is improved, consistent with established behaviors of the various models. Our calculations do not reproduce earlier molecular cluster or crystal orbital results, which gave E < 330 GPa. The most rigorous MP2/6-31G* model is known to overestimate force constants by ∼ 11%; the scaled modulus E = 299 GPa is in good accord with E = 306 GPa from recent calculations based on experimental vibration frequencies. © 1996 John Wiley & Sons, Inc.     

A comprehensive study of the rotational energy profiles of organic systems by ab initio MO theory,forming a basis for peptide torsional parameters

Ab initio molecular orbital calculations have been carried out on over 50 model organic molecules and ions to provide the data necessary in the determination of torsional parameters for a force field involving polypeptides. The rotational energy profiles were obtained at the HF/6-31G*//HF/6-31G* level. The results were supported, in many cases, by full geometry optimizations and with consideration of correlation corrections at the MP2 level. With the exception of the dihedral angle being studied, all of the molecules were fully optimized with C₁ symmetry. © 1995 by John Wiley & Sons, Inc.     

含苯并五元环方酸衍生物的二阶非线性光学性质的理论研究

利用从头算HF/6-31G*方法对用苯并五元环取代的方酸衍生物体系SQ1～SQ15进行几何构型优化和电子结构计算. 以优化后的构型为基础, 应用ZINDO方法计算电子光谱. 同时应用从头算CPHF/6-31G*方法和半经验FF/AM1, FF/PM3, FF/MNDO等有限场方法计算了分子的二阶非线性光学系数β_μ. 研究几何结构、电子结构和前线分子轨道能与β_μ之间的关系, 为设计性能优良的有机非线性光学材料提供理论指导, 并对这四种计算方法的结果进行了比较.     

苯乙烯基吡啶类化合物几何结构与光谱的理论研究

采用量子化学密度泛函理论(DFT)方法分别在B3LYP/6-31G*, 6-31G**, 6-31＋G*水平上对苯乙烯基吡啶类化合物进行计算研究. 通过在相同水平下的振动频率分析发现苯乙烯基吡啶类化合物具有C₁对称性, 酯基的碳氧原子与苯环形成不同的离域大π键, 空间位阻和共轭效应使得两苯环处于两个不同平面, 二面角在60°与62°之间. 使用含时密度泛函理论(TD-DFT)方法计算第一激发态的电子垂直跃迁能, 得到最大吸收波长λ_max. 计算结果表明末端烷基链的长度对该类化合物的几何结构与振动光谱、电子光谱无影响.     

Proton affinities of molecules containing nitrogen and oxygen: comparing ab initio and semiempirical results to experiments

We report a comparison of theoretical and experimental proton affinities at nitrogen and oxygen sites within a series of small molecules. The calculated proton affinities are determined using the semiempirical methods AM 1, MNDO , and PM 3; the ab initio Hartree–Fock method at the following basis levels: 3-21G //3-21G , 3-21+G //3-21G , 6-31G *//6-31G *, and 6-31+G (d, p)//6-31G *; and Møller–Plesset perturbation calculations: MP 2/6-31G *//6-31G *, MP 3/6-31G *//6-31G *, MP 2/6-31G +(d, p)//6-31G *, MP 3/6-31G +(d, p)//6-31G *, and MP 4(SDTQ )/6-31G +G (d, p)//6-31G *. The semiempirical methods have more nonsystematic scatter from the experimental values, compared to even the minimal 3-21G level ab initio calculations. The thermodynamically corrected 6-31G *//6-31G * proton affinities provide acceptable results compared to experiment, and we see no significant improvement over 6-31G *//6-31G * in the proton affinities with any of the higher-level calculations. © 1992 John Wiley & Sons, Inc.     

(HBNH)n团簇的结构和稳定性

Structures and thermodynamic properties of the imidoboranes （HBNH）n （n=1-16） have been investigated theoretically at the B3LYP/6-31G^＊ level of theory. Needle-shaped oligomers that violate the isolated square rule were found to be more stable than cage isomers. The needle-shaped oligomer with n=16 was predicted to be exceptionally stable at low temperature, hexamer and octamer clusters dominated the gas phase at higher temperature. The highest oligomerization degree of the spontaneous cluster fomation has been estimated. It was concluded that generation of the gas phase （HBNH）n clusters with oligomerization degree n ≥24 was viable, making these species possible intermediates involved in the gas phase generation of BN nanoparticles.     

Dependence of collision‐induced dissociation energy on molecular degrees of freedom as a means to assess relative binding affinity in multivalent complexes

In collision‐induced dissociation mass spectrometry experiments, the collision energy required for dissociation linearly depends on the degrees of freedom in the precursor ion. The magnitude of the slope of this relationship previously has been shown to qualitatively correlate to the relative binding strength of a noncovalently bound, monovalent complex. The goal of the work presented here is to determine if a similar methodology can be applied for assessing relative binding strengths in multivalent species. We have tested the method on complexes formed from 18‐crown‐6 and a variety of protonated, primary alkylamines, [C_nH_2n+1NH₃]⁺ (n = 9, 12, 14, 16 and 18) and alkyldiamines, [H₃NC_nH_2nNH₃]²⁺ (n = 3, 5, 6, 9 and 12), and compared our results with dissociation energies calculated using density functional theory at the B3LYP/6‐31G* level. We found that the method correctly assessed the stronger crown ether/headgroup interaction in the two divalent species (1:1 and 2:1 complexes formed from the diaminoalkanes) compared with the weaker interaction in the monovalent species (1:1 complexes formed from mono‐aminoalkanes). However, the experimental method could not distinguish between the binding strengths of the two divalent complexes, perhaps because their calculated dissociation energies were quite similar. Our preliminary results suggest that this method could potentially be used for a quick and simple analysis of binding strengths in multivalent species if the binding strengths of the species are significantly different from one another. Copyright © 2011 John Wiley & Sons, Ltd.     

The photoelectron spectrum and conformation of phenylphosphine and phenylarsine

He(I) and He(II) photoelectron spectra of phenylphosphine and phenylarsine have been investigated and assigned. The rotational barrier of the phosphino group has been investigated at the MP2/6-31G(d,p)//MP2/6-31G(d,p) and HF/6-31G(d,p)//HF/6-31G(d,p) levels of theory, and that of the arsino group at the HF/6-31G(d,p)//6-31G(d,p) levels of theory. The rotational barrier of the two molecules is nearly the same. The energy difference between the two possible conformers of the molecules is low (1.5 kJ/mol at the MP2/6-31G(d,p) level of theory), allowing nearly free rotation about the P-C bond. The photoelectron spectrum cannot be interpreted by considering the most stable rotamer, but all possible conformers should be taken into account. The present interpretation is consistent with the smalln _p - interaction concluded from other investigations. The rotational barrier ofo-phosphinophenol is significantly larger than for phenylphosphine, and the photoelectron spectrum of this compound can be interpreted by considering a single conformer, and no appreciable interaction between the -system of the ring and the phosphorus lone pair.     

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     

Ab initio study of reactions of hydroxyl radicals with chloro‐ and fluoro‐substituted methanes

Ab initio calculations at the unrestricted Hartree–Fock (UHF) level have been performed to investigate the hydrogen abstraction reactions of ^? OH radicals with methane and nine halogen‐substituted methanes (F, Cl). Geometry optimization and vibrational frequency calculations have been performed on all reactants, adducts, products, and transition states at the UHF/6‐31G* level. Single‐point energy calculations at the MP2/6‐31++G* level using the UHF/6‐31G* optimized geometries have also been carried out on all species. Pre‐ and postreaction adducts have been detected on the UHF/6‐31G* potential energy surfaces of the studied reactions. Energy barriers, ΔE^?, reaction energies, ΔE_r, reaction enthalpies, ΔH_r, and activation energies, E_a, have been determined for all reactions and corrected for zero‐point energy effects. Both E_a and ΔH_r come into reasonable agreement with the experiment when correlation energy is taken into account and when more polarized and diffuse basis sets are used. The E_a values, estimated at the PMP2/6‐31++G* level, are found to be in good agreement with the experimental ones and correctly reproduce the experimentally observed trends in fluorine and chlorine substitution effects. A linear correlation between E_a and ΔH_r is obtained, suggesting the presence of an Evans–Polanyi type of relationship. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 426–440, 2001     

Theoretical study of [Li(H2O)n]+ and [K(H2O)n]+ (n = 1−4) complexes

The geometries, successive binding energies, vibrational frequencies, and infrared intensities are calculated for the [Li(H₂O)_n]⁺ and [K(H₂O)_n]⁺ (n = 1?4) complexes. The basis sets used are 6-31G* and LANL 1DZ (Los Alamos ECP +DZ ) at the SCF and MP 2 levels. There is an agreement for calculated structures and frequencies between the MP 2/6-31G* and MP 2/LANL 1DZ basis sets, which indicates that the latter can be used for calculations of water complexes with heavier ions. Our results are in a reasonable agreement with available experimental data and facilitate experimental study of these complexes. © 1995 John Wiley & Sons, Inc.