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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Molecular Simulations of Ultrafast Radiation Induced Melting at Metal-Semiconductor Interfaces

Understanding radiation induced ultrafast melting at material interfaces is essential in designing robust electronic devices for aviation/space applications and in laser machining. While it is difficult to achieve the spatial and temporal resolution required to quantify the phenomenon experimentally, simulations can provide the detailed mechanisms of the structural changes that happen during phase transition. In this work, we use molecular simulations to study the effect of radiation damage on silicon carbide (SiC) - tungsten (W) interfaces which is of interest in high power electronics. A multi-scale approach is involved wherein the reactions at the interfaces are quantified using ab-initio molecular dynamics (MD) simulations and classical MD simulations are employed to understand the structural and diffusional changes across the material interface. Finally, coarse-grained Lennard-Jones type models are used to study the larger scale mechanisms and structures obtained due to the induced damages. We show that the response of the material to radiation damage depends on factors such as energy of the incident radiation, thermal properties, and molecular structure of the material.

Ravichandran, Ashwin↗
NASA NTRS1986

Simulation studies for surfaces and materials strength

During this reporting period three investigations were carried out. The first area of research concerned the analysis of the structure-energy relationship in small clusters. This study is very closely related to the improvement of the potential energy functions which are suitable and simple enough to be used in atomistic simulation studies. Parameters obtained from ab initio calculations for dimers and trimers of Al were used to estimate energetics and global minimum energy structures of clusters continuing up to 15 Al atoms. The second research topic addressed modeling of the collision process for atoms impinging on surfaces. In this simulation study qualitative aspects of the O atom collision with a graphite surface were analyzed. Four different O/graphite systems were considered and the aftermath of the impact was analyzed. The final area of investigation was related to the simulation of thin amorphous Si films on crystalline Si substrates. Parameters obtained in an earlier study were used to model an exposed amorphous Si surface and an a-Si/c-Si interface. Structural details for various film thicknesses were investigated at an atomistic level.

Halicioglu, T.↗
NASA NTRS2016

Computing Highly Accurate Spectroscopic Line Lists for Characterization of Planetary Atmospheres: CO2 and SO2 Line Lists Needed for Modeling Venus

Over the last decade, it has become apparent that the most effective approach for determining highly accurate rotational and rovibrational line lists for molecules of interest in planetary atmospheres and other astrophysical environments is through a combination of highresolution laboratory experiments coupled with state-of-the art ab initio quantum chemistry methods. The approach involves computing the most accurate potential energy surface (PES) possible using state-of-the art electronic structure methods, followed by computing rotational and rovibrational energy levels using an exact variational method to solve the nuclear Schrödinger equation. Then, reliable experimental data from high-resolution experiments is used to refine the ab initio PES in order to improve the accuracy of the computed energy levels and transition energies. From the refinement step, we have been able to achieve an accuracy of approximately 0.015 cm-1 for rovibrational transition energies, and even better for purely rotational transitions. This combined "experiment / theory" approach allows for determination of essentially a complete line list, with hundreds of millions of transitions, and having the transition energies and intensities be highly accurate. Our group has successfully applied this approach to determine highly accurate line lists for NH3, CO2 and isotopologues, and SO2 and isotopologues. Here I will report our latest results for CO2 and SO2 including all isotopologues. Comparisons to the available data in HITRAN2012 and other available databases will be shown, though we note that our line lists for SO2 are significantly more complete than any other databases. Since it is important to span a large temperature range in order to model the spectral signature of Venus as well as exoplanets, we will demonstrate how the spectra change on going from low temperatures (100 K) to higher temperatures (500 K to 1500 K).

potential energy surface (PES)5405 PLANETARY SCIEN↗
NASA NTRS2002

Towards accurate ab initio predictions of the vibrational spectrum of methane

We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born-Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.

Spectrophotometry, Infrared/methods↗
NASA NTRS1998

Simulated infrared emission spectra of highly excited polyatomic molecules: a detailed model of the PAH-UIR hypothesis

