Distribution Of Baryons Research Articles

Cosmological simulations of the same spiral galaxy: connecting the dark matter distribution of the host halo with the subgrid baryonic physics

The role of baryonic physics, star formation and stellar feedback, in shaping the galaxies and their host halos is an evolving topic. The dark matter aspects are illustrated in this work by showing distribution features in a Milky Way sized halo.We focus on the halo morphology, geometry, and profile as well as the phase space distribution using one dark matter only and five hydrodynamical cosmological high-resolution simulations of the same halo with different subgrid prescriptions for the baryonic physics (Kennicut versus multi-freefall star formation and delayed cooling versus mechanical supernovae feedback). If some general properties like the relative halo-galaxy orientation are similar, the modifications of the gravitational potential due to the presence of baryons are found to induce different dark matter distributions (rounder and more concentrated halo). The mass density profile as well as the velocity distribution are modified distinctively according to the specific resulting baryonic distribution highlighting the variability of those properties (e.g inner power index from 1.3 to 1.8, broader speed distribution). The uncertainties on those features are of paramount importance for dark matter phenomenology, particularly when dealing with dark matter dynamics or direct and indirect detection searches. As a consequence, dark matter properties and prospects using cosmological simulations require improvement on baryonic physics description. Modeling such processes is a key issue not only for galaxy formation but also for dark matter investigations.

Reconstructing the extended structure of multiple sources strongly lensed by the ultra-massive elliptical galaxy SDSS J0100+1818

We study the total and baryonic mass distributions of the deflector SDSS J0100+1818 through a full strong lensing analysis. The system is composed of an ultra-massive early-type galaxy at z = 0.581, with a total stellar mass of (1.5 ± 0.3)×1012 M⊙ and a stellar velocity dispersion of (450 ± 40) km s−1, surrounded by ten multiple images of three background sources, two of which are spectroscopically confirmed at z = 1.880. We took advantage of high-resolution HST photometry and VLT/X-shooter spectroscopy to measure the positions of the multiple images and performed a strong lensing study with the software GLEE. We tested different total mass profiles for the lens and modeled the background sources first as point-like and then as extended objects. We successfully predict the positions of the observed multiple images and reconstruct over approximately 7200 HST pixels the complex surface brightness distributions of the sources. We measured the cumulative total mass profile of the lens and find a total mass value of (9.1 ± 0.1)×1012 M⊙, within the Einstein radius of approximately 42 kpc, and stellar-over-total mass fractions ranging from (49 ± 12)%, at the half-light radius (Re = 9.3 kpc) of the lens galaxy, to (10 ± 2)%, in the outer regions (R = 70 kpc). These results suggest that the baryonic mass component of SDSS J0100+1818 is very concentrated in its core and that the lens early-type galaxy (or group) is immersed in a massive dark matter halo, which allows it to act as a powerful gravitational lens, creating multiple images with exceptional angular separations. This is consistent with what has been found in other ultra-high-mass candidates at intermediate redshift. We also measured the physical sizes of the distant sources, resolving them down to a few hundred parsecs. Finally, we quantify and discuss a relevant source of systematic uncertainties on the reconstructed sizes of background galaxies, associated with the adopted lens total mass model.

Phase transition amplification of proton number fluctuations in nuclear collisions from a transport model approach

The time evolution of particle number fluctuations in nuclear collisions at intermediate energies ($E_{\rm lab} = 1.23-10A$ GeV) is studied by means of the UrQMD-3.5 transport model. The transport description incorporates baryonic interactions through a density-dependent potential. This allows for an implementation of a first order phase transition including a mechanically unstable region at large baryon density. The scaled variance of the baryon and proton number distributions is calculated in the central cubic spatial volume of the collisions at different times. A significant enhancement of fluctuations associated with the unstable region is observed. This enhancement persists to late times reflecting a memory effect for the fluctuations. The presence of the phase transition has a much smaller influence on the observable event-by-event fluctuations of protons in momentum space.

Baryonic post-processing of N-body simulations, with application to fast radio bursts

