Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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.

162 records · Page 9

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.↗

A correlated ab initio study of the A2 pi <-- X2 sigma+ transition in MgCCH

The A2 pi <-- X2 sigma+ transition in MgCCH was studied with correlation consistent basis sets and single- and multireference correlation methods. The A2 pi excited state was characterized in detail; the x2 sigma+ ground state has been described elsewhere recently. The estimated complete basis set (CBS) limits for valence correlation, including zero-point energy corrections, are 22668, 23191, and 22795 for the RCCSD(T), MRCI, and MRCI + Q methods, respectively. A core-valence correction of +162 cm-1 shifts the RCCSD(T) value to 22830 cm-1, in good agreement with the experimental result of 22807 cm-1.

Non-NASA Center↗

Calculations of ion-molecule deuterium fractionation reactions involving HD

Gas-phase chemical models of deuterium fractionation in dense interstellar clouds utilize a small number of exothermic reactions to achieve fractionation. Although HD is a major repository of deuterium, it appears not to exchange deuterium with many molecular ions. Useful semiquantitative reasons have been given for the unusual lack of reactivity of exothermic ion-HD deuterium exchange systems, but quantum chemical studies are needed to understand these ideas in more detail and to determine if the lack of reactivity pertains at very low temperatures not studied in the laboratory, or whether tunneling can drive the reactions. Accordingly, the potential energy surfaces of three representative ion-molecule exchange reactions involving protonated ions and HD have been investigated with ab initio quantum chemical techniques. Our results generally confirm the semiquantitative picture as to which reactions are likely to occur and show that tunneling at low temperatures is unlikely to alter this picture.

Maluendes, Sergio A.↗

Charge-transfer energy in closed-shell ion-atom interactions

The importance of charge-transfer energy in the interactions between closed-shell ions and atoms is investigated. Ab initio calculations on H(plus)-He and Li(plus)-He are used as a guide for the construction of approximate methods for the estimation of the charge-transfer energy for more complicated systems. For many alkali ion-rate gas systems the charge-transfer energy is comparable to the induction energy in the region of the potential minimum, although for doubly charged alkaline-earth ions in rare gases the induction energy always dominates. Surprisingly, an empirical combination of repulsion energy plus asymptotic induction energy plus asymptotic dispersion energy seems to give a fair representation of the total interaction, especially if the repulsion energy is parameterized, despite the omission of any explicit charge-transfer contribution. More refined interaction models should consider the charge-transfer energy contribution.

Alvarez-Rizzatti, M.↗

A correlated ab initio study of linear carbon-chain radicals CnH (n = 2-7)

Linear carbon-chain radicals CnH for n = 2-7 have been studied with correlation consistent valence and core-valence basis sets and the coupled cluster method RCCSD(T). Equilibrium structures, rotational constants, and dipole moments are reported and compared with available experimental data. The ground state of the even-n series changes from 2 sigma+ to 2 pi as the chain is extended. For C4H, the 2 sigma+ state was found to lie only 72 cm-1 below the 2 pi state in the estimated complete basis set limit for valence correlation. The C2H- and C3H- anions have also been characterized.

NASA Discipline Exobiology↗

The C4H7+ cation. A theoretical investigation

The potential energy surface of the C4H7+ cation has been investigated with ab initio quantum chemical theory. Extended basis set calculations, including electronic correlation, show that cyclobutyl and cyclopropylcarbinyl cation are equally stable isomers. The saddle point connecting these isomers lies 0.6 kcal/mol above the minima. The global C4H7+ minimum corresponds to the 1-methylallyl cation, which is 9.0 kcal/mol more stable than the cyclobutyl and the cyclopropylcarbinyl cation and 9.5 kcal/mol below the 2-methylallyl cation. These results are in excellent agreement with experimental data.

NASA Discipline Number 52-10↗

adwTools Developed: New Bulk Alloy and Surface Analysis Software for the Alloy Design Workbench

