Engineering Papers⌕ Search

Engineering topics

Liu, Ying

Publications and source records attributed to Liu, Ying.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Observation of odd-parity superconductivity in UTe 2

Symmetry properties of the order parameter are among the most fundamental characteristics of a superconductor. UTe 2 , which was found to feature an exceedingly large upper critical field and striking reentrant behavior at low temperatures, is widely believed to possess a spin-triplet pairing symmetry. However, unambiguous evidence for such a pairing symmetry is still lacking, especially at zero and low magnetic fields. The presence of an inversion crystalline symmetry in UTe 2 requires that, if it is indeed a spin-triplet superconductor, the order parameter must be of odd parity. We report here phase-sensitive measurements of the symmetry of the orbital part of the order parameter using the Josephson effect. The selection rule in the orientation dependence of the Josephson coupling between In, an s -wave superconductor, and UTe 2 suggests strongly that UTe 2 possesses the odd-parity pairing state of B 1 u symmetry near zero magnetic field, making it a spin-triplet superconductor. We also report the apparent formation of Andreev surface bound states on the (1−10) surface of UTe 2 .

Li, Zixuan (ORCID:0009000572025666)↗
DOE OSTI2024

An efficient construction of nano-interfaces for excellent coking tolerance of cermet anodes

Solid oxide fuel cells (SOFCs) are promising energy conversion devices for the effective and convenient utilization of hydrocarbons (for example, methane) to electricity. However, the development of direct methane SOFCs is primarily hindered by the poor coking tolerance of the state-of-the-art Ni-based cermet anodes. Herein, we efficiently construct nano-interfaces in the anode by infiltrating a Ni 0.6 Y 0.064 Zr 0.336 O 2-δ (NYZ) catalyst onto the traditional Ni-based cermet anode to effectively enhance the coking tolerance. After being reduced in H 2 , Ni and Y 0.16 Zr 0.84 O 2-δ (YSZ) nanoparticles (NPs) are in situ formed on the surface of the Ni-YSZ substrate. The roughened anode demonstrates significantly improved fuel oxidation activity and coking tolerance, due likely to the formation of nano-interfaces. Specifically, when applied in the Ni-YSZ-based anode-supported SOFCs, a high peak power density of 1.785 W cm –2 and a stable operation of ~ 240 h with no observable degradation is achieved at 750 °C in nearly dry methane (3% H 2 O). Finally, a density functional theory study suggests that the excellent coking tolerance is attributed to the formation of OH species on Ni/YSZ nano-interfaces, which would further interact with intermediate carbon species to generate COH intermediates.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Intense formation of secondary ultrafine particles from Amazonian vegetation fires and their invigoration of deep clouds and precipitation

New particle formation (NPF) in fire smoke is thought to be unlikely due to large condensation and coagulation sinks that scavenge molecular clusters. We analyze aircraft measurements over the Amazon and find that fires significantly enhance NPF and ultrafine particle (UFP < 50 nm diameter) numbers compared to background conditions, contrary to previous understanding. We identify that the nucleation of dimethylamine with sulfuric acid, which is aided by the formation of extremely low volatility organics in biomass-burning smoke, can overcome the large condensation and coagulation sinks and explain aircraft observations. We show that freshly formed clusters rapidly grow to UFP sizes through biomass-burning secondary organic aerosol formation, leading to a 10-fold increase in UFP number concentrations. Here, we find a contrasting effect of UFPs on deep convective clouds compared to the larger particles from primary emissions for the case investigated here. UFPs intensify the deep convective clouds and precipitation due to increased condensational heating, while larger particles delay and reduce precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Phase segregation of a composite air electrode unlocks the high performance of reversible protonic ceramic electrochemical cells

