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 289 records · Page 16

Computational Study of Electron-Molecule Collisions Related to Low-Temperature Plasmas

Computational study of electron-molecule collisions not only complements experimental measurements, but can also be used to investigate processes not readily accessible experimentally. A number of ab initio computational methods are available for these types of calculations. Here we describe a recently developed technique, the finite element Z-matrix method, Analogous to the R-matrix, method, it partitions the space into regions and employs real matrix elements. However, unlike the implementation of the R-matrix method commonly used in atomic and molecular physics, the Z-matrix method is fully variational. In the present implementation, a mixed basis of finite elements and Gaussians is used to represent the continuum electron, thus offering full flexibility without imposing fixed boundary conditions. Numerical examples include the electron-impact dissociation of N2 via the metastable A3Su+ state, a process which may be important in the lower thermosphere, and the dissociation of the CF radical, a process of interest to plasma etching. To understand the dissociation pathways, large scale quantum chemical calculations have been carried out for all target states which dissociate to the lowest five limits in the case of N2, and to the lowest two limits in the case of CF. For N2, the structural calculations clearly show the preference for predissociation if the initial state is the ground X1 Sg+ state, but direct dissociation appears to be preferable if the initial state is the A3Su+ state. Multi-configuration SCF target functions are used in the collisional calculation.

Huo, Winifred M.↗

Line List for the A 2 Π−X 2 Σ + and the B 2 Σ + −X 2 Σ + Band Systems of CaF

The spectral analysis of molecule-rich asymptotic giant branch (AGB) stars is challenging. Although other calcium and fluorine bearing molecules are observed in the microwave and in the visible spectrum of AGB stars, CaF has never been detected, despite favorable chemical equilibrium predictions. Yet, measuring the CaF abundance could give more insight into the fluorine budget of AGB stars and allow better simulation of stellar spectra. In this work, we present an analysis of the visible spectrum of CaF obtained with a Fourier transform spectrometer. CaF A 2 Π−X 2 Σ + and the B 2 Σ + −X 2 Σ + band systems were excited with a hollow cathode discharge. Using previous fluorescence spectroscopy measurements, highly accurate ground state constants, and reasonable extrapolation schemes, the strongest features of CaF in the visible spectrum can be accurately modeled for both band systems. Spectroscopic constants are determined for A 2 Π with v ≤ 16 and for B 2 Σ + with v ≤ 20. Ab initio transition dipole moment curves of both transitions were calculated and scaled, and we provide a line list with Einstein A coefficients and oscillator strengths. This line list can be used to simulate spectra of CaF at temperatures and pressures relevant to astrophysical environments.

Léo Lavy↗

Theory of collision-induced translation-rotation spectra: H2-He

An adiabatic quantal theory of spectral line shapes in collision-induced absorption and emission is presented which incorporates the induced translation-rotation and translation-vibration spectra. The generalization to account for the anisotropy of the scattering potential is given. Calculations are carried out of the collision-induced absorption spectra of He in collisions with H2 with ab initio electric dipole functions and realistic potentials. The anisotropy of the interaction potential is small and is not included in the calculations. The predicted spectra are in satisfactory agreement with experimental data though some deviations occur which may be significant. The rotational line shapes have exponential wings and are not Lorentzian. The connection between the quantal and classical theories is written out explicitly for the isotropic overlap induction.

Birnbaum, G.↗

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

Metal atom clusters are studied based on the application of ab initio quantum mechanical approaches. Because these large 'molecular' systems pose special practical computational problems in the application of the quantum mechanical methods, there is a special need to find simplifying techniques that do not compromise the reliability of the calculations. Research is therefore directed towards various aspects of the implementation of the effective core potential technique for the removal of the metal atom core electrons from the calculations.

Kahn, L. R.↗

The rate of the reaction between CN and C2H2 at interstellar temperatures

The rate coefficient for the important interstellar reaction between CN and C2H2 has been calculated as a function of temperature between 10 and 300 K. The potential surface for this reaction has been determined through ab initio quantum chemical techniques; the potential exhibits no barrier in the entrance channel but does show a small exit channel barrier, which lies below the energy of reactants. Phase-space calculations for the reaction dynamics, which take the exit channel barrier into account, show the same unusual temperature dependence as determined by experiment, in which the rate coefficient at first increases as the temperature is reduced below room temperature and then starts to decrease as the temperature drops below 50-100 K. The agreement between theory and experiment provides strong confirmation that the reaction occurs appreciably at cool interstellar temperatures.

NASA Discipline Exobiology↗

Curve crossing for low energy elastic scattering of He (plus) by Ne

