Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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 199 records · Page 11

Faster Form of Electron Magnetic Reconnection with a Finite Length X-Line

Observations in Earth’s turbulent magnetosheath downstream of a quasiparallel bow shock reveal a prevalence of electron-scale current sheets favorable for electron-only reconnection where ions are not coupled to the reconnecting magnetic fields. In small-scale turbulence, magnetic structures associated with intense current sheets are limited in all dimensions. And since the coupling of ions are constrained by a minimum length scale, the dynamics of electron reconnection is likely to be 3D. Here, both 2D and 3D kinetic particle-in-cell simulations are used to investigate electron-only reconnection, focusing on the reconnection rate and associated electron flows. A new form of 3D electron-only reconnection spontaneously develops where the magnetic X-line is localized in the out-of-plane (z) direction. The consequence is an enhancement of the reconnection rate compared with two dimensions, which results from differential mass flux out of the diffusion region along z, enabling a faster inflow velocity and thus a larger reconnection rate. This outflow along z is due to the magnetic tension force in z just as the conventional exhaust tension force, allowing particles to leave the diffusion region efficiently along z unlike the 2D configuration.

79 ASTRONOMY AND ASTROPHYSICS↗

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗

Avalanches and many-body resonances in many-body localized systems

Here we numerically study both the avalanche instability and many-body resonances in strongly disordered spin chains exhibiting many-body localization (MBL). Finite-size systems behave like MBL within the MBL regimes, which we divide into the asymptotic MBL phase and the finite-size MBL regime; the latter regime is, however, thermal in the limit of large systems and long times. In both Floquet and Hamiltonian models, we identify some landmarks within the MBL regimes. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition occurs much deeper in the MBL regimes than has been suggested by most previous studies. Our other landmarks involve systemwide many-body resonances: We find that the effective matrix elements producing eigenstates with systemwide many-body resonances are enormously broadly distributed. This broad distribution means that the onset of such resonances in typical samples occurs quite deep in the MBL regimes, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regimes of finite-size systems into three subregimes: (i) at strongest randomness, typical samples do not have any eigenstates that are involved in systemwide many-body resonances; (ii) there is a substantial intermediate subregime where typical samples do have such resonances but the pair of eigenstates with the minimum spectral gap does not, so the size of the minimum gap agrees with expectations from Poisson statistics; and (iii) in the weaker randomness subregime, the minimum gap is larger than predicted by Poisson level statistics because it is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain nonresonant and the system thus still appears MBL in most respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i) and the other subregimes are entirely in the thermal phase, even though they look localized in most respects, so are in the finite-size MBL regime.

36 MATERIALS SCIENCE↗

FREQUENCY REGULATION WITH CONNECTED LIGHTING SYSTEMS

The management of end-use energy resources, e.g., commercial buildings, has been investigated as a promising source of services for the electric power grid. Lighting consumes about 17% electricity of U.S. commercial buildings; however, it may contribute significantly to services that improve the reliability and resilience of the grid due to its rapid speed of response. Connected lighting systems (CLS), which build upon solid-state lamp technology, can change their power demand more quickly than most other building electricity end uses. But the potential of CLS to provide grid services has not been considerably investigated. In this paper, we describe initial research to evaluate the potential of CLS for the grid service of frequency regulation. Frequency regulation is a reliability product that corrects in a matter of seconds for short-term changes in the balance between supply and demand in the balancing area that might affect the stability of the power system. Frequency regulation signals for a medium office building are generated from the normalized test signals for the PJM Interconnection Reg-A and Reg-D regulations services. Control of the CLS locally to follow the Reg-A or Reg-D signal is used to provide contributions to the corresponding frequency regulation service. The performance of CLS for frequency regulation is evaluated using the PJM 40-Minute Performance Score Template. The performance scores obtained for all CLS categories for both Reg-A and Reg-D signals far exceed the minimum score for the qualification, a very promising result for CLS to provide frequency regulation service.

Wang, Peng↗