ABSTRACT Where the cosmic baryons lie in and around galactic dark matter haloes is only weakly constrained. We develop a method to quickly paint on models for their distribution. Our approach uses the statistical advantages of N-body simulations, while painting on the profile of gas around individual haloes in ways that can be motivated by semi-analytic models or zoom-in hydrodynamic simulations of galaxies. Possible applications of the algorithm include extragalactic dispersion measures to fast radio bursts (FRBs), the Sunyaev–Zeldovich effect, baryonic effects on weak lensing, and cosmic metal enrichment. As an initial application, we use this tool to investigate how the baryonic profile of foreground galactic-mass haloes affects the statistics of the dispersion measure (DM) towards cosmological FRBs. We show that the distribution of DM is sensitive to the distribution of baryons in galactic haloes, with viable gas profile models having significantly different probability distributions for DM to a given redshift. We also investigate the requirements to statistically measure the circumgalactic electron profile for FRB analyses that stack DM with impact parameter to foreground galaxies, quantifying the size of the contaminating ‘two-halo’ term from correlated systems and the number of FRBs for a high significance detection. Publicly available python modules implement our CGMBrush algorithm.

Splitting of proton-antiproton directed flow in relativistic heavy-ion collisions

The rapidity dependent directed flow of particles produced in a relativistic heavy ion collision can be generated in the hydrodynamic expansion of a tilted source. The asymmetry of the pressure leads to a build up of a directed flow of matter with respect to the collision axis. The experimentally observed ordering of the directed flow of baryons, pions and antibaryons can be described as resulting from the expansion of a baryon inhomogeneous fireball. An uneven distribution of baryons in the transverse plane leads to an asymmetry in the collective push for protons and antiprotons. Precise measurements of the collective flow as a function of rapidity could serve as a strong constraint on mechanism of baryon stopping in the early phase of the collision.

Understanding galaxy rotation curves with Verlinde’s emergent gravity

We present the results from the analysis of galaxy rotation curves with Verlinde’s emergent gravity. We use the data in the Spitzer Photometry and Accurate Rotation Curves database, which contains a sample of 175 nearby disk galaxies with 3.6 µm surface photometry and rotation curves. We compute the gravitational acceleration at different galactocentric radii expected from the baryon distribution of the galaxies with the emergent gravity, and compare it with the observed gravitational acceleration derived from galactic rotation curves. The predicted and observed accelerations agree well with a mean offset and a scatter by assuming a de Sitter Universe. These offset and scatter become smaller when we assume a more realistic Universe, quasi de Sitter Universe, as and . Our results suggest that Verlinde’s emergent gravity could be a good solution to the missing mass problem without introducing dark matter.

Extended surface brightness modeling of three sources strongly lensed by an ultra-massive elliptical galaxy

AbstractRare systems with multiple sources strongly lensed by a single galaxy provide the most robust way to explore the geometry of the Universe and to study the inner mass structure of lens galaxies. We present here a study of the SDSS J0100+1818 deflector, analyzing its total and baryonic mass distributions. The system comprises an ultra-massive early-type galaxy, surrounded by ten multiple images of three background sources. Exploiting high-resolution HST photometry and VLT/X-shooter spectroscopy we conduct a strong lensing analysis with the software GLEE to reconstruct the complex surface brightness distributions of the background sources over approximately 7200 HST pixels. These results are presented in our recent paper, Bolamperti et al. (2023). Finally, we present some preliminary results from new VLT/MUSE observations, that will allow us to build a new strong lensing+dynamics joint model and measure the values of the total matter density and of the cosmological constant parameters, Ωm and ΩΛ.

Distinguishing Dark Matter, Modified Gravity, and Modified Inertia with the Inner and Outer Parts of Galactic Rotation Curves

The missing gravity in galaxies requires dark matter, or alternatively a modification of gravity or inertia. These theoretical possibilities of fundamental importance may be distinguished by the statistical relation between the observed centripetal acceleration of particles in orbital motion and the expected Newtonian acceleration from the observed distribution of baryons in galaxies. Here predictions of cold dark matter halos, modified gravity, and modified inertia are compared and tested by a statistical sample of galaxy rotation curves from the Spitzer Photometry and Accurate Rotation Curves database. Modified gravity under an estimated mean external field correctly predicts the observed statistical relation of accelerations from both the inner and outer parts of rotation curves. Taken at face value there is a difference between the inner and outer parts on an acceleration plane which would be inconsistent with current proposals of modified inertia. Removing galaxies with possible systematic concerns such as central bulges or special inclinations does not change this trend. Cold dark matter halos predict a systematically deviating relation from the observed one. All aspects of rotation curves appear to be most naturally explained by modified gravity.

