Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional theory”

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 163 records · Page 9

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗

Theoretical Modeling of Formic Acid (HCOOH), Formate (HCOO(-)), and Ammonia (NH(4)) Vibrational Spectra in Astrophysical Ices

Ions embedded in icy grain mantles are thought to account for various observed infrared spectroscopic features, particularly in certain young stellar objects. The dissociation of formic acid (HCOOH) in astrophysical ices to form the formate ion (HCOO(-)) was modeled with density functional theory cluster calculations. Like isocyanic acid (HOCN), HCOOH was found to spontaneously deprotonate when sufficient water is present to stabilize charge transfer complexes. Both ammonia and water can serve as proton acceptors, yielding ammonium (NH4(+)) and hydronium (H3O(+)) counterions. Computed frequencies of weak infrared features produced by stretching and bending modes in both HCOO(-) and HCOOH were compared with experimental and astronomical data. Our results confirm laboratory assignments that a band at 1381 cm(exp -1) can be attributed to the CH bend in either HCOO(-) or HCOOH, but a band at 1349 cm(exp -1) corresponds to CO stretching in HCOO(-). Another feature at 1710 cm(exp -1) (5.85 m) can possibly be assigned to a CO stretching mode in HCOOH, as suggested by experiment, but the agreement is less satisfactory. In addition, we examine and analyze spectroscopic features associated with NH+4, both as a counterion to HCOO(-) or OCN(-) and in isolation, in order to compare with experimental and astronomical data in the 7 m region.

Park, Jin-Young↗

Quantum Chemical Evaluation of the Astrochemical Significance of Reactions between S Atom and Acetylene or Ethylene

