Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio calculations”

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.

At least 199 records · Page 11

Transport properties in the Jovian atmosphere

Transport properties in a Jupiter-like atmosphere (89 mol % hydrogen and 11 mol % helium) are obtained by using the method of the kinetic theory of gases. The transport collision integrals are calculated by fitting various two-body semiempirical interaction potentials for which the collision integrals are tabulated to ab initio quantum mechanical calculations of the two-body interactions. The collision integrals are used to calculate the binary diffusion coefficients, viscosity, and 'total' thermal conductivity of the pure gases and the gas mixtures at 1-atm pressure from 1000 K to 25,000 K.

Biolsi, L.↗

An ab initio investigation of the structure, vibrational frequencies, and intensities of HO2 and HOCl

The infrared spectral intensities for HOCl and HO2 have been calculated using a new ab initio technique. Theoretical results for the geometries, vibrational frequencies, and the dipole moments of these species are also reported. All of the calculations were performed at the SCF level using near Hartree-Fock quality basis sets. The results for the molecular geometries and the vibrational frequencies are in good agreement with available experimental data. It is believed that the computed intensities are accurate to at least 50%. The results should be helpful in attempts to determine the stratospheric abundance of HOCl and HO2 by in situ infrared spectroscopic measurements.

Komornicki, A.↗

Theoretical and Experimental Determination of the Proton Affinity of (CF3CH2)2O

