Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase stability”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Novel Superionic Na Conductors for Solid‐State Na batteries: Materials and Interface

This project aims to advance the fundamental understanding on how the cation or anion doping influence on the crystal structure, phase stability and ion transport of novel Na superionic conductors. Specifically, we focus on the Na 3-x A x SbS 4-y X y family solid electrolytes (A: cation dopant, X: anion dopant), a class of promising materials for all-solid-state sodium batteries owing to their high ionic conductivity and compositional tunability. By integrating materials synthesis, neutron scattering, and electrochemical characterization and theoretical simulations, this project sees to elucidate 1) how does the synthesis approaches (solid-state reaction, low-temperature reaction) on the phase purity and doping process; 2) how chemical doping (A: cation dopants such as Sn; X: anion dopants such as Se, F) modulate lattice dynamics and Na+ migration as well as interface stability toward Na metal. Specifically, this project’s proposed technology involves three objectives: (1) synthesize novel Na-ion conductors (anion or cation doped Na 3 SbS 4 ) with high phase purity and examine their structural features (e.g., phase stability); (2) understand Na-ion conductive properties of new Na-ion conductors through a combination of experimental characterizations and theoretical simulations; (3) investigate the interface stability of novel Na superionic conductors with electrode materials in solid-state Na batteries. These objectives will be achieved through combined experimental synthesis/characterizations and computational approaches. The obtained knowledge will provide a comprehensive synthesis-structure-property correlation, enabling the rational design of novel solid-state conductors with excellent conductivity and stability.

25 ENERGY STORAGE↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2025

Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 ⁢O 3

Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 ⁢O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 ⁢O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 ⁢O 3 and 4.94 eV in 𝛽−Ga 2 ⁢O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4⁢𝑓 multiplets and 4⁢𝑓 −4⁢𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µ⁢m, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.

3-dimensional systems↗
DOE OSTI2026

Combustion synthesis of nanoscale magnesium borides with improved hydrogen uptake and release

The overarching goal of the reported project was to develop reversible metal hydrides that can be dehydrogenated and recharged at relatively low temperatures and pressures. The project complemented the research on reversible metal hydrides for hydrogen storage conducted by the Hydrogen Materials – Advanced Research Consortium (HyMARC). The students involved were trained at both the University of Texas at El Paso (UTEP) and Sandia National Laboratories (SNL). The research objectives were to fabricate nanoscale magnesium boride powders with different B/Mg atomic ratios and to investigate the effects of B/Mg ratio and particle size on the hydrogenation properties of the obtained materials. Magnesium diboride (MgB 2 ) and tetraboride (MgB 4 ) were synthesized by using combustion synthesis and by heating in a tube furnace. The latter method resulted in materials with significantly reduced impurity levels. An effective method for the removal of magnesium oxide impurities was identified. High-energy ball milling produced submicron MgB 2 and MgB 4 powders. However, hydrogenation experiments at pressures up to 700 bar have not shown any uptake of hydrogen. Thermogravimetric analysis of the decomposition of the obtained MgB 2 and MgB 4 demonstrated higher thermodynamic stability of the latter. Thus, the boron-rich nature and phase stability of magnesium borides with high B/Mg ratios make them poor candidates for direct hydrogenation. However, the synthesized MgB 4 is a promising energetic additive to solid fuels for high-speed air-breathing propulsion as it offers a lower oxidation onset temperature and greater energy density than MgB 2 .

08 HYDROGEN↗
DOE OSTI2025

Mesoscale fractal whey protein particles derived from microscale linear-shaped protein assemblies (Part 1): Manufacturing method and particle characteristics

Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.

Antisolvent precipitation↗
DOE OSTI2025

Phonon density of states in uranium Laves phases UNi 2 and UCo 2

Recent computational work evaluated the phase stability of UCo 2 and UNi 2 in the C14, C15, and C36 Laves phases and reported their computationally determined phonon density of states. To corroborate these computational findings, we report the experimentally determined phonon density of states for UCo 2 and UNi 2 using inelastic neutron scattering at the wide angular-range chopper spectrometer at the Spallation Neutron Source at Oak Ridge National Laboratory. We found excellent agreement between experimental and computational results for UNi 2 and noted that some differences were observed for UCo 2 . Furthermore, by comparing the phonon density of states of UCo 2 with that of other U-bearing Laves phases, we found that UCo 2 showed distinct vibrational properties. From the one-phonon density of states, we calculated the lattice contribution to the specific heat of UNi 2 and UCo 2 . Finally, we bound the electronic contribution to the specific heat in UNi 2 .

