Search NASA⌕ Search

Engineering topics

Marinacci, Federico

Publications and source records attributed to Marinacci, Federico.

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Velocity-dependent annihilation radiation from dark matter subhalos in cosmological simulations

Here, we use the suite of Milky Way-like galaxies in the Auriga simulations to determine the contribution to annihilation radiation from dark matter subhalos in three velocity-dependent dark matter annihilation models: Sommerfeld, p-wave, and d-wave models. We compare these to the corresponding distribution in the velocity-independent s-wave annihilation model. For both the hydrodynamical and dark-matter-only simulations, only in the case of the Sommerfeld-enhanced annihilation does the total annihilation flux from subhalos exceed the total annihilation flux from the smooth halo component within the virial radius of the halo. Progressing from Sommerfeld to the s, p, and d-wave models, the contribution from the smooth component of the halo becomes more dominant, implying that for the p-wave and d-wave models the smooth component is by far the dominant contribution to the radiation. Comparing to the Galactic center excess observed by Fermi-LAT, for all simulated halos the emission is dominated by the smooth halo contribution. However, it is possible that for Sommerfeld models, extrapolation down to mass scales below the current resolution limit of the simulation would imply a non-negligible contribution to the gamma-ray emission from the Galactic Center region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Velocity-dependent J-factors for annihilation radiation from cosmological simulations

Abstract We determine the dark matter pair-wise relative velocity distribution in a set of Milky Way-like halos in the Auriga and APOSTLE simulations. Focusing on the smooth halo component, the relative velocity distribution is well-described by a Maxwell-Boltzmann distribution over nearly all radii in the halo. We explore the implications for velocity-dependent dark matter annihilation, focusing on four models which scale as different powers of the relative velocity: Sommerfeld, s-wave, p-wave, and d-wave models. We show that the J-factors scale as the moments of the relative velocity distribution, and that the halo-to-halo scatter is largest for d-wave, and smallest for Sommerfeld models. The J-factor is strongly correlated with the dark matter density in the halo, and is very weakly correlated with the velocity dispersion. This implies that if the dark matter density in the Milky Way can be robustly determined, one can accurately predict the dark matter annihilation signal, without the need to identify the dark matter velocity distribution in the Galaxy.

Astronomy & Astrophysics↗
DOE OSTI2019

First Star-Forming Structures in Fuzzy Cosmic Filaments

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter, we investigate how the first galaxies are assembled in a “fuzzy” dark matter (FDM) cosmology where dark matter is an ultralight ~10 –22 eV boson and the primordial stars are expected to form along dense dark matter filaments. Furthermore, using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical solitonlike cores, which are unstable under gravity and collapse into kiloparsec-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter—a smoking gun signature of FDM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Dust in and around galaxies: dust in cluster environments and its impact on gas cooling

Simulating the dust content of galaxies and their surrounding gas is challenging due to the wide range of physical processes affecting the dust evolution. Here we present cosmological hydrodynamical simulations of a cluster of galaxies, M 200,crit = 6×10 14 M ⊙ , including a novel dust model for the moving mesh code AREPO. This model includes dust production, growth, supernova-shock-driven destruction, ion-collision-driven thermal sputtering, and high-temperature dust cooling through far-infrared reradiation of collisionally deposited electron energies. Adopting a rather low thermal sputtering rate, we find, consistent with observations, a present-day overall dust-to-gas ratio of ~2 × 10 -5 , a total dust mass of ~2×10 9 M ⊙ , and a dust mass fraction of ~3 × 10 -6 . The typical thermal sputtering time-scales within ~100 kpc are around ~10 Myr, and increase towards the outer parts of the cluster to ~10 3 Myr at a cluster-centric distance of 1 Mpc. The condensation of gas-phase metals into dust grains reduces high-temperature metal-line cooling, but also leads to additional dust infrared cooling. The additional infrared cooling changes the overall cooling rate in the outer parts of the cluster, beyond ~1 Mpc, by factors of a few. This results in noticeable changes of the entropy, temperature, and density profiles of cluster gas once dust formation is included. The emitted dust infrared emission due to dust cooling is consistent with observational constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The velocity anisotropy of the Milky Way satellite system

We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, β, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE (A Project Of Simulating The Local Environment) and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from β ~-2 at r ~ 20 kpc to β ~ 0.5 at r ~ 200 kpc, indicating that satellites closer to the Galactic centre have tangentially biased motions while those farther out have radially biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the β( r ) profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the β( r ) profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the MW system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

AREPO-RT : radiation hydrodynamics on a moving mesh

