Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pattern”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,243 records · Page 69

Engineering Assembly Kinetics and Line Roughness in Solvent Vapor-Annealed Block Copolymer/Homopolymer Blends

Block copolymer (BCP) directed self-assembly (DSA) is a promising route to enhance lithography resolution by multiplying nanopattern density and reducing feature roughness. Eliminating kinetically trapped self-assembly defects requires fast self-assembly. However, acceleration strategies like solvent vapor annealing or homopolymer blending broaden domain interfaces, implying a trade-off in increased feature roughness. In this work, we experimentally investigate this apparent dilemma between self-assembly kinetics and line roughness for solvent vapor-annealed thin films of a lamellar poly(styrene-block-2-vinylpyridine) (PS-b-P2VP) BCP blended with PS and P2VP homopolymers. Binary blends with PS or P2VP homopolymers and ternary blends incorporating both in equal weight fractions were solvent vapor annealed using acetone, a near-neutral solvent for PS and P2VP, followed by P2VP-selective vapor-phase infiltration with alumina (AlOx) and polymer etching. Binary blends with P2VP exhibit a modest kinetic enhancement but also higher line-edge and -width roughness due to the increased frequency of P2VP protrusions and bridge defects in the alumina line patterns. In contrast, binary blends with PS self-assemble noticeably faster, while domain asymmetry from the added PS homopolymer reduces roughness by curbing the number of alumina protrusions and bridge defects. Ternary blends maintain DSA line patterns across a wider composition window and, at higher homopolymer loadings, reduce roughness at length scales near the lamellar period, consistent with a reduced impact of intradomain compositional fluctuations. These findings provide important insights for codesigning blend compositions and process flows to achieve high-resolution, defect free patterns with minimal roughness through BCP DSA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wet Brush Homopolymers as “Smart Solvents” for Rapid, Large Period Block Copolymer Thin Film Self-Assembly

The spontaneous formation of periodic nanopatterns through the self-assembly of block copolymer (BCP) thin films offers an enticing path to engineer functional material properties over large areas, but the prohibitively slow ordering of ultrahigh molecular weight BCPs that generate large pattern periods (>100 nm) severely restricts the design space in which BCPs may be used. Here, we demonstrate that blending lamellae-forming polystyrene-block-poly(methyl methacrylate) diblock copolymers with very low molecular weight polystyrene and poly(methyl methacrylate) homopolymers can reduce by an order-of-magnitude the times required to form ordered large period patterns through solvent vapor annealing. Using in situ grazing-incidence small-angle X-ray scattering and ex situ electron microscopy, we show that the added homopolymer plasticizers act as a “smart solvent”: They increase the mobility of the long BCP chains by amplifying the dilution effect during solvent vapor annealing but segregate to like domains during subsequent thermal annealing to produce high-quality patterns with periods that may exceed 200 nm, which are readily transferred to other materials. The dramatic impact of homopolymer plasticizers may help extend the use of BCPs into new areas, particularly for visible light photonic applications that demand larger period nanopatterns.

36 MATERIALS SCIENCE↗

Atom Identification in Bilayer Moiré Materials with Gomb-Net

Moiré patterns in van der Waals bilayer materials complicate the analysis of atomic-resolution images, hindering the atomic-scale insight typically attainable with scanning transmission electron microscopy. Here, we report a method to detect the positions and identities of atoms in each of the individual layers that compose twisted bilayer heterostructures. We developed a deep learning model, Gomb-Net, which identifies the coordinates and atomic species in each layer, deconvoluting the moiré pattern. This enables layer-specific mapping of atomic positions and dopant distributions, unlike other commonly used segmentation models which struggle with moiré-induced complexity. Using this approach, we explored the Se atom substitutional site distribution in a twisted fractional Janus WS 2 -WS 2(1–x) Se 2x heterostructure and found that layer-specific implantation sites are unaffected by the moiré pattern’s local energetic or electronic modulation. In conclusion, this advancement enables atom identification within material regimes where it was not possible before, opening new insights into previously inaccessible material physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sub-10 nm Probing of Ferroelectricity in Heterogeneous Materials by Machine Learning Enabled Contact Kelvin Probe Force Microscopy

