Search NASA⌕ Search

SEARCH · Search NASA

Results for “periodic potential”

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 127 records · Page 7

Continued Development and Advanced Testing of DPC Filler Cements (on FY22 R&D and Demonstration Activities) (Progress Report)

Commercial generation of energy by nuclear power plants in the United States (U.S.) has produced thousands of metric tons of spent nuclear fuel (SNF), the disposal of which is the responsibility of the U.S. Department of Energy (DOE). Utilities typically utilize the practice of storing this SNF in dual-purpose canisters (DPCs). DPCs were designed, licensed, and loaded to meet Nuclear Regulatory Commission (NRC) requirements that preclude the possibility of a criticality event during SNF storage and transport, but were not designed or loaded to preclude the possibility of a criticality event during the regulated post-closure period following disposal, which could be up to 1,000,000 years (Price, 2019). There are several options being investigated that could facilitate the disposal of SNF stored in DPCs in a geologic repository (Hardin et al., 2015; SNL 2020b; SNL 2021b). These include: (1) repackage the SNF into canisters that are designed to prevent criticality during the regulated post-closure period following disposal, but with an increased disposal cost estimated at approximately $\$$20B in United States dollars (USD) (Freeze et al., 2019); (2) analysis of the probability and consequences of criticality from the direct disposal of DPCs during a 1,000,000-year post-closure period in several geologic disposal media (Price, 2019); and (3) filling the void space of a DPC with a material before its disposal that significantly limits the potential for criticality over the post-closure regulatory period. This report further investigates the third option, filling DPC already containing SNF with a material to limit the potential for criticality over the post-closure regulatory period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Impurity screening and Friedel oscillations in Floquet-driven two-dimensional metals

We develop a theory for the non-equilibrium screening of a charged impurity in a two-dimensional electron system under a strong time-periodic drive. Our analysis of the time-averaged polarization function and dielectric function reveals that Floquet driving modifies the screened impurity potential in two main regimes. In the weak drive regime, the time-averaged screened potential exhibits unconventional Friedel oscillations with multiple spatial periods contributed by a principal period modulated by higher-order periods, which are due to the emergence of additional Kohn anomalies in the polarization function. In the strong drive regime, the time-averaged impurity potential becomes almost unscreened and does not exhibit Friedel oscillations. Furthermore, this tunable Friedel oscillations is a result of the dynamic gating effect of the time-dependent driving field on the two-dimensional electron system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Occupant-driven end use load models for demand response and flexibility service participation of residential grid-interactive buildings

As demand response becomes increasingly used as a tool to support improved grid flexibility, it is important to consider that there are many potential types of energy end uses that may be used to support such flexibility. Residential appliances, often accounting for 30 % or more of residential energy use, are a currently untapped source of demand flexibility, particularly when aggregated together across homes. To date there has been very limited analysis of residential appliances for use as grid-interactive loads. As such, this research uses disaggregated energy end use data for 564 households, to model the electricity demand flexibility potential of the use of residential dishwashers, clothes washers, clothes dryers, ovens, and ranges (oven + stovetop) on both weekdays and weekends. This includes both at the building level, as well as aggregated to the grid level, specifically the Midcontinent Independent System Operator (MISO) region. This study was divided into two parts. Part 1 focuses on determining appliance-level loads, and Part 2, which involves aggregation to the grid. Findings suggest that among the studied appliances, clothes dryers provide the greatest demand reduction potential for most times of the day, followed by dishwashers and clothes washers. The maximum potential reduction for clothes dryers is found to be approximately at 11:00 a.m. and this potential sustains throughout most of the daytime period. When considering the willingness of households to participate, based on a survey of households in the Midwest region, clothes dryers still have the most potential for demand reduction. The availability of appliances for load modulation on weekdays and weekends indicates similar load reduction potential for all appliances. Overall, the results of this study suggest that there is an opportunity for shifting appliance usage to optimize grid efficiency and enhance demand response strategies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Spatio-Temporal Assessment of Heavy-Duty Truck Incident and Inspection Data

