Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

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.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Observation of large Rashba spin–orbit coupling at room temperature in compositionally engineered perovskite single crystals and application in high performance photodetectors

Indirect absorption extended below the direct transition edge and increase in carrier lifetime derived from Rashba spin–orbit coupling may advance the optoelectronic applications of metal halide perovskites. Spin-orbit coupling in halide perovskites is due to the presence of heavy elements in their structure. However, when these materials lack an inversion symmetry, for example by the application of strain, spin–orbit coupling becomes odd in the electron’s momentum giving rise to a splitting in the electronic energy bands. Here we report on the observation of a large Rashba splitting of 117 meV at room temperature, as predicted by relativistic first-principles calculations, in halide perovskite single crystals through a facile compositional engineering approach. Partial substitution of organic cations by rubidium in single crystals induces significant indirect absorption and dual peak photoluminescence as a result of a large Rashba splitting. We measured circularly polarized photoluminescence and magneto-photoluminescence in perovskite films printed by single crystals as well as magneto-electroluminescence and magneto-photocurrent in spin-LEDs based on perovskite single crystals. They indicated significant spin-momentum locking due to the large Rashba effect. A hybrid perovskite single crystal photodetector achieved record figures of merit, including detectivity of more than 1.3 × 10 18 Jones which represents a three orders of magnitude improvement compared to the to date record. Furthermore, these findings show that facile compositional engineering of perovskite single crystals holds great promise for further advancing the optoelectronic properties of existing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gamma-ray spectrum variations for surface measurements of uranium hexafluoride cylinders

Here, this work aims to characterize features that contribute to the spectral variations that appear in gamma-23 ray spectra of filled uranium hexafluoride (UF 6 ) cylinders. While the enrichment of the UF 6 in a 24 cylinder can greatly affect gamma-ray spectra at peaks such as the 185.7 keV peak and be used to 25 verify enrichment, many other characteristics of a filled UF 6 cylinder can affect the spectra. These 26 characteristics vary on a per-cylinder basis and can include: the time between when a cylinder was 27 filled and when it was measured; the previous fill history of a cylinder; how the cylinder was most 28 recently filled; the UF 6 geometric distribution within the cylinder; and the measurement location 29 along the cylinder. This work explores spectral variations seen in measurement data and uses 30 validated simulation models to find trends in these uncertainties. These trends are analyzed to observe 31 how each parameter can contribute to the overall uncertainty of UF 6 spectral measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Full-polarization millimeter wavelength variability of Sagittarius A * during the 2018 EHT campaign

