Scaled quantum chemical force fields for 1,1-difluorocyclopropane and the influence of vibrational anharmonicity

Potential functions and harmonic (omega(i)) and anharmonic (nu(i)) fundamental frequencies have been calculated for 1,1-difluorocyclopropane (DFCP) and its d4 and d2 isotopomers using the program Gaussian 03. B3LYP and MP2 models were employed, each with the bases 6-311++G** and cc-pVTZ. Anharmonicity corrections Delta(i) = omega(i) - nu(i) are listed and shown to be different for symmetric and antisymmetric CH stretching modes in situations where Fermi resonance appears to be absent. The same effect is missing in C2H4, for which similar calculations were made. The quadratic force fields for DFCP have been scaled in symmetry coordinate space with 15 scale factors both to observed frequencies nu(obsd)and also to omega (obsd), where omega(obsd) = nu(obsd) + Delta. With nu(obsd) especially, different scale factors are needed for the symmetric and antisymmetric CH stretching force constants due to their differing anharmonicities. The source of the latter in the quartic and cubic force field is explored. MP2 calculations of valence interaction force constants involving the stretching of bonds on a common carbon atom are preferred to those from a B3LYP model. In either model, scaling to omega(obsd) rather than to nu(obsd) does not remove the necessity of varying scale factors for differing types of motion in the same group. Theoretical values of the five quartic centrifugal distortion constants are listed for the normal species and compared with new experimental data. The predictions are sufficiently good to be useful in fitting pure rotational transitions. A weakness is identified in the current Gaussian 03 code for the calculation of vibration-rotation quantities, and limitations are noted in the manner in which Fermi resonance is handled.     

Acetamide, a challenge to theory and experiment? On the molecular structure, conformation, potential to internal rotation of the methyl group and force fields of free acetamide as studied by quantum chemical calculations

The molecular geometry has been optimized without any constraints using different basis sets and levels of theory as: Hartree-Fock with basis sets 6–31+G**, 6–311++G**, cc-pVTZ and aug-cc-pVTZ, MP2 with basis sets 6–311++G** and cc-pVTZ, MP3 with basis set 6–311++G**, and density functional theory with basis sets 6–311++G** and cc-pVTZ. Small basis sets up to 6-31G predict the syn conformation of the methyl group to be the most stable conformation. Larger basis sets predict an unsymmetrical conformation with one of the H atoms perpendicular to the amide skeleton or an anti-like conformation. Dunnings correlation consistent polarized valence triple zeta, cc-pVTZ, basis set including MP2 predict two conformations, one perpendicular and one anti to be the most stable. The DFT calculations predict anti-like conformations. The most accurate calculations predict anti-like conformations which have not been predicted previously. The vibrational frequencies have been calculated for several basis sets and compared to the observed frequencies. The wagging frequency of the NH₂ is very dependent on the basis sets and levels of theory. Most calculations predict a planar NH₂ group in agreement with experiment. A scaled molecular force field has been determined by fitting the calculated frequencies to the observed ones for the perpendicular conformation using MP2/cc-pVTZ. The barrier heights for the methyl group have been calculated. The rotational constants, I_A + I_B − I_C values and dipole moments are compared with experimental values.     

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

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

Anharmonic Effect of the Unimolecular Dissociation of Ethyl Radical

The anharmonic and harmonic rate constants have been calculated for the unimolecular dissociation of ethyl radical using the method proposed by Yao and Lin (YL method) at both B3LYP/6‐311++G** and MP2/6‐311++G** levels. The different rate constants indicate that the results obtained from B3LYP and MP2 method are very close. The anharmonic and tunneling effect of the title reaction has also been examined. The comparison shows that, both in microcanonical and canonical systems, the anharmonic rate constants are higher than those for harmonic cases, especially in the case of high total energies and temperatures, which indicates that anharmonic effect of the unimolecular dissociation of ethyl into C₂H₄ and H is so significant that cannot be neglected. The tunneling effect is very small for the decomposition of C₂H₅ radical.     

Structural studies on ethyl isovalerate by microwave spectroscopy and quantum chemical calculations

