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.
Review of some results pertinent to the development of nonlinear theory on two-way coherent tracking systems. In particular, the model of cascaded systems is presented, and approximations to steady state probability density functions of the two-way system phase and Doppler error are developed. From these, certain numerical results required in the design and planning of these systems are derived.
The presence of point defects, such as vacancies, plays an important role in materials design. Here, we explore the extrapolative power of a graph neural network (GNN) to predict vacancy formation energies. We show that a model trained only on perfect materials can also be used to predict vacancy formation energies (E vac ) of defect structures without the need for additional training data. Such GNN-based predictions are considerably faster than density functional theory (DFT) calculations and show potential as a quick pre-screening tool for defect systems. To test this strategy, we developed a DFT dataset of 530 E vac consisting of 3D elemental solids, alloys, oxides, semiconductors, and 2D monolayer materials. We analyzed and discussed the applicability of such direct and fast predictions. We applied the model to predict 192 494 E vac for 55 723 materials in the JARVIS-DFT database. Our work demonstrates how a GNN-model performs on unseen data.
A fully self-consistent density-functional theory (DFT) with improved functionals is used herein to provide a comprehensive account of structural, electronic, and optical properties of C 3 N 4 polymorphs. Using our recently developed van Leeuwen-Baerends (vLB) corrected local-density approximation (LDA), we implemented LDA + vLB within full-potential N th -order muffin-tin orbital (FP-NMTO) method and show that it improves structural properties and band gaps compared to semi-local functionals (LDA/GGA). We demonstrate that the LDA + vLB predicts band-structure and work-function for well-studied 2D-graphene and bulk-Si in very good agreement with experiments, and more exact hybrid functional (HSE) calculations as implemented in the Quantum-Espresso (QE) package. The structural and electronic-structure (band gap) properties of C 3 N 4 polymorphs calculated using FP-NMTO-LDA + vLB is compared with more sophisticated hybrid-functional calculations. We also perform detailed investigation of photocatalytic behavior using QE-HSE method of C 3 N 4 polymorphs through work-function, band (valence and conduction) position with respect to water reduction and oxidation potential. Our results show $γ$-C 3 N 4 as the best candidate for photocatalysis among all the C 3 N 4 polymorphs but it is dynamically unstable at ‘zero’ pressure. We show that $γ$-C 3 N 4 can be stabilized under hydrostatic-pressure, which improves its photocatalytic behavior relative to water reduction and oxidation potentials.
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.
Here, we have studied the surface structure of a single crystal β-Ga 2 O 3 (010) using quantitative Low Energy Electron Diffraction (LEED) and X-ray photoelectron spectroscopy (XPS). The XPS measurements show spectra typical of stoichiometric Ga 2 O 3 with a clean surface. LEED consistently shows a p(1x1) pattern, free of surface reconstruction. Quantitative LEED I(V) curves are acquired for 41 distinct diffraction spots. The experimental I(V) curves are compared to simulations over the first five layers. The best fits to the experimental LEED I(V) curves acquired at all diffraction spots are then used to calculate the interplanar relaxation and atomic rumpling. Significant atomic rumpling and interplanar relaxation are found over the first 5 atomic layers. As a result of rumpling a polarization of ~ 2 μC/cm 2 develops in the topmost surface layer. The structural results are in good agreement with previous density functional theory calculations and experimental X-ray photoelectron diffraction.
Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode–electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. Here we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.
Ternary intermetallic compound \({\text {Pr}}_2 {\text {Co}}_{0.86} {\text {Si}}_{2.88}\) has been synthesized in single phase and characterized by x-ray diffraction, scanning electron microscopy with energy dispersive x-ray spectroscopy (SEM-EDX) analysis, magnetization, heat capacity, neutron diffraction and muon spin rotation/relaxation ( \(\mu\) SR) measurements. The polycrystalline compound was synthesized in single phase by introducing necessary vacancies in Co/Si sites. Magnetic, heat capacity, and zero-field neutron diffraction studies reveal that the system undergoes magnetic transition below \(\sim\) 4 K. Neutron diffraction measurement further reveals that the magnetic ordering is antiferromagnetic in nature with an weak ordered moment. The high temperature magnetic phase has been attributed to glassy in nature consisting of ferromagnetic clusters of itinerant (3 d ) Co moments as evident by the development of internal field in zero-field \(\mu\) SR below 50 K. The density-functional theory (DFT) calculations suggest that the low temperature magnetic transition is associated with antiferromagnetic coupling between Pr 4 f and Co 3 d spins. Pr moments show spin fluctuation along with unconventional orbital moment quenching due to crystal field. The evolution of the symmetry and the crystalline electric field environment of Pr-ions are also studied and compared theoretically between the elemental Pr and when it is coupled with other elements such as Co. The localized moment of Pr 4 f and itinerant moment of Co 3 d compete with each other below \(\sim\) 20 K resulting in an unusual temperature dependence of magnetic coercivity in the system.
Molecular pair distribution function in perfect crystalline solids, using method based on cell cluster theory
Developing active and durable air electrodes for efficient oxygen reactions is challenging for protonic ceramic cells (PCCs), especially at temperatures below 550°C. Here, in this study, we report a rationally designed conformal coating with a high-entropy PrNi 0.2 Mn 0.2 Co 0.2 Fe 0.2 Cu 0.2 O 3−δ (PNMCFC) perovskite structure on the surface of a state-of-the-art PrBaCo 2 O 5+δ (PBC) air electrode. The formed hybrid air electrode (PNMCFC-PBC) shows faster surface oxygen kinetics and a more stable phase structure in high-humidity air than the bare PBC electrode. Further density functional theory calculations suggest that the conformal coating mitigates Ba segregation at the interface and improves oxygen-related reactions, enhancing overall stability and electrocatalytic performance. The cells with the developed hybrid electrodes show encouraging electrochemical performance at 550°C: a polarization resistance of 0.72 Ω cm 2 , a peak power density of 1.30 W cm −2 , an electrolysis current density of −1.36 A cm −2 at 1.3 V, and reasonable operating stabilities (∼200 h at 550°C).
Among the most significant developments in heterogeneous catalysis in the last 20 years is the emergence of microkinetic models, often parameterized from density functional theory (DFT) calculations, used to quantify observed catalyst performance and to guide the discovery of new catalytic materials. While DFT directly reports binding energies and elementary step activation energies, the free energies that enter into microkinetic models must be computed from additional approximations. Frequently these free energy approximations assume ideal behavior—that adsorbates do not interact with one another, or that adsorbates are immobile or vibrate harmonically about a binding site. These assumptions can and do have an impact on predicted catalyst performance and potentially even on predicted trends. Our work relates to three categories of non-idealities, explored in the context of nitrogen oxidation and reduction catalysis on metal surfaces. One component relates to adsorbate translational free energy, the contribution least well described by conventional models, in which we develop modeling approaches that improve accuracy with limited increase in computational expense. A second component relates adsorbate-adsorbate interactions, in which we develop benchmark on-lattice interaction models and compare kinetic predictions with conventional mean-field models, in an effort to develop more robust coverage-dependent mean field modeling approaches. The last relates the most fundamental assumption of all—that of energy equipartition—in an effort to rationalize and guide plasma-enhanced catalysis.
Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.
In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.
The development of next-generation quaternary chalcogenides, such as Cu2ZnSnS4 (CZTS) and Cu2ZnGeS4 (CZGS), for solar energy and thermoelectric applications hinges upon both careful experimentation and accurate quantum mechanical modeling. To address the latter, many have turned to density functional theory (DFT), which offers several choices for the approximate treatment of electron exchange and correlation (XC). Popular XC functionals include the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) and the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. Extensions of DFT functionals, such as adding a Hubbard U correction and introducing a fraction of the Fock exchange (hybrid functionals), have been used widely to model systems containing 3d metal ions. However, no studies yet have compared comprehensively PBE(+U) and SCAN(+U) in the quality of their predictions of the bulk and defect thermodynamics of quaternary chalcogenides, which play a critical role in device fabrication and performance. Hence, here we calculate the (i) 0 K formation energies of bulk Ge compounds and (ii) neutral defect formation energies including charge-balanced (e.g., CuZn + ZnCu) and charge-imbalanced (e.g., CuSn) combinations of antisites and vacancies in CZTS and CZGS using the PBE, PBE +U, SCAN, SCAN +U, and the hybrid Heyd-Scuseria-Ernzerhof XC frameworks. We find that the formation energies of charge-imbalanced defects are more sensitive to the choice of the XC functional than those of charge-balanced defects, which can be explained by the differences in the extent of penalization of defect-generated delocalized electrons/holes by PBE, PBE +U, SCAN, and SCAN +U. Additionally, our results show that SCAN systematically underbinds Ge-containing compounds, thus highlighting the need for even further improvement of XC functionals. Based on our findings, we recommend the use of SCAN for modeling quaternary chalcogenides because its errors are systematic, and it has the firmest theoretical underpinning. Our work provides guidance for future modeling of quaternary chalcogenides.
Chemical looping with oxygen uncoupling (CLOU) materials is actively sought for combustion of carbonaceous materials to achieve complete conversion and capture of carbon dioxide. These materials may play a vital role in reducing atmospheric carbon via negative carbon output. However, there is no one‐size‐fits‐all approach as different operating conditions and feedstocks may require different CLOU materials. As a result, the exploration and discovery of high‐performance CLOU materials can be a slow process. To address this challenge, a high‐throughput inverse machine learning workflow that identifies optimum materials from perovskite oxides for a given set of targets is developed—temperature and Gibbs free energy of oxygen formation. The model is trained on high‐throughput density functional theory calculations of CLOU materials and inverts the materials design process using a genetic algorithm to produce realistic substituted SrFeO 3‐δ compositions as output. Using the inverse model, it is able to identify several interesting new families of CLOU materials: Sr 1‐ x A x Fe 1‐ y B y O 3‐δ (e.g., A = Ca or K; B = Mg, Bi, Mn, Ni, Co, Cu, or Zn). These materials have shown promising properties, and some of them even outperform the benchmark material in terms of oxygen release kinetics under relevant CLOU operating conditions.
On the one hand, multi-principal element alloys (MPEAs) have created a paradigm shift in alloy design due to large compositional space, whereas on the other, they have presented enormous computational challenges for theory-based materials design, especially density functional theory (DFT), which is inherently computationally expensive even for traditional dilute alloys. In this project, we developed a machine learning framework, namely PREDICT ( PR edict properties from E xisting D atabase I n C omplex alloys T erritory), that opens a pathway to predict elastic constants in large compositional space with little computational expense. The framework only relies on the DFT database of binary alloys and predicts Voigt–Reuss–Hill Young’s modulus, shear modulus, bulk modulus, elastic constants, and Poisson’s ratio in MPEAs. We show that the key descriptors of elastic constants are the A–B bond length and cohesive energy. The framework can predict elastic constants in hypothetical compositions as long as the constituent elements are present in the database, thereby enabling property exploration in multi-compositional systems. We illustrate predictions in a FCC Ni-Cu-Au-Pd-Pt system.
Predicting crystal structure has been a challenging problem in physics and materials science for a long time. A reliable energy calculation engine combined with an efficient global search algorithm, such as particle swarm optimization algorithm or genetic algorithm, is needed to conduct crystal structure prediction. In recent years, machine learning-based interatomic potential energy surface models have been proposed, potentially allowing us to perform crystal structure prediction for systems with the accuracy of density functional theory (DFT) and the speed of empirical force fields. In this paper, we employ a previously developed Deep Potential model to predict the intermetallic compound of the aluminum–magnesium system, and find six meta-stable phases with negative or nearly zero formation energy. In particular, Mg 12 Al 8 shows excellent ductility and Mg 5 Al 27 has a high Young's modulus. Based on our benchmark results, we propose a relatively robust structure screening criterion that selects potentially stable structures from the Deep Potential-based convex hull and performs DFT refinement. By using this criterion, the computational cost needed to construct the convex hull with ab initio accuracy can be dramatically reduced.
A model for the altitude dependence of the hot plasma parameters responsible for the electrostatic charging of spacecraft was developed. Based upon plasma orbit theory, the directed velocity is a function of the ambient magnetic field flux density. A consequence of this approach is that while the thermal velocity distributions (assumed to be Maxwellian) of the plasma particles are independent of the magnetic field strength (and hence altitude), the particle densities increase with magnetic field strength. Thus, according to this model, while the equilibrium voltage is independent of altitude, the charging current density increases with decreasing altitude. However, the probability of such spacecraft charging decreases with decreasing altitude.
With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.