Actinides↗
DOE OSTI2026

High-throughput validation of phase formability and simulation accuracy of Cantor alloys

High-throughput methods enable accelerated discovery of novel materials in complex systems such as high-entropy alloys, which exhibit intricate phase stability across vast compositional spaces. Computational approaches, including Density Functional Theory (DFT) and calculation of phase diagrams (CALPHAD), facilitate screening of phase formability as a function of composition and temperature. However, the integration of computational predictions with experimental validation remains challenging in high-throughput studies. In this work, we introduce a quantitative confidence metric to assess the agreement between predictions and experimental observations, providing a quantitative measure of the confidence of machine learning models trained on either DFT or CALPHAD input in accounting for experimental evidence. The experimental dataset was generated via high-throughput in-situ synchrotron X-ray diffraction on compositionally varied FeNiMnCr alloy libraries, heated from room temperature to ~1000 °C. Agreement between the observed and predicted phases was evaluated using either temperature-independent phase classification or a model that incorporates a temperature-dependent probability of phase formation. This integrated approach demonstrates where strong overall agreement between computation and experiment exists, while also identifying key discrepancies, particularly in FCC/BCC predictions at Mn-rich regions to inform future model refinement.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Back to functional hydrides: Effects of neutron-irradiated microstructure on hydrogen retention in yttrium hydride

Functional hydrides are promising candidates for advanced nuclear reactors, particularly in portable or transportable applications, due to their high hydrogen-retention capabilities, enabling efficient neutron moderation, and compact reactor design. However, hydrogen mobility in hydrides at elevated irradiation temperatures poses significant technological challenges, necessitating a comprehensive understanding of their irradiation behavior. Furthermore, this study investigated the microstructural and chemical stability of neutron-irradiated yttrium hydrides to assess their hydrogen-retention capacity. A targeted literature review was also conducted to contextualize neutron-irradiation effects on functional hydrides with regards to structural stability and hydrogen retention. Experimental characterizations revealed that, at high temperatures, irradiated hydrides retained their phase stability, which was likely enhanced by irradiation-induced microstructure evolution. Notably, an amorphous yttrium and oxygen -rich surface layer was present at the free surface of the hydride. Its thickness decreased while a continuous crystalline Y-O-rich layer was formed with increasing neutron damage. Additionally, the number density of dislocation loops and cavities generally increased as a function of neutron dose. First-principles calculations of hydrogen behavior within yttrium vacancy clusters in yttrium hydrides demonstrated vacancy-size-dependent hydrogen stability and configuration, highlighting the role of vacancy geometry in regulating hydrogen retention. Thus, the presence of irradiation-induced dislocation loops and cavities were hypothesized to improve hydrogen retention. Collectively, these findings advance the understanding of hydride behavior under neutron irradiation as well as their technological readiness for portable or transportable nuclear reactors.

Defect clusters↗
DOE OSTI2025

Real-Time High-Accuracy Digital Wireless Time, Frequency, and Phase Calibration for Coherent Distributed Antenna Arrays

his work presents a fully-digital high-accuracy real-time calibration procedure for frequency and time alignment of open-loop wirelessly coordinated coherent distributed antenna array (CDA) modems, enabling radio frequency (RF) phase coherence of spatially separated commercial off-the-shelf (COTS) software-defined radios (SDRs) without cables or external references such as the global navigation satellite system (GNSS). Building on previous work using high-accuracy spectrally-sparse time of arrival (ToA) waveforms and a multistep ToA refinement process, a high-accuracy two-way time transfer (TWTT)-based time–frequency coordination approach is demonstrated. Due to the two-way nature of the high-accuracy TWTT approach, the time and frequency estimates are Doppler and multipath tolerant, so long as the channel is reciprocal over the synchronization epoch. This technique is experimentally verified using COTS SDRs in a lab environment in static and dynamic scenarios and with significant multipath scatterers. Time, frequency, and phase stability were evaluated by beamforming over coaxial cables to an oscilloscope which achieved time and phase precisions of ~60– 70 ps , with median coherent gains above 99% using optimized coordination parameters, and a beamforming frequency root-mean-square error (RMSE) of 3.73 ppb in a dynamic scenario. Finally, experiments were conducted to compare the performance of this technique with previous works using an analog continuous-wave two-tone (CWTT) frequency reference technique in both static and dynamic settings.

