Ab initio study of spectroscopic constants and anharmonic force field of 74GeCl2

The equilibrium structure, spectroscopy constants, and anharmonic force field of germanium dichloride have been calculated at MP2, B3LYP, and CCSD(T) levels of theory employing two basis sets, cc-pVDZ and cc-pVTZ, respectively. The computed geometries, rotational constants, and vibration-rotation interaction constants, and quartic centrifugal distortion constants are compared with the available experimental data. The harmonic frequencies, anharmonic constants, and cubic and quartic force constants are predicted. The calculated results show that the MP2 results are in excellent agreement with experiment and represent a substantial improvement over the results obtained from B3LYP. The CCSD(T) method is also an advisable choice to study anharmonic force field of molecules.     

Ab initio SCF computation of force constants for CO2

Force constants for CO₂ have been evaluated using SCF wave functions. The effect ofd basis functions and geometry are investigated. Comparison with experimental values shows that a large error, due to neglect of electron correlation, occurs for theK ₁₂ interaction stretch force constant.     

Ab initio simulation of benzene Raman intensities

C₆H₆ Raman scattering activities calculated from harmonic model ab initio Hartree–Fock 6–311 ++ G(d, p) polarizability derivatives (and harmonic force fields calculated at the same level) accurately simulate experiment (to within 1% for the a_1g modes). Accurate predictions are also made for the e_2g modes (to within 5% for ν₇ and ν₉, and more poorly for ν₆ and ν₈ [in Fermi resonance with ν₆ + ν₁]) and for the e_1g out-of-plane mode, ν₁₀. Only the ν₆ in-plane CCC bending mode scattering activity is found to be anomalous. Systematic variation of the basis set indicates that the benzene scattering activities and depolarization ratios are strongly dependent on inclusion of both carbon and hydrogen atom diffuse functions in the basis set. Predictions are also made for ¹²C₆D₆ and for unmeasured intensities in ¹³C₆H₆. Measurements of a_1g mode scattering activities in the latter molecule are predicted to be useful in testing the harmonic Hartree–Fock Raman intensity model.     

Ab initio study of the NMR shielding constants and spin-spin coupling constants in cyclopropene

Summary Ab initio calculations of parameters which characterize the NMR spectrum are presented for the cyclopropene molecule. The London orbitals CHF (or GIAO-CHF, Gauge-Independent Atomic Orbital Coupled Hartree-Fock) results for the shielding constants are in good agreement with the experimental data, accurately determined, and with otherab initio values. The calculations of the NMR spin-spin coupling constants have been performed using the Multiconfiguration Time-Dependent Hartree-Fock (MC TDHF) approach. Different basis sets and MC SCF wavefunctions were used to estimate the accuracy of the results. Good agreement is obtained with the coupling constants estimated using the available experimental data.Dedicated to Professor Werner Kutzelnigg on the occasion of his 60th birthday     

Ab initio calculations of the vibrational force field and IR intensities of the ammonia dimers