A detailed description of the polycyclic aromatic hydrocarbon (PAH)/unidentified infrared band (UIR) mechanism is presented in which experimental spectral bandshape functions are used to simulate IR emission spectra for individual molecules. These spectra are additively superimposed to produce a conglomerate spectrum representative of a family of PAH molecules. Ab initio vibrational frequencies and intensities for nine PAHs (neutral and cationic) as large as ovalene are used in conjunction with measured bandshape and temperature-dependent redshift data to simulate the UIR bands. The calculated spectra of cations provide a closer match to the UIRs than do those of the neutrals. However, the PAH cations used in the simulations fail to reproduce the details of the UIR emission spectra. The discrepancies are potentially alleviated if both larger PAHs and a greater number of PAHs were included in the simulation.

Non-NASA Center↗
NASA NTRS1998

Modeling Intermolecular Interactions in Nanotubes, Fullerenes and Graphite using a New Long-Range Potential

The cohesive energy and compressibility of strands of a single-wall nanotube rope has been computed using a new long-range potential energy function derived from accurate ab initio quantum chemistry calculations of the benzene dimer and calibrated for energetic and mechanical properties of graphite (at pressures up to 12 GPa). We also use this potential to calculate a variety of properties of carbon nanotubes (both single- and multi-wall) and fullerenes. Extensive comparisons are made with previously published potentials.

Jaffe, Richard↗
NASA NTRS1997

Rate Coefficients for Reactions of Ethynyl Radical (C2H) With HCN and CH3CN: Implications for the Formation of Comples Nitriles on Titan

Rate coefficients for the reactions of C2H + HCN yields products and C2H + CH3CN yields products have been measured over the temperature range 262-360 K. These experiments represent an ongoing effort to accurately measure reaction rate coefficients of the ethynyl radical, C2H, relevant to planetary atmospheres such as those of Jupiter and Saturn and its satellite Titan. Laser photolysis of C2H2 is used to produce C2H, and transient infrared laser absorption is employed to measure the decay of C2H to obtain the subsequent reaction rates in a transverse flow cell. Rate constants for the reaction C2H + HCN yields products are found to increase significantly with increasing temperature and are measured to be (3.9-6.2) x 10(exp 13) cm(exp 3) molecules(exp -1) s(exp -1) over the temperature range of 297-360 K. The rate constants for the reaction C2H + CH3CN yields products are also found to increase substantially with increasing temperature and are measured to be (1.0-2.1) x 10(exp -12) cm(exp 3) molecules(exp -1) s(exp -1) over the temperature range of 262-360 K. For the reaction C2H + HCN yields products, ab initio calculations of transition state structures are used to infer that the major products form via an addition/elimination pathway. The measured rate constants for the reaction of C2H + HCN yields products are significantly smaller than values currently employed in photochemical models of Titan, which will affect the HC3N distribution.

Hoobler, Ray J.↗
NASA NTRS2012

Development and Application of Interatomic Potentials for Ultra High Temperature Ceramics (UHTC): ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are characterized by high melting point, good strength, and reasonable oxidation resistance. These materials are of interest for use as sharp leading edges for hypersonic vehicles among other applications. Progress in computational modeling of UHTCs has been limited in part due to the absence of suitable interatomic potentials. We present a Tersoff style parameterization of such potentials for ZrB2 and HfB2 appropriate for atomistic simulations. Parameters are fit to data generated from ab initio calculations. The accuracy of the potentials is assessed against further ab initio data. As a first non ]trivial application, molecular dynamics simulations are performed to evaluate the thermal conductivity of single crystals and the thermal resistance of high symmetry grain boundaries.

Lawson, John W.↗
NASA NTRS2004

MGGPOD: a Monte Carlo Suite for Modeling Instrumental Line and Continuum Backgrounds in Gamma-Ray Astronomy

Intense and complex instrumental backgrounds, against which the much smaller signals from celestial sources have to be discerned, are a notorious problem for low and intermediate energy gamma-ray astronomy (approximately 50 keV - 10 MeV). Therefore a detailed qualitative and quantitative understanding of instrumental line and continuum backgrounds is crucial for most stages of gamma-ray astronomy missions, ranging from the design and development of new instrumentation through performance prediction to data reduction. We have developed MGGPOD, a user-friendly suite of Monte Carlo codes built around the widely used GEANT (Version 3.21) package, to simulate ab initio the physical processes relevant for the production of instrumental backgrounds. These include the build-up and delayed decay of radioactive isotopes as well as the prompt de-excitation of excited nuclei, both of which give rise to a plethora of instrumental gamma-ray background lines in addition t o continuum backgrounds. The MGGPOD package and documentation are publicly available for download. We demonstrate the capabilities of the MGGPOD suite by modeling high resolution gamma-ray spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board Wind during 1995. The TGRS is a Ge spectrometer operating in the 40 keV to 8 MeV range. Due to its fine energy resolution, these spectra reveal the complex instrumental background in formidable detail, particularly the many prompt and delayed gamma-ray lines. We evaluate the successes and failures of the MGGPOD package in reproducing TGRS data, and provide identifications for the numerous instrumental lines.

