Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium constant”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Probing the thermodynamics and kinetics of ethylene carbonate reduction at the electrode–electrolyte interface with molecular simulations

Understanding the formation of the solid-electrolyte interphase (SEI) in lithium-ion batteries (LIBs) is an ongoing area of research due to its high degree of complexity and the difficulties encountered by experimental studies. Herein, we investigate the initial stage of SEI growth, the reduction reaction of ethylene carbonate (EC), from both a thermodynamic and kinetic approach with theory and molecular simulation. We employed both the potential distribution theorem (PDT) and the SMD implicit solvent model to EC solvation for the estimation of reduction potentials of Li$^+$, EC, and Li$^+$-solvating EC (s-EC), as well as reduction rate constants of EC and s-EC. We find that solvation effects greatly influence these quantities of interest, particularly the Li$^+$/Li reference electrode potential in EC solvent. Further, we also compute the inner- and outer-sphere reorganization energies for both EC and s-EC at the interface of liquid EC and a hydroxyl-terminated graphite surface, where total reorganization energies are predicted to be 76.6 and 88.9 kcal/mol, respectively. With the computed reorganization energies, we estimate reduction rate constants across a range of overpotentials and show that EC has a larger electron transfer rate constant than s-EC at equilibrium, despite s-EC being more thermodynamically favorable. Furthermore, this manuscript demonstrates how ion solvation effects largely govern the prediction of reduction potentials and electron transfer rate constants at the electrode-electrolyte interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrochemical Study of OH - in Molten MgCl 2 -NaCl-KCl Salts

Here the present study determines the equilibrium concentration of OH- in molten MgCl 2 -NaCl-KCl salts and measures the kinetic properties data of OH- (diffusion coefficient, exchange current density, charge transfer coefficient, limiting current density, and standard rate constant) accordingly at different temperatures. The equilibrium concentration of OH- descends as the temperature increases, which ranges from 4.11 to 2.84 × 10 -6 mol cm -3 in the temperature range of 600 °C–800 °C. The diffusion coefficient follows Arrhenius law with a value between 2.46 and 4.48 × 10 -4 cm 2 s -1 . The exchange and limit current densities, and standard rate constant increase with temperature. The calculated standard rate constant κ 0 follows Arrhenius law, which is in the range from 1.23 to 2.50 × 10 -4 cm s -1 . The measured charge transfer coefficient is independent of temperature, which validates the present measurements.

14 SOLAR ENERGY↗

Establishing structure-sensitivity of ceria reducibility: real-time observations of surface-hydrogen interactions

The first layer of atoms on an oxide catalyst provides the first sites for adsorption of reactants and the last sites before products or oxygen are desorbed. We employ a unique combination of morphological, structural, and chemical analyses of a model ceria catalyst with different surface terminations under an H 2 environment to unequivocally establish the effect of the last layer of atoms on surface reduction. (111) and (100) terminated epitaxial islands of ceria are simultaneously studied in situ allowing for a direct investigation of the structure–reducibility relationship under identical conditions. Kinetic rate constants of Ce 4+ to Ce 3+ transformation and equilibrium concentrations are extracted for both surface terminations. Unlike the kinetic rate constants, which are practically the same for both types of islands, more pronounced oxygen release, and overall higher reducibility were observed for (100) islands compared to (111) ones. The findings are in agreement with coordination-limited oxygen vacancy formation energies calculated by density functional theory. The results point out the important aspect of surface terminations in redox processes, with particular impact on the catalytic reactions of a variety of catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic, magnetic, and optical properties of bulk and (1 1 1)-surfaces of CoMnZnSi quaternary Heusler alloy