Context. Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims. We investigated the variability and polarization properties of Sgr A* using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results. We find low variability in total intensity (σ/μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions. Our results offer new constraints on the physical mechanisms driving variability in Sgr A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

Galaxy: center↗
DOE OSTI2022

The Optical Two- and Three-dimensional Fundamental Plane Correlations for Nearly 180 Gamma-Ray Burst Afterglows with Swift/UVOT, RATIR, and the Subaru Telescope

Gamma-ray bursts (GRBs) are fascinating events due to their panchromatic nature. We study optical plateaus in GRB afterglows via an extended search into archival data. We comprehensively analyze all published GRBs with known redshifts and optical plateaus observed by many ground-based telescopes (e.g., Subaru Telescope, RATIR) around the world and several space-based observatories such as the Neil Gehrels Swift Observatory. We fit 500 optical light curves, showing the existence of the plateau in 179 cases. This sample is 75% larger than the previous one, and it is the largest compilation so far of optical plateaus. We discover the 3D fundamental plane relation at optical wavelengths using this sample. This correlation is between the rest-frame time at the end of the plateau emission, ${T}_{\mathrm{opt}}^{* }$, its optical luminosity, L opt , and the peak in the optical prompt emission, L peak,opt , thus resembling the three-dimensional (3D) X-ray fundamental plane (the so-called 3D Dainotti relation). We correct our sample for redshift evolution and selection effects, discovering that this correlation is indeed intrinsic to GRB physics. We investigate the rest-frame end-time distributions in X-rays and optical (${T}_{\mathrm{opt}}^{* }$, ${T}_{{\rm{X}}}^{* }$), and conclude that the plateau is achromatic only when selection biases are not considered. We also investigate if the 3D optical correlation may be a new discriminant between optical GRB classes and find that there is no significant separation between the classes compared to the Gold sample plane after correcting for evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗
DOE OSTI2020

Detailed Abundances in the Ultra-faint Magellanic Satellites Carina II and III

We present the first detailed elemental abundances in the ultra-faint Magellanic satellite galaxies Carina II (Car II) and Carina III (Car III). With high-resolution Magellan/MIKE spectroscopy, we determined the abundances of nine stars in Car II, including the first abundances of an RR Lyrae star in an ultra-faint dwarf galaxy (UFD), and two stars in Car III. The chemical abundances demonstrate that both systems are clearly galaxies and not globular clusters. The stars in these galaxies mostly display abundance trends matching those of other similarly faint dwarf galaxies: enhanced but declining [α/Fe] ratios, iron-peak elements matching the stellar halo, and unusually low neutron-capture element abundances. One star displays a low outlying [Sc/Fe] = −1.0. We detect a large Ba scatter in Car II, likely due to inhomogeneous enrichment by low-mass asymptotic giant branch star winds. The most striking abundance trend is for [Mg/Ca] in Car II, which decreases from +0.4 to −0.4 and indicates clear variation in the initial progenitor masses of enriching core-collapse supernovae. So far, the only UFDs displaying a similar [Mg/Ca] trend are likely satellites of the Large Magellanic Cloud. We find two stars with [Fe/H] ≤ −3.5 whose abundances likely trace the first generation of metal-free Population III stars and are well fit by Population III core-collapse supernova yields. An appendix describes our new abundance uncertainty analysis that propagates line-by-line stellar parameter uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Nature of Jets in Double-peaked Emission-line AGN in the KISSR Sample

Double-peaked emission-line AGN (DPAGN) have been regarded as binary black hole candidates. We present here results from parsec-scale radio observations with the Very Long Baseline Array (VLBA) of five DPAGN belonging to the KISSR sample of emission-line galaxies. This work concludes our pilot study of nine type 2 Seyfert and LINER DPAGN from the KISSR sample. In the nine sources, dual compact cores are only detected in the “offset AGN,” KISSR 102. However, the overall incidence of jets in the eight sources that were detected with the VLBA is ≥ 60%. We find a difference in the “missing flux density” going from the Very Large Array to VLBA scales between Seyferts and LINERs, with LINERs showing less missing flux density on parsec-scales. Using the emission-line modeling code, MAPPINGS III, we find that the emission lines are likely to be influenced by jets in 5/9 sources. Furthermore, jet-medium interaction is the likely cause of the emission-line splitting observed in the SDSS spectra of these sources. Jets in radio-quiet AGN are therefore energetically capable of influencing their parsec- and kpcscale environments, which makes them agents of “radio AGN feedback” (similar to radio-loud AGN).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Correlated energy uncertainties in reaction rate calculations

Monte Carlo methods can be used to evaluate the uncertainty of a reaction rate that arises from many uncertain nuclear inputs. However, until now no attempt has been made to find the effect of correlated energy uncertainties in input resonance parameters. Our goal is to investigate the impact of correlated resonance energy uncertainties on reaction rates. Using a combination of numerical and Monte Carlo variation of resonance energies, the effect of correlations are investigated. Five reactions are considered: two fictional, illustrative cases and three reactions whose rates are of current interest. The effect of correlations in resonance energies depends on the specific reaction cross section and temperatures considered. When several resonances contribute equally to a reaction rate, and when they are located on either side of the Gamow peak, correlations between their energies dilute their effect on reaction rate uncertainties. If they are both located above or below the maximum of the Gamow peak, however, correlations between their resonance energies can increase the reaction rate uncertainties. This effect can be hard to predict for complex reactions with wide and narrow resonances contributing to the reaction rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A multi-MeV alpha particle source via proton-boron fusion driven by a 10-GW tabletop laser

Nuclear fusion between protons and boron-11 nuclei has undergone a revival of interest thanks to the rapid progress in pulsed laser technology. Potential applications of such reaction range from controlled nuclear fusion to radiobiology and cancer therapy. A laser-driven fusion approach consists in the interaction of high-power, high-intensity pulses with H- and B-rich targets. We report on an experiment exploiting proton-boron fusion in CN-BN targets to obtain high-energy alpha particle beams (up to 5 MeV) using a very compact approach and a tabletop laser system with a peak power of ~10 GW, which can operate at high-repetition rate (up to 1 kHz). The secondary resonance in the cross section of proton-boron fusion (~150 keV in the center-of-mass frame) is exploited using a laser-based approach. The generated alpha particles are characterized in terms of energy, flux, and angular distribution using solid-state nuclear-track detectors, demonstrating a flux of ~10 5 particles per second at 10 Hz, and ~10 6 per second at 1 kHz. Hydrodynamic and particle-in-cell numerical simulations support our experimental findings. Potential impact of our approach on future spread of ultra-compact, multi-MeV alpha particle sources driven by moderate intensity (10 16 -10 17 W/cm 2 ) laser pulses is anticipated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Smoke from 2020 United States wildfires responsible for substantial solar energy forecast errors

Abstract The 2020 wildfire season (May through December) in the United States was exceptionally active, with the National Interagency Fire Center reporting over 10 million acres ( > 40 000 km 2 ) burned. During the September 2020 wildfire events, large concentrations of smoke particulates were emitted into the atmosphere. As a result, smoke was responsible for ∼10%–30% reduction in solar power production during peak hours as recorded by the California Independent System Operator (CAISO) sites. In this study, we focus on a 9 d period in September when wildfire smoke had a profound impact on solar energy production. During the smoke episodes, hour-ahead forecasts utilized by CAISO did not include the effects of smoke and therefore overestimated the expected power production by ∼10%–50%. Here we use multiple observational networks and a numerical weather prediction (NWP) model to show that the wildfire events of 2020 had a significantly detrimental influence on solar energy production due to high aerosol loading. We find that including the contribution of biomass burning particles greatly improves the day-ahead solar energy bias forecast of both global horizontal irradiance and direct normal irradiance by nearly ∼50%. Our results suggest that a more comprehensive treatment of aerosols, including biomass burning aerosols, in NWP models may be an important consideration for energy grid balancing, in addition to solar resource assessment, as solar power reliance increases.

14 SOLAR ENERGY↗
DOE OSTI2024

Band structure and Fermi surface nesting in LaSb 2

Here we use high-resolution angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) to investigate the electronic structure of the charge density wave (CDW) system LaSb 2 . This compound is among an interesting group of materials that manifests both a CDW transition and lower temperature superconductivity. We find the DFT calculations to be in good agreement with our ARPES data. The Fermi surface of LaSb 2 consists of two small hole pockets close to Γ and four larger pockets near the Brillouin zone boundary. The key features of the Fermi surface do not vary significantly with temperature. A saddle point is present at -0.19 eV below the Fermi level at Γ. Saddle points in the band structure have more pronounced effects on materials properties when they are located closer to the Fermi level, making doped LaSb 2 compounds a potential interesting subject of future research. Multiple peaks are present in the generalized, electronic susceptibility calculations, indicating the presence of three possible nesting vectors. We were not able to detect any signatures of the CDW transition at 355 K down to the lowest temperature of 7 K, pointing to the subtle nature of this transition. This is unusual, given that such a high transition temperature is expected to be associated with the presence of a large CDW gap. This is confirmed through investigation of the Fermi surface and through analysis of momentum distribution curves. It is possible that changes are subtle and occur below current sensitivity of our measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Detecting Long-period Variability in the SDSS Stripe 82 Standards Catalog