We introduce AREPO-RT , a novel radiation hydrodynamic (RHD) solver for the unstructured moving-mesh code AREPO . Our method solves the moment-based radiative transfer equations using the M1 closure relation. We achieve second-order convergence by using a slope-limited linear spatial extrapolation and a first-order time prediction step to obtain the values of the primitive variables on both sides of the cell interface. A Harten–Lax–van Leer flux function, suitably modified for moving meshes, is then used to solve the Riemann problem at the interface. The implementation is fully conservative and compatible with the individual time-stepping scheme of AREPO . It incorporates atomic hydrogen (H) and helium (He) thermochemistry, which is used to couple the ultraviolet radiation field to the gas. Additionally, infrared (IR) radiation is coupled to the gas under the assumption of local thermodynamic equilibrium between the gas and the dust. We successfully apply our code to a large number of test problems, including applications such as the expansion of HII regions, radiation pressure-driven outflows, and the levitation of optically thick layer of gas by trapped IR radiation. The new implementation is suitable for studying various important astrophysical phenomena, such as the effect of radiative feedback in driving galactic scale outflows, radiation-driven dusty winds in high-redshift quasars, or simulating the reionization history of the Universe in a self-consistent manner.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Similar star formation rate and metallicity variability time-scales drive the fundamental metallicity relation

The fundamental metallicity relation (FMR) is a postulated correlation between galaxy stellar mass, star formation rate (SFR), and gas-phase metallicity. At its core, this relation posits that offsets from the mass–metallicity relation (MZR) at a fixed stellar mass are correlated with galactic SFR. In this Letter, we use hydrodynamical simulations to quantify the time-scales over which populations of galaxies oscillate about the average SFR and metallicity values at fixed stellar mass. We find that Illustris and IllustrisTNG predict that galaxy offsets from the star formation main sequence and MZR oscillate over similar time-scales, are often anticorrelated in their evolution, evolve with the halo dynamical time, and produce a pronounced FMR. Our models indicate that galaxies oscillate about equilibrium SFR and metallicity values – set by the galaxy’s stellar mass – and that SFR and metallicity offsets evolve in an anticorrelated fashion. This anticorrelated variability of the metallicity and SFR offsets drives the existence of the FMR in our models. In contrast to Illustris and IllustrisTNG, we speculate that the SFR and metallicity evolution tracks may become decoupled in galaxy formation models dominated by feedback-driven globally bursty SFR histories, which could weaken the FMR residual correlation strength. This opens the possibility of discriminating between bursty and non-bursty feedback models based on the strength and persistence of the FMR – especially at high redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Is There a Disk of Satellites around the Milky Way?

The “disk of satellites” (DoS) around the Milky Way is a highly debated topic with conflicting interpretations of observations and their theoretical models. We perform a comprehensive analysis of all of the dwarfs detected in the Milky Way and find that the DoS structure depends strongly on the plane identification method and the sample size. In particular, we demonstrate that a small sample size can artificially produce a highly anisotropic spatial distribution and a strong clustering of the angular momentum of the satellites. Moreover, we calculate the evolution of the 11 classical satellites with proper motion measurements and find that the thin DoS in which they currently reside is transient. Furthermore, we analyze two cosmological simulations using the same initial conditions of a Milky-Way-sized galaxy, an N -body run with dark matter only, and a hydrodynamic one with both baryonic and dark matter, and find that the hydrodynamic simulation produces a more anisotropic distribution of satellites than the N -body one. Our results suggest that an anisotropic distribution of satellites in galaxies can originate from baryonic processes in the hierarchical structure formation model, but the claimed highly flattened, coherently rotating DoS of the Milky Way may be biased by the small-number selection effect. These findings may help resolve the contradictory claims of DoS in galaxies and the discrepancy among numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Simulating the dust content of galaxies: successes and failures

We present full-volume cosmological simulations, using the moving-mesh code AREPO to study the coevolution of dust and galaxies. We extend the dust model in AREPO to include thermal sputtering of grains and investigate the evolution of the dust mass function, the cosmic distribution of dust beyond the interstellar medium and the dependence of dust-to-stellar mass ratio on galactic properties. The simulated dust mass function is well described by a Schechter fit and lies closest to observations at z = 0. The radial scaling of projected dust surface density out to distances of 10 Mpc around galaxies with magnitudes 17 < i < 21 is similar to that seen in Sloan Digital Sky Survey data, albeit with a lower normalization. At z = 0, the predicted dust density of Ωdust ≈ 1.3 × 10 -6 lies in the range of Ωdust values seen in low-redshift observations. We find that the dust-to-stellar mass ratio anticorrelates with stellar mass for galaxies living along the star formation main sequence. Moreover, we estimate the 850 μm number density functions for simulated galaxies and analyse the relation between dust-to-stellar flux and mass ratios at z = 0. At high redshift, our model fails to produce enough dust-rich galaxies, and this tension is not alleviated by adopting a top-heavy initial mass function. We do not capture a decline in Ωdust from z = 2 to 0, which suggests that dust production mechanisms more strongly dependent on star formation may help to produce the observed number of dusty galaxies near the peak of cosmic star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE ROLE OF COSMIC-RAY PRESSURE IN ACCELERATING GALACTIC OUTFLOWS

We study the formation of galactic outflows from supernova (SN) explosions with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SN placement and energy feedback, including cosmic rays (CRs), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