Vibrational assignment and structure of trifluorobenzoylacetone. A density functional theoretical study

Molecular structure and vibrational frequencies of 4,4,4-trifluoro-1-phenyl-1,3-butanedione, known as trifluorobenzoylacetone (TFBA), have been investigated by means of density functional theory (DFT) calculations. The results were compared with those of benzoylacetone (BA), acetylacetone (AA), and trifluoroacetylacetone (TFAA). Comparing the calculated and experimental band frequencies and intensities suggests coexisting of both stable cis-enol conformers in comparable proportions in the sample. The energy difference between the two stable chelated enol forms is negligible, 0.96 kcal/mol, calculated at B3LYP/6-311++G** level of theory. The molecular stability and the hydrogen bond strength were investigated by applying the natural bond orbital (NBO) theory and geometry calculations. The theoretical calculations and spectroscopic results indicate that the hydrogen bond strength of TFBA is between those of TFAA and AA, considerably weaker than that of BA.     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

去氢抗坏血酸分子振动光谱的理论研究

采用RHF, MP2, DFT(B3LYP)方法, 以6-311++G**为基组研究了去氢抗坏血酸分子(DHA)的平衡几何构型和振动光谱. 计算结果表明, 采用RHF, B3LYP以及MP2 方法优化得到的几何结构以及频率值是一致的. 采用B3LYP/6-311++G**计算了DHA分子平衡构型下的谐振动力场﹑振动频率和振动强度. 使用Wilson的GF矩阵方法对DHA分子进行了简正坐标分析, 依据所得的势能分布对DHA分子的振动基频进行了合理的理论归属.     

Theoretical investigation of the structure and vibrational spectra of carbamoyl azide

Molecular structure and vibrational frequencies of carbamoyl azide NH2CO-NNN have been investigated with ab initio and density functional theory (DFT) methods. The molecular geometries for all the possible conformers of the molecule were optimized using DFT-B3LYP, DFT-BLYP and MP2 applying the standard 6-311++G** basis set. From the calculations, the molecule was predicted to exist predominantly in cis conformation with the cis-trans rotational barrier of about 7.91-9.10 kcal/mol depending on the level of theory applied. The vibrational frequencies and the corresponding vibrational assignments of carbamoyl azide in Cs symmetry were examined theoretically and the calculated Infrared and Raman spectra of the molecule in the cis conformation were plotted. Observed frequencies for normal modes were compare with those calculated from normal mode coordinate analysis carried out on the basis of ab initio and DFT force fields using the standard 6-311++G** basis set of the theoretical optimized geometry. Theoretical IR intensities and Raman activities are reported.     

Density functional theory study on the structure and vibrational frequencies of glycylglycine

Eleven possible conformers of glycylglycine have been studied by using the BLYP, B3LYP methods of density functional theory and the HF method at the basis set of 6-311++G**. BLYP (using Becke's and Lee-Yang-Parr's correlation functionals), ab initio Hartree-Fock (HF) and hybrid DFT/HF B3LYP calculations have been carried out to study the structure and vibrational spectra of glycylglycine. Glycylglycine crystal structure has been determined by X-ray diffraction analysis. The title compound has been crystallizes in the orthorhombic space group C1, with Z=4. And the unit cell parameters are: a=8.1184(12)A, b=9.5542(14)A, c=7.8192(11)A and V=577.95(15)A(3). Molecular conformation calculations have got 11 possible conformers. In these possible conformers, the most stable one has been selected. The BLYP/6-311++G** and scaled HF/6-311++G** frequencies correspond well with available experimental assignments of the normal vibrational modes. Comparison of the observed fundamental vibrational frequencies of glycylglycine and calculated results by density functional B3LYP and Hartree-Fock (HF) methods indicates that B3LYP is superior to the scaled Hartree-Fock (HF) for molecular vibrational issues.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

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.     

Vibrational assignment and structure of dibenzoylmethane. A density functional theoretical study

