乙烯醇锂的从头算研究

用限制的HF/3-21G和HF/6-31G^*优化乙烯醇锂的几种可能构型,比较了它们的稳定性.用限制的HF/3-1G,从乙醛开始,探讨了气相反应生成乙烯醇锂的机理,并在MP2水平上用6-31G^*基组计算了反应热.     

DFT Calculation of Glycine with Methanols Clusters

Density functional theory （DFT） calculations are reported for the structures of neutral and zwitterionic glycine-（CHaOH）n where n=1-6. Initial geometries of the clusters of neutral and zwitterionic glycine with 1-6 methanol molecules are fully optimized at B3LYP/6-31＋G^＊ level of theory. The lowest energy configurations are located and their hydrogen bond structures are analyzed. Theoretical prediction reveals that the methanols prefer to locate near the carboxylic acid group for the small clusters （n_〈3） with the neutral form while the configurations with the methanols bridging the acid and the amino group are favorite in the zwitterionic form clusters. When the number of the methanol molecules in the clusters reaches five and six, the two forms tend to be isoenergetic.     

Investigation of the use of B3LYP zero-point energies and geometries in the calculation of enthalpies of formation

The use of B3LYP/6–31G^* zero-point energies and geometries in the calculation of enthalpies of formation has been investigated for the enlarged G2 test set of 148 molecules [J. Chem. Phys. 106 (1997) 1063]. A scale factor of 0.96 for the B3LYP zero-point energies gives an average absolute deviation nearly the same as scaled HF/6–31G^* zero-point energies for G2, G2(MP2), and B3LYP/6–311 + G(3df,2p) enthalpies. A scale factor of 0.98, which has been recommended in some studies, increases the average absolute deviation by about 0.2 kcal/mol. Geometries from B3LYP/6–31G^* are found to do as well as MP2/6–31G^* geometries in the calculation of the enthalpies of formation.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

A semi-empirical and ab initio combined approach for the full conformational searches of gaseous lysine and lysine–H2O complex

A full structural search of the canonical, zwitterionic, protonated and deprotonated lysine conformers in gas phase is presented. A total of 17,496 canonical, 972 zwitterionic, 11,664 protonated and 1458 trial deprotonated structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were initially optimized at the AM1 level, and the resulting structures were determined at the B3LYP/6-311G^* level. A total of 927 canonical, 730 protonated and 193 deprotonated conformers were found, but there were no stable zwitterionic structures in the gas phase. The most stable conformers of the canonical, protonated and deprotonated lysine were further optimized at the B3LYP/6-311++G^** level. The energies of the most stable structures were determined at the MP2/6-311G(2df,p) level and the vibrational frequencies were calculated at the B3LYP/6-311++G^** level. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, vertical ionization energies, enthalpies, Gibbs free energies and conformational distributions of gaseous lysine were presented. Numerous new structures are found and the lowest-energy lysine conformer is more stable than the existing one by 1.1 kcal/mol. Hydrogen bonds are classified and may cause significant red-shifts to the associated vibrational frequencies. The calculated proton affinity/dissociation energy and gas-phase basicity/acidity are in good agreement with the experiments. Calculations are also presented for the canonical lysine–H₂O and zwitterionic lysine–H₂O clusters. Interaction between lysine and H₂O significantly affects the relative conformational stabilities. Only one water molecule is sufficient to produce the stable zwitterionic structures in gas phase. The lowest-energy structure is found to be zwitterions when applying the conductor-like polarized continuum solvent model (CPCM) to the lysine–H₂O complexes.     

Molecular structure and conformational stability of allylisocyanate—post-Hartree–Fock and density functional theory studies

