Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 541 records · Page 30

An Integrated Study of a Novel Thermal Coating for Nb-Based High Temperature Alloy

This report summarizes our recent works of ab initio density functional theory (DFT) method and molecular dynamics (MD) simulation on the interfaces between niobium substrate and coatings at atomic level. Potential oxidation barrier bond coat, Nb₂AlC and high entropy alloys, and top coat candidates were synthesized, characterized, and evaluated in our labs. The simulation methods, experimental validation techniques, achievements already reached, students and postdoc training, and future improvement are briefly introduced.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nuclear and neutron-star matter from local chiral interactions [Slides]

We use quantum Monte Carlo (QMC) methods to perform exact ab-initio calcula1ons of the nuclear equation of state (EOS) and the symmetry energy. We employ local chiral interactions up to next-to-next-to leading (N 2 LO) fit to few-body observables only, and provide a comprehensive uncertainty quantifica1on.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Photodissociation of CH2. I - Potential energy surfaces of the dissociation into CH and H

The possible photodissociation pathways of the CH2 radical are studied using ab initio multireference configuration-interaction methods, and accurate photodissociation cross sections and branching ratios for the production of CH + H and C + H2 are obtained. Potential energy surfaces were calculated using the Wuppertal-Bonn self-consistent field plus a multireference single and double-excitation configuration interaction package of programs. Two-dimensional potential energy surfaces of the ten lowest triplet states correlating with the seven lowest states of CH were calculated as functions of bond angle and one C-H bond distance, keeping the other C-H bond distance fixed at the equilibrium CH2 value.

Bearda, Robert A.↗

Potential energy curves and collision integrals of air components

Collision integrals are fundamental quantities required to determine the transport properties of the environment surrounding aerospace vehicles in the upper atmosphere. These collision integrals can be determined as a function of temperature from the potential energy curves describing the atomic and molecular collisions. Ab initio calculations provide a practical method of computing the required interaction potentials. In this work we will discuss recent advances with an emphasis on the accuracy that is obtainable. Results for interactions, e.g. N+N, N+O, O+O, and H+N2 will be reviewed and their application to the determination of transport properties, such as diffusion and viscosity coefficients, will be examined.

Partridge, Harry↗

Potential Energy Curves and Collisions Integrals of Air Components: Interactions Involving Ionized Atoms - 2

Collision integrals are fundamental quantities required to determine the transport properties of the environment surrounding aerospace vehicles in the upper atmosphere. These collision integrals can be determined as a function of temperature from the potential energy curves describing the atomic and molecular collisions. Ab initio calculations provide a practical method of computing the required interaction potentials. In this work we will discuss recent advances in scattering calculations with an emphasis on the accuracy that is obtainable. Results for interactions of the atoms and ionized atoms of nitrogen and oxygen will be reviewed and their application to the determination of transport properties, such as diffusion and viscosity coefficients, will be examined.

Stallcop, James R.↗

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗

Magnetic properties of CrFeCoNi based high entropy alloys

Monte Carlo simulations are performed on three high entropy alloys: Cr0.25Fe0.25Co0.25Ni0.25, Cr0.2Fe0.2Co0.2Ni0.2Pd0.2, and Cr0.2Mn0.2Fe0.2Co0.2Ni0.2, with exchange interactions extracted from The ab initio Korringa-Kohn-Rostoker method combined with the coherent potential approximation calculations. Using finite size scaling analyses, we estimate the magnetic phase transition temperature for the four component alloy to be 108(2) K, and although the individual critical exponents are different from 3D Heisenberg universality class, the reduced exponent follows Suzuki weak universality. With the additional Palladium component, the transition temperature elevates to about 200 K. In contrast, we find no magnetic order for the five component alloy with Manganese at any finite temperatures.

Yin, Junqi↗

Thermal scattering law for structure-dependent-Doppler broadening in FLASSH

