Search NASA⌕ Search

Engineering topics

Li, Zhao

Publications and source records attributed to Li, Zhao.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2024

The Impact of the 8–10 March 2012 Geomagnetic Storm on Inner Zone Protons as Measured by Van Allen Probes

The Relativistic Electron Proton Telescope (REPT) instrument on the Van Allen Probes observed a double-peaked inner zone proton population throughout the 7 year lifetime of the mission. M. Hudson et al. (2023) showed that a strong SEP event accompanied by a CME-shock in early March 2012 provided the Solar Energetic Proton (SEP) source for the higher L trapped proton population, which then diffused radially inward to be observed by REPT at L = 2. The study followed trajectories of SEP protons launched isotropically from a sphere at 7 Re for 2.5 hr in fields calculated by the LFM-RCM global MHD model, which includes electric fields needed to model the transport and trapping of the protons by the shock, and then a radial diffusion simulation was run for 2 years using the result from the test-particle simulation as the initial condition. The simulation result was compared with REPT measurement in November 2013 and showed reasonable agreement. However, the simulation overestimated the Phase Space Density by a factor of four due to lack of field line curvature scattering during the storm in the model. In this study, a test-particle simulation is performed for 2 days following the injection and trapping of protons in March 2012 using TS05 fields to simulate the field line curvature scattering of the trapped SEP due to the buildup of the ring current during the geomagnetic storm. The resulting sample distribution was then weighted using the flux at the end of the two-hour MHD-test particle simulation. A radial diffusion simulation is then run using the initial profile that included the loss effect, with improved comparison with REPT measurements after 2 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamic Bubbling Balanced Proactive CO 2 Capture and Reduction on a Triple-Phase Interface Nanoporous Electrocatalyst

The formation and preservation of the active phase of the catalysts at the triple-phase interface during CO 2 capture and reduction is essential for improving the conversion efficiency of CO 2 electroreduction toward value-added chemicals and fuels under operational conditions. Designing such ideal catalysts that can mitigate parasitic hydrogen generation and prevent active phase degradation during the CO 2 reduction reaction (CO 2 RR), however, remains a significant challenge. Herein, we developed an interfacial engineering strategy to build a new SnO x catalyst by invoking multiscale approaches. This catalyst features a hierarchically nanoporous structure coated with an organic F-monolayer that modifies the triple-phase interface in aqueous electrolytes, substantially reducing competing hydrogen generation (less than 5%) and enhancing CO 2 RR selectivity (~90%). This rationally designed triple-phase interface overcomes the issue of limited CO 2 solubility in aqueous electrolytes via proactive CO 2 capture and reduction. Concurrently, we utilized pulsed square-wave potentials to dynamically recover the active phase for the CO 2 RR to regulate the production of C1 products such as formate and carbon monoxide (CO). This protocol ensures profoundly enhanced CO 2 RR selectivity (~90%) compared with constant potential (~70%) applied at -0.8 V (V vs RHE). We further achieved a mechanistic understanding of the CO 2 capture and reduction processes under pulsed square-wave potentials via in situ Raman spectroscopy, thereby observing the potential-dependent intensity of Raman vibrational modes of the active phase and CO 2 RR intermediates. Finally, this work will inspire material design strategies by leveraging triple-phase interface engineering for emerging electrochemical processes, as technology moves toward electrification and decarbonization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Intermetallic Single-Atom Alloy In–Pd Bimetallene for Neutral Electrosynthesis of Ammonia from Nitrate

Harvesting recyclable ammonia (NH 3 ) from the electrocatalytic reduction of nitrate (NO 3 RR) offers a sustainable strategy to close the ecological nitrogen cycle from nitration contamination in an energy-efficient and environmentally friendly manner. The emerging intermetallic single-atom alloys (ISAAs) are recognized to achieve the highest site density of single atoms by isolating contiguous metalatoms into single sites stabilized by another metal within the intermetallic structure, which holds promise to couple the catalytic benefits from intermetallic nanocrystals and single-atom catalysts for promoting NO 3 RR. Herein, ISAA In-Pd bimetallene, in which the Pd single atoms are isolated by surrounding In atoms, is reported to boost neutral NO 3 RR with a NH 3 Faradaic efficiency(FE) of 87.2%, a yield rate of 28.06 mg h -1 mg Pd -1 , and an exceptional electrocatalytic stability with increased activity/selectivity over 100 h and 20 cycles. The ISAA structure induces substantially diminished overlap of Pdd-orbitals and narrowed p-d hybridization of In-p and Pd-d states around the Fermi level, resulting in a stronger NO 3 - adsorption and a depressed energy barrier of the potential-determining step for NO 3 RR. Finally, further integrating the NO 3 RR catalyst into a Zn-NO 3 - flow battery as the cathode delivers a power density of 12.64 mW cm -2 and a FE of 93.4% for NH 3 production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Two-Dimensional Energy Histograms as Features for Machine Learning to Predict Adsorption in Diverse Nanoporous Materials