The influence of LiH and TiH 2 on hydrogen storage in MgB 2 I: Promotion of bulk hydrogenation at reduced temperature

Mg(BH 4 ) 2 is an attractive hydrogen storage material, owing to its high gravimetric capacity of 14.9 wt %. However, the dehydrogenated material MgB 2 is very difficult to rehydrogenate, requiring excessive pressures and temperatures. Here we report the influence of LiH and TiH 2 on hydrogen storage reactions involving Bulk MgB 2 using XRD, XAS, FTIR and NMR. In ball-milled mixtures of LiH/MgB 2 , the LiH loses crystallinity but remains undissociated, forming a weakly bound complex with MgB 2 . The weak interactions produce minor variations in the local electronic structure at B and Mg, but do not markedly affect the underlying MgB 2 hexagonal crystal structure. No evidence is found for a mixed-metal boride Mg 1-x Li x B 2 in the as-prepared LiH/MgB 2 materials. The presence of LiH dramatically improves the hydrogenation of MgB 2 at 700 bar, forming borohydride 100 °C below the minimum hydrogenation temperature of pure MgB 2 and without the formation of undesirable intermediates such as [B 3 H 8 ] - , [B 10 H 10 ] 2- or [B 12 H 12 ] 2- . Evidence is reported for a mixed-metal borohydride of the type Mg (3-x)/2 Li x (BH 4 ) 3 produced by the hydrogenation. Subsequent desorption is also improved compared to pure Mg(BH 4 ) 2 and LiBH 4 , showing single-step hydrogen release up to ~8 wt% by 380 °C, whereas Mg(BH 4 ) 2 and LiBH 4 still retain significant amounts of hydrogen at this temperature. The material produced by desorption contains both MgB 2 and Mg metal, revealing the original LiH/MgB 2 system is not fully reversible. In contrast to LiH, TiH 2 is essentially inert when ball-milled with MgB 2 , and high-pressure hydrogenation leaves only unreacted TiH 2 and MgB 2 . Thus, added TiH 2 provides no benefit to MgB 2 hydrogenation.

08 HYDROGEN↗

Quasi-linear toroidal simulations of resonant magnetic perturbations in eight ITER H-mode scenarios

Abstract Both linear and quasi-linear aspects of the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for various H-mode scenarios in ITER, covering the pre-fusion power operation and the fusion power operation phases. Linear response computations for eight ITER scenarios, with varying plasma current and toroidal magnetic field, reveal that the best coil current phasing for controlling the type-I edge localized modes (ELMs) scales roughly linearly with the edge safety factor. The coil phasing is defined as the relative toroidal phase of the coil currents between different rows, for a given toroidal harmonic. Quasi-linear initial value simulation, which is the focus of the present study, shows that application of the n = 3 ( n is the toroidal mode number) RMP field has a minimum side effect on the plasma core momentum confinement but potentially a large effect on the global particle transport. Generally, the RMP field with the best (worst) coil phasing for ELM control produces the strongest (weakest) effect on the plasma edge flow and the overall density. This robustly holds for all eight ITER scenarios. Consequently, in order to minimize the RMP induced side effects while achieving ELM control (suppression) in ITER, a compromise is necessary in choosing the coil current configuration.

Physics↗

Volt-VAR Optimization in Distribution Networks Using Twin Delayed Deep Reinforcement Learning

Modern distribution grids are undergoing new challenges due to the stochastic nature of distributed energy resources (DERs). High penetration of DERs has a significant impact on Volt-VAR profile and system power losses. This work proposes a deep reinforcement learning (DRL)-based Volt-VAR optimization approach for improving voltage profile and reducing system power loss under high penetration of distributed energy resources, such as battery energy storage and solar photovoltaic units in distribution grids. The twin delayed deep deterministic policy gradient (TD3) method-based DRL agent is proposed to configure optimal set-points of reactive power outputs of fast responding smart inverters. The agent schedules the reactive power of inverters according to their physical capabilities, such as minimum allowed power factor, e.g., 0.9 leading/lagging. The reward function of the proposed DRL scheme is designed carefully to ensure a proper voltage profile of the grids with effective scheduling of reactive power outputs from inverters. The performance of the proposed model is verified on modified IEEE 34- and 123-bus systems and compared with base case with no reactive supply by inverters, and local droop Volt-VAR control approach. The results show that the proposed method performs better than the local droop control and deep deterministic policy gradient (DDPG)-based DRL method for reducing voltage fluctuation and minimizing power loss.

