Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SeS”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Secondary-electron imaging of bulk crystalline specimens in an aberration corrected $\text{STEM}$

Atomic-scale electron microscopy traditionally probes thin specimens, with thickness below 100 nm, and its feasibility for bulk samples has not been documented. Here in this study, we explore the practicality of scanning transmission electron microscope (STEM) imaging with secondary electrons (SE), using a silicon-wedge specimen having a maximum thickness of 18 μm. We find that the atomic structure is present in the entire thickness range of the SE images although the background intensity increases moderately with thickness. The consistent intensity of secondary electron (SE) images at atomic positions and the modest increase in intensity in between (background) observed in silicon suggest a limited contribution from SEs generated by backscattered electrons, a conclusion supported by our multislice calculations. We conclude that achieving atomic resolution in SE imaging for bulk specimens is indeed attainable using aberration-corrected STEM and that an aberration-corrected scanning electron microscope (SEM) may have the capacity for atomic-level resolution, holding great promise for future strides in materials research.

25 ENERGY STORAGE↗
DOE OSTI2025

Digital light processing of porous LLZTO scaffolds for Li-garnet solid-state batteries

Li 7 La 3 Zr 2 O 12 (LLZO)-based solid-state electrolytes (SEs) are promising materials for next-generation solid-state batteries. In this work, digital light processing (DLP), an emerging additive manufacturing technology, is employed to produce porous Ta-doped LLZO (LLZTO) scaffolds. The self-standing scaffolds are 100 μm thick and have 40% porosity. The scaffolds demonstrate symmetric cell cycling stability exceeding 1,500 h at 0.1 mA/cm2 current density, with a capacity of 0.1 mAh/cm 2 (1 h for each half cycle). At higher current densities, reversible soft shorts frequently happen, while immediate hard shorts are prevented due to Li dendrite growth being hindered by the tortuous pore network. In addition to the cycling stability, the phase stability of LLZTO is investigated during the post-printing thermal process for printing resin removal. We discovered that the LLZTO partially decomposes into Li 2 Zr 2 O 7 and other impurity phases from 400°C to 800°C, but the pure LLZTO phase is restored upon the completion of resin removal beyond 800°C.

Li-metal anode↗
DOE OSTI2025

Recent Improvements to the NWChem COSMO Module

This work presents recent developments in the COSMO solvation model implementation in NWChem. A new cavity construction approach, based on the solvent-excluding surface (SES) and utilizing the well-established GEPOL algorithm, has been introduced. Additionally, a straightforward procedure to merge surface segments that are too close─often a source of numerical artifacts─has been implemented. The available methods for correcting outlying surface charges have also been reviewed and improved. To validate the new implementation, we computed dielectric solvation energies for a chemically diverse set of approximately 100 molecules, including neutral species, small ions, and common ionic liquid components. Results were compared to those from GAMESS using the double-cavity method as a reference. Although the double-cavity approach can be regarded as more accurate, the simpler correction schemes available in NWChem─based on scaling factors or Lagrange multipliers─can achieve excellent agreement if the potential is also properly corrected, with mean unsigned deviations of around 0.15 kcal/mol. Predictions of typical vapor–liquid and liquid–liquid equilibria using a COSMO-SAC variant based on NWChem also showed very promising results.

Cavities↗
DOE OSTI2023

Highly Conductive Iodine and Fluorine Dual-Doped Argyrodite Solid Electrolyte for Lithium Metal Batteries