Molecular structure and vibrational frequencies of 1,3-diphenyl-1,3-propanedione, known as dibenzoylmethane (DBM), have been investigated by means of density functional theory (DFT) calculations. The results were compared with those of benzoylacetone (BA) and acetylacetone (AA), the parent molecule. IR and Raman spectra of DBM and its deuterated analogue were clearly assigned.The calculated hydrogen bond energy of DBM is 16.15 kcal/mol, calculated at B3LYP/6-311++G** level of theory, which is 0.28 kcal/mol more than that of AA. This result is in agreement with the vibrational and NMR spectroscopy results. The molecular stability and the hydrogen bond strength were investigated by applying the Natural Bond Orbital analysis (NBO) and geometry calculations. The theoretical calculations indicate that the hydrogen bond in DBM is relatively stronger than that in BA and AA.     

Theoretical determinations of ionization potentials of dopamine

Adiabatic and vertical ionization potentials (IPs) of nine conformers of dopamine in the gas phase are determined using density functional theory (DFT) B3LYP, B3P86, B3PW91 methods and high level ab initio HF method with 6-311++G** basis set, respectively. And the nine stable cationic states have been found in the ionization process of dopamine. Vertical ionization potentials of nine conformers of dopamine are calculated using the older outer-valence Green’s function (OVGF) calculations at 6-311++G** basis set. Vibrational frequencies and infrared spectrum intensities of G1b and G1b⁺ at B3LYP/6-311++G** level are discussed.     

Effect of methyl groups substitution on the strength of intramolecular hydrogen bonding of naphthazarin: DFT and NBO studies

A density functional theory (DFT) study-based method B3LYP/6-311++G** was carried out to investigate the methyl groups substitution effect on the structure and the strength of intramolecular hydrogen bonding in naphthazarin (NZ) (5,8-dihydroxy-1,4-naphthoquinone). The full geometry optimization of molecular structures, the difference between the energies of hydrogen-bonded and non-hydrogen-bonded rotamers, and the proton chemical shift of the hydroxyl groups in NZ and its methyl substituents obtained at the B3LYP/6-311++G** level. The vibrational frequencies of all samples and their deuterated analogues were calculated at the same theoretical level. The ¹H chemical shifts for NZ and its methyl substituents were computed at the B3LYP/6-311++G** level using the gauge-including atomic orbital method. Furthermore, in order to investigate the changes in bond order, electron density, electron delocalization, and steric effects caused by methyl substituents, natural bond orbital analysis were carried out at the B3LYP/6-311++G** level. After comparing these effective parameters in methyl substituents with those of their parent, NZ, we concluded that, in general, intramolecular hydrogen bonding strength increases by substituting methyl groups in the different positions of NZ.     

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.     

Investigation of conformational stability and vibrational spectra of halomethylsulfonyl isocyanates

The conformational behavior and structural stability of chloro- and fluoromethylsulfonyl isocyanates were investigated by quantum mechanical DFT and ab initio MP2 calculations. The 6-311++G^** basis set was employed to include polarization and diffuse functions in the calculations. The molecules were found to exist in a mixture of two stable gauche conformations. The potential scans were calculated from which the rotational barriers could be estimated. The vibrational frequencies and spectra were computed at B3LYP/6-311++G^** level. The potential energy distributions were then calculated to provide tentative vibrational assignment for the normal modes of the stable conformers of both molecules.     

FT-IR and FT-Raman spectroscopic investigation, computed vibrational frequency analysis and IR intensity and Raman activity peak resemblance analysis on 2-nitroanisole using HF and DFT (B3LYP and B3PW91) calculations

Fourier-transform Raman and infrared spectra of 2-nitroanisole are recorded (4000-100 cm(-1)) and interpreted by comparison with respective theoretical spectra calculated using HF and DFT method. The geometrical parameters with C(S) symmetry, harmonic vibrational frequencies, infrared and Raman scattering intensities are determined using HF/6-311++G (d, p), B3LYP/6-311+G (d, p), B3LYP/6-311++G (d, p) and B3PW91/6-311++G (d, p) level of theories. A detailed vibrational spectral analysis has been carried out and assignments of the observed fundamental bands have been proposed on the basis of peak positions and relative intensities. The results of the calculations have been used to simulate IR and Raman spectra for the molecule that showed good agreement with the observed spectra. The SQM method, which implies multiple scaling of the DFT force fields has been shown superior to the uniform scaling approach. The vibrational frequencies and the infrared intensities of the C-H modes involved in back-donation and conjugation are also investigated.     