Hossain, Rakib↗

Refining seasonal performance metrics for room air-conditioning in emerging markets: Integrating building simulations with real-world equipment performance data

Buildings significantly impact worldwide energy consumption, emphasizing the need to reduce the cooling energy demand, especially in warm climates. Minimum energy performance standards (MEPS) and seasonal performance metrics such as the Cooling Seasonal Performance Factor (CSPF) are crucial for improving room air conditioning (RAC) efficiency. However, challenges remain, particularly in emerging markets like Brazil, where seasonal performance metrics have recently been introduced. This study assesses the factors influencing country-level seasonal efficiency metrics and proposes a framework to refine these calculations by considering local climates and expected RAC usage in real-world households via building simulations. Key considerations include outdoor air temperature binning for different climates, RAC usage patterns (i.e., daytime and nighttime usages), envelope thermal performance of households, and urban heat island (UHI) effects. The results reveal that CSPF values can vary significantly based on climate conditions, with observed CSPF ranging from 4.10 to 11.59 Wh/Wh across 577 Brazilian climates. The inclusion of UHI effects led to a reduction in CSPF values by up to 29% during nighttime operations in hot urban areas. Additionally, building envelope efficiency showed contrasting impacts on RAC performance, with CSPFs reaching up to 15.35 Wh/Wh under specific optimized conditions. These findings highlight the need for transparent policymaking in RAC performance databases, facilitating the application of approaches like those proposed in this study and supporting diverse stakeholders in decision-making.

Bavaresco, Mateus↗

A Nanoscale Ternary Amide‐rGO Composite with Boosted Kinetics for Reversible H 2 Storage

Abstract Metal amides are attractive candidates for hydrogen storage due to their high volumetric and gravimetric hydrogen densities. However, the sluggish kinetics and competing side reactions during hydrogen uptake and release limit their practical use. Here, a novel nanoconfined Li 2 Mg(NH) 2 @reduced graphene oxide (rGO) composite is presented, which is fabricated using a melt‐infiltration method with a minimum weight penalty of only 2 wt.%. The presence of rGO ensures close contact between the active phases, effectively preventing aggregation during cycling process. As a result, the reversible capacity of Li 2 Mg(NH) 2 @rGO reaches 4.42 wt.%, with no capacity degradation observed after multiple cycling. Theoretical calculations show that rGO catalyzes the hydrogen bond cleavage at the Mg‐amide/Li hydride interface, leading to local dehydrogenation hotspots and significantly improves kinetics of dehydrogenation compared to the bulk counterpart. This study provides a promising strategy for designing metal imide‐based composites to overcome the kinetic limitations and improve their reversible hydrogen storage performance.

08 HYDROGEN↗

Fabrication of silicon-tipped fiber-optic temperature sensors using aerogel-assisted glass soldering with precise laser heating

We demonstrate the fabrication of fiber-optic Fabry–Perot interferometer (FPI) temperature sensors by bonding a small silicon diaphragm to the tip of an optical fiber using low melting point glass powders heated by a 980 nm laser on an aerogel substrate. The heating laser is delivered to the silicon FPI using an optical fiber, while the silicon temperature is being monitored using a 1550 nm white-light system, providing localized heating with precise temperature control. The use of an aerogel substrate greatly improves the heating efficiency by reducing the thermal loss of the bonding parts to the ambient environment. A desirable temperature for bonding can be achieved with relatively small heating laser power. The bonding process is carried out in an open space at room temperature for convenient optical alignment. The precise temperature control ensures minimum perturbation to the optical alignment and no induced thermal damage to the optical parts during the bonding process. For demonstration, we fabricated a low-finesse and high-finesse silicon FPI sensor and characterized their measurement resolution and temperature capability. The results show that the fabrication method has a good potential for high-precision fabrication of fiber-optic sensors.

