Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Large‐Eddy Simulations of a Convection Cloud Chamber: Sensitivity to Bin Microphysics and Advection

Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Nonlinear Dynamic Models in Advanced Life Support

To facilitate analysis, ALS systems are often assumed to be linear and time invariant, but they usually have important nonlinear and dynamic aspects. Nonlinear dynamic behavior can be caused by time varying inputs, changes in system parameters, nonlinear system functions, closed loop feedback delays, and limits on buffer storage or processing rates. Dynamic models are usually cataloged according to the number of state variables. The simplest dynamic models are linear, using only integration, multiplication, addition, and subtraction of the state variables. A general linear model with only two state variables can produce all the possible dynamic behavior of linear systems with many state variables, including stability, oscillation, or exponential growth and decay. Linear systems can be described using mathematical analysis. Nonlinear dynamics can be fully explored only by computer simulations of models. Unexpected behavior is produced by simple models having only two or three state variables with simple mathematical relations between them. Closed loop feedback delays are a major source of system instability. Exceeding limits on buffer storage or processing rates forces systems to change operating mode. Different equilibrium points may be reached from different initial conditions. Instead of one stable equilibrium point, the system may have several equilibrium points, oscillate at different frequencies, or even behave chaotically, depending on the system inputs and initial conditions. The frequency spectrum of an output oscillation may contain harmonics and the sums and differences of input frequencies, but it may also contain a stable limit cycle oscillation not related to input frequencies. We must investigate the nonlinear dynamic aspects of advanced life support systems to understand and counter undesirable behavior.

Jones, Harry↗
NASA NTRS1990

Spectral behavior of Hercules X-1 - Its long-term variability and pulse phase spectroscopy

Results are presented of the analysis of the broad-energy-band Her X-1 spectra obtained by HEAO 1. Data from 13 pointed observations covered three 35-day cycles in February, August, and September of 1978. The pulse phase-averaged spectrum did not show any apparent variation during the main-on state, except for the overall intensity. The pulse phase-resolved spectrum, on the other hand, showed a moderately strong phase dependence in each of the main-on state pointings. This may imply that the orientation of the emission region is not changing on a monthly time frame, but the geometry of that region may evolve with time.

Soong, Y.↗
DOE OSTI2025

Longitudinal Patterns in Nitrate Revealed Through Drone‐Based Measurements in an Agriculturally Influenced Midwestern U.S. River

Nitrate concentrations in streams and rivers in the Midwestern United States are often elevated, reflecting the predominance of agriculture in the surrounding landscape. Recent advances in technology, including surface water drones and more precise sensors, provide opportunities to investigate nitrate dynamics with high spatial and temporal resolution. We deployed an aquatic drone, the AquaBOT, in a sixth-order, agriculturally influenced river to examine longitudinal patterns in water quality. Our goal was to measure the spatial and temporal heterogeneity in nitrate and nitrate removal processes and determine the influence of tributary inputs on main stem chemistry. We navigated the drone along a 12-km reach of the Des Moines River (Iowa, USA) nine times between June 2021 and August 2022. Across the deployments, mean nitrate concentration was positively related to discharge and was nearly two orders of magnitude higher in spring than summer. We observed contrasting patterns in main stem nitrate, which decreased downstream during some runs (e.g., 3.1–2.7 mg N L −1 in June 2021), demonstrating net nitrate uptake along the reach, and remained constant on other dates. Similarly, tributaries to the Des Moines had a varied influence on riverine nitrate. Tributaries either increased or decreased main stem nitrate concentrations depending on the tributary and the date. Nitrate removal rates were spatially and temporally variable but showed some consistency at the subreach (2 km) scale, with two subreaches often showing elevated rates of nitrate removal across dates. In conclusion, our study reveals nuanced heterogeneity in nitrate dynamics of the Des Moines River despite the homogeneity of agricultural land cover in the watershed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Comparing the Degree of Land-Atmosphere Interaction in Four Atmospheric General Circulation Models

