Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional calculation”

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 37 records · Page 2

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

On the use of the noncentral chi-square density function for the distribution of helicopter spectral estimates

A probability density function for the variability of ensemble averaged spectral estimates from helicopter acoustic signals in Gaussian background noise was evaluated. Numerical methods for calculating the density function and for determining confidence limits were explored. Density functions were predicted for both synthesized and experimental data and compared with observed spectral estimate variability.

Garber, Donald P.↗

An Ab Initio Study of CuCO

Modified coupled-pair functional (MCPF) calculations and coupled cluster singles and doubles calculations, which include a perturbational estimate of the connected triples [CCSD(T)], yield a bent structure for CuCO, thus, supporting the prediction of a nonlinear structure based on density functional (DF) calculations. Our best estimate for the binding energy is 4.9 +/- 1.4 kcal/mol; this is in better agreement with experiment (6.0 +/- 1.2 kcal/mol) than the DF approach which yields a value (19.6 kcal/mol) significantly larger than experiment.

Bauschlicher, Charles W., Jr.↗

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↗

An analysis of spanwise gust gradient data

Plots of winds encountered during flight are presented for severe turbulence cases from flights 6 and 7 of the Joint Airport Weather Studies (JAWS) Project (McCarthy et al., 1982). During run 10 of flight 7, analyzed here in considerable detail, the B-57B Cambera aircraft showed a sharp increase in airspeed of 15 m/sec over a distance of about 130 m and then a more gradual decrease of 20 to 25 m/sec over a distance of about 2 km. This suspected outflow feature was associated with a downdraft exceeding 10 m/sec. The horizontal wind direction changed nearly 180 deg during the pass through the feature. Unexpectedly, the intensity of the velocity differences decreased for all three components within the downdraft. Calculated probability density functions for u, v, and w show a jagged character, whereas delta-u, delta-v, and delta-w show a much smoother, unimodal behavior. The distributions for delta-u, delta-v, and delta-w are distinctly non-Gaussian.

Campbell, W.↗

NASA B-57B Severe Storms Flight Program

Plots of winds encountered in-flight are presented for a severe turbulence case from JAWS flight 7 (near Denver on the afternoon of July 15, 1982). During the flight the B-57B showed a 30-knot increase in airspeed over a distance of about 426 ft and then a more gradual decrease of 40 to 50 knots over a distance of about 3.2 miles. This suspected outflow feature was associated with downdraft in excess of 20 knots. The horizontal wind direction changed almost 180 deg during the pass through the feature, and the intensity of velocity differences decreased for all three components within the downdraft. Calculated probability density functions for u, v, and w showed a jagged character; and distributions for Delta-u, Delta-v, and Delta-w were distinctly non-Gaussian. Based on skewness and kurtosis values, the data could probably best be modeled by a Pearson type VII (Student's t) distribution.

Painter, W. D.↗

Probabilistic structural analysis computer code (NESSUS)

Probabilistic structural analysis has been developed to analyze the effects of fluctuating loads, variable material properties, and uncertain analytical models especially for high performance structures such as SSME turbopump blades. The computer code NESSUS (Numerical Evaluation of Stochastic Structure Under Stress) was developed to serve as a primary computation tool for the characterization of the probabilistic structural response due to the stochastic environments by statistical description. The code consists of three major modules NESSUS/PRE, NESSUS/FEM, and NESSUS/FPI. NESSUS/PRE is a preprocessor which decomposes the spatially correlated random variables into a set of uncorrelated random variables using a modal analysis method. NESSUS/FEM is a finite element module which provides structural sensitivities to all the random variables considered. NESSUS/FPI is Fast Probability Integration method by which a cumulative distribution function or a probability density function is calculated.

Shiao, Michael C.↗

Electronic structure Fermi liquid theory of high Tc superconductors: Comparison of predictions with experiments

Predictions of local density functional (LDF) calculations of the electronic structure and transport properties of high T(sub c) superconductors are presented. As evidenced by the excellent agreement with both photoemission and positron annihilation experiments, a Fermi liquid nature of the 'normal' state of the high T(sub c) superconductors become clear for the metallic phase of these oxides. In addition, LDF predictions on the normal state transport properties are qualitatively in agreement with experiments on single crystals. It is emphasized that the signs of the Hall coefficients for the high T(sub c) superconductors are not consistent with the types of dopants (e.g., electron-doped or hole-doped) but are determined by the topology of the Fermi surfaces obtained from the LDF calculations.

Yu, Jaejun↗

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

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

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

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

Mid-Infrared Spectroscopy of Polycyclic Aromatic Nitrogen Heterocycles (PANHS) and their Ions

In recent years, polycyclic aromatic nitrogen heterocycles (PANHs) have attracted a good deal of attention because of their potent carcinogenic and mutagenic properties, and their prevalence in our environment. Such species also play a prominent role in the chemistry of life up to and including the very nucleobases from which our DNA is constructed. Surprisingly, these compounds may even be common outside of our terrestrial environment. To wit, it is now widely accepted that polycyclic aromatic materials are abundant in space and represent a major reservoir of organic carbon in the interstellar medium and developing planetary systems. Given that nitrogen is the fourth most abundant chemically reactive element in space (surpassed only by hydrogen, carbon, and oxygen), it is entirely reasonable to suspect that PANHs may represent an important component of that organic reservoir. Motivated by their intrinsic merit and with special attention toward evaluating their exobiological significance, we have initiated a program to study the spectroscopic and chemical properties of P A " s under conditions relevant to extraterrestrial environments. Here we present the first results of that program-infrared spectroscopic measurements on a series of PANH"s in neutral and cationic forms, isolated in inert matrices at cryogenic temperatures.temperatures. The species studied include: 1 -, and 2-azabenz[a]anthracene, 1-, 2-, and 4- azachrysene, dibenz[a,h]acridine, and dibenz[a,J)acridine. The experimental measurements are also compared with theoretical spectra calculated using density functional theory. General spectroscopic trends observed in this series of compounds are discussed and the implications of these results for Astrophysics and Exobiology are considered.

Mattioda, Andrew L.↗

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