Clock synchronization↗
DOE OSTI2026

Process-driven roadmap for depositing super duplex stainless steel via wire Arc additive manufacturing

Here, this study systematically investigates the effects of shielding gas, bead spacing, weld mode, and travel speed on the phase balance, porosity, and hardness of wire arc additively manufactured (WAAM) super duplex stainless steel (ER2594). Robotic WAAM was employed to fabricate multilayer walls under systematically varied process conditions, followed by phase transformation simulations, X-ray computed tomography (XCT), electron backscatter diffraction (EBSD), and microhardness evaluation. Thermodynamic simulations predicted rapid cooling of AM process can suppress the potential formation of deleterious precipitates which was later validated via cross-sectional microstructure analyses of printed samples. XCT revealed porosity levels below 0.2% for all deposits, with reduced travel speed significantly lowering defect volume. Microstructural analyses revealed the evolution of various austenite precipitates, including grain boundary austenite (GBA), Widmanstätten austenite (WA), and intergranular austenite (IGA), sequentially upon cooling of the ferrite phase. Among all process parameters, weld transfer mode exhibited the strongest influence on phase balance; pulsed mode promoted higher ferrite retention (~ 36%) compared to RapidX mode. No consistent relationship between stabilized phase fraction and captured microhardness was observed. This work provides critical insights for optimizing WAAM parameters to control phase balance and mechanical performance, which is essential for producing high-integrity super duplex stainless-steel components for nuclear and marine applications.

Grain orientation↗
DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Thermally driven surface phase separation in intermetallic alloys

Intermetallic compounds are widely recognized for their high-temperature phase stability and resistance to composition and structural changes. However, we reveal athermally activated bulk-to-surface mass exchange mechanism that drives surface phase separation, resulting in the formation of surface precipitates with distinct composition andstructure from the bulk matrix. Using the archetypal β-NiAl system, we show that asymmetries in vacancy formation energies between Ni and Al atoms induce preferential Ni segregation to the surface, forming Ni-rich γ'-Ni 3 Al precipitates. By integrating in-situ electron microscopy, synchrotron X-ray absorption spectroscopy and first-principles computational modeling, we establish a direct mechanistic connection between bulk thermal defect dynamics, surface compositional evolution, and phase segregation behavior. This bulk-surface coupling mechanism can be a driver of surface phase separation in multicomponent alloys under thermal stress. In conclusion, these results refine the thermodynamic boundaries of intermetallic stability and provide insights into managingthe performance and durability of intermetallic alloys for demanding high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Epitaxial growth and physical properties of Bi 2 Ru 2 O 7 thin films on YSZ(111) substrates

We systematically investigated the growth of Bi 2 Ru 2 O 7 thin films on a Y-stabilized ZrO 2 (111) substrate using pulsed laser deposition by mapping the influence of growth temperature and oxygen partial pressure on phase stability, lattice parameters, and cation ratio. The results show that the epitaxial stabilization requires a minimum growth temperature, which is rather insensitive to the pressure. Meanwhile, the Bi:Ru ratio decreases when increasing growth temperature or decreasing pressure. By constructing the temperature–pressure phase diagram, an optimal growth window within the epitaxial phase was established. On the other hand, the electrical resistivity remains at a similar level within the epitaxial phase with only subtle changes to the temperature dependence, indicative of the robustness of the conductivity against composition variation. Our study provides a foundation for future investigations on thin films and heterostructures that utilize Bi 2 Ru 2 O 7 .

Crystallography↗
DOE OSTI2025

Stability, growth, and doping of In 2 (Si, Ge) 2 O 7 : Promising n -type wide-bandgap semiconductors

In this paper, we investigate, computationally and experimentally, the phase stability, electronic structure properties, and the propensity for n-type doping of In 2 X 2 O 7 (X = Si, Ge) ternary oxides. This family of materials contains promising novel wide-gap semiconductors based on their estimated high n-type Baliga figures of merit and acceptable thermal conductivity for power electronics applications. Here, we predict that both In 2 Si 2 O 7 and In 2 Ge 2 O 7 are n-type dopable, with Zr providing between 10 16 and above 10 21 cm −3 net donor concentrations under O-poor conditions, depending on the chemistry, structure (ground-state thortveitite or high-pressure pyrochlore), and synthesis temperature. To verify our predictions, we synthesize Zr-doped In 2 Ge 2 O 7 in the thortveitite structure and measure its electrical properties. Initial thin-film growth and annealing lead to polycrystalline thin films with bandgaps over 4 eV and confirm Zr doping predictions by achieving electron concentrations at 10 14 –10 16 cm −3 even under O-rich conditions. While future epitaxial growth development is still needed, this study establishes In 2 X 2 O 7 as promising n-type wide-gap semiconductors for power electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Enhanced stability and ordering of 2D MnBi2Se2Te2