Testing modified gravity theories with numerical solutions of the external field effect in rotationally supported galaxies

The strong equivalence principle is violated by gravity theories of Milgromian dynamics (MOND) through the action of the external field effect. We test two different Lagrangian theories AQUAL and QUMOND based on their numerical solutions of the external field effect, by comparing two independent estimates of the mean external field strength of the nearby universe: a theory-deduced value from fitting the outer rotation curves of 114 galaxies and an empirical value from the large-scale distribution of cosmic baryons. The AQUAL-deduced external field strength from rotation curves agrees with that from the large-scale cosmic environment, while QUMOND-deduced value is somewhat higher. This suggests that AQUAL is likely to be preferred over QUMOND as an effective non-relativistic limit of a potential relativistic modified gravity theory.

Can ultralight dark matter explain the age–velocity dispersion relation of the Milky Way disc: A revised and improved treatment

ABSTRACT Ultralight axion-like particles ma ∼ 10−22 eV, or Fuzzy Dark Matter (FDM), behave comparably to cold dark matter (CDM) on cosmological scales and exhibit a kpc-size de Broglie wavelength capable of alleviating established (sub-)galactic-scale problems of CDM. Substructures inside an FDM halo incur gravitational potential perturbations, resulting in stellar heating sufficient to account for the Galactic disc thickening over a Hubble time, as first demonstrated by Church et al. We present a more sophisticated treatment that incorporates the full baryon and dark matter distributions of the Milky Way and adopts stellar disc kinematics inferred from recent Gaia, APOGEE, and LAMOST surveys. Ubiquitous density granulation and subhalo passages, respectively, drive inner disc thickening and flaring of the outer disc, resulting in an observationally consistent ‘U-shaped’ disc vertical velocity dispersion profile with the global minimum located near the solar radius. The observed age–velocity dispersion relation in the solar vicinity can be explained by the FDM-substructure-induced heating and places an exclusion bound ma ≳ 0.4 × 10−22 eV. We assess non-trivial uncertainties in the empirical core–halo relation, FDM subhalo mass function and tidal stripping, and stellar heating estimate. The mass range ma ≃ 0.5–0.7 × 10−22 eV favoured by the observed thick disc kinematics is in tension with several exclusion bounds inferred from dwarf density profiles, stellar streams, and Milky Way satellite populations, which could be significantly relaxed due to the aforesaid uncertainties. Additionally, strongly anisotropic heating could help explain the formation of ultra-thin disc galaxies.

Statistics for dark matter subhalo searches in gamma rays from a kinematically constrained population model: Fermi-LAT-like telescopes

Cold dark matter subhalos are expected to populate galaxies in numbers. If dark matter self-annihilates, these objects turn into prime targets for indirect searches, in particular with gamma-ray telescopes. Incidentally, the Fermi-LAT catalog already contains many unidentified sources that might be associated with subhalos. In this paper, we determine the probability for subhalos to be identified as gamma-ray pointlike sources from their predicted distribution properties. We use a semi-analytical model for the Galactic subhalo population, which, in contrast to cosmological simulations, can be made fully consistent with current kinematic constraints in the Milky Way and has no resolution limit. The model incorporates tidal stripping effects from a realistic distribution of baryons in the Milky Way. The same baryonic distribution contributes a diffuse gamma-ray foreground which adds up to that, often neglected in subhalo searches, generated by the smooth dark matter and the unresolved subhalos. This configuration implies a correlation between pointlike subhalo signals and diffuse background. Based on this semi-analytical modeling, we generate mock gamma-ray data assuming an idealized telescope resembling Fermi-LAT and perform a likelihood analysis to estimate the current and future sensitivity to subhalos in the relevant parameter space. We find a number of detectable subhalos of order ${\cal O}(<1)$ for optimistic model parameters and a WIMP mass of 100~GeV, maximized for a cored host halo. This barely provides support to the current interpretation of several Fermi unidentified sources as subhalos. We also find it more likely to detect the smooth Galactic halo itself before subhalos, should dark matter in the GeV-TeV mass range self-annihilate through $s$-wave processes.