Reducing the dimensions of ferroelectric materials down to the nanoscale has strong implications on the ferroelectric polarization pattern and on the ability to switch the polarization. As the size of ferroelectric domains shrinks to the nanometer scale, the heterogeneity of the polarization pattern becomes increasingly pronounced, enabling a large variety of possible polar textures in nanocrystalline and nanocomposite materials. Critical to the understanding of fundamental physics of such materials and hence their applications in electronic nanodevices is the ability to investigate their ferroelectric polarization at the nanoscale in a nondestructive way. We show that contact Kelvin probe force microscopy (cKPFM) combined with a k-means response clustering algorithm enables to measure the ferroelectric response at a mapping resolution of 8 nm. In a BaTiO 3 thin film on silicon composed of tetragonal and hexagonal nanocrystals, we determine a nanoscale lateral distribution of discrete ferroelectric response clusters, fully consistent with the nanostructure determined by transmission electron microscopy. Moreover, we apply this data clustering method to the cKPFM responses measured at different temperatures, which allows us to follow the corresponding change in the polarization pattern as the Curie temperature is approached and across the phase transition. This work opens up perspectives for mapping complex ferroelectric polarization textures such as curled/swirled polar textures that can be stabilized in epitaxial heterostructures and more generally for mapping the polar domain distribution of any spatially highly heterogeneous ferroelectric materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Genesis and Propagation of Fractal Structures During Photoelectrochemical Etching of n-Silicon

The genesis, propagation, and dimensions of fractal-etch patterns that form anodically on front- or back-illuminated n-Si(100) photoelectrodes in contact with 11.9 M NH 4 F (aqueous) have been investigated during either a linear potential sweep or a constant potential hold ( E = +6.0 V versus Ag/AgCl). Optical images collected in situ during electrochemical experiments reflected the location and underlying mechanism of initiation and propagation of the structures on the surface. X-ray photoelectron spectroscopic (XPS) data collected for samples emersed from the electrolyte at varied times provided detailed information about the chemistry of the surface during fractal etching. The fractal structure was strongly influenced by the orientation of the crystalline Si sample. Furthermore, the etch patterns were initially generated at points along the circumference of bubbles that formed upon immersion of n-Si(100) samples in the electrolyte, most likely due to the electrochemical and electronic isolation of areas beneath bubbles. XPS data showed the presence of a tensile-stressed silicon surface throughout the etching process as well as the presence of SiO x F y on the surface. The two-dimensional fractal dimension, D f,2D , of the patterns increased with etching time to a maximum observed value of D f,2D = 1.82. Promotion of fractal etching near etch masks that electrochemically and electronically isolated areas of the photoelectrode surface enabled the selective placement of highly branched structures at desired locations on an electrode surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-Dimensional Molecular Charge Density Waves in Single-Layer-Thick Islands of a Dirac Fermion System

Charge density waves have been intensely studied in inorganic materials such as transition metal dichalcogenides however its counterpart in organic materials has yet to be explored in detail. Here we report the finding of a robust two dimensional charge density waves in molecular layers formed by α-(BEDT-TTF) 2 -I 3 on a Ag(111) surface. Low temperature scanning tunnelling microscopy images of a multi-layer thick α-(BEDT-TTF) 2 -I 3 on Ag(111) substrate reveal coexistence of 5a0 x 5a0 and √31 a0 x √31 a0 R9° charge density wave patterns commensurate with the underlying molecular lattice at 80 K. Both charge density wave patterns remain in nano-size molecular islands with just a single constituent molecular-layer thickness at 80 K and 5 K. Local tunneling spectroscopy measurements reveal the variation of the gap from 244 meV to 288 meV between the maximum and minimum charge density wave locations. Density functional theory calculations further confirm a vertical positioning of BEDT-TTF molecules in the molecular layer. While the observed charge density wave patterns are stable for the defect sites, they can be reversely switched for one molecular lattice site by means of inelastic tunnelling electron energy transfer with the electron energies exceeding 400 meV using a scanning tunneling microscope manipulation scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extrinsic Nonlinear Kerr Rotation in Topological Materials under a Magnetic Field