We observed the microwave spectrum of ethyl isovalerate by molecular beam Fourier transform microwave spectroscopy. The rotational and centrifugal distortion constants of the most abundant conformer were determined. Its structure was investigated by comparison of the experimental rotational constants with those obtained by ab initio methods. In a first step, the rotational constants of various conformers were calculated at the MP2/6-311++G** level of theory. Surprisingly, no agreement with the experimental results was found. Therefore, we concluded that in the case of ethyl isovalerate more advanced quantum chemical methods are required to obtain a reliable molecular geometry. Ab initio calculations carried out at MP3/6-311++G**, MP4/6-311++G**, and CCSD/6-311++G** levels and also density functional theory calculations using the B3LYP/6-311++G** method gave similar results for the rotational constants, but they were clearly distinct from those obtained at the MP2/6-311++G** level. With use of these more advanced methods, the rotational constants of the lowest energy conformer were in good agreement with those obtained from the microwave spectrum.     

Systematic ab initio studies of the conformers and conformational distribution of gas-phase tyrosine

A full structural assignment of the conformers of gaseous tyrosine is presented. A total of 1296 unique trial structures were generated by allowing for all combinations of internal single-bond rotamers and optimized at the B3LYP6-311G* level of theory and then subjected to further optimization at the B3LYP6-311++G** level. A total of 76 conformers are found and their dipole moments, rotational constants, and harmonic frequencies are determined. Accurate relative energies are given at the MP26-311G(2df,p)B3LYP6-311++G** level of theory. Characteristic H-bonding types are classified and listed for all the conformers. The four most stable conformers display an intramolecular H bond, COOH...NH(2), and an additional H-bonding interaction between the amino group and pi electron of the aromatic ring. The results further confirm that the global minimum conformations of the aromatic amino acids have the same H-bonding type. Combined with statistical mechanics principles, conformational distributions at various temperatures are computed and the temperatures with which the theoretical results match that of experiments are indicated.     

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.     

Ab initio studies of the conformers and conformational distribution of the gaseous hydroxyamino acid threonine

Systematic and extensive conformational search has been performed to characterize the gas-phase threonine structures. A total of 1296 unique trial 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 of the theory and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 71 conformers were found and their rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies and vertical ionization energies of all the conformers were determined. Single-point energies were also calculated at the MP2/6-311G(2df,p) and B3LYP/6-311G(2df,p) levels. Characteristic H-bonding types were classified and listed for all the conformers. The conformational distributions of gaseous threonine at various temperatures were calculated.     

Electronic structure and vibrational spectra of X20H20 (X=C, Si, Ge, Sn): A theoretical study

Optimized geometries, HOMO–LUMO gaps, vertical ionization potentials and electron affinities are obtained using HF, and B3LYP methods with 6-311G** basis set for C₂₀H₂₀, Si₂₀H₂₀ and Ge₂₀H₂₀. For germanium and tin analogues, B3LYP calculations are performed with LANL2DZ effective core potential. Electron correlation is included by doing MP2 calculation. The harmonic frequencies of all the compounds are obtained using B3LYP with 6-311G** and/or LANL2DZ basis sets. The force field and vibrational spectra are analyzed and 74 symmetry unique non-redundant local force constants are evaluated. Probable assignments are proposed for all the fundamentals based on the potential energy distribution.     

Theoretical study of light-metal tetrahydroborates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroborates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the structural parameters, the energies of different decomposition pathways, barriers to internal rotation of BH₄ groups, normal mode frequencies, magnetic shielding constants, and spin density distribution (in radicals) have been analyzed in various related series of these compounds. The results obtained by ab initio methods and at the DFT level are compared. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** adequately reproduces the results obtained at the higher level of theory CCSD(T)/6-311++G*s*//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

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.     

Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.     

s-trans-1,3-butadiene and isotopomers: vibrational spectra, scaled quantum-chemical force fields, fermi resonances, and C-H bond properties