Sulfide-type argyrodite solid electrolytes (SEs) with halide doping have attracted serious interest. Although other halides such as Cl and Br have been found to enhance Li-ion transport in argyrodites, the direct synthesis of efficiently conductive Li 6 PS 5 I without postprocessing has rarely been investigated. In this work, we report the one-step synthesis of highly conductive Li 6 PS 5 I with an impressive ionic conductivity of 2.5 × 10 –4 S cm –1 at room temperature through a solvent-based method. Moreover, by introducing F – to partially replace I – , hybrid-doped argyrodites Li 6 PS 5 F x I 1–x (x = 0.25, 0.5, 0.75) have been synthesized. Li 6 PS 5 F 0.25 I 0.75 achieves the highest conductivity of 3.5 × 10 –4 S cm –1 due to the energetic preference for anion-disordering among F – /I – and S 2– , which facilitate faster Li transport as supported by density functional theory (DFT) calculations. With higher F content in argyrodites, Li 6 PS 5 F 0.75 I 0.25 displays the best electrochemical stability toward Li metal, as evidenced by long-term stable cycling in Li symmetric cells up to 1100 h. Solid-state Li metal batteries with an active cathode of Li 4 Ti 5 O 12 (LTO) display an initial specific capacity of 140 mAh g –1 and remain at 105 mAh g –1 after 200 cycles, suggesting great battery cycling performance. Furthermore, this research has developed new compositions in the argyrodite SE family which could lead to advancements in the development of solid-state Li metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Lithium Deposition-Induced Fracture of Carbon Nanotubes and Its Implication to Solid-State Batteries

The increasing demand for safe and dense energy storage has shifted research focus from liquid electrolyte-based Li-ion batteries toward solid-state batteries (SSBs). However, the application of SSBs is impeded by uncontrollable Li dendrite growth and short circuiting, the mechanism of which remains elusive. Herein, we conceptualize a scheme to visualize Li deposition in the confined space inside carbon nanotubes (CNTs) to mimic Li deposition dynamics inside solid electrolyte (SE) cracks, where the high-strength CNT walls mimic the mechanically strong SEs. We observed that the deposited Li propagates as a creeping solid in the CNTs, presenting an effective pathway for stress relaxation. When the stress-relaxation pathway is blocked, the Li deposition-induced stress reaches the gigapascal level and causes CNT fracture. Mechanics analysis suggests that interfacial lithiophilicity critically governs Li deposition dynamics and stress relaxation. Furthermore, our study offers critical strategies for suppressing Li dendritic growth and constructing high-energy-density, electrochemically and mechanically robust SSBs.

25 ENERGY STORAGE↗
DOE OSTI2024

Eliminating chemo-mechanical degradation of lithium solid-state battery cathodes during >4.5 V cycling using amorphous Nb2O5 coatings

Abstract Lithium solid-state batteries offer improved safety and energy density. However, the limited stability of solid electrolytes (SEs), as well as irreversible structural and chemical changes in the cathode active material, can result in inferior electrochemical performance, particularly during high-voltage cycling (>4.3 V vs Li/Li + ). Therefore, new materials and strategies are needed to stabilize the cathode/SE interface and preserve the cathode material structure during high-voltage cycling. Here, we introduce a thin (~5 nm) conformal coating of amorphous Nb 2 O 5 on single-crystal LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathode particles using rotary-bed atomic layer deposition (ALD). Full cells with Li 4 Ti 5 O 12 anodes and Nb 2 O 5 -coated cathodes demonstrate a higher initial Coulombic efficiency of 91.6% ± 0.5% compared to 82.2% ± 0.3% for the uncoated samples, along with improved rate capability (10x higher accessible capacity at 2C rate) and remarkable capacity retention during extended cycling (99.4% after 500 cycles at 4.7 V vs Li/Li + ). These improvements are associated with reduced cell polarization and interfacial impedance for the coated samples. Post-cycling electron microscopy analysis reveals that the Nb 2 O 5 coating remains intact and prevents the formation of spinel and rock-salt phases, which eliminates intra-particle cracking of the single-crystal cathode material. These findings demonstrate a potential pathway towards stable and high-performance solid-state batteries during high-voltage operation.

Science & Technology - Other Topics↗
DOE OSTI2022

Different interfacial reactivity of lithium metal chloride electrolytes with high voltage cathodes determines solid-state battery performance

