Analytical potentials for triatomic molecules IX. The prediction of anharmonic force constants from potential energy surfaces based on harmonic force fields and dissociation energies for SO2 and O3

Analytical potential energy functions which are valid at all dissociation limits have been derived for the ground states of SO2 and O3. The procedure involves minimizing the errors between the observed vibrational spectra and spectra calculated by a variational procedure. Good agreement is obtained between the observed and calculated spectra for both molecules. Comparisons are made between anharmonic force fields, previously determined from the spectral data, and the force fields obtained by differentiating the derived analytical functions at the equilibrium configurations.

Calorimetry of dehydrogenation and dangling-bond recombination in several hydrogenated amorphous silicon materials

Differential scanning calorimetry (DSC) was used to study the dehydrogenation processes that take place in three hydrogenated amorphous silicon materials: nanoparticles, polymorphous silicon, and conventional device-quality amorphous silicon. Comparison of DSC thermograms with evolved gas analysis (EGA) has led to the identification of four dehydrogenation processes arising from polymeric chains (A), SiH groups at the surfaces of internal voids (A'), SiH groups at interfaces (B), and in the bulk (C). All of them are slightly exothermic with enthalpies below 50 meV/H atoms , indicating that, after dissociation of any SiH group, most dangling bonds recombine. The kinetics of the three low-temperature processes [with DSC peak temperatures at around 320 (A),360 (A'), and 430°C (B)] exhibit a kinetic-compensation effect characterized by a linea relationship between the activation entropy and enthalpy, which constitutes their signature. Their Si-H bond-dissociation energies have been determined to be E (Si-H)0=3.14 (A), 3.19 (A'), and 3.28 eV (B). In these cases it was possible to extract the formation energy E(DB) of the dangling bonds that recombine after Si-H bond breaking [0.97 (A), 1.05 (A'), and 1.12 (B)]. It is concluded that E(DB) increases with the degree of confinement and that E(DB)>1.10 eV for the isolated dangling bond in the bulk. After Si-H dissociation and for the low-temperature processes, hydrogen is transported in molecular form and a low relaxation of the silicon network is promoted. This is in contrast to the high-temperature process for which the diffusion of H in atomic form induces a substantial lattice relaxation that, for the conventional amorphous sample, releases energy of around 600 meV per H atom. It is argued that the density of sites in the Si network for H trapping diminishes during atomic diffusion

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Dft theoretical study of energetic nitrogen-rich C4N6H8-n(NO2)n derivatives

Density functional theory (DFT) calculations at the B3LYP/6-31G** theoretical level were performed for a series of guanidine-fused bicyclic skeleton derivatives C4N6H8-n(NO2)n (n = 1 - 6). The heats of formation (HOFs) were calculated by isodesmic reactions, and the detonation properties were evaluated using the Kamlet - Jacobs equations. The bond dissociation energies were also analyzed to investigate the thermal stability and sensitivity of the compounds. The results show that all of the derivatives have high positive HOFs, compound G has the highest theoretical density, and compound F1 has the highest detonation velocity and detonation pressure. Considering both the detonation properties and thermal stabilities, compounds D1 and D4 (3 nitro substituents), E1 - E6 (4 nitro substituents), and G (6 nitro substituents) can be regarded as potential candidates for high-energy density materials.

An attempt on the calculation of relativistic potential energy curves: noble gas difluorides

Total energy SCF calculations were performed for noble gas difluorides in a relativistic procedure and compared with analogous non-relativistic calculations. The discrete variational method with numerical basis functions was used. Rather smooth potential energy curves could be obtained. The theoretical Kr - F and Xe - F bond distances were calculated to be 3.5 a.u. and 3.6 a.u. which should be compared with the experimental values of 3.54 a.u. and 3.7 a.u. Although the dissociation energies are off by a factor of about five it was found that ArF_2 may be a stable molecule. Theoretical ionization energies for the outer levels reproduce the experimental values for KrF_2 and XeF_2 to within 2 eV.

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Adsorption of organic molecules at the TiO2(110) surface: the effect of van der Waals interactions

