Gas Density Profiles Research Articles

Can the splashback radius be an observable boundary of galaxy clusters?

The splashback radius was proposed as a physically motivated boundary of clusters as it sets the limit between the infalling and the orbitally dominated regions. However, galaxy clusters are complex objects connected to filaments of the cosmic web from which they accrete matter that disturbs them and modifies their morphology. In this context, estimating the splashback radius and the cluster boundary becomes challenging. In this work, we use a constrained hydrodynamical simulation replicating the Virgo cluster embedded in its large-scale structure to investigate the impact of its local environment on the splashback radius estimate. We identify the splashback radius from 3D radial profiles of dark matter density, gas density, and pressure in three regions representative of different dynamical states: accretion from spherical collapse, filaments, and matter outflow. We also identify the splashback radius from 2D-projected radial profiles of observation-like quantities: mass surface density, emission measure, and Compton-y. We show that the splashback radius mainly depends on the dynamics in each region and the physical processes traced by the different probes. We find multiple values for the splashback radius ranging from 3.3 ± 0.2 to 5.5 ± 0.3 Mpc. In particular, in the regions of collapsing and outflowing materials, the splashback radii estimated from gas density and pressure radial profiles overestimate that of the dark matter density profiles, which is considered the reference value given that the splashback radius was originally defined from dark matter simulations in pioneering works. Consequently, caution is required when using the splashback radius as a boundary of clusters, particularly in the case of highly disturbed clusters like Virgo. We conclude with a discussion of the detection of the splashback radius from pressure radial profiles, which could be more related to an accretion shock, and its detection from stacked radial profiles.

Simulations of nozzle gas flow and gas-puff Z-pinch implosions on the Weizmann Z-pinch

We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.

Correction to: X-ray surface brightness and gas density profiles of galaxy clusters up to 3 × R500c with SRG/eROSITA

Correction to: X-ray surface brightness and gas density profiles of galaxy clusters up to 3 × <i>R</i>500c with SRG/eROSITA

Forecast constraints on the baryonic feedback effect from the future kinetic Sunyaev-Zel'dovich effect detection

The baryonic feedback effect is an important systematic error in the weak lensing (WL) analysis. It contributes partly to the S 8 tension in the literature. With the next generations of large scale structure (LSS) and CMB experiments, the high signal-to-noise kinetic Sunyaev-Zel'dovich (kSZ) effect detection can tightly constrain the baryon distribution in and around dark matter halos, and quantify the baryonic effect in the weak lensing statistics. In this work, we apply the Fisher matrix technique to predict the future kSZ constraints on 3 kSZ-sensitive Baryon Correction Model (BCM) parameters. Our calculations show that, in combination with next generation LSS surveys, the 3rd generation CMB experiments such as AdvACT and Simon Observatory can constrain the matter power spectrum damping S(k) to the precision of σ S (k) < 0.8%√(37.8 Gpc3h-3-3/V) at k ≲ 10h/Mpc, where V is the overlapped survey volume between the future LSS and CMB surveys. For the 4th generation CMB surveys such as CMB-S4 and CMB-HD, the constraint will be enhanced to σ S (k) < 0.4%√(37.8 Gpc3h-3/V). If extra-observations, e.g. X-ray detection and thermal SZ observation, can effectively fix the gas density profile slope parameter δ, the constraint on S(k) will be further boosted to σ S (k) < 0.3%√(37.8 Gpc3 h -3/V) and σ S (k) < 0.1%√(37.8 Gpc3 h -3/V) for the 3rd and 4th generation CMB surveys.

Predicting the Scaling Relations between the Dark Matter Halo Mass and Observables from Generalised Profiles II: Intracluster Gas Emission