Cotton gravity and 84 galaxy rotation curves

Recently, as a generalization of general relativity, a gravity theory has been proposed in which gravitational field equations are described by the Cotton tensor. That theory allows an additional contribution to the gravitational potential of a point mass that rises linearly with radius as $\Phi = -GM/r + \gamma r/2$, where $G$ is the Newton constant. The coefficients $M$ and $\gamma$ are the constants of integration and should be determined individually for each physical system. When applied to galaxies, the coefficient $\gamma$, which has the dimension of acceleration, should be determined for each galaxy. This is the same as having to determine the mass $M$ for each galaxy. If $\gamma$ is small enough, the linear potential term is negligible at short distances, but can become significant at large distances. In fact, it may contribute to the extragalactic systems. In this paper, we derive the effective field equation for Cotton gravity applicable to extragalactic systems. We then use the effective field equation to numerically compute the gravitational potential of a sample of 84 rotating galaxies. The 84 galaxies span a wide range, from stellar disk-dominated spirals to gas-dominated dwarf galaxies. We do not assume the radial density profile of the stellar disk, bulge, or gas; we use only the observed data. We find that the rotation curves of 84 galaxies can be explained by the observed distribution of baryons. This is due to the flexibility of Cotton gravity to allow the integration constant $\gamma$ for each galaxy. In the context of Cotton gravity, "dark matter" is in some sense automatically included as a curvature of spacetime. Consequently, even galaxies that have been assumed to be dominated by dark matter do not need dark matter.

A minimal model for massive neutrinos in Newtonian N-body simulations

We present a novel method for including the impact of massive neutrinos in cold dark matter N-body simulations. Our approach is compatible with widely employed Newtonian N-body codes and relies on only three simple modifications. First, we use commonly employed backscaling initial conditions, based on the cold dark matter plus baryon power spectrum instead of the total matter power spectrum. Second, the accurate Hubble rate is employed in both the backscaling and the evolution of particles in the N-body code. Finally, we shift the final particle positions in a post-processing step to account for the integrated effect of neutrinos on the particles in the simulation. However, we show that the first two modifications already capture most of the relevant neutrino physics for a large range of observationally interesting redshifts and scales. The output of the simulations are the cold dark matter and baryon distributions and can be analysed using standard methods. All modifications are simple to implement and do not generate any computational overhead. By implementing our methods in the N-body codes gadget-4 and gevolution, we show that any state-of-the-art Newtonian N-body code can be utilised out of the box. Our method is also compatible with higher order Lagrangian perturbation theory initial conditions and accurate for masses up to at least ∑ mν = 0.3 eV. Being formulated in relativistic gauge theory, in addition to including the impact of massive neutrinos, our method further includes relativistic corrections relevant on the large scales for free.

Kinematics and mass distributions for non-spherical deprojected Sérsic density profiles and applications to multi-component galactic systems

Using kinematics to decompose the mass profiles of galaxies, including the dark matter contribution, often requires parameterization of the baryonic mass distribution based on ancillary information. One such model choice is a deprojected Sérsic profile with an assumed intrinsic geometry. The case of flattened, deprojected Sérsic models has previously been applied to flattened bulges in local star-forming galaxies (SFGs), but can also be used to describe the thick, turbulent disks in distant SFGs. Here, we extend this previous work that derived density (ρ) and circular velocity (vcirc) curves by additionally calculating the spherically-enclosed 3D mass profiles (Msph). Using these profiles, we compared the projected and 3D mass distributions, quantified the differences between the projected and 3D half-mass radii (Re; r1/2, mass, 3D), and compiled virial coefficients relating vcirc(R) and Msph(&lt; r = R) or Mtot. We quantified the differences between mass fraction estimators for multi-component systems, particularly for dark matter fractions (ratio of squared circular velocities versus ratio of spherically enclosed masses), and we considered the compound effects of measuring dark matter fractions at the projected versus 3D half-mass radii. While the fraction estimators produce only minor differences, using different aperture radius definitions can strongly impact the inferred dark matter fraction. As pressure support is important in analyses of gas kinematics (particularly, at high redshifts), we also calculated the self-consistent pressure support correction profiles, which generally predict less pressure support than for the self-gravitating disk case. These results have implications for comparisons between simulation and observational measurements, as well as for the interpretation of SFG kinematics at high redshifts. We have made a set of precomputed tables and the code to calculate the profiles publicly available.