Addition-elimination reactions of S atom in its P-3 ground state with acetylene (C2H2) and ethylene (C2H4) were characterized with both molecular orbital and density functional theory calculations employing correlation consistent basis sets in order to assess the likelihood either reaction might play a general role in astrochemistry or a specific role in the formation of S2 (X (sup 3 SIGMA (sub g) (sup -)) via a mechanism proposed by Saxena and Misra (Mon. Not. R. Astron. Soc. 1995, 272, 89). The acetylene and ethylene reactions proceed through C2H2S ((sup 3)A")) and C2H4S ((sup 3)A")) intermediates, respectively, to yield HCCS ((sup 2)II)) and C2H3S ((sup 2)A')). Substantial barriers were found in the exit channels for every combination of method and basis set considered in this work, which effectively precludes hydrogen elimination pathways for both S + C2H2 and S + C2H4 in the ultracold interstellar medium where only very modest barriers can be surmounted and processes without barriers tend to predominate. However, if one or both intermediates is formed and stabilized efficiently under cometary or dense interstellar cloud conditions, they could serve as temporary reservoirs for S atom and participate in reactions such as S + C2H2S (right arrow) S2 = C2H2 or S + C2H4S (right arrow) S2 + C2H4. For formation and stabilization to be efficient, the reaction must possess a barrier height small enough to be surmountable at low temperatures yet large enough to prevent redissociation to reactants. Barrier heights computed with B3LYP and large basis sets are very low, but more rigorous QCISD(T) and RCCSD(T) results indicate that the barrier heights are closer to 3-4 kcal/mol. The calculations therefore indicate that S + C2H2 or S + C2H4 could contribute to the formation of S2 in comets and may serve as a means to gauge coma temperature. The energetics of the ethylene reaction are more favorable.

Woon, David E.↗

Computational Confirmation of the Carrier for the "XCN" Interstellar Ice Bank: OCN(-) Charge Transfer Complexes

Recent experimental studies provide evidence that carrier for the so-called XCN feature at 2165 cm(exp -1) (4.62 micron) in young stellar objects is an OCN(-)/NH4(+) charge transfer (CT) complex that forms in energetically processed interstellar icy grain mantles. Although other RCN nitriles and RCN iosonitriles have been considered, Greenberg's conjecture that OCN(-) is associated with the XCN feature has persisted for over 15 years. In this work we report a computational investigation that thoroughly confirms the hypothesis that the XCN feature observed in laboratory studies can result from OCN(-)/NH4(+) CT complexes arising from HNCO and NH3, in a water ice environment. Density functional theory calculations with theory calculations with HNCO, NH3, and up to 12 waters reproduce seven spectroscopic measurements associated with XCN: the band origin of the asymmetric stretching mode of OCN(-), shifts due to isotopic substitutions of C, N, O, and H, and two weak features. However, very similar values are also found for the OCN(-)/NH4(+) CT complex arising from HOCN and NH3. In both cases, the complex forms by barrierless proton transfer from HNCO or HOCN to NH3 during the optimization of the solvated system. Scaled B3LYP/6-31+G** harmonic frequencies for HNCO and HOCN cases are 2181 and 2202 cm(exp -1), respectively.

Park, J.-Y.↗

Modeling Chemical Growth Processes in Titan's Atmosphere: 1. Theoretical Rates for Reactions between Benzene and the Ethynyl (C2H) and Cyano (CN) Radicals at Low Temperature and Pressure

Density functional theory calculations at the B3LYP/6-31+G** level were employed to characterize the critical points for adducts, isomers, products, and intervening transition states for the reactions between benzene and the ethynyl (C2H) or cyano (CN) radicals. Both addition reactions were found to have no barriers in their entrance channels, making them efficient at the low temperature and pressure conditions that prevail in the haze-forming region of Titan's atmosphere as well as in the dense interstellar medium (ISM). The dominant products are ethynylbenzene (C6H5C2H) and cyanobenzene (C6H5CN). Hydrogen abstraction reactions were also characterized but found to be non-competitive. Trajectory calculations based on potentials fit to about 600 points calculated at the ROMP2/6-31+G** level for each interaction surface were used to determine reaction rates. The rates incorporated any necessary corrections for back reactions as ascertained from a multiwell treatment used to determine outcome distributions over the range of temperatures and pressures pertinent to Titan and the ISM and are in good agreement with the limited available experimental data.

Woon, David E.↗

Conformational Structure of Tyrosine, Tyrosyl-Glycine, and Tyrosyl-Glycyl-Glycine by Double Resonance Spectroscopy

We investigated the variation in conformation for the amino acid tyrosine (Y), alone and in the small peptides tyrosine-glycine (YC) and tyrosine-glycine-glycine (YGG), in the gas phase by using UV-UV and IR-UV double resonance spectroscopy and density functional theory calculations. For tyrosine we found seven different conformations, for YG we found four different conformations, and for YGG we found three different conformations. As the peptides get larger, we observe fewer stable conformers, despite the increasing complexity and number of degrees of freedom. We find structural trends similar to those in phenylalanine-glycine glycine (FGG) and tryptophan-glycine-glycine (WGG)j however) the effect of dispersive forces in FGG for stabilizing a folded structure is replaced by that of hydrogen bonding in YGG.

Abo-Riziq, Ali↗

Ab initio Computations of the Electronic, Mechanical, and Thermal Properties of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2

Refractory materials such as metallic borides, often considered as ultra high temperature ceramics (UHTC), are characterized by high melting point, high hardness, and good chemical inertness. These materials have many applications which require high temperature materials that can operate with no or limited oxidation. Ab initio, first principles methods are the most accurate modeling approaches available and represent a parameter free description of the material based on the quantum mechanical equations. Using these methods, many of the intrinsic properties of these material can be obtained. We performed ab initio calculations based on density functional theory for the UHTC materials ZrB2 and HfB2. Computational results are presented for structural information (lattice constants, bond lengths, etc), electronic structure (bonding motifs, densities of states, band structure, etc), thermal quantities (phonon spectra, phonon densities of states, specific heat), as well as information about point defects such as vacancy and antisite formation energies.

Lawson, John W.↗

Association Mechanisms of Unsaturated C2 Hydrocarbons with Their Cations: Acetylene and Ethylene