The perturbation seen in the experimental differential elastic scattering cross section for the 40 eV He(+) + Ne system was attributed to a single crossing of two intermolecular potential energy curves. A new theoretical treatment of the curve crossing problem, that of Delos and Thorson, is employed to obtain the crossing probabilities and phases associated with the crossing. These are determined by utilizing ab initio potentials involved in the crossing and are further used in a partial wave calculation of the cross section, which is compared with our experiment. The origin of the oscillatory structure observed in the differential cross section is discussed in semiclassical terms by defining the problem in terms of two pseudo-deflection functions. A rainbow effect is shown to be related to a particular feature (a maximum rather than a minimum) of these deflection functions.

Bobbio, S. M.↗

A reinvestigation of the rate of the C/+/ + H2 radiative association reaction

It is noted that new experimental results and statistical theories have prompted a reinvestigation of the rate coefficient of the interstellar reaction in which C(+) + H2 yields CH2(+) + h(nu) in the 10-100 K temperature range. The results presented here indicate a rate coefficient between 10 to the -16 and 10 to the -15 cu cm/s at all temperatures studied. In applying the modified thermal and phase space theories, it is expected that they will be as accurate as for radiative association as for three-body association, provided the ab initio value for the radiative decay rate is correct. It is expected that the calculated values of the rate coefficient will be accurate to within an order of magnitude and will have the correct temperature dependence.

Herbst, E.↗

Dissociation and Internal Excitation of Molecular Nitrogen Due to N + N2 Collisions Using Direct Molecular Simulation

In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.

Grover, Maninder S.↗

Curve crossing for low-energy elastic scattering of He/+/ by Ne.

The perturbation seen in the experimental differential elastic-scattering cross section for the 40-eV He/+/ + Ne system has been attributed to a single crossing of two intermolecular potential-energy curves. A new theoretical treatment of the curve-crossing problem, namely, that of Delos and Thorson, is employed to obtain the crossing probabilities and phases associated with the crossing. These are determined by utilizing ab initio potentials involved in the crossing and are further used in a partial-wave calculation of the cross section, which is compared with our experiment. The origin of the oscillatory structure observed in the differential cross section is discussed in semiclassical terms by defining the problem in terms of two pseudo-deflection-functions. A rainbow effect is shown to be related to a particular feature (a maximum rather than a minimum) of these deflection functions.

Bobbio, S. M.↗

Accurate collision-induced line-coupling parameters for the fundamental band of CO in He - Close coupling and coupled states scattering calculations

The first accurate theoretical values for off-diagonal (i.e., line-coupling) pressure-broadening cross sections are presented. Calculations were done for CO perturbed by He at thermal collision energies using an accurate ab initio potential energy surface. Converged close coupling, i.e., numerically exact values, were obtained for coupling to the R(0) and R(2) lines. These were used to test the coupled states (CS) and infinite order sudden (IOS) approximate scattering methods. CS was found to be of quantitative accuracy (a few percent) and has been used to obtain coupling values for lines to R(10). IOS values are less accurate, but, owing to their simplicity, may nonetheless prove useful as has been recently demonstrated.

Green, Sheldon↗

A New Global Potential Energy Surface for the Hydroperoxyl Radical, HO2: Reaction Coefficients for H + O2 and Vibrational Splittings for H Atom Transfer

A new analytic global potential energy surface describing the hydroperoxyl radical system H((sup 2)S) + O2(X (sup 3)Sigma((sup -)(sub g))) (reversible reaction) HO2 ((X-tilde) (sup 2)A'') (reversible reaction) O((sup 3)P) + O H (X (sup 2)Pi) has been fitted using the ab initio complete active space SCF (self-consistent-field)/externally contracted configuration interaction (CASSCF/CCI) energy calculations of Walch and Duchovic. Results of quasiclassical trajectory studies to determine the rate coefficients of the forward and reverse reactions at combustion temperatures will be presented. In addition, vibrational energy levels were calculated using the quantum DVR-DGB (discrete variable representation-distributed Gaussian basis) method and the splitting due to H atom migration is investigated. The material of the proposed presentation was reviewed and the technical content will not reveal any information not already in the public domain and will not give any foreign industry or government a competitive advantage.

Dateo, Christopher E.↗

Theoretical calculation of low-lying states of NaAr and NaXe

Potential curves as well as dipole moments and linking transition moments are calculated for the ground X 2 Sigma + and low lying excited A 2 Pi, B 2 Sigma +, C 2 Sigma +, (4) 2 Sigma +, (2) 2 Pi and (1) 2 Delta states of NaAr and NaXe. Calculations are performed using a self-consistent field plus configuration-interaction procedure with the core electrons replaced by an ab initio effective core potential. The potential curves obtained are found to be considerably less repulsive than the semiempirical curves of Pascale and Vandeplanque (1974) and to agree well with existing experimental data, although the binding energies of those states having potential minima due to van der Waals interactions are underestimated. Emission bands are also calculated for the X 2 Sigma + - C 2 Sigma + excimer transitions of NaAr and NaXe using the calculated transition moments and potential curves, and shown to agree well with experiment on the short-wavelength side of the maximum.

