Ab initio vibrational state calculations with a quartic force field: applications to H2CO, C2H4, CH3OH, CH3CCH, and C6H6

For polyatomic molecules, n-mode coupling representations of the quartic force field (nMR-QFF) are presented, which include terms up to n normal coordinate couplings in a fourth-order polynomial potential energy function. The computational scheme to evaluate third-and fourth-order derivatives by finite differentiations of the energy is fully described. The code to generate the nMR-QFF has been implemented into GAMESS program package and interfaced with the vibrational self-consistent field (VSCF) and correlation corrected VSCF (cc-VSCF) methods. As a demonstration, fundamental frequencies have been calculated by the cc-VSCF method based on 2MR-QFF for formaldehyde, ethylene, methanol, propyne, and benzene. The applications show that 2MR-QFF is a highly accurate potential energy function, with errors of 1.0-1.9% relative to the experimental value in fundamental frequencies. This approach will help quantitative evaluations of vibrational energies of a general molecule with a reasonable computational cost.     

Vibrational spectral analysis and first hyperpolarizability studies of 1-bromonaphthalene based on ab initio and DFT methods

In this work, the experimental and theoretical vibrational spectra of 1-bromonaphthalene (1-BN) were studied. FTIR and FT Raman spectra were recorded in the region 4000-400 cm(-1) and 3500-50 cm(-1), respectively. The structural and spectroscopic data of the molecule in the ground state were calculated by using ab initio Hartree-Fock and density functional method (B3LYP) with the 6-311++G(d,p) basis set. The vibrational frequencies were calculated and scaled values were compared with experimental FTIR and FT Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The optimized geometric parameters were calculated. The predicted first hyperpolarizability also shows that the molecule might have a reasonably good nonlinear optical (NLO) behaviour. The calculated HOMO-LUMO energy gap reveals that charge transfer occurs within the molecule.     

FTIR and FT-Raman spectra, assignments, ab initio HF and DFT analysis of xanthine

The molecular vibrations of xanthine were investigated in polycrystalline sample, at room temperature by Fourier transform infrared (FTIR) and FT-Raman spectroscopies. The spectra of the molecule have been recorded in the regions 4000-50 cm(-1) and 3500-100 cm(-1), respectively. Theoretical information on the optimized geometry, harmonic vibrational frequencies, infrared and Raman intensities were obtained by means of ab initio Hartree-Fock (HF) and density functional theory (DFT) gradient calculations with complete relaxation in the potential energy surface using 6-311++G(d,p) basis set. The vibrational frequencies which were determined experimentally from the spectral data are compared with those obtained theoretically from ab initio and DFT calculations. A close agreement was achieved between the observed and calculated frequencies by refinement of the scale factors. The infrared and Raman spectra were also predicted from the calculated intensities. Thermodynamic properties like entropy, heat capacity, zero point energy have been calculated for the molecule. Unambiguous vibrational assignment of all the fundamentals was made using the potential energy distribution (PED).     

Ab-initio molecular geometry and normal coordinate analysis of tetrahydrothiophene molecule.

The molecular geometry of tetrahydrothiophene (THT) was quantum mechanically calculated using the split valence 6-31G** basis set. Electron correlation energy has been computed employing MP2 method. The molecule showed a twist form puckered structure with a twist torsion angle of 13 degrees and has a total energy of -347,877.514 kcal/mol of which a 436.715 kcal/mol electron correlation energy. The envelope form of the molecule showed an inter-plane angle of 22 degrees and has a total energy of -347,874.430 kcal/mol involving -436.558 kcal/mol electron correlation energy. The normal coordinates of the molecule were theoretically analyzed and the fundamental vibrational frequencies were calculated. The IR and laser Raman spectra of THT molecule was measured. All the observed vibrational bands including combination bands and overtones were assigned to normal modes with the aid of the potential energy distribution values obtained from normal coordinate calculations. The molecular force field was determined by refining the initial set of force constants using the least square fit method instead of using the less accurate scaling factor methods. The determined molecular force field has produced simulated frequencies which best match the observed values. The lowest-energy modes of vibration were two molecular out-of-plane deformations, observed at 114 and 166 cm(-1). The barrier of ring twisting estimated from the observed ring out-of-plane vibrational mode at 114 cm(-1) was estimated.     

Ab initio and DFT studies of the molecular structures and vibrational spectra of succinonitrile

The Molecular structure, conformational stability and vibrational frequencies of succinonitrile NCCH2CH2CN have been investigated with ab initio and density functional theory (DFT) methods implementing the standard 6-311++G* basis set. The potential energy surfaces (PES) have been explored at DFT-B3LYP, HF and MP2 levels of theory. In agreements with previous experimental results, the molecule was predicted to exist in equilibrium mixture of trans and gauche conforms with the trans form being slightly lower in energy. The vibrational frequencies and the corresponding vibrational assignments of succinonitrile in both C2h and C2 symmetry were examined theoretically and the calculated Infrared and Raman spectra of the molecule were plotted. Observed frequencies for normal modes were compared 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.     

