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 Nuclear Data Interface 2.2.0 has been released. It includes the addition of charged particle dE/dx (stopping power) data, the ability to read in NDI-formatted binary data, corrections to TN data, and other minor changes and bug fixes.
Generating high-fidelity spin-spin entanglement is an essential task of quantum repeater networks for the distribution of quantum information across long distances. Solid-state-based spin-photon interfaces are promising candidates to realize nodes of a quantum network, but are often limited by spectral diffusion of the optical transition, which results in phase errors on the entangled states. Here, we introduce a method to correct phase errors from quasi-static frequency fluctuations after the entangled state is generated, by shelving the emitters in the excited state to refocus the unknown phase. For quasi-static frequency fluctuations, the fidelity is determined only by the lifetime of the excited state used for shelving, making it particularly suitable for systems with a long-lived shelving state with correlated spectral diffusion. Such a shelving state may be found in Kramers doublet systems such as rare-earth emitters and color centers in Si or SiC interfaced with nanophotonic cavities with a strongly frequency-dependent Purcell enhancement. Furthermore, the protocol can be used to generate high-fidelity entangled spin pairs without reducing the rate of entanglement generation.
Bispecific antibodies represent an increasingly large fraction of biologics in therapeutic development due to their expanded scope in functional capabilities. Asymmetric monovalent bispecific IgGs (bsIgGs) have the additional advantage of maintaining a native antibody-like structure, which can provide favorable pharmacology and pharmacokinetic profiles. The production of correctly assembled asymmetric monovalent bsIgGs, however, is a complex engineering endeavor due to the propensity for non-cognate heavy and light chains to mis-pair. Previously, we introduced the DuetMab platform as a general solution for the production of bsIgGs, which utilizes an engineered interchain disulfide bond in one of the C H 1-C L domains to promote orthogonal chain pairing between heavy and light chains. While highly effective in promoting cognate heavy and light chain pairing, residual chain mispairing could be detected for specific combinations of Fv pairs. Here, we present enhancements to the DuetMab design that improve chain pairing and production through the introduction of novel electrostatic steering mutations at the C H 1-C L interface with lambda light chains (C H 1-C λ ). These mutations work together with previously established charge-pair mutations at the C H 1-C L interface with kappa light chains (C H 1-C κ ) and Fab disulfide engineering to promote cognate heavy and light chain pairing and enable the reliable production of bsIgGs. Importantly, these enhanced DuetMabs do not require engineering of the variable domains and are robust when applied to a panel of bsIgGs with diverse Fv sequences. We present a comprehensive biochemical, biophysical, and functional characterization of the resulting DuetMabs to demonstrate compatibility with industrial production benchmarks. Overall, this enhanced DuetMab platform substantially streamlines process development of these disruptive biotherapeutics.
The chemical potential of metal atoms, μM, in supported metal nanoparticles is an important descriptor related to both the catalytic activity and the stability of the nanoparticles. Here, we derive an expression relating μM to the radius of the particle’s contact area with the support and the adhesion energy at the metal/support interface (E adh ) that assumes the particles have the shape of spherical caps but of arbitrary contact angle with the support (θ c ) and includes an empirical correction for the increase in metal surface energy and adhesion energy with decreasing radius of curvature. We then show that, at any assumed contact angle, we can simultaneously fit previously reported measurements of both calorimetric μ M (from heats of metal vapor adsorption during nanoparticle growth by vapor deposition) versus metal coverage data and the He + low-energy ion scattering (LEIS) intensities for the metal and/or support versus metal coverage (using our recently developed spherical cap model for quantitative LEIS intensities), to determine the particle size versus coverage and E adh . Only one choice of contact angle gives a pair of values for contact angle and E adh , which is consistent with the Young–Dupré equation for the equilibrium shape of a spherical particle. At this equilibrium shape, we then applied this spherical cap model (SCM) to reanalyze microcalorimetric metal chemical potentials and LEIS signals versus coverage data for nine metal/support combinations that were previously analyzed by assuming that the particles had the shape of hemispherical caps, i.e., with a contact angle of 90°. We show that this revised approach gives close agreement with the calorimetric and LEIS data; the best-fit contact angles vary from 64 to 84°, correcting the earlier assumption of 90°. Furthermore, these results provide significant accuracy improvements in particle size versus coverage, metal chemical potential versus size and coverage, metal/support adhesion energies and contact angles for Cu, Ag and Au on CeO 2 (111), Ni on MgO(100), Ag on Fe 3 O 4 (111) and TiO 2 (100), and Ag, Ni and Pd on Ni-supported graphene. This revised approach is much more broadly applicable than the earlier hemispherical cap model (HCM).
The Nuclear Data Interface 2.2.0alpha has been released. It includes the addition of charged particle dE/dx (stopping power) data, the ability to read in NDI-formatted binary data, corrections to TN data, CP 2011 and CP 2020 data, and other minor changes and bug fixes.
Multipole aberration correctors have given scanning transmission electron microscopes (STEM) the ability to produce high-quality, atomic-resolution images, enabling STEM to be a key tool in material sciences for characterizing the structure and composition of materials. However, the process of correcting these aberrations typically requires human input and is accomplished using scanned STEM or Ronchigram images at different focii and beam tilts. In this work, we demonstrate an automated on-sample aberration correction system using Bayesian Optimization. We have developed a Python-based server able to communicate with the CEOS DCOR aberration corrector and the Thermo Fischer microscope scripting interface. This server allows us to change aberrations, acquire images and perform basic image analysis.
Here, the quest for a deeper understanding of the mechanical properties of HCP alloys has been an ongoing effort. Due to the lack of a sufficient number of slip systems, twin nucleation and migration originate as a substantial plastic deformation mechanism in HCP metals. The paper uncovers the energetics and atomistic mechanism of the {$11\overline{2}1$} < $11\overline{2}6$ > which is one of the most prominent twin modes in HCP materials that has not been fully understood. In view of corrugated {$11\overline{2}1$} interfaces, the positioning of lattice/motif atoms needs to be clarified. First, we establish the initial twin boundary (TB) structure unambiguously via the lattice offsets resulting in volume invariancy and correct relative positioning of the twin variants. Secondly, the TB migration mechanism is determined from a crystallographic (zero net shuffle, $Σ\vec{s}$ = 0) and energetic standpoints (minimum energy path). The Generalized Planar Fault Energy (GPFE) curves are calculated for Ti, Hf, Mg, and Zr. Such GPFE curves for this twin mode are calculated for the first time, and point to unusually high unstable energies (300-850 mJ/m 2 ), and low migration barriers (<60 mJ/m 2 ) for twinning accompanied with large shuffle magnitudes, |$\vec{s}$|, comparable to the Burger’s vector (|$\vec{s}$| ≅ $b_e$ > 0.44Å)
DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.
Reaction rate coefficients for electron-transfer processes at the electrode–electrolyte interface are commonly estimated by using the Butler–Volmer equation, but their values are inaccurate beyond a few tenths of volts of overpotential. The Marcus–Hush–Chidsey (MHC) formalism yields correct asymptotic behavior of the rate coefficients vs applied overpotential but has complex dependencies on the redox system’s intrinsic parameters, which can be difficult to model or measure. In this work, we bridge the two kinetics formalisms to estimate the reorganization energy, one of the important parameters for the MHC formalism, and investigate its dependence on other intrinsic parameters such as activation barriers, electronic coupling strength, and the density of states of the electrode surface. We examine the sensitivity of the reorganization energy to these parameters, establish some general relationships for accurately predicting rate coefficients using the MHC formalism over a wide range of applied overpotentials, and compare this approach to calculating MHC rate constants with other empirical approaches for the mechanisms of CO 2 reduction on different metal electrode surfaces.
Conductivity modulation via ferroelectric polarization coupling with LaNiO 3 (LNO) is demonstrated at an epitaxial ferroelectric–LNO interface. Conductivity measurements, varying the thickness of the LNO channel, show that this phenomenon is confined to a few atomic layers at the interface. Combining in situ biasing and off-axis holography, we mapped out electrostatic potentials at the PbZr 0.2 Ti 0.8 O 3 (PZT)/LNO/SrTiO 3 (STO) heterostructure upon polarization switching. Using aberration-corrected STEM, the interfacial atomic structures were investigated for the two different PZT polarization states. Polarization in PZT induces a significant change in the in-plane O–Ni–O bond angles, with a 37° modulation in the topmost 1 or 2 LNO unit cells, driven by strain in the oxygen sublattice for the two opposite polarization directions in PZT. Both oxygen and cation sublattices exhibit strain responses exceeding 10% upon switching. This atomic-layer structural modulation highlights a mechanism for functional oxide heterostructure development, offering pathways for advancements in nonvolatile memory, sensors, and energy-efficient transistors.
In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.
The US Department of Energy’s SunShot goals look to reduce the cost of Concentrating Solar Power (CSP) technology to 5¢/kWh for baseload plants. This is about a 50% reduction from current costs. To achieve this cost target, a significant reduction in operation and maintenance (O&M) costs of 40 to 50% is likely needed. Advances are needed in the O&M practices of CSP plants if the technology is to achieve the SunShot cost goals. Digitization of plant performance and O&M data has become a new best practice in the world of renewable energy asset management. Owners and operators of large photovoltaic and wind power plants are working to digitize performance and O&M data at their existing assets, to improve their management of the facilities, to increase performance, reduce O&M costs, and lower the overall life cycle cost of ownership. CSP power plants are behind the curve of other technologies on the digitization of plant information to aid in the plant asset management. This project directly addresses the objective of digitizing the O&M data of the solar field, focusing on three areas: 1) creating a framework for sharing data and information, 2) creating a system for monitoring and managing the maintenance of the solar field collectors, and 3) developing analytic tools to identify issues in the solar field. According to the NREL CSP Best Practices Study, the current practice at many CSP plants is to rely on paper lists, spreadsheets, and email for monitoring and managing problems and maintenance in the solar field. The key element to digitize solar field O&M is the creation of a centralized data archive that all users and systems can interface with. This project developed a centralized relational database framework that allows users and applications to access and share data. Conventional power plants utilize Computerized Maintenance Management Systems (CMMSs) to track the corrective, preventive (scheduled), and predictive maintenance of equipment and subsystems in the power plant. CSP plants use these systems in the power block, but while these systems specialize at tracking maintenance on up to thousands of pieces of equipment, they are not well suited for tracking the tens or hundreds of thousands of components in large commercial CSP or photovoltaic solar fields. In this project we developed a new software application referred to as FieldStatus (TM). This is a specialized database program that is used to track the status of each collector and its components. This application is designed to complement the existing CMMS to enable improved tracking and management of maintenance activities in the solar field. One of the major maintenance tasks for solar fields is maintaining the cleanliness of the mirrors. Although seemingly a relatively straight forward task, it has often proven challenging to maintain high levels of cleanliness in an efficient and cost-effective manner. This project developed new tools and metrics for monitoring and optimizing solar field cleaning resources and overall solar field cleanliness.
Abstract We present a simple criterion to predict the explodability of massive stars based on the density and entropy profiles before collapse. If a pronounced density jump is present near the Si/Si–O interface, the star will likely explode. We develop a quantitative criterion by using ∼1300 1D simulations where ν -driven turbulence is included via time-dependent mixing-length theory. This criterion correctly identifies the outcome of the supernova more than 90% of the time. We also find no difference in how this criterion performs on two different sets of progenitors, evolved using two different stellar evolution codes: FRANEC and KEPLER. The explodability as a function of mass of the two sets of progenitors is very different, showing: (i) that uncertainties in the stellar evolution prescriptions influence the predictions of supernova explosions; (ii) the most important properties of the pre-collapse progenitor that influence the explodability are its density and entropy profiles. We highlight the importance that ν -driven turbulence plays in the explosion by comparing our results to previous works.
We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.
In this project, we demonstrated stable nanoscale fracture in single-crystal silicon using an in-situ wedge-loaded double cantilever beam (DCB) specimen. The fracture toughness KIC was calculated directly from instrumented measurement of force and displacement via finite element analysis with frictional corrections. Measurements on multiple test specimens were used to show K IC = 0.72 ± 0.07 MPa m 1/2 on {111} planes and observe the crack-growth resistance curve in <500 nm increments. The exquisite stability of crack growth, instrumented measurement of material response, and direct visual access to observe nanoscale fracture processes in an ideally brittle material differentiate this approach from prior DCB methods.
In this project, the code for the fast-running Freitas wildfire plume model was successfully ported into the NARAC system, creating an operations-ready capability for predicting the transport of hazardous smoke from large wildfires. “Porting” refers to copying a section of code from one software package to another while accounting for unique dependencies in each software package. In this case, the code was ported from the Weather Research and Forecasting (WRF) Model, which is used extensively by the project team. To verify correct implementation, model results from the updated NARAC system were compared to those from WRF and the High-Resolution Rapid Refresh (HRRR) forecast model, which also employs the Freitas scheme. The NARAC User Interface (UI) and code documentation were also updated to include the new wildfire plume option, and a wildfire case was added to the NARAC test suite. The code was used in a real-time, semi-operational setup to predict smoke transport during the Cerro Pelado fire near Los Alamos National Laboratory in May 2022. The Freitas plume model was also tested for LLNL climate applications. Although the code was not ported into the Energy Exascale Earth System Model (E3SM), it was tested offline using E3SM atmospheric inputs, and results were compared to existing high-resolution plume predictions from large-eddy simulations in WRF. The code implementation and testing completed during this project will enable future work predicting smoke transport from large wildfires in LLNL national security and climate mission spaces.
As the power system continues to be flooded with intermittent resources, it becomes more important to accurately assess the role of hydro and its impact on the power grid. While hydropower generation has been studied for decades, dependency of power generation on water availability and constraints in hydro operation are not well represented in power system models used in the planning and operation of large-scale interconnection studies. There are still multiple modeling gaps that need to be addressed; if not, they can lead to inaccurate operation and planning reliability studies, and consequently to unintentional load shedding or even blackouts. As a result, it is very important thathydropower is represented correctly in both steady-state and dynamic power system studies. In this paper, we discuss the development and use of the Hydrological Dispatch and Analysis Tool (Hy-DAT) as an interactive graphical user interface (GUI), that uses a novel methodology to address the hydropower modeling gaps like water availability and inter-dependency using a database, and algorithms, to generate accurate representative models for power system simulation.
Niobium is commonly used for superconducting quantum systems as readout resonators, capacitors, and interconnects. Structural defects at the Nb/Si and air/Nb interface may be a major source of two-level systems (TLS), which are detrimental to the device's coherence time. Thus, identifying and understanding the microscopic origin of possible TLS in Nb-based devices and their relationship to processing is key to superconducting qubit performance improvement. Here this work studied the structure and thermal stability of surface oxide on physical vapor deposited Nb films on Si wafers, using aberration-corrected (scanning) transmission electron microscopy and spectroscopy. Here, all Nb films exhibit columnar growth with strong [110] textures. After in-situ heating of the heterostructure at 360 °C inside the microscope, the initial amorphous niobium surface oxides decompose into face-centered cubic Nb nanograins in the amorphous Nb-O matrix, which may reduce microwave dissipation. Despite changes in the microstructure and chemistry of the niobium oxide surface layer due to heat treatment, the interface between the Nb and the surface oxide layer remains almost unchanged. Our comprehensive study of the Nb surface oxide decomposition mechanism may guide future superconducting qubit device optimization through interfacial scattering center and TLS minimization.