A major obstacle for machine learning (ML) in chemical science is the lack of physically informed feature representations that provide both accurate prediction and easy interpretability of the ML model. In this work, we describe adsorption systems using novel two-dimensional energy histogram (2D-EH) features, which are obtained from the probe-adsorbent energies and energy gradients at grid points located throughout the adsorbent. The 2D-EH features encode both energetic and structural information of the material and lead to highly accurate ML models (coefficient of determination R2 ~ 0.94–0.99) for predicting single-component adsorption capacity in metal–organic frameworks (MOFs). Here, we consider the adsorption of spherical molecules (Kr and Xe), linear alkanes with a wide range of aspect ratios (ethane, propane, n-butane, and n-hexane), and a branched alkane (2,2-dimethylbutane) over a wide range of temperatures and pressures. The interpretable 2D-EH features enable the ML model to learn the basic physics of adsorption in pores from the training data. We show that these MOF-data-trained ML models are transferrable to different families of amorphous nanoporous materials. We also identify several adsorption systems where capillary condensation occurs, and ML predictions are more challenging. Nevertheless, our 2D-EH features still outperform structural features including those derived from persistent homology. The novel 2D-EH features may help accelerate the discovery and design of advanced nanoporous materials using ML for gas storage and separation in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

MOFX-DB: An Online Database of Computational Adsorption Data for Nanoporous Materials

Machine learning and data mining coupled with molecular modeling have become powerful tools for materials discovery. Metal-organic frameworks (MOFs) are a rich area for this due to their modular construction and numerous applications. Here, we make data from several previous large-scale studies in MOFs and zeolites from our groups (and new data for N 2 and Ar adsorption in MOFs) easily accessible in one place. The database includes over 3 million simulated adsorption data points for H 2 , CH 4 , CO 2 , Xe, Kr, Ar, and N 2 in over 160 000 MOFs and zeolites, textural properties like pore sizes and surface areas, and the structure file for each material. We include metadata about the Monte Carlo simulations to enable reproducibility. The database is searchable by MOF properties, and the data are stored in a standardized JSON format that that is interoperable with the NIST adsorption database. We also identify several MOFs that meet high performance targets for multiple applications, such as high storage capacity for both hydrogen and methane or high CO 2 capacity plus good Xe/Kr selectivity. Here, by providing this data publicly, we hope to facilitate machine learning studies on these materials, leading to new insights on adsorption in MOFs and zeolites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Three-dimensional Zn-based alloys for dendrite-free aqueous Zn battery in dual-cation electrolytes

Aqueous zinc-ion batteries, in terms of integration with high safety, environmental benignity, and low cost, have attracted much attention for powering electronic devices and storage systems. However, the interface instability issues at the Zn anode caused by detrimental side reactions such as dendrite growth, hydrogen evolution, and metal corrosion at the solid (anode)/liquid (electrolyte) interface impede their practical applications in the fields requiring long-term performance persistence. Despite the rapid progress in suppressing the side reactions at the materials interface, the mechanism of ion storage and dendrite formation in practical aqueous zinc-ion batteries with dual-cation aqueous electrolytes is still unclear. Herein, we design an interface material consisting of forest-like three-dimensional zinc-copper alloy with engineered surfaces to explore the Zn plating/stripping mode in dual-cation electrolytes. The three-dimensional nanostructured surface of zinc-copper alloy is demonstrated to be in favor of effectively regulating the reaction kinetics of Zn plating/stripping processes. The developed interface materials suppress the dendrite growth on the anode surface towards high-performance persistent aqueous zinc-ion batteries in the aqueous electrolytes containing single and dual cations. This work remarkably enhances the fundamental understanding of dual-cation intercalation chemistry in aqueous electrochemical systems and provides a guide for exploring high-performance aqueous zinc-ion batteries and beyond.

25 ENERGY STORAGE↗
DOE OSTI2022

Engineered disorder in CO 2 photocatalysis

Light harvesting, separation of charge carriers, and surface reactions are three fundamental steps that are essential for an efficient photocatalyst. Here we show that these steps in the TiO 2 can be boosted simultaneously by disorder engineering. A solid-state reduction reaction between sodium and TiO 2 forms a core-shell c-TiO 2 @a-TiO 2-x (OH) y heterostructure, comprised of HO-Ti-[O]-Ti surface frustrated Lewis pairs (SFLPs) embedded in an amorphous shell surrounding a crystalline core, which enables a new genre of chemical reactivity. Specifically, these SFLPs heterolytically dissociate dihydrogen at room temperature to form charge-balancing protonated hydroxyl groups and hydrides at unsaturated titanium surface sites, which display high reactivity towards CO 2 reduction. This crystalline-amorphous heterostructure also boosts light absorption, charge carrier separation and transfer to SFLPs, while prolonged carrier lifetimes and photothermal heat generation further enhance reactivity. The collective results of this study motivate a general approach for catalytically generating sustainable chemicals and fuels through engineered disorder in heterogeneous CO 2 photocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

New black indium oxide—tandem photothermal CO 2 -H 2 methanol selective catalyst