Simulations of vibrational spectra from classical trajectories: calibration with ab initio force fields

An algorithm allowing simulating vibrational spectra from classical time-dependent trajectories was applied for infrared absorption, vibrational circular dichroism, Raman, and Raman optical activity of model harmonic systems. The implementation of the theory within the TINKER molecular dynamics (MD) program package was tested with ab initio harmonic force fields in order to determine the feasibility for more extended MD simulations. The results suggest that sufficiently accurate frequencies can be simulated with integration time steps shorter than about 0.5 fs. For a given integration time step, lower vibrational frequencies ( approximately 0-2000 cm(-1)) could be reproduced with a higher accuracy than higher-frequency vibrational modes (e.g., O-H and C-H stretching). In principle, the algorithm also provides correct intensities for ideal systems. In applied simulations, however, the intensity profiles are affected by an unrealistic energy distribution between normal modes and a slow energy relaxation. Additionally, the energy fluctuations may cause weakening of the intensities on average. For ab initio force fields, these obstacles could be overcome by an arbitrary normal mode energy correction. For general MD simulations, averaging of many shorter MD trajectories started with randomly distributed atomic velocities provided the best spectral shapes. alpha-pinene, D-gluconic acid, formaldehyde dimer, and the acetylprolineamide molecule were used in the tests.     

2-氟-5-溴吡啶分子振动光谱的密度泛涵理论研究

用密度泛函理论方法B3LYP以及6-311++G(2df,2pd)基组对2-氟-5-溴吡啶分子 的平衡几何构型进行了优化并计算了该分子的振动谐力场。使用Pulay的标度方法 对理论力场进行了标度。采用Wilson的GF矩阵方法,根据标度后的理论力场进行了 简正坐标分析,对2-氟-5-溴吡啶分子的振动基频进行了理论研究,得到了势能分 布和红外振动频率。与红外频率的实验值相比较,理论频率的均方差为24 cm~(-1) 。此外,根据振动模式的势能分布对此分子的振动基频进行了理论归属,并对前人 的指认进行了修正和补充。     

Si2分子基态和低激发态的电子结构

本文用abinitioSOCI／MCSCF／／HF／DH＋d方法计算得到了Si2分子的基态和八个激发态的平衡几何构型、总能量、谐振动频率、零点能以及各个态的势能曲线．讨论了各态键长随电子组态的变化关系，指出了成健轨道中。电子数的增加，是引起键长及谐振动频率变化的重要因素·其中3Πg，5∑；3△u，3∑，未见有关文献报导由于硅藻合物的重要应用价值，人们已对其进行了一些实验和理论研究h，’」．以前的研究主要是讨论其成键性质．由于St。分子结构的特殊性，要想全面了解其成键性质，不仅要研究其基态的电子结构，而且还要研究其激发态的电…     

A theoretical investigation on the geometry and vibrational spectra of 10,10,2,6,5-pentamethyl-1-hydroxychroman: a model of alpha-tocopherol

In the present study, density functional theory calculations with the combined Becke's three-parameter exchange functional in combination with the Lee, Yang, and Parr correlation functional (B3LYP) exchange-correlation energy functions were performed by using the 6-311G** basis set to study the structure and vibrational spectra of 10,10,2,6,5-pentamethyl-1-hydroxychroman (a model of alpha-tocopherol). The fully optimized geometry of the molecule was found to be very consistent with the X-ray crystal structure. The predicted vibrational frequencies made it possible to give a reliable assignment of the IR spectrum of the molecule according to the potential energy distributions (PEDs).     

Quantum chemical studies of FT-IR and Raman spectra of methyl 2,5-dichlorobenzoate

In this paper, experimental and theoretical studies on the molecular structure and vibrational spectra of methyl 2,5-dichlorobenzoate (MDCB) are presented. Fourier transform infrared and Raman spectra of the title molecule in the solid phase were recorded and analyzed. The geometrical parameters were calculated using DFT (B3LYP) with 6-311G(d,p) and 6-311++G(d,p) basis sets, and compared with the experimental data. The vibrational frequencies, infrared intensities and Raman scattering activities were also reported. The detailed assignments were given based on the total energy distribution of the vibrational modes, calculated with scaled quantum mechanics method. The observed and calculated frequencies are found to be in good agreement.     

Molecular structure and vibrational assignment of 2-,4-,6-methylquinoline by density functional theory (DFT) and ab initio Hartree-Fock (HF) calculations