Abstract We investigate the connection between a cluster’s structural configuration and observable measures of its gas emission that can be obtained in X-ray and Sunyaev-Zeldovich (SZ) surveys. We present an analytic model for the intracluster gas density profile: parameterised by the dark matter halo’s inner logarithmic density slope, α, the concentration, c, the gas profile’s inner logarithmic density slope, ε, the dilution, d, and the gas fraction, η, normalised to cosmological content. We predict four probes of the gas emission: the emission-weighted, TX, and mean gas mass-weighted, Tmg, temperatures, and the spherically, Ysph, and cylindrically, Ycyl, integrated Compton parameters. Over a parameter space of clusters, we constrain the X-ray temperature scaling relations, M200 – TX and M500 – TX, within 57.3% and 41.6%, and M200 – Tmg and M500 – Tmg, within 25.7% and 7.0%, all respectively. When excising the cluster’s core, the M200 – TX and M500 – TX relations are further constrained, to within 31.3% and 17.1%, respectively. Similarly, we constrain the SZ scaling relations, M200 – Ysph and M500 – Ysph, within 31.1% and 17.7%, and M200 – Ycyl and M500 – Ycyl, within 25.2% and 22.0%, all respectively. The temperature observable Tmg places the strongest constraint on the halo mass, whilst TX is more sensitive to the parameter space. The SZ constraints are sensitive to the gas fraction, whilst insensitive to the form of the gas profile itself. In all cases, the halo mass is recovered with an uncertainty that suggests the cluster’s structural profiles only contribute a minor uncertainty in its scaling relations.

Inferring the dark matter splashback radius from cluster gas and observable profiles in the FLAMINGO simulations

ABSTRACT The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev–Zel’dovich (SZ) observations of the hot intracluster gas. Here, we use the FLAMINGO simulations to investigate whether an extremum gradient in a similar position to the splashback radius is predicted to occur in the cluster gas profiles. We find that the minimum in the gradient of the stacked 3D gas density and pressure profiles, and the maximum in the gradient of the entropy profile, broadly align with the splashback feature though there are significant differences. While the dark matter splashback radius varies with specific mass accretion rate, in agreement with previous work, the radial position of the deepest minimum in the log-slope of the gas density is more sensitive to halo mass. In addition, we show that a similar minimum is also present in projected 2D pseudo-observable profiles: emission measure (X-ray), Compton-y (SZ), and surface mass density (weak lensing). We find that the latter traces the dark matter results reasonably well albeit the minimum occurs at a slightly smaller radius. While results for the gas profiles are largely insensitive to accretion rate and various observable proxies for dynamical state, they do depend on the strength of the feedback processes.

The cosmic baryon partition between the IGM and CGM in the SIMBA simulations

ABSTRACT We use the simba suite of cosmological hydrodynamical simulations to investigate the importance of various stellar and active galactic nuclei (AGN) feedback mechanisms in partitioning the cosmic baryons between the intergalactic (IGM) and circumgalactic (CGM) media in the z ≤ 1 Universe. We identify the AGN jets as the most prominent mechanism for the redistribution of baryons between the IGM and CGM. In contrast to the full feedback models, deactivating AGN jets results in ≈20 per cent drop in fraction of baryons residing in the IGM and a consequent increase of CGM baryon fraction by ≈50 per cent. We find that stellar feedback modifies the partition of baryons on a 10 per cent level. We further examine the physical properties of simulated haloes in different mass bins, and their response to various feedback models. On average, a sixfold decrease in the CGM mass fraction due to the inclusion of feedback from AGN jets is detected in $10^{12}\, {\rm M}_{\odot } \le M_{\rm 200} \le 10^{14}\, {\rm M}_{\odot }$ haloes. Examination of the average radial gas density profiles of $M_{200} \gt 10^{12}\, {\rm M}_{\odot }$ haloes reveals up to an order of magnitude decrease in gas densities due to the AGN jet feedback. We compare gas density profiles from simba simulations to the predictions of the modified Navarro–Frenk–White model, and show that the latter provides a reasonable approximation within the virial radii of the full range of halo masses, but only when rescaled by the appropriate mass-dependent CGM fraction of the halo. The relative partitioning of cosmic baryons and, subsequently, the feedback models can be constrained observationally with fast radio bursts in upcoming surveys.