Land-atmosphere feedback, by which (for example) precipitation-induced moisture anomalies at the land surface affect the overlying atmosphere and thereby the subsequent generation of precipitation, has been examined and quantified with many atmospheric general circulation models (AGCMs). Generally missing from such studies, however, is an indication of the extent to which the simulated feedback strength is model dependent. Four modeling groups have recently performed a highly controlled numerical experiment that allows an objective inter-model comparison of land-atmosphere feedback strength. The experiment essentially consists of an ensemble of simulations in which each member simulation artificially maintains the same time series of surface prognostic variables. Differences in atmospheric behavior between the ensemble members then indicates the degree to which the state of the land surface controls atmospheric processes in that model. A comparison of the four sets of experimental results shows that feedback strength does indeed vary significantly between the AGCMs.

Koster, Randal D.↗
NASA NTRS2013

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗
DOE OSTI2026

Measurement of charged-current quasi elastic-like and inclusive muon (anti-) neutrino interactions using NuMI beam the ICARUS detector

The ICARUS experiment, utilizing Liquid Argon Time Projection Chamber (LArTPC) technology, has been successfully taking physics data at Fermilab since June 2022. The experiment's primary objective is to function as the far detector of the Short Baseline Neutrino program (SBN), searching for hints of physics beyond three-flavour PMNS neutrino oscillations. ICARUS also offers other diverse physics capabilities, including searches beyond the standard model and measurements of cross-sections. In addition to being exposed to the Booster Neutrino (BNB) beamline of the SBN experiment, ICARUS sits 795 m downstream and 5.75 degrees off-axis from the NuMI beam line. Due to the off-axis angle between NuMI and ICARUS, coupled with contributions from both pion and kaon decays to neutrino fluxes, interactions of NuMI neutrinos within ICARUS can be detected over a range of several GeV in energy. Measurements of these interactions present unique opportunities to infer neutrino interaction cross sections on an argon nuclear target within an energy range that overlaps both the SBN oscillation search and a significant portion of the DUNE spectrum. This poster will summarize the current status of muon-neutrino and muon-antineutrino charged current inclusive cross section measurements in ICARUS using the FHC and RHC NuMI sample respectively. Such inclusive measurements provide robust benchmarks for the overall interaction rate and muon kinematics, with reduced dependence on hadronic final-state modeling. The measurement of the muon-antineutrino Argon charged current inclusive cross section is particularly exciting and highly relevant for DUNE, given the scarcity of antineutrino-Argon cross section measurements in the literature. In addition, this poster will highlight ICARUS's first measurement of muon-neutrino Argon charged current mesonless final states using both transverse kinematic imbalance (TKI) variables and muon-proton kinematic variables. These measurements provide sensitivity to nuclear effects, final-state interactions, and nucleon correlations in argon, offering important constraints for improving neutrino interaction modeling for upcoming DUNE.

Roy, Promita [Fermilab; Virginia Tech.]↗
DOE OSTI2021

Differential Organic Carbon Mineralization Responses to Soil Moisture in Three Different Soil Orders Under Mixed Forested System

Soil microbial respiration is one of the largest sources of carbon (C) emissions to the atmosphere in terrestrial ecosystems, which is strongly dependent on multiple environmental variables including soil moisture. Soil moisture content is strongly dependent on soil texture, and the combined effects of texture and moisture on microbial respiration are complex and less explored. Therefore, this study examines the effects of soil moisture on the mineralization of soil organic C Soil organic carbon in three different soils, Ultisol, Alfisol and Vertisol, collected from mixed forests of Georgia, Missouri, and Texas, United States , respectively. A laboratory microcosm experiment was conducted for 90 days under different moisture regimes. Soil respiration was measured weekly, and destructive harvests were conducted at 1, 15, 60, and 90 days after incubation to determine extractable organic C (EOC), phospholipid fatty acid based microbial community, and C-acquiring hydrolytic extracellular enzyme activities (EEA). The highest cumulative respiration in Ultisol was observed at 50% water holding capacity (WHC), in Alfisol at 100% water holding capacity, and in Vertisol at 175% WHC. The trends in Extractable Organic Carbon were opposite to that of cumulative microbial respiration as the moisture levels showing the highest respiration showed the lowest EOC concentration in all soil types. Also, extracellular enzyme activities increased with increase in soil moisture in all soils, however, respiration and EEA showed a decoupled relationship in Ultisol and Alfisol soils. Soil moisture differences did not influence microbial community composition.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Pilot Biofeedback Training in the Cognitive Awareness Training Study (CATS)