It has long been known that the thermal catalyst Cu/ZnO/Al 2 O 3 (CZA) can enable remarkable catalytic performance towards CO 2 hydrogenation for the reverse water-gas shift (RWGS) and methanol synthesis reactions. However, owing to the direct competition between these reactions, high pressure and high hydrogen concentration (≥75%) are required to shift the thermodynamic equilibrium towards methanol synthesis. Herein, a new black indium oxide with photothermal catalytic activity is successfully prepared, and it facilitates a tandem synthesis of methanol at a low hydrogen concentration (50%) and ambient pressure by directly using by-product CO as feedstock. The methanol selectivities achieve 33.24% and 49.23% at low and high hydrogen concentrations, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Surface oxygenation induced strong interaction between Pd catalyst and functional support for zinc–air batteries

Employing the strong metal-support interaction (SMSI) effect for promoting the catalyst's activity toward the oxygen reduction reaction (ORR) is promising due to the electronic structure optimization and high utilization efficiency of platinum group metal (PGM) catalysts. Metal oxides as alternative supports for PGMs facilitate intrinsic activity and improve durability as compared to conventional carbon supports. However, the restricted mass and electron transfer at the metal/support interface need to be addressed. Herein, to strengthen the interaction at the metal/support interfaces and improve the utilization efficiency of PGM, an ultralow loading of Pd was embedded in a surface-oxygenated PdNiMnO porous film. The Mn-doping was designed to promote surface oxygenation using a facile anodization process that created sufficiently exposed interfaces between Pd and the support, strengthening the SMSI effects at the Pd/oxygenated support interface for enhancing ORR performance. Furthermore, the Ni-containing oxygenated catalyst served as both the active component for the oxygen evolution reaction (OER) and the functional support for stabilizing Pd, making PdNiMnO a bifunctional catalyst for zinc–air flow batteries (ZAFB). As a proof-of-concept, the ZAFB (PdNiMnO) shows a maximal power density of 211.6 mW cm –2 and outstanding cycling stability for over 2000 h with a minimal voltage gap of 0.69 V at a current density of 10 mA cm –2 , superior to the state-of-the-art catalysts.

25 ENERGY STORAGE↗
DOE OSTI2021

Probe the Localized Electrochemical Environment Effects and Electrode Reaction Dynamics for Metal Batteries using In Situ 3D Microscopy

Uncontrollable dendrite growth is closely related to non-uniform reaction environments. However, there is a lack of understanding and analysis methods to probe the localized electrochemical environment (LEE). Here the effects of the LEE are investigated, including localized ion concentrations, current density, and electric potential, on metal plating/stripping dynamics and dendrite minimization. A novel in situ 3D microscopy technique is developed to image the morphology dynamics and deposition rate of Zn plating/stripping processes on 3D Zn–Mn anodes. Using the in situ 3D microscope, the electrode morphology changes during the reactions are directly imaged and Zn deposition rate maps at different time points are obtained. It is found that reaction kinetics are highly correlated to LEE and electrode morphology. To further quantify the LEE effects, the digital twin technique is employed that allows the accurate calculation of the electrochemical environments, such as localized ion concentrations, current density, and electric potential, which cannot be directly measured from experiments. We found that the curvature of the 3D electrode surface determines the LEE and significantly influences reaction kinetics. This provides a new strategy to minimize the dendrite formation by designing and optimizing the 3D geometry of the electrode to control the LEE.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Doping-modulated strain control of bifunctional electrocatalysis for rechargeable zinc–air batteries

Changes in the local atomic arrangement in a crystal caused by lattice-mismatch-induced strain can efficiently regulate the performance of electrocatalysts for zinc–air batteries (ZABs) in many manners, mainly due to modulated electronic structure configurations that affect the adsorption energies for oxygen-intermediates formed during oxygen reduction and evolution reactions (ORR and OER). However, the application of strain engineering in electrocatalysis has been limited by the strain relaxation caused by structural instability such as dissolution and destruction, leading to insufficient durability towards the ORR/OER. In this work, we propose a doping strategy to modulate the phase transition and formation of self-supported cobalt fluoride–sulfide (CoFS) nanoporous films using a low amount of copper (Cu) as a dopant. This well-defined Cu–CoFS heterostructure overcomes the obstacle of structural instability. Our study of the proposed Cu–CoFS also helps establish the structure–property relationship of strained electrocatalysts by unraveling the role of local strain in regulating the electronic structure of the catalyst. As a proof-of-concept, the Cu–CoFS electrocatalyst with doping-modulated strain exhibited superior onset potentials of 0.91 V and 1.49 V for the ORR and OER, respectively, surpassing commercial Pt/C@RuO 2 and benchmarking non-platinum group metal (non-PGM) catalysts. ZABs with the Cu–CoFS catalyst delivered excellent charge/discharge cycling performance with an extremely low voltage gap of 0.5 V at a current density of 10 mA cm –2 and successively 0.93 V at a high current density of 100 mA cm –2 and afforded an outstanding peak power density of 255 mW cm –2 .

25 ENERGY STORAGE↗
DOE OSTI2021

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