The molecular structure and conformational stability of allylisocyanate (CH₂CHCH₂NCO) molecule was studied using the ab initio and DFT methods. The geometries of possible conformers, C-gauche (δ=120°, θ=0°) (δ=C=C–C–N and θ=C–C–N=C) and C-cis N-trans (δ=0° and θ=180°) were optimized employing HF/6-31G^*, MP2/6-31G^* levels of theory of ab initio and BLYP, B3LYP, BPW91 and B3PW91 methods of DFT implementing the atomic basis set 6-311+G(d,p). The structural and physical parameters of the above conformers were discussed with the experimental and theoretical values of the related molecules, methylisocyanate and 3-fluoropropene. It has been found that the N=C=O bond angle is not linear as the experimental result for both the conformers and the theoretical bond angle is 173°. The rotational potential energy surfaces have been performed at the HF/6-31G^*, and MP2/6-31G^* levels of theory. The Fourier decomposition potentials were analysed at the HF/6-31G^*, and MP2/6-31G^* levels of theory. The HF/6-31G^* level of theory predicted that the C-gauche conformer is more stable than the C-cis N-trans conformer by 0.41 kJ/mol, but the MP2 and DFT methods predicted the C-cis N-trans conformer is found to be more stable than the C-gauche conformer. The calculated chemical hardness value at the HF/6-31G^* level of theory predicted the C-cis N-trans form is more stable than C-gauche form, whereas the chemical hardness value at the MP2/6-31G^* level of theory favours the slight preference towards the C-gauge conformer.     

An ab initio theoretical study of electronic structure and properties of 2'-deoxyguanosine in gas phase and aqueous media

Molecular geometries of two structural forms of 2'-deoxyguanosine (keto-N9R and keto-N7R, R = the sugar moiety) considering both the C2'-endo and C3'-endo conformations of the sugar ring and those of the complexes of these species with two water molecules each were optimized employing the ab initio RHF procedure. A mixed basis set consisting of the 6-311+G* basis set for the nitrogen atom of the amino group and the 4-31G basis set for all the other atoms was used. The RHF calculations were followed by correlation correction of the total energy at the MP2 level. Both the structural forms of 2'-deoxyguanosine were solvated using the polarized continuum model (PCM) of the self-consistent reaction field (SCRF) theory and the corresponding RHF optimized geometries at the RHF and MP2 levels. Geometry optimization was also performed in aqueous media using the Onsager model at the RHF level using the above-mentioned mixed basis set, and subsequently, using the reoptimized geometries, single-point MP2 calculations were performed. It is found that both the keto-N9R and keto-N7R forms of 2'-deoxyguanosine as well as their complexes with two water molecules each would occur, particularly at the water-air interface. Though the normal Watson-Crick-type base pairing would not be possible with the keto-N7R form of 2'-deoxyguanosine(G*), two other (G*-C and G*-T) base pairing schemes may occur with this form of the nucleoside, which may cause mutation. The present calculated geometry of the keto-N9R form of the anti-conformation of 2'-deoxyguanosine including the dihedral angle chi(CN) agree satisfactorily with the available crystallographic results. The present results also agree satisfactorily with those obtained by other authors earlier for the keto-N9R form of 2'-deoxyguanosine using B3LYP and MP2 methods employing the 6-31G* basis set. Using transition state calculations, it is shown that tautomerism of guanine and other similar molecules where the tautomers would coexist would be facilitated by the occurrence of the H(+) and OH(-) fragments of water molecules. Further, this coexistence of the two tautomers appears to make the C8 carbon atom located between the N7 and N9 nitrogen atoms susceptible to attack by the OH(-) group. Thus, an explanation is obtained for the efficient formation of the reaction product 8-hydroxy-2'-deoxyguanosine, which serves as a biomarker for oxidative damage to DNA in biological systems.     

Conformational stability of 1-bromo-2-fluoroethane from temperature dependent FT-IR spectra of rare gas solutions