Quadratic quantum-chemical force fields have been determined for s-trans-1,3-butadiene using B3LYP and MP2 methods. Basis sets included 6-311++G, cc-pVTZ, and aug-cc-pVTZ. Scaling of the force fields was based on frequency data for up to 11 isotopomers, some of these data being original. A total of 18 scale factors were employed, with, in addition, an alteration to one off-diagonal force constant in the A(u) species. MP2 calculations without f functions in the basis perform badly in respect of out-of-plane bending mode frequencies. Centrifugal distortion constants and harmonic contributions to vibration-rotation constants (alphas) have been calculated. Existing experimental frequency data for all isotopomers are scrutinized, and a number of reassignments and diagnoses of Fermi resonance made, particularly in the nu(CH) region. The three types of CH bond in butadiene were characterized in terms of bond length and isolated CH stretching frequency, the latter reflecting data in the nu(CD) region. Broad agreement was achieved with earlier results from local mode studies. Differences in CH bond properties resemble similar differences in propene. A simplified sample setup for recording FT-Raman spectra of gases was applied to four isotopomers of butadiene.     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

Detailed ab initio studies of the conformers and conformational distributions of gaseous tryptophan

A systematic and extensive conformational search has been performed to characterize the gas-phase tryptophan structures. A total of 648 unique trial structures were generated by allowing for all combinations of internal single-bond rotamers and were optimized at the DFT/B3LYP/6-311G* level of theory. A total of 45 local minima conformers were found. Further optimization of the 45 conformers with B3LYP and MP2/6-311++G** did not produce meaningful structural change, and accurate geometries, dipole moments, rotational constants, harmonic frequencies, and relative energies were then determined. Combined with statistical mechanics principles, the conformational distributions of gas-phase tryptophan at different temperatures are shown. The results clearly support the conclusion drawn by Compagnon et al. that only one dominant isomer existed in the molecular beam at 85 K and add further evidence that the supersonic jet expansion or embedding helium droplets did not produce an equilibrium distribution.     

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.     

An exploration of conformational search of leucine molecule and their vibrational spectra in gas phase using ab initio methods

An extensive exploration of the conformational space has been carried out to characterize all possible gas phase structures of leucine. A total of 324 unique trial structures for canonical leucine were generated by considering all possible combinations of single bond rotamers. All trial structures were optimized at the B3LYP/6-311G* level of the DFT method. A total of 77 unique and stationary canonical conformers were found. Further, 15 most stable conformers were reoptimized at B3LYP/6-311++G** level and their respective relative energies, vertical ionization energies, hydrogen bonding patterns, rotational constants and dipole moments were calculated. A single point energy calculations for leucine conformers have also been done at both B3LYP/6-311++G(2df, p) and MP2/6-311++G(2df, p) levels. The good agreement between our estimates of rotational constants for two most stable conformers and available experimental measurements supports the reliability of the B3LYP/6-311++G** level of theory for describing the conformational behavior of leucine molecule. The proton affinity and gas phase basicity were also determined. Using the statistical approach, conformational distributions at various temperatures have also been performed and analyzed. Vibrational spectra were also calculated. It is also observed that zwitterions of leucine are not stable in gas phase.     

A coupled-cluster study of the structure and vibrational spectra of pyrazole and imidazole

Harmonic force fields were calculated at the corresponding optimized geometries for pyrazole and imidazole at the HF, B3LYP, MP2, CCSD and CCSD(T) levels using the 6-31G* basis set and at the HF and B3LYP levels using the cc-pVTZ basis set. The agreement between the calculated and experimental geometries by the CCSD and CCSD(T) methods was generally similar to that obtained with the B3LYP and MP2 methods. The force fields were scaled using one-scale-factor (1SF), 3SF and 7SF scaling schemes. The scale factors were varied with respect to the experimental frequencies. Using 7SF scaling, the root-mean-square (RMS) deviation of the calculated frequencies from the experimental frequencies by the HF, B3LYP, MP2, CCSD and CCSD(T) methods and the 6-31G* basis set was 16, 7, 13, 11 and 11 cm(-1), respectively. This shows that the B3LYP method is preferred for force field calculations over the perturbative MP2, CCSD and CCSD(T) methods. Using 1SF scaling, the CCSD(T) scale factor was 0.931, the highest among the five methods used but close to that obtained with the B3LYP method and the cc-pVTZ basis set with lower RMS deviation.     

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.     

气相中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)分析了三种复合物中氢键的电子密度拓扑性质.