SEARCH · Engineering Papers
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.
The goal of this lecture is to cover the basics of atomic electrons: Definitions of shells, subshells, and orbitals; Max number of electrons per subshell; Order of filling shells; The Periodic Table and electron configuration; and Covalent and Ionic Bonding.
Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.
Room-temperature realization of macroscopic quantum phases is one of the major pursuits in fundamental physics. The quantum spin Hall phase is a topological quantum phase that features a two-dimensional insulating bulk and a helical edge state. Here we use vector magnetic field and variable temperature based scanning tunnelling microscopy to provide micro-spectroscopic evidence for a room-temperature quantum spin Hall edge state on the surface of the higher-order topological insulator Bi 4 Br 4 . Here, we find that the atomically resolved lattice exhibits a large insulating gap of over 200 meV, and an atomically sharp monolayer step edge hosts an in-gap gapless state, suggesting topological bulk–boundary correspondence. An external magnetic field can gap the edge state, consistent with the time-reversal symmetry protection inherent in the underlying band topology. We further identify the geometrical hybridization of such edge states, which not only supports the Z 2 topology of the quantum spin Hall state but also visualizes the building blocks of the higher-order topological insulator phase. Our results further encourage the exploration of high-temperature transport quantization of the putative topological phase reported here.
The rate of heat transfer to the heat shield of a Jupiter probe has been estimated to be one order of magnitude higher than any previously experienced in an outer space exploration program. More than one-third of this heat load is due to an emission of continuum radiation from atoms and ions. The existing computer code for calculating the continuum contribution to the total load utilizes a modified version of Biberman's approximate method. The continuum radiation absorption cross sections of a C - H - O - N ablation system were examined in detail. The present computer code was evaluated and updated by being compared with available exact and approximate calculations and correlations of experimental data. A detailed calculation procedure, which can be applied to other atomic species, is presented. The approximate correlations can be made to agree with the available exact and experimental data.
Chemical ordering is a common phenomenon and highly correlated with the properties of solid materials. By means of the redistribution of atoms and chemical bonds, it invokes an effective lattice adjustment and tailors corresponding physical properties. To date, however, directly probing the 3D interfacial interactions of chemical ordering remains a big challenge. In this work, we deciphered the interlaced distribution of nanosized domains with chemical order/disorder in Fe 3 Pt bulk alloy. HAADF-STEM images evidence the existence of such nanodomains. The reverse Monte Carlo method with the X-ray pair distribution function data reveal the 3D distribution of local structures and the tensile effect in the disordered domains at the single-atomic level. The chemical bonding around the domain boundary changes the bonding feature in the disordered side and reduces the local magnetic moment of Fe atoms. This results in a suppressed negative thermal expansion and extended temperature range in Fe 3 Pt bulk alloy with nanodomains. Our study demonstrates a local revelation for the chemical order/disorder nanodomains in bulk alloy. The understanding gained from atomic short-range interactions within the domain boundaries provides useful insights with regard to designing new functional compounds.
In compositionally complex materials, there is controversy on the effect of enthalpy versus entropy on the structure and short-range ordering in so-called high-entropy materials. Here, to help address this controversy, we synthesized and probed 40 M 4 AlC 3 layered carbide phases containing two to nine metals and found that short-range ordering from enthalpy was present until the entropy increased enough to achieve complete disordering of the transition metals in their atomic planes. We transformed all of these layered carbide phases into two-dimensional (2D) sheets and showed the effects of the order versus disorder on their surface properties and electronic behavior. This study suggests the key effect that the competition between enthalpy and entropy has on short-range order in multicompositional materials.
Abstract Iron diantimonide is a material with the highest known thermoelectric power. By combining scanning transmission electron microscopic study with electronic transport neutron, X-ray scattering, and first principle calculation, we identify atomic defects that control colossal thermopower magnitude and nanoprecipitate clusters with Sb vacancy ordering, which induce additional phonon scattering and substantially reduce thermal conductivity. Defects are found to cause rather weak but important monoclinic distortion of the unit cell P n n m → P m . The absence of Sb along [010] for high defect concentration forms conducting path due to Fe d orbital overlap. The connection between atomic defect anisotropy and colossal thermopower in FeSb 2 paves the way for the understanding and tailoring of giant thermopower in related materials.
Revealing structure–activity relationships of graphene-supported atomically dispersed transition metal catalysts is of key importance in catalysis chemistry and materials science. In this work, we report a first-principles theoretical study on O 2 activation and CO oxidation catalyzed by atomically dispersed Pt 1 /Pt 2 anchored on a raised graphene (Gr) platform. The unique sharp-tip structure of the protuberant graphene platform (carbon-tip) can collect extra polarization electrons on the Pt-tip, promoting electron transfer between the catalyst and adsorbates, and greatly enhancing the chemical activity. Higher carbon-tips on defective graphene oxide with a carbon vacancy (Gr-V) produce more polarization electrons, a more localized electric field, and a larger up-shift of the d-orbital center of the Pt-tip, resulting in better catalytic performance. The carbon-tip enhancement reduces the energy barrier for the CO oxidation on Pt 1 O 2 /Gr-V from 1.19 eV to 0.56 eV compared with Pt 1 on planar Gr. In addition, CO poisoning is significantly alleviated, with the CO poisoning rate Γ θ reduced from 2.00 eV to 0.52 eV. Moreover, a linear relationship between the activation barrier and the binding energy of adsorbates is found for various atomically dispersed Pt-tip catalysts on the raised graphene oxide platform. Importantly, the order of catalytic activity is consistent with the carbon-tip enhancement, i.e., Pt 1 O 2 /Gr-V > Pt 2 O 4 /Gr-V > Pt 1 O 2 /Gr > Pt 2 O 4 /Gr. These findings provide important guidance for rational design of atomically dispersed catalysts.
A nuclear reaction depth profiling technique previously described by Filleux et al. (1977) has been used to measure the depth distribution of C on grain surfaces for Apollo 11, 15, 16 and 17 soil breccias. The surface C concentration of all samples studied lies between 2 and 8 times 10 to the 15th atoms per sq cm, showing no correlation with the volume C, which varies over an order of magnitude. If the observed variation represents the presence of unexposed grains on the surfaces studied, these results indicate a steady state surface C concentration of 5 to 10 times 10 to the 15th atoms per sq cm, accumulated over a time scale short compared with that required for the formation of volume-related C and with the mean lifetime of grains at the lunar surface. About one-third to one-half of the total C in lunar soil seems to be surface-correlated.
The transition probability values for a number of neutral atomic nitrogen (NI) lines in the visible wavelength range are determined in order to augment those given in the National Bureau of Standards Tables. These values are determined from experimentation as well as by using the published results of other investigators. The experimental determination of the lines in the 410 to 430 nm range was made from the observation of the emission from the arc column of an arc-heated wind tunnel. The transition probability values of these NI lines are determined to an accuracy of +/- 30% by comparison of their measured intensities with those of the atomic oxygen (OI) multiplet at around 615 nm. The temperature of the emitting medium is determined both using a multiple-layer model, based on a theoretical model of the flow in the arc column, and an empirical single-layer model. The results show that the two models lead to the same values of transition probabilities for the NI lines.
NASA's Atmospheric Tomography Mission (ATom) deployed in each of the four seasons during 2016-2018, the DC-8 aircraft in order to establish global-scale datasets intended to improve the representation of chemically reactive gases in global atmospheric chemistry models (ACMs). The Global Modeling Initiative (GMI) executed simulations for each ATom flight using the GMI Chemistry Transport Model (GMI-CTM) to provide species concentrations of chemical gases along the DC-8 flight transects. To solve the problem of translating the GMI-CTM simulation data to the unique spatial resolutions of each ATom flight, the GMI ICARTT Processing Software (GMI-IPS) was developed.The GMI-IPS is written in Python and provides data processing, flight extraction, and visualization support for aircraft research projects using ICARTT format, which is a standard format for airborne instrument data. Additionally, the GMI-IPS interpolates global gridded model data from Hierarchical Data Format (HDF) to ICARTT compatible flight transects. Software classes for instruments and collections provided by the ATom DC-8 aircraft such as MER10, MMS, etc. are derived from a common base class. Other functionality provided by the GMI-IPS are: deriving missing flight entries along a transect, reading ICARTT entries from file, and providing Python data structures for storing flight and model information, and more.The GMI-IPS is GIT source controlled, has approximately 30,000 lines of code, and supports parallelization across data collections. It delivered GMI-CTM data for more than forty distinct DC-8 aircraft flights that took place under ATom. The output ICARTT files adhere to format standard V1.1, and pass the scan utility provided by NASA LaRC Airborne Science Data for Atmospheric Composition. This presentation will include a software and methods overview, and results from ATom, including assessments using the GMI-CTM showing how well observations from ATom flight transects represent a broader region.
The crystal and magnetic structures of the nitride antiperovskite Mn3GeN reveals ferrimagnetic order stemming from a distorted kagome-derived lattice of the Mn atoms. Polycrystalline Mn3GeN was synthesized via a solid-state reaction and characterized using neutron powder diffraction, dc magnetometry, and first-principles calculations. Rietveld refinement reveals near-stoichiometric composition (Mn3GeN0.992(7)) adopting a tetragonal 𝐼4/𝑚𝑐𝑚 structure at 𝑇=500K and below, featuring axially distorted and tilted [NMn6] octahedra that result in a buckled Mn kagome lattice. On heating, the tetragonal distortion and octahedral tilt angle decrease continuously before transitioning to the cubic 𝑃𝑚3𝑚 antiperovskite phase at 𝑇≈527K. Neutron diffraction and magnetometry together reveal noncollinear ferrimagnetic ordering. For 30K≤𝑇≤500K, the magnetic structure is described by magnetic space group 𝐼𝑏𝑎′𝑚′ (72.544), with inequivalent Mn1 and Mn2 sublattices that couple antiferromagnetically to yield a net moment. Density-functional theory-based calculations show that the different local moments originate from the bandwidths associated with distinct Mn–N bond lengths. Temperature-dependent refinements reveal distinct differences in the thermal disordering profiles of the Mn1 and Mn2 sublattices. These findings reveal a subtlety in the magnetic and structural behavior of Mn3GeN, highlighting the interplay between structural distortions, magnetic ordering, and electronic structure in kagome-derived antiperovskite materials.
Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to leading order (–0.221 Z 5/3 and –0.021 Z ln Z respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E x (N,Z) ≈ –0.354 N 2/3 Z (as known before only for Z >> N >> 1) and E c ≈ –0.02 N ln N . These asymptotic estimates are most correct for atomic ions with large N and Z >> N, but we find that they are qualitatively and semi- quantitatively correct even for small N and for N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. Here, it is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.
CoAl and FeAl compounds are developed along two directions. Magnetic susceptibility and specific heat at low temperature on (NiCo)Al and (CoFe)Al ternary alloys are in good agreement with band calculations. Results on magnetization and specific heat under field at low temperature on nonstoichiometric compounds show clearly the importance of the nearest neighbor effects. In the case of CoAl, the isolated cobalt atoms substituting aluminum are characterized by a Kondo behavior, and, for FeAl, the isolated extra iron atoms are magnetic and polarize the matrix. Moreover, for the two compounds, clusters of higher order play a considerable part in the magnetic properties for CoAl, these clusters also seem to be characterized by a Kondo behavior, for FeAl, these clusters whose moment is higher than in the case of isolated atoms, could be constituted of excess parts of iron atoms.
The complex interplay between energetic and kinetic factors that governs the phase and morphology selections can originate at the earliest stage of crystallization in the amorphous parent phases. Because of the extreme difficulties in capturing the microscopic nucleation process, a detailed picture of how initial disordered structures affect the transformation pathway remains unclear. Here, we report the experimental observation of widely varying phase selection and grain size evolution during the devitrification of a homogeneous melt-spun glassy ribbon. Two different crystalline phases, θ-Al 5 Sm and ε-Al 60 Sm 11 , are found to form in the different regions of the same metallic glass ribbon during the devitrification. The grain size of ε-Al 60 Sm 11 phase shows a strong spatial heterogeneity. Coarse-grained ε-Al 60 Sm 11 phase coupled with the small volume fraction of θ-Al 5 Sm phase is preferably formed close to wheel side of the melt-spun ribbon. Combining experimental characterization and computational simulations, we show that phase selection and microstructure evolution can be traced back to different types and populations of atomic clusters that serve as precursors for the nucleation of different crystalline phases. Inhomogeneous cooling rates cause different structure orders across the glass sample during the quenching process. Our findings provide direct insight into the effect of structural order on the crystallization pathways during the devitrification of metallic glass. It also opens an avenue to study the detailed nucleation process at the atomic level using the metallic glass as a platform and suggests the opportunity of microstructure and property design via controlling the cooling process.
The reaction of hydrogen (H) atoms with butadiene (C4H2) was studied at room temperature in a pressure range between w mbar and 10 mbar. The primary step was an addition of H to C4H2 which is in its high pressure range at p 1 mbar. Under these conditions the following addition of a second H atom lies in the transition region between low and high pressure range. Vibrationally excited C4H4 can be deactivated to form buten-(1)-yne-(3)(C4H4) or decomposes into two C2H2 molecules. The rate constant at room temperature for primary step is given. The second order rate constant for the consumption of buten-(1)-yne-(3) is an H atom excess at room temperature is given.
In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.
In a paper with a similar title Yamanaka has calculated the mass polarization effect (to first order in mu/M for several low-lying states of the two-electron atoms and ions with atomic number Z from 2 to 10. Here we improve the previous results by using Hylleraas variational wave functions with up to 560 terms and extend the calculation to include some additional states and the Z=1ground state. In addition, we compute the second-order effect using the method of pseudostate summation. In one appendix another method of computation is discussed and used as a check, while the energies of the Z=1 ions are presented and discussed in a second appendix.