Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic dynamics”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Abstract We present a way to dramatically accelerate Gaussian process models for interatomic force fields based on many-body kernels by mapping both forces and uncertainties onto functions of low-dimensional features. This allows for automated active learning of models combining near-quantum accuracy, built-in uncertainty, and constant cost of evaluation that is comparable to classical analytical models, capable of simulating millions of atoms. Using this approach, we perform large-scale molecular dynamics simulations of the stability of the stanene monolayer. We discover an unusual phase transformation mechanism of 2D stanene, where ripples lead to nucleation of bilayer defects, densification into a disordered multilayer structure, followed by formation of bulk liquid at high temperature or nucleation and growth of the 3D bcc crystal at low temperature. The presented method opens possibilities for rapid development of fast accurate uncertainty-aware models for simulating long-time large-scale dynamics of complex materials.

Chemistry↗
NASA NTRSDate not supplied

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗
DOE OSTI2021

Site-specific interrogation of an ionic chiral fragment during photolysis using an X-ray free-electron laser

Abstract Short-wavelength free-electron lasers with their ultrashort pulses at high intensities have originated new approaches for tracking molecular dynamics from the vista of specific sites. X-ray pump X-ray probe schemes even allow to address individual atomic constituents with a ‘trigger’-event that preludes the subsequent molecular dynamics while being able to selectively probe the evolving structure with a time-delayed second X-ray pulse. Here, we use a linearly polarized X-ray photon to trigger the photolysis of a prototypical chiral molecule, namely trifluoromethyloxirane (C 3 H 3 F 3 O), at the fluorine K-edge at around 700 eV. The created fluorine-containing fragments are then probed by a second, circularly polarized X-ray pulse of higher photon energy in order to investigate the chemically shifted inner-shell electrons of the ionic mother-fragment for their stereochemical sensitivity. We experimentally demonstrate and theoretically support how two-color X-ray pump X-ray probe experiments with polarization control enable XFELs as tools for chiral recognition.

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

Quasiaperiodic grain boundary phases of Σ⁢5 tilt grain boundaries in refractory metals

We report ground-state structures and phase transitions in Σ⁢5⁢[001] tilt grain boundaries (GBs) in body-centered-cubic (bcc) refractory metals Nb, Ta, Mo, and W. Σ⁢5 tilt GBs have been extensively investigated over the past several decades, with their ground-state structure—composed of kite-shaped structural units—previously thought to be well understood. By performing a rigorous GB structure search that optimizes the number of atoms in the boundary core, we predict different quasiaperiodic “split kite” phases analogous to those previously found in GBs in face-centered-cubic metals. Furthermore, our results suggest that complex aperiodic phases of GBs appear to be a general phenomenon, as validated through density functional theory calculations. Moreover, the atoms in the split kite phase demonstrate distinct collective diffusion dynamics. Phase-contrast image simulations of split kites show better agreement with experimental observations, offering an alternative explanation for previous microscopy results and motivating future atomically resolved imaging of the GB structure.

Body-centered cubic↗
DOE OSTI2024

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reduction-Controlled Atomic Migration for Single Atom Alloy Library

Here, picturing the atomic migration pathways of catalysts in a reactive atmosphere is of central significance for uncovering the underlying catalytic mechanisms and directing the design of high-performance catalysts. Here, we describe a reduction-controlled atomic migration pathway that converts nanoparticles to single atom alloys (SAAs), which has remained synthetically challenging in prior attempts due to the elusive mechanism. We achieved this by thermally treating the noble-metal nanoparticles M (M = Ru, Rh, Pd, Ag, Ir, Pt, and Au) on metal oxide (CuO) supports with H2/Ar. Atomic-level characterization revealed such conversion as the synergistic consequence of noble metal-promoted H2 dissociation and concomitant CuO reduction. The observed atomic migration pathway offers an understanding of the dynamic mechanisms study of nanomaterials formation and catalyst design.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Effect of dispersion forces on squeezing with Rydberg atoms