Weidenspointner, G.↗
NASA NTRS1985

The lower electronic states of ClOO - A computational investigation

Eight doublet and eight quartet states of ClOO were investigated by ab initio CI techniques. The potential energy surfaces of the four lowest energy doublet states of both A-double-prime and A-prime symmetry indicate that only the 1 2A-double-prime state is bound. In contrast to the model provided by the HO2 radical, all of the excited doublet states investigated were repulsive with respect to dissociation to Cl + O2 and metastable or bound with respect to dissociation to ClO + O. The transitions to the excited states investigated span the visible and near UV spectral regions, but the transition moments indicate that they are very weak. Since the photolysis products are the same as those of the rapid thermal dissociation, photolysis is not expected to be an important atmospheric process. The soft bending potential for the 1 2A-prime state and the shape of the 1 4A-double-prime state in the entrance channel of the ClO + O yields Cl + O2 reaction provide a qualitative explanation for the underprediction of the low temperature reaction rate by previous trajectory calculations.

Jafri, J. A.↗
NASA NTRS1993

Theoretical infrared spectra of some model polycyclic aromatic hydrocarbons - Effect of ionization

In order to test the hypothesis of ionized PAHs as possible carriers of the UIR bands, we realized a computational exploration on selected PAHs of small dimension in order to identify which changes ionization would induce on their IR spectra. In this study we performed ab initio calculations of the spectra of neutral and positively ionized naphthalene, anthracene, and pyrene. The results are significantly important. The frequencies in the cations are slightly shifted with respect to the neutral species, but no general conclusion can be reached from the three molecules considered. By contrast, the relative intensities of most vibrations are strongly affected by ionization, leading to a much better agreement between the calculated CH/CC vibration intensity ratios and those deduced from observations.

De Frees, D. J.↗
NASA NTRSDate not supplied

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗
NASA NTRS2005

[Ti II] and [Ni II] Emission from the Strontium Filament of eta Carinae

We study the nature of the [Ti II] and [Ni II] emission from the so-called strontium filament found in the ejecta of eta Carinae. To this purpose we employ multilevel models of the Ti II and Ni II systems which are used to investigate the physical condition of the filament and the excitation mechanisms of the observed lines. For the Ti II ion, for which no atomic data was previously available, we carry out ab initio calculations of radiative transition rates and electron impact excitation rate coefficients. It is found that the observed spectrum is consistent with the lines being excited in a mostly neutral region with electron density of the order of 10(exp 7) cm(exp -3) and a temperature around 6000 K. In analyzing three observations with different slit orientations recorded between March 2000 and November 2001 we find line ratios that change among various observations, in a way consistent with changes of up to an order of magnitude in the strength of the continuum radiation field. These changes result from different samplings of the extended filament, due to the different slit orientations used for each observation, and yield clues on the spatial extent and optical depth of the filament. The observed emission indicates a large Ti/Ni abundance ratio relative to solar abundances. It is suggested that the observed high Ti/Ni ratio in gas is caused dust-gas fractionation processes and does not reflect the absolute Ti/Ni ratio in the ejecta of eta Carinae. The condensation chemistry shows that if dust condensed in a sequence of layers according to decreasing temperature and increasing distance from the central star, the most refractory dust could be selectively affected by photoevaporation. Thus, Ti would be released back to the gas and the Ti/Ni ratio in the gas would increase to the observed super-solar ratio.

Bautista, M. A.↗
NASA NTRS2007

[Ti II] and [Ni II] Emission from the Strontium Filament of eta Carinae

