Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional calculation”

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

Uncovering an Interfacial Band Resulting from Orbital Hybridization in Nickelate Heterostructures

The interaction of atomic orbitals at the interface of perovskite oxide heterostructures has been investigated for its profound impact on the band structures and electronic properties, giving rise to unique electronic states and a variety of tunable functionalities. In this study, we conducted an extensive investigation of the optical and electronic properties of epitaxial NdNiO 3 thin films grown on a series of single crystal substrates. Unlike films synthesized on other substrates, NdNiO 3 on SrTiO 3 (NNO/STO) gives rise to a unique band structure which features an additional unoccupied band situated above the Fermi level. Our comprehensive investigation, which incorporated a wide array of experimental techniques and density functional theory calculations, revealed that the emergence of the interfacial band structure is primarily driven by the orbital hybridization between Ti 3d orbitals of the STO substrate and O 2p orbitals of the NNO thin film. Furthermore, exciton peaks have been detected in the optical spectra of the NNO/STO film, attributable to the pronounced electron-electron (e-e) and electron-hole (e-h) interactions propagating from the STO substrate into the NNO film. Further, these findings underscore the substantial influence of interfacial orbital hybridization on the electronic structure of oxide thin-films, thereby offering key insights into tuning their interfacial properties.

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

Lithium's low-temperature phase transitions: Insights into quantum lattice dynamics and superconductivity

The large lattice dynamics of lithium, driven by its low atomic mass, results in energetically similar structures and significant isotope effects under pressure, posing challenges to current theoretical models. Above 20 GPa and at low temperatures, lithium's electronic properties deviate from simple metallic behavior, with superconductivity emerging in a complex, pressure-dependent manner, alongside an unusual isotope effect. The structural phases of 7 Li reported under these conditions are inconsistent across studies, and the structures of 6 Li remain unexamined. These gaps limit our understanding of the effects of pressure on lithium's electronic properties and the role of quantum lattice effects on its structural behavior under pressure. Here, we integrate experimental and theoretical approaches to investigate the low-temperature structural phase boundaries in lithium isotopes. We map the structural phase diagram of 7 Li from 5 to 55 GPa and 15–75 K, identifying the sequence fcc → ℎ⁡R⁢1 → cI16. A pronounced isotope effect is observed, with 6 Li shifting the fcc → ℎ⁡R⁢1 phase boundary to lower pressures at 15 K. Density functional theory calculations further clarify how these structural changes affect superconducting properties, particularly emphasizing the role of the fcc → ℎ⁡R⁢1 transition in lithium's superconductivity. Furthermore, our findings offer insights into the unique behavior of lithium isotopes under pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Intercalation induced quasi-freestanding layer in TiSe 2

Angle-resolved photoemission spectroscopy is employed to study the electronic structure of bulk TiSe 2 before and after doping with potassium impurities. A splitting in the conduction band into two branches is observed after room-temperature deposition. The splitting energy increases to approximately 130 meV when the sample is cooled to 40 K. One branch exhibits a nondispersive two-dimensional feature, while the other shows the characteristics of three-dimensional bulk band dispersion. Core-level spectroscopy suggests that the K impurities predominantly occupy the intercalated sites within the van der Waals gap. Furthermore, the results indicate the formation of a quasi-freestanding TiSe 2 layer. Additionally, doping completely suppresses the periodic lattice distortion in the surface region. These findings are further supported by density functional theory calculations, which compare the band structure of monolayer and bulk TiSe 2 with experimental data. Thus, the dimensional and intrinsic electronic properties of 1⁢𝑇−TiSe 2 can be controlled through the intercalation procedure used in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2026

First-principles insights into chromium-induced oxide phases in NiO