The molecular geometry, the normal mode frequencies and corresponding vibrational assignments of 2-,4-,6-methylquinoline (2-,4-,6-mq) in the ground state were performed by HF and DFT/B3LYP levels of theory using the 6-31++G(d,p) basis set. Harmonic and anharmonic vibrational frequencies were calculated. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method by using parallel quantum mechanic solutions program. The general agreements between the observed and calculated frequencies are shown.     

A Dodecahedrane-like Molecule C_（12）H_（12）B_8 with Uncommon T_h Symmetry

The molecule with Th symmetry is rare.A dodecahedrane-like molecule C12H12B8 with uncommon Th symmetry has been reported here.Density functional calculations and minimization techniques have been employed to characterize its structural and electronic properties.Its geometry,electronic properties,vibrational frequencies and heat of formation have been calculated at the B3LYP/6-311+G(d,p) level of theory.The absence of imaginary vibrational frequency confirms that it corresponds to true minimum on the potential energy hypersurface.Its vibrational bands in the IR intensity have been discussed and compared with future experimental identification.At the B3LYP/6-311+G(d,p) level,the heat of formation has been calculated to be 720.9 kJ mol-1 using the isodesmic reaction.According to this value,it is a potential high energy density molecule.     

FT-IR,FT-Raman,NMR spectra and DFT calculations on 4-chloro-N-methylaniline

In this work, the vibrational spectral analysis was carried out by using FT-IR and FT-Raman spectroscopy in the range 400–4000 and 50–3500 cm^?1 respectively, for the title molecule. The structural and spectroscopic data of the molecule in the ground state were calculated by using density functional method using 6-311++G(d,p) basis set. The vibrational frequencies were calculated and scaled values were compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments of all the vibrational mode were performed on the basis of the total energy distributions (TED). ¹³C and ¹H NMR chemical shifts results were given and are in agreement with the corresponding experimental values. The theoretically constructed FT-IR and FT-Raman spectra exactly coincides with experimental one.     

5,6]-Heterofullerene-like C58Ge: odd atoms assembling a cage

Density functional calculations and structural minimization techniques have been employed to characterize the structural and electronic properties of [5,6]-heterofullerene-like C(58)Ge. The heterofullerene molecule has an odd number of atoms on the skeleton of the cage and is a novel molecule. The vibrational frequencies of the molecule have been calculated at the B3LYP/3-21G level of theory. The absence of imaginary vibrational frequency confirms that the molecule corresponds to a true minimum on the potential energy hypersurface. Its heat of formation was estimated in this study. Since the symmetry of the molecule of [5,6]-heterofullerene-like C(58)Ge is C(2), it is a chiral molecule.     

Ab-initio molecular geometry and normal coordinate analysis of pyrrolidine molecule.

The molecular geometry of pyrrolidine was quantum mechanically calculated using the split valence 6-31G** basis set. Electron correlation energy has been computed employing MP2 method. The molecule showed an envelope form puckered structure with inter-plane angle of 36.4 degrees and has a total energy of -132976.80 kcal mol(-1) of which a -464.86 kcal mol(-1) electron correlation energy. The twist form of the molecule showed a twist angle of 10.2 degrees from planarity and has a total energy of -132976.05 kcal mol(-1) involving -464.097 kcal mol(-1) electron correlation energy. The normal coordinates of the molecule were theoretically analyzed on the basis of the Cs point symmetry of the envelope form. Using initial set of force constants obtained from the ab-initio calculations the fundamental vibrational frequencies were computed. The IR and laser Raman spectra of Pyrrolidine molecules were measured. All the observed vibrational bands including combination bands and overtones were assigned to normal modes with the aid of the potential energy distribution values obtained from normal coordinate calculations. The molecular force field was obtained by refining the initial set of force constants using the least square fit method. The molecular force field was determined by refining the initial set of force constants using the least square fit method instead of using the less accurate scaling factor methods. The determined molecular force field has produced simulated frequencies best match to the observed values. The low frequency molecular out-of-plane deformation modes were observed in both infrared and Raman spectra at 298 and 163 cm(-1). The barrier of ring twisting estimated from the observed ring out-of-plane vibrational mode at 163 cm(-1) was found 3.1 kcal mol(-1).     

RDX分子内振动能量重分配的动力学研究

基于从头算分子动力学(Born-oppenheimer molecular dynamics, BOMD)模拟, 构建了环硝胺六氢-1,3,5-三硝基-1,3,5-三嗪(RDX)单分子不同振动模式之间的耦合矩阵, 并计算了在不同加载能量下从低频振动模式到高频振动模式的最优能量传输路径. 结果表明, RDX单分子中—NNO₂基团更有利于能量局域化, 振动模式v₃和v₄在从低频振动模式到高频振动模式的能量传输过程中扮演着重要角色. 通过对v₃和v₄两个振动模式的进一步分析发现, 加载能量的不同会导致RDX单分子能量传输路径的不同. 当加载能量较低时, RDX单分子倾向于从低频振动模式到中频振动模式再到高频振动模式的能量传输路径; 当加载能量较高时, 能量更倾向于从低频振动模式直接传输到高频振动模式上. 揭示了RDX分子内振动耦合能量转移的微观机制, 为进一步探索RDX将“机械能”转化为“化学能”的微观过程提供了理论基础.     