Sheng, Qiwen↗

A Solution Method for the Filtered Lifting Line Theory

The filtered lifting line theory presents a continuous form of the inviscid momentum equations of flow over a lifting device, such as a wing or rotor blade, using body forces without mathematical singularities. This theory is also consistent with an actuator line representation of a lifting device. In this work, we present a reformulation of the equations in terms of the local flow angle along the line, which allows solving the stand-alone equations using multivariate root-finding algorithms. This approach can be used to obtain a fast, computationally inexpensive solution of the loading distribution along a wing without the need to perform computational fluid dynamic simulations. We study the requirements in terms of resolution in the spanwise direction and establish the criteria for spacing and minimum amount of points required along the blade to obtain converged solutions. The solutions are compared to results from large-eddy simulations, and we observed excellent agreement with less than a percent difference in quantities along the blade between the methods.

17 WIND ENERGY↗

Transition Core Planning and Safety Analyses in Support of LEU Fuel Conversion of the University of Missouri Research Reactor (MURR)

The University of Missouri Research Reactor (MURR®) is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is working with the National Nuclear Security Administration (NNSA) Office of Material Management and Minimization (M3) Reactor Conversion Program to convert from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. The M3 Reactor Conversion USHPRR Project objectives include the development of LEU fuel element designs that will ensure safe reactor operations and to maintain the existing experimental performance of each facility. The work is being conducted through many inter-related activities being completed by four Project Pillars: Fuel Qualification (FQ), Fuel Fabrication (FF), Reactor Conversion (RC), and Cross Cutting (CC). A new type of LEU fuel based on an alloy of uranium-10 wt% molybdenum (U-10Mo) is expected to allow the conversion of those USHPRR, like MURR, requiring higher density fuels. The very-high-density LEU U-10Mo monolithic fuel is currently undergoing irradiation testing and post-irradiation examination under a planned and documented fuel qualification effort. The FQ Pillar will document fuel property and fuel performance data and qualify the fuel for use in these reactors. The FF Pillar is fabricating fuel for ongoing and future irradiation tests, as well as conducting fabrication demonstrations to validate or update preliminary fabrication assumptions. The FF Pillar is also working to develop and install commercial manufacturing capacity with the U-10Mo monolithic fuel to produce prototypic fuel. Working with the RC Pillar at Argonne, MURR has progressed through a preliminary fuel element design using preliminary data for the proposed monolithic alloy of U-10Mo. Analyses were completed in previous work that found for typical equilibrium operations with the preliminary LEU fuel element design, in conjunction with a power uprate to 12 MW and appropriate changes to the MURR Limiting safety system settings (LSSS), MURR will have adequate margins to safety for steady-state operations and postulated transient accidents and will have experimental performance in key locations that meets or exceeds current operations with HEU fuel. The purpose of this work is to develop a sequence of transition cycles that will enable MURR to transition from operation with the reactor core loaded with fresh LEU fuel elements only to typical equilibrium operations with mixed-burnup cores following conversion while meeting operational requirements on safety and experimental performance. It is expected that the use of fresh LEU fuel at conversion and subsequent low burnup of the LEU fuel elements that will initially be available for use following conversion will result in critical control blade positions that will substantially change the axial power distribution in the core and the neutron flux available in key experimental locations relative to equilibrium LEU operations. Given the constraints of MURR safety margins, operational practices, and production and research, a novel method has been developed to identify a transition sequence that minimizes the time MURR operates atypically compared to the current prototypic cycles using HEU fuel. The proposed transition sequence moves quickly to the same sort of equilibrium cycles for the LEU fuel that have already been evaluated in documented preliminary safety analyses. Although shifting the neutron flux peak to the lower half of the core during initial cycles with LEU at 12 MW reduces the experiment performance in some key locations relative to current HEU operations at 10 MW, all LEU cores provide an average performance that meets or exceeds that of HEU. An LEU cycle is reached that meets or exceeds the level of experimental performance predicted for current HEU and equilibrium LEU operations in more than 450 key locations identified by a reactor specialist at MURR by the 23rd cycle following conversion and that afterwards will enable MURR to consistently meet its experimental performance requirements. The proposed transition sequence only requires the fabrication of 34 fresh LEU elements in the first year of operation and does not exceed the anticipated availability of fresh elements that can be produced by the fuel fabricator. By the third year after conversion, 22 fresh LEU elements will be required each year, which is the same as expected for equilibrium LEU operations and the same as current operations with HEU fuel. The proposed transition sequence thus combines a relatively short time period before equilibrium burnup is achieved, a temporary increase of fuel elements needed annually relative to typical operations that are within the production capabilities of the fuel fabricator, and demonstrates comparable experimental performance of the LEU cores relative to current HEU operations. Further measures may be taken to reduce any initial experimental performance penalty even further, where possible, by repositioning certain experiments to leverage the increased performance in the lower axial experimental positions in the initial cycles following conversion or leaving the experiments in the irradiation facilities longer in order to achieve the required neutron fluence. This analysis may require refinement depending on the experimental facilities in use at the time of conversion. Nonetheless, the results presented here, including the experimental performance, core burnup, and critical control blade positions throughout the transition cycles, show that the proposed transition cycle fuel management patterns are consistent with what is expected and desired for MURR operation with LEU U-10Mo fuel. Detailed core power distributions from the neutronics models were also used to evaluate safety margins during steady-state operations for the selected transition cycles and the equilibrium LEU core. It is shown that there are adequate safety margins for both steady-state operations and postulated accident scenarios. For the steady-state operations with the preliminary LEU fuel element design the analysis predicts at least 2.49 MW margin to the onset of flow instability at the LSSS power of 15 MW. Considering the LSSS power is 125% of full license power, the margin to OFI is sufficient. In addition, the critical heat flux ratio at LSSS power is well above the requirement of CHFR > 2.0 from NUREG-1537 for all considered cases. For postulated transient accidents, the minimum margin to the fuel temperature safety limit is at least 109 °C. In summary, the proposed sequence of core loadings for MURR operations following conversion to LEU fuel and a power uprate to 12 MW provides sufficient safety margins for both steady-state operations and postulated transient accidents during a proposed sequence of transition cycles to equilibrium operations. Analysis has shown that there are some local experimental performance penalties during the initial cycles. Although there are local shifts in the experimental performance, on average all LEU cores at 12 MW have equal or higher performance than HEU at 10 MW. Temporary adjustments are being planned that will produce suitable experimental performance during these cycles. The results indicate that for the equilibrium LEU core the experimental performance exceeds that of current HEU operations in all key locations while also demonstrating sufficient safety margins.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase diagram of infinite layer praseodymium nickelate Pr 1– x Sr x NiO 2 thin films