The structural stability, electronic structure, and magnetic and optical properties of CoMnZnSi quaternary Heusler alloy (QHA) in the bulk and (1 1 1)-slab forms have been investigated by performing density functional theory (DFT) calculations. Formation and cohesive energies, and elastic constants confirm that the bulk CoMnZnSi is chemically and mechanically stable at an equilibrium lattice parameter of 5.81 Å. The calculated elastic constants also indicate that this QHA has ductile and anisotropic features. We investigate and discuss the bonding behavior from charge density distribution and density of states. Furthermore, these calculations in the bulk phase show a perfect half-metallic behavior with an integer value of magnetic moment (4 µB) and a large spin-flip gap of 0.39 eV. On the other hand, the Co (1 1 1) and Si (1 1 1)-slabs exhibit semi-metallic nature at the Fermi level. These results indicate that CoMnZnSi qualifies for spintronic applications. In the bulk and Co (1 1 1), Zn (1 1 1), and Si (1 1 1)-surfaces, the magnetic moment of Co and Mn atoms are ferromagnetically aligned, while the ferrimagnetic alignment of Co with Mn has been found in Mn (1 1 1)-surface. We also investigate the optimal conductivity, dielectric functions, reflectivity, absorption coefficient, refractive index, and loss function to understand the underlying optical properties of this alloy.

36 MATERIALS SCIENCE↗

Evaluating Novel Polymer-Composite Resins for the Ability to Derive Elastic Properties for Explicit Formulations, or to Develop Nonlinear Models Thereof

This project evaluates the ability to derive elastic material properties for novel additively manufactured carbon-polymer matrix materials. Initially stress-strain equilibrium equations were thought to be sufficient to characterize the elastic region of material performance. However, initial test data shows the cured material behaves nonlinearly in the elastic region. To properly determine the elastic performance, additional strain data must be collected in future testing. This paper discusses the initial test data, methods and models used to describe nonlinear elasticity, and the determination of mean elastic material properties, such as Lame’s constants, using the stress-strain equilibrium equations.

36 MATERIALS SCIENCE↗

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assemblies of Polyacrylonitrile-Derived Photoactive Polymers as Blue and Green Light Photo-Cocatalysts for Cu-Catalyzed ATRP in Water and Organic Solvents

Photoluminescent nanosized quasi-spherical polymeric assemblies prepared by the hydrothermal reaction of polyacrylonitrile (PAN), ht-PLP PAN , were demonstrated to have the ability to photo-induce atom transfer radical polymerization (ATRP) catalyzed by low, parts per million concentrations of Cu II complex with tris(2-pyridylmethyl)amine (TPMA). Such photo induced ATRP reactions of acrylate and methacrylate monomers were performed in water or organic solvents, using ht-PLP PAN as the photo-cocatalyst under blue or green light irradiation. Mechanistic studies indicate that ht-PLP PAN helps to sustain the polymerization by facilitating the activation of alkyl bromide species by two modes: 1) green or blue light-driven photoreduction of the Cu II catalyst to the activating Cu I form, and 2) direct activation of dormant alkyl bromide species which occurs only under blue light. The photoreduction of the Cu II complex by ht-PLP PAN was confirmed by linear sweep voltammetry performed under illumination. Analysis of the polymerization kinetics in aqueous media indicated even though Cu I complexes comprised only 1–1.4% of all Cu species at equilibrium, they exhibited high activation rate constant and activated the alkyl bromide initiators five to six orders of magnitude faster than ht-PLP PAN .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗

Evidence for Photoinduced Polaron Generation in a High Persistence Length Low Bandgap Conjugated Polymer in Solution

Using ultrafast time-resolved infrared (TRIR) spectroscopy, we studied the solution-phase excited-state structural evolution of an indacenodithiophene-co-benzothiadiazole polymer (C8-IDTBT). Following band gap excitation, the TRIR spectra reveal vibrational features that develop within 10 ps and decay over 4 ns. Using pulse radiolysis measurements, charge-modulation spectroscopy, and quantum-chemical calculations, the IR features are assigned to polaron pairs. Interestingly, these features appear on an evolving broad mid-IR electronic absorption background, with kinetics correlating with the formation and decay of the cationradical vibrational bands. A three-state kinetic model successfully reproduces the spectral evolution, revealing that the polaron and exciton populations exist in dynamic equilibrium on picosecond time scales, with time constants for exciton dissociation in the range of 3− 5 ps and polaron-to-exciton reformation between 20 and 100 ps, while both species decay to the ground state on much slower nanosecond time scales (∼1 ns), yielding a remarkably high polaron-generation efficiency, higher than 50%. These findings provide fundamental insights into intramolecular charge photogeneration mechanisms in conjugated polymers, demonstrating efficient bound-polaron formation in isolated polymer chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗

A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)

Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solubility Product of Vanadinite Pb 5 (VO 4 ) 3 Cl at 25 °C—A Comprehensive Approach to Incongruent Dissolution Modeling

Although vanadinite (Pb 5 (VO 4 ) 3 Cl occurs in abundance in various terrestrial geochemical systems of natural and anthropogenic origin and is seriously considered as a potential nuclear waste sequestering agent, its actual application is severely limited by a lack of understanding of its basic thermodynamic parameters. In this regard, the greatest challenge is posed by its incongruent dissolution, which is a pivotal hurdle for effective geochemical modeling. Our paper presents an universal approach for geochemical computing of systems undergoing incongruent dissolution which, along with unique, long-term experiments on vanadinites’ stability, allowed us to determine the mineral solubility constant. The dissolution experiments were carried out at pH = 3.5 for 12 years. Vanadinite has dissolved incongruently, continuously re-precipitating into chervetite (Pb 2 V 2 O 7 ) with the two minerals remaining in mutual equilibrium until termination of the experiments. The empirically derived solubility constant K sp,V,298 = 10 –91.89 ± 0.05 of vanadinite was determined for the first time. The proposed modeling method is versatile and can be adopted to other mineral systems undergoing incongruent dissolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗

Reversible Ligand Protonation in Noninnocent Constrained-Geometry-Like Group 4 Complexes

The proligands C 5 H 5 CMe 2 CHPhOx R (R = Me 2 , CHMe 2 , and CMe 3 ) react with M(NMe 2 ) 4 (M = Ti, Zr, and Hf) to give monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 , doubly deprotonated (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 , or a mixture of both. The observed products depend on reaction conditions, the oxazoline substituent, and the metal center, with 4-t-butyl-oxazoline or Ti giving only the monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 . Amine elimination from (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 to give (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 is reversible for 4,4-dimethyl- and 4-isopropyl-oxazoline-based ligands in Zr or Hf complexes. Temperature-dependent kinetic studies of the equilibration of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 , HNMe 2 and (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 provide the experimental thermodynamic parameters ΔS° = -17.4 ± 2.6 cal·mol –1 K –1 and ΔH° = -6.8 ± 0.8 kcal·mol –1 . An Eyring plot of the rate constants, determined from the system as it approaches equilibrium, gives the activation entropy and activation enthalpy for the addition of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 and HNMe 2 as -36 ± 4 cal·mol –1 K –1 and 8.9 ± 1.3 kcal·mol –1 , respectively; the elimination of HNMe 2 from (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 is characterized by ΔS ‡ = -19 ± 5 cal·mol –1 K –1 and ΔH ‡ = 15.7 ± 1.5 kcal·mol –1 . DFT computational models indicate a single-step, nonlinear transfer of the H between the benzylic position of the noninnocent, oxazoline-coordinated ligand and NMe 2 . Computations also confirm the negative activation entropy and the trends in the barriers support the experimental results. Together, these studies indicate the importance of steric effects from the oxazoline ligand, metal center, ancillary ligands, and leaving group on the shuttling of the proton between HNMe 2 and the noninnocent ligand. Finally, these effects suggest that coordination of oxazoline to the metal center is a key part of the benzylic deprotonation and noninnocent behavior of the cyclopentadienyl-oxazoline ligand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation of the oxyl’s potential energy surface by the spectral kinetics of a vibrational mode

One of the most reactive intermediates for oxidative reactions is the oxyl radical, an electron-deficient oxygen atom. The discovery of a new vibration upon photoexcitation of the oxygen evolution catalysis detected the oxyl radical at the SrTiO3 surface. The vibration was assigned to a motion of the sub-surface oxygen underneath the titanium oxyl (Ti–O●−) created upon hole transfer to (or electron extraction from) a hydroxylated surface site. Evidence for such an interfacial mode is derived from its spectral shape, which exhibited a Fano resonance—a coupling of a sharp normal mode to continuum excitations. Here, this Fano resonance is utilized to derive precise formation kinetics of the oxyl radical and its associated potential energy surface (PES). From the Fano lineshape, the formation kinetics are obtained from the anti-resonance (the kinetics of the coupling factor), the resonance (the kinetics of the coupled continuum excitations), and the frequency integrated spectrum (the kinetics of the normal mode’s cross-section). All three perspectives yield logistic function growth with a half-rise of 2.3 ± 0.3 ps and a time constant of 0.48 ± 0.09 ps. A non-equilibrium transient associated with photoexcitation is separated from the rise of the equilibrated PES. The logistic function characterizes the oxyl coverage at the very initial stages (t ∼ 0) to have an exponential growth rate that quickly decreases toward zero as a limiting coverage is reached. Such time-dependent reaction kinetics identify a dynamic activation barrier associated with the formation of a PES and quantify it for oxyl radical coverage.

Chemistry↗

Mirror at the edge of the universe: Reflections on an accelerated boundary correspondence with de Sitter cosmology

An accelerated boundary correspondence (ABC) is solved for the de Sitter moving mirror cosmology. The beta Bogoliubov coefficients reveal the particle spectrum is a Planck distribution with temperature inversely proportional to horizon radius. The quantum stress tensor indicates a constant emission of energy flux consistent with eternal equilibrium, while the total energy carried by the particles remains finite. In conclusion, the curved spacetime transformation to flat spacetime with an accelerated boundary is illustrated, and also shown for anti-de Sitter (AdS) spacetime.

79 ASTRONOMY AND ASTROPHYSICS↗

Manuscript model outputs, model source code, and figure scripts: The role of geomorphology in mediating biomass allocation impacts on salt-marsh resilience and carbon accumulation

This data package provides model source code (C++), model outputs, and figure generation scripts (R) needed to reproduce the following manuscript: Bruns, Nicholas E., Genevieve L. Noyce, and Matthew L. Kirwan. "The Role of Geomorphology in Mediating Biomass Allocation Impacts on Salt-Marsh Resilience and Carbon Accumulation." Estuarine, Coastal and Shelf Science 327 (December 2025): 109549. https://doi.org/10.1016/j.ecss.2025.109549. This manuscript investigates how geomorphology mediates the impact of biomass allocation shifts on salt marsh persistence and carbon (C) sequestration under sea level rise. We use a 1-D soil-column model (Kirwan and Mudd 2012) to perform experiments across a range of static root:shoot ratios (RSR = 1-4) spanning observed values. The model explicitly simulates interactions between tidal inundation, productivity, inorganic sediment deposition, and organic matter accumulation. Experiments determine whether geomorphic feedbacks amplify, dampen, or leave unchanged the ecosystem response to biomass allocation shifts. A first experiment uses constant sea level rise (2.5 mm/yr) to examine equilibrium dynamics and their influence on carbon accumulation rates across different suspended sediment concentrations (SSC). A second set of experiments calculates threshold sea level rise rates for marsh drowning across RSR and SSC combinations. Final experiments apply accelerating sea level rise scenarios (2000-2200) derived from NOAA projections (Sweet et al. 2022) to generate an envelope of expected responses, quantifying the importance of biomass allocation shifts on marsh carbon accumulation and persistence. All experiments are repeated across SSC ranging from 5-50 mg/L to investigate how these interactions vary in micro-tidal marshes with different sediment supplies. Package contents: * README.txt with detailed description of package contents and instructions for reproducing manuscript figures and model outputs * R scripts for generating all manuscript figures * C++ baseline model code from Kirwan and Mudd (2012) * C++ source code for 4 experimental model variants used in the manuscript, extending above baseline code * Model input files (.csv, .txt) including sea level rise scenarios * Model output files (.txt) used in manuscript analyses Temporal coverage: Model simulations span years 1900-2200, with accelerating sea level rise scenarios for 2000-2200. Key variables: Root:shoot ratio, suspended sediment concentration, carbon accumulation rate, vertical accretion rate, marsh elevation, inundation depth, threshold sea level rise rate.

54 ENVIRONMENTAL SCIENCES↗