One of the ongoing problems that pilots face today is a diminished state of awareness such as boredom, sleepiness, or fatigue during cruise conditions that could result in various pilot errors. This study utilized a cognitive training exercise to sharpen the pilot's awareness during simulated flight thereby providing them with a means to overcome these diminished states of awareness. This study utilizes psychophysiological methods in an attempt to assess a pilot's state of awareness more directly. In turn, the pilots will be able to train themselves to recognize these states of awareness and be more mentally sharp during mundane tasks such as those experienced in cruise conditions. The use of these measurement tools may be beneficial for researchers working within the NASA Aviation Safety Program. This paper will provide the reader with some background information concerning the motivation for the study, a brief description of the experimental setup and design matrix, the dependent and independent variables that were employed, and some preliminary findings based on some of the subjective and objective data that was collected. These preliminary findings are of part of an ongoing study being conducted at the NASA Langley Research Center in Hampton, Virginia.

Uenking, M.↗
NASA NTRS2004

Diels-Alder Trapping of Photochemically Generated o-Quinodimethane Intermediates: An Alternative Route to Photocured Polymer Film Development

Photolysis of o-methylphenyl ketones generates bis-o-quinodimethane intermediates that can be trapped in situ by dienophiles through Diels-Alder cycloadditions. This well-known photochemical process is applied to a series of six new photoreactive monomers containing bis-(o-methylphenyl ketone) functionalities combined with diacrylate and triacrylate ester monomers for the development of acrylic ester copolymer blends. Irradiation of cyclohexanone solutions of the bis-(o-methylphenyl ketone)s and acrylate esters produce thin polymer films. Solid state 13C NMR data indicated 47- 100% reaction of the bis-(o-methylphenyl ketone)s, depending on experimental conditions, to yield the desired products. DSC and TGA analyses were performed to determine the glass transition temperature, T,, and onset of decomposition, Td, of the resulting polymer films. A statistical Design of Experiments approach was used to obtain a systematic understanding of the effects of experimental variables on the extent of polymerization and the final polymer properties.

Tyson, Daniel S.↗
NASA NTRS2019

Analyzing Machine Learning Predictions of Passive Microwave Brightness Temperature Spectral Difference Over Snow-Covered Terrain in High Mountain Asia

Snow is an important component of the terrestrial freshwater budget in high mountainAsia (HMA) and contributes to the runoff in Himalayan rivers through snowmelt. Despitethe importance of snow in HMA, considerable spatiotemporal uncertainty exists across the different estimates of snow water equivalent for this region. In order to better estimate snow water equivalent, radiative transfer models are often used in conjunction with microwave brightness temperature measurements. In this study, the efficacy of support vector machines (SVMs), a machine learning technique, to predict passive microwave brightness temperature spectral difference (1Tb) as a function of geophysical variables (snow water equivalent, snow depth, snow temperature, and snow density) is explored through a sensitivity analysis. The use of machine learning (as opposed to radiative transfer models) is a relatively new and novel approach for improving snow water equivalent estimates. The Noah-MP land surface model within the NASALand Information System framework is used to simulate the hydrologic cycle over HMA and model geophysical variables that are then used for SVM training. The SVMsserve as a nonlinear map between the geophysical space (modeled in Noah-MP) andthe observation space (1Tb as measured by the radiometer). Advanced MicrowaveScanning Radiometer-Earth Observing System measured passive microwave brightness temperatures over snow-covered locations in the HMA region are used as training data during the SVM training phase. Sensitivity of well-trained SVMs to each Noah-MP modeled state variable is assessed by computing normalized sensitivity coefficients. Sensitivity analysis results generally conform with the known first-order physics. Input states that increase volume scattering of microwave radiation, such as snow density and snow water equivalent, exhibit a plurality of positive normalized sensitivity coefficients. In general, snow temperature was the most sensitive input to the SVM predictions. The sensitivity of each state is location and time dependent. The signs of normalized sensitivity coefficients that indicate physical irrationality are ascribed to significant cross-correlation between Noah-MP simulated states and decreased SVM prediction capability at specific locations due to insufficient training data. SVM prediction pitfalls do exist that serve to highlight the limitations of this particular machine learning algorithm.

high mountain Asia↗
NASA NTRSDate not supplied

Improving Computational Efficiency of Prognostics Algorithms in Resource-Constrained Settings