Understanding the interaction of organic molecules with TiO2 surfaces is important for a wide range of technological applications. While density functional theory (DFT) calculations can provide valuable insight about these interactions, traditional DFT approaches with local exchange-correlation functionals suffer from a poor description of non-bonding van der Waals (vdW) interactions. We examine here the contribution of vdW forces to the interaction of small organic molecules (methane, methanol, formic acid and glycine) with the TiO2 (110) surface, based on DFT calculations with the optB88-vdW functional. The adsorption geometries and energies at different configurations were also obtained in the standard generalized gradient approximation (GGA-PBE) for comparison. We find that the optB88-vdW consistently gives shorter surface adsorbate-to-surface distances and slightly stronger interactions than PBE for the weak (physisorbed) modes of adsorption. In the case of strongly adsorbed (chemisorbed) molecules both functionals give similar results for the adsorption geometries, and also similar values of the relative energies between different chemisorption modes for each molecule. In particular both functionals predict that dissociative adsorption is more favourable than molecular adsorption for methanol, formic acid and glycine, in general agreement with experiment. The dissociation energies obtained from both functionals are also very similar, indicating that vdW interactions do not affect the thermodynamics of surface deprotonation. However, the optB88-vdW always predicts stronger adsorption than PBE. The comparison of the methanol adsorption energies with values obtained from a Redhead analysis of temperature programmed desorption data suggests that optB88-vdW significantly overestimates the adsorption strength, although we warn about the uncertainties involved in such comparisons.

Excited Electronic States, Transition Probabilities, and Radiative Lifetimes of CAs: A Theoretical Contribution and Challenge to Experimentalists

High-level CASSCF/MRCI calculations with a quintuple-zeta quality basis set are reported by characterizing for the first time a manifold of electronic states of the CAs radical yet to be investigated experimentally. Along with the potential energy curves and the associated spectroscopic constants, the dipole moment functions for selected electronic states as well as the transition dipole moment functions for the most relevant electronic transitions are also presented. Estimates of radiative transition probabilities and lifetimes complement this investigation, which also assesses the effect of spin-orbit interaction on the A (2)Pi state. Whenever pertinent, comparisons of similarities and differences with the isovalent CN and CP radicals are made.

Theoretical Characterization of Hydrogen Polyoxides: HOOH, HOOOH, HOOOOH, and HOOO

We have investigated the polyoxides HOOH, HOOOH, HOOOOH, and HOOO employing the CCSD(T) methodology, and the correlation consistent basis sets. For all molecules, we have computed fundamental vibrational frequencies, structural parameters, rotational constants, and rotation-vibration corrections. For HOOOH, we have obtained a good agreement between our results and microwave and infrared spectra measurements, although for the symmetric OO stretch some important differences were found. Heats of formation were computed using atomization energies, and our recommendation is as follows: Delta H degrees(f,298)(HOOOH) = -21.50 kcal/mol and Delta H degrees(f,298)(HOOOOH) = -10.61 kcal/mol. In the case of HOOO, to estimate the heat of formation, we have constructed three isodesmic reactions to cancel high order correlation effects. The results obtained confirmed that the latter effects are very important for HOOO. The new Delta H degrees(f,298)(HOOO) obtained is 5.5 kcal/mol. We have also calculated the zero-point energies of DO and DOOO to correct the experimental lower limit determined for the Delta H degrees(f,298)(HOOO). The Delta(Delta ZPE) decreases the binding energy of HOOO by 0.56 kcal/mol. Employing the latter value, the new experimental lower limit for Delta H degrees(f,298)(HOOO) is 3.07 kcal/mol, just 2.4 kcal/mol lower than our determination. We expect that the fundamental vibrational frequencies and rotational constants determined for HOOOOH and DOOOOD contribute to its identification in the gas phase. The vibrational spectrum of HOOOOH shows some overlapping with that of HOOOH thus indicating that one may encounter some difficulties in its characterization. We discuss the consequences of the thermochemical properties determined in this work, and suggest that the amount of HOOO present in the atmosphere is smaller than that proposed recently in this journal (J. Phys. Chem A 2007, 111, 4727).

Electronic structure and spectra of a new molecular species: SI. A theoretical contribution

A high level theoretical approach is used to characterize for the first time a manifold of doublet and quartet A + S and Omega states correlating with the first two dissociation channels of an as yet experimentally unknown molecular species, SI, sulfur monoidide. A set of spectroscopic constants is determined, including vibrationally averaged spin-orbit coupling constants, vibrationally averaged dipole moments, and dissociation energies. The transition dipole moment function for the spin-forbidden transition a (4)Sigma -X (2)Pi, and the associated radiative lifetimes were also evaluated. Two possibilities to detect transitions experimentally and to derive spectroscopic constants are suggested. (C) 2011 Elsevier B. V. All rights reserved.

Spin contamination in XOO radicals X = F, Cl, Br, HO: How is the investigation of the HOOO radical affected

