Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic simulations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Irradiation-Induced Structural Disorder and Its Influence on the Mechanical Response of Polycrystalline MoS2

Molybdenum disulfide (MoS2) thin films are widely used as dry-film lubricants and protective coatings in aerospace and other radiation-exposed environments. Conventional synthesis routes produce polycrystalline films whose grain boundaries and other native defects cause their mechanical and tribological behavior to differ substantially from that of ideal single crystals. Under irradiation, these films progressively evolve from polycrystalline structures, composed of layered MoS2 grains, into highly disordered and eventually amorphous structures, altering both their tribological performance and mechanical integrity. Here, we employ reactive atomistic simulations to investigate irradiation-driven structural evolution in bulk polycrystalline MoS2. Using controlled primary knock-on atom (PKA) events, we characterize the progressive transition from a polycrystalline microstructure to an amorphous network by tracking defect accumulation and structural disorder. We then establish how this transition modifies the dominant deformation mechanisms and the temperature-dependent tensile response. Specifically, irradiation suppresses interlayer sliding and delamination, mechanisms which facilitate the deformation of the pristine polycrystal, resulting in defect-induced hardening. Broadly, our results establish direct process–structure–property relationships linking irradiation-induced defect accumulation, microstructural evolution, deformation mechanisms, and mechanical behavior, providing an atomistic framework for understanding the structural integrity and long-term reliability of irradiated MoS2 coatings.

Moore, Daniel [Sandia National Laboratories (SNL)]↗
DOE OSTI2024

Stabilizing Nickel‐Rich Cathodes in Aqueous Process through Nanocellulose as Water Barrier

Nickel-rich LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC 811) cathode offers high voltage and high specific capacity, making it promising for high energy density batteries. However, large-scale manufacturing of aqueous-processed NMC 811 electrodes remains challenging due to proton exchange causing material decomposition and capacity loss. This work addresses this issue by constructing an in situ nanocellulose protective layer for NMC 811 particles via electrostatic interactions during the slurry preparation. For the first time, the interatomic spacing between inter-chains of nanocellulose is measured through wide-angle X-ray scattering and demonstrate the ability to effectively confine interlayer water using atomistic simulations. Moreover, this nanocellulose coverage simultaneously minimizes Li + surface segregation and mitigates water infiltration. Owing to less material decomposition during the aqueous processing, nanocellulose-protected NMC electrodes exhibit higher initial coulombic efficiency (83% vs 62% at 0.1C) and capacity (133 vs 59 mAh g −1 at 6C) than unprotected electrodes. Additionally, optimized aqueous-processed NMC electrodes offer comparable or even superior electrochemical properties compared to the electrodes fabricated using the conventional toxic organic solvent, N-methyl-2-pyrrolidone. Consequently, the developed approach enables affordable, sustainable aqueous processing for Nickel-rich NMC 811 cathodes with excellent electrochemical performances.

25 ENERGY STORAGE↗
DOE OSTI2025

On the Structure–Property Relationship of Semi‐Coherent FeCr 2 O 4 /Cr 2 O 3 Spinel/Corundum Interfaces

Abstract Oxide heterointerfaces are extremely common in both natural and artificial composite structures, including corroded structural materials. Often, key properties such as segregation and atomic transport are dictated by the structure of these interfaces. However, despite this critical link, very few heterointerfaces have been studied in any detail at the atomic scale. Here, one important oxide heterointerface is examined, between spinel and corundum, using the chemical system FeCr 2 O 4 /Cr 2 O 3 as a representative and technologically important case. Using atomistic simulation techniques, it is found that the structure, particularly the local chemistry, of the interface depends on the crystal chemistry at the interface. This atomic and chemical structure further impacts important properties such as defect segregation and mass transport. It is found that defects can nucleate at some regions of these interfaces and migrate back and forth across the corundum layer, suggesting high atomic mobility that may be important for the evolution of spinel/corundum composite structures in extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Direct Observation of Elusive (DTBM‐SEGPHOS)CuH Monomer Enables Mechanistic Insights Into Hydrocupration, Aggregation, and Dynamics of Alkene Functionalization Catalysis

