Search NASA⌕ Search

SEARCH · Search NASA

Results for “dominant balance”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shear Stress Dependence of Force Networks in 3D Dense Suspensions

The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicate a uniformly spreading network that rapidly forms across empty domains, however the growth patterns have structure that are significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening.

Edens, Lance E.↗
DOE OSTI2024

Hysteresis between gas breakdown and plasma discharge

In direct-current (DC) discharge, it is well known that hysteresis is observed between the Townsend (gas breakdown) and glow regimes. Forward and backward voltage sweep is performed using a one-dimensional particle-in-cell Monte Carlo collision (PIC-MCC) model considering a ballast resistor. When increasing the applied voltage after reaching the breakdown voltage (V b ), transition from Townsend to glow discharges is observed. When decreasing the applied voltage from the glow regime, the discharge voltage (V d ) between the anode–cathode gap can be smaller than the breakdown voltage, resulting in a hysteresis, which is consistent with experimental observations. Next, the PIC-MCC model is used to investigate the self-sustaining voltage (V s ) in the presence of finite initial plasma densities between the anode and cathode gap. It is observed that the self-sustaining voltage coincides with the discharge voltage obtained from the backward voltage sweep. In addition, the self-sustaining voltage decreases with increased initial plasma density and saturates above a certain initial plasma density, which indicates a change in plasma resistivity. The decrease in self-sustaining voltage is associated with the electron heat loss at the anode for the low pd (rarefied) regime. In the high pd (collisional) regime, the ion energy loss toward the cathode due to the cathode fall and the inelastic collision loss of electrons in the bulk discharge balance out. Finally, it is demonstrated that the self-sustaining voltage collapses to a singular value, despite the presence of a initial plasma, for microgaps when field emission is dominant, which is also consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resolving Salt-Induced Agglomeration of Laponite Suspensions Using X-ray Photon Correlation Spectroscopy and Molecular Dynamics Simulations

Linking the physics of the relaxation behavior of viscoelastic fluids as they form arrested gel states to the underlying chemical changes is essential for developing predictive controls on the properties of the suspensions. In this study, 3 wt.% laponite suspensions are studied as model systems to probe the influence of salt-induced relaxation behavior arising from the assembly of laponite nanodisks. X-ray Photon Correlation Spectroscopy (XPCS) measurements show that laponite suspensions prepared in the presence of 5 mM concentrations of CaCl2, MgCl2 and CsCl salts accelerate the formation of arrested gel states, with CaCl2 having a significant impact followed by CsCl and MgCl2 salts. The competing effects of ion size and charge on relaxation behavior are noted. For example, the relaxation times of laponite suspensions in the presence of Mg2+ ions are slower compared to Cs+ ions despite the higher charge, suggesting that cation size dominates in this scenario. The faster relaxation behavior of laponite suspensions in the presence of Ca2+ ions compared to Cs+ ions shows that a higher charge dominates the size of the ion. The trends in relaxation behavior are consistent with the cluster formation behavior of laponite suspensions and the electrostatic interactions predicted from MD simulations. Charge balance is achieved by the intercalation of the cations at the negatively charged surfaces of laponite suspensions. These studies show that the arrested gel state of laponite suspensions is accelerated in the presence of salts, with ion sizes and charges having a competing effect on relaxation behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

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

Arbuscular mycorrhizal fungi equalize differences in plant fitness and facilitate plant species coexistence through niche differentiation

Mycorrhizal fungi are essential to the establishment of the vast majority of plant species but are often conceptualized with contradictory roles in plant community assembly. On the one hand, host-specific mycorrhizal fungi may allow a plant to be competitively dominant by enhancing growth. On the other hand, host-specific mycorrhizal fungi with different functional capabilities may increase nutrient niche partitioning, allowing plant species to coexist. Here, to resolve the balance of these two contradictory forces, we used a controlled greenhouse study to manipulate the presence of two main types of mycorrhizal fungus, ectomycorrhizal fungi and arbuscular mycorrhizal fungi, and used a range of conspecific and heterospecific competitor densities to investigate the role of mycorrhizal fungi in plant competition and coexistence. We find that the presence of arbuscular mycorrhizal fungi equalizes fitness differences between plants and stabilizes competition to create conditions for host species coexistence. Furthermore, our results show how below-ground mutualisms can shift outcomes of plant competition and that a holistic view of plant communities that incorporates their mycorrhizal partners is important in predicting plant community dynamics.

09 BIOMASS FUELS↗
DOE OSTI2020

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗
DOE OSTI2021