A deep understanding of the interaction of the surface of cathode materials with solid electrolytes is crucial to design advanced solid-state batteries (SSBs). This is especially true for the new class of lithium metal chloride (Li-M-Cl) solid electrolytes which are receiving rapidly growing attention due to their very high oxidative stability (>4 V) in combination with good ionic conductivity that can enable long cell cycle life. While Li-M-Cl electrolytes typically contain resource-limited metals (M) such as indium or rare earths, work has focused on substituting M with more abundant elements such as zirconium. Via operando resistance measurements using intermittent current interruption we explore the dynamic evolution of the interphase at the surface of Ni-rich NCM85 or NCM111 cathode particles inside a working SSB with three different Li-(M 1 ,M 2 )-Cl catholytes (Li 3 InCl 6 , Li 2 Sc 1/3 In 1/3 Cl 4 and Li 5/2 Y 1/2 Zr 1/2 Cl 6 ) to reveal the impact of the cationic metal substitution on the interfacial chemistry. We show that the metal plays a critical role in determining high voltage stability, contrary to prior assumptions. Using a combination of cyclic voltammetry and ultraviolet photoelectron spectroscopy measurements of the electronic band structure to assess oxidative stability; coupled with DFT calculations and ToF-SIMS to evaluate products formed at the interface at different upper cutoff potentials and degrees of delithiation, we are able to differentiate between electrochemical and chemical degradation. We find that Li 2 Sc 1/3 In 1/3 Cl 4 yields the highest (and Li 3 InCl 6 the lowest) stability against electrochemical oxidation, while Li 5/2 Y 1/2 Zr 1/2 Cl 6 undergoes a detrimental chemical reaction with oxygen released from Ni-rich NCM85 at high potentials, resulting in fast capacity fading. Altogether, our work establishes a platform for the metrics and an approach that can be utilized to efficiently evaluate the stability of new halide SEs in SSB cells.

25 ENERGY STORAGE↗
DOE OSTI2020

Sub-system quantum dynamics using coupled cluster downfolding techniques

In this paper, we discuss extending the sub-system embedding sub-algebra coupled cluster (SES-CC) formalism and the double unitary coupled cluster (DUCC) ansatz to the time domain. As we demonstrated in earlier studies, it is possible, using these formalisms, to calculate the energy of the entire system as an eigenvalue of downfolded/effective Hamiltonian in the active space, that is identifiable with the sub-system of the composite system. In these studies, we demonstrated that downfolded Hamiltonians integrate out Fermionic degrees of freedom that do not correspond to the physics encapsulated by the active space. We extend these results to the time-dependent Schrödinger equation, showing that a similar construct is possible to partition a system into a sub-system that varies slowly in time and a remaining subsystem that corresponds to fast oscillations. This time dependent formalism allows coupled cluster quantum dynamics to be extended to larger systems and for the formulation of novel quantum algorithms based on the quantum Lanczos approach, which have recently been considered in the literature.

coupled cluster, Electron correlation, quantum dyn↗
DOE OSTI2025

A Diffuse-Interface Model for Predicting the Evolution of Metallic Negative Electrodes and Interfacial Voids in Solid-State Batteries with Homogeneous and Polycrystalline Solid Electrolyte Separators

Here, this paper presents a novel diffuse-interface electrochemical model that simultaneously simulates the evolution of the metallic negative electrode and interfacial voids during the stripping and plating processes in solid-state batteries. The utility and validity of this model are demonstrated for the first time on a cell with a sodium (Na) negative electrode and a Na-β″-alumina ceramic solid electrolyte (SE) separator. Three examples are simulated. First, stripping and plating with a perfect electrode/electrolyte interface; second, stripping and plating with a single interfacial void at the electrode/electrolyte interface; third, stripping with multiple interfacial voids. Both homogeneous and polycrystalline SEs with low and high-conductivity grain boundaries (GBs) are considered for all three examples. Heterogeneous GB conductivity marginally impacts the thickness of the Na electrode in cases with a perfect electrode/electrolyte interface. Moreover, it results in local changes to void growth due to the interactions between the void edge and the GBs. The void growth rate is a linear function of the flux of Na atoms at the void edge, which in turn depends on the applied current density. We also show that the void coalescence rate increases with applied current density and can be marginally influenced by GB conductivity.

batteries↗
DOE OSTI2025

Observed and Imputed Volumetric Soil Water Content Timeseries for the New Mexico Elevation Gradient