A suite of atomistic modeling software, called the Alloy Design Workbench, has been developed by the Computational Materials Group at the NASA Glenn Research Center and the Ohio Aerospace Institute (OAI). The main goal of this software is to guide and augment experimental materials research and development efforts by creating powerful, yet intuitive, software that combines a graphical user interface with an operating code suitable for real-time atomistic simulations of multicomponent alloy systems. Targeted for experimentalists, the interface is straightforward and requires minimum knowledge of the underlying theory, allowing researchers to focus on the scientific aspects of the work. The centerpiece of the Alloy Design Workbench suite is the adwTools module, which concentrates on the atomistic analysis of surfaces and bulk alloys containing an arbitrary number of elements. An additional module, adwParams, handles ab initio input for the parameterization used in adwTools. Future modules planned for the suite include adwSeg, which will provide numerical predictions for segregation profiles to alloy surfaces and interfaces, and adwReport, which will serve as a window into the database, providing public access to the parameterization data and a repository where users can submit their own findings from the rest of the suite. The entire suite is designed to run on desktop-scale computers. The adwTools module incorporates a custom OAI/Glenn-developed Fortran code based on the BFS (Bozzolo- Ferrante-Smith) method for alloys, ref. 1). The heart of the suite, this code is used to calculate the energetics of different compositions and configurations of atoms.

Bozzolo, Guillermo↗

Solar flare gamma-ray line shapes

A computer code has been developed which is used to calculate ab initio the laboratory shapes and energy shifts of gamma-ray lines from (C-12)(p, gamma/4.438/)p-prime(C-12) and (O-16)(p, gamma/6.129/)p-prime(O-16) reactions and to calculate the expected shapes of these lines from solar flares. The sensitivity of observable solar flare gamma-ray line shapes to the directionality of the incident particles is investigated for several projectile angular distributions. Shapes of the carbon and oxygen lines are calculated assuming realistic proton energy spectra for particles in circular orbits at the mirror points of magnetic loops, for particle beams directed downward into the photosphere, and for isotropic particle distributions. Line shapes for flare sites near the center of the sun and on the limb are shown for both thin-target and thick-target interaction models.

Werntz, C.↗

Electron-nitrogen molecular collisions in high-temperature nonequilibrium air

Ab initio calculation of vibrational excitation and deexcitation cross sections of N2 by low energy resonant electron impact have been carried out. The calculation includes initial target state in v = 0-12 and final target state in v = 0-17. The calculated cross sections are found to agree with available experimental data to within reported experimental error. Effect of target rotation has been investigated. Comparison of vibrational excitation cross sections at J = 0, 50, and 150 shows that the major effect of target rotation is to lower the resonance energy. Vibrational excitation and deexcitation rate coefficients are tabulated for electron temperatures of 1100-55,000 K. The present tabulation supercedes an earlier compilation by the Naval Research Laboratory that was derived from experimental data with a factor-of-two error. The tabulated rate coefficients are suited for flowfield modeling of aeroassisted orbital transfer vehicles.

Huo, W. M.↗

A priori predictions of the rotational constants for HC13N, HC15N, C5O

Ab initio molecular orbital theory is used to estimate the rotational constant for several carbon-chain molecules that are candidates for discovery in interstellar space. These estimated rotational constants can be used in laboratory or astronomical searches for the molecules. The rotational constant for HC13N is estimated to be 0.1073 +/- 0.0002 GHz and its dipole moment 5.4 D. The rotational constant for HC15N is estimated to be 0.0724 GHz, with a somewhat larger uncertainty. The rotational constant of C5O is estimated to be 1.360 +/- 2% GHz and its dipole moment 4.4. D.

NASA Discipline Exobiology↗

Ab initio predictions on the rotational spectra of carbon-chain carbene molecules

We predict rotational constants for the carbon-chain molecules H2C=(C=)nC, n=3-8, using ab initio computations, observed values for the earlier members in the series, H2CCC and H2CCCC with n=1 and 2, and empirical geometry corrections derived from comparison of computation and experiment on related molecules. H2CCC and H2CCCC have already been observed by radioastronomy; higher members in the series, because of their large dipole moments, which we have calculated, are candidates for astronomical searches. Our predictions can guide searches and assist in both astronomical and laboratory detection.

NASA Discipline Number 52-10↗

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↗

Theoretical study of the dissociation energy and the red and violet band systems of CN

The dissociation energy (D sub O) of CN is determined to be 7.65 + or - 0.06 eV. This corresponds to delta H sub f (CN) = 105.3 + or - 1.5 kcal/mole, in excellent agreement with Engleman and Rouse (1975), but considerably larger than the recent value deduced from shock-tube studies by Colket (1984). The result is obtained not only from extensive ab initio MRCI calculations using a very large Gaussian basis set, but also from extrapolation of the directly computed value by comparison of computed and experimental results fo NO, C2, and N2. As an additional calibration of the methods, the D sub O value for CN was computed from the corresponding value for CN(-) using the experimental electron affinity data. The lifetime of the nu prime = 0 level of the violet (B 2 sigma + yields X 2 sigma +) system was computed to be 62.4 ns, in good agreement with both experiment and previous calculations. Lifetimes for the red (A 2 pi yields X 2 sigma +) system decrease with increasing nu prime, which is consistent both with the recent experiment and calculations. While the computed lifetimes are significantly longer that those obtained from the experiment, they are shorter than those deduced from an analysis of the solar spectrum. However the D sub O and f (sub OO) are consistent with Lambert's model for the solar spectrum.