Laskowski, B. C.↗

AB Initio Characterization of MgCCH, MgCCH(+), and MgC2, and Pathways to their Formation in the Interstellar Medium

A study of Mg-bearing compounds has been performed in order to determine molecular properties which are critical for planning new astronomical searches and laboratory studies. The primary focus of the work is on MgCCH, MgCCH(+), and the isomers of MgC2. Only MgCCH has been identified in laboratory studies. Additional calculations have been carried out on MgH, MgNC, MgCN, and their cations in an effort to evaluate pathways to the formation of MgCCH and MgCCH(+) in the InterStellar Medium (ISM) or in circumstellar envelopes. Correlated ab initio methods and correlation-consistent basis sets have been employed. Properties including structures, rotational constants, dipole moments, and harmonic frequencies are reported. A transition state between linear MgCC and cyclic MgC2 has been characterized and was found to yield a minimal barrier (approx. 0.5 kcal/mole), indicating easy interconversion to the cyclic form. Direct reactions in the ISM between Mg or Mg(+) and HCCH are precluded by energetic considerations, but a number of ion- molecule or neutral-neutral exchange reactions between CCH and various Mg-containing species offer plausible pathways to MgCCH or MgCCH(+). Weakly bound MgH may react with CCH to form MgCCH, but MgH has not been detected. Both MgNC and MgCN have been observed, but reactions with CCH are slightly endothermic by 1-3 kcal/mole. Although MgH(+), MgNC(+), and MgCN(+) have not been detected, their reactions with CCH to form MgCCH(+) are all exothermic. With only a small barrier separating linear MgCC and cyclic MgC2, the dissociative recombination of MgCCH(+) with an electron is expected to yield cyclic MgC2, and regenerate Mg and CCH. New astronomical searches for MgCCH, MgCCH(+), cyclic MgC2, MgNC(+), and MgCN(+) will provide further insight into organo-magnesium astrochemistry.

Woon, David E.↗

An Extended Ab Initio and Theoretical Thermodynamics Studies of the Bergman Reaction and the Energy Splitting of the Singlet Ortho-, Meta-, and Para-Benzynes

The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Lindh, Roland↗

Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes

Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.

Jaffe, Richard↗

Rotational and vibrational transitions for Li + H2 collisions

Close coupling calculations for integral and differential cross sections have been carried out for Li + H2 collisions with an ab initio Hartree-Fock potential energy surface. Rotational, vibrational, and vib-rotational excitation cross sections are reported at 0.4336 eV, 0.7 eV, and 0.8673 eV in the center of mass system. For pure rotational excitations, which dominate the inelastic scattering, coupling with vibrational states is not very important. For vibrational transitions, the influence of large multiquantum rotational transitions is far less than that found for Li(+) + H2 collisions.

Choi, B. H.↗

Predicted star counts in selected fields and photometric bands Applications to galactic structure, the disk luminosity function, and the detection of a massive halo

Tables are presented of predicted star counts in 17 selected directions on the sky as a function of apparent B, V, R, and I magnitudes and absolute visual magnitudes calculated from the Bahcall-Soneira Galaxy model. It is noted that eight of the fields were selected ab initio for their usefulness in determining from star counts either parameters of galactic structure, the faint end of the disk luminosity function, or the detection of a massive stellar halo. Star counts in specified fields to m sub I = 19 mag are able to reveal the faint end of the disk luminosity function down to the end of the hydrogen-burning main sequence. It is thought that star counts at high galactic latitudes to m sub I = 22 mag should reveal the stellar constituents of a massive halo if they are massive enough to be on the hydrogen-burning main sequence.

Bahcall, J. N.↗

Ab initio molecular orbital studies of low-energy, metastable isomers of the ubiquitous cyclopropenylidene

The discovery of cyclopropenylidene in space suggests that other C3H2 isomers may be present, and a tentative detection of one such isomer, propargylene (HCCCH), has been reported. Ab initio molecular orbital theory has been used to characterize five low-lying, metastable isomers of cyclopropenylidene. Extended calculations including the electron correlation energy, show that the lowest in energy is singlet propadienylidene, followed by propargylene; the singlet and triplet of the latter are too close in energy to allow an assignment of the ground state; triplet propadienylidene is at a significantly higher energy. Rotational frequencies computed to an expected accuracy of + or - 1 -2 percent do not confirm the tentative detection of propargylene in space, although the discrepancy between theory and the observation is not so great as to unequivocally rule out this possibility.

Defrees, D. J.↗