Exploring chromium (VI) dioxodihalides chemistry: is density functional theory the most suitable tool?

A comparative systematic study of the CrO2F2 compound has been performed using different conventional ab initio methodologies and density functional procedures. Two points have been analyzed: first, the accuracy of results yielded by each method under study, and second, the computational cost required to reach such results. Weighing up both aspects, density functional theory has been found to be more appropriate than the Hartree-Fock (HF) and the analyzed post-HF methods. Hence, the structural characterization and spectroscopic elucidation of the full CrO2X2 series (X=F,Cl,Br,I) has been done at this level of theory. Emphasis has been given to the unknown CrO2I2 species, and specially to the UV/visible spectra of all four compounds. Furthermore, a topological analysis in terms of charge density distributions has revealed why the valence shell electron pair repulsion model fails in predicting the molecular shape of such CrO2X2 complexes

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

Algorithm for generating defective graphene sheets

An algorithm is presented for the generation of molecular models of defective graphene fragments, containing a majority of 6-membered rings with a small number of 5- and 7-membered rings as defects. The structures are generated from an initial random array of points in 2D space, which are then subject to Delaunay triangulation. The dual of the triangulation forms a Voronoi tessellation of polygons with a range of ring sizes. An iterative cycle of refinement, involving deletion and addition of points followed by further triangulation, is performed until the user-defined criteria for the number of defects are met. The array of points and connectivities are then converted to a molecular structure and subject to geometry optimization using a standard molecular modeling package to generate final atomic coordinates. On the basis of molecular mechanics with minimization, this automated method can generate structures, which conform to user-supplied criteria and avoid the potential bias associated with the manual building of structures. One application of the algorithm is the generation of structures for the evaluation of the reactivity of different defect sites. Ab initio electronic structure calculations on a representative structure indicate preferential fluorination close to 5-ring defects.

Harmonic force fields from scaled SCF calculations: Program ASYM40

We report an extended version of our normal coordinate program ASYM40, which may be used to transform Cartesian force constants from ab initio calculations to a force field in nonredundant internal (symmetry) coordinates. When experimental data are available, scale factors for the theoretical force field may then be optimized by least-squares refinement. The alternative of refining an empirical force field to fit a wide variety of data, as with the previous version ASYM20, has been retained. We compare the results of least-squares refinement of the full harmonic force field with least-squares refinement of only the scale factors for an SCF calculated force field and conclude that the latter approach may be useful for large molecules where more sophisticated calculations are impractical. The refinement of scale factors for a theoretical force field is also useful when there are only limited spectroscopic data. The program will accept ab initio calculated force fields from any program that presents Cartesian force constants as output. The program is available through Quantum Chemistry Program Exchange.

Analytical potentials for triatomic molecules from spectroscopic data V. Application to HOX (X=F, Cl, Br, I)

Analytical potential energy functions are reported for HOX (X=F, Cl, Br, I). The surface for HOF predicts two metastable minima as well as the equilibrium configuration. These correspond to HFO (bent) and OHF (linear). Ab initio calculations performed for the HOF surface confirm these predictions. Comparisons are drawn between the two sets of results, and a vibrational analysis is undertaken for the hydrogen bonded OHF species. For HOCl, one further minimum is predicted, corresponding to HClO (bent), the parameters for which compare favourably with those reported from ab initio studies. In contrast, only the equilibrium configurations are predicted to be stable for HOBr and HOI.

A potential energy surface for the ground state of acetylene, C2H2

A potential energy function has been derived for the ground state surface of C2H2 as a many-body expansion. The 2- and 3-body terms have been obtained by preliminary investigation of the ground state surfaces of CH2( 3B1) and C2H( 2Σ+). A 4-body term has been derived which reproduces the energy, geometry and harmonic force field of C2H2. The potential has a secondary minimum corresponding to the vinylidene structure and the geometry and energy of this are in close agreement with predictions from ab initio calculations. The saddle point for the HCCH-H2CC rearrangement is predicted to lie 2•530 eV above the acetylene minimum.

Nitrous acid: vibrational frequency shifts due to 15N substitution, and the harmonic force field

