Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

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 271 records · Page 15

A Spatiotemporal-Aware Weighting Scheme for Improving Climate Model Ensemble Predictions

Multimodel ensembling has been widely used to improve climate model predictions, and the improvement strongly depends on the ensembling scheme. In this work, we propose a Bayesian neural network (BNN) ensembling method, which combines climate models within a Bayesian model averaging framework, to improve the predictive capability of model ensembles. Our proposed BNN approach calculates spatiotemporally varying model weights and biases by leveraging individual models' simulation skill, calibrates the ensemble prediction against observations by considering observation data uncertainty, and quantifies epistemic uncertainty when extrapolating to new conditions. More importantly, the BNN method provides interpretability about which climate model contributes more to the ensemble prediction at which locations and times. Thus, beyond its predictive capability, the method also brings insights and understanding of the models to guide further model and data development. In this study, we design experiments using an ensemble of CMIP6 climate model simulations to illustrate the BNN ensembling method's capability with respect to prediction accuracy, interpretability, and uncertainty quantification (UQ). We demonstrate that BNN can correctly assign larger weights to the regions and seasons where the individual model fits the observation better. Moreover, its offered interpretability is consistent with our understanding of localized climate model performance. Additionally, BNN shows an increasing uncertainty when the prediction is farther away from the period with constrained data, which appropriately reflects our trustworthiness of the models in the changing climate.

54 ENVIRONMENTAL SCIENCES↗

Characterizing and mitigating coherent errors in a trapped ion quantum processor using hidden inverses

Quantum computing testbeds exhibit high-fidelity quantum control over small collections of qubits, enabling performance of precise, repeatable operations followed by measurements. Currently, these noisy intermediate-scale devices can support a sufficient number of sequential operations prior to decoherence such that near term algorithms can be performed with proximate accuracy (like chemical accuracy for quantum chemistry problems). While the results of these algorithms are imperfect, these imperfections can help bootstrap quantum computer testbed development. Demonstrations of these algorithms over the past few years, coupled with the idea that imperfect algorithm performance can be caused by several dominant noise sources in the quantum processor, which can be measured and calibrated during algorithm execution or in post-processing, has led to the use of noise mitigation to improve typical computational results. Conversely, benchmark algorithms coupled with noise mitigation can help diagnose the nature of the noise, whether systematic or purely random. Here, we outline the use of coherent noise mitigation techniques as a characterization tool in trapped-ion testbeds. We perform model-fitting of the noisy data to determine the noise source based on realistic physics focused noise models and demonstrate that systematic noise amplification coupled with error mitigation schemes provides useful data for noise model deduction. Further, in order to connect lower level noise model details with application specific performance of near term algorithms, we experimentally construct the loss landscape of a variational algorithm under various injected noise sources coupled with error mitigation techniques. This type of connection enables application-aware hardware codesign, in which the most important noise sources in specific applications, like quantum chemistry, become foci of improvement in subsequent hardware generations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Conductivity and dissociation in liquid metallic hydrogen and implications for planetary interiors

Liquid metallic hydrogen (LMH) is the most abundant form of condensed matter in our solar planetary structure. The electronic and thermal transport properties of this metallic fluid are of fundamental interest to understanding hydrogen’s mechanism of conduction, atomic or pairing structure, as well as the key input for the magnetic dynamo action and thermal models of gas giants. Here, we report spectrally resolved measurements of the optical reflectance of LMH in the pressure region of 1.4–1.7 Mbar. We analyze the data, as well as previously reported measurements, using the free-electron model. Fitting the energy dependence of the reflectance data yields a dissociation fraction of 65 ± 15%, supporting theoretical models that LMH is an atomic metallic liquid. We determine the optical conductivity of LMH and find metallic hydrogen’s static electrical conductivity to be 11,000–15,000 S/cm, substantially higher than the only earlier reported experimental values. Here, the higher electrical conductivity implies that the Jovian and Saturnian dynamo are likely to operate out to shallower depths than previously assumed, while the inferred thermal conductivity should provide a crucial experimental constraint to heat transport models.

08 HYDROGEN↗

The effect of relative humidity and temperature on the response of stomatal conductance to vapor pressure deficit in tropical trees

