SEARCH · Search NASA

Results for “Dynamic Consider Parameter”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Burton-Cabrera-Frank theory for the influence of off-cut direction on the growth of hcp crystals

The precise direction of the off-cut of a vicinal crystal surface determines the azimuth (in-plane angle) of the surface steps and thus their kink density. Here we develop a Burton-Cabrera-Frank (BCF) theory for step-flow crystal growth that includes the effects of step azimuth. We focus on the basal-plane surfaces of hexagonal close-packed (hcp) and related crystal structures, for which surface steps have alternating properties that depend on azimuth. We build on a previous treatment that considered only a single-step azimuth, 𝜙 = 0, that gives alternating steps of pure A and B type. As 𝜙 increases from zero, the density of kinks on the steps increases. For hcp-type crystals, in addition the structural difference between alternating steps decreases, until they become equivalent at 𝜙 = 30°. We consider the azimuth dependence of the step transparency, attachment, and repulsion parameters of the BCF model, including the effects of both geometrically required and thermally generated kinks. Example calculations are presented that demonstrate the overall effect of step azimuth on BCF predictions for step spacings during steady-state growth and for their dynamics during growth transients. Effects of symmetric and asymmetric step-step repulsion are discussed.

Ju, Guangxu↗

DOE OSTI2024

Impact of prescribed fire on soil microbial communities in a Southern Appalachian Forest clear-cut

Escalating wildfire frequency and severity, exacerbated by shifting climate patterns, pose significant ecological and economic challenges. Prescribed burns, a common forest management tool, aim to mitigate wildfire risks and protect biodiversity. Nevertheless, understanding the impact of prescribed burns on soil and microbial communities in temperate mixed forests, considering temporal dynamics and slash fuel types, remains crucial. Our study, conducted at the University of Tennessee Forest Resources AgResearch and Education Center in Oak Ridge, TN, employed controlled burns across various treatments, and the findings indicate that low-intensity prescribed burns have none or minimal short-term effects on soil parameters but may alter soil nutrient concentrations, as evidenced by significant changes in porewater acetate, formate, and nitrate concentrations. These burns also induce shifts in microbial community structure and diversity, with Proteobacteria and Acidobacteria increasing significantly post-fire, possibly aiding soil recovery. In contrast, Verrucomicrobia showed a notable decrease over time, and other specific microbial taxa correlated with soil pH, porewater nitrate, ammonium , and phosphate concentrations. Our research contributes to understanding the intricate relationships between prescribed fire, soil dynamics, and microbial responses in temperate mixed forests in the Southern Appalachian Region, which is valuable for informed land management practices in the face of evolving environmental challenges.

Microbiology↗

DOE OSTI2025

Catalytic resonance theory for parametric uncertainty of programmable catalysis

Microkinetic models are useful tools for screening catalytic materials; however, errors in their input parameters can lead to significant uncertainty in model predictions of catalyst performance. Here, in this work, we investigate the impact of linear scaling and Brønsted-Evans-Polanyi relation parametric uncertainty on microkinetic predictions of programmable-catalyst performance. Two case studies are considered: a generic A-to-B prototype reaction and the oxygen evolution reaction (OER). The results show that error-unaware models can accurately predict trends and, for the prototype reaction, values of optimal waveform parameters. The specific model parameters driving output uncertainty are identified via variance-based global sensitivity analysis. However, predictions of dynamic rate enhancement can decrease when uncertainty is propagated into the models. In both cases, we identify operating conditions where the programmable catalyst achieves a rate enhancement of at least one order of magnitude despite parametric uncertainty in the model, supporting programmable catalysis as a viable strategy for exceeding the Sabatier limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Models of SIV rebound after treatment interruption that involve multiple reactivation events

