Search NASA⌕ Search

SEARCH · Search NASA

Results for “activity coefficient”

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

Developing predictive models for µ opioid receptor binding using machine learning and deep learning techniques

Opioids exert their analgesic effect by binding to the µ opioid receptor (MOR), which initiates a downstream signaling pathway, eventually inhibiting pain transmission in the spinal cord. However, current opioids are addictive, often leading to overdose contributing to the opioid crisis in the United States. Therefore, understanding the structure-activity relationship between MOR and its ligands is essential for predicting MOR binding of chemicals, which could assist in the development of non-addictive or less-addictive opioid analgesics. This study aimed to develop machine learning and deep learning models for predicting MOR binding activity of chemicals. Chemicals with MOR binding activity data were first curated from public databases and the literature. Molecular descriptors of the curated chemicals were calculated using software Mold2. The chemicals were then split into training and external validation datasets. Random forest, k-nearest neighbors, support vector machine, multi-layer perceptron, and long short-term memory models were developed and evaluated using 5-fold cross-validations and external validations, resulting in Matthews correlation coefficients of 0.528–0.654 and 0.408, respectively. Furthermore, prediction confidence and applicability domain analyses highlighted their importance to the models’ applicability. Our results suggest that the developed models could be useful for identifying MOR binders, potentially aiding in the development of non-addictive or less-addictive drugs targeting MOR.

Research & Experimental Medicine↗

Simplified calculation of the area specific impedance for solid-state battery design

Simplified algebraic area specific impedance (ASI) correlations have been developed for solid-state composite battery electrodes made of a single ion conducting electrolyte, conductive additive, and intercalation active material. Two ASI expressions were developed, one for short times ( i.e. , pulsed power operation) and another for the pseudo steady state operation (i.e., sustained discharge for energy estimation). A full electrochemical model based on porous electrode theory was developed to examine the accuracy of the simplified ASI expressions. The simplified expressions agree favorably with full model results over a wide range of parameters (i.e. , electrode thicknesses, electrolyte conductivities, solid-state diffusion coefficients, specific surface areas, etc.) and conditions (i.e. , C-rates, states of charge, and pulse times). Under most conditions, the error between the full model and the correlations is well below 7 %. Higher errors were observed for the pseudo steady state expression at high/ low states of charge where the assumption of uniform reaction distributions loses validity. Here, the short time ASI has higher error at low states of charge due to the nonlinearity of the open circuit voltage equation, which is assumed linear in the formulation of the simplified algebraic expression.

25 ENERGY STORAGE↗

Diffusion coefficients predicting facet-dependent crystallization in amorphous silicon nitride

Amorphous silicon nitride is a common material in microelectronics devices, which acts as an insulating barrier. Extended annealing times at elevated temperature can initiate crystallization of α-Si 3 ⁢ N 4 , which does not possess the same barrier properties. Molecular dynamics can resolve the fundamental mechanism for α-Si 3 ⁢N 4 crystallization and the influence of local environments. Here, we compare two interatomic potentials and conclude that these models predict structural features (e.g., angular distributions and densities) which span the range of experimental measurements. We confirmed these models reproduce experimental estimates of activation energy and leveraged these models to identify crystallization drivers. We conclude that near-T g , facet-dependent silicon nitride crystal growth rates can be predicted directly by either bulk or interfacial diffusion properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CTRL-STEER: Closed-Loop Neuron Activation Control in Vision-Language-Action Models

Vision-Language-Action (VLA) models enable test-time behavioral steering via neuron-level interventions, but existing methods use fixed strengths and operate in open loop. This static modulation fails under evolving task dynamics, leading to overcorrection, oscillations, and reduced task success—especially for temporal attributes like speed. We propose CTRL-STEER, a control-theoretic framework that casts activation steering as closed-loop feedback with adaptive, time-varying interventions. Instead of assuming neurons encode temporal concepts, we steer along motion-aligned residual directions and regulate intervention magnitude via feedback. We instantiate this with both PID and reinforcement learning controllers that jointly optimize concept adherence and task success. Experiments on fine-tuned OpenVLA policies across four LIBERO suites show improved stability and a better steering–success trade-off over fixed-coefficient baselines, without retraining the base model.