Understanding how leaf gas exchange responds to changes in vapor pressure deficit (VPD) is key to predicting tropical forest resilience to climate change. Stomata regulate leaf water and CO2 diffusion, and respond to changes in temperature and relative humidity (RH), two drivers of VPD. At high temperatures, the cuticular pathway may also become significant and participate in the overall leaf conductance. Here, we measured gas exchange under light and dark conditions to investigate the stomatal and cuticular responses to temperature and RH on detached branches in five tropical tree species. Leaf conductance in the dark, when stomata are essentially closed, was not markedly impacted by temperature and RH, suggesting a minimal response of the cuticular pathway to these conditions. We compared six steady-state conductance models incorporating different effects of photosynthesis and evaporative demand on stomatal control. All models performed well (residual standard deviation, σ, < 0.025 mol m-2 s-1), but the best-fitting model (σ = 0.017 mol m-2 s-1) used a nonlinear relationship between photosynthesis and stomatal conductance and incorporated RH rather than vapor pressure difference. All models overestimated the steady-state conductance at high RH. Furthermore, leaf conductance immediately increased after a decrease in RH (wrong-way response), but not after an increase in leaf temperature. This suggests that the mechanisms underlying stomatal response to VPD need further investigation. The non-linear coupling between photosynthetic rate and stomatal conductance indicates a sharper physiological response than previously acknowledged.

Lamour, Julien↗

Estimating Regional Methane Emissions Through Atmospheric Measurements and Inverse Modeling

In this report we describe an enhanced methodology for performing stochastic Bayesian inversions of atmospheric trace gas inversions that allows the time variation of model parameters to be inferred. We use measurements of methane atmospheric mixing ratio made in Livermore, California along with atmospheric transport modeling and published prior estimates of emissions to estimate the regional emissions of methane and the temporal variations in inferred bias parameters. We compute Bayesian model evidence and continuous rank probability score to optimize the model with respect to temporal resolution. Using two different emissions inventories, we perform inversions for a series of models with increasing temporal resolution in the model bias representation. We show that temporal variation in the model bias can improve the model fit and can also increase the likelihood that the parameterization is appropriate, as measured by the Bayesian model evidence.

54 ENVIRONMENTAL SCIENCES↗

Acoustic Tomography of the Atmosphere: A Large-Eddy Simulation Sensitivity Study

Accurate measurement of atmospheric turbulent fluctuations is critical for understanding environmental dynamics and improving models in applications such as wind energy. Advanced remote sensing technologies are essential for capturing instantaneous velocity and temperature fluctuations. Acoustic tomography (AT) offers a promising approach that utilizes sound travel times between an array of transducers to reconstruct turbulence fields. This study presents a systematic evaluation of the time-dependent stochastic inversion (TDSI) algorithm for AT using synthetic travel-time measurements derived from large-eddy simulation (LES) fields under both neutral and convective atmospheric boundary-layer conditions. Unlike prior work that relied on field observations or idealized fields, the LES framework provides a ground-truth atmospheric state, enabling quantitative assessment of TDSI retrieval reliability, sensitivity to travel-time measurement noise, and dependence on covariance model parameters and temporal data integration. A detailed sensitivity analysis was conducted to determine the best-fit model parameters, identify the tolerance thresholds for parameter mismatch, and establish a maximum spatial resolution. The TDSI algorithm successfully reconstructed large-scale velocity and temperature fluctuations with root mean square errors ( RMSE s) below 0.35 m/s and 0.12 K, respectively. Spectral analysis established a maximum spatial resolution of approximately 1.4 m, and reconstructions remained robust for travel-time measurement uncertainties up to 0.002 s. These findings provide critical insights into the operational limits of TDSI and inform future applications of AT for atmospheric turbulence characterization and system design.

17 WIND ENERGY↗

Traffic safety analysis and model updating for freeways using Bayesian method