Reliable soil water content (SWC) data are essential for understanding dryland ecosystem dynamics, but high-frequency SWC sensors often fail, creating gaps in critical datasets. To address this, we developed a Bayesian mixture model that imputes missing SWC using both linear interpolation and an ecosystem water balance model (SOILWAT2), tested across six AmeriFlux eddy covariance tower sites in the New Mexico Elevation Gradient, demonstrating its effectiveness in reconstructing SWC patterns while providing insights into the factors driving SWC variability. Daily volumetric soil water content (SWC) data are provided as csv-formatted spreadsheets for the six AmeriFlux sites (US-Seg, US-Ses, US-Wjs, US-Mpi, US-Vcp, and US-Vcs). For each site there is an observed SWC file (site_SWC_gapfill.csv) and a file that contains imputed SWC (imputed_SWC_site.csv). The observed SWC files contain temperature corrected sensor values, tower precipitation data, as well as outputs from SOILWAT2 simulations that were used to impute SWC. The imputed files contain the original observed SWC values and the imputed missing SWC values. When SWC was missing from the original data, the missing value was imputed based on the Bayesian imputation mixture model. The posterior mean of all imputed values is reported as "mean_X". When the observed SWC was NOT missing, mean_X = observed SWC value (original data). The standard deviation, 2.5th percentile and the 97.5th percentile for the imputed values are also reported in the imputed files. There are readme text files for each file type explaining the contents of each column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Design and test of a special magnet for the Mu2e experiment

The Mu2e experiment (Muon-to-Electron Conversion Experiment) is looking for evidence that a muon can change into an electron without the emission of neutrinos. The experiment will search for the charged-lepton flavor violating (CLFV) neutrino-less conversion of a negative muon into an electron in the field of a nucleus. The conversion process results in a monochromatic electron with an energy of 104.97 MeV, slightly below the muon rest mass. The goal of the experiment is to improve the previous upper limit by four orders of magnitude and reach a SES (single event sensitivity) of 3 × 10 -17 on the conversion rate, a 90% CL of 8 × 10 -17 , and a 5σ discovery reach at 2 × 10 -16 . The experiment will use an intense pulsed negative muon beam, which is essential to reducing backgrounds. The other essential element is a sophisticated magnetic system composed of three consecutive solenoids (figure 1) that form the muon beam.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Evaluating the Accuracy of Machine Learning Forecasts

To improve the accuracy of forecasting in machine learning, we must investigate multiple machine learning models and see how accurately they can predict values after training. We used seven machine learning models to try and get more accurate predictions. The models that were used were ARIMA, SES, MLP, CART, LightGBM, and XGBoost. We used a processed dataset from a Terminal at LAX that had the number of people traveling through terminal X every hour in March from 2015-2019. We trained our models with the dates March 6 - March 19 to predict the value for March 20th and the hours 6:00 am to 6:00 pm since those are the most popular traveling hours. By using the different models, we had varying results of accuracy when estimating the amount of people traveling through terminal X on March 20th. We know that machine learning models are helpful for forecasting and by seeing how accurately these models can predict, we can see how forecasting can be helpful for other issues. Using these methods, airports can use forecasting to predict the amount of people coming in and out and can use these predictions to prepare their resource management, operational efficiency, and overall passenger experience.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Time-Evolved Hetero-Alkali Interphases Enable Long-Life Sulfide-Based Anode-Free Solid-State Batteries

Sulfide-based anode-free solid-state batteries (AFSSBs) offer compelling advantages in terms of energy density and safety, yet their practical implementation is severely hindered by undesirable interfacial reactions between sulfide solid electrolytes (SEs) and freshly plated lithium (Li), as well as non-uniform Li plating/stripping behavior. Herein, an effective interfacial stabilization strategy by incorporating sodium bis(fluorosulfonyl)imide (NaFSI) additive into the Li5.4PS4.4Cl1.6 (LPSC) is investigated. Unlike conventional Li-based additives that form static passivation layers, NaFSI introduces a transient hetero-alkali chemistry that kinetically governs interphase evolution during fresh Li plating. NaFSI induces a timesequenced interphase evolution: an initial NaF/LiF-rich layer that suppresses early sulfide reduction, followed by a LiF/Li3N-rich layer that optimizes Li⁺ transport during repeated anode-free cycling. This evolved robust and fast ion conducting layer mitigates interfacial impedance growth, enhances Li + transport kinetics, and suppresses localized Li growth and filamentary shorting. As a result, the anode-free full cell with NaFSI modified LPSC as the interlayer exhibits an excellent cycling stability over 500 cycles at 0.2 C with a capacity retention of 77.6%, whereas the cell with bare LPSC suffers from rapid capacity decay after 100 cycles, retaining only 32.1% of its initial capacity. This work establishes dynamic heteroalkali additive chemistry as a general strategy to kinetically program solid-solid interphases, guiding the interface design in anode-free solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