Topological properties in quantum materials are often governed by symmetry and tuned by crystal structure and external fields, and hence, symmetry-sensitive nonlinear optical measurements in a magnetic field are a valuable probe. Here, we report nonlinear magneto-optical second harmonic generation (SHG) studies of nonmagnetic topological materials including bilayer WTe 2 , monolayer WSe 2 , and bulk TaAs. The polarization-resolved patterns of optical SHG under a magnetic field show nonlinear Kerr rotation in these time-reversal symmetric materials. For materials with 3-fold rotational symmetric lattice structure, the SHG polarization pattern rotates just slightly in a magnetic field, whereas in those with mirror or 2-fold rotational symmetry, the SHG polarization pattern rotates greatly and distorts. Further, these different magneto-SHG characters can be understood by considering the superposition of the magnetic field-induced time-noninvariant nonlinear optical tensor and the crystal-structure-based time-invariant counterpart. The situation is further clarified by scrutinizing the Faraday rotation, whose subtle interplay with crystal symmetry accounts for the diverse behavior of the extrinsic nonlinear Kerr rotation in different materials. Our work illustrates the application of magneto-SHG techniques to directly probe nontrivial topological properties, and underlines the importance of minimizing extrinsic nonlinear Kerr rotation in polarization-resolved magneto-optical studies.

36 MATERIALS SCIENCE↗

Solid-State Lighting Using Side-by-Side White Phosphorescent Organic Light-Emitting Diodes

White organic light-emitting diodes (WOLEDs) have become increasingly popular for use in solid-state illumination, where diffuse, large area light sources that achieve a high color rendering index and luminous power efficiency are desirable. Color-tunable emission, where the light source combines emission from multiple, separately addressed color elements, is conveniently provided by WOLEDs for the purpose of adapting the lighting source to a particular illumination requirement. In this work, we demonstrate a method for side-by-side positioning of monochromatic blue and yellow phosphorescent OLED stripes that are combined to create tunable white light, using a high resolution mechanical peel-off patterning method. We achieve a peak luminous power efficiency of 17.1 ± 0.3 lm W –1 and an external quantum efficiency of 11.8 ± 0.2% and demonstrate color tunability of the 1960 Commission Internationale d’Eclairage chromaticity coordinates from (u,v) = (0.33,0.36) to (0.12,0.32). This corresponds to a tuning rnge of the color rendering index from 74 ± 1 to 86 ± 1 and the correlated color temperature from 2000 to 8000 ± 500 K. Here, due to the nondestructive nature of the peel-off technique, patterned devices exhibit a lifetime comparable to conventional, shadow mask-patterned devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

J -Resolved Molecular Fingerprinting by Parahydrogen Hyperpolarized Low-Field NMR

A J-resolved spectroscopy that depends on homonuclear scalar coupling in the strong-coupling regime and heteronuclear coupling in the weak regime expands complex peak patterns to a second axis. Hyperpolarization by Signal Amplification by Reversible Exchange (SABRE) enables the spectroscopy at a low magnetic field of 0.82 mT. Overlapping peaks of molecules such as 3-fluoropyridine and 3,5-difluoropyridine are resolved. Density matrix simulations of the 1 H and 19 F spins indicate a strong dependence on the signs and values of the J-coupling constants, including the homonuclear couplings that are not directly observable. The best matching peak positions and intensities predict coupling constants, including couplings between chemically equivalent nuclear spins, ranging in magnitude from 0.4 to 9.0 Hz for the two molecules. Simulations of other spin systems show unique patterns for molecules containing 1 H and 19 F or 13 C. The dependence of the J-resolved peak patterns on all coupling constants in a spin system presents a new modality for portable and inexpensive identification of molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relative Contributions of Atmospheric, Oceanic, and Coupled Processes to North Pacific and North Atlantic Variability

Patterns of sea surface temperature (SST) variability over the northern oceans arise from a combination of atmospheric, oceanic, and coupled processes. In this work, we use a novel methodology and a suite of observations to quantify the processes contributing to the dominant patterns of interannual SST variability over these regions. We decompose the upper ocean heat content tendency associated with such dominant patterns into contributions from different heat fluxes: (a) atmospherically driven, (b) surface feedbacks, and (c) oceanic. We find that in the subtropics, cloud radiative flux, turbulent heat flux, and residual oceanic processes each contributes substantially to North Pacific SST variability, whereas turbulent heat flux primarily induces North Atlantic SST variability. Cloud radiative fluxes therefore provide a major source of interannual SST variability in the North Pacific but not in the North Atlantic. In midlatitudes, SST fluctuations over the northern oceans are driven by the combination of turbulent and oceanic heat fluxes.