We report the results of a search for long-period (100 < P < 600 days) periodic variability in the SDSS Stripe 82 standards catalog. The SDSS coverage of Stripe 82 enables such a search because there are on average 20 observations per band in ugriz bands for about one million sources, collected over about 6 yr, with a faint limit of r ~ 22 mag and precisely calibrated 1%–2% photometry. We calculated the periods of variable source candidates in this sample using the Lomb–Scargle periodogram and considered the three highest periodogram peaks in each of the gri filters as relevant. Only those sources with gri periods consistent within 0.1% were later studied. We use the Kuiper statistic to ensure uniform distribution of data points in phased light curves. We present five sources with the spectra consistent with quasar spectra and plausible periodic variability. This SDSS-based search bodes well for future sensitive large-area surveys, such as the Rubin Observatory Legacy Survey of Space and Time, which, due to its larger sky coverage (about a factor of 60) and improved sensitivity (~2 mag), will be more powerful for finding such sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

InAs sidewall tunnel diodes enabled by surface states

Negative differential resistance (NDR), where the device current decreases with increasing bias voltage, is a representative phenomenon where quantum mechanics induces counterintuitive physical behavior and offers promising applications such as high-frequency oscillators, amplifiers, and multilevel logic circuits. While the NDR behavior has been extensively studied in various materials and devices, the role of surface properties in NDR, particularly in InAs-based diodes, remains underexplored. In this work, we report the observation of NDR in vertically structured InAs p + n diodes that exhibit a peak-valley current ratio of ∼6, which is suitably high for applications. Circumference-normalized current–voltage characterization revealed that the NDR originates from band-to-band tunneling between the valence band of p + -InAs and the conduction band of an n + -InAs surface, where the n + surface is due to surface states on otherwise n − -InAs. In addition, by comparing devices with various surface passivation methods (without intentional passivation, benzocyclobutene polymer, and silicon nitride), we found that the surface termination significantly affects the NDR characteristics. We present an equivalent circuit model to explain the observed device behavior. These findings offer insights into surface-enabled NDR phenomena and present new knobs for engineering NDR devices.