Variable temperature (−55 to −135°C) studies of the infrared spectra (3500–400 cm⁻¹) of 1-bromo-2-fluoroethane, BrCH₂CH₂F, dissolved in liquid krypton and xenon have been recorded. From these data, the enthalpy difference has been determined to be 108±9 cm⁻¹ (1.296±0.113 kJ/mol) and 112±8 cm⁻¹ (1.346±0.098 kJ/mol) from the krypton and xenon solutions, respectively, with the trans conformer the more stable rotamer. Complete vibrational assignments are presented for both conformers which are consistent with the predicted frequencies obtained from the ab initio MP2/6-31G^* calculations. The optimized geometries, conformational stabilities, harmonic force fields, infrared intensities, Raman activities, and depolarization ratios have been obtained from RHF/6-31G^* and/or MP2/6-31G^* ab initio calculations. These quantities are compared to the corresponding experimental quantities when appropriate. Structural parameters and conformational stability have also been obtained from MP2/6-311+G^** calculations. Combining the ab initio predicted structural parameters with the microwave rotational constants, r_o parameters have been obtained for the gauche conformer.     

异硫氰酸与甲亚胺环加成反应机理的理论研究

对异硫氰酸与甲亚胺形成四元或六元环产物的环加成反应进行了理论研究,结果表明,一分子异硫氰酸与一分子甲亚胺形成四元环(1:1)产物P1的反应(1)为经过一个两性离子中间体的分步反应,其中第二步为速控步骤,其活化热垒为107.86kJ/mol.此外,反应(1)的中间体还可与另一甲亚胺或异硫氰酸分子继续反应形成两个不同的六元环(1:2或2:1)产物P2或P3;这两个反应均为协同反应,其活化势垒分别为15.88和21.82kJ/mol.这些结果与当异硫氰酸酯与亚胺发生环加成反应时,只有类似于P2和P3的两种六元环产物生成的实验事实一致。     

Ab initio studies of aspartic acid conformers in gas phase and in solution

Systematic and extensive conformational searches of aspartic acid in gas phase and in solution have been performed. For the gaseous aspartic acid, a total of 1296 trial canonical structures and 216 trial zwitterionic structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 139 canonical conformers were found, but no stable zwitterionic structure was found. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, and vertical ionization energies of the canonical conformers were determined. Single-point energies were also calculated at the MP2/6-311++G** and CCSD/6-311++G** levels. The equilibrium distributions of the gaseous conformers at various temperatures were calculated. The proton affinity and gas phase basicity were calculated and the results are in excellent agreement with the experiments. The conformations in the solution were studied with different solvation models. The 216 trial zwitterionic structures were first optimized at the B3LYP/6-311G* level using the Onsager self-consistent reaction field model (SCRF) and then optimized at the B3LYP/6-311++G** level using the conductorlike polarized continuum model (CPCM) SCRF theory. A total of 22 zwitterions conformers were found. The gaseous canonical conformers were combined with the CPCM model and optimized at the B3LYP/6-311++G** level. The solvated zwitterionic and canonical structures were further examined by the discrete/SCRF model with one and two water molecules. The incremental solvation of the canonical and zwitterionic structures with up to six water molecules in gas phase was systematically examined. The studies show that combining aspartic acid with at least six water molecules in the gas phase or two water molecules and a SCRF solution model is required to provide qualitatively correct results in the solution.     

A combined effective fragment potential-fragment molecular orbital method. II. Analytic gradient and application to the geometry optimization of solvated tetraglycine and chignolin

The gradient for the fragment molecular orbital (FMO) method interfaced with effective fragment potentials (EFP), denoted by FMO∕EFP, was developed and applied to polypeptides solvated in water. The structures of neutral and zwitterionic tetraglycine immersed in water layers of 2.0, 2.5, 3.0, 3.5, 4.0, and 4.5 A? are investigated by performing FMO∕EFP geometry optimizations at the RHF∕cc-pVDZ level of theory for the solutes. The geometries optimized with FMO-RHF∕EFP are compared to those from the conventional RHF∕EFP and are found to be in very close agreement. Using the optimized geometries, the stability of the hydrated zwitterionic and neutral structures is discussed structurally and in terms of energetics at the second-order M?ller-Plesset theory (MP2)∕cc-pVDZ level. To demonstrate the potential of the method for proteins, the geometry of hydrated chignolin (protein data bank ID: 1UAO) was optimized, and the importance of the inclusion of water was examined by comparing the solvated and gas phase structures of chignolin with the experimental NMR structure.     

