Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory 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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Time Dependent Density Functional Theory Calculations of Large Compact PAH Cations: Implications for the Diffuse Interstellar Bands

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

Weisman, Jennifer L.↗
NASA NTRS2002

Electronic Absorption Spectra of Neutral Perylene (C20H12), Terrylene (C30H16), and Quaterrylene (C40H20) and their Positive and Negative Ions: Ne Matrix-Isolation Spectroscopy and Time Dependent Density Functional Theory Calculations

We present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.

Halasinski, Thomas M.↗
NASA NTRSDate not supplied

Does Chirality Influence the Stability of Amino Acid –Cu Complexes in the Salt-Induced Peptide Formation Reaction? Insights from Density Functional Theory Calculations

The polymerization of amino acids into peptides and ultimately proteins is critical to evolution of life. Ribosomally-synthesized proteins are homochiral, composed entirely of L-amino acids. Abiotic syntheses of amino acids, however, produce racemic mixtures(equal amounts of L and D enantiomers), meaning that the pool of monomers available for abiotic protein synthesis would be racemic or nearly so. There are no prebiotically plausible mechanisms known to yield homochiral peptides from racemic pools of amino acids. While some mechanisms capable of producing or amplifying enantiomeric excesses have been observed, these excesses are small or not applicable to all biologically relevant amino acids. As a result, a prebiotically-feasible mechanism to create enantiopure pools of amino acids has not been identified. However, it is not clear that these enantio-enriched or enantiopure pools are necessary to produce homochiral peptides. Steric and electronic effects, or a combination of both, may influence how different enantiomers interact with each other, potentially leading to homochiral peptides from racemic mixtures. Here, we use density functional theory (DFT) calculations to assess the stability of homochiral and heterochiral reactive complexes produced during polymerization via the salt-induced peptide formation (SIPF) reaction. Experimental studies of the SIPF reaction have shown that it enables the formation of peptides under diverse environmental conditions, making it a plausible pathway to peptide formation on early Earth. In the SIPF reaction, NaCl acts as a condensation reagent while a divalent metal cation, primarily Cu2+, facilitates polymerization by complexing and activating amino acids forming a monochlorocuprate complex. To assess if hetero-or homochiral monochlorocuprate complexes were energetically favored, we compared the stability of LL, DD, and LD enantiomers of Cu2+–(alanine)2and Cu2+–(valine) 2 complexes. Gaussian 09 density DFT energy-minimization calculations were made for each Cu –amino acid complex in the cis and trans configuration. Models were energy minimized using the BH and HLYP/6-31++G(d,p) level of theory and free energies were compared to determine the most stable configuration.

A C Fox↗
NASA NTRS1995

Using Density Functional Theory (DFT) for the Calculation of Atomization Energies

The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.

Bauschlicher, Charles W., Jr.↗
NASA NTRS2021

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗
NASA NTRS1998

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerene and Their Correlation with Three Geometric Properties: Symmetry, Aromaticity, and Size

The static third-order polarizabilities (gamma) of C60, C70, five isomers of C78 and two isomers of C84 were analyzed in terms of three properties, from a geometric point of view: symmetry, aromaticity and size. The polarizability values were based on the finite field approximation using a semiempirical Hamiltonian (AM1) and applied to molecular structures obtained from density functional theory calculations. Symmetry was characterized by the molecular group order. The selection of 6-member rings as aromatic was determined from an analysis of bond lengths. Maximum interatomic distance and surface area were the parameters considered with respect to size. Based on triple linear regression analysis, it was found that the static linear polarizability (alpha) and gamma in these molecules respond differently to geometrical properties: alpha depends almost exclusively on surface area while gamma is affected by a combination of number of aromatic rings, length and group order, in decreasing importance. In the case of alpha, valence electron contributions provide the same information as all-electron estimates. For gamma, the best correlation coefficients are obtained when all-electron estimates are used and when the dependent parameter is ln(gamma) instead of gamma.

Moore, C. E.↗
NASA NTRS2000

Theoretical Study of Indium Compounds of Interest for Organometallic Chemical Vapor Deposition