Traditional Doppler broadening methods assume that target nuclei velocities follow a Maxwell-Boltzmann distribution which inherently assumes that the atoms are a free gas or that their velocities are independent of each other. This approximation is reasonable at high target temperatures and high neutron energies. However, nuclei are not independent: the lattice temperature and chemical binding will define the atomic motion. The thermal scattering law (i.e., TSL or S(α,β)) is a material property which describes the chemical binding and temperature response in terms of available momentum and energy states. In the thermal energy range, this TSL will define the thermal scattering cross sections. In the resonance region, the TSL offers a first-principles description of the probability distribution function for the velocity of the target nuclei. Using the TSL in Doppler broadening captures the structure of the material consistently from the thermal range into the resonance cross sections. In this work, both traditional free-gas and precise TSL Doppler broadening methods for resonance cross section evaluation have been implemented within the Full Law Analysis Scattering System Hub (FLASSH). This framework provides a generalized formulation for evaluating low-lying resonance data and streamlines nuclear data processing methods. TSL data generated using ab initio lattice dynamics (AILD) can be used to broaden cross section libraries in ENDF format for high-fidelity input into reactor physics calculations. These capabilities are demonstrated for the lowest absorption resonances of {sup 238}U in UO{sub 2}, UC and UN. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bridging a gap between continuum-QCD and ab initio predictions of hadron observables

Within contemporary hadron physics there are two common methods for determining the momentum-dependence of the interaction between quarks: the top-down approach, which works toward an ab initio computation of the interaction via direct analysis of the gauge-sector gap equations; and the bottom-up scheme, which aims to infer the interaction by fitting data within a well-defined truncation of those equations in the matter sector that are relevant to bound-state properties. We unite these two approaches by demonstrating that the renormalisation-group-invariant running-interaction predicted by contemporary analyses of QCD's gauge sector coincides with that required in order to describe ground-state hadron observables using a nonperturbative truncation of QCD's Dyson–Schwinger equations in the matter sector. This bridges a gap that had lain between nonperturbative continuum-QCD and the ab initio prediction of bound-state properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Iron corrosion in the “inert” supercritical CO 2 , ab initio dynamics insights: How impurities matter

Metal corrosion from exposure to supercritical fluids is a long-standing challenge lacking atomistic understanding for effective corrosion prevention in diverse important applications. Here, in this work, reactions relevant to corrosion in supercritical CO 2 (sCO 2 ), with H 2 O and NO 2 as impurities, are revealed using ab initio molecular dynamics and enhanced sampling methods. While little reactivity is observed in bulk sCO 2 , the iron surface is found to play a crucial almost-catalytic role in its own corrosion, activating the inert CO2 and hence facilitating $\ast$HCO 3 – and $\ast$CO 3 2– formation at the interface. Moreover, the hydrogen bond network of $\ast$H 2 O assists the H shuttling to the activated CO 2 and impurities, forming reactive species observed in experiments. Such atomistic insight elucidates the reaction mechanisms of metal corrosion in aqueous and supercritical fluids and can point to possible remedies.

36 MATERIALS SCIENCE↗

An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl 2

Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl 2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl 2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. Here a wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl 2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl 2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl 2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl 2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl 2 . Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.

36 MATERIALS SCIENCE↗

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48 Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. In this work, we use our scheme to compute the matrix element that governs the neutrinoless double beta decay of 48 Ca to 48 Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76 Ge, 130 Te, and 136 Xe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Second spectrum of chromium (Cr II), Part III: Radiative lifetimes and transition probabilities from highly excited 3d{sup 4}5s levels

In this study we have extended our previous investigations of Cr II spectrum to levels of the 3d{sup 4}5s configuration. Once again an overall good agreement is observed for these levels between experimental oscillator strength values and our calculated data using two different approaches, namely the oscillator strength parameterization and the HFR+CPOL methods. Recurring to the latter model we have computed radiative lifetime values for 20 3d{sup 4}5s levels, confirming the well founded basis of recent experimental data, given in literature. From the use of the former method and with the help of a least squares fitting procedure to available experimental gf-values we have taken advantage of extracting for the first time the transition radial integral value for 3d{sup 4}4p−3d{sup 4}5s, 〈4p|r{sup 1}|5s〉=1.962±0.005, which is confirmed favorably by ab initio results, obtained by the HFR method scaled by a factor 0.92 or directly by the HFR+CPOL approach. Finally a long list of semi-empirical oscillator strength, transition probability and branching fraction values is generated for around 500 Cr II transitions depopulating the 3d{sup 4}5s levels. Our new gf data are compared successfully with those given by Kurucz.

74 ATOMIC AND MOLECULAR PHYSICS↗