Freeway crash prediction models are the basic of traffic safety research, yet crash occurrence and the influencing factors change over time. In order to make sure the implemented safety models fit the current traffic environment, this study conducts a comparative analysis of 2017 and 2020 datasets collected from freeways in Suzhou, China. Herein, considering the spatial correlation among analysis units and the hierarchical data structure, a Bayesian conditional autoregressive negative binomial (CAR-NB) model and a Bayesian hierarchical CAR-NB (HCAR-NB) model were used to explore the safety influencing factors, and a traditional NB model was developed for further comparison. To update the HCAR-NB model from 2017 to 2020, Bayesian inference with informative priors was used to improve its goodness of fit and efficiency. Preliminary results showed that 1) the HCAR-NB model outperformed the NB model and CAR-NB model in prediction accuracy, and 2) the number of crashes was significantly correlated with average speed, speed variance, road segment length, number of lanes, and presence of ramps. The potential for safety improvement (PSI) method was applied to the modeling results to identify hotspots for the two years. The results confirmed that the hotspots spatiotemporally shifted among the freeways. The proposed crash prediction model and updating method are expected to assist implementation of informed countermeasures for freeway safety improvement.

97 MATHEMATICS AND COMPUTING↗

The active CGCG 077-102 NED02 galaxy within the Abell 2063 galaxy cluster

Context.Within the framework of investigating the link between the central super massive black holes in the cores of galaxies and the galaxies themselves, we detected a variable X-ray source in the center of CGCG 077-102 NED02, which is a member of the CGCG 077-102 galaxy pair within the Abell 2063 cluster of galaxies. Aims.Our goal is to combine X-ray and optical data to demonstrate that this object harbors an active super massive black hole in its core, and to relate this to the dynamical status of the galaxy pair within the Abell 2063 cluster. Methods.We usedChandraandXMM-Newtonarchival data to derive the X-ray spectral shape and variability. We also obtained optical spectroscopy to detect the expected emission lines that are typically found in active galactic nuclei. Finally, we used public ZTF imaging data to investigate the optical variability. Results.There is no evidence of multiple X-ray sources or extended components within CGCG 077-102 NED02. Single X-ray spectral models fit the source well. We detect significant, nonrandom inter-observation 0.5–10 keV X-ray flux variabilities, for observations separated by ∼4 days for short-term variations and by up to ∼700 days for long-term variations. Optical spectroscopy points toward a passive galaxy for CGCG 077-102 NED01 and a Seyfert for CGCG 077-102 NED02. The classification of CGCG 077-102 NED02 is also consistent with its X-ray luminosity of over 10 42 erg s −1 . We do not detect short-term variability in the optical ZTF light curves. However, we find a significant long-term stochastic variability in theg-band that can be well described by the damped random walk model with a best-fit characteristic damping timescale ofτ DRW = 30 −12 +28 days. Finally, the CGCG 077-102 galaxy pair is deeply embedded within the Abell 2063 potential, with a long enough history within this massive structure to have been affected by the influence of this cluster for a long time. Conclusions.Our observations point toward a moderately massive black hole in the center of CGCG 077-102 NED02 of ∼10 6 M ⊙ . As compared to another similar pair in the literature, CGCG 077-102 NED02 is not heavily obscured, perhaps because of the surrounding intracluster medium ram-pressure stripping.

Astronomy & Astrophysics↗

On the relaxation behaviors of slow and classical glitches: Observational biases and their opposite recovery trends

We study the pulsar timing properties and the data analysis methods of glitch recoveries. In some cases one first fits the times of arrival (TOAs) to obtain the 'time-averaged' frequency ν and its first derivative ν-dot , and then fits models to them. However, our simulations show that ν and ν-dot obtained in this way are systematically biased unless the time intervals between the nearby data points of TOAs are smaller than about 10{sup 4} s, which is much shorter than typical observation intervals. Alternatively, glitch parameters can be obtained by fitting the phases directly with relatively smaller biases, but the initial recovery time scale is usually chosen by eye, which may introduce a strong bias. We also construct a phenomenological model by assuming a pulsar spin-down law of ν-dot ν{sup −3}=−H{sub 0}G(t) with G(t) = 1 + κe {sup –t/τ} for a glitch recovery, where H {sub 0} is a constant and κ and τ are the glitch parameters to be found. This model can reproduce the observed data of slow glitches from B1822–09 and the giant classical glitch of B2334+61, with κ < 0 or κ > 0, respectively. We then use this model to simulate TOA data and test several fitting procedures for a glitch recovery. The best procedure is (1) to use a very high order polynomial (e.g., to 50th order) to precisely describe the phase, (2) obtain ν(t) and ν-dot (t) from the polynomial, then (3) obtain the glitch parameters from ν(t) or ν-dot (t). Finally, the uncertainty in the starting time t {sub 0} of a classical glitch causes uncertainties in some glitch parameters, but less so for a slow glitch, t {sub 0}, of which can be determined from data.