One breakthrough in developing highly efficient air electrodes for reversible protonic ceramic electrochemical cells (R-PCECs) is optimizing the sluggish oxygen reduction and water oxidation reactions. Here, we present a novel composite material with a nominal formula of high-entropy Ce 0.2 Ba 0.2 Sr 0.2 La 0.2 Ca 0.2 CoO 3–δ (CBSLCC) that spontaneously self-assembles to three-phase electrocatalysts composed of deficient Ce 0.2–y Ba 0.2 Sr 0.2–x La 0.2–x Ca 0.2 CoO 3–δ (CD-CBSLCC), CeO 2 , and La 0.5 Sr 0.5 CoO 3–δ (LSC). Mechanistic studies corroborate that oxygen reduction may occur on entire air electrode surfaces, followed by water formation preferentially at or near CD-CBSLCC. The CeO 2 phase could provide or consume protons to facilitate the oxygen evolution/reduction kinetics in R-PCECs. The developed electrodes demonstrate a record-high electrochemical performance in dual modes of fuel cells and electrolysis cells, delivering a peak power density of 1.66 W cm –2 at 600 °C and a current density of –1.76 A cm –2 at 1.3 V and 600 °C. Finally, excellent operational stabilities of the fuel cell (200 h at 600 °C), electrolysis cell (200 h at 600 °C), and reversible cycling (548 h at 550 °C) provide a promising and reliable step towards realizing the commercialization of R-PCECs.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Bending-induced isostructural transitions in ultrathin layers of van der Waals ferrielectrics

Using Landau-Ginzburg-Devonshire (LGD) phenomenological approach we analyze the bending-induced re-distribution of electric polarization and field, elastic stresses and strains inside ultrathin layers of van der Waals ferrielectrics. We consider a CuInP2S6 (CIPS) thin layer with fixed edges and suspended central part, the bending of which is induced by external forces. The unique aspect of CIPS is the existence of two ferrielectric states, FI1 and FI2, corresponding to big and small polarization values, which arise due to the specific four-well potential of the eighth-order LGD functional. When the CIPS layer is flat, the single-domain FI1 state is stable in the central part of the layer, and the FI2 states are stable near the fixed edges. With an increase of the layer bending below the critical value, the sizes of the FI2 states near the fixed edges decreases, and the size of the FI1 region increases. When the bending exceeds the critical value, the edge FI2 states disappear being substituted by the FI1 state, but they appear abruptly near the inflection regions and expand as the bending increases. The bending-induced isostructural FI1-FI2 transition is specific for the bended van der Waals ferrielectrics described by the eighth (or higher) order LGD functional with consideration of linear and nonlinear electrostriction couplings. Furthermore, the isostructural transition, which is revealed in the vicinity of room temperature, can significantly reduce the coercive voltage of ferroelectric polarization reversal in CIPS nanoflakes, allowing for the curvature-engineering control of various flexible nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stress and Curvature Effects in Layered 2D Ferroelectric CuInP 2 S 6

Nanoscale ferroelectric 2D materials offer the opportunity to investigate curvature and strain effects on materials functionalities. Among these, CuInP 2 S 6 (CIPS) has attracted tremendous research interest in recent years due to combination of room temperature ferroelectricity, scalability to a few layers thickness, and ferrielectric properties due to coexistence of 2 polar sublattices. Here, we explore the local curvature and strain effect on polarization in CIPS via piezoresponse force microscopy and spectroscopy. To explain the observed behaviors and decouple the curvature and strain effects in 2D CIPS, we introduce the finite element Landau–Ginzburg–Devonshire model, revealing strong changes in hysteresis characteristics in regions subjected to tensile and compressive strain. The piezoresponse force microscopy (PFM) results show that bending induces ferrielectric domains in CIPS, and the polarization-voltage hysteresis loops differ in bending and nonbending regions. Furthermore, these studies offer insights into the fabrication of curvature-engineered nanoelectronic devices.

2D materials↗
DOE OSTI2023

Ferroelectric Schottky diodes of CuInP 2 S 6 nanosheet