The bulky diphosphine DTBM-SEGPHOS is widely employed in CuH-catalyzed transformations as it provides remarkably active catalyst systems. The transient (DTBM-SEGPHOS)CuH monomer (LCuH) is the often-invoked active species. However, its instability has prevented spectroscopic characterization and mechanistic elucidation, hindering mechanistic understanding. We report low-temperature NMR spectroscopic characterization of LCuH, enabling quantitative kinetic analysis of the stoichiometric hydrocupration and catalytic hydroboration of cyclopentene, as well as the structural identification of two CuH clusters. LCuH inserts cyclopentene at −43°C, reaffirming its high reactivity toward olefins. LCuH deactivates to form L 2 Cu 3 H 3 and L 2 Cu 4 H 4 clusters, in which LCuH dimerization initiates aggregation. Kinetic analysis of reactions of unactivated alkenes indicates that competing on-cycle alkene hydrocupration and LCuH dimerization impact performance, as catalyst deactivation and turnover occur on comparable timescales. Structure–activity analysis using atomistic simulations shows that the steric profile of DTBM-SEGPHOS increases the CuH dimerization barrier by ∼7.7 kcal mol−1 compared to that of SEGPHOS, rationalizing the unique ability of DTBM-SEGPHOS to stabilize a reactive monomer for hydrocupration of broader alkene substrates. These findings illustrate the fundamental design principle that steric control of aggregation governs CuH catalyst performance, explaining both the exceptional activity of (DTBM-SEGPHOS)CuH and the limitations imposed by competing deactivation.

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

Unlocking High Capacity and Reversible Alkaline Iron Redox Using Silicate-Sodium Hydroxide Hybrid Electrolytes

Alkaline iron (Fe) batteries are attractive due to the high abundance, low cost, and multiple valent states of Fe but show limited columbic efficiency and storage capacity when forming electrochemically inert Fe 3 O 4 on discharging and parasitic H 2 on charging. Herein, sodium silicate is found to promote Fe(OH) 2 /FeOOH against Fe(OH) 2 /Fe 3 O 4 conversions. Electrochemical experiments, operando X-ray characterization, and atomistic simulations reveal that improved Fe(OH) 2 /FeOOH conversion originates from (i) strong interaction between sodium silicate and iron oxide and (ii) silicate-induced strengthening of hydrogen-bond networks in electrolytes that inhibits water transport. Furthermore, the silicate additive suppresses hydrogen evolution by impairing energetics of water dissociation and hydroxyl de-sorption on iron surfaces. In conclusion, this new silicate-assisted redox chemistry mitigates H 2 and Fe 3 O 4 formation, improving storage capacity (199 mAh g -1 in half-cells) and coulombic efficiency (94 % after 400 full-cell cycles), paving a path to realizing green battery systems built from earth-abundant materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2024

Ultimate compressive strength and severe plastic deformation of equilibrated single-crystalline copper nanoparticles

Mechanical properties and deformation mechanisms of defect-free copper nanoparticles are investigated by combining experiments with atomistic simulations. The compressive strength of the particles increases with decreasing size and tends to saturate near the theoretical strength in the small-size limit. In this limit, the intrinsic size dependence of the strength is governed by the stochastic nature of dislocation nucleation near the particle surface. The particle deformation process evolves from the initial strain softening to strain hardening as the particle accumulates residual damage. The normalized strength-size relation for Cu is compared with those for Au, Ni, and Pt. The lack of universal behavior among the four FCC metals is discussed. Heavily deformed Cu nanoparticles develop polycrystalline structures and change the lattice orientation from [111] to [110]. The experiments and simulations reveal the twinning mechanism of the lattice rotation leading to the new grain formation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Formation of carbon homonuclear bonds in β-SiC under neutron irradiation at various temperatures and neutron doses