Vibration rotation spectra of HO15 NO and DO15 NO have been measured at a resolution of 0•04 cm-1 to determine the isotopic shifts in the vibrational band origins. These have been used together with recently determined data on the vibrational band origins, Coriolis constants, and centrifugal distorition constants, to determine the harmonic force field of both cis and trans nitrous acid in least squares refinement calculations. The results are discussed in relation to recent ab initio calculations, the inertia defects, and the torsional potential function.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

Analytical functions for the ground state potential surfaces of C3 (X 1 Σ g+) and HCN (X 1 Σ +)

Analytic functions have been obtained to represent the potential energy surfaces of C3 and HCN in their ground electronic states. These functions closely reproduce the available data on the energy, geometry, and force constants in all stable conformations, as well as data on the various dissociation products, and ab initio calculations of the energy at other conformations. The form of the resulting surfaces are portrayed in various ways and discussed briefly.

Calculated rotational and rovibrational spectra of D2S and HDS

Rovibrational energy levels, transition frequencies, and linestrengths are computed variationally for the sulfur hydrides D2S and HDS, using ab initio potential energy and dipole surfaces. Wave-numbers for the pure rotational transitions agree to within 0.2 cm−1 of the experimental lines. For the fundamental vibrational transitions, the band origins for D2S are 860.4, 1900.6, and 1912.0 cm−1 for ν2, ν1, and ν3, respectively, compared with the corresponding experimental values of 855.4, 1896.4, and 1910.2 cm−1. For HDS, we compute ν2 to be 1039.4 cm−1, compared with the experimental value of 1032.7 cm−1. The relative merits of local and normal mode descriptions for the overtone stretching band origins are discussed. Our results confirm the local mode nature of the H2S, D2S, and HDS system.

Analytical potentials for triatomic molecules VIII. A two valued surface for the lowest 1A¢ surfaces of H2O

Ab initio calculations of the energy have been made at approximately 150 points on the two lowest singlet A' potential energy surfaces of the water molecule, 1A' and 1A', covering structures having D∞h, C∞v, C2v and Cs symmetries. The object was to obtain an ab initio surface of uniform accuracy over the whole three-dimensional coordinate space. Molecular orbitals were constructed from a double zeta plus Rydberg basis, and correlation was introduced by single and double excitations from multiconfiguration states which gave the correct dissociation behaviour. A two-valued analytical potential function has been constructed to fit these ab initio energy calculations. The adiabatic energies are given in our analytical function as the eigenvalues of a 2 2 matrix, whose diagonal elements define two diabatic surfaces. The off-diagonal element goes to zero for those configurations corresponding to surface intersections, so that our adiabatic surface exhibits the correct Σ/II conical intersections for linear configurations, and singlet/triplet intersections of the O + H2 dissociation fragments. The agreement between our analytical surface and experiment has been improved by using empirical diatomic potential curves in place of those derived from ab initio calculations.

The equilibrium structure of monofluoroacetylene

The equilibrium rotational constants Be of HCCF and DCCF have been determined from the ground state rotational constants B0, by determining the αr constants for all five fundamentals from the high-resolution vibrational—rotation spectrum making appropriate corrections for the effects of Fermi resonance. By combination with results from the 13C isotopomers and the recent ab initio calculations by Botschwina (Chem. Phys. Lett., 209 (1993) 117), the equilibrium structure is deduced to be: re(CH) = 1.0555(15) Å, re(CC) = 1.1955(8) Å and re(CF) = 1.2781(8) Å.

The general harmonic force field of methyl fluoride

The complete general harmonic force field of methyl flouride was recalculated using the most recent literature frequency, Coriolis ζ, and centrifugal distortion data for 12CH3F, 13CH3F, 12CD3F, 12CHD2F and 12CH2DF. The anharmonic corrections applied to the observed frequency data and the adopted molecular geometry are considered to be more realistic than those used hitherto. There is excellent overall agreement between the fitted force constants and the highest quality ab initio force field currently available.

Potential energy functions for ground state surfaces of HCO and HNO

Analytical potential functions are reported for the ground state surfaces of HCO and HNO, the functions being derived from spectroscopic and ab initio data. Harmonized force fields have been deduced for the stable configurations of both molecules and vibration frequencies predicted for the metastable species COH and NOH.