In this paper, we report the phase diagram of infinite layer Pr 1-x Sr x NiO 2 thin films synthesized via topotactic reduction from the perovskite precursor phase using CaH 2 . Based on the electrical transport properties, we find a doping-dependent superconducting dome extending between x = 0.12 and 0.28, with a maximum superconducting transition temperature T c of 14 K at x = 0.18, bounded by weakly insulating behavior on both sides. In contrast to the narrower dome observed in Nd 1–x Sr x NiO 2 , a local T c suppression near x = 0.2 was not observed for the Pr 1-x Sr x NiO 2 system. Normal state Hall effect measurements indicate mixed carrier contributions of both electrons and holes, and show a sign change in the Hall coefficient as functions of temperature and x, quite similar to that in Nd 1-x Sr x NiO 2 . Also similar is the observation of a minimum in the normal state resistivity associated with the superconducting compositions. These findings indicate an infinite layer nickelate phase diagram that is relatively insensitive to the rare-earth element, but suggest that disorder arising from the variations of the ionic radii on the rare-earth site affects the superconducting dome.

36 MATERIALS SCIENCE↗

An airfoil-based synthetic actuator disk model for wind turbine aerodynamic and structural analysis

Here, this study introduces an airfoil-based refinement technique to enhance the Actuator Disk Model (ADM) for improved wind turbine aerodynamic load prediction and structural simulation in conjunction with Large Eddy Simulations of the wind flow. While ADM offers higher computational efficiency than the more detailed but resource-intensive Actuator Line Model (ALM), it traditionally lacks the resolution needed to capture the localized blade forces accurately. To address this limitation, we introduce a refinement technique that uses airfoil-specific data and employs interpolation-based grid point refinement, achieving ALM-comparable accuracy while preserving ADM's efficiency. Unlike conventional ADM that provides only rotor-disk averaged forces, our synthetic method tracks transient aerodynamic load variations over multiple blade revolutions, allowing us to calculate the distributions of maximum and minimum loads during typical cycles. Applied to the NREL 5 MW reference turbine, our enhanced ADM accurately predicts key aerodynamic parameters (angle of attack, axial velocity, lift, drag, axial and tangential forces along the blades) as well as structural responses (blade tip deflection, maximum stress, and stress concentration). Our results show that the tip deflection ranges from 2.33m (3.69 % of blade length) to 4.28m (6.79 %), with maximum stress concentration occurring near the blade root. This research demonstrates that a refined synthetic ADM approach can serve as a computationally efficient alternative for both aerodynamic analysis and structural simulation of wind turbine blades subjected to realistic wind fields.