To elucidate radiation defect processes in SiC, Raman spectroscopy was systematically applied to high-purity, polycrystalline β-SiC that was neutron irradiated at a range of temperature and dose conditions. The analysis specifically focused on formation of carbon homonuclear bonds by irradiation; these bonds were indicated by D and G bands and amorphous carbon peaks. Intensity of the carbon peaks relative to SiC peaks significantly decreased in the case of high temperature and/or high neutron dose of 500 °C to 29 displacements per atom (dpa) and about 800 °C to 1.38 and 29 dpa. The absence of carbon bond peaks under those conditions was explained by growth of stoichiometric defect clusters, consistent with previous atomistic simulations on SiC defect stability. The lack of Raman bands associated with carbon clusters under high-temperature and high-dose radiation conditions accounts for the resistance of SiC to phase separation under irradiation. The findings further suggest that material compositions and chemical properties that are inherently resistant to chemical disordering under high-dose radiation conditions are indicative of the long-term durability of ceramic compounds in radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.

Critical benchmarks↗
DOE OSTI2024

A high-throughput workflow to analyze sequence-conformation relationships and explore hydrophobic patterning in disordered peptoids

Understanding how a macromolecule’s primary sequence governs its conformational landscape is crucial for elucidating its function, yet these design principles are still emerging for macromolecules with intrinsic disorder. Herein, we introduce a high-throughput workflow that implements a practical colorimetric conformational assay, introduces a semi-automated sequencing protocol using matrix-assisted laser desorption/ionization and tandem mass spectrometry (MALDI-MS/MS), and develops a generalizable sequence-structure algorithm. Using a model system of 20mer peptidomimetics containing polar glycine and hydrophobic N-butylglycine residues, we identified nine classifications of conformational disorder and isolated 122 unique sequences across varied compositions and conformations. Conformational distributions of three compositionally identical library sequences were corroborated through atomistic simulations and ion mobility spectrometry coupled with liquid chromatography. A data-driven strategy was developed using existing sequence variables and data-derived “motifs” to inform a machine-learning algorithm toward conformation prediction. Here, this multifaceted approach enhances our understanding of sequence-conformation relationships and offers a powerful tool for accelerating the discovery of materials with conformational control.

data-driven analysis↗
DOE OSTI2025

Machine learning enabled discovery of superhard and ultrahard carbon polymorphs

The demand for multifunctional materials has motivated the move from near-equilibrium materials to metastable i.e. out-of-equilibrium phases that can meet several desired target properties. The search for such metastable phases with exotic properties is non-trivial and often serendipitous. Inverse design approaches based on evolutionary search have been powerful tools, but such traditional searches have focused on identifying primarily stable and metastable materials with the lowest enthalpy. The inverse design of materials, with a focus on a desired property such as, for example, hardness is a challenging task because of the expensive computational cost involved in sampling multiple structures. The recent advances in machine learning have brought new powerful AI techniques to the forefront which can potentially revolutionize the inverse design and discovery of materials, especially metastable phases capable of meeting multifunctionality. Here, in this work, we develop and apply an automated reinforcement learning workflow for inverse design that integrates first principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the superhard and ultrahard metastable phases of Carbon. We demonstrate an automatic machine learning based inverse design workflow to map new undiscovered metastable states ranging from near equilibrium to those far-from-equilibrium that satisfy multiple property objectives, specifically bulk moduli, shear moduli and hardness. We create a comprehensive library of carbon stable and metastable phases with varying hardness and subsequently shortlist 10 top performing candidate carbon structures, including two newly reported phases, based on their hardness and characterize their temperature dependent mechanical properties. A neural network model is built using featurization of allotropes of carbon to predict the quasi-harmonic Gibbs free energies. The Gibbs free energies of the top performing phases are analyzed to get an estimate of the experimental synthesizability of these superhard and ultrahard carbon phases. In general, we show using machine learning based inverse design approaches how hitherto inaccessible metastable states can be identified and potentially synthesized to meet the demand for multifunctional materials.