The high-temperature oxidation of Ni–Cr alloys leads to complex oxide scales comprising Ni(Cr)O solid solutions, NiCr 2 O 4 spinel, and Cr 2 O 3 corundum within the NiO matrix. Understanding the atomic-scale mechanisms of Cr segregation and precipitate formation is crucial for enhancing oxidation resistance. Here, we employ density functional theory calculations to investigate Cr behavior on NiO(100), (110), and (111) surfaces and in the bulk. Our results reveal that isolated Cr atoms preferentially segregate to the surfaces, stabilizing Ni(Cr)O solid solutions via strong Cr–O bonding, whereas Cr pairs and clusters favor subsurface migration and bulk aggregation, promoting nucleation of NiCr 2 O 4 and Cr 2 O 3 phases. These findings elucidate a size-dependent segregation mechanism linking Cr coordination environments to oxide phase evolution. This atomic-scale insight informs strategies to tailor oxide microstructures and enhance the high-temperature oxidation resistance of Ni–Cr alloys.

Chromium Segregation↗
DOE OSTI2024

Evolution of structure and spectroscopic properties of a new 1,3-diacetylpyrene polymorph with temperature and pressure

A new polymorph of 1,3-diacetylpyrene has been obtained from its melt and thoroughly characterized using single-crystal X-ray diffraction, steady-state UV–Vis spectroscopy and periodic density functional theory calculations. Experimental studies covered the temperature range from 90 to 390 K and the pressure range from atmospheric to 4.08 GPa. Optimal sample placement in a diamond anvil cell according to our previously presented methodology ensured over 80% data coverage up to 0.8 Å for a monoclinic sample. Unrestrained Hirshfeld atom refinement of the high-pressure crystal structures was successful and anharmonic behavior of carbonyl oxygen atoms was observed. Unlike the previously characterized polymorph, the structure of 2°AP-β is based on infinite π-stacks of antiparallel 2°AP molecules. 2°AP-β displays piezochromism and piezofluorochromism which are directly related to the variation in interplanar distances within the π-stacking. The importance of weak intermolecular interactions is reflected in the substantial negative thermal expansion coefficient of -55.8 (57) MK -1 in the direction of C—H∙∙∙O interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

4-Electron Oxygen Reduction Reaction (ORR) with Iron Phthalocyanine (FePc) Functionalized Nanowire Templated-3D Fuzzy Graphene (NT-3DFG)

Iron phthalocyanine (FePc) is a promising alternative to platinum-based catalysts for sustainable energy devices; however, the plane-symmetry of Fe-N 4 sites, random aggregation, and poor conductivity of FePc present major barriers for their application as oxygen reduction reaction (ORR) electrocatalysts. Here, we report the synergistic effects of FePc electrocatalysts supported by a nanowire-templated three-dimensional fuzzy graphene (FePc@NT-3DFG) substrate. The in situ functionalized oxygen groups (iFOGs) at the edge of NT-3DFG localize Fe active sites in FePc under alkaline ORR conditions. With a uniform FePc distribution through many single layers of graphene, the NT-3DFG substrates improve O 2 adsorption and catalytic activity while stabilizing the electrochemical activity during reactions. The FePc@NT-3DFG catalyst exhibits fast ORR kinetics with an extremely low Tafel slope of 28.3 ± 2.7 mV dec −1 , a higher half-wave potential of 0.911 ± 0.004 V (vs RHE), and notable long-term stability at 0.5 V (vs RHE) of 96.0 ± 0.4% retention after 30 h. Surface chemistry spectra validate electronic configuration modification of Fe at the iFOGs. Density functional theory calculations indicate that the extra layers of graphene improve oxygen adsorption. Moreover, additional exploration of other transition metal phthalocyanines supports the effects of iFOGs through the transition toward 4e − ORR. This work offers an expanded strategy for active site modification through edge-based graphene substrates for 4e − ORR.

X-ray absorption spectroscopy↗
DOE OSTI2024

Destabilizing high-capacity high entropy hydrides via earth abundant substitutions: From predictions to experimental validation