Ferroelectricity in van der Waals (vdW) layered material has attracted a great deal of interest recently. CuInP 2 S 6 (CIPS), the only vdW layered material whose ferroelectricity in the bulk was demonstrated by direct polarization measurements, was shown to remain ferroelectric down to a thickness of a few nanometers. However, its ferroelectric properties have just started to be explored in the context of potential device applications. Here, we report here the preparation and measurements of metal-ferroelectric semiconductor-metal heterostructures using nanosheets of CIPS obtained by mechanical exfoliation. Four bias voltage and polarization dependent resistive states were observed in the current–voltage characteristics, which we attribute to the formation of ferroelectric Schottky diode, along with switching behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sensitivity of Precipitation Displacement of a Simulated MCS to Changes in Land Surface Conditions

Abstract This study investigates the role of the land surface on the precipitation produced by an elevated mesoscale convective system (MCS) in Iowa between 24–25 June 2015 during the Plains Elevated Convection at Night (PECAN) field campaign. Previous studies have shown a strong effect of low‐level atmospheric moisture on the location of this MCS. A series of semi‐idealized and realistic simulations with irrigation are conducted to understand the effect of moisture perturbations on the MCS precipitation displacement. In general, numerical simulations place the MCS east of the observed location. Adding moisture directly in the low‐level atmosphere in the semi‐idealized experiments reduces this displacement error. However, experiments with perturbed soil moisture result in drying over Iowa induced by moisture flux divergence from cooler low‐level temperatures and higher surface pressure, causing the MCS to move further to the east. The irrigation impact on low‐level moisture is highly dependent on the length of simulation period. Shorter simulations on the order of days generate similar drying over Iowa but the opposite is found for month‐long simulations. Despite the lack of low‐level moistening in the perturbed soil moisture and short‐term irrigation experiments, the sensitivity to low‐level moisture is similar in all runs. More low‐level moisture generates a more convectively unstable environment with less inhibition and a lower level of free convection that leads to more rapid MCS development and a change in MCS location.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling the Impact of the Organic Aerosol Phase State on Multiphase OH Reactive Uptake Kinetics and the Resultant Heterogeneous Oxidation Timescale of Organic Aerosol in the Amazon Rainforest

Accurate predictions of chemical lifetime, i.e., oxidation timescale and change in the mass of organic aerosols (OAs) in atmospheric models, are critical to quantify the impacts of OA on aerosol–cloud interactions, radiative forcing, and air quality. The heterogeneous oxidation of OA by hydroxyl (OH) radicals is a key process governing OA mass changes during their oxidation. Recently, laboratory observations of the OH uptake coefficient (γ OH ) and heterogeneous reaction rate (k het,OH ) for different OA systems with varying phase states have been combined to develop a new parameterization of γ OH and k het,OH as a function of the OA phase state, i.e., viscosity (η OA ). Here, in this work, we use a recently published viscosity prediction framework to analyze the new γ OH and k het,OH parameterization with a box model. Subsequently, we implement this box model to predict γ OH and k het,OH over the Amazon rainforest in the dry-to-wet transition season with both significant biomass burning and biogenic influences. Relevant parameters within our box model are specified based on detailed regional model simulations over the Amazon using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). η OA is predicted as a function of species volatility, OA composition (including water uptake) along with ambient relative humidity (RH), and temperature. Based on ambient conditions simulated by WRF-Chem over the entire atmospheric column, we use the box model to predict the upper bounds of the heterogeneous oxidation timescale of OA. Based on previous laboratory measurements, we assume that this heterogeneous OH oxidation causes fragmentation (carbon loss) of OA, followed by the evaporation of the fragmented molecules that results in the exponential decay of OA. We predict that the oxidation timescale of OA is ~1 month near the Earth’s surface in the pristine Amazonian background because of low OH concentrations. But OA is oxidized more rapidly within urban and wildfire plumes near the Earth’s surface with an order of magnitude higher OH concentrations compared to the pristine background, causing the simulated oxidation timescale of OA to be much shorter, ~3–4 days. At 3–5 km altitudes where biomass burning OA is predicted to be semisolid, the heterogeneous oxidation timescale is estimated as ~6 days to 3 weeks and decreases with increasing OH concentrations within plumes. We show that the simulated mass loss of OA is strongly size-dependent, where smaller particles oxidize more rapidly compared to larger particles due to their greater per-particle surface area-to-volume ratio. Increasing urbanization and deforestation in the Amazon in the future might increase OH concentrations in the background Amazon, causing faster oxidation of OA. At higher altitudes above liquid clouds, especially in the upper troposphere where temperatures approach ~250 K, future measurements are needed to reduce the uncertainties related to the mass loss of OA at colder temperatures and low relative humidity conditions due to its heterogeneous OH oxidation.

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