In order to assess the efficacy of novel HIV-1 treatments leading to a functional cure, the time to viral rebound is frequently used as a surrogate endpoint. The longer the time to viral rebound, the more efficacious the therapy. In support of such an approach, mathematical models serve as a connection between the size of the latent reservoir and the time to HIV-1 rebound after treatment interruption. The simplest of such models assumes that a single successful latent cell reactivation event leads to observable viremia after a period of exponential viral growth. Here we consider a generalization developed by Pinkevych et al. and Hill et al. of this simple model in which multiple reactivation events can occur, each contributing to the exponential growth of the viral load. We formalize and improve the previous derivation of the dynamics predicted by this model, and use the model to estimate relevant biological parameters from SIV rebound data. We confirm a previously described effect of very early antiretroviral therapy (ART) initiation on the rate of recrudescence and the viral load growth rate after treatment interruption. We find that every day ART initiation is delayed results in a 39% increase in the recrudescence rate (95% credible interval: [18%, 62%]), and a 11% decrease of the viral growth rate (95% credible interval: [4%, 20%]). We show that when viral rebound occurs early relative to the viral load doubling time, a model with multiple successful reactivation events fits the data better than a model with only a single successful reactivation event.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Pressure Exchanger for Energy Recovery in a Trans-Critical CO2 Refrigeration System

Trans-critical CO2 vapor compression (VC) refrigeration cycles require a high compression ratio, which is associated with high expansion losses. To recover these expansion losses, a pressure exchange process between the low- and high-pressure sides of the VC cycle is proposed and examined in this study. The proposed pressure exchange system is an open type constant volume process where the high- and low-pressure flows mix inside the system. This prototype is inspired by the pressure exchangers used in reverse-osmosis (RO) desalination systems. In this system, a 2D model was generated and modeled using the computational fluid dynamics (CFD) technique. The numerical model ignored any losses due to leakage or hydraulic friction and the process is considered adiabatic. For the modeling, it was assumed that the inlet conditions for the two pressure exchanger flows are similar to the flow conditions at the evaporator and gas cooler outlets in a VC cycle. Two parameters are examined to test the validity of the system and understand their effect on the performance, including the inlet flow rate represented by the inlet velocity and the process time represented by the speed of rotation. A total of nine cases were simulated and analyzed in this study.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2025

DESI 2024: Constraints on physics-focused aspects of dark energy using DESI DR1 BAO data

Baryon acoustic oscillation data from the first year of the Dark Energy Spectroscopic Instrument (DESI) provide near percent-level precision of cosmic distances in seven bins over the redshift range z=0.1–4.2. Here, this paper is the follow-up to the original DESI BAO cosmology paper [A. G. Adame et al. (DESI Collaboration), arXiv:2404.03002], which considered the conventional w 0 w a cold dark matter (CDM) model. We use the novel DESI data, together with other cosmic probes, to constrain the background expansion history using some well-motivated physical classes of dark energy. In particular, we explore three physics-focused behaviors of dark energy from the equation of state and energy density perspectives: the thawing class (matching many simple quintessence potentials), emergent class (where dark energy comes into being recently, as in phase transition models), and mirage class [where phenomenologically the distance to cosmic microwave background (CMB) last scattering is close to that from a cosmological constant Λ despite dark energy dynamics]. All three classes fit the data at least as well as Λ ⁢CDM, and indeed can improve on it by Δ⁢χ 2 ≈ –5 to –17 for the combination of DESI BAO with CMB and supernova data while having one more parameter. The mirage class does essentially as well as w 0 ⁢w a CDM and exhibits moderate to strong Bayesian evidence preference with respect to Λ⁢ CDM. These classes of dynamical behaviors highlight worthwhile avenues for further exploration into the nature of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Mapping dynamical ejecta and disk masses from numerical relativity simulations of neutron star mergers

Abstract We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions, that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear EOS lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

DOE OSTI2025

Resistive Wall Simulations for the DARHT Multi-Pulse Test Line

The multi-pulse test line (MPTL) presently under T development at the Dual Axis Radiography Hydrodynamic Test (DARHT) facility may have several meters of drifting beam transport without external magnetic focusing field. The absence of external magnetic focusing is the worst case for growth of the resistive-wall instability. This instability is usually considered to be strictly a long-pulse (low-frequency) problem, but has been theoretically shown to couple the multiple pulses in a long pulse train, such as are expected to be the subject of experimentation on the MPTL. In this note, we explore the resistive wall instability in the MPTL parameter range using our LAMDA beam dynamics code. (LAMDA stands for Linear Accelerator Model for DARHT)

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Role of Intra-Domain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: An Approach for Spatially Resolving Dynamics and Ion Transport