The vast chemical space of high entropy alloys (HEAs) makes trial-and-error experimental approaches for materials discovery intractable and often necessitates data-driven and/or first principles computational insights to successfully target materials with desired properties. In the context of materials discovery for hydrogen storage applications, a theoretical prediction-experimental validation approach can vastly accelerate the search for substitution strategies to destabilize high-capacity hydrides based on benchmark HEAs, e.g. TiVNbCr alloys. Here, in this study, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb) 75 Cr 25-x Fe x system. The as-prepared alloys crystallize in a single-phase bcc lattice for limited Fe content x < 7, while larger Fe content favors the formation of a secondary C14 Laves phase intermetallic. Short range order for alloys with x < 7 can be well described by a random distribution of atoms within the bcc lattice without lattice distortion. Hydrogen absorption experiments performed on selected alloys validate the predicted thermodynamic destabilization of the corresponding fcc hydrides and demonstrate promising lifecycle performance through reversible absorption/desorption. This demonstrates the potential of computationally expedited hydride discovery and points to further opportunities for optimizing bcc alloy ↔ fcc hydrides for practical hydrogen storage applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deep chalcogen donors and electron localization in aluminum nitride

Abstract We show with hybrid density functional theory calculations that chalcogen donors other than oxygen (i.e. S N , Se N , and Te N ) give rise to deep donor states in aluminum nitride. These donors trap a localized electron in their neutral charge state, leading to deep ( + / 0 ) donor levels that are 0.45 eV or more from the conduction-band edge. As such, this behavior is distinct from the DX behavior leads to deep (+/−) levels which affects other donors such as O N and Si Al . We highlight how these results hint at the formation of small electron polarons in AlN, which are found to be unstable in the bulk, but metastable when bound to donor dopants like Si Al and the chalocogens, with activation energies on the order of 0.2–0.3 eV. These results indicate that S, Se, and Te are not shallow donor dopants in aluminum nitride and identify origins of the experimentally observed ∼200–300 meV activation energies for dopant activation in donor-doped samples.

aluminum nitride↗
DOE OSTI2024

Enhancing Nitrogen Activation in Electrochemical Reduction: The Role of Rare Earth Oxide Surface Configurations

The pursuit of sustainable ammonia synthesis has prompted the exploration of ambient electrochemical nitrogen reduction reaction (e-NRR) as an alternative to the energy-intensive Haber-Bosch process. Here, this study conducted a theoretical investigation into the use of rare earth oxide materials, specifically dysprosium oxide (Dy 2 O 3 ), as potential electrocatalysts for NRR. Utilizing spin-polarized density functional theory calculations, we explored the interaction between Dy 2 O 3 surfaces and nitrogen (N 2 ) molecules, examining the capability of Dy 2 O 3 to adsorb and activate N 2 under ambient conditions. The results indicate that Dy 2 O 3 surfaces exhibit diverse configurations and bonding environments, providing a variety of reactive sites that display different behaviors in N 2 adsorption and activation. The distinctive electronic structure and surface chemistry of a particular Dy 2 O 3 surface configuration were found to significantly enhance the activation of N 2 by promoting charge transfer, which facilitates the NRR process. This research provides deep insights into the mechanistic pathways of N 2 reduction over Dy 2 O 3 , highlighting the surface properties as pivotal in catalysis. These theoretical insights serve as a foundation for the development of novel rare earth-based electrocatalytic materials for efficient ambient e-NRR, potentially transforming ammonia production into a greener and more energy-efficient process.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Soft x-ray spectroscopy of ammonium sulfate aerosols at varied pH conditions