17 WIND ENERGY↗

Deep reinforcement learning for predicting kinetic pathways to surface reconstruction in a ternary alloy

Abstract The majority of computational catalyst design focuses on the screening of material components and alloy composition to optimize selectivity and activity for a given reaction. However, predicting the metastability of the alloy catalyst surface at realistic operating conditions requires an extensive sampling of possible surface reconstructions and their associated kinetic pathways. We present CatGym, a deep reinforcement learning (DRL) environment for predicting the thermal surface reconstruction pathways and their associated kinetic barriers in crystalline solids under reaction conditions. The DRL agent iteratively changes the positions of atoms in the near-surface region to generate kinetic pathways to accessible local minima involving changes in the surface compositions. We showcase our agent by predicting the surface reconstruction pathways of a ternary Ni 3 Pd 3 Au 2 (111) alloy catalyst. Our results show that the DRL agent can not only explore more diverse surface compositions than the conventional minima hopping method, but also generate the kinetic surface reconstruction pathways. We further demonstrate that the kinetic pathway to a global minimum energy surface composition and its associated transition state predicted by our agent is in good agreement with the minimum energy path predicted by nudged elastic band calculations.

36 MATERIALS SCIENCE↗

Superconvergence of Online Optimization for Model Predictive Control

We develop a one-Newton-step-per-horizon, online, lag-L, model predictive control (MPC) algorithm for solving discrete-time, equality-constrained, nonlinear dynamic programs. Based on recent sensitivity analysis results for the target problems class, we prove that the approach exhibits a behavior that we call superconvergence; that is, the tracking error with respect to the full horizon solution is not only stable for successive horizon shifts, but also decreases with increasing shift order to a minimum value that decays exponentially in the length of the receding horizon. The key analytical step is the decomposition of the one-step error recursion of our algorithm into algorithmic error and perturbation error. We show that the perturbation error decays exponentially with the lag between two consecutive receding horizons, while the algorithmic error, determined by Newton’s method, achieves quadratic convergence instead. Overall this approach induces our local exponential convergence result in terms of the receding horizon length for suitable values of L. In conclusion, numerical experiments validate our theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