The design of safe and high-performance, nanostructured, block polymer (BP) electrolytes for lithium-ion batteries requires a thorough understanding of the key parameters that govern local structure and dynamics. Yet, the interfaces between microphase separated domains can introduce complexities in this local behavior that can be challenging to quantify. Herein, the local polymer, cation (Li + ), and anion dynamics were described in salt-doped polystyrene block poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) through a quantitative framework that considered the effects of polymer architecture, segmental mixing, chain stretching, and confinement on polymer mobility and ion transport. This framework was validated through nuclear magnetic resonance (NMR) spectroscopy measurements on solid (dry) polymer electrolyte samples. Notably, a mobility transition temperature (T mobility ) was identified through NMR spectroscopy that captured the local dynamics more accurately than the thermal glass transition temperature. Additionally, the approach quantitatively described the mobility gradient across a domain when segmental mixing effects were combined with chain stretching and confinement information, especially at higher segregation strengths – facilitating the assessment of local ion diffusion and conductivity. As a result, spatially averaged local ion diffusion predictions quantitatively matched NMR-measured ion diffusivities in the BP samples, while spatially summed ionic conductivity predictions across a domain qualitatively captured trends in the measured ionic conductivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2024

KiDS-1000 cosmology: Combined second- and third-order shear statistics

Aims.In this work, we perform the first cosmological parameter analysis of the fourth release of Kilo Degree Survey (KiDS-1000) data with second- and third-order shear statistics. This paper builds on a series of studies aimed at describing the roadmap to third-order shear statistics. Methods.We derived and tested a combined model of the second-order shear statistic, namely, the COSEBIs and the third-order aperture mass statistics 〈ℳ ap 3 〉 in a tomographic set-up. We validated our pipeline withN-body mock simulations of the KiDS-1000 data release. To model the second- and third-order statistics, we used the latest version of HMCODE2020 for the power spectrum and BIHALOFITfor the bispectrum. Furthermore, we used an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering only equal smoothing radii for the 〈ℳ ap 3 〉 part of the data vector. This makes it possible to carry out a data analysis of the KiDS-1000 data release using a combined analysis of COSEBIs and third-order shear statistics. Results.We first validated the accuracy of our modelling by analysing a noise-free mock data vector, assuming the KiDS-1000 error budget, finding a shift in the maximum of the posterior distribution of the matter density parameter, ΔΩ m < 0.02 σ Ω m , and of the structure growth parameter, ΔS 8 < 0.05 σ S 8 . Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained Ω m = 0.248 −0.055 +0.062 andS 8 = σ 8 √(Ω m /0.3 )= 0.772 ± 0.022. The geometric average on the errors of Ω m andS 8 of the combined statistics decreases, compared to the second-order statistic, by a factor of 2.2.

Astronomy & Astrophysics↗

DOE OSTI2024

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al . J. Am. Chem. Soc. 142, 10459–10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. Here, in this work, we explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck–Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

The Ecosystem as Super-Organ/ism, Revisited: Scaling Hydraulics to Forests under Climate Change

Synopsis Classic debates in community ecology focused on the complexities of considering an ecosystem as a super-organ or organism. New consideration of such perspectives could clarify mechanisms underlying the dynamics of forest carbon dioxide (CO2) uptake and water vapor loss, important for predicting and managing the future of Earth’s ecosystems and climate system. Here, we provide a rubric for considering ecosystem traits as aggregated, systemic, or emergent, i.e., representing the ecosystem as an aggregate of its individuals or as a metaphorical or literal super-organ or organism. We review recent approaches to scaling-up plant water relations (hydraulics) concepts developed for organs and organisms to enable and interpret measurements at ecosystem-level. We focus on three community-scale versions of water relations traits that have potential to provide mechanistic insight into climate change responses of forest CO2 and H2O gas exchange and productivity: leaf water potential (Ψcanopy), pressure volume curves (eco-PV), and hydraulic conductance (Keco). These analyses can reveal additional ecosystem-scale parameters analogous to those typically quantified for leaves or plants (e.g., wilting point and hydraulic vulnerability) that may act as thresholds in forest responses to drought, including growth cessation, mortality, and flammability. We unite these concepts in a novel framework to predict Ψcanopy and its approaching of critical thresholds during drought, using measurements of Keco and eco-PV curves. We thus delineate how the extension of water relations concepts from organ- and organism-scales can reveal the hydraulic constraints on the interaction of vegetation and climate and provide new mechanistic understanding and prediction of forest water use and productivity.

Zoology↗