Search NASA⌕ Search

SEARCH · Search NASA

Results for “load changes”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

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

Flattop-HEU Benchmark Reevaluation Summary

The Flattop critical assembly was first constructed in the 1950’s at Los Alamos National Laboratory as a follow-on to the Topsy experiment. Flattop is composed of a sphere of special nuclear material (SNM) surrounded by a thick spherical reflector made of natural uranium (NU). Two SNM cores currently exist: a highly-enriched-uranium (HEU) core and a plutonium core. The reflector is composed of three parts: a stationary hemisphere and two movable quarter spheres. For fine reactivity control, there are three control rods of NU that are inserted into the stationary hemisphere from underneath the assembly. The final components that allow for reactivity adjustment are the glory hole pieces, mass adjustment buttons, and hemispherical caps. These pieces can be loaded in various configurations to change the available reactivity loaded in the system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Influence of Phase-Change Materials and Operating Periods on Load Shifting for Thermal Storage Heat Pumps in European Climate Zones

As the reliance on heat pumps (HPs) for space conditioning increases worldwide, ex-tensive grid strains during utility peak hours are expected. Significant research efforts focus on investigating the use of integrated thermal energy storage (TES) systems to reduce peak energy demands. Conventionally, TES systems are often integrated to store and provide either cooling or heating, but not both. In this work, a single room-temperature phase-change material (PCM)-TES was integrated for dual-mode opera-tion with an R290 air-to-water 10.5 kW (3-Ton) HP. Two PCM melting temperatures (17°C, 22°C) were assessed to investigate their impact on system performance and cost reductions for Milan, Italy (4A); Barcelona, Spain (3A); and Oslo, Norway (6A). HP and HP-TES performance maps were generated using Modelica-based transient models and co-simulated with a prototype residential building using the Spawn of En-ergyPlusTM with simple rule-based controls. Our simulations indicate that PCM-17°C provided annual cost savings at higher peak-to-off-peak cost ratios, as the demand re-duction dominated the savings potential. The recharge energy costs and requirements became more dominant at lower cost ratios, making PCM-22°C more beneficial at lower cost ratios. Finally, despite Oslo's heating-dominated climate, the utility peak hours favor cooling demand reductions, making PCM-17°C more appealing. In all cases, the cost savings potential stabilized between 4 – 10% regardless of cost ratio, suggesting that system-level component and operation control optimization is required to maximize cost savings further.

25 ENERGY STORAGE↗
DOE OSTI2024

Reviewing Flexibility in Industrial Electrification: U.S. Green Ammonia and Steel Industries [Slides]

The renewable energy transition in the power sector involves a paradigm shift for flexibility. Supply flexibility faces new constraints due to the increased share of variable renewable resources. Increased demand flexibility can allow less use of peaking power plants and delay need for additional capacity and transmission. Industrial customers are larger on average than residential and commercial consumers and have typically provided the largest share of demand response in the United States. We consider industrial demand, studying characteristics of flexible industrial loads. We examine the dynamics of change occurring around industrial load flexibility by focusing on two case studies: green ammonia and steel production via electric arc furnaces. Electric arc furnace steel production is an important component of current demand response programs, whereas green ammonia and green fuels offer new paradigms for flexibility. We analyze the structure and functions of the technological innovation systems of load flexibility in those two industries via interviews with twenty-two stakeholders. We conclude that the technological innovation systems are not well-functioning for flexibility in EAF based steelmaking, but are in the green ammonia space. Explicit connections between scope two greenhouse gas emissions reporting and flexibility are lacking, and industry stakeholders do not appear to make a connection between decarbonization and load flexibility.

25 ENERGY STORAGE↗
DOE OSTI2025

Hydraulic Conductivity Measurements, Utqiagvik (Barrow), Alaska, 2014