Coordination and control – limits in standard representations of multi-reservoir operations in hydrological modeling

Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Influence of interfacial responses of Berea Sandstone in low-salinity waterflooding environments

Wetting conditions and fluid-fluid interfacial properties are essential to understanding oil recovery in the context of low-salinity and smart waterflooding. Many have deemed wettability alteration to be the dominant interaction, while some highlight the importance of fluid-fluid interactions. In some rock-fluid systems, benefits from the use of additives can arise to favor hydrocarbon production. Understanding the balance between the fluid-fluid and rock-fluid interactions when utilizing additives is crucial to evaluate benefits in improved efficiency. Ultimately, this understanding will shed light on recovery mechanisms. Regarding oil recovery, previous studies indicate that the injection of a naphthenic acid (NA) blend dissolved in an injection brine improves oil recovery in Berea Sandstone, particularly under low-salinity conditions. Interfacial rheology results for the brine-oil interface correlated well with the recovery response. This observation was speculated to serve as a potential improved recovery mechanism through oil-trapping suppression. However, the said studies did not evaluate potentially concurrent rock wettability alteration during recovery or the initial wetting condition. The purpose of this work is to determine the production behavior in systems with NAs, namely NA-10 (3-Cyclohexanepropinoic acid), NA-14 (Cyclopentane carboxylic acid), NA-18 (Cyclohexane carboxylic acid), and NA-mix (Naphthenic acid mixture), at different initial wetting conditions under spontaneous imbibition. In addition, we tested the NAs ability to alter wettability in Berea Sandstone and determine each “additive’s” capacity to improve recovery. Here, time-domain nuclear magnetic resonance (TD-NMR), contact angle, and Amott cell tests, i.e. spontaneous imbibition, were used to analyze the influence of the selected naphthenic acids on oil displacement under different initial wetting conditions. Our results indicate that individual naphthenic acids indeed have potential as additives in low-salinity waterflooding for oil-wet systems while negatively impacting recovery in water-wet systems. However, results suggest that the evolution of the interfacial properties and wettability conditions are equally important on the recovery behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An accelerated approach to inline non-LTE modeling

A goal across many physics and engineering disciplines is to use radiation-hydrodynamics codes to accurately simulate systems that are expected to deviate from local thermodynamic equilibrium (LTE), such as laser-driven hohlraums, stellar coronae, and supernova ejecta. The computational expense of traditional non-LTE (NLTE) opacity calculations can dominate the cost of simulating these systems and render NLTE physics untenable in most cases. In this paper, we present a new parallel implementation for NLTE opacity calculations that employs advanced load balancing, vectorization, and memory management to significantly reduce the computational costs of these calculations. Using the radiation-hydrodynamics code xRAGE , we demonstrate this new implementation on three test problems and present order-of-magnitude speed improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Doped GaInAs/GaP Quantum Well Superlattice Solar Cells With 27.5% Efficiency

Quantum wells can extend the absorption range of a solar cell and are typically placed in the intrinsic region of the device to enable efficient carrier collection via drift. Thick intrinsic regions are needed for significant absorption in the low bandgap quantum wells, resulting in a large depletion region recombination and ultimately a solar cell with a low fill factor (FF). However, in quantum well superlattice solar cells where tunneling plays a dominant role in carrier transport, collection by carrier diffusion may be possible, relieving the requirement for quantum wells to be placed in the intrinsic region. Here, we investigate doping in stress-balanced quantum well superlattice solar cells using thin 2 nm GaP barriers. Doping reduces J02 depletion region recombination and improves the solar cell FF up to 86.7%, but very high doping eventually reduces the carrier collection, leading to a tradeoff in efficiency. Finally, we show that the barrier thickness also plays an important role in carrier collection, and we demonstrate high efficiency 27.5% single-junction devices with a doped superlattice.

13 HYDRO ENERGY↗
DOE OSTI2024

The need for carbon-emissions-driven climate projections in CMIP7

