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.
This erratum corrects a factor 2 error in previous work on current-dependent TDDFT.
Spin-1/2 chains with alternating antiferromagnetic (AFM) and ferromagnetic (FM) couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density-matrix renormalization-group (DMRG) methods, we have systematically studied the dimerized chain system Na 2 Cu 2 TeO 6 with a d 9 electronic configuration. Near the Fermi level, in the nonmagnetic phase the dominant states are mainly contributed by the Cu 3d x 2 -y 2 orbitals highly hybridized with the O 2p orbitals, leading to an “effective” single-orbital low-energy model. By calculating the relevant hoping amplitudes, we explain the size and sign of the exchange interactions in Na 2 Cu 2 TeO 6 . In addition, a single-orbital Hubbard model is constructed for this dimerized chain system where the quantum fluctuations are taken into account. Both AFM and FM couplings (leading to an ↑|-|↓-|↓-↑| state) along the chain were found in our DMRG and Lanczos calculations, in agreement with density functional theory and neutron-scattering results. The hole pairing binding energy ΔE is predicted to be negative at Hubbard U~ 11eV, suggesting incipient pairing tendencies.
This paper introduces a novel track-length extension fitting algorithm for measuring the kinetic energies of inelastically interacting particles in liquid argon time projection chambers (LArTPCs). The algorithm finds the most probable offset in track length for a track-like object by comparing the measured ionization density as a function of position with a theoretical prediction of the energy loss as a function of the energy, including models of electron recombination and detector response. The algorithm can be used to measure the energies of particles that interact before they stop, such as charged pions that are absorbed by argon nuclei. The algorithm's energy measurement resolutions and fractional biases are presented as functions of particle kinetic energy and number of track hits using samples of stopping secondary charged pions in data collected by the ProtoDUNE-SP detector, and also in a detailed simulation. Additional studies describe the impact of the dE/dx model on energy measurement performance. The method described in this paper to characterize the energy measurement performance can be repeated in any LArTPC experiment using stopping secondary charged pions.
Wurtzite LaN (wz-LaN) is a semiconducting nitride that has piezoelectric and ferroelectric properties, making it promising for applications in electronics, either as a binary compound or in alloys such as LaAlN. The prospects for wz-LaN in devices are influenced by the properties of point defects and impurities; here, we use first-principles density functional theory with a hybrid functional to calculate their formation energies, as well as their atomic and electronic structures. Among native point defects, we find that nitrogen-related defects, both vacancies ($V^+_N$) and interstitials ($N^-_i$), are energetically most favorable under most relevant chemical potentials and positions of the Fermi level; $V^0_N$ may additionally be observed under N-poor conditions, and $N^0_i$ may be prominent under N-rich conditions. We also investigate the incorporation of oxygen and hydrogen, which will likely be present as unintentional impurities. We find that the $O^+_N$ substitutional species readily forms, but oxygen will not lead to n-type conductivity due to formation of DX centers and compensation by interstitial defects. Similarly, substitutional HN and interstitial H i can compensate both p- and n-type dopants. Our results provide detailed, microscopic guidance for the development of electronic devices based on wz-LaN.
Grading of bandgap by alloying CdTe with selenium to form a CdSe x Te 1–x /CdTe-graded bilayer device has led to a device efficiency over 19%. A CdSe x Te 1–x absorber would increase the short-circuit current due to its lower bandgap but at the expense of open-circuit voltage. It has been demonstrated that adding a CdTe layer at the back of such a CdSe x Te 1–x film reduces the voltage deficit caused by the lower bandgap of absorber from selenium alloying while maintaining the higher short-circuit current. This leads to a photovoltaic device that draws advantage from both materials with an efficiency greater than either of them. In this report there is a detailed account using device data, ultraviolet photoelectron spectroscopy, electron microscopy, and first-principles density functional theory modeling is provided, which shows that CdTe acts as an electron reflector for CdSe x Te 1–x .
Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.
Surface wettability is determined by the attraction of a liquid phase to a solid surface. It is typically quantified by using contact angle measurements at mineral surfaces in the case of the flotation of mineral particles. Contact angle research to describe wettability has been investigated at different scales by sessile drop measurements, molecular dynamic simulation, and atomic force microscopy. In this study, the density functional theory (DFT) was employed for predicting the surface free energy and contact angles of a well-known hydrophobic phyllosilicate mineral talc and a well-known hydrophilic phyllosilicate mineral muscovite based on the calculated interfacial energy and surface charge. The results revealed that the predicted contact angle at the atomic scale was larger than the experimental value, and identified two interactions: electrostatic interaction and hydrogen bonding, between the hydrophilic muscovite surface and the water layer, while a water-exclusion zone of 3.346 Å was found between the hydrophobic talc surface and the first water layer. This investigation gives a new perspective for wettability determination at the atomic scale.
Potassium (K) is an essential nutrient for plant growth, and despite its abundance in soil, most of the K is structurally bound in minerals, limiting its bioavailability and making this soil K reservoir largely inaccessible to plants. Microbial biochemical weathering has been shown to be a promising pathway to sustainably increase plant available K. However, the mechanisms underpinning microbial K uptake, transformation, storage, and sharing are poorly resolved. Here, to better understand the controls on microbial K transformations, we performed K K-edge x-ray absorption near-edge structure (XANES) spectroscopy on K-organic salts, including acetate, citrate, nitrate, oxalate, and tartrate, which are frequently observed as low molecular weight organic acids secreted by soil microbes, as well as humic acid, which acts as a proxy for higher molecular weight organic acids. The organic salts display feature-rich K XANES spectra, each demonstrating numerous unique features spanning ~13 eV range across the absorption edge. In contrast, the spectra for humic acid have one broad, wide feature across the same energy range. We used a combination of time-dependent density functional theory and the Bethe–Salpeter equation based approach within the OCEAN code to simulate the experimental spectra for K-nitrate (KNO 3 ) and K-citrate [K 3 (C 6 H 5 O 7 )·H 2 O] to identify the electronic transitions that give rise to some of the outlying and unique spectral features in the organic salts. KNO 3 has both the lowest and highest lying energy features, and K 3 (C 6 H 5 O 7 )·H 2 O is produced by several soil microbes and is effective at mineral weathering. Our results analyze the K-organic salt bonding in detail to elucidate why the spectral shapes differ and indicate that the K K-edge XANES spectra are associated with the entire ligand despite similar first-shell bonding environments around the K center. The improved understanding of K bonding environments with organic ligands and their use for interpretation of the K-XANES spectra provides an important toolkit to understand how K is transformed by microbial processes and made bioavailable for plant uptake.
Ammonia affords dense storage for renewable energy as a fungible liquid fuel, provided it can be efficiently synthesized from hydrogen and nitrogen. In this work, the catalysis of ammonia synthesis was computationally explored beyond the Sabatier limit by dynamically straining a ruthenium crystal (±4%) at the resonant frequencies (10 2 to 10 5+ Hz) of N 2 surface dissociation and hydrogenation. Density functional theory calculations at different strain conditions indicated that the energies of NH x surface intermediates and transition states scale linearly, allowing the description of ammonia synthesis at a continuum of strain conditions. A microkinetic model including multiple sites and surface diffusion between step and Ru(0001) terrace sites of varying ratios for nanoparticles of differing size revealed that dynamic strain yields catalytic ammonia synthesis conversion and turnover frequency comparable to industrial reactors (400°C, 200 atm) but at lower temperature (320°C) and an order of magnitude lower pressure (20 atm).
We detail the estimation of activation energies and quantum nuclear vibrational tunneling effects for hydrogen diffusion in PuO 2 based on Density Functional Theory calculations and a quantum double well approximation. We find that results are relatively insensitive to choice of exchange correlation functional. In addition, the representation of spin in the system and use of an extended Hubbard U correction has only a small effect on hydrogen point defect formation energies when the PuO 2 lattice is held fixed at the experimental density. We then compute approximate activation energies for transitions between hydrogen interstitial sites seeded by a semi-empirical quantum model and determine the quantum tunneling enhancement relative to classical kinetic rates. Our model indicates that diffusion rates in H/PuO 2 systems could be enhanced by more than one order of magnitude at ambient conditions and that these effects persist at high temperature. The method we propose here can be used as a fast screening tool for assessing possible quantum nuclear vibrational effects in any number of condensed phase materials and surfaces, where hydrogen hopping tends to follow well defined minimum energy pathways.
Surface energies of metal-based systems are important for determining the Wulff-constructed shapes of metal nanoparticles and understanding the stability. We have developed a physics informed method to predict the total energy of metal-based systems across a wide range of configurations. Our method has been tested against density functional theory (DFT) calculations for late transition metals. This method enables on-the-fly surface energy predictions based on the fundamental of materials property and allows for the Wulff construction of metal particles for a random number of elemental atoms and without the need for DFT calculations. By making a division between atoms in the different layers of the model system, we can improve the accuracy of the model, suggesting a dissimilarity between the electronic structure due to an alternating compression and expansion of atomic layers. We find that our model accurately and effectively provides valuable insights into the distribution and stability of nanoparticles.
In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal–organic frameworks (MOFs) known as Cu–BTT (BTT–3 = 1,3,5-benzenetristetrazolate) and Cu–BTTri (BTTri–3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu–BTTri compared to its tetrazole counterpart, Cu–BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu–BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu–BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn–Teller effect in Cu–BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials’ potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu–BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.
Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.
Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.
This paper reports the effect of mechanically applied elastic strain on the hydrogen evolution reaction (HER) activity of graphene under acidic conditions. Here, applied tensile strain of 0.2% on a graphene electrode is shown to lead to a 1–3% increase of HER current. Tensile strain increases HER activity while compressive strain decreases it. Density functional theory (DFT) calculations using a periodic graphene slab model show increasing in the adsorption energy of H atom with growing tensile strain, leading to a corresponding enhancement of current density in HER, similar to that observed experimentally.
Equations of state for molecular and metallic hydrogen are calculated in an accurate and comparable manner by using the density-functional method applied to static lattices. By isolating the electronic energies in this way, attention can be drawn to the relative importance of the protonic degrees of freedom. From the equation-of-state results it can be concluded that a remnant molecular pairing is preferred in a band-overlap metallic state, complete dissociation occurring only at very high densities. An accurate determination of the pressure needed to achieve complete dissociation is shown to require a self-consistent treatment of electron and proton degrees of freedom.
An overarching challenge of the electrochemical carbon dioxide reduction reaction (eCO 2 RR) is finding an earth-abundant, highly active catalyst that selectively produces hydrocarbons at relatively low overpotentials. Here, we report the eCO 2 RR performance of two-dimensional transition metal carbide class of materials. Our results indicate a maximum methane (CH 4 ) current density of -421.63 mA/cm 2 and a CH 4 faradic efficiency of 82.7% ± 2% for di-tungsten carbide (W 2 C) nanoflakes in a hybrid electrolyte of 3 M potassium hydroxide and 2 M choline-chloride. Powered by a triple junction photovoltaic cell, we demonstrate a flow electrolyzer that uses humidified CO 2 to produce CH 4 in a 700-h process under one sun illumination with a CO 2 RR energy efficiency of about 62.3% and a solar-to-fuel efficiency of 20.7%. Density functional theory calculations reveal that dissociation of water, chemisorption of CO 2 and cleavage of the C-O bond—the most energy consuming elementary steps in other catalysts such as copper—become nearly spontaneous at the W2C surface. This results in instantaneous formation of adsorbed CO—an important reaction intermediate—and an unlimited source of protons near the tungsten surface sites that are the main reasons for the observed superior activity, selectivity, and small potential.
In this work, the heat capacities of nanocrystalline Ni 3 Fe and control materials with larger crystallites were measured from 0.4–300 K. The heat capacities were integrated to obtain the enthalpy, entropy, and Gibbs free energy and to quantify how these thermodynamic functions are altered by nanocrystallinity. From the phonon density of states (DOS) measured by inelastic neutron scattering, we find that the Gibbs free energy is dominated by phonons and that the larger heat capacity of the nanomaterial below 100 K is attributable to its enhanced phonon DOS at low energies. Besides electronic and magnetic contributions, the nanocrystalline material has an additional contribution at higher temperatures, consistent with phonon anharmonicity. The nanocrystalline material shows a stronger increase with temperature of both the enthalpy and entropy compared to the bulk sample. Its entropy exceeds that of the bulk material by 0.4 $k_B$/atom at 300 K. This is insufficient to overcome the enthalpy of grain boundaries and defects in the nanocrystalline material, making it thermodynamically unstable with respect to the bulk control material.