Balasubramanian, Karthik [Univ. of Illinois, Chica↗
DOE OSTI2025

Computational materials assessment of the D/Li-stripping neutron source as a prototypical facility for fusion materials testing

As the US fusion materials community awaits the selection and design of a fusion prototypical neutron source (FPNS), a risk reduction exercise has been conducted to (i) provide an updated materials performance evaluation using state-of-the-art computational materials modeling, (ii) expand on legacy analysis based on pure Fe to other relevant fusion structural materials types, and (iii) ensure that materials response under FPNS operational conditions is consistent with referential fusion reactor conditions. The current paper describes the efforts undertaken to assemble a comprehensive computational methodology that includes neutronics, primary damage calculations, atomistic simulations of displacement cascades, chemical inventory evolution calculations, and a computational thermodynamic analysis of emerging phases during irradiation. Our work extends existing studies in pure Fe to reduced-activation ferritic/martensitic steels, tungsten, silicon carbide, and vanadium alloys. We focus on the single-beam deuteron/lithium-stripping neutron source behind the IFMIF-DONES concept, which we assess against ITER, two DEMO designs, and an ideal pure 14-MeV flux. Our analysis indicates that, within standard uncertainties inherent to the models employed, the DONES concept adequately captures fusion conditions in the four materials analyzed. Our work is intended as a comprehensive irradiation damage analysis of fusion-representative neutron sources, to be used for further neutron source evaluation and fusion facility operation.

Marian, Jaime [Univ. of California, Los Angeles, C↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2025

Interplay between hydrogen, temperature, and character angle on dissociated dislocation energies in Fe–Ni–Cr austenitic stainless steels

Dislocation energy has an important role in the mechanical performance of structural metals. While dislocation energies cannot be fully obtained from continuum theories due to the contribution of the dislocation core, they have been calculated via atomistic simulations in elemental metals. However, constraints on the local atomic environments have prevented the use of such approaches in systems that incorporate alloying or interstitial solutes. In this work, we develop robust molecular dynamics methods to resolve these issues through a geometric construction of dislocation dipoles and the calculation of time-averaged energies. Furthermore, we apply these methods to calculate dislocation energies (including core energies) in an Fe 70 Ni 11 Cr 19 austenitic steel at a variety of character angles, hydrogen concentrations, temperatures, and dipole spacings. The resulting highly converged energies show an excellent agreement with continuum expressions. Overall, hydrogen concentrations up to 1.0 % do not have a significant effect on the elastic parameters and dislocation energy. The methods and insights derived in this work have the potential to facilitate the calculation of dislocation energies in a wide range of systems, and to guide our understanding of hydrogen embrittlement.

Alloyed systems↗
DOE OSTI2024

Hydrogen effects on the deformation and slip localization in a single crystal austenitic stainless steel

Hydrogen is known to embrittle austenitic stainless steels, which are widely used in high-pressure hydrogen storage and delivery systems, but the mechanisms that lead to such material degradation are still being elucidated. The current work investigates the deformation behavior of single crystal austenitic stainless steel 316L through combined uniaxial tensile testing, characterization and atomistic simulations. Thermally precharged hydrogen is shown to increase the critical resolved shear stress (CRSS) without previously reported deviations from Schmid’s law. Molecular dynamics simulations further expose the statistical nature of the hydrogen and vacancy contributions to the CRSS in the presence of alloying. Slip distribution quantification over large in-plane distances (> 1 mm), achieved via atomic force microscopy (AFM), highlights the role of hydrogen increasing the degree of slip localization in both single and multiple slip configurations. The most active slip bands accumulate significantly more deformation in hydrogen precharged specimens, with potential implications for damage nucleation. For $\langle$110$\rangle$ tensile loading, slip localization further enhances the activity of secondary slip, increases the density of geometrically necessary dislocations and leads to a distinct lattice rotation behavior compared to hydrogen-free specimens, as evidenced by electron backscatter diffraction (EBSD) maps. Finally, the results of this study provide a more comprehensive picture of the deformation aspect of hydrogen embrittlement in austenitic stainless steels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructurally informed synchrotron x-ray analysis revealing helium defect transitions in ultrafine grained tungsten

The formation of insoluble gaseous defects in materials due to nuclear transmutation or ion implantation involves the diffusion of impurity atoms to form atomic defect clusters that coalesce into bubbles or cavities and ultimately degrade the material properties. Transmission electron microscopy (TEM) is limited in its ability to resolve sub-nanometer gas clusters whereas X-ray diffraction (XRD) provides information pertaining to local atomic changes. Here, in this study, helium (He) implanted ultrafine grained tungsten is explored through a multimodal defect characterization campaign combining TEM-informed Small Angle X-ray Scattering (SAXS) analysis, XRD lattice parameter measurements, and nanoscale He cluster quantification from a region of reciprocal space accessible via Wide Angle X-ray Scattering (WAXS). Moderate elevated temperature implantations are shown to produce high concentrations of sub-nanoscale He clusters and small, homogeneously distributed cavities, which collectively are linked to lattice expansion and further substantiated through complementary atomistic simulations. Increased implantation temperatures encourage the diffusion of these defects to the grain boundaries (GBs), leading to lattice relaxation and the growth of large GB cavities manifesting as bimodal size distributions in the SAXS analysis. Overall, our results demonstrate the utility of multimodal synchrotron X-ray analysis in bridging the gap between microscale He cavity quantification and atomic-scale defect analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning-enabled multiscale modeling of mechanical deformation of aluminum and Al-SiC nanocomposites

A machine learning-enabled multiscale framework is developed for modeling the mechanical response of both pure metal and nanoparticle-reinforced metal matrix nanocomposites (MMNCs). Using aluminum–silicon carbide (Al-SiC) as an example MMNC, atomistic simulations reveal three distinct deformation mechanisms (i.e., defect-free, dislocation-based, and interface separation) governed by the interfaces between the Al matrix and SiC nanoparticles. As compared with single crystal Al, the lattice undergoes a more abrupt failure once the dislocation network becomes extensive and void nucleation initiates, whereas in Al-SiC, nanoparticle interfaces enable a more gradual progression of damage. These mechanisms are captured through a combined classification-regression neural network surrogate model that bridges atomic-scale insights with continuum-scale finite element analysis. Machine learning-enabled multiscale modeling of pure Al accurately predicted strain localization and confirmed by in-situ scanning electron microscopic tensile testing on perforated Al specimens. This study underscores the promise of integrating physics-informed machine learning with hierarchical modeling to capture the interface dominated phenomena and guide the design of advanced MMNCs.

Al-SiC↗
DOE OSTI2025

Pressure-driven density match nucleates metastable r8 phases from amorphous Si and Ge

The pressure–temperature phase behavior of covalent disordered solids such as amorphous silicon and germanium is complex. Questions remain on possible glass transitions, on polyamorphism via amorphous–amorphous transitions, on connections with liquid–liquid transitions, on structure-behavior relationships, and on their potential as precursor for novel methods for material discovery. Here we demonstrate experimentally the nucleation of a metastable, four-fold coordinated rhombohedral r8 phase from pure amorphous silicon and germanium upon room temperature compression at pressures below 10 GPa. Accompanying theory reveals a strong pressure-driven distortion of the bond angle transforming the starting tetrahedral low-density amorphous network to a distorted four-fold coordinated medium-density state. This state is of lower density than metallic high-density networks, resembles the crystalline r8 phase and initiates its nucleation. Our finding shows that polyamorphism is not the only possible transformation mode for these amorphous solids and that instead nucleation of interesting functional phases at potentially useful pressures is possible. Such novel access modes to metastable structures are critical for future exploitability and could be useful for other tetrahedral materials including carbon, where the related (bc8) post-diamond phase remains elusive. Our observed density match between an amorphous and a metastable crystalline phase clearly allows for a new phase transition pathway, while corresponding theory demonstrates how carefully validated atomistic simulations can guide prediction, discovery and synthesis of novel material structures.

Materials discovery↗