Six individual ice cores were collected from the Barrow Environmental Observatory in Barrow, Alaska, in May of 2013 as part of the Next Generation Ecosystem Experiment (NGEE). Each core was drilled at a different location to varying depths. After drilling, the cores were stored in coolers packed with dry ice and flown to Lawrence Berkeley National Laboratory (LBNL) in Berkeley, CA. 3-dimensional images of the cores were constructed using medical X-ray computed tomography (CT) scanner at 120kV. Hydraulic conductivity samples were extracted from these cores at LBNL Richmond Field Station in Richmond, CA, in February 2014 by cutting 5 to 8 inch segments using a chop saw. Samples were packed individually and stored at -20C freezing temperatures to minimize any changes in structure or loss of ice content prior to analysis. Hydraulic conductivity was determined through falling head tests using a permeameter [ELE International, Model #: K-770B] (Appendix A). Samples were placed in a latex membrane via a membrane stretcher while frozen. Use of a membrane stretcher made the membranes easier to secure and minimized contact with the sample. A clear polycarbonate sleeve, fabricated with a stainless steel ring at the bottom to keep the sleeve from floating, was placed around the sample inside the permeameter to minimize deformation during analysis. The permeameter was filled with water and 1.0 PSI of air was applied for confining pressure during sample defrost. Outflow valves were left open to allow for incremental thawing and samples were left to thaw for approximately 12 hours. After approximately 12 hours of thaw, initial falling head tests were performed. When the flow was significantly too fast or too slow, the analysis was stopped and the burette size was adjusted accordingly (i.e. a larger diameter burette was used for flows that were faster than desired or a smaller diameter burette was used for flows that were slower than desired). Two to four measurements were collected on each sample and collection stopped when the applied head load exceeded 25% change from the original load. Analyses were performed between 2 to 3 times for each sample. The final hydraulic conductivity calculations were computed using methodology of Das et al., 1985.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Projecting Electric Vehicle Electricity Demands and Charging Loads

After over a century of petroleum dominance, electric vehicles are rapidly disrupting the transportation energy landscape. At the same time, the electric power systems are undergoing profound changes as variable renewables replace dispatchable fossil generators. It is critically important to understand how transportation electrification will impact electricity demand, including changes in the load shapes that characterize the system, and the value of flexible electric vehicle charging to better balance electricity demand and supply. This talk focuses on methods to project electricity demand for EV charging with high spatiotemporal fidelity to enable electricity system modeling and analysis and summarizes results from recent NREL studies in this area.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2024

A Framework for Identifying Building Energy Models of Localized Utility Service Areas Using Smart Meter Data

Bottom-up load modeling of buildings offers a versatile approach to simulating baseline demand and scenarios of future technology evolution and adoption at the individual building level. This capability is essential to understanding how future load shapes may change with the adoption of electric equipment and vehicles, particularly as it relates to grid planning and infrastructure investments. Traditionally, grid planning techniques have used historical load data to predict future load and infrastructure needs. However, with the anticipated rise in adoption of electrification technologies such as heat pumps and electric vehicles, historical data become less reliable predictors of the future. By employing ResStock, a high-fidelity building stock modeling tool, we can fine-tune electrification scenarios and aggregate models to represent varying geographic resolutions of the grid system, while considering the underlying features of homes. This may enable a more accurate and responsive approach to anticipate and plan for the evolving landscape of energy demands. We present a new framework that leverages building stock energy modeling to identify building models that align with the load shapes and housing attributes of buildings with AMI data. This approach applies two model layers: (1) a classification step that identifies the presence of air conditioning, electric heating, and electric water heating, and (2) an optimization routine that identifies building energy models aligning with load profile data from advanced metering infrastructure meters. This report demonstrates one approach to deploying this framework, and presents results for three test cases that use both modeled and AMI data to assess performance. For a test case using AMI data in Fort Collins, Colorado, we observed a median monthly electricity load CV-RMSE of 16.6%, and a top ten daily heating and cooling median absolute percent error of 7.7% and 8.3%, respectively. For each AMI meter, we identify a set of potential energy models so that downstream use-cases can account for uncertainty driven by variability of baseline technologies and occupant behavior, which impact the response to electrification and energy efficiency scenarios. Our results indicate that ResStock has potential as a scalable solution for modeling residential energy demand at local grid resolutions. Its performance depends on location-specific factors, underlying building characteristics, and the level of aggregation, offering a path towards more precise and adaptive distribution grid planning for the evolving energy landscape.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Radioisotope Power System Cask Trailer Chiller Low Heat Load Study

This report documents the study of the Radioisotope Thermoelectric Generator Transport System (RTGTS) chillers capability to transport a Radioisotope Power System (RPS) with low heat loads in the cask. Tests of the RTGTS chillers included setting the glycol/water (coolant) temperature to a control set point of 40° and 60°F while applying heat loads ranging from 2500 to 0 W heat load. The tests performed showed that the chiller systems could handle the entire range of heat loads. However, a noise/vibration started to occur when the set point was 60°F with low heat loads. Three recommended changes to the RTGTS chillers are recommended with varying implementation windows. The first is an operational change that can be implemented immediately. The RTGTS chillers should have a set point of 40°F when transporting an RPS less than 1500 W. The second is a short-term change (within the next two years) augmenting RPS heat by utilizing one RTGTS onboard heater (2500 W) while transporting an RPS with less than 1500 W which will also require a change to the Instrument and Data Acquisition System (IDAS) to allow monitoring and possible control of power to the onboard heaters. However, if the heater does fail during transport, no emergency actions or immediate repairs will be required due to the chillers’ ability to handle low heat loads. Running a 2500 W heater when transporting low heat loads aids in prolonging the life of the chillers. The remaining recommended long-term change is to replace the current chillers with a new portable or permanently installed chiller system in the next five to ten years. This time frame takes advantage of the proven reliability and durability of the current chillers while the new chillers are phased in and proven to the RPS program.

42 ENGINEERING↗
DOE OSTI2026

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

On the effect of strain rate during the cyclic compressive loading of liquid crystal elastomers and their 3D printed lattices

Nematic liquid crystal elastomers (LCEs) are a unique class of network polymers with the potential for enhanced mechanical energy absorption and dissipation capacity over conventional network polymers because they exhibit both conventional viscoelastic behavior and soft-elastic behavior (nematic director changes under shear loading). This additional inelastic mechanism makes them appealing as candidate damping materials in a variety of applications from vibration to impact. The lattice structures made from the LCEs provide further mechanical energy absorption and dissipation capacity associated with packing out the porosity under compressive loading. Understanding the extent of mechanical energy absorption, which is the work per unit mass (or volume) absorbed during loading, versus dissipation, which is the work per unit mass (or volume) dissipated during a loading cycle, requires measurement of both loading and unloading response. Here, in this study, a bench-top linear actuator was employed to characterize the loading-unloading compressive response of polydomain and monodomain LCE polymers and polydomain LCE lattice structures with two different porosities (nominally, 62% and 85%) at both low and intermediate strain rates at room temperature. As a reference material, a bisphenol-A (BPA) polymer with a similar glass transition temperature (9 °C) as the nematic LCE (4 °C) was also characterized at the same conditions for comparing to the LCE polymers. Based on the loading-unloading stress-strain curves, the energy absorption and dissipation for each material at different strain rates (0.001, 0.1, 1, 10 and 90 s -1 ) were calculated with considerations of maximum stress and material mass/density. The strain-rate effect on the mechanical response and energy absorption and dissipation behaviors was determined. The energy dissipation ratio was also calculated from the resultant loading and unloading stress-strain curves. All five materials showed significant but different strain rate effects on energy dissipation ratio. The solid LCE and BPA materials showed greater energy dissipation capabilities at both low (0.001 s -1 ) and high (above 1 s -1 ) strain rates, but not at the strain rates in between. The polydomain LCE lattice structure showed superior energy dissipation performance compared with the solid polymers especially at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Single-Atom Nickel Dispersed on Mo 2 TiC 2 T x MXene for Upcycling of Polyolefins via Catalytic Hydropyrolysis for Fuels and Lubricants

Plastic waste has been a major social and environmental problem. Chemical upcycling of plastic waste into transportation fuels or lubricant oil has emerged as an attractive and promising approach to convert plastic waste into valuable products. However, conventional catalyst systems often require high-pressure hydrogen, prolonged residence time, and the use of noble metal catalysts. Here, we present zero-valent nickel single atoms on Mo-based MXenes (Ni/Mo 2 TiC 2 T x ) for effective hydrogenation of polyolefins. In situ spectroscopic and microscopic characterizations demonstrate the formation of Ni–Mo and Ni–C bonds at a low Ni loading, resulting in the formation of dispersed intercalated Ni atoms. The catalysts were then employed to convert polyethylene into fuels and lubricant molecules under atmospheric hydrogen pressure, short vapor residence time (τ < 1 s), and mild reaction temperature (300 °C). We further demonstrated the possibility of tuning the Ni structure by changing the Ni loading, attributed to the strong metal–support interactions (MSIs) with the Mo 2 TiC 2 T x support. Under optimized conditions, 22.5 C% gasoline, 12.0 C% diesel, 19.2 C% jet fuel (JP-8), and 37.4 C% lubricant were produced over 0.5% Ni/Mo 2 TiC 2 T x . This work highlights the potential of utilizing MSIs of Mo 2 TiC 2 T x MXenes to synthesize single-atom catalysts (SACs) for upcycling plastic waste into higher-value products.

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

Fuel Fabrication Specification Impact Analysis for NBSR LEU Conversion

As part of a national initiative to enhance nuclear security and reduce proliferation risks, significant efforts have been undertaken by the National Nuclear Security Administration Material Management and Minimization Office of Reactor Conversion Program to convert U.S. high performance research reactors (USHPRRs) from the use of highly enriched uranium (HEU) to low-enriched uranium (LEU), including the National Bureau of Standards Reactor (NBSR). The current plan is to procure LEU fuel assemblies from commercial fabricators according to fuel specifications tailored for each USHPRR. The analysis conducted at Brookhaven National Laboratory was part of an effort to identify the sources of uncertainty in the fuel specifications that may impact the performance of the NBSR core after its conversion and, in particular, to assess the range of acceptable tolerance limits from the perspective of core safety and reactor performance. Using the stochastic neutronics code MCNP 6.2, the variations in important NBSR neutronics characteristics were analyzed as a function of the specification parameters independently and in combination. The important NBSR specification parameters analyzed were the fuel isotopic composition, the amount of impurity content in cladding, the fuel plate thickness, and the fuel element 235U mass loading. The range of variation of each specification parameter was based on the technical specification limit or available as-fabricated assay data and uncertainties. The NBSR neutronics characteristics selected for analysis were the reactor reactivity characteristics at equilibrium and the equilibrium fuel cycle length. Results show that with variations in the fabrication parameters of the as-fabricated U-10Mo fuel within the specification limitations, the excess reactivity of the NBSR LEU core remains well below the 15% Δk/k technical specification limit, and the shutdown margin is always significantly greater than the required 0.68% Δk/k. This ensures that the NBSR can be operated safely and reliably shut down for all analyzed cases within the specified fabrication limits after the LEU conversion. In the prototypic case, the fuel cycle length was 1.5 days longer than the targeted 38.5 days. In a credible worst-case scenario, where all low-reactivity parameters were combined, the fuel cycle length was reduced to 35.5 days, which is still considered manageable for reactor operations. Variations in cycle length are primarily driven by changes in 235U loading, with other parameters having secondary effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dispatching Grid-Forming Inverters in Grid-Connected and Islanded Mode

This paper explores the dispatchability of grid-forming (GFM) inverters in grid-connected and islanded mode. GFM inverters usually use droop control to automatically share power with other GFM sources (inverters and synchronous generators) and follow the change in the load demand; however, they can be dispatched like their grid-following (GFL) counterparts to output the target active and reactive power. This will help grid operators better manage their inverter-based resources (IBRs) to improve operation efficiency and reliability; therefore, this paper proposes an innovative concept of dispatching GFM sources (inverters and synchronous generators) to output the target power in both grid-connected and islanded mode by adjusting their droop intercepts. The fundamental principle is that the GFM inverter's active and reactive power is dictated by its frequency and voltage, and thus dispatching the active and reactve power of a GFM inverter can be achieved through dispatching its frequency and voltage. Moreover, the concept distinguishes the dispatch rules for grid-connected and islanded mode. Finally, the concept is validated with an example microgrid system with two GFM inverters, one diesel generator, one GFL inverter, and the load in both grid-connected and islanded mode. This pioneering work results in practical guidance for the development of energy management systems for future electric grids with GFM and GFL inverters.

dispatching↗
DOE OSTI2024

Energy Storage Improves Power Plant Flexibility and Economic Performance

Most existing coal-fired power plants were designed for sustained operation at full load to maximize efficiency, reliability, and revenue, as well as to operate air pollution control devices at design conditions. Depending on plant type and design, these plants can adjust output within a fixed range in response to plant operating or market conditions. The need for flexibility driven by increased penetration of variable and non-dispatchable power generation, such as wind and solar, is shifting the traditional mission profile of thermoelectric power plants in three ways: more frequent shutdowns when market or grid conditions warrant, more aggressive load ramp rates (rate of output change), and a lower minimum sustainable load, which provides a wider operating range and helps avoid costly plant shutdowns. Recent studies have shown that the flexibility of a coal-fired power plant can be improved by energy storage. The objective of this work was to analyze a set of energy storage options and determine their impact on the flexibility and economics of a representative coal-fired power plant. The effect of three energy storage systems integrated with a coal power plant on plant flexibility and economics was investigated. The results obtained in this project show that energy storage systems integrated with a thermal power plant improve plant flexibility and participation in the energy and ancillary services markets, which improves plant financial performance.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Dispatching Grid-Forming Inverters in Grid-Connected and Islanded Mode: Preprint

This paper explores the dispatch-ability of grid-forming (GFM) inverters in grid-connected and islanded mode. Grid-forming (GFM) inverters usually use droop control to automatically share power with other GFM sources (inverters and synchronous generator (SG)) and follow the change of the load demand. However, they can be dispatched like their grid-following (GFL) counterparts to output the target active and reactive power. This will help the grid operator better manage their inverter-based resources (IBRs) for improved operation efficiency and reliability. Therefore, this paper proposes an innovative concept of dispatching GFM sources (inverters and SG) to output the target power for both grid-connected and islanded mode by adjusting their droop intercepts. The fundamental principle of doing so is that the GFM inverter's active and reactive power is dictated by its frequency and voltage, thus, dispatching active and reactve power of a GFM inverter can be achieved through dispatching its frequency and voltage. Moreover, the concept distinguishes the dispatch rules for grid-connected and islanded mode. Finally, the concept is validated with an example microgrid system with two GFM inverters, one diesel generator, one grid-following inverter and loads in both grid-connected and islanded mode. This pioneering work results in a practical guidance for power system energy management system (EMS) development to control a future grid with GFM and GFL inverters.

droop control↗
DOE OSTI2025

Wood template-supported phase change material composites for durable and form-stable thermal energy storage in buildings

Here, to reduce and shift peak energy loads in buildings, phase change materials (PCMs) with high transition enthalpies and transition temperatures near human thermal comfort are desirable for thermal energy storage (TES). Traditional solid/liquid PCMs suffer from leakage during thermal cycling, requiring encapsulation that lowers heat storage capacity. Wood templates (WTs), with porous and hierarchical structures, provide natural encapsulation scaffolds for PCM containment. Polyethylene glycol (PEG) is a compatible PCM with WTs, but when infiltrated alone, ~30 wt% of PEG binds to wood with no detectable phase change, limiting TES efficiency. To address this, we developed a method to (1) improve the form stability of balsa and pine-based composites (BWT + PCMs and PWT + PCMs) and (2) reduce inactive PCM within the composites to ~10 wt%. BWT + PCMs exhibit transition properties upwards of 114.2 J/g at 25.4 degrees C, with no degradation after 1000 thermal cycles, and similar stiffness compared to raw balsa. Meanwhile, PWT + PCMs exhibit 25 % higher storage efficiency with the addition of poly(ethylene glycol) diacrylate compared to solely PEG-infiltrated PWT. PWT + PCM retains 84 % of raw pine's mechanical stiffness, sufficient for light-duty construction. Our shape-stabilized WT + PCM composites enrich the functionality of wood materials as both ideal TES material candidates and light-duty building construction applications.

25 ENERGY STORAGE↗
DOE OSTI2024

Infinitely rugged intra-cage potential energy landscape in metallic glasses caused by many-body interaction

The absence of translational symmetry in glassy materials poses a significant challenge in establishing effective structure-property relationships in real space. Consequently, the potential energy landscape (PEL) in phase space is widely utilized to comprehend the complex phenomena in glasses. The classical PEL features a two-scale profile comprising mega-basins and sub-basins, corresponding to α-relaxations (e.g. glass transition) and β-relaxations (e.g. local cage-breaking atomic rearrangements), respectively. Recent studies, however, reveal that sub-basins are not smooth and contain finer structures, the origins of which remain elusive. Here we probe the smoothness of sub-basin bottoms in glasses' PEL by introducing small intra-cage cyclic loading and then measuring the net changes in atomic-level stresses. Compared to glasses with pair interaction, glasses with many-body interaction exhibit orders-of-magnitude larger and loading-dependent stress changes even before the first cage-breaking event takes place, which reflect much more feature-rich sub-basins. We further demonstrate this stark contrast stems from the spatial distribution of individual atom's constraining force field. Specifically, at vanishing perturbations, many-body interactions disrupt the positive-definite synchrony in energy variations of the perturbed atom and the whole system, causing inherently less confined atomic responses and infinitely rugged sub-basins. The implications of these findings for the selective addition or removal of fine structures in the PEL and the subsequent tuning of glassy materials' responses to external stimuli are also explored.

36 MATERIALS SCIENCE↗