Vehicular incidents, especially those involving tractor trailers, are increasing in number every year. These events are extremely costly for fleets, in terms of damage or loss of property, loss of efficiency, and certainly in terms of loss of life. Although the U.S. Department of Transportation (DOT) is responsible for performing inspections, and fleet managers are encouraged to maintain their fleet and participate in regular inspections, it is uncertain whether these inspections are occurring at a frequency that is necessary to prevent incidents. The Federal Motor Carrier Safety Administration (FMCSA) of the DOT manages and maintains the Motor Carrier Management Information System (MCMIS) dataset, which contains all incident and inspection data regarding commercial vehicles in the U.S. The purpose of this preliminary analysis was to explore the MCMIS dataset through spatiotemporal analyses, to uncover findings that may hint at potential improvements in the DOT inspection process and highlight location-specific trends in the dataset. These analyses are novel, as previous research using the MCMIS dataset only examined the data at the state or county level, not at a national scale. The results from the analyses pinpointed specific major metropolitan areas, namely Harris County (Houston), Texas, and three of the New York boroughs (Kings, Queens, and the Bronx), which were found to have increasing incident rates during the study period (2016–2020). An overview of potential causal factors contributing to this increase are provided as well as an overview of the inspection process, and suggestions for improvement relative to the highlighted locations in Texas and New York are also provided. Ultimately, it is suggested that the incorporation of advanced technology and automation may prove beneficial in reducing the occurrence of events that lead to incidents and may also help in the inspection process.

99 GENERAL AND MISCELLANEOUS↗

Thermal and hydraulic performance of volumetrically heated triply periodic minimal surface heaters

Advanced heat transfer surfaces, mainly Triply Periodic Minimal Surfaces (TPMSs) have great potential to increase heat transfer performance over traditional heat transfer surfaces due to their tortuous flow path and higher surface area-volume ratio. The emergence of additive manufacturing of nuclear fuel motivates experimental investigation of such advanced geometries in the context of a nuclear reactor core. Here, this study employed a novel method to explore the performance of gyroid and diamond TPMS geometries. Volumetrically heated heaters were additively manufactured from conductive polymer filament, and airflow was used to simulate the convective heat transfer from the core within a nuclear reactor. The mass flow and volumetric generation rates were varied, local Nusselt number correlations were obtained using embedded thermocouples to measure the solid and fluid temperature, and friction factor correlations were developed using differential pressure measurements. The gyroid TPMS presented a friction factor 1.5 times higher than the diamond TPMS and 90 times that of the rod bundle. TPMS geometries showed Nusselt numbers 8–10 times higher than that of the rod bundle, also maintaining colder internal temperatures, which suggested that they could sustain 250–275 % higher power density than the rod bundle for a given centerline temperature. Hence, having the potential to greatly increase the safety margins and power output of a nuclear reactor core with similar volume.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulations of Interacting Binary Systems—Pathways to Radio-bright Gamma-Ray Burst Progenitors

Although the association of gamma-ray bursts (GRBs) with massive stellar death is on firm footing, the nature of the progenitor system and the key ingredients required for a massive star to produce a GRB remain open questions. Here, we investigate the evolution of a 15–25M ⊙ massive star with a 10–15M ⊙ black hole (BH) using Modules for Experiments in Stellar Astrophysics. We quantify companion-influenced angular-momentum evolution over a stellar lifetime for orbital periods where tides are significant, varying stellar and BH masses, initial stellar spin, and accretion and dynamo prescriptions while tracking mass loss and angular momentum. The final spin is set by tidal torques versus stellar winds. For binaries that initially avoid Roche-lobe overflow, tides can spin up the star, but a late-stage expansion can drive tidal stripping; associated mass and angular-momentum loss can suppress spin-up. We find that massive-star BH binaries at comparable mass ratios may be potential GRB progenitors for short orbital periods (∼20–5 × 10 2 days) and long orbital periods (∼2 × 10 3 –4 × 10 3 days), where our suite of lifetime simulations reveals a favored parameter space with negligible mass loss and enough spin angular momentum to power a GRB jet. For initially nonrotating stars, this provides a lower limit on the final spin above a threshold estimate consistent with forming a postcollapse BH mass of 5–10M ⊙ with a spin parameter ≥ 0.5. For initially rapidly rotating stars, tidal interactions may sustain high spin when mass loss is negligible because the binary is not tidally synchronized.