In engineering and aerospace applications, it is vital to operational success to have insight into the expected performance and health of physical systems. The field of prognostics and health management provides quantitative methods for monitoring, predicting, and managing system health. Prognostics algorithms can be employed to assess the current state of a system, propagate the system throughout time, and predict potential anomalies or failures that may occur. While they can provide accurate prediction results, effective prognostics algorithms can be challenging to use in resource-constrained settings due to computational limitations and high computational latency, leading to obsolete predictions. Thus, computationally efficient and accurate algorithms are necessary for future remaining useful life predictions. In this work, we implement new algorithmic approaches for prediction, quantitatively compare them via a battery degradation use-case, and provide recommendations of potential improvements to a prognostics framework. One approach to prediction is through sampling, whereby the current state of a physical system is sampled many times and each sample is propagated forward until failure is reached, resulting in a distribution of failure values. To improve the efficiency of this process, we implemented five new algorithmic approaches to prediction, including three distinct sampling methods (standard Monte Carlo, Quasi-Monte Carlo, and Latin Hypercube Sampling), a variable time step algorithm, and a variable sample size algorithm. To compare the algorithms, we employ a variety of metrics designed specifically to analyze both computational efficiency and model accuracy. Our metrics include accuracy to compare the average predicted value to ground truth, mean absolute deviation to illustrate dispersion, specific percentile error to describe accuracy within a user-defined risk tolerance, and code run-time. To quantitatively analyze our results, we employ a use-case of degradation of a Lithium-ion battery. We use an electrochemistry-based model to describe the current health state of the battery, and implement our prediction algorithms to propagate forward in time until end-of-discharge (EOD) is reached. Notably, through this work it was found that none of our sampling approaches had a significant impact on computational efficiency or model accuracy in predicting EOD of the battery. We find that while the sampling methods are unique, the distributions they generate are similar, ultimately producing final predictions that are nearly identical. In exploring the effect of the time step within the prediction algorithm, we found that prediction accuracy was highly dependent on the time step used, and that implementing a variable time step within a particular prediction may provide an increase in computational efficiency while also maintaining prediction accuracy. Finally, implementing a variable sample size also affected prediction, and our results show that tuning both the magnitude and timing of the sample size adjustment can result in improved computation speed and maintained prediction accuracy. Taken together, our findings highlight the challenge of performing prognostics in resource-constrained settings, and illustrate the potential of developing new prediction algorithms to improve computational efficiency.

prognostics↗
DOE OSTI2022

Role of ripples in altering the electronic and chemical properties of graphene

Ripples of graphene are known to manipulate electronic and hydrogenation properties of graphitic materials. More detailed work is needed to elucidate the structure–property relationship of these systems. In this work, the density functional theory is used to compute the energy and electronic structure of the graphene models with respect to variable curvatures and hydrogen adsorption sites. The magnitude of finite bandgap opening depends on the orientation of ripples, and the hydrogen adsorption energy depends on the local curvature of graphene. An adsorbed hydrogen alters the local curvature, resulting in relatively weakened adsorption on the neighboring three sites, which gives a rationale to experimentally observed dynamic equilibrium stoichiometry (H:C = 1:4) of hydrogenated graphene. The surface diffusion transition state energy of adsorbed hydrogen is computed, which indicates that the Eley–Rideal surface recombination mechanism may be important to establish the dynamic equilibrium, instead of the commonly assumed Langmuir–Hinshelwood mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

A model for predicting damage induced fatigue life of laminated composite structural components

This paper presents a model for predicting the life of laminated composite structural components subjected to fatigue induced microstructural damage. The model uses the concept of continuum damage mechanics, wherein the effects of microcracks are incorporated into a damage dependent lamination theory instead of treating each crack as an internal boundary. Internal variables are formulated to account for the effects of both matrix cracks and internal delaminations. Evolution laws for determining the damage variables as functions of ply stresses are proposed, and comparisons of predicted damage evolution are made to experiment. In addition, predicted stiffness losses, as well as ply stresses are shown as functions of damage state for a variety of stacking sequences.

Allen, David H.↗
NASA NTRS2018

Modeling Principles Using the Relation Between Lagrange Multipliers and Bond Graphs