Incorporating baryon-driven contraction of dark matter halos in rotation curve fits

The condensation of baryons within a dark matter (DM) halo during galaxy formation should result in some contraction of the halo as the combined system settles into equilibrium. We quantify this effect on the cuspy primordial halos predicted by DM-only simulations for the baryon distributions observed in the galaxies of the SPARC database. We find that the DM halos of high surface brightness galaxies (with Σeff ≳ 100 L⊙ pc−2 at 3.6 μm) experience strong contraction. Halos become more cuspy as a result of compression: the inner DM density slope increases with the baryonic surface mass density. We iteratively fit rotation curves to find the balance between initial halo parameters (constrained by abundance matching), compression, and stellar mass-to-light ratio. The resulting fits often require lower stellar masses than expected for stellar populations, particularly in galaxies with bulges: stellar mass must be reduced to make room for the DM it compresses. This trade off between dark and luminous mass is reminiscent of the cusp-core problem in dwarf galaxies, but occurs in more massive systems: the present-epoch DM halos cannot follow from cuspy primordial halos unless (1) the stellar mass-to-light ratios are systematically smaller than expected from standard stellar population synthesis models, and/or (2) there is a net outward mass redistribution from the initial cusp, even in massive galaxies widely considered to be immune from such effects.

In-medium hadronization of heavy quarks and its effect on charmed meson and baryon distributions in heavy-ion collisions

We present a new model for the description of heavy-quark hadronization in relativistic heavy-ion collisions in the presence of a reservoir of lighter thermal particles with which recombination can occur leading to the formation of color-singlet clusters. Color neutralization is assumed to occur locally, within the same fluid cell occupied by the heavy quark and it proceeds via the recombination of the latter with light antiquarks (quarks) or diquarks (anti-diquarks), which are assumed to be present in the medium with thermal abundance as effective degrees of freedom around the QCD crossover temperature T_c. Typically the resulting color-singlet clusters have quite low invariant mass, in most of the cases below 4 GeV, and in this case they are assumed to undergo an isotropic 2-body decay in their local rest-frame. Heavier clusters are instead fragmented as Lund strings. The possibility of recombination with light diquarks enhances the yields of charmed baryons, in qualitative agreement with recent measurements. The assumption of local color neutralization leads to a strong space-momentum correlation, which provides a substantial enhancement of the collective flow of the final-state charmed hadrons, affecting both their momentum and their angular distributions.

Probing the shape of the Milky Way dark matter halo with hypervelocity stars: A new method

We propose a new method for determining the shape of the gravitational potential of the dark matter (DM) halo of the Milky Way (MW) with the galactocentric tangential velocities of a sample of hypervelocity stars (HVSs). We compute the trajectories of different samples of HVSs in a MW where the baryon distribution is axisymmetric and the DM potential either is spherical or is spheroidal or triaxial with radial-dependent axis ratios. We create ideal observed samples of HVSs with known latitudinal components of the tangential velocity, vϑ, and azimuthal component of the tangential velocity, vφ. We determine the shape of the DM potential with the distribution of |vϑ| when the Galactic potential is axisymmetric, or with the distribution of |vϑ| and of a function, $ \bar{v}_{\varphi} $, of vφ when the Galactic potential is non-axisymmetric. We recover the correct shape of the DM potential by comparing the distribution of |vϑ| and $ \bar{v}_{\varphi} $ of the ideal observed sample against the corresponding distributions of mock samples of HVSs that traveled in DM halos of different shapes. We use ideal observed samples of ∼800 HVSs, which are the largest samples of 4 M⊙ HVSs ejected with the Hills mechanism at a rate ∼10−4 yr−1, currently outgoing, and located at more than 10 kpc from the Galactic Center. In our ideal case of galactocentric velocities with null uncertainties and no observational limitations, the method recovers the correct shape of the DM potential with a success rate S ≳ 89% when the Galactic potential is axisymmetric, and S &gt; 96% in the explored non-axisymmetric cases. The unsuccessful cases yield axis ratios of the DM potential that are off by ±0.1. The success rate decreases with decreasing size of the HVS sample: for example, for a spherical DM halo, S drops from ∼98% to ∼38% when the sample size decreases from ∼800 to ∼40 HVSs. Accurate estimates of the success rate of our method applied to real data require more realistic samples of mock observed HVSs. Nevertheless, our analysis suggests that a robust determination of the shape of the DM potential requires the measure of the galactocentric velocity of a few hundred HVSs of robustly confirmed galactocentric origin.