Babu, Abhijith [Florida International University, ↗

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗

Fast Growth of Ga2O3 on Highly Offcut (100)-Oriented Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. More importantly, the impact ionization coefficients along the [100] direction are low, leading to the highest possible critical fields. This is advantageous compared to the current state of the art, (001), due to reduced surface defects and increased possible breakdown voltage. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption, and on-axis (100) growth easily forms twin domains. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction. These films show step-flow growth from (20-1) step-edges and high electron mobility due to suppressed twins. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (>5 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3 and critical breakdown field in Schottky barrier diodes comparable with state-of-the-art (001) Ga2O3 without edge termination.

36 MATERIALS SCIENCE↗

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗

Day-to-day reliability of basal heart rate and short-term and ultra short-term heart rate variability assessment by the Equivital eq02+ LifeMonitor in US Army soldiers

Introduction The present study determined the (1) day-to-day reliability of basal heart rate (HR) and HR variability (HRV) measured by the Equivital eq02+ LifeMonitor and (2) agreement of ultra short-term HRV compared with short-term HRV. Methods Twenty-three active-duty US Army Soldiers (5 females, 18 males) completed two experimental visits separated by >48 hours with restrictions consistent with basal monitoring (eg, exercise, dietary), with measurements after supine rest at minutes 20–21 (ultra short-term) and minutes 20–25 (short-term). HRV was assessed as the SD of R–R intervals (SDNN) and the square root of the mean squared differences between consecutive R–R intervals (RMSSD). Results The day-to-day reliability (intraclass correlation coefficient (ICC)) using linear-mixed model approach was good for HR (0.849, 95% CI: 0.689 to 0.933) and RMSSD (ICC: 0.823, 95% CI: 0.623 to 0.920). SDNN had moderate day-to-day reliability with greater variation (ICC: 0.689, 95% CI: 0.428 to 0.858). The reliability of RMSSD was slightly improved when considering the effect of respiration (ICC: 0.821, 95% CI: 0.672 to 0.944). There was no bias for HR measured for 1 min versus 5 min (p=0.511). For 1 min measurements versus 5 min, there was a very modest mean bias of −4 ms for SDNN and −1 ms for RMSSD (p≤0.023). Conclusion When preceded by a 20 min stabilisation period using restrictions consistent with basal monitoring and measuring respiration, military personnel can rely on the eq02+ for basal HR and RMSSD monitoring but should be more cautious using SDNN. These data also support using ultra short-term measurements when following these procedures.

General & Internal Medicine↗

A critical review of electrochemical heat pump technologies: Status, challenges, and perspectives

The development of advanced heat pump technologies is critical for reducing global energy consumption in the building sector, where space heating and cooling account for nearly 50% of energy use. Electrochemical heat pumps (EHPs) offer a promising alternative to vapor compression systems by enabling direct electrochemical-to-thermal energy conversion, often with environmentally benign working fluids that exhibit low or zero global warming potential (GWP). Prior literature has predominantly focused on chemically reactive heat pumps, while comprehensive assessments of electrochemical mechanisms remain limited. Here, this review addresses this gap by systematically evaluating the underlying principles, architectures, and performance metrics of EHP systems. Compared to conventional vapor compression systems, EHPs can achieve 10%-30% higher energy efficiency, with reported cooling coefficients of performance (COP c ) ranging from 3.5 to 14.3 under standard operating conditions. Despite these advantages, widespread adoption is hindered by challenges including membrane degradation, electrode fouling, sluggish redox kinetics, and elevated system-level capital costs. To address these limitations, the review outlines three research priorities: (i) the development of advanced membranes, catalysts, and electrode materials with enhanced chemical and mechanical stability; (ii) the application of molecular-level simulations for the rational design of high-performance redox-active working fluids; and (iii) the integration of advanced diagnostic techniques for real-time monitoring and sustained operation of EHPs. By consolidating recent advances and explicitly identifying technological and scientific gaps, this work uniquely contributes a comprehensive framework for guiding future electrochemical heat pump research and facilitating the transition to sustainable thermal management technologies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Cs X Si 15 P 21 ( X = Sn or Pb): Polar Noncentrosymmetric Si–P Frameworks Stabilized by Covalent X –P Bonding

Metal silicon phosphides composed of earth-abundant Si and P tend to exhibit semiconducting properties and adopt diverse crystal structures with relatively small additions of structure-directing elements. The potential of silicon phosphide materials in nonlinear optical applications has been hindered by the inability to systematically produce noncentrosymmetric structures with such a flexible framework. Here, in this work, two isostructural compounds with a novel noncentrosymmetric structure were made possible by the inclusion of elements with stereochemically active lone pairs (Sn 2+ and Pb 2+ ). The structures were determined through single-crystal and synchrotron powder X-ray diffraction. Analysis of chemical bonding in real space through the electron localization function revealed stereochemically active Pb 2+ and Sn 2+ species in a trigonal pyramidal coordination with {Pb/Sn}–P bonds. Such covalent bonding between Pb and P is quite uncommon in extended solids and has been reported in a few rare instances. Band structure calculations and linear optical measurements confirm the semiconducting nature of Cs X Si 15 P 21 ( X = Sn or Pb). The synthesis was optimized to yield high-purity polycrystalline samples. The nonlinear optical properties show promising second-harmonic generation (SHG) coefficients from the Kurtz–Perry method. First-principles calculations of the nonlinear optical properties support the experimentally determined SHG values and provide moderate values of birefringence, suggesting Cs X Si 15 P 21 could be phase-matchable and practical nonlinear optical materials in the mid-IR region.

crystal structure↗

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

Design of an alumina forming coating for Nb-base refractory alloys

Refractory multi-principal element alloys (RMPEAs) promise to significantly enhance gas turbine engine efficiency, but their poor oxidation performance inhibits their implementation. Alumina-forming bond coat alloys can provide oxidation protection, but discovering suitable chemistries remains a challenge. We employed a design methodology that screens for alumina-formation capability using Al activity and phase constitution predictions from CalPhaD (Thermo-Calc). Alloy down-selection from approximately 7,800 alloys in the Nb-Si-Ti-Al-Hf system was conducted via analysis of calculated thermodynamic properties with number-density topology style maps. This approach is validated by creating and testing the composition Nb 12 Si 23 Ti 24 Al 36 Hf 5 , which forms protective alumina scales up to 1400 °C and resists pesting at 800 °C. Further, the alloy has an average coefficient of thermal expansion of ~10.1 ppm/K, making it well matched to Nb-based refractory alloys. The methodology will be useful for the design of coatings for RMPEAs, enabling their implementation and significant efficiency benefits sooner.

36 MATERIALS SCIENCE↗

Hanford Tank Waste Matrix Impact on Ion Exchange Performance Using Crystalline Silicotitanate

The removal of radiocesium from Hanford tank waste supernate is a critical step in preparing feed for low-activity waste immobilization. This study evaluated cesium ion exchange performance using crystalline silicotitanate (CST) media in a series of tests designed to evaluate the influence of waste matrix variability on capacity and kinetics. Tank waste supernate subsampled from five Hanford double-shell tanks encompassed a range of sodium, hydroxide, nitrate, and nitrite concentrations in order to assess the impact of feed variability on the performance of the ion exchange system. Both equilibrium and dynamic ion exchange tests were conducted to quantify cesium distribution coefficients and breakthrough behavior under prototypic operating conditions. Results indicated that effective cesium capacity varied by up to a factor of five across the matrices tested, with higher sodium concentrations significantly reducing uptake. Kinetic behavior was similarly matrix-dependent, with solution viscosity contributing to a twofold variation in mass-transfer rates. These results demonstrate the strong dependence of CST ion exchange performance on waste composition and must be incorporated into predictive models for future treatment system design and optimization.

Westesen, Amy M.↗

Structure and Dynamics of Water Confined in Transition Metal Carbide MXenes: Implications for Electrochemical Applications

Two-dimensional transition metal carbides and nitrides (MXenes) are an important family of electrochemically active 2D materials. MXenes combine high conductivity with hydrophilicity, making them attractive materials for many applications, including electrochemical energy storage, sensing, desalination, and others. In order to better understand the role of structure on MXene properties, here, we investigated the vibrational properties and diffusion of water in MXenes with differing layer thicknesses and transition metal compositions using inelastic, quasi-elastic, and small-angle neutron scattering. We found that all of the Mo-containing MXenes studied here exhibited comparable vibrational dynamics and diffusion coefficients to each other and to previously studied Ti 3 C 2 T x . However, Ti 2 CT x was distinguished by its faster diffusion and more hydroxyl groups compared to the other MXenes studied. These results can help guide the selection of appropriate MXenes for energy storage and electrochemical water purification applications.

MXenes↗

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗

Mechanism, Thermochemistry, and Kinetics for the CH + N 2 Reaction Leading to Prompt NO Formation in Combustion

Here, the reaction of CH + N 2 forming H + NCN is a remarkable example of activation of the nitrogen triple bond and is an important source of prompt NO in combustion. The reaction pathway is complex and proceeds through two competing mechanisms: a cyclic addition channel initiated by c-HC(NN) and a chain-addition channel initiated by HCNN, both of which eventually form HNCN prior to dissociation to H + NCN. This work reinvestigates this reaction with composite coupled cluster protocols, including a novel spin-flip equation of motion coupled cluster scheme, combined with pragmatic two-dimensional master equation simulations of the resulting rate coefficients. These improved calculations predict the CH + N 2 rate coefficient between the two more recent previous theoretical results and reduce the uncertainties of the best theoretical models of this reaction to less than a factor of 1.3. Additionally, we provide a closer theoretical investigation of the simultaneous dependence of the CH + N 2 rate coefficient on pressure and temperature, and affirm that collisionally stabilized HNCN, another potential source of prompt NO, emerges as an appreciable product of this reaction under conditions relevant to automotive internal combustion engines and aircraft gas turbine engines.

Nguyen, Thanh Lam [Univ. of Florida, Gainesville, ↗

How Does Feedback Affect the Star Formation Histories of Galaxies?

Star formation in galaxies is regulated by the interplay of a range of processes that shape the multiphase gas in the interstellar and circumgalactic media. Using the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) suite of cosmological simulations, we study the effects of varying feedback and cosmology on the average star formation histories (SFHs) of galaxies at z ∼ 0 across the IllustrisTNG, SIMBA, and ASTRID galaxy formation models. We find that galaxy SFHs in all three models are sensitive to changes in stellar feedback, which affect the efficiency of baryon cycling and the rates at which central black holes grow, whereas the effects of varying active galactic nucleus (AGN) feedback depend on model-specific implementations of black hole seeding, accretion, and feedback. We also find strong interaction terms that couple stellar and AGN feedback, usually by regulating the amount of gas available for the central black hole to accrete. Using a double power law to describe the average SFHs, we derive a general set of equations relating the shape of the SFHs to physical quantities like baryon fraction and black hole mass across all three models. We find that a single set of equations (albeit with different coefficients) can describe the SFHs across all three CAMELS models, with cosmology dominating the SFH at early times, followed by halo accretion, and feedback and baryon cycling at late times. Galaxy SFHs provide a novel, complementary probe to constrain cosmology and feedback, and can connect the observational constraints from current and upcoming galaxy surveys with the physical mechanisms responsible for regulating galaxy growth and quenching.

Iyer, Kartheik G. [Columbia Univ., New York, NY (U↗