Showing 1–50 of 159 results for author: Hoffman, Y

The large scale velocity field from the Cosmicflows-4 data

Authors: Yehuda Hoffman, Aurélien Valade, Noam I. Libeskind, Jenny G. Sorce, R. Brent Tully, Simon Pfeifer, Stefan Gottlöber, Daniel Pomaréde

Abstract: The reconstruction of the large scale velocity field from the grouped Cosmicflows-4 (CF4) database is presented. The lognormal bias of the inferred distances and velocities data is corrected by the Bias Gaussianization correction (BGc) scheme, and the linear density and velocity fields are reconstructed by means of the Wiener filter (WF) and constrained realizations (CRs) algorithm. These tools ar… ▽ More The reconstruction of the large scale velocity field from the grouped Cosmicflows-4 (CF4) database is presented. The lognormal bias of the inferred distances and velocities data is corrected by the Bias Gaussianization correction (BGc) scheme, and the linear density and velocity fields are reconstructed by means of the Wiener filter (WF) and constrained realizations (CRs) algorithm. These tools are tested against a suite of random and constrained Cosmicflows-3-like mock data. The CF4 data consists of 3 main subsamples - the 6dFGS and the SDSS data - and the `others'. The individual contributions of the subsamples have been studied. The quantitative analysis of the velocity field is done mostly by the mean overdensity ($Δ_L(R)$) and the bulk velocity ($V_{\mathrm{bulk}}(R)$) profiles of the velocity field out to $300\, h^{-1}{\rm Mpc}$. The $V_{\mathrm{bulk}}(R)$ and $Δ_{\mathrm L}(R)$ profiles of the CF4 data without its 6dFGS component are consistent with the cosmic variance to within $1σ$. The 6dFGS sample dominates the $V_{\mathrm{bulk}}$ ($Δ_{\mathrm L}$) profile beyond $\sim120\, h^{-1}{\rm Mpc}$, and drives it to roughly a $3.4σ$ ($-1.9σ$) excess (deficiency) relative to the cosmic variance at $R\sim250\ (190)\ \, h^{-1}{\rm Mpc}$. The excess in the amplitude of $V_{\mathrm{bulk}}$ is dominated by its Supergalactic X component, roughly in the direction of the Shapley Concentration. The amplitude and alignment of the inferred velocity field from the CF4 data is at $\sim(2\,-\,3)\,σ$ discrepancy with respect to the $Λ$CDM model. Namely, it is somewhat atypical but yet there is no compelling tension with the model. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 17 pages, 14 figures, accepted for publication to MNRAS

Exploring the centre of mass properties of LG-like galaxies

Authors: Jean-Baptiste Salomon, Noam Libeskind, Yehuda Hoffman

Abstract: From high resolution cosmological simulations of the Local Group in realistic environment, namely HESTIA simulations, we study the position and kinematic deviations that may arise between the disc of a Milky Way (or Andromeda)-like galaxy and its halo. We focus on the 3-dimensional analysis of the centres of mass (COM). The study presents two parts. We first consider individual particles to track… ▽ More From high resolution cosmological simulations of the Local Group in realistic environment, namely HESTIA simulations, we study the position and kinematic deviations that may arise between the disc of a Milky Way (or Andromeda)-like galaxy and its halo. We focus on the 3-dimensional analysis of the centres of mass (COM). The study presents two parts. We first consider individual particles to track down the very nature and amplitude of the physical deviations of the COM with respect to the distance from the disc centre. Dark matter dominates the behaviour of the COM of all particles at all distances. But the total COM is also very close to the COM of stars. In the absence of a significant merger, the velocity offsets are marginal (10 km/s) but the positional shifts can be important compared to the disc characteristics (>10 kpc). In the event of a massive accretion, discrepancies are of the same order as the recent finding for the MW under the Magellanic Clouds influence. In a second part, the accent is put on the study of various populations of subhaloes and satellites. We show that satellites properly represent the entire subhalo population. There exists strong mismatch in phase space between the satellites' COM and the host disc. Moreover, the results are highly inhomogeneous between the simulations, and thus between the accretion histories. Finally, we point out that these shifts are mainly due to a few of the most massive objects. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS

A Local Universe model for constrained simulations

Authors: Simon Pfeifer, Aurélien Valade, Stefan Gottlöber, Yehuda Hoffman, Noam I. Libeskind, Wojciech A. Hellwing

Abstract: The aim of cosmological simulations is to reproduce the properties of the observed Universe, serving as tools to test structure and galaxy formation models. Constrained simulations of our local cosmological region up to a few hundred Mpc/h , the local Universe, are designed to reproduce the actual cosmic web of structures as observed. A question that often arises is how to judge the quality of con… ▽ More The aim of cosmological simulations is to reproduce the properties of the observed Universe, serving as tools to test structure and galaxy formation models. Constrained simulations of our local cosmological region up to a few hundred Mpc/h , the local Universe, are designed to reproduce the actual cosmic web of structures as observed. A question that often arises is how to judge the quality of constrained simulations against the observations of the Local Universe. Here we introduce the Local Universe model (LUM), a new methodology, whereby many constrained simulations can be judged and the ''best'' initial conditions can be identified. By characterising the Local Universe as a set of rich clusters, the model identifies haloes that serve as simulated counterparts to the observed clusters. Their merit is determined against a null hypothesis, the probability that such a counterpart could be identified in a random, unconstrained simulation. This model is applied to 100 constrained simulations using the Cosmicflows-3 data. Cluster counterparts are found for all constrained simulations, their distribution of separation from the true observed cluster position and their mass distribution are investigated. Lastly, the ''best'' constrained simulation is selected using the LUM and discussed in more detail. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 10 pages, 8 figures, submitted to MNRAS

The undiscovered ultra-diffuse galaxies of the Local Group

Authors: Oliver Newton, Arianna Di Cintio, Salvador Cardona-Barrero, Noam I. Libeskind, Yehuda Hoffman, Alexander Knebe, Jenny Sorce, Matthias Steinmetz, Elmo Tempel

Abstract: Ultra-diffuse galaxies (UDGs) are attractive candidates to probe cosmological models and test theories of galaxy formation at low masses; however, they are difficult to detect because of their low surface brightness. In the Local Group a handful of UDGs have been found to date, most of which are satellites of the Milky Way and M31, and only two are isolated galaxies. It is unclear whether so few U… ▽ More Ultra-diffuse galaxies (UDGs) are attractive candidates to probe cosmological models and test theories of galaxy formation at low masses; however, they are difficult to detect because of their low surface brightness. In the Local Group a handful of UDGs have been found to date, most of which are satellites of the Milky Way and M31, and only two are isolated galaxies. It is unclear whether so few UDGs are expected. We address this by studying the population of UDGs formed in hydrodynamic constrained simulations of the Local Group from the HESTIA suite. For a Local Group with a total enclosed mass $M_{\rm LG}\!\left(< 2.5\, {\rm Mpc}\right)=8\times10^{12}\, {\rm M_\odot}$, we predict that there are $12\pm3$ isolated UDGs (68% confidence) with stellar masses $10^6 \leq M_\ast\, /\, {\rm M_\odot} < 10^9,$ and effective radii $R_{\rm e} \geq 1.5\, {\rm kpc}$, within 2.5 Mpc of the Local Group, of which $2^{+2}_{-1}$ (68% confidence) are detectable in the footprint of the Sloan Digital Sky Survey (SDSS). Accounting for survey incompleteness, we find that almost the entire population of UDGs in the Local Group field would be observable in a future all-sky survey with a depth similar to the SDSS, the Dark Energy Survey, or the Legacy Survey of Space and Time. Our results suggest that there is a population of UDGs in the Local Group awaiting discovery. △ Less

Submitted 28 March, 2023; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 9 pages including acknowledgements and references. Contains 5 figures and two tables. A repository of reduced data and scripts to reproduce the results in this manuscript is available at https://github.com/Musical-Neutron/lg_field_udgs and archived at https://doi.org/10.5281/zenodo.7757685 . Accepted by ApJL

Constrained simulations of the local Universe with Modified Gravity

Authors: Krishna Naidoo, Wojciech A. Hellwing, Maciej Bilicki, Noam Libeskind, Simon Pfeifer, Yehuda Hoffman

Abstract: We present a methodology for constructing modified gravity (MG) constrained simulations of the local Universe using positions and peculiar velocities from the CosmicFlows data set. Our analysis focuses on the following MG models: the normal branch of the Dvali-Gabadadze-Porrati (nDGP) model and Hu-Sawicki $f(R)$ model. We develop a model independent methodology for constructing constrained simulat… ▽ More We present a methodology for constructing modified gravity (MG) constrained simulations of the local Universe using positions and peculiar velocities from the CosmicFlows data set. Our analysis focuses on the following MG models: the normal branch of the Dvali-Gabadadze-Porrati (nDGP) model and Hu-Sawicki $f(R)$ model. We develop a model independent methodology for constructing constrained simulations with any given power spectra and numerically calculated linear growth functions. Initial conditions (ICs) for a set of constrained simulations are constructed for the standard cosmological model $Λ$CDM and the MG models. Differences between the model's reconstructed Wiener filtered density and the resultant simulation density are presented showing the importance for the generation of MG constrained ICs to study the subtle effects of MG in the local Universe. These are the first MG constrained simulations ever produced. The current work paves the way to improved approximate methods for models with scale-dependent growth functions, such as $f(R)$, and for high-resolution hydrodynamical MG zoom-in simulations of the local Universe. △ Less

Submitted 24 February, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 10 pages (+references), 6 figures, accepted in Phys. Rev. D

Testing Bayesian reconstruction methods from peculiar velocities

Authors: Aurélien Valade, Noam I Libeskind, Yehuda Hoffman, Simon Pfeifer

Abstract: Reconstructing the large scale density and velocity fields from surveys of galaxy distances, is a major challenge for cosmography. The data is very noisy and sparse. Estimated distances, and thereby peculiar velocities, are strongly affected by the Malmquist-like lognormal bias. Two algorithms have been recently introduced to perform reconstructions from such data: the Bias Gaussian correction cou… ▽ More Reconstructing the large scale density and velocity fields from surveys of galaxy distances, is a major challenge for cosmography. The data is very noisy and sparse. Estimated distances, and thereby peculiar velocities, are strongly affected by the Malmquist-like lognormal bias. Two algorithms have been recently introduced to perform reconstructions from such data: the Bias Gaussian correction coupled with the Wiener filter (BGc/WF) and the HAMLET implementation of the Hamiltonian Monte Carlo forward modelling. The two methods are tested here against mock catalogs that mimic the Cosmicflows-3 data. Specifically the reconstructed cosmography and moments of the velocity field (monopole, dipole) are examined. A comparison is made to the ``exact'' wiener filter as well - namely the Wiener Filter in the unrealistic case of zero observational errors. This is to understand the limits of the WF method. The following is found. In the nearby regime ($d \lesssim 40 {\rm Mpc}/h$) the two methods perform roughly equally well. HAMLET does slightly better in the intermediate regime ($ 40 \lesssim d \lesssim 120 {\rm Mpc}/h$). The main differences between the two appear in the most distant regime ($d \gtrsim 120 {\rm Mpc}/h$), close to the edge of the data. The HAMLET outperforms the BGc/WF in terms of better and tighter correlations, yet in the distant regime the HAMLET yields a somewhat biased reconstruction. Such biases are missing from the BGc/WF reconstruction. In sum, both methods perform well and create reliable reconstructions with significant differences apparent when details are examined. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 13 pages, 12 figures

Anisotropic satellite accretion onto the Local Group with HESTIA

Authors: Alexandra Dupuy, Noam I. Libeskind, Yehuda Hoffman, Hélène M. Courtois, Stefan Gottlöber, Robert J. J. Grand, Alexander Knebe, Jenny G. Sorce, Elmo Tempel, R. Brent Tully, Mark Vogelsberger, Peng Wang

Abstract: How the cosmic web feeds halos, and fuels galaxy formation is an open question with wide implications. This study explores the mass assembly in the Local Group within the context of the local cosmography by employing simulations whose initial conditions have been constrained to reproduce the local environment. The goal of this study is to inspect whether the direction of accretion of satellites on… ▽ More How the cosmic web feeds halos, and fuels galaxy formation is an open question with wide implications. This study explores the mass assembly in the Local Group within the context of the local cosmography by employing simulations whose initial conditions have been constrained to reproduce the local environment. The goal of this study is to inspect whether the direction of accretion of satellites on to the Milky Way and Andromeda galaxies, is related to the cosmic web. The analysis considers the three high-resolution simulations available in the HESTIA simulation suite, as well as the derived velocity shear and tidal tensors. We notice two eras in the Local Group accretion history, delimited by an epoch around $z \approx 0.7$. We also find that satellites can travel up to $\sim 4$ Mpc, relative to their parent halo before crossing its viral radius $R_{200}$. Finally, we observe a strong alignment of the infall direction with the axis of slowest collapse $\vec{e_3}$ of both tidal and shear tensors, implying satellites of the Local Group originated from one particular region of the cosmic web and were channeled towards us via the process of accretion.This alignment is dominated by the satellites that enter during the early infall era, i.e $z>0.7$. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 10 pages, 11 figures, accepted for publication in MNRAS

The stellar halo in Local Group Hestia simulations III. Chemical abundance relations for accreted and in-situ stars

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We fou… ▽ More Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that the merger debris are chemically distinct from the survived dwarf galaxies. The mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal abundance gradients in the ELz, where the most metal-rich stars have formed in the inner parts of the disrupted systems before the merger and mainly contribute to the central regions of the hosts. Therefore, we suggest that abundance measurements in the inner MW will allow constraining better the parameters of building blocks of the MW stellar halo. The MDFs of the individual debris show several peaks and the majority of debris have lower metallicity than the in-situ stars for Lz>0, while non-rotating and retrograde accreted stars are similar to the in-situ. Prograde accreted stars show a prominent knee in the [Fe/H]-[Mg/Fe] plane while the retrograde stars typically deposit to a high-[Mg/Fe] sequence. We found that the metal-poor stars ([Fe/H]<-1) of the HESTIA galaxies exhibit between zero to 80 km/s net rotation which is consistent with the Aurora population. At higher metallicities, we detect a sharp transition (spin-up) from the turbulent phase to a disk-like rotation. Mergers debris are similar in the [Fe/H]-[Mg/Fe] plane. However, combining a set of abundances allows to capture chemical patterns corresponding to different debris, which are the most prominent as a function of stellar age. △ Less

Submitted 12 September, 2023; v1 submitted 11 June, 2022; originally announced June 2022.

Comments: Paper III in series of III, 21 pages, 17 figures, A&A in press

Journal ref: A&A 677, A91 (2023)

The stellar halo in Local Group Hestia simulations II. The accreted component

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Misha Haywood, Paola Di Matteo, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological… ▽ More In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience a few dozen mergers but only 1-4 mergers have the stellar mass ratio >0.2 where, depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo. Individual merger remnants show diverse density distributions at z=0, significantly overlapping with each other and with the in-situ stars in the ELz, UV and RVphi coordinates. The mergers debris often change their position in the ELz with time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinct ELz features. In the UV plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features. We found an age gradient in Elz space for the individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host. The bulk of these stars is being formed during the last stages of accretion, making it possible to date the merger. In actions space (Jr, Jz, Jφ), the mergers debris do not appear as isolated substructures but are instead scattered over a large parameters area and overlapping with the in-situ stars. We also introduce a purely kinematic space (Jz/Jr-eccentricity), where different merger debris can be disentangled better from each other and from the in-situ stars. △ Less

Submitted 12 September, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Paper II in a series of III, 22 pages, 16 figures, A&A in press

Journal ref: A&A 677, A90 (2023)

The stellar halo in Local Group Hestia simulations I. The in-situ component and the effect of mergers

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Misha Haywood, Paola Di Matteo, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the… ▽ More Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7-11Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vphi of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30-40%, while in the outer parts it typically does not exceed ~5% beyond 15 kpc. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax-Zmax space are mainly populated by the stars born between significant mergers. △ Less

Submitted 12 September, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Paper I in a series of III, 17 pages, 15 figures, A&A in press

Journal ref: A&A 677, A89 (2023)

Estimation of the masses in the Local Group by Gradient Boosted Decision Trees

Authors: Edoardo Carlesi, Yehuda Hoffman, Noam Libeskind

Abstract: Our goal is to estimate the mass of the Local Group (LG) and the individual masses of its primary galaxies, the M31 and the Milky Way (MW). We do this by means of a supervised machine learning algorithm, the gradient boosted decision trees (GBDT) and using the observed distance and relative velocity of the two as input parameters. The GBDT is applied to a sample of 2148 mock LGs drawn from a set o… ▽ More Our goal is to estimate the mass of the Local Group (LG) and the individual masses of its primary galaxies, the M31 and the Milky Way (MW). We do this by means of a supervised machine learning algorithm, the gradient boosted decision trees (GBDT) and using the observed distance and relative velocity of the two as input parameters. The GBDT is applied to a sample of 2148 mock LGs drawn from a set of 5 dark matter (DM)-only simulations, ran within the standard $Λ$CDM\ cosmological model. The selection of the mock LGs is guided by a LG model, which defines such objects. The role of the observational uncertainties of the input parameters is gauged by applying the model to an ensemble of mock LGs pairs whose observables are these input parameters perturbed by their corresponding observational errors. Finally the observational data of the actual LG is used to infer its relevant masses. Our main results are the sum and the individual masses of the MW and M31: $M_{tot} = 3.31 ^{+0.79}_{-0.67} $, $M_{MW}=1.15^{+0.25}_{-0.22}$ and $M_{M31}=2.01^{+0.65}_{-0.39} \ \ \times 10^{12}M_{\odot}$ (corresponding to the median and the 1st and 3rd quartiles). The ratio of the masses is $M_{M31}/M_{MW}=1.75^{+0.54}_{-0.28}$, where by convention the M31 is defined here to be the more massive of the two halos. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS

Hamiltonian Monte Carlo reconstruction from peculiar velocities

Authors: Aurélien Valade, Yehuda Hoffman, Noam I Libeskind, Romain Graziani

Abstract: The problem of the reconstruction of the large scale density and velocity fields from peculiar velocities surveys is addressed here within a Bayesian framework by means of Hamiltonian Monte Carlo (HMC) sampling. The HAmiltonian Monte carlo reconstruction of the Local EnvironmenT (Hamlet) algorithm is designed to reconstruct the linear large scale density and velocity fields in conjunction with the… ▽ More The problem of the reconstruction of the large scale density and velocity fields from peculiar velocities surveys is addressed here within a Bayesian framework by means of Hamiltonian Monte Carlo (HMC) sampling. The HAmiltonian Monte carlo reconstruction of the Local EnvironmenT (Hamlet) algorithm is designed to reconstruct the linear large scale density and velocity fields in conjunction with the undoing of lognormal bias in the derived distances and velocities of peculiar velocities surveys such as the Cosmicflows data. The Hamlet code has been tested against Cosmicflows mock catalogs consisting of up to 30 000 data points with mock errors akin to those of the Cosmicflows-3 data, within the framework of the LCDM standard model of cosmology. The Hamlet code outperforms previous applications of Gibbs sampling MCMC reconstruction from the Cosmicflows-3 data by two to four orders of magnitude in CPU time. The gain in performance is due to the inherent higher efficiency of the HMC algorithm and due to parallel computing on GPUs rather than CPUs. This gain will enable an increase in the reconstruction of the large scale structure from the upcoming Cosmicfows-4 data and the setting of constrained initial conditions for cosmological high resolution simulations. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 14 pages, 10 figures

Cold and Hot gas distribution around the Milky-Way-M31 system in the HESTIA simulations

Authors: Mitali Damle, Martin Sparre, Philipp Richter, Maan H. Hani, Sebastián E. Nuza, Christoph Pfrommer, Robert J. J. Grand, Yehuda Hoffman, Noam Libeskind, Jenny G. Sorce, Matthias Steinmetz, Elmo Tempel, Mark Vogelsberger, Peng Wang

Abstract: Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three hig… ▽ More Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three high-resolution {\sc arepo}-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the $z = 0$ simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (HI and SiIII) are more clumpy in comparison to warmer gas tracers (OVI, OVII and OVIII). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. HESTIA under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way. △ Less

Submitted 16 March, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 22 pages, 9 figures, revised version, Accepted in MNRAS

COWS: A filament finder for Hessian cosmic web identifiers

Authors: Simon Pfeifer, Noam I. Libeskind, Yehuda Hoffman, Wojciech A. Hellwing, Maciej Bilicki, Krishna Naidoo

Abstract: The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells… ▽ More The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells into an ensemble of individual discreet filaments. Specifically, a thinning algorithm is applied to velocity shear tensor based cosmic web (V-web) to identify the spine of the filaments. This results in a set of filaments with well defined end-point and length. It is confirmed that these sit at local density ridges and align with the appropriate direction defined by the underlying velocity field. The radial density profile of these curved cylindrical filaments, as well as the distribution of their lengths is also examined. The robustness of all results is checked against resolution and the V-web threshold. The code for the COWS method has been made publicly available. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 10 pages, 10 figures, submitted to MNRAS

NIHAO-LG: The uniqueness of Local Group dwarf galaxies

Authors: Nikhil Arora, Andrea V. Macció, Stéphane Courteau, Tobias Buck, Noam I. Libeskind, Jenny G. Sorce, Chris B. Brook, Yehuda Hoffman, Gustavo Yepes, Eduardo Carlesi, Connor Stone

Abstract: Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulat… ▽ More Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulations, performed with NIHAO-like hydrodynamics which track the evolution of the Milky Way (MW) and M31 along with ~100 dwarfs in the LG, reveal the total gas mass and stellar properties (velocity dispersion, evolution history, etc.) of present-day LG dwarfs to be similar to field systems. However, relative to field galaxies, LG dwarfs have more cold gas in their central parts and more metal-rich gas in the halo stemming from interactions with other dwarfs living in a high-density environment like the LG. Interestingly, the direct impact of massive MW/M31 analogues on the metallicity evolution of LG dwarfs is minimal; LG dwarfs accrete high-metallicity gas mostly from other dwarfs at late times. We have also tested for the impact of metal diffusion on the chemical evolution of LG dwarfs, and found that it does not affect the stellar or gaseous content of LG dwarfs. Our simulations suggest that the stellar components of LG dwarfs offer a unique and unbiased local laboratory for galaxy-formation tests and comparisons, especially against the overall dwarf population in the Universe. △ Less

Submitted 14 April, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: 16 pages, 12 figures, Accepted for publication in MNRAS

From Cosmicflows distance moduli to unbiased distances and peculiar velocities

Authors: Yehuda Hoffman, Adi Nusser, Aurelien Valade, Noam I. Libeskind, R. Brent Tully

Abstract: Surveys of galaxy distances and radial peculiar velocities can be used to reconstruct the large scale structure. Other than systematic errors in the zero-point calibration of the galaxy distances the main source of uncertainties of such data are errors on the distance moduli, assumed here to be Gaussian and thus turn into lognormal errors on distances and velocities. Naively treated, it leads to s… ▽ More Surveys of galaxy distances and radial peculiar velocities can be used to reconstruct the large scale structure. Other than systematic errors in the zero-point calibration of the galaxy distances the main source of uncertainties of such data are errors on the distance moduli, assumed here to be Gaussian and thus turn into lognormal errors on distances and velocities. Naively treated, it leads to spurious nearby outflow and strong infall at larger distances. The lognormal bias is corrected here and tested against mock data extracted from a $Λ$CDM simulation, designed to statistically follow the grouped Cosmicflows-3 (CF3) data. Considering a subsample of data points, all of which have the same true distances or same redshifts, the lognormal bias arises because the means of the distributions of observed distances and velocities are skewed off the means of the true distances and velocities. Yet, the medians are invariant under the lognormal transformation. That invariance allows the Gaussianization of the distances and velocities and the removal of the lognormal bias. This Bias Gaussianization correction (BGc) algorithm is tested against mock CF3 catalogs. The test consists of a comparison of the BGC estimated with the simulated distances and velocities and of an examination of the Wiener filter reconstruction from the BGc data. Indeed, the BGc eliminates the lognormal bias. The estimation of Hubble's ($H_{0}$) constant is also tested. The residual of the BGc estimated $H_{0}$ from the simulated values is $0.6 \pm 0.7 {\rm kms}^{-1}{\rm Mpc}^{-1}$ and is dominated by the cosmic variance. The BGc correction of the actual CF3 data yields $H_{0} = 75.8 \pm 1.1 {\rm kms}^{-1}{\rm Mpc}^{-1}$ . △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 14 pages, 16 figures, accepted for publication in MNRAS

Hermeian haloes: Field haloes that interacted with both the Milky Way and M31

Authors: Oliver Newton, Noam I. Libeskind, Alexander Knebe, Miguel A. Sánchez-Conde, Jenny G. Sorce, Sergey Pilipenko, Matthias Steinmetz, Ruediger Pakmor, Elmo Tempel, Yehuda Hoffman, Mark Vogelsberger

Abstract: The Local Group is a unique environment in which to study the astrophysics of galaxy formation. The proximity of the Milky Way and M31 enhances the frequency of interactions of the low-mass halo population with more massive dark matter haloes, which increases their concentrations and strips them of gas and other material. Some low-mass haloes pass through the haloes of the Milky Way or M31 and are… ▽ More The Local Group is a unique environment in which to study the astrophysics of galaxy formation. The proximity of the Milky Way and M31 enhances the frequency of interactions of the low-mass halo population with more massive dark matter haloes, which increases their concentrations and strips them of gas and other material. Some low-mass haloes pass through the haloes of the Milky Way or M31 and are either ejected into the field or exchanged between the two primary hosts. We use high resolution gas-dynamical simulations to describe a new class of field haloes that passed through the haloes of both the Milky Way and M31 at early times and are almost twice as concentrated as field haloes that do not interact with the primary pair. These 'Hermeian' haloes are distributed anisotropically at larger distances from the Local Group barycentre than the primary haloes and appear to cluster along the line connecting the Milky Way and M31. Hermeian haloes facilitate the exchange of dark matter, gas, and stars between the Milky Way and M31 and can enhance the star formation rate of the gas in the primary haloes during their interactions with them. We also show that some Hermeian haloes can host galaxies that, because they are embedded in haloes that are more concentrated than regular field haloes, are promising targets for indirect dark matter searches beyond the Milky Way virial radius and can produce signals that are competitive with those of some dwarf galaxies. Hermeian galaxies in the Local Group should be detectable by forthcoming wide-field imaging surveys. △ Less

Submitted 5 July, 2022; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 14 pages, 7 figures, 1 appendix. Accepted by MNRAS. Plotting scripts and data available from: http://doi.org/10.5281/zenodo.6629724

Journal ref: MNRAS 514 (2022) 3612-3625

The proper motion of Andromeda from Gaia eDR3: confirming a nearly radial orbit

Authors: J. -B. Salomon, R. Ibata, C. Reylé, B. Famaey, N. I. Libeskind, A. W. McConnachie, Y. Hoffman

Abstract: We present an analysis of the proper motion of the Andromeda galaxy (M31), based on the Early Third Data Release of the Gaia mission. We use the Gaia photometry to select young blue main sequence stars, and apply several quality cuts to obtain clean samples of these tracers. After correcting the proper motion measurements for the internal rotation of the M31 disk motion, we derive an apparent moti… ▽ More We present an analysis of the proper motion of the Andromeda galaxy (M31), based on the Early Third Data Release of the Gaia mission. We use the Gaia photometry to select young blue main sequence stars, and apply several quality cuts to obtain clean samples of these tracers. After correcting the proper motion measurements for the internal rotation of the M31 disk motion, we derive an apparent motion of 52.5 +/- 5.8 muas/yr with respect to the Gaia reference frame, or 61.9 +/- 9.7 muas/yr after applying a zero-point correction determined from quasars within 20 degrees from M31 and a correction from systemic biases. Accounting for the Solar reflex motion we deduce a relative velocity between Andromeda and the Milky way (in a non-rotating frame at the current location of the Sun) of 42.2 +/- 39.3 km/s along right ascension (40.0 +/- 39.3 km/s along galactic longitude) and -59.4 +/- 30.3 km/s along declination (-60.9 +/- 30.3 km/s along galactic latitude), with a total transverse velocity of V_trans = 82.4 +/- 31.2 km/s. These values are consistent with (but more accurate than) earlier Hubble Space Telescope measurements that predict a future merger between the two galaxies. We also note a surprisingly large difference in the derived proper motion between the blue stars in M31 and samples of red stars that appear to lie in that galaxy. We propose several hypotheses to explain the discrepancy but found no clear evidence with the current data to privilege any one of them. △ Less

Submitted 3 September, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

Comments: 10 pages, 7 figures, Accepted for publication in MNRAS

The sum of the masses of the Milky Way and M31: a likelihood-free inference approach

Authors: Pablo Lemos, Niall Jeffrey, Lorne Whiteway, Ofer Lahav, Niam I Libeskind, Yehuda Hoffman

Abstract: We use Density Estimation Likelihood-Free Inference, $Λ$ Cold Dark Matter simulations of $\sim 2M$ galaxy pairs, and data from Gaia and the Hubble Space Telescope to infer the sum of the masses of the Milky Way and Andromeda (M31) galaxies, the two main components of the Local Group. This method overcomes most of the approximations of the traditional timing argument, makes the writing of a theoret… ▽ More We use Density Estimation Likelihood-Free Inference, $Λ$ Cold Dark Matter simulations of $\sim 2M$ galaxy pairs, and data from Gaia and the Hubble Space Telescope to infer the sum of the masses of the Milky Way and Andromeda (M31) galaxies, the two main components of the Local Group. This method overcomes most of the approximations of the traditional timing argument, makes the writing of a theoretical likelihood unnecessary, and allows the non-linear modelling of observational errors that take into account correlations in the data and non-Gaussian distributions. We obtain an $M_{200}$ mass estimate $M_{\rm MW+M31} = 4.6^{+2.3}_{-1.8} \times 10^{12} M_{\odot}$ ($68 \%$ C.L.), in agreement with previous estimates both for the sum of the two masses and for the individual masses. This result is not only one of the most reliable estimates of the sum of the two masses to date, but is also an illustration of likelihood-free inference in a problem with only one parameter and only three data points. △ Less

Submitted 21 December, 2020; v1 submitted 16 October, 2020; originally announced October 2020.

Comments: 14 pages, 9 figures. Accepted by Physical Review D

Journal ref: Phys. Rev. D 103, 023009 (2021)

The Hestia project: simulations of the Local Group

Authors: Noam I. Libeskind, Edoardo Carlesi, Rob J. J. Grand, Arman Khalatyan, Alexander Knebe, Ruediger Pakmor, Sergey Pilipenko, Marcel S. Pawlowski, Martin Sparre, Elmo Tempel, Peng Wang, Helene M. Courtois, Stefan Gottloeber, Yehuda Hoffman, Ivan Minchev, Christoph Pfrommer, Jenny G. Sorce, Volker Springel, Matthias Steinmetz, R. Brent Tully, Mark Vogelsberger, Gustavo Yepes

Abstract: We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter on… ▽ More We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: 20 pages, 13 figures, submitted to MNRAS

Cosmicflows-3: The South Pole Wall

Authors: Daniel Pomarede, R. Brent Tully, Romain Graziani, Helene M. Courtois, Y. Hoffman, Jeremy Lezmy

Abstract: Velocity and density field reconstructions of the volume of the universe within 0.05c derived from the Cosmicflows-3 catalog of galaxy distances has revealed the presence of a filamentary structure extending across ~ 0.11c. The structure, at a characteristic redshift of 12,000 km/s, has a density peak coincident with the celestial South Pole. This structure, the largest contiguous feature in the l… ▽ More Velocity and density field reconstructions of the volume of the universe within 0.05c derived from the Cosmicflows-3 catalog of galaxy distances has revealed the presence of a filamentary structure extending across ~ 0.11c. The structure, at a characteristic redshift of 12,000 km/s, has a density peak coincident with the celestial South Pole. This structure, the largest contiguous feature in the local volume and comparable to the Sloan Great Wall at half the distance, is given the name the South Pole Wall. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Accepted for publication by The Astrophysical Journal. Links to a video and 4 interactive Sketchfab visualizations. High quality figures, video, and interactive visualizations are available at http://irfu.cea.fr/southpolewall

Cosmicflows-3: Two Distance$-$Velocity Calculators

Authors: Ehsan Kourkchi, Helene M. Courtois, Romain Graziani, Yehuda Hoffman, Daniel Pomarede, Edward J. Shaya, R. Brent Tully

Abstract: Tools are provided at the Extragalactic Distance Database website that provide relationships between the distances and velocities of galaxies based on smoothed versions of the velocity fields derived by the Cosmicflows program. Tools are provided at the Extragalactic Distance Database website that provide relationships between the distances and velocities of galaxies based on smoothed versions of the velocity fields derived by the Cosmicflows program. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Comments: 10 pages, 3 figures, accepted for publication in The Astronomical Journal

Journal ref: 2020, AJ, 159, 67

On the Mass Assembly History of the Local Group

Authors: Edoardo Carlesi, Yehuda Hoffman, Stefan Gottlöber, Noam I. Libeskind, Alexander Knebe, Gustavo Yepes, Sergey V. Pilipenko

Abstract: In this work an ensemble of simulated Local Group analogues is used to constrain the properties of the mass assembly history of the Milky Way (MW) and Andromeda (M31) galaxies. These objects have been obtained using the constrained simulation technique, which ensures that simulated LGs live within a large scale environment akin to the observed one. Our results are compared against a standard $Λ$ C… ▽ More In this work an ensemble of simulated Local Group analogues is used to constrain the properties of the mass assembly history of the Milky Way (MW) and Andromeda (M31) galaxies. These objects have been obtained using the constrained simulation technique, which ensures that simulated LGs live within a large scale environment akin to the observed one. Our results are compared against a standard $Λ$ Cold Dark Matter ($Λ$CDM) series of simulations which use the same cosmological parameters. This allows us to single out the effects of the constraints on the results. We find that (a) the median constrained merging histories for M31 and MW live above the standard ones at the 1-$σ$ level, (b) the median formation time takes place $\approx$ 0.5 Gyr earlier than unconstrained values, while the latest major merger happens on average 1.5 Gyr earlier and (c) the probability for both LG haloes to have experienced their last major merger in the first half of the history of the Universe is $\approx$ 50% higher for the constrained pairs. These results have been estimated to be significant at the 99% confidence level by means of a Kolmogorov-Simirnov test. These results suggest that the particular environment in which the Milky Way and Andromeda form plays a role in shaping their properties, and favours earlier formation and last major merger time values in agreement with other observational and theoretical considerations. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: 10 pages, 3 figures. Accepted for publication in MNRAS

The orientation of planes of dwarf galaxies in the quasi-linear Universe

Authors: Noam I Libeskind, Edoardo Carlesi, Oliver Muller, Marcel S Pawlowski, Yehuda Hoffman, Daniel Pomarede, Helene M Courtois, R. Brent Tully, Stefan Gottloeber, Matthias Steinmetz, Jenny Sorce, Alexander Knebe

Abstract: To date at least 10 highly flattened planes of dwarf galaxies are claimed to have been discovered in the Local Universe. The origin of these planes of galaxies remains unknown. One suggestion is that they are related to the large-scale structure of the cosmic web. A recent study found that the normal of a number of these dwarf galaxy planes are very closely aligned with the eigenvector of the shea… ▽ More To date at least 10 highly flattened planes of dwarf galaxies are claimed to have been discovered in the Local Universe. The origin of these planes of galaxies remains unknown. One suggestion is that they are related to the large-scale structure of the cosmic web. A recent study found that the normal of a number of these dwarf galaxy planes are very closely aligned with the eigenvector of the shear tensor corresponding to the direction of greatest collapse obtained by reconstructing the full velocity field in the linear regime. Here we extend that work by both considering an additional 5 planes beyond the 5 examined previously and by examining the alignment with respect to the quasi-linear field, a more sophisticated reconstruction, which is a better approximation on smaller (quasi-linear) scales. Our analysis recovers the previous result while not finding a significantly tight alignment with the additional 5 planes. However, the additional 5 planes normals also do not appear to be randomly oriented. We conclude that this could either be due to the normals of the new planes being poorly defined and described; the quasi-linear field at those locations being poorly constrained; or different formation mechanisms for the orientation of planes of dwarf galaxies. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 9 pages, 3 figures, 1 table

Cosmicflows-3: Cosmography of the Local Void

Authors: R Brent Tully, Daniel Pomarede, Romain Graziani, Helene M Courtois, Yehuda Hoffman, Edward J Shaya

Abstract: Cosmicflows-3 distances and inferred peculiar velocities of galaxies have permitted the reconstruction of the structure of over and under densities within the volume extending to 0.05c. This study focuses on the under dense regions, particularly the Local Void that lies largely in the zone of obscuration and consequently has received limited attention. Major over dense structures that bound the Lo… ▽ More Cosmicflows-3 distances and inferred peculiar velocities of galaxies have permitted the reconstruction of the structure of over and under densities within the volume extending to 0.05c. This study focuses on the under dense regions, particularly the Local Void that lies largely in the zone of obscuration and consequently has received limited attention. Major over dense structures that bound the Local Void are the Perseus-Pisces and Norma-Pavo-Indus filaments separated by 8,500 km/s. The void network of the universe is interconnected and void passages are found from the Local Void to the adjacent very large Hercules and Sculptor voids. Minor filaments course through voids. A particularly interesting example connects the Virgo and Perseus clusters, with several substantial galaxies found along the chain in the depths of the Local Void. The Local Void has a substantial dynamical effect, causing a deviant motion of the Local Group of 200-250 km/s. The combined perturbations due to repulsion from the Local Void and attraction toward the Virgo Cluster account for about 50% of the motion of the Local Group in the rest frame given by the cosmic microwave background. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: 18 pages; 9 figures; links to 1 video and 2 interactive models

The peculiar velocity field up to $z \sim 0.05$ by forward-modeling Cosmicflows-3 data

Authors: R. Graziani, H. M. Courtois, G. Lavaux, Y. Hoffman, R. B. Tully, Y. Copin, D. Pomarède

Abstract: A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $Λ$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate bia… ▽ More A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $Λ$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate biases inherent to peculiar velocity analyses, such as the Homogeneous Malmquist bias or the log-normal distribution of peculiar velocities. The statistical uncertainty on the recovered velocity field is about 150 km/s depending on the location, and we study systematics coming from the selection function and calibration of distance indicators. The resulting velocity field and related density fields recover the cosmography of the Local Universe which is presented in an unprecedented volume of universe 10 times larger than previously reached. This methodology open the doors to reconstruction of initial conditions for larger and more accurate constrained cosmological simulations. This work is also preparatory to larger peculiar velocity datasets coming from Wallaby, TAIPAN or LSST. △ Less

Submitted 7 January, 2019; originally announced January 2019.

Comments: 15 pages, 14 figures, 1 table. Accepted by MNRAS

Cosmic Dawn II (CoDa II): a new radiation-hydrodynamics simulation of the self-consistent coupling of galaxy formation and reionization

Authors: Pierre Ocvirk, Dominique Aubert, Jenny G. Sorce, Paul R. Shapiro, Nicolas Deparis, Taha Dawoodbhoy, Joseph Lewis, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Kyungjin Ahn, Ilian T. Iliev, Yehuda Hoffman

Abstract: Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RA… ▽ More Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $Λ$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other. △ Less

Submitted 14 October, 2020; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Cosmic-Ray Anisotropy from Large Scale Structure and the effect of magnetic horizons

Authors: Noémie Globus, Tsvi Piran, Yehuda Hoffman, Edoardo Carlesi, Daniel Pomarède

Abstract: Motivated by the ~7 % dipole anisotropy in the distribution of ultra-high energy cosmic-rays (UHECRs) above 8 EeV, we explore the anisotropy induced by the large scale structure, using constrained simulations of the local Universe and taking into account the effect of magnetic fields. The value of the intergalactic magnetic field (IGMF) is critical as it determines the UHECR cosmic horizon. We cal… ▽ More Motivated by the ~7 % dipole anisotropy in the distribution of ultra-high energy cosmic-rays (UHECRs) above 8 EeV, we explore the anisotropy induced by the large scale structure, using constrained simulations of the local Universe and taking into account the effect of magnetic fields. The value of the intergalactic magnetic field (IGMF) is critical as it determines the UHECR cosmic horizon. We calculate the UHECR sky maps for different values of the IGMF variance and show the effect of the UHECR horizon on the observed anisotropy. The footprint of the local ($\lesssim350$ Mpc) Universe on the UHECR background, a small angular scale enhancement in the Northern Hemisphere, is seen. At 11.5 EeV (the median value of the energy bin at which the dipole has been reported), the LSS-induced dipole amplitude is $A_1\sim10\%$, for IGMF in the range [0.3-3] nG for protons, helium and nitrogen, compatible with the rms value derived from the cosmic power spectrum. However at these energies the UHECRs are also influenced by the Galactic Magnetic Field (GMF) and we discuss its effect on the LSS-induced anisotropy. △ Less

Submitted 19 November, 2018; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 7 pages, 6 figures, submitted to MNRAS (revised version)

The quasi-linear nearby Universe

Authors: Yehuda Hoffman, Edoardo Carlesi, Daniel Pomarède, Hélène M. Courtois, Stefan Gottlöber, Noam I. Libeskind, Jenny G. Sorce, Gustavo Yepes

Abstract: The local Universe provides a unique opportunity for testing cosmology and theories of structure formation. To facilitate this opportunity we present a new method for the reconstruction of the quasi-linear matter density and velocity fields from galaxy peculiar velocities and apply it to the Cosmicflows-2 data. The method consists of constructing an ensemble of cosmological simulations, constraine… ▽ More The local Universe provides a unique opportunity for testing cosmology and theories of structure formation. To facilitate this opportunity we present a new method for the reconstruction of the quasi-linear matter density and velocity fields from galaxy peculiar velocities and apply it to the Cosmicflows-2 data. The method consists of constructing an ensemble of cosmological simulations, constrained by the standard cosmological model and the observational data. The quasi-linear density field is the geometric mean and variance of the fully non-linear density fields of the simulations. The main nearby clusters (Virgo, Centaurus, Coma), superclusters (Shapley, Perseus-Pisces) and voids (Dipole Repeller) are robustly reconstructed. Galaxies are born biased with respect to the underlying dark matter distribution. Using our quasi-linear framework we demonstrate that the luminosity-weighted density field derived from the 2M++ redshift compilations is non-linearly biased with respect to the matter density field. The bias diminishes in the linear regime. △ Less

Submitted 10 July, 2018; originally announced July 2018.

Comments: Accepted for publication in Nature Astronomy. Full article available here https://www.nature.com/articles/s41550-018-0502-4. A version of the arxiv paper including higher resolution versions of the paper images can be obtained from https://edoardocarlesi.github.io/files/CF2_QL_NatAst_arxiv.pdf

Journal ref: Nature Astronomy 2018

Suppression of Star Formation in Low-Mass Galaxies Caused by the Reionization of their Local Neighborhood

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Dominique Aubert, Nicolas Gillet, Jun-Hwan Choi, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman

Abstract: Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical s… ▽ More Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later. △ Less

Submitted 25 July, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 15 pages, 15 figures, accepted by MNRAS 07/17/18, comments welcome

Journal ref: MNRAS, 480, 1740-1753 (2018)

The Inhomogeneous Reionization Times of Present-day Galaxies

Authors: Dominique Aubert, Nicolas Deparis, Pierre Ocvirk, Paul R. Shapiro, Ilian T. Iliev, Gustavo Yepes, Stefan Gottloeber, Yehuda Hoffman, Romain Teyssier

Abstract: Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen t… ▽ More Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both). △ Less

Submitted 20 February, 2018; v1 submitted 5 February, 2018; originally announced February 2018.

Comments: 8 pages, 5 figures, updated version accepted for publication in ApJL

Action Dynamics of the Local Supercluster

Authors: Edward J. Shaya, R. Brent Tully, Yehuda Hoffman, Daniel Pomarède

Abstract: The fully nonlinear gravitationally induced trajectories of a nearly complete set of galaxies, groups, and clusters in the Local Supercluster are constructed in a Numerical Action Method (NAM) model constrained by data from the CosmicFlows survey and various distance indicators. We add the gravity field due to inhomogeneities external to the sample sphere by making use of larger scale peculiar flo… ▽ More The fully nonlinear gravitationally induced trajectories of a nearly complete set of galaxies, groups, and clusters in the Local Supercluster are constructed in a Numerical Action Method (NAM) model constrained by data from the CosmicFlows survey and various distance indicators. We add the gravity field due to inhomogeneities external to the sample sphere by making use of larger scale peculiar flow measurements. Assignments of total masses were made to find the best overall set of mutual attractions, as determined by a goodness criterion based on present day radial velocities, individually for the Virgo Cluster, M31, and the Milky Way (MW), and via a mass-to-light ratio relationship for other masses. The low median chi-square found indicates the model fits the present day velocity flow well, but a slightly high mean chi-square may indicate that some masses underwent complex orbits. The best fit, when setting the value of $H_0$ to the CosmicFlows value of 75 km/s/Mpc and the WMAP value for $Ω_m=0.244$ consistent with that $H_0$, occurs with the following parameters: $Ω_{orphan} = 0.077\pm0.016$, $M/L_K $= $40 \pm 2$ $L_{10}^{0.15}$ $M_\odot/L_\odot$ ($L_{10}$ is K-band luminosity in units of $10^{10} L_\odot$), a Virgo mass of $6.3\pm0.8 \times 10^{14} M_\odot$ ($M/L_K = 113\pm15 M_\odot/L_\odot$), and a mass for the MW plus M31 of $5.15\pm0.35 \times 10^{12} M_\odot$. The best constant mass-to-light ratio is $M/L_K = 58\pm3 M_\odot/L_\odot$. The Virgocentric turnaround radius is $7.3\pm0.3$ Mpc. We explain several interesting trends in peculiar motions for various regions now that we can construct the 3D orbital histories. △ Less

Submitted 26 October, 2017; v1 submitted 24 October, 2017; originally announced October 2017.

Comments: 18 pages, 17 Figures, links to online video, interactive plots and a large table

Universal subhalo accretion in cold and warm dark matter cosmologies

Authors: Bogna Kubik, Noam I. Libeskind, Alexander Knebe, Hélène Courtois, Gustavo Yepes, Stefan Gottlöber, Yehuda Hoffman

Abstract: The influence of the large scale structure on host halos may be studied by examining the angular infall pattern of subhalos. In particular, since warm and cold dark matter cosmologies predict different abundances and internal properties for halos at the low mass end of the mass function, it is interesting to examine if there are differences in how these low mass halos are accreted. The accretion e… ▽ More The influence of the large scale structure on host halos may be studied by examining the angular infall pattern of subhalos. In particular, since warm and cold dark matter cosmologies predict different abundances and internal properties for halos at the low mass end of the mass function, it is interesting to examine if there are differences in how these low mass halos are accreted. The accretion events are defined as the moment a halo becomes a substructure, namely when it crosses its host's virial radius. We quantify the cosmic web at each point by the shear tensor and examine where, with respect to its eigenvectors, such accretion events occur in cold ($Λ$CDM) and warm (1keV sterile neutrino WDM) dark matter cosmological models. We find that the CDM and WDM subhalos are preferentially accreted along the principal axis of the shear tensor corresponding to the direction of weakest collapse. The beaming strength is modulated by the host and subhalo masses and by the redshift at which the accretion event occurs. Although strongest for the most massive hosts and subhalos at high redshift, the preferential infall is found to be always aligned with the axis of weakest collapse, thus we say that it has universal nature. We compare the strength of beaming in the WDM cosmology with the one found in the $Λ$CDM scenario. While the main findings remain the same, the accretion in the WDM model for the most massive host halos appears more beamed than in $Λ$CDM cosmology across all the redshifts. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: 12 pages, 9 figures

Cosmicflows-3: Cold Spot Repeller?

Authors: Helene M. Courtois, R. Brent Tully, Yehuda Hoffman, Daniel Pomarede, Romain Graziani, Alexandra Dupuy

Abstract: The three-dimensional gravitational velocity field within z~0.1 has been modeled with the Wiener filter methodology applied to the Cosmicflows-3 compilation of galaxy distances. The dominant features are a basin of attraction and two basins of repulsion. The major basin of attraction is an extension of the Shapley concentration of galaxies. One basin of repulsion, the Dipole Repeller, is located n… ▽ More The three-dimensional gravitational velocity field within z~0.1 has been modeled with the Wiener filter methodology applied to the Cosmicflows-3 compilation of galaxy distances. The dominant features are a basin of attraction and two basins of repulsion. The major basin of attraction is an extension of the Shapley concentration of galaxies. One basin of repulsion, the Dipole Repeller, is located near the anti-apex of the cosmic microwave background dipole. The other basin of repulsion is in the proximate direction toward the 'Cold Spot' irregularity in the cosmic microwave background. It has been speculated that a vast void might contribute to the amplitude of the Cold Spot from the integrated Sachs-Wolfe effect. △ Less

Submitted 24 August, 2017; originally announced August 2017.

Comments: Accepted for publication in The Astrophysical Journal Letters

The Cosmic V-Web

Authors: Daniel Pomarede, Yehuda Hoffman, Helene M Courtois, R. Brent Tully

Abstract: The network of filaments with embedded clusters surrounding voids seen in maps derived from redshift surveys and reproduced in simulations has been referred to as the cosmic web. A complementary description is provided by considering the shear in the velocity field of galaxies. The eigenvalues of the shear provide information on whether a region is collapsing in three dimensions, the condition for… ▽ More The network of filaments with embedded clusters surrounding voids seen in maps derived from redshift surveys and reproduced in simulations has been referred to as the cosmic web. A complementary description is provided by considering the shear in the velocity field of galaxies. The eigenvalues of the shear provide information on whether a region is collapsing in three dimensions, the condition for a knot, expanding in three-dimensions, the condition for a void, or in the intermediate condition of a filament or sheet. The structures that are quantitatively defined by the eigenvalues can be approximated by iso-contours that provide a visual representation of the cosmic velocity (V) web. The current application is based on radial peculiar velocities from the Cosmicflows-2 collection of distances. The three-dimensional velocity field is constructed using the Wiener filter methodology in the linear approximation. Eigenvalues of the velocity shear are calculated at each point on a grid. Here, knots and filaments are visualized across a local domain of diameter ~0.1c. △ Less

Submitted 11 June, 2017; originally announced June 2017.

Comments: Accepted for publication in The Astrophysical Journal. Accompanied by a video and an interactive figure

Tracing the cosmic web

Authors: Noam I Libeskind, Rien van de Weygaert, Marius Cautun, Bridget Falck, Elmo Tempel, Tom Abel, Mehmet Alpaslan, Miguel A. Aragoon-Calvo, Jaime E. Forero-Romero, Roberto Gonzalez, Stefan Gottloober, Oliver Hahn, Wojciech A. Hellwing, Yehuda Hoffman, Bernard J. T. Jones, Francisco Kitaura, Alexander Knebe, Serena Manti, Mark Neyrinck, Sebastiaan E. Nuza, Nelson Padilla, Erwin Platen, Nesar Ramachandra, Aaron Robotham, Enn Saar , et al. (5 additional authors not shown)

Abstract: The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effe… ▽ More The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public. △ Less

Submitted 8 May, 2017; originally announced May 2017.

Comments: 24 pages, 8 figures, 2 tables. Submitted to MN. Comments Welcome

Journal ref: MNRAS, 2018, volume 473, pages 1195-1217

Reionization of the Milky Way, M31, and their satellites I: reionization history and star formation

Authors: Keri L. Dixon, Ilian T. Iliev, Stefan Gottlöber, Gustavo Yepes, Alexander Knebe, Noam Libeskind, Yehuda Hoffman

Abstract: Observations of the Milky Way (MW), M31, and their vicinity, known as the Local Group (LG), can provide clues about the sources of reionization. We present a suite of radiative transfer simulations based on initial conditions provided by the Constrained Local UniversE Simulations (CLUES) project that are designed to recreate the Local Universe, including a realistic MW-M31 pair and a nearby Virgo.… ▽ More Observations of the Milky Way (MW), M31, and their vicinity, known as the Local Group (LG), can provide clues about the sources of reionization. We present a suite of radiative transfer simulations based on initial conditions provided by the Constrained Local UniversE Simulations (CLUES) project that are designed to recreate the Local Universe, including a realistic MW-M31 pair and a nearby Virgo. Our box size (91 Mpc) is large enough to incorporate the relevant sources of ionizing photons for the LG. We employ a range of source models, mimicking the potential effects of radiative feedback for dark matter haloes between $10^{8}-10^{9}$ M$_{\odot}$. Although the LG mostly reionizes in an inside-out fashion, the final 40 per cent of its ionization shows some outside influence. For the LG satellites, we find no evidence that their redshift of reionization is related to the present-day mass of the satellite or the distance from the central galaxy. We find that less than 20 per cent of present-day satellites for MW and M31 have undergone any star formation prior to the end of global reionization. Approximately five per cent of these satellites could be classified as fossils, meaning the majority of star formation occurred at these early times. The more massive satellites have more cumulative star formation prior to the end of global reionization, but the scatter is significant, especially at the low-mass end. Present-day mass and distance from the central galaxy are poor predictors for the presence of ancient stellar populations in satellite galaxies. △ Less

Submitted 17 March, 2017; originally announced March 2017.

Comments: 14 pages, 12 figures, 2 tables. Submitted to MNRAS

Toward an Optimal Sampling of Peculiar Velocity Surveys For Wiener Filter Reconstructions

Authors: Jenny G. Sorce, Yehuda Hoffman, Stefan Gottlöber

Abstract: The Wiener Filter (WF) technique enables the reconstruction of density and velocity fields from observed radial peculiar velocities. This paper aims at identifying the optimal design of peculiar velocity surveys within the WF framework. The prime goal is to test the dependence of the quality of the reconstruction on the distribution and nature of data points. Mock datasets, extending to 250 Mpc/h,… ▽ More The Wiener Filter (WF) technique enables the reconstruction of density and velocity fields from observed radial peculiar velocities. This paper aims at identifying the optimal design of peculiar velocity surveys within the WF framework. The prime goal is to test the dependence of the quality of the reconstruction on the distribution and nature of data points. Mock datasets, extending to 250 Mpc/h, are drawn from a constrained simulation that mimics the local Universe to produce realistic mock catalogs. Reconstructed fields obtained with these mocks are compared to the reference simulation. Comparisons, including residual distributions, cell-to-cell and bulk velocities, imply that the presence of field data points is essential to properly measure the flows. The fields reconstructed from mocks that consist only of galaxy cluster data points exhibit poor quality bulk velocities. In addition, the quality of the reconstruction depends strongly on the grouping of individual data points into single points to suppress virial motions in high density regions. Conversely, the presence of a Zone of Avoidance hardly affects the reconstruction. For a given number of data points, a uniform sample does not score any better than a sample with decreasing number of data points with the distance. The best reconstructions are obtained with a grouped survey containing field galaxies: Assuming no error, they differ from the simulated field by less than 100 km/s up to the extreme edge of the catalogs or up to a distance of three times the mean distance of data points for non-uniform catalogs. The overall conclusions hold when errors are added. △ Less

Submitted 6 March, 2017; originally announced March 2017.

Comments: Accepted for publication in MNRAS, 12 pages, 7 figures, 1 table

The Dipole Repeller

Authors: Yehuda Hoffman, Daniel Pomarede, R. Brent Tully, Helene Courtois

Abstract: In the standard (LCDM) model of cosmology the universe has emerged out of an early homogeneous and isotropic phase. Structure formation is associated with the growth of density irregularities and peculiar velocities. Our Local Group is moving with respect to the cosmic microwave background (CMB) with a velocity 631+/-20 km s-1 and participates in a bulk flow that extends out to distances of at lea… ▽ More In the standard (LCDM) model of cosmology the universe has emerged out of an early homogeneous and isotropic phase. Structure formation is associated with the growth of density irregularities and peculiar velocities. Our Local Group is moving with respect to the cosmic microwave background (CMB) with a velocity 631+/-20 km s-1 and participates in a bulk flow that extends out to distances of at least 20,000 km s-1. Since the discovery of the CMB dipole, the implicit assumption was that excesses in the abundance of galaxies induce the Local Group motion. Yet, underdense regions push as much as overdensities attract but they are deficient of light and consequently difficult to chart. It was suggested a decade ago that an underdensity in the northern hemisphere roughly 15,000 km s-1 away is a significant actor in the local flow. Here we report on kinematic evidence for such an underdensity. We map the large scale 3D velocity field using a Wiener filter reconstruction from the Cosmicflows-2 dataset of peculiar velocities, and identify the attractors and repellers that dominate the local dynamics. We show here that the local flow is dominated by a single attractor -associated with the Shapley Concentration- and a single previously unidentified repeller. Multipole expansion of the local flow provides further support for the existence and role played by the attractor and repeller. The bulk flow (i.e. dipole moment) is closely (anti)aligned with the repeller at a distance of 16,000+/-4,500 km s-1. The expansion eigenvector of the shear tensor (quadrupole moment) is closely aligned with the Shapley Attractor out to 7,000 km s-1. The close alignment of the local bulk flow with the repeller provides further support for its dominant role in shaping the local flow. This Dipole Repeller is predicted to be associated with a void in the distribution of galaxies. △ Less

Submitted 8 February, 2017; originally announced February 2017.

Comments: published in Nature Astronomy Volume 1, 36

Journal ref: 2017NatAs...1E..36H

Cosmography and Data Visualization

Authors: Daniel Pomarede, Helene M. Courtois, Yehuda Hoffman, R. Brent Tully

Abstract: Cosmography, the study and making of maps of the universe or cosmos, is a field where visual representation benefits from modern three-dimensional visualization techniques and media. At the extragalactic distance scales, visualization is contributing in understanding the complex structure of the local universe, in terms of spatial distribution and flows of galaxies and dark matter. In this paper,… ▽ More Cosmography, the study and making of maps of the universe or cosmos, is a field where visual representation benefits from modern three-dimensional visualization techniques and media. At the extragalactic distance scales, visualization is contributing in understanding the complex structure of the local universe, in terms of spatial distribution and flows of galaxies and dark matter. In this paper, we report advances in the field of extragalactic cosmography obtained using the SDvision visualization software in the context of the Cosmicflows Project. Here, multiple visualization techniques are applied to a variety of data products: catalogs of galaxy positions and galaxy peculiar velocities, reconstructed velocity field, density field, gravitational potential field, velocity shear tensor viewed in terms of its eigenvalues and eigenvectors, envelope surfaces enclosing basins of attraction. These visualizations, implemented as high-resolution images, videos, and interactive viewers, have contributed to a number of studies: the cosmography of the local part of the universe, the nature of the Great Attractor, the discovery of the boundaries of our home supercluster of galaxies Laniakea, the mapping of the cosmic web, the study of attractors and repellers. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: Accepted for publication in the PASP Special Focus Issue: Techniques and Methods for Astrophysical Data Visualization

Constraining the mass of the Local Group

Authors: Edoardo Carlesi, Yehuda Hoffman, Jenny G. Sorce, Stefan Gottlöber

Abstract: The mass of the Local Group (LG) is a crucial parameter for galaxy formation theories. However, its observational determination is challenging - its mass budget is dominated by dark matter which cannot be directly observed. To meet this end the posterior distributions of the LG and its massive constituents have been constructed by means of constrained and random cosmological simulations. Two prior… ▽ More The mass of the Local Group (LG) is a crucial parameter for galaxy formation theories. However, its observational determination is challenging - its mass budget is dominated by dark matter which cannot be directly observed. To meet this end the posterior distributions of the LG and its massive constituents have been constructed by means of constrained and random cosmological simulations. Two priors are assumed - the LCDM model that is used to set up the simulations and a LG model,which encodes the observational knowledge of the LG and is used to select LG-like objects from the simulations. The constrained simulations are designed to reproduce the local cosmography as it is imprinted onto the Cosmicflows-2 database of velocities. Several prescriptions are used to define the LG model, focusing in particular on different recent estimates of the tangential velocity of M31. It is found that (a) different $v_{tan}$ choices affect the peak mass values up to a factor of 2, and change mass ratios of $M_{M31}$ to $M_{MW}$ by up to 20%; (b) constrained simulations yield more sharply peaked posterior distributions compared with the random ones; (c) LG mass estimates are found to be smaller than those found using the timing argument; (d) preferred MW masses lie in the range of $(0.6 - 0.8)\times10^{12} M_{\odot}$ whereas (e) $M_{M31}$ is found to vary between $(1.0 - 2.0)\times10^{12} M_{\odot}$, with a strong dependence on the $v_{tan}$ values used. △ Less

Submitted 24 November, 2016; originally announced November 2016.

Comments: 10 pages, 1 figure, 6 tables. Accepted for publication in MNRAS

Estimating the Mass of the Local Group using Machine Learning Applied to Numerical Simulations

Authors: Michael McLeod, Noam Libeskind, Ofer Lahav, Yehuda Hoffman

Abstract: We revisit the estimation of the combined mass of the Milky Way and Andromeda (M31), which dominate the mass of the Local Group. We make use of an ensemble of 30,190 halo pairs from the Small MultiDark simulation, assuming a $Λ$CDM (Cosmological Constant with Cold Dark Matter) cosmology, to investigate the relationship between the bound mass and parameters characterising the orbit of the binary an… ▽ More We revisit the estimation of the combined mass of the Milky Way and Andromeda (M31), which dominate the mass of the Local Group. We make use of an ensemble of 30,190 halo pairs from the Small MultiDark simulation, assuming a $Λ$CDM (Cosmological Constant with Cold Dark Matter) cosmology, to investigate the relationship between the bound mass and parameters characterising the orbit of the binary and their local environment with the aid of machine learning methods (artificial neural networks, ANN). Results from the ANN are most successful when information about the velocity shear is provided, which demonstrates the flexibility of machine learning to model physical phenomena and readily incorporate new information as it becomes available. The resulting estimate for the Local Group mass, when shear information is included, is $4.9 \times 10^{12} M_\odot$, with an error of $\pm0.8 \times 10^{12} M_\odot$ from the 68% uncertainty in observables, and a 68% confidence interval of $^{+1.3}_{-1.4} \times 10^{12}M_\odot$ from the intrinsic scatter from the differences between the model and simulation masses. We also consider a recently reported large transverse velocity of M31 relative to the Milky Way, and produce an alternative mass estimate of $3.6\pm0.3\pm1.4 \times 10^{12}M_\odot$. Although different methods predict similar values for the most likely mass of the LG, application of ANN compared to the Timing Argument reduces the scatter in the log mass by over half when tested on samples from the simulation. △ Less

Submitted 25 October, 2017; v1 submitted 8 June, 2016; originally announced June 2016.

Comments: 20 pages, 6 figures, 5 tables

The abundance and environment of dark matter haloes

Authors: Ofer Metuki, Noam I Libeskind, Yehuda Hoffman

Abstract: An open question in cosmology and the theory of structure formation is to what extent does environment affect the properties of galaxies and haloes. The present paper aims at shedding light on this problem. The paper focuses on the analysis of a dark matter only simulation and it addresses the issue of how the environment affects the abundance of haloes, which are are assigned four attributes: the… ▽ More An open question in cosmology and the theory of structure formation is to what extent does environment affect the properties of galaxies and haloes. The present paper aims at shedding light on this problem. The paper focuses on the analysis of a dark matter only simulation and it addresses the issue of how the environment affects the abundance of haloes, which are are assigned four attributes: their virial mass, an ambient density calculated with an aperture that scales with $R_{vir}$ ($Δ_M$), a fixed-aperture ($Δ_R$) ambient density, and a cosmic web classification (i.e. voids, sheets, filaments, and knots, as defined by the V--web algorithm). $Δ_M$ is the mean density around a halo evaluated within a sphere of a radius of $5$\rvir, where \rvir\ is the virial radius. $Δ_R$ is the density field Gaussian smoothed with $R=4h^{-1}$Mpc, evaluated at the center of the halo. The main result of the paper is that the difference between haloes in different web elements stems from the difference in their mass functions, and does not depend on their adaptive-aperture ambient density. A dependence on the fixed-aperture ambient density is induced by the cross correlation between the mass of a halo and its fixed-aperture ambient density. △ Less

Submitted 5 June, 2016; originally announced June 2016.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS

How did the Virgo cluster form?

Authors: Jenny G. Sorce, Stefan Gottloeber, Yehuda Hoffman, Gustavo Yepes

Abstract: While the Virgo cluster is the nearest galaxy cluster and therefore the best observed one, little is known about its formation history. In this paper, a set of cosmological simulations that resemble the Local Universe is used to shed the first light on this mystery. The initial conditions for these simulations are constrained with galaxy peculiar velocities of the second catalog of the Cosmicflows… ▽ More While the Virgo cluster is the nearest galaxy cluster and therefore the best observed one, little is known about its formation history. In this paper, a set of cosmological simulations that resemble the Local Universe is used to shed the first light on this mystery. The initial conditions for these simulations are constrained with galaxy peculiar velocities of the second catalog of the Cosmicflows project using algorithms developed within the Constrained Local UniversE Simulation project. Boxes of 500 Mpc/h on a side are set to run a series of dark matter only constrained simulations. In each simulation, a unique dark matter halo can be reliably identified as Virgo's counterpart. The properties of these Virgo halos are in agreement at a 10-20% level with the global properties of the observed Virgo cluster. Their zero-velocity masses agree at one-sigma with the observational mass estimate. In all the simulations, the matter falls onto the Virgo objects along a preferential direction that corresponds to the observational filament and the slowest direction of collapse. A study of the mass accretion history of the Virgo candidates reveals the most likely formation history of the Virgo cluster, namely a quiet accretion over the last 7 Gigayears. △ Less

Submitted 22 May, 2016; originally announced May 2016.

Comments: Accepted for publication in MNRAS, 10 pages, 10 figures

Goodness-of-fit analysis of the Cosmicflows-2 database of velocities

Authors: Yehuda Hoffman, Adi Nusser, Helene M. Courtois, R. Brent Tully

Abstract: The goodness-of-fit (GoF) of the Cosmicflows-2 (CF2) database of peculiar velocities with the LCDM standard model of cosmology is presented. Standard application of the Chi^2 statistics of the full database, of its 4,838 data points, is hampered by the small scale nonlinear dynamics which is not accounted for by the (linear regime) velocity power spectrum. The bulk velocity constitutes a highly co… ▽ More The goodness-of-fit (GoF) of the Cosmicflows-2 (CF2) database of peculiar velocities with the LCDM standard model of cosmology is presented. Standard application of the Chi^2 statistics of the full database, of its 4,838 data points, is hampered by the small scale nonlinear dynamics which is not accounted for by the (linear regime) velocity power spectrum. The bulk velocity constitutes a highly compressed representation of the data which filters out the small scales non-linear modes. Hence the statistics of the bulk flow provides an efficient tool for assessing the GoF of the data given a model. The particular approach introduced here is to use the (spherical top-hat window) bulk velocity extracted from the Wiener filter reconstruction of the 3D velocity field as a linear low pass filtered highly compressed representation of the CF2 data. An ensemble 2250 random linear realizations of the WMAP/LCDM model has been used to calculate the bulk velocity auto-covariance matrix. We find that the CF2 data is consistent with the WMAP/LCDM model to better than the 2 sigma confidence limits. This provides a further validation that the CF2 database is consistent with the standard model of cosmology. △ Less

Submitted 11 May, 2016; v1 submitted 8 May, 2016; originally announced May 2016.

Comments: submitted to MNRAS, V2 : solved page sizing problem

The tangential velocity of M31: CLUES from constrained simulations

Authors: Edoardo Carlesi, Yehuda Hoffman, Jenny G. Sorce, Stefan Gottloeber, Gustavo Yepes, Helene Courtois, R. Brent Tully

Abstract: Determining the precise value of the tangential component of the velocity of M31 is a non trivial astrophysical issue, that relies on complicated modeling. This has recently lead to con- flicting estimates, obtained by several groups that used different methodologies and assump- tions. This letter addresses the issue by computing a Bayesian posterior distribution function of this quantity, in orde… ▽ More Determining the precise value of the tangential component of the velocity of M31 is a non trivial astrophysical issue, that relies on complicated modeling. This has recently lead to con- flicting estimates, obtained by several groups that used different methodologies and assump- tions. This letter addresses the issue by computing a Bayesian posterior distribution function of this quantity, in order to measure the compatibility of those estimates with LambdaCDM. This is achieved using an ensemble of local group (LG) look-alikes collected from a set of Con- strained Simulations (CSs) of the local Universe, and a standard unconstrained LambdaCDM. The latter allows us to build a control sample of LG-like pairs and to single out the influence of the environment in our results. We find that neither estimate is at odds with LambdaCDM; how- ever, whereas CSs favour higher values of vtan , the reverse is true for estimates based on LG samples gathered from unconstrained simulations, overlooking the environmental element △ Less

Submitted 31 March, 2016; originally announced March 2016.

Comments: 5 pages, 2 figures, accepted for publications in MNRAS Letters

Constrained Local UniversE Simulations: A Local Group Factory

Authors: Edoardo Carlesi, Jenny G. Sorce, Yehuda Hoffman, Stefan Gottlöber, Gustavo Yepes, Noam I. Libeskind, Sergey V. Pilipenko, Alexander Knebe, Helene Courtois, R. Brent Tully, Matthias Steinmetz

Abstract: Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations… ▽ More Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations of our cosmic environment are performed. The aim of these simulations is to reproduce the LG and its local environment. Our main result is that the LG is likely a robust outcome of the LCDM scenario when subjected to the constraint derived from CF2 data, emerging in an environment akin to the observed one. Three levels of criteria are used to define the simulated LGs. At the base level, pairs of halos must obey specific isolation, mass and separation criteria. At the second level the orbital angular momentum and energy are constrained and on the third one the phase of the orbit is constrained. Out of the 300 constrained simulations 146 LGs obey the first set of criteria, 51 the second and 6 the third. The robustness of our LG factory enables the construction of a large ensemble of simulated LGs. Suitable candidates for high resolution hydrodynamical simulations of the LG can be drawn from this ensemble, which can be used to perform comprehensive studies of the formation of the LG △ Less

Submitted 11 February, 2016; originally announced February 2016.

Comments: 13 pages, accepted for publication in MNRAS

The alignment of galaxy spin with the shear field in observations

Authors: Isha Pahwa, Noam I. Libeskind, Elmo Tempel, Yehuda Hoffman, R. Brent Tully, Helene M. Courtois, Stefan Gottlöber, Matthias Steinmetz, Jenny G. Sorce

Abstract: Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities comb… ▽ More Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction. △ Less

Submitted 14 December, 2015; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: 9 pages, 3 figures, accepted in MNRAS

Journal ref: MNRAS 457 (2016) 695

Mergers and the outside-in formation of dwarf spheroidals

Authors: Alejandro Benítez-Llambay, Julio F. Navarro, Mario G. Abadi, Stefan Gottloeber, Gustavo Yepes, Yehuda Hoffman, Matthias Steinmetz

Abstract: We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few example… ▽ More We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals. △ Less

Submitted 19 November, 2015; originally announced November 2015.

Comments: 10 pages, 8 figures, Accepted for publication in MNRAS

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Authors: Pierre Ocvirk, Nicolas Gillet, Paul R. Shapiro, Dominique Aubert, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Jun-Hwan Choi, David Sullivan, Alexander Knebe, Stefan Gottloeber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman, Timothy Stranex

Abstract: Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while reso… ▽ More Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged) △ Less

Submitted 6 September, 2016; v1 submitted 30 October, 2015; originally announced November 2015.

Comments: 26 pages, 16 figures, accepted in MNRAS