Aqueous aerosols are useful model systems for understanding the physical chemistry that takes place in the Earth’s atmosphere, as well as the chemistry of other solar system objects. Ammonium sulfate aerosols are impacted by anthropogenic sources such as farming and shipping and are understood to be important seeds of cloud nucleation in the atmosphere, affecting climate. X-ray photoelectron spectroscopy and near edge x-ray absorption fine structure spectroscopy were used in tandem to probe the surface and bulk properties of aqueous (NH 4 ) 2 SO 4 , respectively. Aerosolized solutions of (NH 4 ) 2 SO 4 , some altered with NaOH to modify pH, were introduced via an aerodynamic lens to a velocity map imaging instrument for detection. The results show that as pH and sodium cation concentration increase, sulfate anion and ammonia become more prevalent at the surface, and the aerosol surface appears to become drier. However, most of the aerosol remains aqueous, with ion concentrations insensitive to the addition of NaOH. Density functional theory calculations were also performed to probe the coordination environment at the aerosol surfaces to understand the underlying phenomena. Finally, the changes seen on the surface of the aerosol underscore the importance of pH in regulating the structure and chemical properties of atmospherically relevant aerosols.

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

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Coexistence of Interacting Charge Density Waves in a Layered Semiconductor

Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe 4 , a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect FS nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that mono- and bi-layers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Finally, our findings not only offer microscopic insights into the correlated ground state of EuTe 4 but also provide a general non-equilibrium approach to understand coexisting, layer-dependent orders in a complex system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗
DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electron–electrophile coupled dinitrogen reduction in a cerium– meta -tetraphenolate system: a computational study

The use of lanthanide complexes for catalytic dinitrogen reduction is a new development in homogeneous catalysis. Density functional theory calculations on our recently reported cerium phenolate catalyst [K 2 Ce 2 (sol) 4 (mTP) 2 ] (mTP = {(OC 6 H 2 - 2- t Bu-4-Me) 2 CH} 2 - 1,3-C 6 H 4 ; sol = OMe 2 here; THF in the experiment) have been undertaken to elucidate the reduction, activation and silylation steps at the bound dinitrogen molecule, in the presence of the reductant, potassium metal (K 0 ) and the electrophile Me 3 SiCl (TMSCl). Out of the total of six electron reductions required to cleave the N 2 , the first two-electron reduction step was found to be highly disfavoured unless potassium cations (K + ) are included, upon which the step is rendered strongly exergonic; N–Si bond formation at the two-electron stage is predicted to be unfavourable. The three-electron-reduced N 2 -adduct is found to be at the reductive activation limit in the absence of added electrophiles, which can form N-element bonds and lower the overall charge. Added electron density beyond three-electron reduction no longer localises on N 2 , preventing formal N 2 4− formation. A pathway in which both K 0 and Me 3 SiCl work in concert was modelled, and six sequential reduction–silylation steps were calculated, showing how the N–N bond is cleaved after the third reduction, eventually releasing two equivalents of N(SiMe 3 ) 3 , and regenerating the starting complex with the highest barrier of any step being 22 kcal mol −1 . We establish alkali metal coordination and coupled electron–electrophile transfer as key factors in the design of rare-earth-mediated dinitrogen functionalisation.

Ahmad, Shahbaz [Univ. of Manchester (United Kingdo↗
DOE OSTI2025

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗
DOE OSTI2024

Generation of out-of-plane ferroelectric behavior in a one-atom-thick monolayer

Ferroelectricity with out-of-plane polarization has so far been found in several two-dimensional (2D) materials, including monolayers comprising three to five planes of atoms, e.g. α-In 2 Se 3 and MoTe 2 . Here, we explore the generation of out-of-plane polarization within a one-atom-thick monolayer material, namely hexagonal boron nitride. We performed density-functional-theory calculations to explore inducing ferroelectric-like distortions through incorporation of isovalent substitutional impurities that are larger than the host atoms. This disparity in bond lengths causes a buckling of the h-BN, either up or down, which amounts to a dipole with two equivalent energies and opposing orientations. We tested several impurities to explore the magnitude of the induced dipole and the switching energy barrier for dipole inversion. The effects of strain, dipole–dipole interactions, and vertical heterostructures with graphene are further explored. Our results suggest a highly-tunable system with ground state antiferroelectricity and metastable ferroelectricity. Finally, we expect that this work will help foster new ways to include functionality in layered 2D-material-based applications.

36 MATERIALS SCIENCE↗