Ab initio SCF calculations using the 4-31G basis set have been carried out to determine the equilibrium geometry, force constants and dipole moment derivatives of the linear (C_s and cylic (C_2h) ammonia dimers. The results are compared with monomer calculations and experimental data.     

Ab initio MO calculation of force constants and dipole derivatives for the formamide dimer. An estimation of hydrogen-bond force constants

Ab initio SCF MO calculations (with the 4-31G basis set) have been carried out to determine the equilibrium geometry, vibrational frequencies, dipole-moment derivatives, and force constants for intermolecular modes of the formamide dimer and its d₄ and d₆ derivatives. The results are correlated with monomer calculations and experimental data for crystalline formamide.     

Ab initio vibrational transition dipole moments and intensities of formaldehyde

Vibrational transition dipole moments and absorption band intensities for the ground state of formaldehyde, including the deuterated isotopic forms, are calculated. The analysis is based on ab initio SCF and CI potential energy and dipole moment surfaces. The formalism derives from second-order perturbation theory and involves the expansion of the dipole moment in terms of normal coordinates, as well as the incorporation of point group symmetry in the selection of the dipole moment components for the allowed transitions. Dipole moment expansion coefficients for the three molecule-fixed Cartesian coordinates of formaldehyde are calculated for internal and normal coordinate representations. Transition dipole moments and absorption band intensities of the fundamental, first overtone, combination, and second overtone transitions are reported. The calculated intensities and dipole moment derivatives are compared to experiment and discussed in the context of molecular orbital and bond polarization theory.     

Ab initio computation of force constants: Part VI. Applications of the force relaxation method for geometry optimization

Techniques to improve the computational efficiency of the force relaxation method are discussed. Force constants for fragments in previously computed smaller molecules can be transferred to construct a guess force constant matrix. Additional force constants that may be needed can be computed by a procedure which uses only one additional force calculation per diagonal force constant required. A scaling technique to improve convergence on the optimized geometry is discussed.     

Ab initio and DFT predictions of infrared intensities and Raman activities

Relative infrared (IR) intensities and relative Raman activities have been computed for vibrations of test molecules, including from two to nine heavy atoms, using second-order Moller-Plesset perturbation theory (MP2), and three hybrid density functionals (B3LYP, M05, and M05-2X). The basis set convergence of vibrational properties is discussed. Our results demonstrate that B3LYP offers the most cost-effective choice for the prediction of molecular vibrational properties, but the predictions of another two tested hybrid functionals are very similar and in very good agreement with experimental data. MP2 shows good performance for the IR intensities, whereas the quality of prediction of the relative Raman activities should be characterized as only moderate. B3LYP calculations of the relative IR intensities using highly compact Sadlej's Z3PolX basis set retain the high accuracy of the more CPU expensive Sadlej's pVTZ and much more expensive aug-cc-pVTZ calculations. Relative Raman activities are more sensitive to basis set effects and require at least Sadlej's pVTZ to obtain quantitative results.     

Ab initio based polarizable force field parametrization

Experimental and simulation studies of anion-water systems have pointed out the importance of molecular polarization for many phenomena ranging from hydrogen-bond dynamics to water interfaces structure. The study of such systems at molecular level is usually made with classical molecular dynamics simulations. Structural and dynamical features are deeply influenced by molecular and ionic polarizability, which parametrization in classical force field has been an object of long-standing efforts. Although when classical models are compared to ab initio calculations at condensed phase, it is found that the water dipole moments are underestimated by approximately 30%, while the anion shows an overpolarization at short distances. A model for chloride-water polarizable interaction is parametrized here, making use of Car-Parrinello simulations at condensed phase. The results hint to an innovative approach in polarizable force fields development, based on ab initio simulations, which do not suffer for the mentioned drawbacks. The method is general and can be applied to the modeling of different systems ranging from biomolecular to solid state simulations.     

Ab initio study of interaction-induced NMR shielding constants in mixed rare gas dimers

The interaction-induced contribution to the NMR shielding constants in homonuclear A2 and heteronuclear AB (A,B=He,Ne,Ar) dimers is obtained ab initio by employing a coupled cluster singles and doubles with perturbative treatment of triples wave function model and extended correlation-consistent basis sets. The second virial coefficients entering the expansion of the property with the density are then computed in a fully quantum mechanical approach, for temperatures ranging from the limit of dissociation of the dimer to well above standard conditions. The results can be used to describe the density and temperature dependence of the shielding constants in binary mixtures of helium, neon, and argon. The predicted effects should be observable for the interaction of 21Ne with other rare gases.     

Reaction of triphenyl borate with 1,3,5-trioxane

Triphenyl borate was prepared by reaction of boric acid with phenol in xylene. Its reaction with 1,3,5-trioxane involved replacement of protons in the para positions of the benzene rings by methylene group.     

Ab initio study of KN

The potential energy curves for the lowest (3)Sigma(-), (3)Pi, and (5)Sigma(-) states of the KN molecule have been calculated by the multireference singles and doubles configuration interaction method, including Davidson's corrections for quadruple excitations [MRCI(+Q)]. It is shown that the former two are bound, while the last one is repulsive. The electronic ground state of KN is predicted as (3)Sigma(-) state, although the term energy of the (3)Pi state is very small, 177.3 cm(-1). The binding energy for the (3)Sigma(-) state is evaluated as 0.838 eV, the rotational constant B(0) as 0.250 63 cm(-1), and harmonic frequency as 324.4 cm(-1). The spin-orbit coupling effects between the (3)Sigma(-) and (3)Pi states of KN are evaluated and discussed. The same MRCI(+Q) computational procedures are applied to the isovalent LiN, KC, KO, and KCl to confirm the accuracy of present calculations. Theoretical spectroscopic constants presented here will inspire experimental studies of KN.     

Anharmonic force field and spectroscopic constants of silene: an ab initio study

High-level ab initio calculations with large basis sets are reported for silene, H₂C=SiH₂. Correlated harmonic force fields are obtained from coupled cluster CCSD(T) calculations with the cc-pVQZ basis (cc-pVTZ for H) while the anharmonic force fields are computed at the MP2/TZ2Pf level. There is excellent agreement with the available experimental data, in particular the equilibrium geometry and the fundamental vibrational frequencies. Many other spectroscopic constants are predicted for the C _{2 v} isotopomers of silene. Received: 27 May 1998 / Accepted: 23 July 1998 / Published online: 9 October 1998     

Ab initio rotational constants of interstellar species: Cyanoacetylene hydrogenated derivatives

By means of ab initio calculations, the rotational constants and dipole moments of H_nC₃N (n = 1, 3, 5, 7, and 9) species have been calculated at the HF / 6-31G * level of theory. Selected cases have been also calculated at the MP 2/6-31G * level and the influence of calculation level on rotational constant values is briefly discussed. Some of these species were discovered in the interstellar medium, while others have still not been detected there, although their existence is very probable. The results given here could help in their detection. © 1993 John Wiley & Sons, Inc.     

Ab initio study of hydrazoic acid

SCF and CI calculations were carried out on the ground^1A state of HN₃. The equilibrium geometry and vibration frequencies were computed. The results point to a planar structure (groupC _s) but to a non-linear (170 °) N-N-N conformation. The calculated vibration frequencies are in fair agreement with experimental assignments.The dissociation path of the molecule to NH and N₂ products was investigated and compared to the isoelectronic reaction of diazomethane. The dissociation energy of hydrazoic acid is estimated to be about –8 kcal/mole, with a potential barrier to dissociation of about 30 kcal/mole.Boursier IRSIA     

Ab initio study of cyclic hexahalotriphosphazenes

Using standard ab initio methods, the electronic structure and optimal geometries of cyclo-(NPX₂)₃ (X = F, Cl) are investigated at the DZ+P basis set level. Out-of-plane π overlap populations for P-N bonds (evaluated as the contribution of a“₂ and e” molecular orbitais) are roughly two times greater than the in-plane gp ones evaluated as half the difference between the gross and net d orbital Mulliken populations on phosphorus belonging to é molecular orbitais.     

Ab initio force field for simulations of proteins and nucleic acids

We present a new self-consistent set of ab initio analytical pair potential to predict specific nonbonded interactions of protein with nucleic acid, of protein with protein, and of nucleic acid with nucleic acid. The purpose of this study is to represent the interaction between biological molecules with an accuracy equivalent to the ab initio molecular orbital calculations, which are used as reference data to obtain the pair potentials. Atoms in nucleic acids and proteins are classified according to their chemical environments. An “effective charge,” a modification of a charge obtained from the Mulliken population analysis, is introduced and used to represent the electrostatic energy. More than 30,000 SCF interaction energies have been calculated to provide the reference data for the fitting procedure that we have adopted in the parameterization of the potentials. The standard deviation is 1.61 kcal/mol for interaction energies spanning the range from about ?220 kcal/mol to +20 kcal/mol. Molecular dynamics simulations, using the above new set of force field, have been performed successfully for the systems where adequate treatments of specific interactions are required: The stability of α-helix of C-peptide and the interaction of spermine with oligonucleotide are examined as preliminary examples.