SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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.

SourceResearch connectionFrom yearThrough year

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

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Extended Barton–Bandis model for rock joints under cyclic loading: Formulation and implicit algorithm

In this paper, the Barton–Bandis model for rock joints is extended to cyclic loading conditions, without any new material parameter. Also developed herein is an implicit solution algorithm for the extended Barton–Bandis model, which can also be used for the original Barton–Bandis model for which an implicit algorithm has been unavailable. To this end, we first cast the Barton–Bandis model into an incremental elasto-plastic framework, deriving an expression for the elastic shear stiffness being consistent with the original model formulation. We then extend the model formulation to cyclic loading conditions, incorporating the dependence of shear stress and dilation on the joint position and the shearing direction. The extension is achieved by introducing a few state-dependent variables which can be calculated with the existing material parameters. For robust and accurate utilization of the model, we also develop an implicit algorithm based on return mapping, which is unconditionally stable and guarantees the satisfaction of the strength criterion. We verify that the proposed model formulation and algorithm produce virtually the same results as the original Barton–Bandis model under monotonic shearing conditions. Here, we then validate the extended Barton–Bandis model against experimental data on natural rock joints under cycling loading conditions. The present work thus enables the Barton–Bandis model, which has been exceptionally popular in research and practice, to be applicable to a wider range of problems in rock mechanics and rock engineering.

58 GEOSCIENCES↗

DOE OSTI2023

NMR relaxometry characterization of water adsorption in corn stover anatomical fractions

Nuclear magnetic resonance (NMR) relaxometry is applied to provide direct measurement of water adsorption in anatomical fractions of corn stover. NMR transverse T 2 relaxation time distribution measurements indicate multiple water populations, which vary with anatomical fraction and water adsorption. Measured T 2 data are used to calculate thermodynamic properties of Brunauer-Emmet-Teller adsorption theory using a model to estimate mono and bilayer relaxation. T 2 data are used directly to determine rotational diffusion correlation times indicating adsorption interaction strength. T 1 -T 2 longitudinal-transverse relaxation time correlation measurements quantify differences in the molecular level structural order of the adsorbate surface water as a function of water activity, i.e. relative humidity or water vapor partial pressure. The T 1 /T 2 ratio provides a measure of the surface energy related to the adsorption strength and surface diffusive mobility of the water adsorbate and differentiates the anatomical fractions. Furthermore, the results indicate that direct measurement of NMR relaxation times can be used to characterize corn stover biomass water adsorption, data relevant to processing and handling considerations.

09 BIOMASS FUELS↗

DOE OSTI2022

Increasing the composition range of a novel τ 11 -Al 4 Fe 1.7 Si alloy with additions of Mn

Automotive applications need low-cost, lightweight, high-temperature alloys to increase vehicle efficiency. The Al–Fe–Si system provides an opportunity to develop such a material, as it consists of three low-cost elements that are all abundant in nature. Specifically, the τ 11 -Al 4 Fe 1.7 Si ternary intermetallic phase is a high-temperature, lightweight phase with high strength and good corrosion resistance. However, this phase exhibits a narrow compositional range of stability, resulting in undesirable microstructures forming during solidification and processing, limiting its use in potential applications. Density functional theory (DFT) calculations and a thermodynamically-driven experimental approach utilizing diffusion couples were employed to study the effect of Mn on the stability and composition range of τ 11 -Al 4 Fe 1.7 Si. The DFT calculations showed a decrease in the energy of the structure when alloying with Mn. Experimental results confirmed the predictions from the DFT calculations, indicating that alloying with Mn increases the compositional range, and thus the processability of this phase. New phase diagrams and equilibria are proposed by exploring and determining phase boundaries for the τ 11 -Al 4 Fe 1.7 Si phase with Mn.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Identification of Active Metal Carbide and Nitride Catalytic Facets for Hydrodeoxygenation Reactions