Relative distribution of dark matter, gas, and stars around cosmic filaments in the IllustrisTNG simulation

We present a comprehensive study of the distribution of matter around different populations of large-scale cosmic filaments, using the IllustrisTNG simulation atz = 0. We computed the dark matter (DM), gas, and stellar radial density profiles of filaments, and we characterise the distribution of the baryon fraction in these structures. We find that baryons exactly follow the underlying DM distribution only down tor ∼ 7 Mpc to the filament spines. At shorter distances (r &lt; 7 Mpc), the baryon fraction profile of filaments departs from the cosmic value Ωb/Ωm. While in ther ∼ 0.7−7 Mpc radial domain this departure is due to the radial accretion of the warm-hot intergalactic medium (WHIM) towards the filament cores (creating an excess of baryons with respect to the cosmic fraction), the cores of filaments (r &lt; 0.7 Mpc) show a clear baryon depletion instead. The analysis of the efficiency of active galactic nuclei (AGN) feedback events in filaments reveals that they are potentially powerful enough to eject gas outside of the gravitational potential wells of filaments. We show that the large-scale environment (i.e. denser versus less dense, hotter versus colder regions) has a non-negligible effect on the absolute values of the DM, gas, and stellar densities around filaments. Nevertheless, the relative distribution of baryons with respect to the underlying DM density field is found to be independent of the filament population. Finally, we provide scaling relations between the gas density, temperature, and pressure for the different populations of cosmic filaments. We compare these relations to those pertaining to clusters of galaxies, and find that these cosmic structures occupy separate regions of the density-temperature and density-pressure planes.

Mapping the Circumgalactic Medium in Emission

Mapping the circumgalactic medium (CGM) in emission is a crucial next step to understand both the baryonic distribution of the universe and the multi-scale multi-phase medium that is a primary mediator of galaxy formation and evolution. Here we discuss the current state of the field and next steps.

Model of hybrid star with baryonic and strange quark matter in Tolman–Kuchowicz spacetime

The purpose of our work is to investigate some new features of a static anisotropic relativistic hybrid compact star composed of strange quark matter (SQM) in the inner core and normal baryonic matter distribution in the crust. Here we apply the simplest form of the phenomenological MIT bag model equation of state [Formula: see text] to correlate the density and pressure of strange quark matter within the stellar interior, whereas radial pressure and matter density due to baryonic matter are connected by the simple linear equation of state [Formula: see text]. In order to obtain the solution of the Einstein field equations, we have used the Tolman–Kuchowicz ansatz [R. C. Tolman, Static solutions of Einstein’s field equations for spheres of fluid, Phys. Rev. 55 (1939) 364–373; B. Kuchowicz, Acta Phys. Pol. 33 (1968) 541] and further derivation of the arbitrary constants from some physical conditions. Here, we examine our proposed model graphically and analytically in detail for physically plausible conditions. In particular, for this investigation, we have reported on the compact object Her [Formula: see text] [[Formula: see text]; [Formula: see text] km] in our paper as a strange quark star candidate. In order to check the physical validity and stability of our suggested model, we have performed various physical tests both analytically and graphically, namely, dynamical equilibrium of applied forces, energy conditions, compactness factor and surface redshift. Finally, we have found that our present model meets all the necessary physical requirements for a realistic model and can be studied for strange quark stars (SQS).