Circuit theorems↗
DOE OSTI2025

Measurements of the temperature and $E$-mode polarization of the cosmic microwave background from the full 500-square-degree SPTpol dataset

Using the full four-year SPTpol 500 deg 2 dataset in both the 95 and 150 GHz frequency bands, we present measurements of the temperature and E-mode polarization of the cosmic microwave background (CMB), as well as the E-mode polarization autopower spectrum (EE) and temperature-E-mode cross-power spectrum (TE) in the angular multipole range 50 < ℓ < 8000. We find the SPTpol dataset to be self-consistent, passing several internal consistency tests based on maps, frequency bands, bandpowers, and cosmological parameters. The full SPTpol dataset is well-fit by the ΛCDM model, for which we find H 0 = 70.48 ± 2.16 km s -1 Mpc -1 and Ω m = 0.271 ± 0.026, when using only the SPTpol data and a Planck-based prior on the optical depth to reionization. The ΛCDM parameter constraints are consistent across the 95 GHz-only, 150 GHz-only, TE-only, and EE-only data splits. Between the ℓ < 1000 and ℓ > 1000 data splits, the ΛCDM parameter constraints are borderline consistent at the ∼2σ level. This consistency improves when including a parameter A L , the degree of lensing of the CMB inferred from the smearing of acoustic peaks. When marginalized over A L , the ΛCDM parameter constraints from SPTpol are consistent with those from Planck. In conclusion, the power spectra presented here are the most sensitive measurements of the lensed CMB damping tail to date for roughly ℓ > 1700 in TE and ℓ > 2000 in EE.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗
DOE OSTI2025

Puzzling Variation of Gamma Rays from the Sun over the Solar Cycle Revealed with Fermi-LAT

The steady-state gamma-ray emission from the Sun is thought to consist of two emission components due to interactions with Galactic cosmic rays: (1) a hadronic disk component, and (2) a leptonic extended component peaking at the solar edge and extending into the heliosphere. The flux of these components is expected to vary with the 11 yr solar cycle, being highest during solar minimum and lowest during solar maximum, as it varies with the cosmic-ray flux. No study has yet analyzed the flux variation of each component over solar cycles. In this work, we measure the temporal variations of the flux of each component over 15 yr of Fermi Large Area Telescope observations and compare them with the sunspot number and Galactic cosmic-ray flux from AMS-02 near Earth. We find that the flux variation of the disk anticorrelates with the sunspot number and correlates with cosmic-ray protons, as expected, confirming its emission mechanism. In contrast, the extended component exhibits a more complex variation: despite an initial anticorrelation with the sunspot number, we find neither anticorrelation with the sunspot number nor correlation with cosmic-ray electrons over the full 15 yr period. This most likely suggests that cosmic-ray transport and modulation in the inner heliosphere are unexpectedly complex and may differ for electrons and protons or, alternatively, that there is an additional, unknown component of gamma rays or cosmic rays. These findings impact space weather research and emphasize the need for close monitoring of Cycle 25 and the ongoing polarity reversal.

Acharyya, A. [Univ. of Southern Denmark, Odense (D↗
DOE OSTI2023

Decomposition of liquid-liquid extraction organic phase structure into critical and pre-peak contributions

Solution nanostructure induced by amphiphile self-association plays an important role in various chemical and physical processes, including liquid-liquid extraction (LLE). Small angle scattering techniques are an essential means of probing this structure, but, due to non-uniqueness of scattering patterns, their interpretation is not trivial. Here, in this study, we decompose small angle x-ray scattering (SAXS) patterns for a range of binary LLE organic phases into two components: Ornstein-Zernike and pre-peak contributions resulting from, respectively, critical concentration fluctuations and packing of the amphiphilic extractant molecules in the nonpolar aliphatic diluent. While we previously applied Ornstein-Zernike scattering models to similar organic phases, including the pre-peak explicitly in the fit allows us to measure weak fluctuations at high extractant concentration and simultaneously obtain information on the position and intensity of the correlation peak related to amphiphile packing. Scattering patterns and partial structure factors calculated from molecular dynamics simulations support this interpretation. We demonstrate how this minimal scattering model describes nanostructuring for a large number of extractant types over their entire extractant/diluent composition ranges, suggesting simple and universal behavior. The straightforward assignment of these structural contributions facilitates the comparison of these features between different classes of extractant molecules and will enable further studies of organic phase aggregation. Applicability to more complex, process-relevant organic phases is illustrated with post-contact organic phases containing significant quantities of extracted lanthanide nitrate salt, which we find are readily described by this model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