Abstract. Previous phases of the Coupled Model Intercomparison Project (CMIP) have primarily focused on simulations driven by atmospheric concentrations of greenhouse gases (GHGs), for both idealized model experiments and climate projections of different emissions scenarios. We argue that although this approach was practical to allow parallel development of Earth system model simulations and detailed socioeconomic futures, carbon cycle uncertainty as represented by diverse, process-resolving Earth system models (ESMs) is not manifested in the scenario outcomes, thus omitting a dominant source of uncertainty in meeting the Paris Agreement. Mitigation policy is defined in terms of human activity (including emissions), with strategies varying in their timing of net-zero emissions, the balance of mitigation effort between short-lived and long-lived climate forcers, their reliance on land use strategy, and the extent and timing of carbon removals. To explore the response to these drivers, ESMs need to explicitly represent complete cycles of major GHGs, including natural processes and anthropogenic influences. Carbon removal and sequestration strategies, which rely on proposed human management of natural systems, are currently calculated in integrated assessment models (IAMs) during scenario development with only the net carbon emissions passed to the ESM. However, proper accounting of the coupled system impacts of and feedback on such interventions requires explicit process representation in ESMs to build self-consistent physical representations of their potential effectiveness and risks under climate change. We propose that CMIP7 efforts prioritize simulations driven by CO2 emissions from fossil fuel use and projected deployment of carbon dioxide removal technologies, as well as land use and management, using the process resolution allowed by state-of-the-art ESMs to resolve carbon–climate feedbacks. Post-CMIP7 ambitions should aim to incorporate modeling of non-CO2 GHGs (in particular, sources and sinks of methane and nitrous oxide) and process-based representation of carbon removal options. These developments will allow three primary benefits: (1) resources to be allocated to policy-relevant climate projections and better real-time information related to the detectability and verification of emissions reductions and their relationship to expected near-term climate impacts, (2) scenario modeling of the range of possible future climate states including Earth system processes and feedbacks that are increasingly well-represented in ESMs, and (3) optimal utilization of the strengths of ESMs in the wider context of climate modeling infrastructure (which includes simple climate models, machine learning approaches and kilometer-scale climate models).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

MULTIPHYSICS-MODELING OF FIRE-INDUCED URANIUM AEROSOL FORMATION – A-POSTERIORI BENCHMARKING OF EXPERIMENTS

The formulation and solution of a model that properly described the temperature profiles of uranium in fire conditions was prepared, verified and validated, to support a-posteriori benchmarking of historical experiments. Data from the multi-physics model combined with data regression from the experiments provides a useful tool for exploring sources of potential bias in historical experiments. The approaches used are augmented by visual observations and photographic evidence allowing for correlations between various phenomena, independent of whether or not that specific phenomena are part of the validated model. In this manner, the validated model is a tool for a-posteriori benchmarking, and not a model designed to replicate the experimental output exactly for any specific experiment. For experiments like the Clark (2015) experiments, this integrated model (including chemistry) demonstrated that the largest bias was associated in the metallurgy of the coupons. The model developed does serve as a basis for defending which of the historical alloys tested are the most appropriate for any other alloy being considered (e.g., one that has not yet been tested). The model identifies key reactions that impact uranium chemistry and can significantly bias experimental results. Through use of the model and energy balances, data from historical experiments by Elder and Tinkle (1980) were regressed and benchmarked. This study identifies that all uranium alloys are not equal and that there are distinct differences between the dominant alloys representing distinct metallurgical phases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Mining for Gluon Saturation at Colliders

Quantum chromodynamics (QCD) is the theory of strong interactions of quarks and gluons collectively called partons, the basic constituents of all nuclear matter. Its non-abelian character manifests in nature in the form of two remarkable properties: color confinement and asymptotic freedom. At high energies, perturbation theory can result in the growth and dominance of very gluon densities at small-x. If left uncontrolled, this growth can result in gluons eternally growing violating a number of mathematical bounds. The resolution to this problem lies by balancing gluon emissions by recombinating gluons at high energies: phenomena of gluon saturation. High energy nuclear and particle physics experiments have spent the past decades quantifying the structure of protons and nuclei in terms of their fundamental constituents confirming predicted extraordinary behavior of matter at extreme density and pressure conditions. In the process they have also measured seemingly unexpected phenomena. We will give a state of the art review of the underlying theoretical and experimental tools and measurements pertinent to gluon saturation physics. We will argue for the need of high energy electron-proton/ion colliders such as the proposed EIC (USA) and LHeC (Europe) to consolidate our knowledge of QCD knowledge in the small x kinematic domains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The advancing wave front on a sloping channel covered by a rod canopy following an instantaneous dam break

The drag coefficient Cd for a rigid and uniformly distributed rod canopy covering a sloping channel following the instantaneous collapse of a dam was examined using flume experiments. The measurements included space x and time t high resolution images of the water surface h(x, t) for multiple channel bed slopes So and water depths behind the dam Ho along with drag estimates provided by sequential load cells. Using these data, an analysis of the Saint-Venant equation (SVE) for the front speed was conducted using the diffusive wave approximation. An inferred Cd=0.4 from the h(x, t) data near the advancing front region, also confirmed by load cell measurements, is much reduced relative to its independently measured steady-uniform flow case. This finding suggests that drag reduction mechanisms associated with transients and flow disturbances are more likely to play a dominant role when compared to conventional sheltering or blocking effects on Cd examined in uniform flow. The increased air volume entrained into the advancing wave front region as determined from an inflow–outflow volume balance partly explains the Cd reduction from unity.

