SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Ionization, excitation and total energies of oxygen and nitrogen, using self-consistent field molecular orbital method
Using the optimized valence configurations (OVC) multiconfiguration self-consistent-field (MCSCF) method, the dipole moment function for the ground state of CO in the vicinity of the equilibrium internuclear distance has been calculated. The OVC MCSCF calculation results are compared with existing Hartree-Fock and configuration interaction treatments of this molecule at single points and also the dipole moment function deduced from experimental infrared intensities. A general prescription for constructing OVC wavefunctions for diatomic molecules is also presented.
Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.
Thioperoxy (XSO or XOS) radicals exist in a variety of chemical environments, and they have as a consequence drawn some interest. HSO, an important species in the chemistry of the troposphere, has been examined both experimentally. The halogenated (X = F, Cl or Br) peroxy species and isovalent thioperoxy species have been studied less, but they too are potentially interesting because oxidized sulfur species and halogen sources are present in the atmosphere. Learning the fate of XSO and XOS radicals is important to understanding the atmospheric oxidation chemistry of sulfur compounds. Of these, FSO and ClSO are particularly interesting because they have been directly detected spectroscopically. Recent studies in our laboratory on the photochemistry of thionyl halides (X2SO; where X = F or Cl) have suggested new ways to generate XSO species. The laser-induced photodissociation of thionyl fluoride, F2SO, at 193 nm and thionyl chloride, ClSO, at 248 nm is characterized by a radical mechanism, X2SO -> XSO + X. The structure of FSO has been characterized experimentally by Endo et cd. employing microwave spectroscopy. Using the unrestricted Hartree-Fock (UHF) self-consistent field (SCF) method, Sakai and Morokuma computed the electronic structure of the ground (sup 2)A" and the first excited (sup 2)A' states of FSO. Electron correlation was not taken into account in their study. In a laser photodissociation experiment, Huber et al. identified ClSO mass spectromctrically. ClSO has also been detected in low temperature matrices by EPR and in the gas phase by far IR laser magnetic resonance. Although the structure of FSO is known in detail, the only study, experimental or theoretical, of CISO has been an ab initio HFSCF study by Hinchliffe. Electron correlation corrections were also excluded from this study. In order to better understand the isomerization and dissociation dynamics of the radical species, we have performed ab initio correlated studies of the potential energy surfaces (PES) of ClSO and its isomer ClOS at the QCISD(T)/6-31 G* level of theory. For FSO and FOS, more extensive QCISD/6-31 1G(2df) calculations have been possible, and the results are summarized here.
Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman intensities, are calculated for MgSO4.nH20 (n=1-3). Electronic structure theory at the second order Moller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set is used to determine the geometry, properties, and vibrational spectra of pure and hydrated MgSO4 salts. The direct vibrational self-consistent field (VSCF) method and its correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to vibrational frequencies and intensities for the pure MgSO4 and its complex with one water molecule. Very significant differences are found between vibrational of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies are shifted to the red very significantly (by up to 1500-2000/cm) upon complexation with magnesium sulfate. They should be observed between 1700 and 3000/cm in a region very different from the corresponding O-H stretch frequency region of pure water (3700-3800/cm). In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. They can serve as unique identifiers for the presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable help in the design of future missions and analysis of observed data from the ice surface of Jupiter's moon Europa that possibly contains hydrated MgSO4 salts.
We report a theoretical study of the dissociative ionization (DI) and product distributions of benzene (C6H6) and pyridine (C5H5N) from their low-lying ionization channels. Our approach makes use of the fact that electronic motion is much faster than nuclear motion allowing DI to be treated as a two-step process. The first step is the electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step, the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model [W.M. Huo, Phys. Rev. A64,042719-I (2001)]. For the unimolecular dissociation, we use multiconfigurational self-consistent field (MCSCF) methods to determine the steepest descent pathways to the possible product channels. More accurate methods are then used to obtain better energetics of the paths which are used to determine unimolecular dissociation probabilities and product distributions. Our analysis of the dissociation products and the thresholds of their productions for benzene are compared with the recent dissociative photoionization meausurements of benzene by Feng et al. [R. Feng, G. Cooper, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,211 (2002)] and the dissociative photoionization measurements of pyridine by Tixier et al. [S. Tixier, G. Cooper, R. Feng, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,185 (2002)] using dipole (e,e+ion) coincidence spectroscopy.
Self consistent field molecular orbital method in LCAO /linear combination of atomic orbitals/ approximation applied to LiH ground state for potential energy curve
HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.
HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.
The potential energy surface (PES) for NH2+NO is characterized using complete active space self-consistent-field (CASSCF)/derivative methods to locate the stationary points, followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. Production of N2+H2O is found to involve a complex mechanism, which, however, has no barrier with respect to NH2+NO. This pathway is exothermic by 124.5 kcal/mol. Production of NH2+OH can occur with no barrier other than the exothermicity from any of three isomers of HNNOH.
The reaction of vinylidene (CH2C) with acetylene may be an initiating reaction in soot formation. We report minimum energy paths and accurate energetics for a pathway leading to vinyl-acetylene and for a number of isomers of C4H4. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contacted configuration interaction (ICCI) and/or coupled cluster singles and doubles with a perturbational estimate of triple excitations (CCSD(T)) to determine the energetics. We find an entrance channel barrier of about 5 kcal/mol for the addition of vinylidene to acetylene, but no barriers above reactants for the reaction pathway leading to vinyl-acetylene.
The reaction of vinylidene (CH2C) with acetylene may be an initiating reaction in soot formation. We report minimum energy paths and accurate energetics for a pathway leading to vinylacetylene and for a number of isomers Of C4H4. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contacted configuration interaction (ICCI) and/or coupled cluster singles and doubles with a perturbational estimate of triple excitations (CCSD(T)) to determine the energetics. We find an entrance channel barrier of about 5 kcal/mol for the addition of vinylidene to acetylene, but no barriers above reactants for the reaction pathway leading to vinylacetylene.
We report calculations of the minimum energy pathways connecting CH2 + N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use Complete Active Space Self-Consistent Field (CASSCF) derivative methods to characterize the stationary points and Internally Contracted Configuration Interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.
We report calculations of the minimum energy pathways connecting (1)CH2+N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contracted configuration interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction of (1)CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 Surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.
A minimum energy pathway for interchange of the CH groups in ketene via a C2v structure has been obtained using complete active space self consistent field (CASSCF) derivative methods with a polarized valence double zeta basis set to define the reaction pathway followed by multi-reference internally contracted configuration interaction (ICCI) calculations with a [3s3p2d/3s2p] basis set to determine the energetics. Qualitatively, the C2v structure is found to be a shallow minimum on the potential energy surface separated from ketene by a small barrier (0.2 kcal/mol), a second minimum, and a larger barrier (3.0 kcal/mol). The minimum energy pathway leading from the C2v minimum to ketene starts by simultaneous rotation of the farther CH group out of the plane and away from the oxygen followed by increase of the CCO angle and subsequent 1,2-migration of the H of the nearer CH group toward the carbon of the farther CH group.
Eleven stationary geometries of quartet O4(+) have been studied by ab initio methods. The geometries were optimized at the Complete Active Space Self-Consistent Field (CASSCF) level of theory and the energies were calculated by the multiconfigurational second order pertubation method (CASPT2), using Double-Zeta Plus polarization (DZP), Triple-Zeta Plus Double Polarization (TZ2P), average Atomic Natural Orbital (ANO) (5s4p2d) and average ANO (6s5p3d2f) basis sets. The rectangular and trans-planar structures are found to be the most stable, with an energy barrier to conversion between the two at the threshold of dissociation. Both have a delocalized hole and are stable relative to separated 02 and 02(+)by 11.0 and 11.5 kcal/mol for the rectangular and the 2 trans-planar structure, respectively, compared with the experimentally deduced energy in the range of 9.2 to 10.8 kcal/mol. The adiabatic ionization potentials of 04 and 02 are computed to be 11.67 and 12.21 eV, while experimental values are 11.66 and 12.07 eV, respectively. The vibrational frequencies have been computed for all degrees of freedom at the CASSCF level of theory. Symmetry breaking is found to be a particular problem in the computation of the antisymmetric stretch frequency for the delocalized structures at the CASSCF level of theory. Attempts to rectify these problems using the restricted active space self-consistent field (RASSCF) method leads to additional difficulties, but further analysis yields insight into the symmetry breaking and problems with earlier calculations. Finally, a nonorthogonal configuration interaction (CI) calculation based on the interaction of localized CASSCF wave functions using the Complete Active Space State Interation (CASSI) method leads to a balanced treatment of the antisymmetric stretch which is free from symmetry breaking. The study explains the four most prominent absorption frequencies observed in the partially unassigned IR spectrum of O4(+) isolated in solid neon as the antisymmetric OO stretch, and the combination band of the symmetric and antisymmetric 00 stretch of both the rectangular and trans-planar structures.
The expectation of the solution process in a stochastic operator equation can be obtained from averaged equations only under very special circumstances. Conditions for validity are given and the significance and validity of the approximation in widely used hierarchy methods and the ?self-consistent field' approximation in nonequilibrium statistical mechanics are clarified. The error at any level of the hierarchy can be given and can be avoided by the use of the iterative method.
The second-order Moller-Plesset ab initio electronic structure method is used to compute points for the anharmonic mode-coupled potential energy surface of N-methylacetamide (NMA) in the trans(sub ct) configuration, including all degrees of freedom. The vibrational states and the spectroscopy are directly computed from this potential surface using the Correlation Corrected Vibrational Self-Consistent Field (CC-VSCF) method. The results are compared with CC-VSCF calculations using both the standard and improved empirical Amber-like force fields and available low temperature experimental matrix data. Analysis of our calculated spectroscopic results show that: (1) The excellent agreement between the ab initio CC-VSCF calculated frequencies and the experimental data suggest that the computed anharmonic potentials for N-methylacetamide are of a very high quality; (2) For most transitions, the vibrational frequencies obtained from the ab initio CC-VSCF method are superior to those obtained using the empirical CC-VSCF methods, when compared with experimental data. However, the improved empirical force field yields better agreement with the experimental frequencies as compared with a standard AMBER-type force field; (3) The empirical force field in particular overestimates anharmonic couplings for the amide-2 mode, the methyl asymmetric bending modes, the out-of-plane methyl bending modes, and the methyl distortions; (4) Disagreement between the ab initio and empirical anharmonic couplings is greater than the disagreement between the frequencies, and thus the anharmonic part of the empirical potential seems to be less accurate than the harmonic contribution;and (5) Both the empirical and ab initio CC-VSCF calculations predict a negligible anharmonic coupling between the amide-1 and other internal modes. The implication of this is that the intramolecular energy flow between the amide-1 and the other internal modes may be smaller than anticipated. These results may have important implications for the anharmonic force fields of peptides, for which N-methylacetamide is a model.