Mass bias in clusters of galaxies: Projection effects on the case study of Virgo replica

When measuring the observed pressure, density, or temperature profiles of the intracluster gas, and hence the mass of clusters of galaxies, projection effects or departures from the spherical symmetry hypothesis may induce biases. To estimate how strongly the cluster’s observed properties depend on the direction of observation, we use a constrained hydrodynamical simulation of the Virgo cluster that replicates the actual cluster of galaxies. In this case study, we analysed Virgo properties when projected in different directions, including along the Milky Way-Virgo axis, which mimics our observation direction. We compared the hydrostatic mass and the hydrostatic mass bias from the projection along the different observation directions to that derived from the 3D simulation. We show that projection effects impact the determination of Virgo mass. We particularly demonstrate that the mass and pressure along the line of sight correlate with the 2D- and 3D-deprojected electron density and pressure profiles intensity and thus impact the derived hydrostatic mass. We also show that the deviations to the hydrostatic equilibrium induced by pressure discontinuities within the cluster are emphasised by the deprojection process and thus make the hydrostatic mass estimation invalid at these radii.

The formation of compact massive relic galaxies in MOND

ABSTRACT Compact massive relic galaxies are a class of galaxies that exhibit characteristics suggesting they have remained largely unchanged since their initial formation, making them ‘relics’ of the early Universe. These galaxies represent a distinct class characterized by strongly peaked high-velocity dispersion profiles with large rotational velocities. This study investigates the formation of such galaxies within the framework of Milgromian dynamics (MOND), offering a unique perspective on their origin without invoking the presence of cold or warm dark matter. Our analysis focuses on the collapse dynamics of isolated non-rotating post-big bang gas clouds, revealing kinematic and density profiles comparable to observed compact massive relic galaxies like NGC 1277, Mrk 1216, and PGC 032873. The findings underscore the natural emergence of compact massive relic galaxies within a MOND-based Universe, providing valuable insights into the interplay between gravitational dynamics and galaxy formation.

Determining the baryon impact on the matter power spectrum with galaxy clusters

ABSTRACT The redistribution of baryonic matter in massive haloes through processes like active galactic nuclei feedback and star formation leads to a suppression of the matter power spectrum on small scales. This redistribution can be measured empirically via the gas and stellar mass fractions in galaxy clusters, and leaves imprints on their electron density profiles. We constrain two semi-analytical baryon correction models with a compilation of recent Bayesian population studies of galaxy groups and clusters sampling a mass range above ∼3 × 1013 M⊙, and with cluster gas density profiles derived from deep, high-resolution X-ray observations. We are able to fit all the considered observational data, but highlight some anomalies in the observations. The constraints allow us to place precise, physically informed priors on the matter power spectrum suppression. At a scale of k = 1 h Mpc−1 we find a suppression of $0.042^{+0.012}_{-0.014}$ ($0.049^{+0.016}_{-0.012}$), while at k = 3 h Mpc−1 we find $0.184^{+0.026}_{-0.031}$ ($0.179^{+0.018}_{-0.020}$), depending on the model used. In our fiducial setting, we also predict at 97.5 per cent credibility, that at scales k &amp;lt; 0.37 h Mpc−1 baryon feedback impacts the matter power less than 1 per cent. This puts into question if baryon feedback is the driving factor for the discrepancy between cosmic shear and primary CMB results. We independently confirm results on this suppression from small-scale cosmic shear studies, while we exclude some hydro-dynamical simulations with too strong and too weak baryonic feedback. Our empirical prediction of the power spectrum suppression shows that studies of galaxy groups and clusters will be instrumental in unlocking the cosmological constraining power of future cosmic shear experiments like Euclid and Rubin-LSST, and invites further investigation of the baryon correction models.