We have analysed the effect of spin contamination in the wavefunction of HOOO. At least, two solutions can be found for the HF wavefunction. One, lower in energy, presents a high spin contamination and gives qualitatively incorrect structural parameters. On the other hand, the less contaminated HF reference gives structural parameters that are in better agreement with experiment, and positive spin densities on all atoms. Some of the problems described during previous investigations of HOOO can now be traced to problems in the HF reference. For the first time we report a CCSD(T) estimation of the structure of HOOO cis employing a HF reference with small spin contamination. (C) 2008 Elsevier B.V. All rights reserved.

Excited states of the CaAl molecule: An MRCI study

Accurate potential energy curves, dissociation energies and spectroscopic constants for several low-lying doublet and quartet electronic states of CaAl were investigated using the CASSCF/MRCI methodology, and the cc-pVQZ basis set. Our results represent an improvement over a previous theoretical description, and also characterizes new higher excited states not previously investigated, thus confirming the assignment of four excited states investigated experimentally. With the theoretical transition moment functions, transition probabilities and radiative lifetimes were estimated via Einstein spontaneous emission coefficients. (c) 2008 Elsevier B. V. All rights reserved.

A CASSCF/MRCI Study of the Low-Lying Electronic States of the BeS Molecule

Accurate potential energy curves, dipole moment functions, dissociation energies, and molecular constants for several low-lying singlet and triplet electronic states of BeS were investigated using the CASSCF/MRCI methodology, and the cc-pV5Z basis set for beryllium, and the aug-cc-pV(5+d) Z set for sulfur. Besides presenting improved results for the three lowest lying states, this study presents the first theoretical characterization of another set of nine excited states so far unknown experimentally. Our results are sufficiently accurate to reliably guide the experimental search and characterization of these states, and also to confirm the experimental assignment of the B(1)Sigma(+) excited state. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1694-1700, 2011

A theoretical study on the XeF(2) molecule

A relativistic four-component study was performed for the XeF(2) molecule by using the Dirac-Coulomb (DC) Hamiltonian and the relativistic adapted Gaussian basis sets (RAGBSs). The comparison of bond lengths obtained showed that relativistic effects on this property are small (increase of only 0.01 angstrom) while the contribution of electron correlation, obtained at CCSD(T) or CCSD-T levels, is more important (increase of 0.05 angstrom). Electron correlation is also dominant over relativistic effects for dissociation energies. Moreover, the correlation-relativity interaction is shown to be negligible for these properties. The electron affinity, the first ionization potential and the double ionization potential are obtained by means of the Fock-space coupled cluster (FSCC) method, resulting in DC-CCSD-T values of 0.3 eV, 12.5 eV and 32.3 eV, respectively. Vibrational frequencies and some anharmonicity constants were also calculated under the four-component formalism by means of standard perturbation equations. All these molecular properties are, in general, ill satisfactory agreement with available experimental results. Finally, a partition in terms of charge-charge flux-dipole flux (CCFDF) contributions derived by means of the quantum theory of atoms in molecules (QTAIM) in non-relativistic QCISD(FC)/3-21G* calculations was carried out for XeF(2) and KrF(2). This analysis showed that the most remarkable difference between both molecules lies on the charge flux contribution to the asymmetric stretching mode, which is negligible in KrF(2) but important in XeF(2). (c) 2008 Elsevier B.V. All rights reserved.

Cellulose pyrolysis and quantum chemistry

Cellulose is the major constituent of most plants of interest as renewable sources of energy and is the most extensively studied form of biomass or biomass constituent. Predicting the mass loss and product yields when cellulose is subjected to increased temperature represents a fundamental problem in the thermal release of biomass energy. Unfortunately, at this time, there is no internally consistent model of cellulose pyrolysis that can organize the varied experimental data now available or provide a guide for additional experiments. Here, we present a model of direct cellulose pyrolysis using a multistage decay scheme that we first presented in the IJQC in 1984. This decay scheme can, with the help of an inverse method of assigning reaction rates, provide a reasonable account of the direct fast pyrolysis yield measurements. The model is suggestive of dissociation states of d-glucose (C6H10O5,), the fundamental cellulose monomer. The model raises the question as to whether quantum chemistry could now provide the dissociation energies for the principal breakup modes of glucose into C-1, C-2, C-3, C-4, and C-5 compounds. These calculations would help in achieving a more fundamental description of volatile generation from cellulose pyrolysis and could serve as a guide for treating hemicellulose and lignin, the other major biomass constituents. Such advances could lead to the development of a predictive science of biomass pyrolysis that would facilitate the design of liquifiers and gasifiers based upon renewable feedstocks. (C) 1998 John Wiley & Sons, Inc.