54 ENVIRONMENTAL SCIENCES↗

Varying Contributions of Drivers to the Relationship Between Canopy Photosynthesis and Far‐Red Sun‐Induced Fluorescence for Two Maize Sites at Different Temporal Scales

Abstract Sun‐induced fluorescence (SIF) has been found to be strongly correlated with gross primary production (GPP) in a quasi‐linear pattern at the scales beyond leaves. However, the causes of the GPP:SIF relationship deviating from a linear pattern remain unclear. In the current study conducted at two maize sites in Nebraska in 2017 summer growing season, we investigated the relationship between GPP and SIF at 760 nm (F 760 ) at two temporal scales and quantified the contributions of incoming photosynthetically active radiation (PAR in ), fraction of absorbed PAR (fPAR), light use efficiency (LUE), and F 760 yield (F 760,y , defined as F 760 /(PAR in ×fPAR)) to GPP and F 760 variabilities to further understand the linearity and deviations in the GPP:F 760 relationship. We found the following: (1) For individual growth stages when canopy structure and chlorophyll content were stable, GPP and F 760 were strongly controlled by PAR in , while LUE and F 760,y had much lower contributions to the GPP:F 760 relationship; during this period, LUE and F 760,y had either a slightly negative or no clear relationship, which explained some deviations in the GPP:SIF relationship. (2) At the seasonal scale, the contribution of LUE to GPP variability as well as the contribution of F 760,y to F 760 variability increased and was comparable to the contribution of PAR in ; the LUE:F 760,y relationship showed a strong linear relationship, which strengthened the linear GPP:F 760 relationship. Both maize sites showed similar patterns. A framework was applied to estimate LUE at individual stages and as a result, significantly improved the GPP estimation, thus enhancing the SIF potential for inferring photosynthesis.

Miao, Guofang↗

Stable Methane Isotopologues From Northern Lakes Suggest That Ebullition Is Dominated by Sub‐Lake Scale Processes

Abstract Stable isotopes have emerged as popular study targets when investigating emission of methane (CH 4 ) from lakes. Yet little is known on how isotopic patterns conform to variations in emission magnitudes—a highly relevant question. Here, we present a large multiyear data set on stable isotopes of CH 4 ebullition (bubbling) from three small adjacent subarctic lakes. The δ 13 C‐CH 4 and δD‐CH 4 range from −78.4‰ to −53.1‰ and from −369.8‰ to −218.8‰, respectively, and vary greatly among the lakes. The signatures suggest dominant hydrogenotrophic methanogenesis, particularly in the deep zones, but there are also signals of seemingly acetoclastic production in some high fluxing shallow areas, possibly fueled by in situ vegetation, but in‐sediment anaerobic CH 4 oxidation cannot be ruled out as an alternative cause. The observed patterns, however, are not consistent across the lakes. Neither do they correspond to the spatiotemporal variations in the measured bubble CH 4 fluxes. Patterns of acetoclastic and hydrogenotrophic production plus oxidation demonstrate that gains and losses of sediment CH 4 are dominated by sub‐lake scale processes. The δD‐CH 4 in the bubbles was significantly different depending on measurement month, likely due to evaporation effects. On a larger scale, our isotopic data, combined with those from other lakes, show a significant difference in bubble δD‐CH 4 between postglacial and thermokarst lakes, an important result for emission inventories. Although this characteristic theoretically assists in source partitioning studies, most hypothetical future shifts in δD‐CH 4 due to high‐latitude lake area or production pathway are too small to lead to atmospheric changes detectable with current technology.

54 ENVIRONMENTAL SCIENCES↗

Midwest US Croplands Determine Model Divergence in North American Carbon Fluxes