Ion acceleration from the interaction of ultrahigh-intensity laser pulses with near-critical density, nonuniform gas targets

Using one-dimensional, long-timescale particle-in-cell simulations, we study the processes of ion acceleration from the interaction of ultraintense (1020 W cm−2), ultrashort (30 fs) laser pulses with near-critical, nonuniform gas targets. The considered initially neutral, nitrogen gas density profiles mimic those delivered by an already developed noncommercial supersonic gas shock nozzle: they have the generic shape of a narrow (20 μm wide) peak superimposed on broad (∼1 mm, ∼180 μm scale length), exponentially decreasing ramps. While keeping its shape constant, we vary its absolute density values to identify the interaction conditions leading to collisionless shock-induced ion acceleration in the gas density ramps. We find that collisionless electrostatic shocks (CES) form when the laser pulse is able to shine through the central density peak and deposit a few 10% of its energy into it. Under our conditions, this occurs for a peak electron density between 0.35 nc and 0.7 nc. Moreover, we show that the ability of the CES to reflect the upstream ions is highly sensitive to their charge state and that the laser-induced electron pressure gradients mainly account for shock generation, thus highlighting the benefit of using sharp gas profiles, such as those produced by shock nozzles.

Impact of the neutral molecule trapping on beam lifetime and beam profile

Recent investigations have shown that neutral molecules may be trapped by the electromagnetic field of a charged particle beam. This effect changes the residual gas distribution according to a characteristic temperature of the system consisting of beam and residual gas. Here, we investigate the effect of the modified gas density profile on the beam lifetime.

Proton and helium ions acceleration in near-critical density gas targets by short-pulse Ti:Sa PW-class laser

The ability to quickly refresh gas-jet targets without cycling the vacuum chamber makes them a promising candidate for laser-accelerated ion experiments at high repetition rate. Here we present results from the first high repetition rate ion acceleration experiment on the VEGA-3 PW-class laser at CLPU. A near-critical density gas-jet target was produced by forcing a 1000 bar H $_2$ and He gas mix through bespoke supersonic shock nozzles. Proton energies up to 2 MeV were measured in the laser forward direction and 2.2 MeV transversally. He $^{2+}$ ions up to 5.8 MeV were also measured in the transverse direction. To help maintain a consistent gas density profile over many shots, nozzles were designed to produce a high-density shock at distances larger than 1 mm from the nozzle exit. We outline a procedure for optimizing the laser–gas interaction by translating the nozzle along the laser axis and using different nozzle materials. Several tens of laser interactions were performed with the same nozzle which demonstrates the potential usefulness of gas-jet targets as high repetition rate particle source.

Impact of Galaxy Clusters on the Propagation of Ultrahigh-energy Cosmic Rays

Galaxy clusters are the largest objects in the Universe kept together by gravity. Most of their baryonic content is made of a magnetized diffuse plasma. We investigate the impact of such a magnetized environment on the propagation of ultrahigh-energy cosmic rays (UHECRs). The intracluster medium (ICM) is described according to the self-similar assumption, in which gas density and pressure profiles are fully determined by the cluster mass and redshift. The magnetic field is scaled to the thermal components of the ICM under different assumptions. We model the propagation of UHECRs in the ICM using a modified version of the Monte Carlo code SimProp, where hadronic processes and diffusion in the turbulent magnetic field are implemented. We provide a universal parameterization that approximates the UHECR fluxes escaping from the environment as a function of the most relevant quantities, such as the mass of the cluster, the position of the source with respect to the center of the cluster, and the nature of the accelerated particles. We show that galaxy clusters are an opaque environment, especially for UHECR nuclei. The role of the most massive nearby clusters in the context of the emerging UHECR astronomy is finally discussed.

Probing the roles of orientation and multiscale gas distributions in shaping the obscuration of active galactic nuclei through cosmic time