The catalytic hydrodeoxygenation (HDO) reaction is of considerable interest for biomass conversion to valuable chemicals and fuels, where one of the critical bottlenecks is the lack of cost-effective and efficient catalysts. To discover cost-efficient catalysts for the HDO reaction, we employed a density functional theory-based hierarchical catalyst design strategy based on catalytic descriptors, reaction energy profiles, and microkinetic modeling (MKM). We focused on the carbide and nitride catalyst space, for which we calculated 121 catalyst surfaces of Mo 2 C, MoC, Mo 2 N, W 2 C, NbC, VC, VN, and NbN catalysts. Based on the computed surface energies, reaction energies of oxygen removal, carbon binding strength, and the surface area of nanoparticles, the likely active facets are the Mo 2 C(111), MoC(011), VN(100), Mo 2 N(001), Mo 2 N(011), and Mo 2 N(100) surfaces. Further, detailed energy profiles were obtained, and MKM was performed for a model reaction (glycolaldehyde + 2H 2 . ethylene + 2H 2 O) on the Mo 2 C(111), VN(100), and MoC(100) surfaces. Based on the computed volcano map obtained from MKM, the predicted active facets for this HDO reaction are the Mo 2 C(111), MoC(011), VN(011), Mo 2 N(001), Mo 2 N(011), and Mo 2 N(100) surfaces. Additionally, none of the carbide and nitride catalyst surfaces are located in the optimal catalytic activity part. Therefore, it is essential to modify the catalyst via adding dopants or alloying to improve the catalytic activity. Catalytic modifications that can destabilize the surface adsorption of O*/H 2 O* and decrease the energy barriers of O-H bond formation are recommended to facilitate the HDO on the carbide and nitride catalysts. These a priori investigations provide guidelines for future low-cost HDO catalyst development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Interference-Limited Absorption in Dense Molecular Nanolayers Near Reflecting Surfaces

We investigate linear resonant absorption by a dense ensemble of molecules confined to a sub-wavelength layer in two geometries: (i) a free-standing film in a homogeneous space and (ii) the same film placed at a controlled distance from a reflecting surface. In both cases, increasing the effective light–matter coupling (via molecular density/oscillator strength) produces a nonmonotonic response: absorption rises to an optimum and then decreases as the film becomes increasingly radiatively bright and reflective. Finite-difference time-domain simulations and analytical transfer-matrix calculations agree quantitatively and yield compact ridge conditions for the optimum. We interpret the trends using a scattering/port picture: the isolated film is a symmetric two-port system (reflection and transmission), which bounds single-sided resonant absorption to ≤50% in the ultrathin limit (reflecting transition saturation), whereas adding a mirror suppresses transmission and converts the structure into an effectively one-port absorber. In the mirror-backed geometry, interference can cancel reflection, and unity absorption is obtained at critical coupling, when radiative leakage is balanced by intrinsic molecular loss. These results clarify fundamental limits and design rules for collective absorption in dense molecular layers near dielectric or metallic boundaries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach

We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al 2020 Phys. Plasmas 27 070701; Loizu and Bonfiglio 2023 J. Plasma Phys. 89 905890507), we make use of Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502), an equilibrium solver based on the variational principle of the multi-region relaxed magnetohydrodynamics (MHDs), featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model J bs = C$\boldsymbol{\nabla}$p to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter Δ' and bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHDs code (Huang and Bhattacharjee 2016 Astrophys. J. 818 20) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical modified Rutherford equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and NTMs, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Shape fluctuations and optical transition of He 2 * excimer tracers in superfluid He 4

Metastable He$_{2}^{*}$ excimer molecules have been utilized as tracer particles of the normal component in superfluid $^4$He (He II) which can be imaged via laser-induced fluorescence. These excimer molecules form tiny bubbles in He II and can bind to quantized vortices at sufficiently low temperatures, thereby allowing for direct visualization of vortex dynamics in an inviscid superfluid. However, the $a^{3}\Sigma^+_{u}$${\rightarrow}$$c^{3}\Sigma^+_{g}$ optical absorption line, which is responsible for the fluorescence imaging of the He$_{2}^{*}$ molecules, is controlled by fluctuations on the bubble shape, and its exact line profile is not known at low temperatures. In this paper, we present a bubble model for evaluating the surface fluctuation eigenmodes of the excimers in He II. The line profile of the $a^{3}\Sigma^+_{u}{\rightarrow}c^{3}\Sigma^+_{g}$ transition is calculated at different temperatures by considering both the zero-point and thermal fluctuations on the bubble shape. We show that, as the temperature drops from 2~K to 20 mK, the peak absorption strength is enhanced by a factor of about five, accompanying a blueshift of the peak location by about 2 nm. A double-peak line profile due to the rotational levels of the molecular core can be resolved. This bubble model also allows us to evaluate the stiffness of the He$_{2}^{*}$ bubbles and hence their diffusion constant in He II due to scattering off thermal phonons. Our results will aid the design of future experiments on imaging quantized vortices in He II using He$_{2}^{*}$ tracers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Exploring the Dirac nature of RbBi 2