Mechanics↗
DOE OSTI2022

Checking nonflow assumptions and results via PHENIX published correlations in p+p, p+Au, d+Au, and 3 He+Au at √ s NN=200 GeV

Recently the PHENIX Collaboration has made available two-particle correlation Fourier coefficients for multiple detector combinations in minimum bias p+p and 0–5% central p+Au, d+Au, and 3 He+Au collisions at √ s NN = 200 GeV [Phys. Rev. C 105, 024901 (2022)]. Using these coefficients for three sets of two-particle correlations, azimuthal anisotropy coefficients v 2 and v 3 are extracted for midrapidity charged hadrons as a function of transverse momentum. In this paper, we use the available coefficients to explore various nonflow hypotheses as well as to compare the results with theoretical model calculations. The nonflow methods fail basic closure tests with ampt and pythia/angantyr, particularly when including correlations with particles in the low multiplicity light-projectile going direction. In data, the nonflow adjusted v 2 results are modestly lower in p+Au and the adjusted v 3 results are more significantly higher in p+Au and d+Au. However, the resulting higher values for the ratio v 3 /v 2 in p+Au at RHIC compared to p+Pb at the LHC is additional evidence for a significant overcorrection. Incorporating these additional checks, the conclusion that these flow coefficients are dominated by initial geometry coupled with final-state interactions (e.g., hydrodynamic expansion of quark-gluon plasma) remains true, and explanations based solely on initial-state glasma are ruled out. The detailed balance between intrinsic and fluctuation-driven geometry and the exact role of weakly versus strongly coupled prehydrodynamic evolution remains an open question for triangular flow, requiring further theoretical and experimental investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

MULTI-FIDELITY MODELING AND UNCERTAINTY QUANTIFICATION OF INVERTER BASED RESOURCES IN INTEGRATED T&D SYSTEMS

Uncertainty quantification plays a pivotal role in improving the accuracy and reliability of inverter operation within modern power systems that are increasingly dominated by inverter-based resources (IBRs). IBRs, especially those operating under grid forming (GFM) control, rely heavily on a complex set of control parameters and system measurements to maintain voltage, frequency, and power balance. Traditional deterministic modeling approaches often fail to capture these parameter deviations, potentially resulting in suboptimal control actions, reduced system stability, or even instability under high penetration of IBRs. In this paper, we demonstrate the application of model calibration and uncertainty quantification (UQ) principles to an integrated transmission and distribution (T&D) model involving a GFM converter and provide a framework for prioritizing control improvements, guiding robust design, and informing adaptive strategies that can accommodate real-time variability in system conditions. The proposed approach could be valuable in enhancing the robustness of current and future power systems under increased IBR penetrations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Role of Surface Hydroxyls and Lattice Oxygen in Governing Selectivity and Restructuring During Acetic Acid Conversion on Fe 3 O 4 (001)

Understanding the reactivity of carboxylic acids on metal oxide surfaces is critical for elucidating ketonization mechanisms relevant to biomass upgrading. Here, we investigate the adsorption and thermal decomposition of acetic acid (CH 3 COOH) on Fe 3 O 4 (001) using scanning tunneling microscopy (STM), temperature-programmed reaction spectroscopy (TPRS), and X-ray photoelectron spectroscopy (XPS). At room temperature, acetic acid adsorbs dissociatively to form ordered bidentate acetate (CH 3 COO) overlayer that lifts the (√2 × √2)R 45° surface reconstruction. TPRS reveals ketene (CH 2 CO) as the dominant product, along with CO, CO 2 , and H 2 O, the latter evolving via a Mars–van Krevelen (MvK) mechanism. Isotopic labeling shows preferential CO 2 formation from the carboxyl carbon and a more balanced CO/CO 2 ratio from the methyl carbon, suggesting distinct oxidation pathways. STM imaging reveals embedded acetate intermediates filling surface oxygen vacancies created in MvK steps. Upon product formation completion (~700 K), extensive surface etching is observed, with pits elongated along the octahedral Fe rows. Approximately 20% of the surface oxygen is removed, consistent with vacancy formation stoichiometry inferred from product distributions. These findings demonstrate that carboxylate-induced restructuring of Fe 3 O 4 (001) involves both surface healing and reduction processes, offering mechanistic insights relevant to ketonization and broader carboxylic acid chemistry on metal oxides.

acetic acid↗