Large uncertainties in North American terrestrial carbon fluxes hinder regional climate projections. Terrestrial biosphere models (TBMs), the essential tools for understanding continental-scale carbon cycle, diverge on whether temperate forests or croplands dominate carbon uptake in North America. Evidence from novel photosynthetic proxies, such as those based on chlorophyll fluorescence, has cast doubt on the “weak cropland, strong forest” carbon uptake patterns simulated by most TBMs. However, no systematic evaluation of TBMs has yet been attempted to pin down space-time patterns that are most consistent with regional CO 2 observational constraints. Here, we leverage atmospheric CO 2 observations and satellite-observed photosynthetic proxies to understand emergent space-time patterns in North American carbon fluxes from a large suite of TBMs and data-driven models. To do so, we evaluate how well the atmospheric signals resulting from carbon flux estimates reproduce the space-time variability in atmospheric CO 2 , as is observed by a network of continuous-monitoring towers over North America. Models with gross or net carbon fluxes that are consistent with the observed CO 2 variability share a salient feature of growing-season carbon uptake in Midwest US croplands. Conversely, the remaining models place most growing-season uptake in boreal or temperate forests. Differences in model explanatory power depend mainly on the simulated annual cycles of cropland uptake—especially, the timing of peak uptake—rather than the distribution of annual mean fluxes across biomes. Our results suggest that improved model representation of cropland phenology is crucial to robust, policy-relevant estimation of North American carbon exchange.

54 ENVIRONMENTAL SCIENCES↗

Identifying Key Drivers of Wildfires in the Contiguous US Using Machine Learning and Game Theory Interpretation

Abstract Understanding the complex interrelationships between wildfire and its environmental and anthropogenic controls is crucial for wildfire modeling and management. Although machine learning (ML) models have yielded significant improvements in wildfire predictions, their limited interpretability has been an obstacle for their use in advancing understanding of wildfires. This study builds an ML model incorporating predictors of local meteorology, land‐surface characteristics, and socioeconomic variables to predict monthly burned area at grid cells of 0.25° × 0.25° resolution over the contiguous United States. Besides these predictors, we construct and include predictors representing the large‐scale circulation patterns conducive to wildfires, which largely improves the temporal correlations in several regions by 14%–44%. The Shapley additive explanation is introduced to quantify the contributions of the predictors to burned area. Results show a key role of longitude and latitude in delineating fire regimes with different temporal patterns of burned area. The model captures the physical relationship between burned area and vapor pressure deficit, relative humidity (RH), and energy release component (ERC), in agreement with the prior findings. Aggregating the contribution of predictor variables of all the grids by region, analyses show that ERC is the major contributor accounting for 14%–27% to large burned areas in the western US. In contrast, there is no leading factor contributing to large burned areas in the eastern US, although large‐scale circulation patterns featuring less active upper‐level ridge‐trough and low RH two months earlier in winter contribute relatively more to large burned areas in spring in the southeastern US.

54 ENVIRONMENTAL SCIENCES↗

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗

Seismic Imaging of a Shale Landscape Under Compression Shows Limited Influence of Topography–Induced Fracturing

We used seismic refraction to image the P-wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill-valley-hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. Here, this contradiction could be explained by the greater importance of infiltration-driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.

58 GEOSCIENCES↗

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗

Spring Regional Sea Surface Temperatures as a Precursor of European Summer Heatwaves

Different spring and early summer North Atlantic sea surface temperature anomalies (SSTAs) have been shown to precede recent European summer heatwaves (EuSHWs). So far, the limited number of observed events associated with several physical mechanisms has prevented a robust identification of SSTAs as precursors. Here, we extend beyond previous studies by combining 100 historical simulations (1850–2005) of the MPI Grand-Ensemble with an explainable neural-network method. We find that the spring tripolar North Atlantic pattern with positive SSTAs in the Subtropical Gyre is a precursor of EuSHWs. In addition, positive SSTAs west of the Iberian Peninsula, and in the North Sea, the Baltic Sea and the Mediterranean Sea relate to distinct early summer soil moisture anomaly patterns and are precursors of western and southeastern EuSHWs, respectively. While the phase of the tripolar North Atlantic pattern indicates whether a EuSHW might emerge, regional SSTAs indicate the spatial characteristics of EuSHWs.

54 ENVIRONMENTAL SCIENCES↗