The association between neighborhood obesogenic factors and prostate cancer risk and mortality: the Southern Community Cohort Study

Background: Prostate cancer is one of the leading causes of cancer-related mortality among men in the United States. We examined the role of neighborhood obesogenic attributes on prostate cancer risk and mortality in the Southern Community Cohort Study (SCCS). Methods: From the total of 34,166 SCCS male participants, 28,356 were included in the analysis. We assessed the relationship between neighborhood obesogenic factors [neighborhood socioeconomic status (nSES) and neighborhood obesogenic environment indices including the restaurant environment index, the retail food environment index, parks, recreational facilities, and businesses] and prostate cancer risk and mortality by controlling for individual-level factors using a multivariable Cox proportional hazards model. We further stratified prostate cancer risk analysis by race and body mass index (BMI). Results: Median follow-up time was 133 months [interquartile range (IQR): 103, 152], and the mean age was 51.62 (SD: ± 8.42) years. There were 1,524 (5.37%) prostate cancer diagnoses and 98 (6.43%) prostate cancer deaths during follow-up. Compared to participants residing in the wealthiest quintile, those residing in the poorest quintile had a higher risk of prostate cancer (aHR = 1.32, 95% CI 1.12–1.57, p = 0.001), particularly among non-obese men with a BMI < 30 (aHR = 1.46, 95% CI 1.07–1.98, p = 0.016). The restaurant environment index was associated with a higher prostate cancer risk in overweight (BMI ≥ 25) White men (aHR = 3.37, 95% CI 1.04–10.94, p = 0.043, quintile 1 vs. None). Obese Black individuals without any neighborhood recreational facilities had a 42% higher risk (aHR = 1.42, 95% CI 1.04–1.94, p = 0.026) compared to those with any access. Compared to residents in the wealthiest quintile and most walkable area, those residing within the poorest quintile (aHR = 3.43, 95% CI 1.54–7.64, p = 0.003) or the least walkable area (aHR = 3.45, 95% CI 1.22–9.78, p = 0.020) had a higher risk of prostate cancer death. Conclusion: Living in a lower-nSES area was associated with a higher prostate cancer risk, particularly among Black men. Restaurant and retail food environment indices were also associated with a higher prostate cancer risk, with stronger associations within overweight White individuals. Finally, residing in a low-SES neighborhood or the least walkable areas were associated with a higher risk of prostate cancer mortality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

SCIMON: Scientific Inspiration Machines Optimized for Novelty

We explore and enhance the ability of neu- ral language models to generate novel scien- tific directions grounded in literature. Work on literature-based hypothesis generation has traditionally focused on binary link prediction— severely limiting the expressivity of hypothe- ses. This line of work also does not focus on optimizing novelty. We take a dramatic depar- ture with a novel setting in which models use as input background contexts (e.g., problems, experimental settings, goals), and output natu- ral language ideas grounded in literature. We present SCIMON, a modeling framework that uses retrieval of “inspirations” from past scien- tific papers, and explicitly optimizes for novelty by iteratively comparing to prior papers and up- dating idea suggestions until sufficient novelty is achieved. Comprehensive evaluations reveal that GPT-4 tends to generate ideas with over- all low technical depth and novelty, while our methods partially mitigate this issue. Our work represents a first step toward evaluating and developing language models that generate new ideas derived from the scientific literature.

Ji, Heng↗