Sharpening of cold-season storms over the western United States

Winter storms are responsible for billion-dollar economic losses in the western US. As storm structures are not well resolved by global climate models, event-scale storm response to warming has not been well established. Here regional storm-resolving simulations are used to investigate climate change impact on western US winter storms. In mid-century under a high emissions scenario, precipitation volume from the top 20% of winter storms is projected to increase by up to 40% across the region. Increased area coverage and storm intensity explain 69% and 60% of this increase, respectively, while a robust sharpening of future storms with larger increase in storm center precipitation than increase in storm area offsets the increase by 31%. Ignoring the changing area-intensity relationship due to storm sharpening could overestimate the watershed design storms by up to 10% or 7-75mm across the western US. Future storm structural changes have broad infrastructure, hydrological, and ecological implications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Two-dimensional ruthenium boride: a Dirac nodal loop quantum electrocatalyst for efficient hydrogen evolution reaction

Catalysts with high carrier mobility, high activity, and an active basal plane have been highly sought for the hydrogen evolution reaction (HER). However, combining these advantages into one single material is a grand challenge. Herein, using first principles computations, we predicted that a two-dimensional (2D) Dirac nodal loop semimetal, namely the RuB 4 monolayer, is promising as a superior catalyst for the HER. Our systematic computations showed that the single layer RuB 4 is thermodynamically, dynamically, mechanically, and thermally stable and presents multiple Ru and B sites for the HER on the basal plane. The estimated Gibbs free energy for hydrogen adsorption at a Ru site is approaching zero (–8.8 meV), suggesting its excellent HER performance. The RuB 4 monolayer is a Dirac nodal loop semimetal with high Fermi velocities, which can accelerate charge transfer between catalysts and reaction intermediates. The RuB 4 monolayer is an auxetic material with an out-of-plane negative Poisson's ratio, implying its novel mechanical properties. Finally, this work provides an example of using a Dirac nodal loop semimetal for high-performance HER catalysts, which is a promising alternative to the known catalysts with trivial metallic properties.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Formation Process of Particles and Cloud Condensation Nuclei Over the Amazon Rainforest: The Role of Local and Remote New-Particle Formation

Understanding the formation processes of particles and cloud condensation nuclei (CCN) in pristine environments is a major challenge in accurately assessing anthropogenic radiative forcing. Here, using a state-of-the-art model that systematically simulates the chemical processes controlling new-particle formation (NPF) and multi-scale chemical transport, we find that NPF contributes ~90% of the particle number and ~80% of the CCN at 0.5% supersaturation (CCN0.5%) in the pristine Amazon boundary layer during the wet season. The corresponding contributions are only ~30% and ~20% during the dry season because of prevalent biomass burning. In both seasons, about half of the NPF-induced particles and ~85% of the NPF-induced CCN0.5% in the boundary layer originate from the long-range transport of new particles formed hundreds to thousands of kilometers away. Moreover, about 50–65% of the NPF-induced particles and 35–50% of the NPF induced CCN0.5% originate from the downward transport of new particles formed in the free troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗
DOE OSTI2022

Revealing the Electrophilic‐Attack Doping Mechanism for Efficient and Universal p‐Doping of Organic Semiconductors