Theoretical Studies on theStructures and Reactivity of Stannylenoides(Ⅱ)——Structures and Isomerization of Stannylenoid H2SnNaF

The structures of stannylenoid HzSnNaF have been studied using RHF/3-21G and RHF/3-21G^* gradient method.     

Reaction between ��-phenyl-��-alanine and 4-hydroxy-4-methylpent-2-ynenitrile: a spectroscopic and quantum chemical study

The sequence of isomeric transformations of intermediate adducts formed upon nuceleophilic addition of ??-phenyl-??-alanine to 4-hydroxy-4-methylpent-2-ynenitrile affording a neutral or two alternative zwitterionic forms (with protonated imino or amino groups) of the product was studied quantum chemically by the B3LYP/6-31G(d,p) method and by UV spectroscopy. In the isolated state, the equilibrium is completely shifted toward the neutral form. Taking into account the polyhydrate environment increases the portion of the zwitterionic form bearing a protonated imino group of the dihydrofuran cycle. Protonation of the amino group of amino acid is thermodynamically unfavorable.     

Computational treatment of the microsolvation of neutral and zwitterionic forms of alanine

B3LYP/6-31++G** and MP2/6-31++G**//B3LYP/6-31++G** calculations are reported for the structures of neutral alanine–(H₂O)_n and zwitterionic alanine–(H₂O)_n clusters where n = 2–10. Optimized geometries and energies were obtained. In general, with an increasing number of water molecules, the hydrated zwitterionic form becomes more thermodynamically stable. In the presence of six or seven water molecules, the energetics indicate that the two forms may coexist. Eight water molecules are sufficient to computationally guarantee the reported experimental observation of zwitterionic dominance in solution.     

双环金属铱(Ⅲ)异腈配合物的二阶非线性光学性质

采用密度泛函理论(DFT)方法对双环金属Ir(III)异腈配合物的非线性光学(NLO)性质进行计算研究。用B3PW91(UB3PW91)(金属原子采用LANL2DZ基组,非金属原子采用6-31G*基组)方法对配合物进行几何结构优化。在优化构型基础上,采用B3PW91(UB3PW91)和B3LYP(UB3LYP)方法计算了配合物的第一超极化率(βtot),并用CAM-B3LYP(UCAM-B3LYP)(金属原子采用LANL2DZ基组,非金属原子采用6-31G**基组)方法模拟配合物的吸收光谱。结果表明,主配体的取代基(R1)和副配体的取代基(R2)对第一超极化率值贡献不大。配合物发生氧化还原反应,电荷转移方式增多,电荷转移程度增大,使βtot值显著增加,其中1a+([(C∧N)2Ir(CNR)2]+(R=CH3))发生氧化反应和还原反应的βtot值分别增大了75倍和144倍。因此,这类双环金属铱(III)异腈配合物的氧化还原反应可以有效地调节其二阶NLO性质。     

Ab initio molecular orbital study of dialuminum-ethylene complexes

The geometries of the various isomeric dialuminum-ethylene complexes Al₂C₂H₄ have been optimized using the STO-3G, 3-21G^(*) and 6-31G^(**) basis sets. Vibrational frequency calculations show that only two out of the six structures examined are thermodynamically stable. Unlike the case of Li₂C₂H₄, 1,2-dialuminoethane is most stable and more ethane-like than 1,2-dilithioethane perhaps due to the larger covalent character of the C---Al bond. Besides, 1,2-dialuminoethane has been found, in contrast to ethane and 1,2-dilithioethane, to have no stable gauche conformer.     