This dataset includes density functional theory and quantum Monte Carlo calculations on MnBi(Se(1-x)Te(x))4 system for its phase stability and defect formation energies. The data supports publication "Increased Defect Resistance and Ordering in MnBi2(Se(1–x)Te(x))4 via Accurate Diffusion Monte Carlo" (https://doi.org/10.1021/acs.jpcc.5c04299).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bromine Incorporation Affects Phase Transformations and Thermal Stability of Lead Halide Perovskites

Mixed-cation and mixed-halide lead halide perovskites show great potential for their application in photovoltaics. Many of the high-performance compositions are made of cesium, formamidinium, lead, iodine, and bromine. However, incorporating bromine in iodine-rich compositions and its effects on the thermal stability of the perovskite structure has not been thoroughly studied. In this work, we study how replacing iodine with bromine in the state-of-the-art Cs 0.17 FA 0.83 PbI 3 perovskite composition leads to different dynamics in the phase transformations as a function of temperature. Through a combination of structural characterization, cathodoluminescence mapping, X-ray photoelectron spectroscopy, and first-principles calculations, we reveal that the incorporation of bromine reduces the thermodynamic phase stability of the films and shifts the products of phase transformations. Our results suggest that bromine-driven vacancy formation during high temperature exposure leads to irreversible transformations into PbI 2 , whereas materials with only iodine go through transformations into hexagonal polytypes, such as the 4H-FAPbI 3 phase. This work sheds light on the structural impacts of adding bromine on thermodynamic phase stability and provides new insights into the importance of understanding the complexity of phase transformations and secondary phases in mixed-cation and mixed-halide systems.

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

The Ponderomotive Effects of Narrow-band, Superconducting Resonators in Open and Closed Loop

In this work, we present measurements of ponderomotive instabilities in narrow-band, coaxial resonators in open and closed control loop systems. We show an analytical scheme we will use in future studies to examine the dependency of open loop stability on mechanical parameters. Analytical and simulation models are used to predict the onset of the oscillatory instability in half wave resonators with active disturbance rejection control for amplitude and phase stabilization. We demonstrate that for high amplitude controller bandwidths, ponderomotive oscillations can couple with the controller frequency response via higher harmonics and lower the threshold for the oscillatory instability. We reaffirm the superiority of in-phase/quadrature $(I/Q)$ component control in regards to preventing the oscillatory instability, and show how the phase controller can amplify the crosstalk of disturbances in amplitude and phase control. In cases with large disturbances, this effect can lead to instabilities not present in $I/Q$ control.

Brown, Jacob [Michigan State U.; Fermilab] (ORCID:↗
DOE OSTI2026

Computational discovery of a tetragonal MBene phase: diamond-shaped boron lattices, viable MAB precursors, and selective CO 2 reduction

MBenes are a class of two-dimensional transition-metal borides derived from layered MAB phases. Here, we use first-principles calculations to predict a new tetragonal MB (1 : 1) MBene phase featuring a diamond-shaped B–B lattice and systematically assess its phase stability, synthetic accessibility, and catalytic performance towards the CO 2 reduction reaction (CO 2 RR). Four tetragonal configurations are benchmarked against known hexagonal and orthorhombic phases, and a multi-tier screening identifies six robust members (CrB, FeB, MoB, WB, IrB, and PtB). Exploration of the I4/mmm MAB space reveals 11 viable precursors, including the experimentally synthesized Ir 2 ZnB 2 , which validates our screening approach. CO 2 RR free-energy diagrams show that CrB favors CH 3 OH while FeB, MoB, WB, IrB, and PtB preferentially yield HCOOH. The corresponding limiting potentials are competitive with, and in some cases superior to, those reported for state-of-the-art orthorhombic and hexagonal MBenes. These findings expand the MBene landscape by establishing a tetragonal phase with structural, synthetic, and catalytic promise.

Lu, Linguo [University of Puerto Rico, San Juan, P↗