Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium constant”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The In Situ Observation of Phase Transformations During Intercritical Annealing of a Medium Manganese Advanced High Strength Steel by High Energy X-Ray Diffraction

Microstructural changes during thermal processing of a medium manganese steel containing (in wt%) 0.19C and 4.39 Mn were evaluated in situ with a high energy X-ray diffraction system (HEXRD). Samples with an initial fully martensitic microstructure were heated to intercritical annealing (IA) temperatures of 600 or 650°C, held for 30 min, and cooled to room temperature. Diffraction data were analyzed to determine the variations in austenite and ferrite phase fractions and phase lattice constants throughout the ICA cycles. On heating, the 2 vol. pct of austenite present in the starting microstructure decomposed, and cementite precipitation then occurred. During isothermal holding, the austenite fraction increased, up to 20% for the sample annealed at 650°C. The measured austenite fractions were less than those calculated by Thermo-Calc for equilibrium conditions, indicating that the 30-min hold time was insufficient to achieve near-equilibrium conditions. Observed changes in lattice parameters during isothermal holding were interpreted to reflect composition changes due to redistribution of the C and Mn between austenite and ferrite. The results are discussed in relation to the potential for controlling austenite stability during ambient temperature deformation.

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

Molecular Weight Dependence of Block Copolymer Micelle Fragmentation Kinetics

The effect of molecular weight (M) on the fragmentation kinetics of micelles formed by 1,2-polybutadiene-block-poly(ethylene oxide) (PB-PEO) copolymers was studied in the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. A series of six samples, with total M ranging from 10 4 to 10 5 g mol –1 and nearly constant composition (f PEO ≈ 0.4), were examined; all six formed spherical micelles with PEO coronas. Nonequilibrium PB-PEO micelles were prepared by direct dissolution, a process that systematically produces nanoparticles with mean aggregation numbers more than twice the equilibrium values. When subjected to high temperature annealing (170 °C), the average micelle radius was found to decrease substantially, as determined by temperature-jump dynamic light scattering (T-jump DLS) and time-resolved small-angle X-ray scattering (TR-SAXS). Furthermore, the characteristic fragmentation times (τ) were found to increase strongly with increasing degree of polymerization N, as τ ~ N 1.8 . This result compares favorably with the prediction of a previously untested model.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Computing Free Energies with Fluctuation Relations on Quantum Computers

One of the most promising applications for quantum computers is the dynamic simulation of quantum materials. Current hardware, however, sets stringent limitations on how long such simulations can run before decoherence begins to corrupt results. The Jarzynski equality, a fluctuation theorem that allows for the computation of equilibrium free energy differences from an ensemble of short, non-equilibrium dynamics simulations, can make use of such short-time simulations on quantum computers. Here, we present a quantum algorithm based on the Jarzynski equality for computing free energies of quantum materials. We demonstrate our algorithm using the transverse field Ising model on both a quantum simulator and real quantum hardware. As the free energy is a central thermodynamic property that allows one to compute virtually any equilibrium property of a physical system, the ability to perform this algorithm for larger quantum systems in the future has implications for a wide range of applications including the construction of phase diagrams, prediction of transport properties and reaction constants, and computer-aided drug design.

Bassman, Lindsay↗
DOE OSTI2022

Short- vs . long-range elastic distortion: structural dynamics of a [2 × 2] tetrairon( II ) spin crossover grid complex observed by time-resolved X-Ray crystallography

