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.
Long range interaction of two H atoms calculated with electrostatic Hellmann-Feynman theorem, determining part of second order molecular wave function
Solid-liquid phase transitions, i.e. solidification processes, have applications in data storage, development of novel thermoelectric materials, cooling of microelectronic substrates and air conditioning condensers. Standard analyses of solidification (Stefan problem) assume constant thermal properties of the solid and liquid sides. It is not known how these properties change across the spatial transition interface, though most studies report a discontinuity in the solid and liquid properties through the transition temperature [1]. When the phase transition releases enthalpy, recent research has shown if the phonon or electron transport time is of the same order of magnitude as the time scale of the atomic transformation, this increases the heat capacity of the solid material at temperatures near the phase transition temperature [2]. A fundamental understanding of the phase transition may help shed light on the molecular origins of supercooling and spontaneous nucleation, which will help with applications involving them. We report studies of supercooling and nucleation of sodium sulfate decahydrate, a salt hydrate which is of recent interest in thermal storage, and of a hydrogel-sodium sulfate complex which shows limited supercooling. We will also report the measurement of their thermal properties during the phase transition process. This will be done with a hot-wire setup which produces small temperature changes of the order of ~1 C.
The 𝑅𝑀 6 𝑋 6 family of materials (𝑅 = rare-earth, 𝑀 = transition metal, 𝑋 = Ga, Si, Ge, Sn) produces an array of emergent phenomena, such as charge density waves, intrinsic Hall effects, and complex magnetic order, due to its kagome net of transition metal atoms, its local-moment magnetic anisotropies, and its extensive chemical tunability. Here, in this study, we report UNb 6 Sn 6 , a new “166” material containing both an actinide and a 4𝑑 transition metal, along with its nonmagnetic analog ThNb6Sn6, to investigate the properties of a 5𝑓−4𝑑 electron 166 system. UNb 6 Sn 6 crystallizes in the hexagonal 𝑃6/𝑚𝑚𝑚 space group with a small degree of disorder due to shifts in the size of the CoSn-like cages along the 𝑐 axis. Upon cooling at zero magnetic field, the material undergoes two magnetic phase transitions at 𝑇 2 = 46K and 𝑇 N = 43K. The low-temperature, zero-field phase is an antiferromagnet with ordered uranium moments and a k = (0,0,1/2) propagation vector determined by neutron diffraction. Remarkably, with a magnetic field applied along the 𝑐 axis, five additional magnetic transitions occur, evidenced by magnetization and resistivity data, before the moment saturates at 2.62 𝜇 B /U at 2 K and ≥ 13.6T. In two magnetic phase regions, the Hall resistivity of UNb 6 Sn 6 significantly deviates from the magnetization, suggesting that the phases have a large Berry curvature or a change in the Fermi surface. The unknown magnetic ordering of the field-dependent phases of UNb 6 Sn 6 demonstrates the complexity of the 5𝑓−4𝑑 166 system and encourages further study of its properties.
While many phenomena in cold atoms and other Bose–Einstein condensate (BEC) systems are often described using the mean-field approaches, understanding the kinetics of BECs requires the inclusion of particle scattering via the collision integral of the quantum Boltzmann equation. A rigorous approach for many problems in the dynamics of the BEC, such as the nucleation of the condensate or the decay of the persistent current, requires, in the presence of factors making a symmetry breaking possible, considering collisions with thermal atoms via the collision integral. These collisions permit the emergence of vorticity or other signatures of long-range order in the nucleation of the BEC or the transfer of angular momentum to thermal atoms in the decay of persistent current, due to corresponding terms in system Hamiltonians. Here, we also discuss the kinetics of spin–orbit-coupled BEC. The kinetic equation for the particle spin density matrix is derived. Numerical simulations demonstrate significant effects of the collision integral on the dynamics of the spin–orbit-coupled BEC upon quenching of the Raman coupling that generates synthetic electric and magnetic fields.
Here, the atomic structure of amorphous (FeCoNi) 70 Ti 10 B 20 alloy synthesized by mechanical alloying was investigated using high energy synchrotron X-ray diffraction and inverse Monte Carlo simulations of pair distribution functions. Empirical potential structure refinement indicates a chemical short-range order at the length scales of 2.1–2.5 Å via local atomic arrangements forming deformed bcc-like clusters. The structural model obtained was described by bond lengths, coordination numbers, and bond angle distribution functions determined for the first neighbor atoms by x-ray scattering supplemented with 3D Monte Carlo simulations.
Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.
We report the synthesis, characterization, and electrocatalytic water oxidation activity of two cobalt complexes, (6-FP)Co(NO 3 ) 2 (1) (6-FP = 8,8'-(1,2-phenylene)diquinoline) and (5-FP)Co(NO 3 ) 2 (2) (5-FP = 1,2-bis(N-7-azaindolyl)benzene), containing "capping arene" bidentate ligands with nitrogen atom donors. The cobalt complexes 1 and 2 were supported on ordered mesoporous carbon (OMC) by π-π stacking, resulting in heterogenized cobalt materials 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, respectively, and studied for electrocatalytic water oxidation. We find that 6-FP-Co-OMC-1 exhibits an overpotential of 355 mV for a current density of 10 mA cm -2 and a turnover frequency (TOF) of ~0.53 s -1 at an overpotential of 400 mV at pH 14. 6-FP-Co-OMC-1 exhibits activity that is ~1.6 times that of 5-FP-Co-OMC-2, which gives a TOF of 0.32 s -1 at 400 mV overpotential. The structural stability of the single-atom Co site was demonstrated for 6-FP-Co-OMC-1 using X-ray absorption spectroscopy for the molecular complex supported on OMC, but slow degradation in catalyst activity can be attributed to eventual formation of Co oxide clusters. DFT computations of electrocatalytic water oxidation using the molecular complexes as models provide a description of the catalytic mechanism. These studies reveal that the mechanism for O-O bond formation involves an intermediate Co IV oxo complex that undergoes an intramolecular reductive O-O coupling to form a Co II -OOH species. Further, the calculations predict that the molecular 6-FP-Co structure is more active for electrocatalytic water oxidation than 5-FP-Co, which is consistent with experimental studies of 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, highlighting the possibility that the ligand structure influences the catalytic activity of the supported molecular catalysts.
We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.
Aberration correction is an important aspect of modern high-resolution scanning transmission electron microscopy. Most methods of aligning aberration correctors require specialized sample regions and are unsuitable for fine-tuning aberrations without interrupting on-going experiments. Here, we present an automated method of correcting first- and second-order aberrations called BEACON, which uses Bayesian optimization of the normalized image variance to efficiently determine the optimal corrector settings. We demonstrate its use on gold nanoparticles and a hafnium dioxide thin film showing its versatility in nano- and atomic-scale experiments. BEACON can correct all first- and second-order aberrations simultaneously to achieve an initial alignment and first- and second-order aberrations independently for fine alignment. Ptychographic reconstructions are used to demonstrate an improvement in probe shape and a reduction in the target aberration.
The proposed research will investigate the stability and cryogenic properties of solid propellants that are critical to NASA s goal of realizing practical propellant designs for future spacecraft. We will determine the stability and thermal properties of a solid hydrogen-liquid helium stabilizer in a laboratory environment in order to design a practical propellant. In particular, we will explore methods of embedding atomic species and metallic nano-particulates in hydrogen matrices suspended in liquid helium. We will also measure the characteristic lifetimes and diffusion of atomic species in these candidate cryofuels. The most promising large-scale advance in rocket propulsion is the use of atomic propellants; most notably atomic hydrogen stabilized in cryogenic environments, and metallized-gelled liquid hydrogen (MGH) or densified gelled hydrogen (DGH). The new propellants offer very significant improvements over classic liquid oxygen/hydrogen fuels because of two factors: (1) the high energy-release, and (ii) the density increase per unit energy release. These two changes can lead to significant reduced mission costs and increased payload to orbit weight ratios. An achievable 5 to 10 percent improvement in specific impulse for the atomic propellants or MGH fuels can result in a doubling or tripling of system payloads. The high-energy atomic propellants must be stored in a stabilizing medium such as solid hydrogen to inhibit or delay their recombination into molecules. The goal of the proposed research is to determine the stability and thermal properties of the solid hydrogen-liquid helium stabilizer. Magnetic resonance techniques will be used to measure the thermal lifetimes and the diffusive motions of atomic species stored in solid hydrogen grains. The properties of metallic nano-particulates embedded in hydrogen matrices will also be studied and analyzed. Dynamic polarization techniques will be developed to enhance signal/noise ratios in order to be able to detect low concentrations of the introduced species. The required lifetimes for atomic hydrogen and other species can only be realized at low temperatures to avoid recombination of atoms before use as a fuel.
As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.
In this work we demonstrate a newly developed capability to capture pitting of carbon fibers in DSMC simulations, specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. State-of-the-art reactive surface models in DSMC compute collision dependent carbon consumption rates (usually through desorption of CO) based on a set of surface reactions that has been derived from molecular beam experiments. The reactivity on each carbon surface element is constant in those models, such that the carbon surface recedes uniformly as a result of ablation. However, it is well known that in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. These defective sites have a much higher reactivity than the average sites (2-3 orders of magnitude) and are first to react during ablation leading to its removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these ”active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. Within this framework, we have implemented the capability of a single surface having multiple site sets with different reactivities. Using this feature, we can simulate the presence of active sites on carbon surfaces, whose reactivity is much greater than an average site as a result of defects. We have implemented the active site fraction as a property of surface elements within SPARTA, which is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface, and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we are able to capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.
The parent compounds of iron-based superconductors are magnetically ordered bad metals, with superconductivity appearing near a putative magnetic quantum critical point. The presence of both Hubbard repulsion and Hund’s coupling leads to rich physics in these multiorbital systems, and motivated descriptions of magnetism in terms of itinerant electrons or localized spins. The NaFe 1–x Cu x As series consists of magnetically ordered bad metal (x = 0), superconducting (x ≈ 0.02) and magnetically ordered semiconducing/insulating (x ≈ 0.5) phases, providing a platform to investigate the connection between superconductivity, magnetism and electronic correlations. Here we use x-ray absorption spectroscopy and resonant inelastic x-ray scattering to study the valence state of Fe and spin dynamics in two NaFe 1–x Cu x As compounds (x = 0 and 0.47). We find that magnetism in both compounds arises from Fe 2+ atoms, and exhibits underdamped dispersive spin waves in their respective ordered states. The dispersion of spin excitations in NaFe 0.53 Cu 0.47 As is consistent with being quasi-one-dimensional. Compared to NaFeAs, the band top of spin waves in NaFe 0.53 Cu 0.47 As is slightly softened with significantly more spectral weight of the spin excitations. Here, our results indicate the spin dynamics in NaFe 0.53 Cu 0.47 As arise from localized magnetic moments and suggest the iron-based superconductors are proximate to a correlated insulating state with localized iron moments.
We report on the experimental observation of non-resonant, second-order optical sum-frequency generation (SFG) in five different atomic and molecular gases. The measured signal is attributed to a SFG process by characterizing its intensity scaling and its polarization behavior. We show that the electric quadrupole mechanism cannot explain the observed trends and suggest a mechanism based on symmetry breaking along the incident beam path arising from laser-induced species ground state number density gradients. Our results demonstrate that the SFG is about four orders of magnitude stronger than the third-harmonic generation (THG) and independent from any externally applied electric fields. In conclusion, these features make this method suitable for gas number density measurements at the picosecond time scale in reactive flows and plasmas.
The electrochemical nitrate reduction reaction (NO 3 RR) to ammonia is an essential step toward restoring the globally disrupted nitrogen cycle. In search of highly efficient electrocatalysts, tailoring catalytic sites with ligand and strain effects in random alloys is a common approach but remains limited due to the ubiquitous energy-scaling relations. With interpretable machine learning, we unravel a mechanism of breaking adsorption-energy scaling relations through the site-specific Pauli repulsion interactions of the metal d-states with adsorbate frontier orbitals. The non-scaling behavior can be realized on (100)-type sites of ordered B2 intermetallics, in which the orbital overlap between the hollow *N and subsurface metal atoms is significant while the bridge-bidentate *NO 3 is not directly affected. Among those intermetallics predicted, we synthesize monodisperse ordered B2 CuPd nanocubes that demonstrate high performance for NO 3 RR to ammonia with a Faradaic efficiency of 92.5% at –0.5 V RHE and a yield rate of 6.25 mol h –1 g –1 at –0.6 V RHE . This study provides machine-learned design rules besides the d-band center metrics, paving the path toward data-driven discovery of catalytic materials beyond linear scaling limitations.
Neutral atomic oxygen is the most abundant component of the ionospheric plasma in the low Earth orbit environment (LEO; 200 to 700 kilometers altitude) and can produce significant degradation of some spacecraft materials. In order to produce a more complete understanding of the materials chemistry of atomic oxygen, the chemistry and physics of O-atom interactions with materials were determined in three radically different environments: (1) The Space Shuttle cargo bay in low Earth orbit (the EOIM-3 space flight experiment), (2) a high-velocity neutral atom beam system (HVAB) at Los Alamos National Laboratory (LANL), and (3) a microwave-plasma flowing-discharge system at JSC. The Space Shuttle and the high velocity atom beam systems produce atom-surface collision energies ranging from 0.1 to 7 eV (hyperthermal atoms) under high-vacuum conditions, while the flowing discharge system produces a 0.065 eV surface collision energy at a total pressure of 2 Torr. Data obtained in the three different O-atom environments referred to above show that the rate of O-atom reaction with polymeric materials is strongly dependent on atom kinetic energy, obeying a reactive scattering law which suggests that atom kinetic energy is directly available for overcoming activation barriers in the reaction. General relationships between polymer reactivity with O atoms and polymer composition and molecular structure have been determined. In addition, vacuum ultraviolet photochemical effects have been shown to dominate the reaction of O atoms with fluorocarbon polymers. Finally, studies of the materials chemistry of O atoms have produced results which may be of interest to technologists outside the aerospace industry. Atomic oxygen 'spin-off' or 'dual use' technologies in the areas of anisotropic etching in microelectronic materials and device processing, as well as surface chemistry engineering of porous solid materials are described.
Recent studies on Ba 2 CoO 4 (BCO) and SrRuO 3 (SRO) have unveiled a variety of intriguing phenomena, such as magnetic polarization on oxygens, unexpected nodes in the spin density profile along bonds, and unusual zigzag spin patterns in triangular lattices. Here, using simple model calculations supplemented by density functional theory we explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the transition metal (TM) to oxygen bonds. Our simple model also allows us to explain the net polarization on oxygen as originated from the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Overall, our results are not limited to BCO and SRO, but they are generic and qualitatively predict the net polarization expected on any ligands, according to the spin order of the closest TM atoms and the number of intermediate ligand atoms. Finally, we propose that a robust easy-axis anisotropy would suppress the competing 120° antiferromagnetic order to stabilize the zigzag pattern order as ground state in a triangular lattice. Our generic predictions should be applicable to any other compound with characteristics similar to those of BCO and SRO.
The kinetics of the Cl + Cl2O reaction, a possible source of ClO(2 Pi) radicals for atmospheric photochemical studies, are investigated at 298 K. The discharge flow/mass spectrometry and discharge flow/resonance fluorescence techniques were used to monitor the decay of C12O in the presence of excess concentrations of atomic chlorine and chlorine monoxide, respectively. The pseudo-first order rate constants obtained from both experiments are found to be in excellent agreement, averaging 9.8 + or - 0.8 x 10 to the -11th cu cm/molecule per sec. Results are consistent with the lower limit obtained by Edgecombe et al. (1957) but differ by a factor of 150 from those of Basco and Dogra (1971). The present value is also noted to be consistent with a lower value for the rate constant of the reaction of oxygen atoms with Cl2O.