A permutationally invariant full-dimensional ab initio potential energy surface for the abstraction and exchange channels of the H + CH{sub 4} system

We report a permutationally invariant global potential energy surface (PES) for the H + CH{sub 4} system based on ∼63 000 data points calculated at a high ab initio level (UCCSD(T)-F12a/AVTZ) using the recently proposed permutation invariant polynomial-neural network method. The small fitting error (5.1 meV) indicates a faithful representation of the ab initio points over a large configuration space. The rate coefficients calculated on the PES using tunneling corrected transition-state theory and quasi-classical trajectory are found to agree well with the available experimental and previous quantum dynamical results. The calculated total reaction probabilities (J{sub tot} = 0) including the abstraction and exchange channels using the new potential by a reduced dimensional quantum dynamic method are essentially the same as those on the Xu-Chen-Zhang PES [Chin. J. Chem. Phys. 27, 373 (2014)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Utilizing Amino Acid Composition and Entropy of Potential Open Reading Frames to Identify Protein-Coding Genes

One of the main steps in gene-finding in prokaryotes is determining which open reading frames encode for a protein, and which occur by chance alone. There are many different methods to differentiate the two; the most prevalent approach is using shared homology with a database of known genes. This method presents many pitfalls, most notably the catch that you only find genes that you have seen before. The four most popular prokaryotic gene-prediction programs (GeneMark, Glimmer, Prodigal, Phanotate) all use a protein-coding training model to predict protein-coding genes, with the latter three allowing for the training model to be created ab initio from the input genome. Different methods are available for creating the training model, and to increase the accuracy of such tools, we present here GOODORFS, a method for identifying protein-coding genes within a set of all possible open reading frames (ORFS). Our workflow begins with taking the amino acid frequencies of each ORF, calculating an entropy density profile (EDP), using KMeans to cluster the EDPs, and then selecting the cluster with the lowest variation as the coding ORFs. To test the efficacy of our method, we ran GOODORFS on 14,179 annotated phage genomes, and compared our results to the initial training-set creation step of four other similar methods (Glimmer, MED2, PHANOTATE, Prodigal). We found that GOODORFS was the most accurate (0.94) and had the best F1-score (0.85), while Glimmer had the highest precision (0.92) and PHANOTATE had the highest recall (0.96).

59 BASIC BIOLOGICAL SCIENCES↗

Ab initio simulations of metal contacts for graphene-based devices

The precise atomic structure of a metal contact significantly affects the performance of nanoscale electronic devices. We use an accurate, DFT-based non-equilibrium Green’s function method to evaluate various metal contacts with graphene or graphene nanoribbons. Furthermore, for surface metal contacts not chemically bound to graphene, Ti contacts have lower resistance than those of Au, Ca, Ir, Pt, and Sr. However, as an edge contact, Ti has larger resistance than Au. Bridging O atoms at Ti and Au edge contacts lowers the transmission by over 30%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessing the Age of an Asteroid's Surface with Data from the International Rosetta Mission

Rosetta is an international mission led by the European Space Agency (ESA) with key support and instrumentation from the National Aeronautics and Space Administration (NASA). Rosetta is currently on a ten-year mission to catch comet 67P/Churyumov-Gerasimenko (C-G); throughout its voyage, the spacecraft has performed flybys of two main belt asteroids (MBA): Steins and Lutetia. Data on the physical, chemical, and geological properties of these asteroids are currently being processed and analyzed. Accurate interpretation of such data is fundamental in the success of Rosetta's mission and overall objectives. Post-flyby data analyses strive to correlate the size, shape, volume, and rotational rate of Lutetia, in addition to interpreting its multi-color imagining, albedo, and spectral mapping. Although advancements in science have contributed to the examination of celestial bodies, methods to analyze asteroids remain largely empirical, not semi-empirical, nor ab initio. This study aims to interpret and document the scientific methods currently utilized in the characterization of asteroid (21) Lutetia in order to render these processes and methods accessible to the public. Examples include a standardized technique for assessing the age of an asteroid surface, complete with clickable reference maps, methodology of grouping surface characteristics together, and a standardized power law equation for the age. Other examples include determining the density of an object. Context for what both density and age mean is a bi-product of this study. Results of the study will aid in the development of pedagogical material on asteroids for public use, and in creation of an academic database for selected targets that might be used as a reference.

asteriods↗