Hernandez, Angel [Univ. of Colorado, Boulder, CO (↗

Staircase formation by resonant and non-resonant transport of potential vorticity

The E × B staircase is a quasi-periodic pattern of pressure profile corrugations. In this work, we present a new mechanism for E × B staircase formation that involves resonant transport versus non-resonant transport. We start from a potential vorticity evolution system and use quasi-linear theory, a model dispersion relation, and a bi-Lorentzian spectrum approximation, to construct the relation between the fluxes and the profiles. With these fluxes, we close the profile evolution equations and the extended turbulence intensity evolution equation, which together constitute a turbulence-profile evolution system. In this system, the Doppler effect from the E × B mean flow can cause resonance between trapped ion precession motion and the trapped ion mode, which drives a resonant transport contribution to the fluxes. The profiles will be flattened where the resonant transport is switched on. In contrast, for the regions of non-resonant transport, profiles are steeper. A quasi-periodic pattern of profile corrugation (the E × B staircase) spontaneously emerges in this system, which is the two states mentioned above, arranged as alternating layers in space. The feedback processes during the staircase pattern formation are identified. An estimate of the critical value of the boundary heat flux is obtained, above which the staircase formation will be triggered. An estimate scaling of the step size in the staircase pattern is obtained. The resonant turbulent transport is also a mechanism for collisionless saturation of zonal flow. Furthermore, this work is related to internal transport barrier formation and suggests some new scenarios, such as an enhanced confined L mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗

Charged-Particle Bound States in Periodic Boxes

We consider the binding energy of a two-body system with a repulsive Coulomb interaction in a finite periodic volume. We define the finite-volume Coulomb potential as the usual Coulomb potential, except that the distance is defined the shortest separation between the two bodies in the periodic volume. We investigate this problem in one and three-dimensional periodic boxes and derive the asymptotic behavior of the volume dependence for bound states with zero angular momentum in terms of Whittaker functions. We benchmark our results against numerical calculations and show how the method can be used to extract asymptotic normalization coefficients for charged-particle bound states. Furthermore, the results we derive here have immediate applications for calculations of atomic nuclei in finite periodic volumes for the case where the leading finite-volume correction is associated with two charged clusters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A simple, validated approach for design of two-dimensional periodic particle patterns via acoustophoresis

Two-dimensional patterning of microparticles enables a wide range of functional materials, including patterned energy storage electrodes, flexible electronics, and sensor arrays. Particle patterning via acoustics offers an attractive path to generate a wide variety of 2D periodic patterns that introduce tailorable hierarchical porosity, useful for controlling surface area, transport distances, and other properties. This method is most effective with micron scale particles and patterns of tens to hundreds of microns. To enable systematic exploration of the broad design space for such patterns, this work develops a model of 2D and 3D assembly of particles at high loadings and validates the obtained patterns against both experiments and more computationally intensive modeling techniques. Using this simple model, connections are mapped between input parameters (like actuation conditions, particle volume fraction, material properties) and output geometrical features (like void size and shape, pattern connectivity, and surface area) so that they can be tailored to given applications. The utility of this simple model is illustrated by predicting and then experimentally demonstrating new hierarchical patterns resulting from multiple waves of different frequencies interacting. These multiscale patterns offer the potential to lift the limits on surface area, diffusion distances, and other features.

42 ENGINEERING↗

DOE Cold Climate Heat Pump Challenge: Development, Metrics, and Early Field Observations

Space heating in residential buildings is a major contributor of Greenhouse Gas (GHG) emissions in the United States. New, advanced electric heat pumps are poised to provide a low carbon alternative to traditional fossil-based heating, especially in colder climates. Widespread deployment of cold climate heat pumps could help address the significant portion of building emissions and primary energy used in American households, but these gains will require broader acceptance from consumers and decision makers. As part of the “Energy, Emissions, and Equity” (E3) initiative, the U.S. Department of Energy launched the Cold Climate Heat Pump Challenge (CCHP) in 2021 to accelerate the deployment of the next generation of air source heat pumps. The Challenge is currently focused on residential, centrally ducted, electric heat pumps, with a nominal cooling capacity greater than or equal to 24,000 Btu/h (2 tons) and less than or equal to 65,000 Btu/h (5 tons). The Challenge specifications represent a best-in-class heat pump product that provides high-efficiency heating performance in cold climates, employs environmentally friendly low-Global Warming Potential (GWP) refrigerants, and is designed to be grid interactive. Spearheaded by DOE in partnership with the US Environmental Protection Agency (EPA) and Natural Resources Canada (NRCan), the CCHP Challenge brings together numerous major heat pump manufacturers and key stakeholders including utilities and state agencies across the country. In addition to performance testing in a laboratory environment, a key component of this research is evaluating the in-field performance of the prototype heat pumps developed as part of the Challenge. This paper will discuss the development of the CCHP Challenge, key performance specifications, and the energy and non-energy metrics that will be evaluated through the field study. The types of data that are being collected from the units installed in homes across North America as part of the field validation will be described, including details about the sensors and data acquisition system. Data cleaning and analysis methodologies will be discussed, along with early observations from the winter 2022-2023, spring and summer performance periods, and an assessment of non-energy metrics through pre- and post-installation homeowner surveys. Challenges uncovered and lessons learned throughout the process will also be discussed. Finally, the paper will discuss next steps for the second winter assessment period, areas needing additional research, and potential applications of the data collected and analyzed through this work.

Mendon, Vrushali V.↗

Empirical Assessment of Interregional Coordination to Support Resource Adequacy [Slides]

This study examines where interregional transmission could most effectively support resource adequacy in the contiguous United States. We use hourly load, renewable generation, and real-time price data from 2016–2023 for 18 planning subregions to identify periods of elevated adequacy risk, defined as the top 100 annual hours of net load and wholesale prices in each region. We then measure the temporal coincidence of these peak periods between adjacent regions and compare price patterns to assess the potential for capacity sharing. Results show that NorthernGrid West, a winter-peaking region, has low coincidence of peak net load with its summer-peaking neighbors, indicating high potential for interregional support. In contrast, regions in the Northeast have highly coincident peak periods, suggesting limited adequacy benefits from additional transmission. Price-based analysis shows peak-hour differences in the Midwest and between ERCOT and neighboring regions, indicating potential economic benefits from increased transfers. The findings provide an empirical screening of where transmission may offer the greatest reliability benefits without adding new generation capacity.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impacts of Climate Change on the Generation Potential of Solar and Wind Energy Systems in India

Low-carbon energy sources like wind and solar are essential for decarbonizing the electricity sector. In addition, the cost of electricity generated from these sources has plummeted over the last decade. Therefore, these energy sources are poised to take a significant share of the total installed capacity soon. However, they are susceptible to the impacts of climate change as their generation potential depends on the weather conditions. Estimating the installed capacity requirements of solar and wind energy to decarbonize the power sector without accounting for these possible changes in generation potential could lead to missing out on the set climate goals and meeting future electricity demand. This study evaluates the effect of climate change on the generation potential of wind and solar energy systems in India for two future periods, 2050 and 2070, under two climate scenarios or Shared Socioeconomic Pathway (SSP): SSP245 and SSP585. Almost all regions show a decrease, and most regions show a significant decline (>5%) in the generation potential of solar Photovoltaic (PV) as compared to 2010 levels under both climate scenarios and future periods. The changes in the generation potential of wind energy are more significant (>10%), and the majority of regions show a decline in generation potential. Southwestern and central regions show an increase in wind generation potential for 2070 as compared to 2050 levels under the SSP245 scenario and the SSP585 scenario, respectively.

climate change↗

Importance of Depth and Artificial Structure as Predictors of Female Red Snapper Reproductive Parameters

Abstract The Red Snapper Lutjanus campechanus is a structure‐associated species occurring across a wide depth range in the northern Gulf of Mexico. We used the random forest machine learning algorithm to understand which habitat and individual fish characteristics could predict reproductive parameters of female Red Snapper. We evaluated fish captured from 2016 to 2018 on three artificial structure types with various structure heights at depths of 100 m or less. Overall, we found that depth and month were important predictors for most reproductive parameters, but the type of structure (artificial reefs, oil platforms, and rigs‐to‐reefs structures) was not important. Maturity was correctly classified in 88.9% of the cases when using the random forest ensemble model, with important predictors including FL, depth, structure height, and month of collection. Spawning seasonality (measured as gonadosomatic index [GSI]) was correctly classified in 59.5% of the cases when using histology reproductive phase, FL, month, and depth variables. Reproductively active or inactive females were correctly classified in 89.3% of the cases using GSI, month, FL, and depth, while females in the developing versus spawning capable phases were correctly classified in 82.2% of the cases using GSI, FL, month, and depth. Histological indicators that show potential spawning within a 36‐h period were correctly classified 61.5% of the time, with the best predictors being depth, FL, GSI, and month. Stepwise regression indicated that month was the only factor that significantly predicted contrasts in relative batch fecundity, with significantly greater values in August compared to all other months. Our findings suggest that female Red Snapper reproductive effort is not consistently or well predicted by artificial structure type or height but that a combination of fish FL, month, and depth can predict reproductive characteristics of female Red Snapper.

Brown‐Peterson, Nancy J.↗

Key predictors of soil organic matter vulnerability to mineralization differ with depth at a continental scale

Abstract Soil organic matter (SOM) is the largest terrestrial pool of organic carbon, and potential carbon-climate feedbacks involving SOM decomposition could exacerbate anthropogenic climate change. However, our understanding of the controls on SOM mineralization is still incomplete, and as such, our ability to predict carbon-climate feedbacks is limited. To improve our understanding of controls on SOM decomposition, A and upper B horizon soil samples from 26 National Ecological Observatory Network (NEON) sites spanning the conterminous U.S. were incubated for 52 weeks under conditions representing site-specific mean summer temperature and sample-specific field capacity (−33 kPa) water potential. Cumulative carbon dioxide respired was periodically measured and normalized by soil organic C content to calculate cumulative specific respiration (CSR), a metric of SOM vulnerability to mineralization. The Boruta algorithm, a feature selection algorithm, was used to select important predictors of CSR from 159 variables. A diverse suite of predictors was selected (12 for A horizons, 7 for B horizons) with predictors falling into three categories corresponding to SOM chemistry, reactive Fe and Al phases, and site moisture availability. The relationship between SOM chemistry predictors and CSR was complex, while sites that had greater concentrations of reactive Fe and Al phases or were wetter had lower CSR. Only three predictors were selected for both horizon types, suggesting dominant controls on SOM decomposition differ by horizon. Our findings contribute to the emerging consensus that a broad array of controls regulates SOM decomposition at large scales and highlight the need to consider changing controls with depth.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of modeled time horizon on quantifying the need for long-duration storage

Long-Duration Energy Storage (LDES) has gained interest due to its key role in attaining a decarbonized, low-cost, and stable grid driven by variable renewable electricity (VRE). Currently, there is a wide range of LDES technologies being developed to provide electricity with 8+ hours of consecutive discharge. However, current capacity expansion models used in long-term planning processes rarely consider low cost LDES as a candidate technology. If they do, the storage balancing horizon (SBH) of the model usually only considers non-consecutive 1-day periods that do not capture the potential of LDES to shift energy across multiple days or even seasons. Addressing these limitations in existing models, this work explores the ways in which the optimal energy storage changes when increasing the number of consecutive days in the SBH and how these changes will impact planners who are determining the future roles of energy storage. Our analysis uses SWITCH, an open-source capacity expansion model with a high spatial resolution for the entire Western Electricity Coordinating Council (WECC) in a zero-carbon scenario in 2050. We find that the number of consecutive days in the SBH changes both the total selected power and energy capacity of LDES when storage energy and power capacity overnight costs are $\$$13 USD/kWh (or less) and $113 USD/kW, respectively. We also find that the amount of required energy in storage to drive a future VRE-driven WECC grid ranges from 2.5 TWh to 16.0 TWh depending on the length of the SBH. The optimal storage duration (energy to power ratio) we obtain ranges from 10 h to 620 h among all the scenarios. Furthermore, depending on the storage cost assumption, we observe different charge/discharge patterns when varying the length of the SBH. Given our results, we anticipate that as more LDES technologies become commercially available, it will be critical to increase the length of the SBH to fully capture the benefits of LDES assets in long-term planning processes of high VRE-driven grids.

42 ENGINEERING↗

Energy savings and life cycle cost analysis of advanced double skin facade system applied to old apartments in South Korea

South Korea is one of the most densely populated countries in the world, and the population density in urban areas is much higher among others. Apartments are the most common form of residential buildings due to their high population density. More than 60% of domestic residential buildings are apartments, with more than 10 million households. The high population density is a problem that has long plagued South Korea because more than 45% of apartments are old apartments that are more than 20 years old. Old apartment buildings have low thermal performance that results in lots of energy consumption. Balconies, which worked as a thermal buffer, are also being removed by residents to expand indoor space. According to the green remodeling project supported by the Korean government, the typical retrofitting method is replacing external windows with the high-efficiency window in old apartments in South Korea. Here, this paper proposes a retrofitting method using an extended double-skin façade (DSF) system that replaces exterior windows and acts as a thermal buffer. The simulation model is developed with the EnergyPlus simulation program to conduct energy performance analyses. The simulation model has been calibrated using actual experimental data collected between October 1, 2019, and January 15, 2020. Results show that the cooling, heating, and lighting energy can be reduced up to 44.1% by fully utilizing electric energy generated by solar PV panels integrated with the DSF system. The payback period is about 15 years based on the energy price provided by the utility company. Although the payback period is long, it has great potential of energy savings and CO 2 emission reductions. The DSF system should be considered as a way of renovation, considering other factors such as increased real estate values and energy cost increases in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