ABSTRACT The origin of obscuration in active galactic nuclei (AGNs) is still an open debate. In particular, it is unclear what drives the relative contributions to the line-of-sight column densities from galaxy-scale and torus-linked obscuration. The latter source is expected to play a significant role in Unification Models, while the former is thought to be relevant in both Unification and Evolutionary models. In this work, we make use of a combination of cosmological semi-analytic models and semi-empirical prescriptions for the properties of galaxies and AGN, to study AGN obscuration. We consider a detailed object-by-object modelling of AGN evolution, including different AGN light curves (LCs), gas density profiles, and also AGN feedback-induced gas cavities. Irrespective of our assumptions on specific AGN LC or galaxy gas fractions, we find that, on the strict assumption of an exponential profile for the gas component, galaxy-scale obscuration alone can hardly reproduce the fraction of log (NH/cm−2) ≥ 24 sources at least at z ≲ 3. This requires an additional torus component with a thickness that decreases with luminosity to match the data. The torus should be present in all evolutionary stages of a visible AGN to be effective, although galaxy-scale gas obscuration may be sufficient to reproduce the obscured fraction with 22 &amp;lt; log (NH/cm−2) &amp;lt; 24 (Compton-thin, CTN) if we assume extremely compact gas disc components. The claimed drop of CTN fractions with increasing luminosity does not appear to be a consequence of AGN feedback, but rather of gas reservoirs becoming more compact with decreasing stellar mass.

Design, manufacturing, evaluation, and performance of a 3D-printed, custom-made nozzle for laser wakefield acceleration experiments.

Laser WakeField Acceleration (LWFA) is extensively used as a high-energy electron source, with electrons achieving energies up to the GeV level. The produced electron beam characteristics depend strongly on the gas density profile. When the gaseous target is a gas jet, the gas density profile is affected by parameters, such as the nozzle geometry, the gas used, and the backing pressure applied to the gas valve. An electron source based on the LWFA mechanism has recently been developed at the Institute of Plasma Physics and Lasers. To improve controllability over the electron source, we developed a set of 3D-printed nozzles suitable for creating different gas density profiles according to the experimental necessities. Here, we present a study of the design, manufacturing, evaluation, and performance of a 3D-printed nozzle intended for LWFA experiments.

Xe Recovery from Nuclear Power Plants Off-Gas Streams: Molecular Simulations of Gas Permeation through DD3R Zeolite Membrane.

Recent experimental work has shown zeolite membrane-based separation as a promising potential technology for Kr/Xe gas mixtures due to its much lower energy requirements in comparison to cryogenic distillation, the conventional separation method for such mixtures. Such a separation is also economically rewarding because Xe is in high demand, as a valuable product for many applications/processes. In this work, we have used Molecular Dynamics (MD) simulations to study the effects of different conditions, i.e., temperature, pressure, and gas feed composition, on Kr/Xe separation performance via DD3R zeolite membranes. We provide a comprehensive study of the permeation of the different gas species, density profiles, and diffusion coefficients. Molecular simulations show that if the feed is changed from pure Kr/Xe to an equimolar mixture, the Kr/Xe separation factor increases, which agrees with experiments. In addition, when Ar is introduced as a sweep gas, the adsorption of both Kr and Xe increases, while the permeation of pure Kr increases. A similar behavior is observed with equimolar mixtures of Kr/Xe with Ar as the sweep gas. High-separation Kr/Xe selectivity is observed at 50 atm and 425 K but with low total permeation rates. Changing pressure and temperature are found to have profound effects on optimizing the separation selectivity and the permeation throughput.

X-ray surface brightness and gas density profiles of galaxy clusters up to 3 × R500c with SRG/eROSITA

