SEARCH · Search NASA
Results for “HYSTERESIS”
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.
Deconvolution of electroosmotic flow in hysteresis ion transport through single asymmetric nanopipettes
This report establishes a universal method to diagnose and deconvolute electroosmotic flow in nanoscale transport processes through current–potential measurements and analysis without simulation.
Energy balance closure and evapotranspiration hysteresis in central Amazon under contrasting conditions during the GoAmazon project in 2014 and 2015
Explore the source record for details and available documents.
A Sorption-Induced Swelling Model for Sorption Hysteresis in Glassy Polyurethane Foam
Explore the source record for details and available documents.
Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)
Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.
Experimental discovery of structure–property relationships in ferroelectric materials via active learning
Emergent functionalities of structural and topological defects in ferroelectric materials underpin an extremely broad spectrum of applications ranging from domain wall electronics to high dielectric and electromechanical responses. Many of these functionalities have been discovered and quantified via local scanning probe microscopy methods. However, the search has until now been based on either trial and error, or using auxiliary information such as the topography or domain wall structure to identify potential objects of interest on the basis of the intuition of operator or pre-existing hypotheses, with subsequent manual exploration. Here we report the development and implementation of a machine learning framework that actively discovers relationships between local domain structure and polarization-switching characteristics in ferroelectric materials encoded in the hysteresis loop. The hysteresis loops and their scalar descriptors such as nucleation bias, coercive bias and the hysteresis loop area (or more complex functionals of hysteresis loop shape) and corresponding uncertainties are used to guide the discovery of these relationships via automated piezoresponse force microscopy and spectroscopy experiments. As such, this approach combines the power of machine learning methods to learn the correlative relationships between high-dimensional data, as well as human-based physics insights encoded into the acquisition function. For ferroelectric materials, this automated workflow demonstrates that the discovery path and sampling points of on- and off-field hysteresis loops are largely different, indicating that on- and off-field hysteresis loops are dominated by different mechanisms. Here, the proposed approach is universal and can be applied to a broad range of modern imaging and spectroscopy methods ranging from other scanning probe microscopy modalities to electron microscopy and chemical imaging.
Seasonal changes in the drivers of water physico-chemistry variability of a small freshwater tidal river
Where rivers meet the sea, tides can exert a physical and chemical influence on the lower reaches of a river. How tidal dynamics in these tidal river reaches interact with upstream hydrological drivers such as storm rainfall, which ultimately determines the quantity and composition of material transferred from watersheds to estuaries, is currently unknown. We monitored a small freshwater tidal river in the Pacific Northwest, USA in high resolution over one year to evaluate the relative importance of tides versus upstream hydrological flows (i.e., base flow and precipitation events) on basic physico-chemical parameters (pH, dissolved oxygen, turbidity, specific conductivity, and temperature), and how these interactions relate to the downstream estuary. Tidal variability and diurnal cycles (i.e. solar radiation) dominated water physico-chemical variability in the summer, but the influence of these drivers was overshadowed by storm-driven sharp pulses in river physico-chemistry during the remainder of the year. Within such events, we found incidences of counterclockwise hysteresis of pH, counterclockwise hysteresis of dissolved oxygen, and clockwise hysteresis of turbidity, although systematic trends were not observed across events. The dominance of storm rainfall in the river’s physico-chemistry dynamics, and similar pulses of decreased pH observed in adjacent estuarine waters, suggest that the linkage between tidal streams and the broader system is variable throughout the year. High-frequency monitoring of tidal river biogeochemistry is therefore crucial to enable the assessment of how the relative strength of these drivers may change with future sea level rise and altered precipitation patterns to modulate biogeochemical dynamics across the land-ocean-atmosphere continuum.
Antiferroelectric Ceramics for Energy–Efficient Capacitors by Theory–Guided Discovery
Antiferroelectric ceramics, via the electric-field-induced antiferroelectric (AFE)–ferroelectric (FE) phase transitions, show great promise for high-energy-density capacitors. Yet, currently, only 70–80% energy release is found during a charge–discharge cycle. Here, for PbZrO 3 -based oxides, geometric nonlinear theory of martensitic phase transitions is applied (first used to guide supercompatible shape-memory alloys) to predict the reversibility of the AFE–FE transition by using density-functional theory to assess AFE/FE interfacial lattice-mismatch strain that assures ultralow electric hysteresis and extended fatigue lifetime. A good correlation of mismatch strain with electric hysteresis, hence, with energy efficiency of AFE capacitors is observed. Here, guided by theory, high-throughput material search is conducted and AFE compositions with a near-perfect charge–discharge energy efficiency (98.2%), i.e., near-zero hysteresis are discovered. And the fatigue life of the capacitor reaches 79.5 million charge–discharge cycles, a factor of 80 enhancement over AFE ceramics with large electric hysteresis.
Effect of Al substitution on the magnetocaloric properties of Ni-Co-Mn-Sn multifunctional alloys
Ni-Mn-based metamagnetic shape memory alloys have been extensively studied due to their multifunctional properties including magnetocaloric effect (MCE), magnetoresistance effect, and magnetic-field-induced strain. However, the large intrinsic thermal hysteresis of these alloys sets a limit on further development for practical applications. Here, with a few atomic percent Al substitution, we greatly reduced the thermal hysteresis and transformation temperature interval while maintaining the large magnetization difference (ΔM) between austenite and martensite. The simultaneous achievement of a large ΔM and low thermal hysteresis leads to a remarkable enhancement in MCE upon Al substitution, and a maximum magnetic entropy change of 31.6 J/kg K is achieved under a field of 5 T in Ni 39 Co 11 Mn 40 Sn 8 Al 2 . Furthermore, our in situ high-energy X-ray diffraction experiments show that the Ni 39 Co 11 Mn 40 Sn 8 Al 2 has good geometric compatibility between austenite and martensite, which accounts for the reduced thermal hysteresis in this alloy. This study is instructive for the development of high-performance magnetocaloric materials in the metamagnetic shape memory alloys.
Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles
Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $n_{crit}\sim Q^{1/3}$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thus preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($b_{st}^{2}$) scaling of the critical density as $n_{crit}\sim(1+b_{st}^{2})^{-5/3}$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.
Effective Landau-type model of a Hf x Zr 1 - x O 2 -graphene nanostructure
Here, to describe the charge-polarization coupling in the nanostructure formed by a thin Hf x Zr 1-x O 2 film with a single-layer graphene as a top electrode, we develop the “effective” Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf 1-x Zr x O 2 films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the Hf 1-x Zr x O 2 film polar properties on the film thickness, grain size and surface energy. The proposed Landau free energy with five “effective” expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf 1-x Zr x O 2 film influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 – 25) nm thick Hf 1-x Zr x O 2 films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 ≤ x ≤ 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs, because the working voltages applied to the Hf 1-x Zr x O 2 film (which acts as a gate) can be relatively low (less than 2 V). These low voltages are sufficient to induce the pronounced hysteresis of ferroelectric polarization in the Hf 1-x Zr x O 2 gate, which, due to the strong electric coupling, induces the hysteresis of the graphene charge.
Anomalous polarization reversal in strained thin films of CuInP 2 S 6
Strain-induced transitions of polarization reversal in thin films of a ferrielectric CuInP 2 S 6 (CIPS) with ideally conductive electrodes are explored using the Landau-Ginzburg-Devonshire approach with an eighth-order free energy expansion in polarization powers. Due to multiple potential wells, the height and position of which are temperature and strain dependent, the energy profiles of CIPS can flatten in the vicinity of the nonzero polarization states. Thereby, we reveal an unusually strong effect of a mismatch strain on the out-of-plane polarization reversal, hysteresis loop shape, dielectric susceptibility, and piezoelectric response of CIPS films. By varying the sign of the mismatch strain and its magnitude in a fairly narrow range, quasistatic hysteresisless paraelectric curves can transform into double, triple, and other types of pinched and single hysteresis loops. The strain effect on the polarization reversal is opposite, i.e., anomalous, in comparison with many other ferroelectric films, for which the out-of-plane remanent polarization and coercive field increase strongly for tensile strains and decrease or vanish for compressive strains. For definite values of temperature and mismatch strain, the low-frequency hysteresis loops of polarization may exhibit negative slope in a relatively narrow range of external field amplitude and frequency. The low-frequency susceptibility hysteresis loops, which correspond to the negative slope of polarization loops, contain only positive values, which can be giant in the entire range of field changes. In conclusion, the corresponding piezoresponse also reaches giant values, being maximal near coercive fields.
Important powertrain dynamics for developing models for control of connected and automated electrified vehicles
Connected and Automated Vehicles (CAV) technology presents significant opportunities for energy saving in the transportation sector. CAV technology forecasts vehicle and powertrain power needs under various terrain, ambient, and traffic conditions. Integration of the CAV technology in Hybrid Electric Vehicles (HEVs) provides the opportunity for optimal vehicle operation. Indeed, Hybrid Electric Vehicle powertrains present high degrees of flexibility and possibility for choosing optimum powertrain modes based on the predicted traction power needs. In modeling complex CAV powertrain dynamics, the modeler needs to consider short-time scale powertrain dynamics, such as engine transients, and hysteresis of mode-switching for a multi-mode HEV. Therefore, the powertrain dynamics essential for developing powertrain controllers for a class of connected HEVs is presented. To this end, control-oriented powertrain dynamic models for a test vehicle consisting of full electric, hybrid, and conventional engine operating modes are developed. The resulting powertrain model can forecast vehicle traction torque and energy consumption for the specified prediction horizon of the test vehicle. The model considers different operating modes and associated energy penalty terms for mode switching. Thus, the vehicle controller can determine the optimum powertrain mode, torque, and speed for forecasted vehicle operation via utilizing connectivity data. The powertrain model is validated against the experimental data and shows prediction error of less than 5% for predicting vehicle energy consumption. The model is used to create energy penalty maps that can be used for CAV control, for example fuel penalty map for engine torque changes (10–40 Nm) at each engine speed. The results of model-based optimization show optimum switching delays ranging from 0.4 to 1.4 s to avoid hysteresis in mode switching.
100th Anniversary of the discovery of Ferroelectricity: How did it impact the Current Day Physics
Ferroelectric phenomenon in materials was identified for the first time about 100 years ago in 1920 (paper presented at the American Physical Society (APS) Meeting in Washington, DC, April 1920), at the University of Minnesota, MN, USA; by then a graduate student Joseph Valasek working under the supervision of Prof. W. F. G. Swann. Before this announcement, some prominent physicists had envisaged and even predicted the possibility of an analogous electrical behavior in some materials to that of the well-established hysteresis behavior of magnetic materials. Debye’s work tried to put the concept of Curie temperature in some solids in 1912 but it Schrodinger in the same year who formalized the term ferroelectricity. The names ferroelectrics, ferroelectricity, etc. are fully symbolic of the hysteresis loop-like behavior originally displayed by the magnetic materials. Later on when the physics of ferroelectrics was formalized, most of the terminology, e.g. Curie constant, Curie temperature, etc. that characterized the hysteresis loop or its related characteristics, was inherited from the magnetic materials. Over the years basic ferroelectricity concepts were significantly developed. The observed ferroelectric hysteresis loop between electric polarization vs applied electric field in the then well-known material, Rochelle salt, was presented at the April meeting of the APS, held at, then known as, the National Bureau of Standards or NBS (currently known as the National Institute for Standards and Technology or NIST). Interestingly, even though Valasek submitted the abstract but it was presented by Swann (and Valasek was not present). It happened to be a unique case of the earlier predictions but it did not catch the attention of many researchers working in the field of electronics. Whatever the reason, there was a little attention paid to this newly identified effect in Rochelle salt in 1920s. In late 1920, Valasek’s paper was submitted to Physical Review and that appeared in April 1921. Both ferroelectricity and ferromagnetism are the observational or experimental characteristics of a special class of materials which were not predicted by any basic theoretical approaches to start with at that time.
Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site
Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.
Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site
Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.
Inkjet-printed SnO x as an effective electron transport layer for planar perovskite solar cells and the effect of Cu doping
Inkjet printing is a more sustainable and scalable fabrication method than spin coating for producing perovskite solar cells (PSCs). Although spin-coated SnO 2 has been intensively studied as an effective electron transport layer (ETL) for PSCs, inkjet-printed SnO 2 ETLs have not been widely reported. Here, we fabricated inkjet-printed, solution-processed SnO x ETLs for planar PSCs. A champion efficiency of 17.55% was achieved for the cell using a low-temperature processed SnO x ETL. The low-temperature SnO x exhibited an amorphous structure and outperformed high-temperature crystalline SnO 2 . The improved performance was attributed to enhanced charge extraction and transport and suppressed charge recombination at ETL/perovskite interfaces, which originated from enhanced electrical and optical properties of SnO x , improved perovskite film quality, and well-matched energy level alignment between the SnO x ETL and the perovskite layer. Furthermore, SnO x was doped with Cu. Cu doping increased surface oxygen defects and upshifted energy levels of SnO x , leading to reduced device performance. A tunable hysteresis was observed for PSCs with Cu-doped SnO x ETLs, decreasing at first and turning into inverted hysteresis afterwards with increasing Cu doping level. This tunable hysteresis was related to the interplay between charge/ion accumulation and recombination at ETL/perovskite interfaces in the case of electron extraction barriers.