We report exact results concerning the effect of dipole-dipole interaction (dispersion forces) on dynamic and steady-state characteristics of squeezing in the emitted fluorescent field from two identical coherently driven two-level atoms. The atomic system is subjected to three different damping baths in particular the normal vacuum, a broad band thermal field and a broad band squeezed vacuum. The atomic model is the Dicke model, hence possible experiments are most likely to agree with theory when performed on systems of Rydberg atoms making microwave transitions. The presence of dipole-dipole interaction can enhance squeezing for realizable values of the various parameters involved.

Ng, S. K.↗
DOE OSTI2022

Mechanisms of Retention and Transport of Fission Products in Virgin and Irradiated Nuclear Graphite

To fulfill the design and licensing requirements for advanced reactors, the U.S. Department of Energy (DOE) is specifically interested in targeted research to quantify the mechanisms of fission product transport and retention in graphitic grades that will likely be used in next-generation very/high temperature reactors, fluoride salt-cooled high temperature reactors and molten salt reactors. We have examined an integrated experimental-computational approach to determine the diffusivities of selected fission products – ruthenium (Ru) and silver (Ag) in five nuclear graphite grades – POCO AXF-5Q, POCO ZXF-5Q, PCEA, IG110, and NBG-18. Experimental investigation for diffusion of other fission products (iodine, cesium, and krypton) in graphite could not be carried out due to extreme evaporation due to high vapor pressure, despite several attempts incorporating various experimental schemes. Experimental investigation with irradiated graphite could not be carried out because necessary approval to transport irradiated graphite could not be obtained due to procedural reasons at ORNL and national policy change in UK. An effort was made to “produce” by ion sputtering, but it was discontinued because the outcome did not adequately represent irradiated graphite. Thin film and “cup-and-cap” methods were employed to deposit fission products, diffusion annealed were carried out in an argon atmosphere, and concentration profiles were determined by dynamic secondary ion mass spectroscopy depth profiling technique. Regardless of graphite type, anomalous diffusion behavior of Ru (i.e., lower than expected based on Arrhenius temperature dependence) at an intermediate temperature, around 600 to 700 °C, was observed. Diffusion behavior of silver on the other hand agreed well to the Arrhenius temperature dependence although the excessive evaporation of silver hampered the accurate determination of diffusivities at temperature higher than 900 °C. In general, silver was observed to diffuse faster than Ru in the temperature range compared from 500 to 900 °C, had consistently had higher pre-exponential factor corresponding to higher jump frequency. To understand the abnormal diffusional behavior of Ru in graphite in intermediate temperature range, a detailed microstructural analysis was carried out. Clustering of Ru into particles on the surface of graphite was observed around 700 °C corresponding to de-wetting and spheroidization. This change would influence the boundary condition of diffusion, although features of clustering was much smaller (<1 mm) than the SIMS raster size (120 x 120 mm). Raman spectroscopy also demonstrated that Raman peak intensity for Ru decreased significantly for samples annealed higher than 700 °C and presence of Ru in graphite altered the crystallinity of graphite. To further elucidate the abnormal diffusion behavior observed in intermediate temperature investigated, the vibrational properties of graphene, bilayer graphene and Kr atoms in bilayer graphene were investigated with molecular dynamics. Krypton, which is a fission product, was introduced between a bilayer graphene system to evaluate (i) whether the phonon dispersion curves change with high temperature in bilayer graphene (ii) whether the addition of Kr atoms between the bilayer introduces any significant change in the vibrational properties, and (iii) whether any change in the vibrational properties can be correlated to the diffusive behavior of Kr atoms. The results show that the Kr atoms indeed exhibit an anomalous diffusive behavior at 550°C. Phonon analysis shows that a few phonon frequencies are enhanced as much as 20 to 30 % with the addition of just 0.2 at.% Kr. So, this may explain the unexpected self-diffusivity of fission products.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

28. Dynamics of Molecular Collisions Conference