ABSTRACT Using the data of the SRG/eROSITA all-sky survey, we stacked a sample of ∼40 galaxy cluster images in the 0.3–2.3 keV band, covering the radial range up to 10 × R500c. The excess emission on top of the Galactic and extragalactic X-ray backgrounds and foregrounds is detected up to ∼3 × R500c. At these distances, the surface brightness of the stacked image drops below ∼1 per cent of the background. The density profile reconstructed from the X-ray surface brightness profile agrees well (within ∼30 per cent) with the mean gas profile found in numerical simulations, which predict the local gas overdensity of ∼ 20–30 at 3 × R500c and the gas fraction close to the universal value of $\frac{\Omega _b}{\Omega _m}\approx 0.15$ in the standard Λ cold dark matter model. Taking at face value, this agreement suggests that up to ∼3 × R500c the X-ray signal is not strongly boosted by the gas clumpiness, although a scenario with a moderately inhomogeneous gas cannot be excluded. A comparison of the derived gas density profile with the electron pressure profile based on the Sunyaev-Zeldovich (SZ) effect measurements suggests that by r ∼ 3 × R500c the gas temperature drops by a factor of ∼ 4–5 below the characteristic temperature of a typical cluster in the sample within R500c, while the entropy keeps growing up to this distance. Better constraints on the gas properties just beyond 3 × R500c should be possible with a sample larger than used for this pilot study.

The Three Hundred Project: the evolution of physical baryon profiles

ABSTRACT The distribution of baryons provides a significant way to understand the formation of galaxy clusters by revealing the details of its internal structure and changes over time. In this paper, we present theoretical studies on the scaled profiles of physical properties associated with the baryonic components, including gas density, temperature, metallicity, pressure and entropy as well as stellar mass, metallicity and satellite galaxy number density in galaxy clusters from z = 4 to z = 0 by tracking their progenitors. These mass-complete simulated galaxy clusters are coming from The Three Hundred with two runs: Gizmo-SIMBA and Gadget-X. Through comparisons between the two simulations, and with observed profiles that are generally available at low redshift, we find that (1) the agreements between the two runs and observations are mostly at outer radii r ≳ 0.3r500, in line with the self-similarity assumption. While Gadget-X shows better agreements with the observed gas profiles in the central regions compared to Gizmo-SIMBA; (2) the evolution trends are generally consistent between the two simulations with slightly better consistency at outer radii. In detail, the gas density profile shows less discrepancy than the temperature and entropy profiles at high redshift. The differences in the cluster centre and gas properties imply different behaviours of the AGN models between Gadget-X and Gizmo-SIMBA, with the latter, maybe too strong for this cluster simulation. The high-redshift difference may be caused by the star formation and feedback models or hydrodynamics treatment, which requires observation constraints and understanding.

A Bayesian chemical evolution model of the DustPedia galaxy M74

ABSTRACT We introduce a new, multizone chemical evolution model of the DustPedia galaxy M74, calibrated by means of Markov Chain Monte Carlo methods. We take into account the observed stellar and gas density profiles and use Bayesian analysis to constrain two fundamental parameters characterizing the gas accretion and star formation time-scale, i.e. the infall time-scale τ and the SF efficiency ν, respectively, as a function of galactocentric radius R. Our analysis supports an infall time-scale increasing with R and a star formation efficiency decreasing with R, thus supporting an ‘Inside-Out’ formation for M74. For both τ and ν, we find a weaker radial dependence than in the Milky Way. We also investigate the dust content of M74, comparing the observed dust density profile with the results of our chemical evolution models. Various prescriptions have been considered for two key parameters, i.e. the typical dust accretion time-scale τ0 and the mass of gas cleared out of dust by a supernova remnant, Mclear, regulating the dust growth and destruction rate, respectively. Two models with a different current balance between destruction and accretion, i.e. with an equilibrium and a dominion of accretion over destruction, can equally reproduce the observed dust profile of M74. This outlines the degeneracy between these parameters in shaping the interstellar dust content in galaxies. Our methods will be extended to more DustPedia galaxies to shed more light on the relative roles of dust production and destruction.