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     

Theoretical study of anharmonic force field and spectroscopic constants for 1-chlorophosphaethene,CH2PCl and CD2PCl

The anharmonic force field and spectroscopy constants for CH₂PCl are determined using CCSD(T), VPT2, and density functional theory employing cc-pVQZ basis sets. The molecule structure, rotational spectroscopic constants, and vibrational wave numbers are compared with the available experimental data. Anharmonicity constants, vibration-rotation interaction constants and cubic force constants are predicted. Vibrational wave numbers and rotational constants for CD₂PCl are also determined using the same levels. The isotopic shifts of vibrational wave numbers are remarkable by D atom substitution for 1-chlorophosphaethene.     

Ab initio study of force constants and intensities of 1,3,5-trioxane

The structure, vibrational force constants and infrared and Raman intensities have been calculated for 1,3,5-trioxane using a 3-21G basis set. These results have been used to identify some possible inaccuracies in experimental diffraction based structures and in vibrational assignments. It is demonstrated that there is a marked contrast between the trends in the vibrational force constants and in the King's effective atomic charges of the axial and equatorial CH bonds in the series cyclohexane, 1,4-dioxane and 1,3,5-trioxane. The axial CH stretching force constant decreases by 0.04 mdyne Å⁻¹ for each adjacent oxygen atom, whereas that of the equatorial CH bond increases by 0.15 mdyne Å⁻¹ per oxygen. In trioxane the effective atomic charge of the axial hydrogen is twice that of the equatorial. Atomic polar tensors are calculated in a bond oriented frame, and the effect of the oxygens on CH stretching and bending mode intensities discussed. Some properties are also calculated using the 4-31G basis.     

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 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.     

Calculation of anharmonic force constants

Non-linear transformations from internal to cartesian displacements are considered for anharmonic calculation of vibrational frequencies. Force constants in normal coordinates, up to quartic terms, are related to Christoffel symbols. The latter are tabulated for valence angles or, for torsion and wagging out of plane, evaluated through expressions which avoid differentiation of explicit functional forms for each type of internal coordinate. The same symbols are used for cartesian tensors and redundancy coefficients up to the third order.     

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.     

Ab initio study of azolides: Energetic and spectroscopic properties

We report the ab initio study of twenty‐four azolides derived from pyrrole, imidazole, pyrazole, both triazoles, tetrazole, pentazole, indole and carbazole bearing at the nitrogen atom the groups COMe, CHO, COCF₃ and CO₂Me. Theoretical values (isomerism, barriers, dipole moments, C=O stretching) are compared with experimental ones, when available, and also internally compared. A special effort has been devoted to the calculation of the absolute shieldings for the different nuclei present in azolides. At the level of calculation used (RHF/6‐311G**) the results are satisfactory. To complete the nmr data from the literature, some ¹H, ¹³C, ¹⁵N, ¹⁷O and 19F chemical shifts have been determined.     

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 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.     

SiH2Cl2: ab initio anharmonic force field, dipole moments, and infrared vibrational transitions

The vibrational spectra of SiH2Cl2 have been recorded in the 1000-13,000 cm(-1) region, utilizing the Fourier-transform spectroscopy and Fourier-transform intracavity laser absorption spectroscopy. Totally 61 band centers and intensities are derived from the infrared spectra. An ab initio quartic force field is obtained by applying the second-order Moller-Plesset perturbation theory and correlation-consistent polarized valence triplet-zeta basis sets [J. Chem. Phys. 90, 1007 (1989); 98, 1358 (1993)]. Most observed bands are assigned by the vibration analysis based on the second-order perturbation theory. Reduced-dimensional ab initio dipole moment functions (two dimensional and three dimensional) have also been calculated to investigate the absolute band intensities of the SiH2 chromophore. The calculated values agree reasonably with the observed ones.     

Ab initio study of the vibrational assignment and force field of thiosemicarbazide-d0 and -d5

The complete harmonic force field and optimized geometry of thiosemicarbazide have been calculated at the ab initio Hartree—Fock level using the 3-21G basis set. On the basis of this, the frequencies of thiosemicarbazide-d₀ and -d₅ and their ¹⁵N isotopic molecules have been calculated. The calculated frequencies and their band assignments are utilized to critically examine our previous experimental assignments which were based on normal coordinate calculations. The theoretical IR and Raman intensities, together with qualitative experimental band intensities, are also presented.     

Ab initio calculations and Franck-Condon simulation of the absorption spectra of GeCl2 including anharmonicity.

Geometrical parameters, vibrational frequencies and relative electronic energies of the X1A1, ?3B1 and A1B1 states of GeCl2 have been calculated at the CCSD(T) and/or CASSCF/MRCI level with basis sets of up to aug-cc-pV5Z quality. Core electron correlation and relativistic contributions were also investigated. RCCSD(T)/ aug-cc-pVQZ potential energy functions (PEFs) of the X1A1 and ?3B, states, and a CASSCF/MRCl/aug-cc-pVQZ PEF of the A1B1 state of GeCl2 are reported. Anharmonic vibrational wavefunctions of these electronic states of GeCl2, obtained variationally using the computed PEFs, are employed to calculate the Franck-Condon factors (FCFs) of the ?-X and A-X transitions of GeCl2. Simulated absorption spectra of these transitions based on the computed FCFs are compared with the corresponding experimental laser-induced fluorescence (LIF) spectra of Karolczak et al. [J. Chem. Phys. 1993, 98, 60-70]. Excellent agreement is obtained between the simulated absorption spectrum and observed LIF spectrum of the ?-X transition of GeCl2, which confirms the molecular carrier, the electronic states involved and the vibrational assignments of the LIF spectrum. However, comparison between the simulated absorption spectrum and experimental LIF spectrum of the A-X transition of GeCl2 leads to a revision of vibrational assignments of the LIF spectrum and suggests that the X1A1 state of GeCl2 was prepared in the experimental work, with a non-Boltzmann vibrational population distribution. The X(0,0,1) level is populated over 4000 times more than expected from a Boltzmann distribution at 60 K, which is appropriate for the relative population of the other low-lying vibrational levels, such as the X(1,0,0) and X(0,1,0) levels.     