The ion-molecule association mechanism of acetylene and ethylene with their cations is investigated by ab initio quantum chemical methods to understand the structures, association energies, and the vibrational and electronic spectra of the products. Stable puckered cyclic isomers are found as the result of first forming less stable linear and bridge isomers. The puckered cyclic complexes are calculated to be strongly bound, by 87, 35 and 56 kcal/mol for acetylene-acetylene cation, ethylene-ethylene cation and acetylene-ethylene cation, respectively. These stable complexes may be intermediates that participate in further association reactions. There are no association barriers, and no significant inter-conversion barriers, so the initial linear and bridge encounter complexes are unlikely to be observable. However, the energy gap between the bridged and cyclic puckered isomers greatly differs from complex to complex: it is 44 kcal/mol in C4H4 +, but only 6 kcal/mol in C4H8 +. The accurate CCSD(T) calculations summarized above are also compared against less computationally expensive MP2 and density functional theory (DFT) calculations for structures, relative energies, and vibrational spectra. Calculated vibrational spectra are compared against available experiments for cyclobutadiene cation. Electronic spectra are also calculated using time-dependent DFT.

ethylene↗

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Silicon-based ceramic components for next-generation jet turbine engines offer potential weight savings, as well as higher operating temperatures, both of which lead to increased efficiency and lower fuel costs. Silicon carbide (SiC), in particular, offers low density, good strength at high temperatures, and good oxidation resistance in dry air. However, reaction of SiC with high-temperature water vapor, as found in the hot section of jet turbine engines in operation, can cause rapid surface recession, which limits the lifetime of such components. Environmental Barrier Coatings (EBCs) are therefore needed if long component lifetime is to be achieved. Rare earth silicates such as Yb2Si2O7 and Yb2SiO5 have been proposed for such applications; in an effort to better understand diffusion in such materials, we have performed kinetic Monte Carlo (kMC) simulations of oxygen diffusion in Ytterbium disilicate, Yb2- Si2O7. The diffusive process is assumed to take place via the thermally activated hopping of oxygen atoms among oxygen vacancy sites or among interstitial sites. Migration barrier energies are computed using density functional theory (DFT).

Computer Simulation↗

Computational Design of Materials: Planetary Entry to Electric Aircraft and Beyond

NASA's projects and missions push the bounds of what is possible. To support the agency's work, materials development must stay on the cutting edge in order to keep pace. Today, researchers at NASA Ames Research Center perform multiscale modeling to aid the development of new materials and provide insight into existing ones. Multiscale modeling enables researchers to determine micro- and macroscale properties by connecting computational methods ranging from the atomic level (density functional theory, molecular dynamics) to the macroscale (finite element method). The output of one level is passed on as input to the next level, creating a powerful predictive model.

Materials Design↗

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Ytterbium disilicate is of interest as a potential environmental barrier coating for aerospace applications, notably for use in next generation jet turbine engines. In such applications, the transport of oxygen and water vapor through these coatings to the ceramic substrate is undesirable if high temperature oxidation is to be avoided. In an effort to understand the diffusion process in these materials, we have performed kinetic Monte Carlo simulations of vacancy-mediated and interstitial oxygen diffusion in Ytterbium disilicate. Oxygen vacancy and interstitial site energies, vacancy and interstitial formation energies, and migration barrier energies were computed using Density Functional Theory. We have found that, in the case of vacancy-mediated diffusion, many potential diffusion paths involve large barrier energies, but some paths have barrier energies smaller than one electron volt. However, computed vacancy formation energies suggest that the intrinsic vacancy concentration is small. In the case of interstitial diffusion, migration barrier energies are typically around one electron volt, but the interstitial defect formation energies are positive, with the result that the disilicate is unlikely to exhibit experience significant oxygen permeability except at very high temperature.

Coatings↗

Application of Ab Initio Methods in the Development of Advanced Technical Ceramics

A significant need exists to develop materials not only capable of providing desired electronic and mechanical properties but also survival at extreme temperatures during service and device fabrication, such as co-firing. Computational methods offers an efficient and systematic manner to design new materials and guide their development. As an example computational -based material approaches can be used to determine the suitability of a given materials as a practical thermoelectric for energy harvesting. In this presentation examples are given of applications to relevant technical ceramics such as thermoelectrics, dielectrics, and magnetic systems. The calculation were based on density functional theory and carried out with norm conserving and projector augmented wave (PAW) methods using commercial codes Materials Studio ( Biovia, Inc) with the Cambridge Serial Total Energy Package (CASTEP) and MedeA (Materials Design Inc.) utilizing the Vienna Ab-initio Simulation Package (VASP) as the respective computational engines. This study makes predictions of relevant technical properties of ceramic materials.

thermoelectric↗