79 ASTRONOMY AND ASTROPHYSICS↗

NINE-YEAR WILKINSON MICROWAVE ANISOTROPY PROBE (WMAP) OBSERVATIONS: FINAL MAPS AND RESULTS

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail. We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground-reduced CMB maps are presented. We now implement an optimal C {sup –1} weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained ΛCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N {sub eff} = 3.84 ± 0.40). The model fit also implies that the age of the universe is t {sub 0} = 13.772 ± 0.059 Gyr, and the fit Hubble constant is H {sub 0} = 69.32 ± 0.80 km s{sup –1} Mpc{sup –1}. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n{sub s} = 0.9608 ± 0.0080); and the universe is close to flat/Euclidean (Ω{sub k} = -0.0027{sup +0.0039}{sub -0.0038}). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter ΛCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness of fit, the inflationary flat ΛCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement uncertainty of rossi-alpha neutron experiments

Rossi-alpha neutron experiments are used to estimate the prompt neutron decay constant of a fissile assembly, a quantity of widespread interest in applications including in nuclear nonproliferation and criticality safety. This work develops a mathematical model to efficiently estimate measurement uncertainty of Rossi-alpha neutron experiments inferred from a two-exponential fit model with histogram binning. The derived uncertainty estimates were validated using repeated Rossi-alpha measurements of a subcritical, 4.5-kg sphere of weapons-grade, alpha-phase plutonium with nickel, copper, tungsten, and polyethylene reflectors. The estimates of uncertainty for the histogram data produced by the model were conservative and agree with the reference uncertainties within noise. The estimates of the prompt neutron decay constant uncertainty agreed with the reference uncertainties within one standard deviation. In conclusion, the proposed model will reduce total measurement times, ultimately reducing operational and procedural costs in application.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Kinetic analysis of the interactions between calcium ferrite and coal char for chemical looping gasification applications: Identifying reduction routes and modes of oxygen transfer