Characteristics of topological semimetals such as a nonsaturating magnetoresistance (MR), a field-induced metal to semiconducting crossover and a robust resistivity plateau are observed under a magnetic field in type-I RbBi 2 bulk superconductor with $\mathrm{T}$ c = 4.15 K. The MR exhibits a notable 3500% increase at 2 K and 9 T and the resistivity follows a power law temperature dependence, while the MR ∝ $\mathrm{H}$ 1.26 , indicating weak carrier compensation. Further, first principles calculations provided insights into the dynamical stability of the cubic structure at 0 K. Both hole and electron pockets are observed at the Fermi surface. The electron-phonon interaction constant indicates weak coupling strength (<1) that leads to a maximum predicted $\mathrm{T}$ c of 2.852 K. Just below the Fermi level, E F , the electronic band structure consists of linear band crossings at the $\mathrm{X}$ points in the Brillouin zone (BZ) corresponding to massless, symmetry-protected Dirac fermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Lowest-energy broad 𝛼-cluster resonances in 19 F

There is strong astrophysical interest in the structure of 19 F states near the α-decay threshold, as their properties are important for understanding the development of α clustering in the 20 Ne region. The emergence of clustered states, and more generally of states that couple strongly to reaction channels near their decay thresholds, is also a topic of current interest in theoretical nuclear physics. In this work, we determine the parameters of broad low-spin states in 19 F near the α-decay threshold and present a theoretical study of these states. The analysis is restricted to ℓ = 0 and 1 resonances in the α + 15 N system close to the α threshold in 19 F. Excitation functions for 15 N(α,α) elastic scattering were measured using the Thick Target Inverse Kinematics (TTIK) method, and these new data, together with older high–energy-resolution measurements, were analyzed within the R-matrix framework. The nuclear structure of 19 F was calculated using configuration-interaction methods with a recently developed effective interaction Hamiltonian. As a result, we identify a sequence of α-cluster resonances in 19 F and map the distribution of clustering strength, which is relevant for astrophysical reaction modeling and for the theoretical understanding of many-body dynamics and the emergence of clustering in loosely bound or unstable nuclei. Furthermore, the work advances theoretical insight into the origins of clustering and highlights open questions for future theoretical and experimental studies.

19F↗

DOE OSTI2021

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Surpassing the Performance of Phenolate-derived Ionic Liquids in CO 2 Chemisorption by Harnessing the Robust Nature of Pyrazolonates

Superbase-derived ionic liquids (SILs) are promising sorbents to tackle the carbon challenge featured by tunable interaction strength with CO 2 via structural engineering, particularly the oxygenate-derived counterparts (e. g., phenolate). However, for the widely deployed phenolate-derived SILs, unsolved stability issues severely limited their applications leading to unfavorable and diminished CO 2 chemisorption performance caused by ylide formation-involved side reactions and the phenolate-quinone transformation via auto-oxidation. Here, in this work, robust pyrazolonate-derived SILs possessing anti-oxidation nature were developed by introducing aza-fused rings in the oxygenate-derived anions, which delivered promising and tunable CO 2 uptake capacity surpassing the phenolate-based SIL via a carbonate formation pathway (O-C bond formation), as illustrated by detailed spectroscopy studies. Further theoretical calculations and experimental comparisons demonstrated the more favorable reaction enthalpy and improved anti-oxidation properties of the pyrazolonate-derived SILs compared with phenolate anions. The achievements being made in this work provides a promising approach to achieve efficient carbon capture by combining the benefits of strong interaction strength of oxygenate species with CO 2 and the stability improvement enabled by aza-fused rings introduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Interlayer pairing in bilayer nickelates

The discovery of T c ~ 80 K superconductivity in pressurized La 3 Ni 2 O 7 has launched a new platform to study high-temperature superconductivity. Using non-perturbative dynamic cluster approximation quantum Monte Carlo calculations, we characterize the magnetic and superconducting pairing behavior of a realistic bilayer two-orbital Hubbard-Hund model of this system that describes the relevant Ni e g states with physically relevant interaction strengths. We find a leading s ± superconducting instability in this model at a temperature T ~ 100 K close to the experimentally observed T c . Analyzing the orbital and spatial structure of the effective pairing interaction giving rise to this state reveals that the interaction predominantly acts between local interlayer pairs of the $d_{3z^2-r^2}$ orbital. By correlating the strength of the interaction with that of the magnetic spin fluctuations we show that it is driven by strong interlayer spin-fluctuations arising from the $d_{3z^2-r^2}$ orbital. These results provide first-time non-perturbative evidence supporting the picture that a simple single-orbital bilayer Hubbard model for the Ni $d_{3z^2-r^2}$ orbital provides an excellent low-energy effective description of the superconducting behavior of La 3 Ni 2 O 7 .