Ab initio based force field and molecular dynamics simulations of crystalline TATB

An all-atom force field for 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is presented. The classical intermolecular interaction potential for TATB is based on single-point energies determined from high-level ab initio calculations of TATB dimers. The newly developed potential function is used to examine bulk crystalline TATB via molecular dynamics simulations. The isobaric thermal expansion and isothermal compression under hydrostatic pressures obtained from the molecular dynamics simulations are in good agreement with experiment. The calculated volume-temperature expansion is almost one dimensional along the c crystallographic axis, whereas under compression, all three unit cell axes participate, albeit unequally.     

Ab initio parameterization of an all-atom polarizable and dissociable force field for water

A novel all-atom, dissociative, and polarizable force field for water is presented. The force field is parameterized based on forces, stresses, and energies obtained form ab initio calculations of liquid water at ambient conditions. The accuracy of the force field is tested by calculating structural and dynamical properties of liquid water and the energetics of small water clusters. The transferability of the force field to dissociated states is studied by considering the solvation of a proton and the ionization of water at extreme conditions of pressure and temperature. In the case of the solvated proton, the force field properly describes the presence of both Eigen and Zundel configurations. In the case of the pressure-induced ice VIII/ice X transition and the temperature-induced transition to a superionic phase, the force field is found to describe accurately the proton symmetrization and the melting of the proton sublattice, respectively.     

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 study of the structure,force field,and vibrational spectrum of the LiClO3 molecule

The geometrical structure and the vibrational spectrum of the LiClO₃ molecule are studied by the Hartree-Fock-Roothaan (HF) and configurational interaction (CI) methods taking into account single and double excitations and Davidson's correction for quadruple excitations. Double-zeta basis sets of Huzinaga-Dunning and McLean-Chandler complemented with polarization and diffuse functions are used. Potential surface sections corresponding to migration of the Li⁺ cation around the ClO ₃ ⁻ fragment are investigated. It was found that the LiClO₃ molecule has a single stable configuration of C_s symmetry with the bidentate coordination of the Li⁺ cation by the ClO ₃ ⁻ anion. The cyclic fragment is nonplanar (the dihedral angle θ(LiO₂Cl)=173°). The tridentate configuration of C_3v symmetry lies higher than the equilibrium configuration by 24.5 (HF) or 18.3 (CI) kJ/mole and is not an isomer. The ab initio force field of the molecule was refined by the scaling method. Some assignments of the IR bands of the matrix-isolated molecular forms existing in vapor over lithium chlorate are corrected. The vibration frequencies (cm⁻¹) and IR intensities (km/mole; in parentheses) are calculated with the refined force field: A′ type 1099(236), 856(81), 630(73), 557(119), 481(87), 156(66); A″ type 887(229), 459(35), 367(23). Ivanovo State Chemical Technological Academy. Translated fromZhurnal Strukturnoi Khimii, Vol. 37, No. 3, pp. 440–449, May–June, 1996.     

Ab initio based polarizable force field generation and application to liquid silica and magnesia

We extend the program potfit, which generates effective atomic interaction potentials from ab initio data, to electrostatic interactions and induced dipoles. The potential parametrization algorithm uses the Wolf direct, pairwise summation method with spherical truncation. The polarizability of oxygen atoms is modeled with the Tangney-Scandolo interatomic force field approach. Due to the Wolf summation, the computational effort in simulation scales linearly in the number of particles, despite the presence of electrostatic interactions. Thus, this model allows to perform large-scale molecular dynamics simulations of metal oxides with realistic potentials. Details of the implementation are given, and the generation of potentials for SiO(2) and MgO is demonstrated. The approach is validated by simulations of microstructural, thermodynamic, and vibrational properties of liquid silica and magnesia.     

Ab initio spectroscopic study for the NaRb molecule in ground and excited states

A wide adiabatic study is performed for NaRb molecule, involving 15¹Σ⁺ electronic states including the ionic state Na^?Rb⁺, as well as 14³Σ⁺, 1–9^1,3Π, and 1–5^1,3Δ states. This investigation is performed using an ab initio approach which involves the effective core potential, the core polarization potential with l‐dependent cut‐off functions. The NaRb system has been treated as a two‐electron system and the full valence configuration interaction is easily achieved. The spectroscopic constants R_e, D_e, T_e, ω_e, ω_ex_e, B_e, and D₀ for all these states are derived. We have also computed the vibrational levels as well their spacing for different values of J. In addition, permanent and transition dipole moments are determined and analyzed. The Dunham coefficients have been used to perform experimental spacing to compare directly with our results. The present calculations on NaRb extend previous theoretical works to numerous electronic excited states in the various symmetries. © 2014 Wiley Periodicals, Inc.