Chemical Looping Gasification (CLG) is an emerging technology that shows promise for efficient coal gasification by eliminating the need for energy intensive gas separations to achieve a non-nitrogen diluted syngas stream. Oxygen from oxygen carriers, such as CaFe 2 O 4 , are used for coal gasification in place of conventionally produced gaseous oxygen from cryogenic separation of air. These oxygen carriers are unique for their ability to selectively oxidize coal to form syngas and show limited reactivity with syngas components (H 2 , CO). To gain a deeper understanding of how these unique oxygen carriers perform and to offer a first attempt at the reaction modeling of solid mediated interactions of this nature, this study was carried out to determine the kinetic parameters associated with the selective oxidation of coal derived char (Wyodak and Illinois #6) with a metal ferrite, CaFe 2 O 4 . Using thermogravimetric analysis (TGA) coupled with mass spectrometry, the selective oxygen release of metal ferrite in the presence of char by proximal contact was examined. The application of combinatory model fitting approaches was used to describe controlling resistances during oxygen release. A combination of the modified shrinking core model (SCM) with planar oxygen ion diffusion control and reaction order based models was used for kinetic parameter determination. CaFe 2 O 4 particle size plays a major role in the prevailing mode of oxygen release. Particle sizes on the order of 40–50 μm tend to favor first order kinetically controlled regimes independent of geometric and diffusion controls. The probability for oxygen ion diffusion controlling regimes increased when the particle size range of the oxygen carrier was increased up to 350 μm. Char type also impacted the prevalence of the controlling regime. Higher ranked chars react in a slower manner, limiting the gradient for oxygen ion release from the oxygen carrier. Activation energies determined for this process range from 120–200kJ/mol and oxygen ion diffusion coefficients are on the order of 10-8 cm 2 /s. It is suggested that oxygen ion movement is regulated by lattice diffusion out of partially reduced phases (Ca 2 Fe 2 O 5 ) and through reduced outer layers composed of CaO and Fe. The controlled movement of oxygen ions influences the rate of carbon oxidation in the char and therefore the selectivity towards partial oxidation products, which are desirable in CLG applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A fusion relevant data-driven engineering void swelling model for 9Cr tempered martensitic steels

The UCSB database on cavity evolution in 9-12Cr tempered martensitic steels (TMS), includes the results for both dual heavy and helium ion (DII), and High Flux Isotope Reactor (HFIR) in situ helium injection (ISHI) neutron irradiations at 500°C. These results were combined with literature single ion and fission neutron irradiation data to derive a model for the void volume fraction, f v , as a function of displacements per atom (dpa) and transmutant helium concentrations in atomic parts per million (appm). The scientific foundation for the paper is described in a companion paper entitled “Cavity Evolution and Void Swelling in Dual Ion Irradiated Tempered Martensitic Steels”. Here, in this study, we show that f v (dpa, He/dpa) is described by the incubation dose, dpa i , for the onset of void growth, and the post-incubation growth rate, f v ’(%/dpa). Both dpa i and f v ’ decrease with increasing He/dpa at > ~ 5. The dpa i is also lower for the ISHI neutron irradiations at the same He/dpa. Single heavy ion and fission reactor neutron irradiations, with low He/dpa ratios, have a much larger dpa i . Based on a combined analysis of DII, single ion, ISHI and fission neutron data, we further show that the post-incubation f v data analyzed here have a common empirical curve shape, with f v ’ reaching up to ~ 0.2%/dpa at very high dpa. We also show that f v ’ can be predicted based on a physical model of defect partitioning between evolving sinks. At 500°C and fusion relevant He/dpa ≈ 10, the best-fit model predicts nominal swelling, S = f v /(1-f v ), of ~ 1.1, 4.9 and 16% at 50, 100 and 200 dpa, respectively. The physically motivated, data-driven model includes estimated uncertainties for both dpa i and f v ’.

36 MATERIALS SCIENCE↗

A novel invariant mass method to isolate resonance backgrounds from the chiral magnetic effect

The Chiral Magnetic Effect (CME) refers to charge separation along a strong magnetic field, due to topological charge fluctuations in QCD. Charge correlation (Δγ) signals consistent with CME have been first observed almost a decade ago. It has also been known since then that the Δγ is contaminated by a major background from resonance decays coupled with elliptic flow. In this contribution, we propose differential Δγ measurements as function of the pair invariant mass (m inv ). The Δγ in the high minv region is essentially free of resonance backgrounds. In the low m inv region, the Δγ backgrounds show resonance peaks. The CME signal, presumably smooth in m inv , may thus be extracted from a two-component model fit. We demonstrate the feasibility and effectiveness of this novel method by using the AMPT and toy-model Monte-Carlo simulations. We also discuss an application of the method in data analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Capturing the 3D motion of an infalling galaxy via fluid dynamics

The Fornax Cluster is the nearest (⩽20 Mpc) galaxy cluster in the southern sky. NGC 1404 is a bright elliptical galaxy falling through the intracluster medium (ICM) of the Fornax Cluster. The sharp leading edge of NGC 1404 forms a classical “cold front” that separates 0.6 keV dense interstellar medium and 1.5 keV diffuse ICM. We measure the angular pressure variation along the cold front using a very deep (670 ks) Chandra X-ray observation. We are taking the classical approach—using stagnation pressure to determine a substructure’s speed—to the next level by not only deriving a general speed but also directionality, which yields the complete velocity field as well as the distance of the substructure directly from the pressure distribution. We find a hydrodynamic model consistent with the pressure jump along NGC 1404's atmosphere measured in multiple directions. The best-fit model gives an inclination of 33° and a Mach number of 1.3 for the infall of NGC 1404, in agreement with complementary measurements of the motion of NGC 1404. Our study demonstrates the successful treatment of a highly ionized ICM as ideal fluid flow, in support of the hypothesis that magnetic pressure is not dynamically important over most of the virial region of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