Spin crossover complexes (SCO) are among the most studied molecular switches due to their potential use in displays, sensors, actuators and memory components. A prerequisite to using these materials is the understanding of the structural changes following the spin transition at out-of-equilibrium conditions. So far, out-of-equilibrium studies in SCO solids have been focused on mononuclear complexes, though a growing number of oligonuclear SCO complexes showing cooperative effects are being reported. Here, we use time-resolved pink Laue crystallography to study the out-of-equilibrium dynamics of a [2 × 2] tetranuclear metallogrid of the form [Fe II 4 L Me 4 ](BF 4 ) 4 ·2MeCN ([L Me ] - = 4-methyl-3,5-bis{6-(2,2'-bipyridyl)}pyrazolate). The out-of-equilibrium spin state switching induced by a ps laser pulse demonstrates that the metallogrid exhibits a multi-step response similar to that reported for mononuclear complexes. Contrary to the mononuclear complexes, the metallogrid shows two types of elastic distortions at different time scales. The first is a short-range distortion that propagates over the entire Fe 4 grid complex during the ps time scale, and it is caused by the rearrangement of the coordination sphere of the photo-switching ion and the constant feedback between strongly linked metal ions. The second is a long-range distortion caused by the anisotropic expansion of the lattice during the ns time scale, observed in mononuclear materials. The structural analysis demonstrates that the long-range prevails over the short-range distortion, inducing the largest deformation of both the entire grid and the coordination sphere of each metal ion. The present study sheds light on the out-of -equilibrium dynamics of a non-cooperative oligonuclear complex.

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

Improved Accuracy in Semi-Experimental Structure Determination by Resolving Problems Associated with Rotation of Principal Inertial Axes of Isotopologues: Structures of 1,3-Oxazole ( c -C 3 H 3 NO)

The rotational spectrum of the normal isotopologue of 1,3-oxazole (c-C 3 H 3 NO) was observed from 43 to 750 GHz. Over 3900 transitions for the ground vibrational state are measured, assigned, and least-squares fit to sextic centrifugally distorted-rotor Hamiltonians. The measured frequencies and resulting spectroscopic constants from this extended spectral range, combined with previous measurements of the nuclear quadrupole coupling constants, will facilitate astronomical searches for oxazole across the majority of the range of modern radiotelescopes. Spectra for a set of 30 oxazole isotopologues, which include multiple isotopic substitutions of each atom, are used to determine the first semi-experimental equilibrium ($r$$^{SE}_{e}$) structure and semi-experimental substitution structure ($r$$^{SE}_{e}$), each using CCSD(T) computed values for the vibration–rotation interaction and electron-mass corrections. The large number of isotopologues, including 21 isotopologues observed for the first time, and the redundant substitutions of each atom provide sufficient spectroscopic information to determine the $r$$^{SE}_{e}$ structure with the expected high level of accuracy and precision (0.0001 or 0.0002 Å in bond distances and 0.013 to 0.025° in bond angles). In the course of this study, we analyzed a known issue for some $r$$^{SE}_{e}$ structure determinations of near-oblate asymmetric tops in which inclusion of individual isotopologues degrades the structure determination. We demonstrate that this problem primarily arises from the difference in the values of the computed vibration–rotation interaction corrections as evaluated at the computed re geometry vs the $r$$^{SE}_{e}$ geometry of the “real” molecule. Our solution to this problem substantially improves the $r$$^{SE}_{e}$ structure of oxazole and likely can be generalized to many other molecules.

Chemical structure↗
DOE OSTI2021

Fluctuation–dissipation relation for a quantum Brownian oscillator in a parametrically squeezed thermal field

In this paper we study the nonequilibrium evolution of a quantum Brownian oscillator, modeling the internal degree of freedom of a harmonic atom or an Unruh–DeWitt detector, coupled to a nonequilibrium and nonstationary quantum field bath and inquire whether a fluctuation–dissipation relation (FDR) can exist after/if it approaches equilibration. This is a nontrivial issue because a squeezed field bath cannot reach equilibration and yet, as this work shows, the system oscillator indeed can, which is a necessary condition for FDRs. We discuss three different settings: (A) The bath field essentially remains in a squeezed thermal state throughout, whose squeeze parameter is a mode- and time-independent constant. This situation is often encountered in quantum optics and quantum thermodynamics. (B) The bath field is initially in a thermal state, but is subjected to a parametric process leading to mode- and time-dependent squeezing. This scenario is encountered in cosmology and dynamical Casimir effects. The squeezing in the bath in both types of processes will affect the oscillator’s nonequilibrium evolution. We show that at late times it approaches equilibration and this stationarity condition warrants the existence of a FDR. The trait of squeezing is marked by the oscillator’s effective equilibrium temperature, and the proportionality factor in the FDR is only related to the stationary component of the noise kernel of the bath field. Setting (C) is more subtle: A finite system–bath coupling strength can set the oscillator in a squeezed state even though the bath field is stationary and does not engage in any parametric process. The squeezing of the system in this case is in general time-dependent but becomes constant when the internal dynamics is fully relaxed. We begin with comments on the broad range of physical processes involving squeezed thermal baths and end with some remarks on the significance of FDRs in capturing the essence of quantum backreaction in nonequilibrium and stochastic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Straintronics with single-layer MoS 2 : A quantum Monte Carlo study

Using state-of-the-art quantum Monte Carlo (QMC) methods, we study straintronic properties of a single MoS 2 monolayer. 2D MoS 2 is a quintessential straintronic material for which many experiments have been performed. First, we determine the equilibrium atomic structure which is not known experimentally and is strictly needed to correctly determine the straintronic properties. That enables us to precisely analyze the quasiparticle band gaps for any applied biaxial strain, which we describe by a bivariate paraboloid function of lattice constant and internal structural parameter. Using the fixed-node QMC calculations fitted by analytical formulas, we localize the following excited state crossings between the direct, K → K , and indirect Γ → X and K → K /2 excitations. Based on this highly accurate many-body treatment, we predict a gauge factor of 136 meV/% for the K → K transition and a fairly narrow window of ≈ 2.8 % from compressive to mildly tensile strains, accounting for only ≈ 0.3 eV band gap change maintaining the direct character of the gap. Consequently, we suggest that, compared to other 2D materials, such as phosphorene, there is only a limited straintronic tuneability in this material often studied for its straintronic properties. QMC results are compared to results of standard DFT modeling, which reveal insights into the corresponding inaccuracies and therefore open a window for educated use of rapid DFT approaches. Published by the American Physical Society 2024

Huang, Y.↗
DOE OSTI2023

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Temperature-Dependent Sn Incorporation and Defect Formation in Pseudomorphic SiSn Layers on Si (001) via Molecular Beam Epitaxy

SiSn alloys have attracted growing interest for group-IV bandgap engineering, although their epitaxial growth remains challenging due to the extremely low equilibrium solubility of Sn in Si. In this work, fully strained (pseudomorphic) SiSn epitaxial layers were grown on Si (001) substrates by means of molecular beam epitaxy. A systematic investigation reveals a strong inverse correlation between growth temperature and Sn incorporation efficiency. Despite a constant Sn flux, the incorporated Sn composition decreases from 5.5% to 3.2% as the growth temperature increases, indicating a pronounced temperature dependence of Sn incorporation. Reflection high-energy electron diffraction indicates a gradual transition of the growth from two-dimensional to three-dimensional with increasing film thickness. Structural characterization by means of X-ray diffraction, atomic force microscopy, and transmission electron microscopy confirms the pseudomorphic growth and smooth surface morphology and reveals twins and stacking faults near the surface region. These results establish a quantitative reference for SiSn growth kinetics and provide guidance for future studies of SiSn and SiGeSn alloys in silicon-compatible electronic and optoelectronic applications.

SiSn alloys↗
DOE OSTI2023

Plasma catalysis: separating plasma and surface contributions for an Ar/N 2 /O 2 atmospheric discharge interacting with a Pt catalyst

Atmospheric pressure non-equilibrium plasmas can form nitrogen oxide (NO x ) compounds directly from nitrogen and oxygen without a catalyst, and at lower catalyst temperatures than would be possible without plasma. In this work, the oxidation of plasma-produced NO from an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet (APPJ) over a platinum-on-alumina powder catalyst was investigated with in-situ infrared spectroscopy. Products downstream of the catalyst bed were analyzed along with catalyst surface species. The catalyst was exposed to plasma at both constant temperature and a cyclic temperature ramp in order to study long-lasting and transient surface changes. Primary incident reactive species to the catalyst were assessed to be NO and O 3 . Pt-Al 2 O 3 at 350 °C increased oxidation of NO relative to Al 2 O 3 or an empty chamber. The surface state of Pt-Al 2 O 3 evolves during plasma-effluent exposure and requires upwards of 20 min exposure for stabilization compared to Al 2 O 3 . Once stable surface conditions are achieved, thermal cycling reveals a repeatable hysteresis pattern in downstream products. At low temperature, oxygen and NO x accumulate on the catalyst surface and react at elevated temperatures to form NO 2 . Increasing plasma power and O 2 :N 2 ratio increases the hysteresis of the heating relative to the cooling curves in the pattern of NO 2 formation. The limitation on NO oxidation at high temperatures was assessed to be Pt-O which is depleted as the catalyst is heated. Once stored species have been depleted, NO oxidation rates are determined by incoming reactants. Two overlapping NO oxidation patterns are identified, one determined by surface reactants formed at low temperature, and the other by reactants arriving at the surface at high temperature. As a result, the plasma is responsible for providing the reactants to the catalyst surface, while the catalyst enables reaction at high temperature or storage at low temperature for subsequent reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Water and Solute Activities Regulate CO 2 Reduction in Gas-Diffusion Electrodes

Electrolysis of CO 2 at gas-diffusion electrodes (GDEs) has typically been limited by the supply of gas to the electrocatalyst, overshadowing the importance of the supply of water. However, at high current densities that approach 1 A cm –2 , where the electrolyte becomes highly concentrated in the catalyst layer of a GDE, the activity of water and solutes deviate from their bulk dilute solution values, potentially slowing reaction rates and changing reaction equilibrium potentials. In addition, as flow plates for the gas stream are introduced to enable larger electrodes and high single pass conversion of CO 2 to product, variations in the gas composition will become important. By drawing upon literature for the oxygen reduction reaction (ORR), here we explain how to account for these effects in future modeling and experimental work, with particular attention to accurate use of the Nernst equation for electrode potentials and the Arrhenius equation for reaction rates. Specifically, using measurements of KOH solvent and solute activity reported in literature, and assuming the second protonation of CO 2 by water as the rate-determining step, we show the Nernst equation dilute-solution approximation of the CO 2 to CO equilibrium potential to be accurate below 5 M KOH, but it has a 74 mV error when increasing the concentration up to 10 M KOH. Finally, a simple one-dimensional model of a serpentine flow-field on a GDE demonstrated that a reactor with constant pressure of 1 bar and 1 A cm –2 at the inlet had only ~0.3 A cm –2 at the outlet for a conversion in CO 2 partial pressure from 0.90 to 0.48 bar, showing the significant practical implications of this work.

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

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

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

Guiding Self-Assembly of Active Colloids by Temporal Modulation of Activity

Self-organization phenomena in ensembles of self-propelled particles open pathways to the synthesis of new dynamic states not accessible by traditional equilibrium processes. The challenge is to develop a set of principles that facilitate the control and manipulation of emergent active states. In this study, we report that dielectric rolling colloids energized by a pulsating electric field self-organize into alternating square lattices with a lattice constant controlled by the parameters of the field. We combine experiments and simulations to examine spatiotemporal properties of the emergent collective patterns and investigate the underlying dynamics of the self-organization.We reveal the resistance of the dynamic lattices to compression and expansion stresses leading to a hysteretic behavior of the lattice constant. The general mechanism of pattern synthesis and control in active ensembles via temporal modulation of activity can be applied to other active colloidal systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Forming mechanism of equilibrium and non-equilibrium metallurgical phases in dissimilar aluminum/steel (Al–Fe) joints

Abstract Forming metallurgical phases has a critical impact on the performance of dissimilar materials joints. Here, we shed light on the forming mechanism of equilibrium and non-equilibrium intermetallic compounds (IMCs) in dissimilar aluminum/steel joints with respect to processing history (e.g., the pressure and temperature profiles) and chemical composition, where the knowledge of free energy and atomic diffusion in the Al–Fe system was taken from first-principles phonon calculations and data available in the literature. We found that the metastable and ductile (judged by the presently predicted elastic constants) Al 6 Fe is a pressure ( P ) favored IMC observed in processes involving high pressures. The MoSi 2 -type Al 2 Fe is brittle and a strong P -favored IMC observed at high pressures. The stable, brittle η-Al 5 Fe 2 is the most observed IMC (followed by θ-Al 13 Fe 4 ) in almost all processes, such as fusion/solid-state welding and additive manufacturing (AM), since η-Al 5 Fe 2 is temperature-favored, possessing high thermodynamic driving force of formation and the fastest atomic diffusivity among all Al–Fe IMCs. Notably, the ductile AlFe 3 , the less ductile AlFe, and most of the other IMCs can be formed during AM, making AM a superior process to achieve desired IMCs in dissimilar materials. In addition, the unknown configurations of Al 2 Fe and Al 5 Fe 2 were also examined by machine learning based datamining together with first-principles verifications and structure predictions. All the IMCs that are not P- favored can be identified using the conventional equilibrium phase diagram and the Scheil-Gulliver non-equilibrium simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of microstructure on elastic properties in the alloy Ti-42Al-8.5 Nb

The microstructure of a γ-TiAl alloy containing niobium undergoes continuous transformations during annealing at 550 °C. These take place within the α 2 -phase of lamellar (α 2 +γ) colonies. We have characterized these changes by transmission electron microscopy and made a correlation to the elastic constants determined by resonant ultrasonic spectroscopy during in-situ annealing at the same temperature. Further, the results show a continuous increase of the E-modulus. This can be attributed to elastic strains acting within α 2 lamellae as a result of lattice transformations. After 5000h, a thermodynamic equilibrium was established with equal fractions of α 2 -phase and O-phase that differ only in lattice symmetry and Nb content. Generally, the presence of O-phase should contribute to lower elastic moduli, however in the present case the evolving microstructural characteristics were responsible for the observed increase in the E-modulus.

36 MATERIALS SCIENCE↗