We report the experimental determination of the proton affinity of the molecule (CF3CH2)2O using chemical ionization mass spectrometry, and we compare it to the theoretical value obtained for protonation at the oxygen atom using the calculational methodology (MP2/6-31G**//MP2/3-21G). The proton affinity for this molecule as measured by bracketing experiments was between 724 kJ/mole and 741 kJ/mole. Ab initio (MP2/6-31G**//MP2/3-21G) calculations yield a value of about 729 kJ/mole, in agreement with the chemical ionization experiments. The results of these and related calculations suggest that the (MP2/6-31G**//MP2/3-21G) methodology is acceptable for estimating the proton affinities of partially-and fully-fluorinated methyl and ethyl ethers. We submit that any conclusions about the chemistry of fluoroether polymer lubricants based on their basicity can also be predicted reliably with such calculations.

Zehe, Michael J.↗

Effective Potential Energies and Transport Cross Sections for Atom-Molecule Interactions of Nitrogen and Nitrogen

The potential energy surfaces for H2-N and N2-N interactions are calculated by accurate ab initio methods and applied to determine transport data. The results confirm that an effective potential energy for accurately determining transport properties can be calculated using a single orientation. A simple method is developed to determine the dispersion coefficients of effective potential energies Effective potential energies required for O2-O collisions are determ=ined. The H2-N, N2-N, O2-H, and O2-O collision integrals are calculated and tabulated for a large range of temperatures. The theoretical values of the N2-N and O2-O diffusion coefficients compare well with measured data available at room temperature.

Stallcop, James R.↗

MCSCF potential energy surface for photodissociation of formaldehyde

The ground state potential energy surface for the dissociation of formaldehyde (H2CO to H2 and CO) is calculated with the ab initio MCSCF method with an extended (4-31G) basis set. The location, barrier height, and force constants of the transition state are determined, and the normal coordinate analysis is carried out. The calculated barrier height is 4.5 eV. Based on the calculated quantities, the detailed mechanism of the photochemical dissociation is discussed.

Jaffe, R. L.↗

A global ab initio potential for HCN/HNC, exact vibrational energies, and comparison to experiment

An ab initio (i.e., from first principles) calculation of vibrational energies of HCN and HNC is reported. The vibrational calculations were done with a new potential derived from a fit to 1124 ab initio electronic energies which were calculated using the highly accurate CCSD(T) coupled-cluster method in conjunction with a large atomic natural orbital basis set. The properties of this potential are presented, and the vibrational calculations are compared to experiment for 54 vibrational transitions, 39 of which are for zero total angular momentum, J = 0, and 15 of which are for J = 1. The level of agreement with experiment is unprecedented for a triatomic with two nonhydrogen atoms, and demonstrates the capability of the latest computational methods to give reliable predictions on a strongly bound triatomic molecule at very high levels of vibrational excitation.

Bentley, Joseph A.↗

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗

How much hydrogen is there in a white dwarf?

Stratified hydrogen/helium envelope models in diffusive equilibrium are calculated for a 0.6-solar-mass white dwarf for effective temperatures between 10,000 and 80,000 K in order to investigate the observational constraints placed on the total hydrogen mass. Convective mixing is included ab initio in the calculations, and synthetic spectra are used for comparing these models with observational materials. It is shown that evolutionary changes in the surface composition of white dwarfs cannot be explained by a model in which a small amount of hydrogen floats to the surface from initially being mixed in the outer parts of a helium envelope. It is pointed out that the shape of the hydrogen lines can be used for constraining theories of convective overshoot.

Macdonald, James↗

Ab initio projected-unrestricted Hartree-Fock calculation of some potential energy curves for carbonyl fluoride

Some potential energy curves for CF2O were calculated using projected-unrestricted Hartree-Fock (PUHF) theory. The calculations employed a contracted (4s 3p) Gaussian-type atomic orbital basis set. Bound states were found for the X-tilde 1A1 and 1,3A2 states while the 1,3B1 and 1,3B2 states were repulsive in the valence representation. The merits of the PUHF treatment for excited states are discussed. The results are discussed in terms of available experimental information and previous calculations with particular emphasis on the question of the photolysis channels open in the solar spectral region.

Brewer, D. A.↗

The electronic Hamiltonian for cuprates

A realistic many-body Hamiltonian for the cuprate superconductors should include both copper d and oxygen p states, hopping matrix elements between them, and Coulomb energies, both on-site and inter-site. We have developed a novel computational scheme for deriving the relevant parameters ab initio from a constrained occupation local density functional. The scheme includes numerical calculation of appropriate Wannier functions for the copper and oxygen states. Explicit parameter values are given for La2CuO4. These parameters are generally consistent with other estimates and with the observed superexchange energy. Secondly, we address whether this complicated multi-band Hamiltonian can be reduced to a simpler one with fewer basis states per unit cell. We propose a mapping onto a new two-band effective Hamiltonian with one copper d and one oxygen p derived state per unit cell. This mapping takes into account the large oxygen-oxygen hopping given by the ab initio calculations.

Annett, James F.↗

Bound-free Spectra for Diatomic Molecules

It is now recognized that prediction of radiative heating of entering space craft requires explicit treatment of the radiation field from the infrared (IR) to the vacuum ultra violet (VUV). While at low temperatures and longer wavelengths, molecular radiation is well described by bound-bound transitions, in the short wavelength, high temperature regime, bound-free transitions can play an important role. In this work we describe first principles calculations we have carried out for bound-bound and bound-free transitions in N2, O2, C2, CO, CN, NO, and N2+. Compared to bound ]bound transitions, bound-free transitions have several particularities that make them different to deal with. These include more complicated line shapes and a dependence of emission intensity on both bound state diatomic and atomic concentrations. These will be discussed in detail below. The general procedure we used was the same for all species. The first step is to generate potential energy curves, transition moments, and coupling matrix elements by carrying out ab initio electronic structure calculations. These calculations are expensive, and thus approximations need to be made in order to make the calculations tractable. The only practical method we have to carry out these calculations is the internally contracted multi-reference configuration interaction (icMRCI) method as implemented in the program suite Molpro. This is a widely used method for these kinds of calculations, and is capable of generating very accurate results. With this method, we must first of choose which electrons to correlate, the one-electron basis to use, and then how to generate the molecular orbitals.

Schwenke, David W.↗

Nucleosynthesis in AGB stars: Observation of Mg-25 and Mg-26 in IRC+10216 and possible detection of Al-26