The DMC meeting, since 1965, has a great history of providing a unique platform and focal point for the gathering of experimentalists and theoreticians in the field of chemical dynamics. At this biannual conference which spans five days, the latest developments in bimolecular reactions, unimolecular reactions, reactions in clusters and solutions, reactions at interfaces, reactions in combustion, atmospheres, and interstellar media, as well as cold chemistry, new spectroscopic techniques and non-adiabatic reactions will be discussed by experts in these fields. The 2021 meeting was postponed due to Covid but is back in 2023 including the traditional awarding of two Herschbach Medals, honoring important contributions to the fields of theoretical and experimental dynamics respectively. The 2023 meeting will be organized by Richard Dawes (Chair) and Amy Mullin (Vice Chair). The DMC format encourages participation of young scientists through invited and contributed talks as well as posters. This proposal seeks support for students, post-docs and invited scientists to attend the meeting.

2021 conference↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2024

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

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

Laser-Induced Thermal Decomposition of Uranium Coordination Compounds with Non-oxidic Ligands to Produce Nitride and Carbide Materials

The production of ceramics from uranium coordination compounds can be achieved through thermal processing if an excess amount of the desired atoms (i.e., C or N), or reactive gaseous products (e.g., methane or nitrogen oxide) is made available to the reactive uranium metal core via decomposition/fragmentation of the surrounding ligand groups. Here, computational thermodynamic approaches were utilized to identify the temperatures necessary to produce uranium metal from some starting compounds—UI 4 (TMEDA) 2 , UCl 4 (TMEDA) 2 , UCl 3 (pyridine) x , and UI 3 (pyridine) 4 . Experimentally, precursors were irradiated by a laser under various gaseous environments (argon, nitrogen, and methane) creating extreme reaction conditions (i.e., fast heating, high temperature profile >2000 °C, and rapid cooling). Despite the fast dynamics associated with laser irradiation, the central uranium atom reacted with the thermal decomposition products of the ligands yielding uranium ceramics. Residual gas analysis identified vaporized products from the laser irradiation, and the final ceramic products were characterized by powder X-ray diffraction. The composition of the uranium precursor as well as the gaseous environment had a direct impact on the production of the final phases.

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

Tuning the Coordination Environment of Rh Single Atoms on Highly Dispersed Reducible Oxides for Enhanced Reverse Water-Gas Shift Performance

Controlling the dynamic mobility of catalyst surface active sites and their interactions with the surrounding environment is critical in generating active surfaces that directly influence the catalytic activity and selectivity. Here, we report a strategy for tailoring the dispersion and electronic environment of single atom Rh catalysts by decorating the alumina support with highly dispersed (HD) cerium and molybdenum oxides. The resulting catalysts exhibit markedly different behavior in the Reverse Water Gas Shift (RWGS) reaction. In particular, Rh/MoOx(HD)/Al2O3 maintains atomically dispersed Rh even at elevated temperatures (up to 400 °C), achieving CO selectivity of up to 100% and resists sintering via the formation of a newly developed structure featuring Rh single atoms embedded in MoOx clusters. In situ spectroscopy and microscopy analyses confirm the stabilization of Rh and the dynamic evolution of Rh–Mo coordination under reaction conditions. Our findings highlight the power of support modification in steering active site structure and activity, offering a pathway toward enhanced and tunable single atom catalysts for CO2 valorization.

CO selectivity↗
DOE OSTI2021

On the Development of Order and Interfaces during the Growth of Ultrathin La 2 CuO 4 Films by Molecular Beam Epitaxy