Maier, Thomas A. [Oak Ridge National Laboratory (O↗

DOE OSTI2024

β -decay feeding intensity distribution of Mn 64

Nuclei around the N=40 “Island of Inversion” exhibit interesting structure features that have been the focus of several experimental and theoretical studies. Here, the present work presents the first complete study of the β-decay feeding intensity distribution and Gamow-Teller distribution for the β decay of 64 Mn to 64 Fe up to ~10 MeV. The β-decay intensity function was extracted from Total Absorption Spectroscopy measurements made at the National Superconducting Cyclotron Laboratory with the Summing NaI(Tl) (SuN) detector. The experimental results are compared to shell model calculations with and without the inclusion of the νg 9/2 orbital. From this comparison it is clear that the νg 9/2 orbital is essential for the accurate description of the 64 Fe β-decay strength above ~3 MeV, emphasizing once again the transitional nature of this nucleus into the N=40 Island of Inversion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Higher-order symmetry plane correlations in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The correlations between event-by-event fluctuations of symmetry planes are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the Large Hadron Collider. This analysis is conducted using the Gaussian estimator technique, which is insensitive to biases from correlations between different flow amplitudes. The study presents, for the first time, the centrality dependence of correlations involving up to five different symmetry planes. The correlation strength varies depending on the harmonic order of the symmetry plane and the collision centrality. Comparisons with measurements from lower energies indicate no significant differences within uncertainties. Additionally, the results are compared with hydrodynamic model calculations. Although the model predictions provide a qualitative explanation of the experimental results, they overestimate the data for some observables. This is particularly true for correlators that are sensitive to the nonlinear response of the medium to initial-state anisotropies in the collision system. As these new correlators provide unique information—independent of flow amplitudes—their usage in future model developments can further constrain the properties of the strongly interacting matter created in ultrarelativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗

DOE OSTI2026

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong↗

DOE OSTI2020

Topological Hall Effect in a Topological Insulator Interfaced with a Magnetic Insulator

A topological insulator (TI) interfaced with a magnetic insulator (MI) may host an anomalous Hall effect (AHE), a quantum AHE, and a topological Hall effect (THE). Recent studies, however, suggest that coexisting magnetic phases in TI/MI heterostructures may result in an AHE-associated response that resembles a THE but in fact is not. Here, we report a genuine THE in a TI/MI structure that has only one magnetic phase. The structure shows a THE in the temperature range of T = 2–3 K and an AHE at T = 80–300 K. Over T = 3–80 K, the two effects coexist but show opposite temperature dependencies. Control measurements, calculations, and simulations together suggest that the observed THE originates from skyrmions, rather than the coexistence of two AHE responses. The skyrmions are formed due to a Dzyaloshinskii–Moriya interaction (DMI) at the interface; the DMI strength estimated is substantially higher than that in heavy metal-based systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Local Chemical Clustering Enabled Ultrahigh Capacitive Energy Storage in Pb-Free Relaxors

Designing Pb-free relaxors with both a high capacitive energy density (W rec ) and high storage efficiency (η) remains a remarkable challenge for cutting-edge pulsed power technologies. Local compositional heterogeneity is crucial for achieving complex polar structure in solid solution relaxors, but its role in optimizing energy storage properties is often overlooked. Here, in this paper, we report that an exceptionally high W rec of 15.2 J cm –3 along with an ultrahigh η of 91% can be achieved through designing local chemical clustering in Bi 0.5 Na 0.5 TiO 3 –BaTiO 3 -based relaxors. A three-dimensional atomistic model derived from neutron/X-ray total scattering combined with reverse Monte Carlo method reveals the presence of subnanometer scale clustering of Bi, Na, and Ba, which host differentiated polar displacements, and confirming the prediction by density functional theory calculations. This leads to a polar state with small polar clusters and strong length and direction fluctuations in unit-cell polar vectors, thus manifesting improved high-field polarizability, steadily reduced hysteresis, and high breakdown strength macroscopically. The favorable polar structure features also result in a unique field-increased η, excellent stability, and superior discharge capacity. Our work demonstrates that the hidden local chemical order exerts a significant impact on the polarization characteristic of relaxors, and can be exploited for accessing superior energy storage performance.

25 ENERGY STORAGE↗