We report the detection in the circumstellar envelope IRC+10216 of millimeter lines of the rare isotopomers (25)MgNC and (26)MgNC, as well as of a line at 234433 MHz, which could be the J= 7-6 transition of (26)AlF (an alternate, although less likely identified would be the J= 9-8 transition of NaF). The derived Mg-24:Mg-25:Mg-26 isotopic abundance ratios (78 : 11+/- 1 : 11 +/-1) are consistent with the solar system values (79.0:10.0:11.0), following Anders & Grevesse 1989). According to new calculations of evolutionary models of 3 solar mass and 5 solar mass asymptotic giant branch (AGB) stars, these ratios and the previously measured N, O and Si isotopic ratios imply that the central star had an initial mass 3 solar mass (less than or equal to M(sub *, ini) less than 5 solar mass and has already experienced many 3rd dredge-up events. From this, it can be predicted that the Al-26/Al-27 isotopics ratio lies between 0.01 and 0.08; in fact, the value derived in the case that U234433 arises from (26)AlF is Al-26/Al-27 = 0.04. The identification of the (25)MgNC and (26)MgNC lines was made possible by ab-initio quantum mechanical calculations of the molecule geometrical structure. It was confirmed through millimeter-wave laboratory measurements. The quantum mechanical calculations are briefly described and the laboratory results presented in some detail. The rotation constants B, D, H and the spin-rotation constant gamma of (25)MgNC and (26)MgNC are determined from a fit of laboratory and astronomical data.

Guelin, M.↗

Anharmonic Rovibrational Calculations of Singlet Cyclic C4 Using a New Ab Initio Potential and a Quartic Force

We report a CCSD(T)/cc-pCV5Z quartic force field (QFF) and a semi-global CCSD(T)-F12b/aug-cc-pVTZ potential energy surface (PES) for singlet, cyclic C4. Vibrational fundamentals, combinations and overtones are obtained using vibrational second-order perturbation theory (VPT2) and the vibrational configurationinteraction (VCI) approach. Agreement is within 10 cm(exp −1) between the VCI calculated fundamentals on the QFF and PES using the MULTIMODE (MM) program, and VPT2 and VCI results agree for the fundamentals. The agreement between VPT2- QFF and MM-QFF results is also good for the C4 combinations and overtones. The J = 1 and J = 2 rovibrational energies are reported from both VCI (MM) on the PES and VPT2 on the QFF calculations. The spectroscopic constants of (12)C4 and two C(sub 2v)-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

quartic force field↗

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum mechanical theories of collision induced electronic transitions. Within the two state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases are considered: the avoided crossing between two potential surfaces; and the spin orbit interaction due to a P-2 halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum-mechanical theories of collision-induced electronic transitions. Within the two-state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases have been considered: (1) the avoided crossing between two potential surfaces, and (2) the spin-orbit interaction due to a 2P halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗

Electron-impact excitation of the low-lying electronic states of HCN

The first study of the low-energy electron-impact excitation of low-lying electronic transitions in the HCN molecule is reported. Measurements were made at incident electron energies of 11.6 and 21.6 eV in the energy-loss range of 3-10 eV, and at scattering angles of 20-130 deg. Inelastic scattering spectra were placed on the absolute cross-section scale by determining first the ratio of inelastic-to-elastic scattering cross sections, and then separately measuring the absolute elastic scattering cross section. Several new electronic transitions are observed which are intrinsically overlapped in the molecule itself. Assignments of these electronic transitions are suggested. These assignments are based on present spectroscopic and cross-sections measurements, high-energy electron scattering spectra, optical absorption spectra, and ab initio molecular orbital calculations.

Chutjian, A.↗

Prediction of molecular infrared band and line intensities

The general problem of predicting infrared intensities is examined. The problem of predicting line intensities is discussed with attention focussed specifically on the problem of predicting vibrational strengths (or total band intensities). Two methods are considered: predictions based on ab initio quantum mechanical calculations, and predictions based on the transfer of atomic polar tensor parameters from small model molecules to the molecule of interest. Comparisons of predicted spectra with experimental spectra are shown for CH3OH, CH3OCH3, C2H5OH, and HOC2. It is suggested that the most important application of this work may be to the prediction of intensities for free radicals, ions, and reactive intermediates so that the infrared spectrum can be used in the quantitative analysis for these interesting species.

Person, W. B.↗