The structural. electronic and therinochemical properties of indium compounds which are of interest in halide transport and organometallic chemical vapor deposition processes have been studied by ab initio and statistical mechanics methods. The compounds reported include: indium halides and hydrides (InF, InCl, InCl3, InH, InH2, InH3); indium clusters (In2, In3); methylindium, dimethylindium, and their hydrogen derivatives [In(CH3), In(CH3)H, In(CH3)H2, In(CH3)2, In(CH3)2H]; dimethyl-indium dimer [In2(CH3)4], trimethyl-indium [In(CH3)3]; dehydrogenated methyl, dimethyl and trimethylindium [In(CH3)2CH2, In(CH3)CH2, In(CH2)], trimethylindium adducts with ammonia, trimethylamine and hydrazine [(CH3)3In:NH3, (CH3)3In:N(CH3)3, (CH3)3In:N(H2)N(H2)]; dimethylamino-indium and methylimino-indium [In(CH3)2(NH2), In(CH3)(NH)]; indium nitride and indium nitride dimer (InN, In2N2), indium phosphide, arsenide and antimonide ([InP, InAs, InSb). The predicted electronic properties are based on density functional theory calculations; the calculated thermodynamic properties are reported following the format of the JANAF (Joint Army, Navy, NASA, Air Force) Tables. Equilibrium compositions at two temperatures (298 and 1000 K) have been analyzed for groups of competing simultaneous reactions.

Cardelino, B. H.↗
NASA NTRS1998

Energetics and Dynamics of GaAs Epitaxial Growth via Quantum Wave Packet Studies

The dynamics of As(sub 2) molecule incorporation into the flat Ga-terminated GaAs(100) surface is studied computationally. The time-dependent Schrodinger equation is solved on a two-dimensional potential energy surface obtained using density functional theory calculations. The probabilities of trapping and subsequent dissociation of the molecular As(sub 2) bond are calculated as a function of beam translational energy and vibrational quantum number of As(sub 2).

Dzegilenko, Fedor N.↗
NASA NTRS2003

The Infrared Spectrum of Matrix Isolated Aminoacetonitrile: A Precursor to the Amino Acid Glycine

We present infrared (IR) spectral data from matrix isolation experiments and density functional theory calculations on the pre-biologically interesting molecule aminoacetonitrile, a precursor to glycine. We find that this nitrile has an unusually weak nitrile (C=N) stretch in the infrared, in contrast to expectations based on measurements and models of other nitriles under astrophysical conditions. The absence of an observable nitrile absorption feature in the infrared will make the IR search for this molecule considerably more difficult, and will raise estimates of upper limits on nitriles in interstellar and outer Solar System ices. This is also of relevance to assessing the formation routes of the amino acid glycine, since aminoacetonitrile is the putative precursor to glycine via the Strecker synthesis, the mechanism postulated to have produced the amino acids in meteorites.

Bernstein, Max P.↗
NASA NTRS2007

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

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

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

Bonding and Microstructural Stability in Ni55Ti45 Studied by Experimental and Theoretical Methods

Spiral orbit tribometry friction tests performed on Ni-rich Ni55Ti45 titanium ball bearings indicate that this alloy is a promising candidate for future aerospace bearing applications. Microstructural characterization of the bearing specimens was performed using transmission electron microscopy and energy dispersive spectroscopy, with NiTi, Ni4Ti3, Ni3Ti, and Ni2Ti4Ox phases identified within the microstructure of the alloy. Density functional theory was applied to predict the electronic structure of the NixTiy phases, including the band structure and site projected density of states. Ultraviolet photoemission spectroscopy was used to verify the density of states results from the density functional theory calculations, with good agreement observed between experiment and theory.

Stott, Amanda C.↗
NASA NTRS2011

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↗
NASA NTRS2018

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↗
NASA NTRS2019

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↗
NASA NTRSDate not supplied

Modeling Ring-Opening Polymerization with Machine-Learned Committors

Many important industrial and biological processes depend on ring-opening polymerization (ROP), including large-scale production polymers, such as nylon 6, and a variety of biodegradable polymers and high-performance resins. Large-scale simulations (>200K atoms) of ROP were performed using a deep neural network-based method that captured the effect of the local chemical environment of the reactive sites using training data from quantum chemical methods. The utility of this method was demonstrated by modeling the cationic ROP of bis-benzoxazine, a promising resin for use in advanced high-temperature application composites (>200°C). Density functional theory calculations were used to map out the various ring opening and polymerization reaction paths that occur during the thermal processing of bis-benzoxazine. The probabilities of these reactions, as a function of local atomic configurations, were used to construct committor functions using a neural network approach. The learned committor functions were then used to determine when a particular reaction occurred within the REACTER framework for reactive molecular dynamics simulations. The experimentally observed exotherm that occurs during bis-benzoxazine ROP was captured in an effective manner in the simulations by imposing a local heating of the system that corresponds to the enthalpy of reaction for each reaction. The dependence of the final polymer morphology on the ROP reaction temperature was analyzed at the molecular scale. The results presented in this work indicate that machine-learned committor functions are a promising approach for incorporating high-fidelity reactivity criteria that capture the effect of local chemical environments into large-scale polymerization simulations.

machine learning, polymer simulations, molecular d↗
NASA NTRS2023

An Ensemble Neural Network Model for Predicting Rare-Earth Oxide and Silicate Heat Capacities at High Temperature

In this work, a neural network model was developed to predict the constant pressure heat capacity for materials in the rare-earth oxide—silica material space. Several model architectures were trained and tested on heat capacity data generated from first-principles density functional theory calculations. Hyperparameter optimization was performed, and the optimal model was selected for heat capacity predictions. The optimal model architecture was found to have a root-mean-squared error of 5.12 ± 3.37 J/mol-K. The optimal model architecture was then used in a bagging ensemble model trained using the leave-one-group-out method to provide error estimates for model predictions. The out-of-bag score for the ensemble model was 0.997. The predicted heat capacities agree well with the DFT and experimental results and were computed orders of magnitude faster than DFT simulations. Machine learning shows the potential to provide a suitable surrogate model for thermochemical property predictions for candidate environmental barrier coating materials but refining of input material features and model architectures could further improve accuracy for these models.

environmental barrier coatings↗