Intramolecular vibrational energy redistribution involving the torsion in CF3CH3: a molecular dynamics study

Classical trajectory calculations on intramolecular vibrational energy redistribution (IVR) involving the torsion in 1,1,1-trifluoroethane (TFE) are reported. Two potential energy functions (PEFs) are used to describe the potential energy surface. The "full" PEF gives excellent agreement with the experimental vibrational frequencies. The "simple" PEF omits nondiagonal interaction terms, but still gives very good agreement with the experimental frequencies. The "simple" PEF is intended to minimize mode-mode coupling. Neither PEF includes the HF elimination reaction. Calculations are carried out both with nominal microcanonical selection of initial coordinates and momenta, and with a modified selection method that places controlled amounts of energy in the torsion. Total (classical) vibrational energies from 0.005 to 140 kcal mol(-1) are investigated. The calculated time constants describing energy flow out of the torsional mode are <10 ps for classical vibrational energies near the classical reaction threshold energy (approximately 75 kcal mol(-1)) and greater. It is found that the rate of decay from the torsion largely depends on the amount of energy in the other vibrational modes. Analysis using power spectra shows that the torsional mode in TFE is strongly coupled to the other vibrational modes. These results strongly suggest that vibrational energy in TFE will not be sequestered in the torsion for time periods greater than a few tens of picoseconds when the molecule has enough energy to react via HF elimination.     

Three rotor potential energy scans, conformational equilibrium constants and vibrational analysis of 3-fluoro-1-propanol CH(2)FCH(2)CH(2)OH

The conformational stability and the three rotor internal rotations in 3-fluoro-1-propanol were investigated by the DFT-B3LYP/6-311+G** and the ab initio MP2/6-311+G** levels of theory. The calculated potential energy curves of the molecule at both levels of theory were consistent with complex conformational equilibria of about 12 minima, all of which were predicted to have real frequencies at both the B3LYP and the MP2 levels. The lowest energy minimum in the potential curves of 3-fluoro-1-propanol was predicted to correspond to the Gauche-gauche-trans (Ggt) conformer in excellent agreement with microwave and electron diffraction results. The equilibrium constants for the conformational interconversion of the molecule were calculated and found to correspond to an equilibrium mixture of about 33% Ggt, 14% Ggg1 and 13% Gg1g and about 43% Ggt, 12% Ggg1 and 10% Gg1g distribution by the B3LYP/6-311+G** and the MP2/6-311+G** calculations, respectively, at 298.15K. The vibrational frequencies of each molecule in its three stable forms were computed at B3LYP level and complete vibrational assignments were made based on normal coordinate calculations and comparison with experimental data of the molecule.     

Optimal internal coordinates, vibrational spectrum, and effective Hamiltonian for ozone

In this paper the authors use the optimal internal vibrational coordinates previously determined for the electronic ground state of the ozone molecule to study the vibrational spectrum of the molecule employing the second empirical potential energy surface calculated by Tyuterev et al. [Chem. Phys. Lett. 316, 271 (2000)]. First, the authors compute variationally all the bound vibrational energy levels of the molecule up to the dissociation limit and state the usefulness of the optimal coordinates in this respect, which allows us to converge all the bound levels using relatively small anharmonic basis sets. By analyzing the expansion coefficients of the wave functions, they show then that a large portion of the vibrational spectrum of O3 can be structured in nearly separable polyadic groups characterized by the polyad quantum number N=n1+n2+n(theta) corresponding to the optimal internal coordinates. Accordingly, they determine an internal effective vibrational Hamiltonian for O3 by fitting the effective Hamiltonian parameters to the experimental vibrational frequencies, using as input parameters in the fit those extracted from an analytical second-order Van Vleck perturbation theory calculation. It is finally shown that the internal effective Hamiltonian thus obtained accurately describes the vibrational spectrum of ozone in the low and medium energy regimes.     

IR and Raman spectra,density functional computations of vibrational spectrum,molecular geometry,atomic charges,and some molecular properties of 3‐aminobenzonitrile molecule

Geometry, vibrational wavenumbers, and several thermodynamic parameters have been calculated using ab initio quantum chemical methods for the 3‐aminobenzonitrile molecule for the first time. The results were compared with experimental values. With the help of specific scaling procedures, the observed vibrational wavenumbers were analyzed and assigned to different normal modes of the molecule. In general, the error obtained was very low. Using potential energy distribution (PED), the contributions of the different modes to each wavenumber were determined. Other general conclusions were also deduced. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006