Si~2Cl~6分子的内旋转能垒和振动基频的理论研究

用从头算方法HF/6-31G^*^*和密度函方法B3LYP/6-31G^*^*,对Si~2Cl~6分子的平衡几何构型进行优化,优化的结果与实验结果吻合得较好.并用上述两种不同的方法计算Si~2Cl~6分子的内旋转能垒,结果分别为8.786和6.694kJ/mol,其中DFT方法的计算结果与实验结果4.18kJ/mol吻合得较好.对Si~2Cl~6分子的振动基频进行计算.用HF/6-31G^*^*SQM力场所计算的频率理论值与实验值的平均误差为7.3cm^-^1,用B3LYP/6-31G^*^*未标度的力场所计算的频率理论值与实验值的平均误差为6.0cm^-^1.该密度泛函方法(B3LYP/~6-31G^*^*)的理论计算值比用HF/6-31G^*^*标度后的SQM力场计算的频率与实验值(除Si--Si键扭转振动基频之外的11条振动基频)吻合得更好.并给出了Si--Si键扭转振动基频的预测值。     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

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.     

Potential energy scans and vibrational assignments of cyclopropanecarboxylic acid and cyclopropanecarboxamide

The structural stability and internal rotations in cyclopropanecarboxylic acid and cyclopropanecarboxamide were investigated by the DFT-B3LYP and the ab initio MP2 calculations using 6-311G** and 6-311+G** basis sets. The computations were extended to the MP4//MP2/6-311G** and CCSD(T)//MP2/6-311G** single-point calculations. From the calculations the molecules were predicted to exist predominantly in the cis (C=O group eclipses the cyclopropane ring) with a cis-trans barrier of about 4-6kcal/mol. The OCOH torsional barrier in the acid was estimated to be about 12-13kcal/mol while the corresponding OCNH torsional barrier in the amide was calculated to be about 20kcal/mol. The equilibrium constant k for the cis<-->trans interconversion in cyclopropanecarboxylic acid was calculated to be 0.1729 at 298.15K that corresponds to an equilibrium mixture of about 85% cis and 15% trans. The vibrational frequencies were computed at the DFT-B3LYP level. Normal coordinate calculations were carried out and potential energy distributions were calculated for the low energy cis conformer of the molecules. Complete vibrational assignments were made on the basis of normal coordinate calculations and comparison with experimental data of the molecules.     

A study of internal rotations and vibrational spectra of oxiranemethanol (glycidol)

The conformational stability and the C–O and O–H internal rotations in oxiranemethanol were investigated at the DFT-B3LYP/6-311G**, MP2/6-311G** and MP4(SDQ)/6-311G** levels of theory. Three minima were predicted in the CCOH potential energy scans of the molecule to have relative energies of about 2 kcal/mol or less and all were calculated to have real frequencies upon full optimization of structural parameters at the DFT and the MP2 levels of calculations. The Cg1 (H bond inner) conformation was predicted to be the lowest energy conformation for oxiranemethanol in excellent agreement with an earlier microwave study. The equilibrium mixture was calculated from Gibb's free-energy changes to be about 79% Cg1, 17% G1g and 3% Gg1 at the B3LYP/6-311G** level and about 87% Cg1, 11% G1g and 2% Gg1 at the MP2/6-311G** level for oxiranemethanol at 298.15 K. No conclusive evidence was obtained for the presence of high-energy form in the liquid phase of oxiranemethanol. The vibrational frequencies of oxiranemethanol in its three stable forms were computed at the B3LYP level and complete vibrational assignments were made for the lowest energy Cg1 form on basis of calculated and experimental data of the molecule.     

儿茶素结构和振动光谱的密度泛函理论研究

本文利用密度泛函理论方法B3LYP在6-311 G(d,p)水平上研究了儿茶素及其异构体表儿茶素分子的几何构型,计算结果与实验所得结构参数一致.计算了儿茶素分子平衡构型下的力常数,使用Wilson的GF矩阵方法计算了振动基频以及相应的势能分布,据此结合理论计算的光谱强度,对儿茶素分子的振动基频进行了完善合理的理论归属.