Final document is attached. Modeling dynamic systems by bond graphs has become state of the art technology since hundreds of researchers around the world have incorporated the technology in many fields of engineering and science. The legacy of its invertor Prof. Henry Paynter at MIT in 1959 is now a fundamental and practical technique to understand reality by building computer models. This paper addresses a particular aspect of this technology when modeling of mechanical systems require relaxation of constraints by means of Lagrange principles. Lagrange's equations are a useful means of describing and solving systems with kinematic constraints. Lagrange multipliers are variables used in equations to find the extremes of multivariate functions. Here we explore the relation of Lagrange multipliers to solve modeling difficulties of a space vehicle with equations with dependent derivatives. Lagrange multipliers were used in conjunction with bond graphs to simulate a system where joints of kinematic linkages produce dependent derivatives. NASA's Morpheus Project lunar lander was used as a case study. The Morpheus Project is a terrestrial test vehicle designed to fly the terminal descent trajectory of a lunar lander to advance the Autonomous Landing Hazard Avoidance Technology (ALHAT). An objective of this study is to apply the modeling approach herein to capture the dynamic movement of the lander as the propellant is sloshed and consumed. This paper expands further the analysis presented by (Granda, J J. Nguyen, L, Carlson, T, Sahragard-Monfared, G., Fornalski, E., Brocker 2016). Using an automated approach bond graph models of state space equations were generated using the Computer Aided Modeling Program (CAMPG). Integral causality models and derivative causality models were considered in order to find the simpler solution for the mathematical dependencies produced in modeling this vehicle.

Granda, Jose J.↗
NASA NTRS2016

An Evaluation of MODIS-Retrieved Aerosol Optical Depth over a Mountainous AERONET Site in the Southeastern US

The literature shows that aerosol optical depth (AOD) derived from the MODIS Collection 5 (C5) dark target algorithm has been extensively validated by spatiotemporal collocation with AERONET sites on both global and regional scales.Although generally comparing well over the eastern US region, poor performance over mountains in other regions indicate the need to evaluate the MODIS product over a mountain site. This study compares MODIS C5 AOD at 550nm to AOD measured at the Appalachian State University AERONET site in Boone, NC over 30 months between August 2010 and September 2013. For the combined Aqua and Terra datasets, although more than 70% of the 500 MODIS AOD measurements agree with collocated AERONET AOD to within error envelope of +/- (0.05 + 15%), MODIS tends to have a low bias (0.02-0.03). The agreement between MODIS and AERONET AOD does not depend on MODIS quality assurance confidence (QAC) value. However, when stratified by satellite, MODIS-Terra data does not perform as well as Aqua, with especially poor correlation (r = 0.39) for low aerosol loading conditions (AERONET AOD less than 0.15).Linear regressions between Terra and AERONET possess statistically-different slopes for AOD < 0.15 and AOD > or = 0.15. AERONET AOD measured only during MODIS overpass hours is highly correlated with daily-averaged AERONET AOD. MODIS monthly-averaged AOD also tracks that of AERONET over the study period. These results indicate that MODIS is sensitive to the day-to-day variability, as well as the annual cycle of AOD over the Appalachian State AERONET site. The complex topography and high seasonality in AOD and vegetation indices allow us to specifically evaluate MODIS dark target algorithm surface albedo and aerosol model assumptions at a regionally-representative SE US mountain site.

air quality↗
DOE OSTI2025

Tuning water dissociation at oxide–electrolyte interfaces with electric fields

Understanding how electric fields influence water dissociation at heterogeneous interfaces is crucial for controlling interfacial chemical reactions and advancing next-generation energy technologies. Herein, ab initio–based machine learning simulations show that even small electric field changes can significantly alter the water dissociation fraction at planar TiO 2 –electrolyte interfaces. The resulting free energy difference between undissociated and dissociated interfacial water exhibits a linear dependence on the field change with a slope of 1.97 eÅ, which far exceeds the dissociation-induced dipole change of a water molecule. Employing a machine-learned collective variable to investigate the reaction statistics of thousands of water dissociation/recombination events, we find that small electric field changes exert minor effects on individual reaction energy barriers but significantly influence the populations of local configurations associated with initial states that are most favorable for reactions. These findings elucidate the pronounced impact of electric fields on interfacial water dissociation and reveal a mechanism for electric-field-controlled chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