While the atomic structure of interfaces in complex oxide heterostructures created by epitaxial growth has been investigated extensively, few studies have been conducted on how interfaces form and restructure at the initial stage of film growth. The dynamic aspects of the growth behavior can strongly influence the final interfacial atomic structure, which may lead to the emergence of interface-specific properties, such as the rise of interfacial superconductivity between certain Ruddlesden-Popper oxide materials. In this report, the structural development of La2CuO 4 thin films grown by molecular beam epitaxy on a LaSrAlO 4 substrate is investigated by X-ray diffraction measurements with rapid scans over a volume of reciprocal space. This method provides far more detailed information on interface formation than traditional fixed-point measurements. The results show that the atomic structure of the interface becomes fully established after just a single unit cell of growth. Interestingly, restructuring continues to occur within the topmost half to one unit cell of the film during the deposition process. However, diffraction intensity oscillations from both reflection high-energy electron and X-ray measurements stabilize only after the growth of two unit cells, indicating that the growth front morphology continues to evolve until the start of the third unit cell. This multimodal investigation provides insights into the atomic processes taking place during layered oxide interface formation, including the dynamical rearrangement of LaO and CuO 2 layers. Such information is not only relevant to the engineering and optimization of functional layer structures but can also be critical for ultrathin films.

36 MATERIALS SCIENCE↗
DOE OSTI2025

X-ray crystallographic and hydrogen deuterium exchange studies confirm alternate kinetic models for homolog insulin monomers

Despite the crucial role of various insulin analogs in achieving satisfactory glycemic control, a comprehensive understanding of their in-solution dynamic mechanisms still holds the potential to further optimize rapid insulin analogs, thus significantly improving the well-being of individuals with Type 1 Diabetes. Here, we employed hydrogen-deuterium exchange mass spectrometry to decipher the molecular dynamics of newly modified and functional insulin analog. A comparative analysis of H/D dynamics demonstrated that the modified insulin exchanges deuterium atoms faster and more extensively than the intact insulin aspart. Additionally, we present new insights derived from our 2.5 Å resolution X-ray crystal structure of modified hexamer insulin analog at ambient temperature. Furthermore, we obtained a distinctive side-chain conformation of the Asn3 residue on the B chain (AsnB3) by operating a comparative analysis with a previously available cryogenic rapid-acting insulin structure (PDB_ID: 4GBN). The experimental conclusions have demonstrated compatibility with modified insulin’s distinct cellular activity, comparably to aspart. Additionally, the hybrid structural approach combined with computational analysis employed in this study provides novel insight into the structural dynamics of newly modified and functional insulin vs insulin aspart monomeric entities. It allows further molecular understanding of intermolecular interrelations driving dissociation kinetics and, therefore, a fast action mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2022

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

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

Nanoscale dynamics of Dynamin 1 helices reveals squeeze-twist deformation mode critical for membrane fission

Dynamin 1 (Dyn1) GTPase, a principal driver of membrane fission during synaptic endocytosis, self-assembles into short mechanoactive helices cleaving the necks of endocytic vesicles. While structural information about Dyn1 helix is abundant, little is known about the nanoscale dynamics of the helical scaffolding at the moment of fission, complicating mechanistic understanding of Dyn1 action. To address the role of the helix dynamics in fission, we used High-Speed Atomic Force Microscopy (HS-AFM) and fluorescence microscopy to track and compare the spatiotemporal characteristics of the helices formed by wild-type Dyn1 and its K44A mutant impaired in GTP hydrolysis on minimal lipid membrane templates. In the absence of nucleotide, membrane-bound WT Dyn1 and K44A Dyn1 self-assembled into tubular protein scaffolding of similar diameter encaging the lipid bilayer. In both cases, the GTP addition caused scaffold constriction coupled with formation of 20 to 30 nm nanogaps in the protein coverage. While both proteins reached scaffold diameters characteristic for membrane superconstriction causing fission, the fission was detected only with WT Dyn1. We associated the fission activity with the dynamic evolution of the nanogaps: K44A Dyn1 gaps were static, while WT Dyn1 gaps actively evolved via repetitive nonaxisymmetric constriction-bending deformations caused by localized GTP hydrolysis. Modeling of the deformations implicated filament twist as an additional deformation mode which combines with superconstriction to facilitate membrane fission. Our results thus show that the dynamics of the Dyn1 helical scaffold goes beyond radial constriction and involves nonaxisymmetric deformations, where filament twist emerges as a critical driver of membrane fission.

59 BASIC BIOLOGICAL SCIENCES↗