锂钠复合超价氧化物的电子结构和稳定性

采用6-31G^*基组,用B3LYP密度泛函方法,对OLi_nNa₂(n=1-4)及OLi₂Na系列分子的所有可能结构进行了几何优化和频率计算,得到各分子的几种能量极小点.除OLi₄Na₂采取D_2d对称性外,该系列其它4种分子都采取C_2v对称作为它们的能量最低点.OLi_nNa₂(n=2,3,4)分子解离出一个Li原子的离解能分别为212.25432,68.49208和245.55896kJ/mol,表明此类分子在热力学上是相当稳定的.     

An ab initio study of water clusters in gas phase and bulk aqueous media: (H2O)n,n=1–12

Geometries of several clusters of water molecules including single minimum energy structures of n‐mers (n=1–5), several hexamers and two structures of each of heptamer to decamer derived from hexamer cage and hexamer prism were optimized. One structural form of each of 11‐mer and 12‐mer were also studied. The geometry optimization calculations were performed at the RHF/6‐311G* level for all the cases and at the MP2/6‐311++G** level for some selected cases. The optimized cluster geometries were used to calculate total energies of the clusters in gas phase employing the B3LYP density functional method and the 6‐311G* basis set. Frequency analysis was carried out in all the cases to ensure that the optimized geometries corresponded to total energy minima. Zero‐point and thermal free energy corrections were applied for comparison of energies of certain hexamers. The optimized cluster geometries were used to solvate the clusters in bulk water using the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory, the 6‐311G* basis set, and the B3LYP density functional method. For the cases for which MP2/6‐311++G** geometry optimization was performed, solvation calculations in water were also carried out using the B3LYP density functional method, the 6‐311++G** basis set, and the PCM model of SCRF theory, besides the corresponding gas‐phase calculations. It is found that the cage form of water hexamer cluster is most stable in gas phase among the different hexamers, which is in agreement with the earlier theoretical and experimental results. Further, use of a newly defined relative population index (RPI) in terms of successive total energy differences per water molecule for different cluster sizes suggests that stabilities of trimers, hexamers, and nonamers in gas phase and those of hexamers and nonamers in bulk water would be favored while those of pentamer and decamer in both the phases would be relatively disfavored. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 90–104, 2001     

An ab initio study on the geometries and vibrational frequencies of hydrazoic acid and methyl azide

The geometric parameters for hydrazoic acid and methyl azide were optimized at the HF/6-31G^** and MP2/6-31G^** levels and the vibrational frequencies of the compounds were calculated by use of these optimized geometries. The experimental frequencies are assigned on the basis of the calculated results. The effects of deutero-substitution and substitution of hydrogen in HN₃ by a methyl group are also discussed.     

A theoretical investigation on the potential energy surface of CN2H rotation of the encapsulated 1-bicyclo[2.2.1]heptyldiazirine

The potential energy surface (PES) of CN₂H rotation of the encapsulated 1-bicyclo[2.2.1]heptyldiazirine (BHD) inside a molecular container: Cram’s hemicarcerand (CH) was explored using two different DFT involved ONIOM methods: B3LYP/6-31G^**//ONIOM(B3LYP/6-31G^*: AM1) and B971/6-31G^**//ONIOM(B971/6-31G^*: AM1). The free-state PES of CN₂H rotation was also calculated, respectively by B3LYP/6-31G^**//B3LYP/6-31G^* and B971/6-31G^**//B971/6-31G^* methods for comparison. The findings in this study have shown that the PES profiles differ from each other notably in the two states. In the encapsulated state the rotation barrier corresponding to the free-state conversion with the largest rotation barrier increases by about 2 kcal/mol, which has exceeded the largest rotation barrier in the free-state. The conformational preference behavior towards certain BHD isomers, which might be in better conformational compatibility with the container, has been demonstrated.