Li(+)-ligand Binding Energies and the Effect of Ligand Fluorination on the Binding Energies

The Li(+)-ligand binding energies are computed for seven ligands and their perfluoro analogs using Density Functional Theory. The bonding is mostly electrostatic in origin. Thus the size of the binding energy tends to correlate with the ligand dipole moment, however, the charge-induced dipole contribution can be sufficiently large to affect the dipole-binding energy correlation. The perfluoro species are significantly less strongly bound than their parents, because the electron withdrawing power of the fluorine reduces the ligand dipole moment.

Charles W Bauschlicher↗

Aladyn - Adaptive Neural Network Molecular Dynamics Simulation Code: Computational Materials Mini-Application

This report provides an overview and commands description of the Computational Materials mini-application, Aladyn. Aladyn is a simple molecular dynamics code written in FORTRAN 2008, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of quantum mechanics-based energy calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. The Aladyn code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing many core architectures based on graphic processing units (GPUs). The effort is related to the High Performance Computing Incubator (HPCI) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗

Segregation and Phase Transformations Along Superlattice Intrinsic Stacking Faults in Ni-Based Superalloys

In this study, local chemical and structural changes along superlattice intrinsic stacking faults combine to represent an atomic-scale phase transformation. In order to elicit stacking fault shear, creep tests of two different single crystal Ni-base superalloys, ME501 and CMSX-4, were performed near 750 degrees Centigrade using stresses of 552 megapascals and 750 megapascals, respectively. Through high resolution Scanning Transmission Electron Microscopy (STEM) and state-of-the-art energy dispersive x-ray spectroscopy, ordered compositional changes were measured along SISFs (Superlattice Intrinsic Stacking Faults) in both alloys. For both instances, the elemental segregation and local crystal structure present along the SISFs are consistent with a nanoscale to D019 (ordering of a hexagonal close-packed crystal) phase transformation. Another notable observation is prominent Cr- and Co-rich Cottrell atmospheres and new evidence of more complex reordering processes responsible for the formation of these faults. These findings are further supported using density functional theory calculations and High Angle Annular Dark Field (HAADF) STEM image simulations.

STEM-EDS↗

The Effect Stacking Fault Segregation and Phase Transformations have on Creep Strength in Ni-based Superalloys

In this study, two commercially similar polycrystalline Ni-based disk superalloys (LSHR and ME3) were creep tested at 760C and 552MPa to approximately 0.3% plastic strain. LSHR consistently displayed superior creep properties at this stress/temperature regime even though the microstructural characteristics between the two alloys were comparable. High resolution structural and chemical analysis, however, revealed significant differences between the two alloys among active gamma prime shearing modes involving superlattice intrinsic and extrinsic stacking faults. In ME3, Co and Cr segregation and Ni and Al depletion were observed along the intrinsic faults - revealing a gamma prime -to- gamma phase transformation. Conversely in LSHR, an alloy with a higher W content, Co and W segregation was observed along the intrinsic faults. This observation combined with scanning transmission electron microscopy (STEM) simulations confirm a gamma prime-to-D019 phase transformation along the intrinsic faults in LSHR. Using experimental observations and density functional theory calculations, a novel local phase transformation strengthening mechanism is proposed that could be further utilized to improve the high temperature creep capabilities of Ni-base disk alloys.

Superalloy↗

Production and Potential Detection of Functionalized Hexamethylene-Tetramine Compounds in Space

Laboratory studies have shown that exposure of mixed ices of astrophysical interest to ionizing radiation such as ultraviolet (UV) photons or energetic particles (electrons, protons) leads to the production of large numbers of new, more complex compounds. A significant portion of these new species appear to belong to a family of molecules that consist of hexamethylenetetramine (HMT; C6N4H12) and HMT to which different chemical side groups have been substituted for a peripheral H atom. This work presents the identification of HMT-methanol (HMT-CH2OH), one of these HMT variants, in organic residues produced from the UV irradiation of astrophysically relevant ice mixtures at < 20 K. We also present the infrared (IR) spectra of HMT, HMTCH2OH, and a number of other HMT variants computed using density functional theory (DFT) computations. These spectra can be compared with each other and show similarities that can be used to search for this family of compounds in space.

UV Irradiation↗