Bauschlicher, Charles W., Jr.↗

A correlated ab initio study of the X2A1 and A2E states of MgCH3

The X2A1 and A2E states of the MgCH3 radical have been studied with correlation consistent basis sets and the coupled cluster method RCCSD(T) in order to compare with two recent experimental efforts [M. A. Anderson and L. M. Ziurys, Astrophys. J. 452, L157 (1995); R. Rubino, J. M. Williamson, and T. A. Miller, J. Chem. Phys. 103, 5964 (1995)]. The best computed values [RCCSD(T)/cc-pCVTZ] for the X2A1 state are (experimental results in parentheses): Ae = 160.433 GHz, Be = 10.948 GHz (B0 = 11.008 GHz), and Mue = 1.011 D. The Mg-CH3 bond is weak, 26.3 kcal/mol. Values for the A2E state are Ae = 154.648 GHz (A0 = 149.666 GHz), Be = 10.87 GHz (B0 = 10.932 GHz), and Mue = 1.022 D. The excitation energy (Te) for the A2E <-- X2A1 transition is 19 999 cm-1 (T00 = 20 030 cm-1). A brief discussion of bonding trends in Mg-containing radials is included.

Non-NASA Center↗

Theoretical IR spectra of ionized naphthalene

We report the results of a theoretical study of the effect of ionization on the IR spectrum of naphthalene, using ab initio molecular orbital theory. For that purpose we determined the structures, band frequencies, and intensities of neutral and positively ionized naphthalene. The calculated frequencies and intensities allowed an assignment of the most important bands appearing in the newly reported experimental spectrum of the positive ion. Agreement with the experimental spectrum is satisfactory enough to take into consideration the unexpected and important result that ionization significantly affects the intensities of most vibrations. A possible consequence on the interpretation of the IR interstellar emission, generally supposed to originate from polycyclic aromatic hydrocarbons (PAHs), is briefly presented.

Non-NASA Center↗

An Experimental and Theoretical Study of Nitrogen-Broadened Acetylene Lines

We present experimental nitrogen-broadening coefficients derived from Voigt profiles of isotropic Raman Q-lines measured in the 2 band of acetylene (C2H2) at 150 K and 298 K, and compare them to theoretical values obtained through calculations that were carried out specifically for this work. Namely, full classical calculations based on Gordon's approach, two kinds of semi-classical calculations based on Robert Bonamy method as well as full quantum dynamical calculations were performed. All the computations employed exactly the same ab initio potential energy surface for the C2H2N2 system which is, to our knowledge, the most realistic, accurate and up-to-date one. The resulting calculated collisional half-widths are in good agreement with the experimental ones only for the full classical and quantum dynamical methods. In addition, we have performed similar calculations for IR absorption lines and compared the results to bibliographic values. Results obtained with the full classical method are again in good agreement with the available room temperature experimental data. The quantum dynamical close-coupling calculations are too time consuming to provide a complete set of values and therefore have been performed only for the R(0) line of C2H2. The broadening coefficient obtained for this line at 173 K and 297 K also compares quite well with the available experimental data. The traditional Robert Bonamy semi-classical formalism, however, strongly overestimates the values of half-width for both Qand R-lines. The refined semi-classical Robert Bonamy method, first proposed for the calculations of pressure broadening coefficients of isotropic Raman lines, is also used for IR lines. By using this improved model that takes into account effects from line coupling, the calculated semi-classical widths are significantly reduced and closer to the measured ones.

coefficients↗

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.

NASA Discipline Exobiology↗

Ab initio calculations of anharmonic vibrational spectroscopy for hydrogen fluoride (HF)n (n = 3, 4) and mixed hydrogen fluoride/water (HF)n(H2O)n (n = 1, 2, 4) clusters

Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n, with n = 3, 4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n = 1, 2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the Moller-Plesset (MP2) potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.

Water/chemistry↗