We study the nature of the [Ti II] and [Ni II] emission from the so-called strontium filament found in the ejecta of eta Carinae. To this purpose we employ multilevel models of the Ti II and Ni II systems which are used to investigate the physical condition of the filament and the excitation mechanisms of the observed lines. For the Ti II ion, for which no atomic data was previously available, we carry out ab initio calculations of radiative transition rates and electron impact excitation rate coefficients. It is found that the observed spectrum is consistent with the lines being excited in a mostly neutral region with an electron density of the order of 10(exp 7) per cubic centimeter and a temperature around 6000 K. In analyzing three observations with different slit orientations recorded between March 2000 and November 2001 we find line ratios that change among various observations, in a way consistent with changes of up to an order of magnitude in the strength of the continuum radiation field. These changes result from different samplings of the extended filament, due to the different slit orientations used for each observation, and yield clues on the spatial extent and optical depth of the filament. The observed emission indicates a large Ti/Ni abundance ratio relative to solar abundances. It is suggested that the observed high Ti/Ni ratio in gas is caused by dust-gas fractionation processes and does not reflect the absolute Ti/Ni ratio in the ejecta of eta Carinae. We study the condensation chemistry of Ti, Ni and Fe within the filament and suggest that the observed gas phase overabundance of Ti is likely the result of selective photo-evaporation of Ti-bearing grains. Some mechanisms for such a scenario are proposed.

Bautista, M. A.↗
NASA NTRSDate not supplied

Utilizing Experimental-based Line Positions for Semi-empirical IR Line Lists

Molecular IR line lists computed from semi-empirically refined ab initio potential energy surface and high quality ab initio dipole moment surface may have line position accuracy typically within 0.01 - 0.05 cm−1. For high resolution spectroscopy databases, it is necessary to integrate the computed theoretical IR intensities with the line positions accurately determined from experiments and (or) reliably derived from Effective Hamiltonian models. However, it is not a trivial task, and the choices in reality are heavily contingent upon the coverage, consistency, and accuracy of the data available for a specific molecule or isotopologue. We will present several approaches applied to recent Ames IR line lists of CO 2 1, N 2 O 2 , and OCS to demonstrate how challenging such integration may become, and what we have learned From these projects. Discussions will focus on the relation between each different scenario and the corresponding choice or solution, including their advantages and limits, how the semi-empirical IR line lists can help, and what else may be needed for future improvements. We will emphasize that the “Theory+Experiment” synergy may still play significant role in the determination of the best line positions.

Effective Hamiltonian models↗
NASA NTRS1976

Electronic states of Zn2 - Ab initio calculations of a prototype for Hg2

The electronic states of Zn2 are investigated by ab initio polarization configuration-interaction calculations. Molecular states dissociating to Zn(1S) + Zn(1S, 3P, 1P) and Zn(3P) + Zn(3P) are treated. Important effects from states arising from Zn(+)(25) + Zn(-)(2P) are found in the potential-energy curves and electronic-transition moments. A model calculation for Hg2 based on the Zn2 curves and including spin-orbit coupling leads to a new interpretation of the emission bands in Hg vapor.

Hay, P. J.↗
NASA NTRS2002

Dynamics of a Finite Liquid Oxygen (LOX) Column in a Pulsed Magnetic Field

It is well known that liquid oxygen has a sufficient paramagnetic susceptibility that a strong magnetic field gradient can lift it in the earth's gravitational field. The movement of liquid oxygen is vital to the space program since it one of the primary oxidizers used for propulsion. Transport of liquid oxygen (LOX) via direct interaction of the magnetic fields (B field) with the fluid is a current topic of research and development at Kennedy Space Center, FL. This method of transporting (i.e. pumping) LOX may have particular advantages on Mars and other reduced gravitational environments, namely safety and reliability. This paper will address transport of a magnetic fluid, LOX, via phased-pulsed electromagnets acting on the edge of the column of fluid. The authors have developed a physical model from first-principles for the motion of a magnetic fluid in a particular U-tube geometry subjected to a pulsed magnetic field from an arbitrary solenoidal electromagnet. Experimental data that have been collected from the analogous geometry correlate well to that of the ab-initio calculations.

Youngquist, Robert↗
NASA NTRS1993

Ab initio method for calculating total cross sections

A method for calculating total cross sections without formally including nonelastic channels is presented. The idea is to use a one channel T-matrix variational principle with a complex correlation function. The derived T matrix is therefore not unitary. Elastic scattering is calculated from T-parallel-squared, but total scattering is derived from the imaginary part of T using the optical theorem. The method is applied to the spherically symmetric model of electron-hydrogen scattering. No spurious structure arises; results for sigma(el) and sigma(total) are in excellent agreement with calculations of Callaway and Oza (1984). The method has wide potential applicability.

Bhatia, A. K.↗