Abstract Doping is of great importance to tailor the electrical properties of semiconductors. However, the present doping methodologies for organic semiconductors (OSCs) are either inefficient or can only apply to some OSCs conditionally, seriously limiting their general applications. Herein, a novel p‐doping mechanism is revealed by investigating the interactions between the dopant trityl tetrakis(pentafluorophenyl) borate (TrTPFB) and poly(3‐hexylthiophene) (P3HT). It is found that electrophilic attack of the trityl cations on thiophenes results in the formation of tritylated thiophenium ions, which subsequently induce electron transfer from neighboring P3HT chains to realize p‐doping. This unique p‐doping mechanism enables TrTPFB to p‐dope various OSCs including those with high ionization energy (IE ≈ 5.8 eV). Moreover, this doping mechanism endows TrTPFB with strong doping capability, leading to doping efficiency of over 80% in P3HT. The discovery and elucidation of this novel doping mechanism not only points out that strong electrophiles are a class of efficient p‐dopants for OSCs, but also provides new opportunities toward highly efficient doping of various OSCs.

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

Calibration of cloud and aerosol related parameters for solar irradiance forecasts in WRF-solar

Model parameters are a major source of uncertainty in numerical weather prediction. Recently, the Weather Research and Forecasting model with Solar extensions (WRF-Solar) has been upgraded by enhancing the treatment of sub-grid scale cloud and aerosols with augmentations of a sub-grid scale cloud scheme (CLD3) and an upgraded aerosol-aware Thompson-Eidhammer scheme (TE14). However, the value of model parameters associated with these parameterizations are assigned based on limited measurements or theoretical calculations. Calibrating the most sensitive parameters has the potential to improve solar irradiance predictions. Here, we adopted a multiobjective surrogate-based optimization (SBO) framework to calibrate nine parameters used in CLD3 and TE14 that lead to the largest sensitivity in simulated irradiance. The normalized mean-absolute-error (NMAE) of global horizontal irradiance (GHI) and direct normal irradiance (DNI) are minimized by calibrating WRF-Solar over two regions including the Southern Great Plains (SGP) and Central California, in order to focus on parameter calibration under cloudy conditions with different aerosol loading. The results show that generalized linear model (GLM)-based surrogate models approximate physical models well, particularly when the third order and three-way interaction terms are considered. The SBO framework efficiently searches the parameter space for optimal solutions with less computational costs than directly calibrating the physical model. We first calibrate CLD3 parameters over the less-polluted SGP region. Optimized CLD3 parameters alone result in NMAE reduction by 14% for the site-mean and up to 33% for individual cases over the SGP region. With further calibration of TE14 parameters over the Central California during active fire periods, the optimized parameters lead to over 20% reductions of NMAE. Our investigation reveals, however, that optimizing TE14 has a limited impact on irradiance simulations under less-polluted conditions in the SGP.

14 SOLAR ENERGY↗
DOE OSTI2022

Critical role of acceptor dopants in designing highly stable and compatible proton-conducting electrolytes for reversible solid oxide cells

Proton-conducting electrolytes are receiving increasing attention due to their high ionic conductivity at intermediate temperatures, enabling the operation of solid oxide cells with high energy efficiency at low cost. However, the effect of B-site dopants on the properties of doped barium hafnate-cerate electrolyte materials, especially in single cells under operating conditions, has not been systematically studied. Here we report our findings in the development of a series of proton-conducting electrolytes with a general formula of BaHf 0.1 Ce 0.7 R 0.2 O 3–δ (BHCR172, R = Yb, Er, Y, Gd, Sm). Here, the results reveal that electrical conductivity, ionic transference number, chemical stability against steam and CO 2 , and compatibility with NiO during sintering are all closely correlated with the dopant size. In particular, the reaction with NiO is found to strongly affect the properties of the electrolytes and hence cell performance. Among all tested compositions, BaHf 0.1 Ce 0.7 Yb 0.2 O 3–δ (BHCYb172) shows excellent chemical stability and minimal reactivity towards NiO, as predicted from density functional theory (DFT)-based calculations and confirmed by experimental results. In addition, proton-conducting reversible solid oxide cells (P-ReSOCs) based on the optimized electrolyte composition, BHCYb172, demonstrate exceptional performance and stability, achieving a remarkable peak power density of 1.74 W cm –2 (O 2 as the oxidant) at 600 °C in the fuel cell mode and a high current density of 2.0 A cm –2 at 1.3 V and 600 °C in the steam electrolysis mode while maintaining excellent durability for over 1000 h.

25 ENERGY STORAGE↗