Showing 1–50 of 71 results for author: Sorce, J G

The Glow of Axion Quark Nugget Dark Matter: (II) Galaxy Clusters

Authors: Julian S. Sommer, Klaus Dolag, Ludwig M. Böss, Ildar Khabibullin, Xunyu Liang, Ludovic Van Waerbeke, Ariel Zhitnitsky, Fereshteh Majidi, Jenny G. Sorce, Benjamin Seidel, Elena Hernández-Martínez

Abstract: (abridged) We analyze the emission of axion quark nuggets in a large sample of 161 simulated galaxy clusters using the SLOW simulation. These clusters are divided into a sub-sample of 150 galaxy clusters, ordered in five mass bins ranging from $0.8$ to $31.7 \times 10^{14} \,M_\odot$, along with 11 cross-identified galaxy clusters from observations. We investigate dark matter-baryonic matter inter… ▽ More (abridged) We analyze the emission of axion quark nuggets in a large sample of 161 simulated galaxy clusters using the SLOW simulation. These clusters are divided into a sub-sample of 150 galaxy clusters, ordered in five mass bins ranging from $0.8$ to $31.7 \times 10^{14} \,M_\odot$, along with 11 cross-identified galaxy clusters from observations. We investigate dark matter-baryonic matter interactions in galaxy clusters in their present stage at redshift $z=0$ by assuming all dark matter consists of axion quark nuggets. The resulting electromagnetic signatures are compared to thermal Bremsstrahlung and non-thermal cosmic ray synchrotron emission in each galaxy cluster. We further investigate individual frequency bands imitating the observable range of the WMAP, Planck, Euclid, and XRISM telescopes for the most promising cross-identified galaxy clusters hosting detectable signatures of axion quark nugget emission. We propose that the Fornax and Virgo clusters represent the most promising candidates to search for axion quark nugget emission signatures. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 22 pages, 14 figures

Euclid: Early Release Observations -- The intracluster light and intracluster globular clusters of the Perseus cluster

Authors: M. Kluge, N. A. Hatch, M. Montes, J. B. Golden-Marx, A. H. Gonzalez, J. -C. Cuillandre, M. Bolzonella, A. Lançon, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, A. Boselli, M. Cantiello, J. G. Sorce, F. R. Marleau, P. -A. Duc, E. Sola, M. Urbano, S. L. Ahad, Y. M. Bahé, S. P. Bamford, C. Bellhouse, F. Buitrago, P. Dimauro , et al. (163 additional authors not shown)

Abstract: We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus clu… ▽ More We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus cluster hosts 70000$\pm$2800 GCs and $1.6\times10^{12}$ L$_\odot$ of diffuse light from the BCG+ICL in the near-infrared H$_E$. This accounts for 37$\pm$6% of the cluster's total stellar luminosity within this radius. The ICL and ICGCs share a coherent spatial distribution, suggesting a common origin or that a common potential governs their distribution. Their contours on the largest scales (>200 kpc) are offset from the BCG's core westwards by 60 kpc towards several luminous cluster galaxies. This offset is opposite to the displacement observed in the gaseous intracluster medium. The radial surface brightness profile of the BCG+ICL is best described by a double Sérsic model, with 68$\pm$4% of the H$_E$ light in the extended, outer component. The transition between these components occurs at ~50 kpc, beyond which the isophotes become increasingly elliptical and off-centred. The radial ICGC number density profile closely follows the BCG+ICL profile only beyond this 50 kpc radius, where we find an average of 60 GCs per $10^9$ M$_\odot$ of diffuse stellar mass. The BCG+ICL colour becomes increasingly blue with radius, consistent with the stellar populations in the ICL having subsolar metallicities [Fe/H]~-0.6. The colour of the ICL, and the specific frequency and luminosity function of the ICGCs suggest that the ICL+ICGCs were tidally stripped from the outskirts of massive satellites with masses of a few $\times10^{10}$ M$_\odot$, with an increasing contribution from dwarf galaxies at large radii. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations. 24 pages, 17 figures, 5 tables. Submitted to A&A

Euclid: Early Release Observations -- Dwarf galaxies in the Perseus galaxy cluster

Authors: F. R. Marleau, J. -C. Cuillandre, M. Cantiello, D. Carollo, P. -A. Duc, R. Habas, L. K. Hunt, P. Jablonka, M. Mirabile, M. Mondelin, M. Poulain, T. Saifollahi, R. Sánchez-Janssen, E. Sola, M. Urbano, R. Zöller, M. Bolzonella, A. Lançon, R. Laureijs, O. Marchal, M. Schirmer, C. Stone, A. Boselli, A. Ferré-Mateu, N. A. Hatch , et al. (171 additional authors not shown)

Abstract: We make use of the unprecedented depth, spatial resolution, and field of view of the Euclid Early Release Observations of the Perseus galaxy cluster to detect and characterise the dwarf galaxy population in this massive system. The Euclid high resolution VIS and combined VIS+NIR colour images were visually inspected and dwarf galaxy candidates were identified. Their morphologies, the presence of n… ▽ More We make use of the unprecedented depth, spatial resolution, and field of view of the Euclid Early Release Observations of the Perseus galaxy cluster to detect and characterise the dwarf galaxy population in this massive system. The Euclid high resolution VIS and combined VIS+NIR colour images were visually inspected and dwarf galaxy candidates were identified. Their morphologies, the presence of nuclei, and their globular cluster (GC) richness were visually assessed, complementing an automatic detection of the GC candidates. Structural and photometric parameters, including Euclid filter colours, were extracted from 2-dimensional fitting. Based on this analysis, a total of 1100 dwarf candidates were found across the image, with 638 appearing to be new identifications. The majority (96%) are classified as dwarf ellipticals, 53% are nucleated, 26% are GC-rich, and 6% show disturbed morphologies. A relatively high fraction of galaxies, 8%, are categorised as ultra-diffuse galaxies. The majority of the dwarfs follow the expected scaling relations. Globally, the GC specific frequency, S_N, of the Perseus dwarfs is intermediate between those measured in the Virgo and Coma clusters. While the dwarfs with the largest GC counts are found throughout the Euclid field of view, those located around the east-west strip, where most of the brightest cluster members are found, exhibit larger S_N values, on average. The spatial distribution of the dwarfs, GCs, and intracluster light show a main iso-density/isophotal centre displaced to the west of the bright galaxy light distribution. The ERO imaging of the Perseus cluster demonstrates the unique capability of Euclid to concurrently detect and characterise large samples of dwarfs, their nuclei, and their GC systems, allowing us to construct a detailed picture of the formation and evolution of galaxies over a wide range of mass scales and environments. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 44 pages, 24 figures, 5 tables, paper submitted to A&A as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Overview of the Perseus cluster and analysis of its luminosity and stellar mass functions

Authors: J. -C. Cuillandre, M. Bolzonella, A. Boselli, F. R. Marleau, M. Mondelin, J. G. Sorce, C. Stone, F. Buitrago, Michele Cantiello, K. George, N. A. Hatch, L. Quilley, F. Mannucci, T. Saifollahi, R. Sánchez-Janssen, F. Tarsitano, C. Tortora, X. Xu, H. Bouy, S. Gwyn, M. Kluge, A. Lançon, R. Laureijs, M. Schirmer, Abdurro'uf , et al. (177 additional authors not shown)

Abstract: The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exception… ▽ More The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 35 figures, Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Programme overview and pipeline for compact- and diffuse-emission photometry

Authors: J. -C. Cuillandre, E. Bertin, M. Bolzonella, H. Bouy, S. Gwyn, S. Isani, M. Kluge, O. Lai, A. Lançon, D. A. Lang, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, Abdurro'uf, N. Aghanim, B. Altieri, F. Annibali, H. Atek, P. Awad, M. Baes, E. Bañados, D. Barrado, S. Belladitta, V. Belokurov , et al. (240 additional authors not shown)

Abstract: The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline t… ▽ More The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

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

Authors: Théo Lebeau, Stefano Ettori, Nabila Aghanim, Jenny G. Sorce

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

Submitted 27 March, 2024; originally announced March 2024.

Comments: Submitted to A&A on March the 27th 2024, 14 pages, 8 figures and 1 table

Simulating the LOcal Web (SLOW) -- II: Properties of local galaxy clusters

Authors: Elena Hernández-Martínez, Klaus Dolag, Benjamin Seidel, Jenny G. Sorce, Nabila Aghanim, Sergey Pilipenko, Stefan Gottloeber, Théo Lebeau, Milena Valentini

Abstract: This is the second paper in a series presenting the results from a 500 $h^{-1}$Mpc large constrained hydro-dynamical simulation of the local Universe (SLOW). The initial conditions are based on peculiar velocities derived from the CosmicFlows-2 catalogue. The inclusion of galaxy formation treatment, allows to directly predict observable properties of the Intra-Cluster Medium (ICM) within galaxy cl… ▽ More This is the second paper in a series presenting the results from a 500 $h^{-1}$Mpc large constrained hydro-dynamical simulation of the local Universe (SLOW). The initial conditions are based on peculiar velocities derived from the CosmicFlows-2 catalogue. The inclusion of galaxy formation treatment, allows to directly predict observable properties of the Intra-Cluster Medium (ICM) within galaxy clusters. Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts, enables us to assess the effectiveness of the initial condition constraints in accurately replicating the non-linear properties of the largest, collapsed objects within the simulation. Based on the combination of several, publicly available surveys, we identified 45 local Universe galaxy clusters in SLOW, including the 13 most massive from the Planck SZ catalog and 70% of those with $M_{500} > 2\times 10^{14}$ M$_{\odot}$. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature and Compton-y by comparing it to a random selection. In relation to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross identified sample allows to independently deduce a hydrostatic mass bias of $(1-b)\approx0.87$. The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of the galaxy clusters within our local neighbourhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: Submitted to Astronomy & Astrophysics

Journal ref: A&A 687, A253 (2024)

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

Simulating the LOcal Web (SLOW) -- III: Synchrotron Emission from the Local Cosmic Web

Authors: Ludwig M. Böss, Klaus Dolag, Ulrich P. Steinwandel, Elena Hernández-Martínez, Benjamin Seidel, Jenny G. Sorce

Abstract: Aims: Detecting diffuse synchrotron emission from the cosmic web is still a challenge for current radio telescopes. We aim to make predictions for the detectability of cosmic web filaments from simulations. Methods: We present the first cosmological MHD simulation of a 500 $h^{-1} c$Mpc volume with an on-the-fly spectral cosmic ray (CR) model. This allows us to follow the evolution of populations… ▽ More Aims: Detecting diffuse synchrotron emission from the cosmic web is still a challenge for current radio telescopes. We aim to make predictions for the detectability of cosmic web filaments from simulations. Methods: We present the first cosmological MHD simulation of a 500 $h^{-1} c$Mpc volume with an on-the-fly spectral cosmic ray (CR) model. This allows us to follow the evolution of populations of CR electrons and protons within every resolution element of the simulation. We model CR injection at shocks, while accounting for adiabatic changes to the CR population and high energy loss processes of electrons. The synchrotron emission is then calculated from the aged electron population, using the simulated magnetic field, as well as different models for origin and amplification of magnetic fields. We use constrained initial conditions, which closely resemble the local Universe and compare the results of the cosmological volume to zoom-in simulation of the Coma cluster, to study the impact of resolution and turbulent re-acceleration of CRs on the results. Results: We find consistent injection of CRs at accretion shocks onto cosmic web filaments and galaxy clusters. This leads to diffuse emission from filaments of the order $S_ν\approx 0.1 \: μ$Jy beam$^{-1}$ for a potential LOFAR observation at 144 MHz, when assuming the most optimistic magnetic field model and the inclusion of an on-the-fly treatment of re-acceleration of electrons by turbulence. The flux can be increased by up-to two orders of magnitude for different choices of CR injection parameters. This can bring the flux within a factor of 10 of the current limits for direct detection. We find a spectral index of the simulated synchrotron emission from filaments of α {\approx} 1.0 - 1.5. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 23 pages, 12 figures, submitted to A&A. Comments welcome!

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

Authors: Théo Lebeau, Jenny G. Sorce, Nabila Aghanim, Elena Hernández-Martínez, Klaus Dolag

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

Submitted 7 March, 2024; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: Published in A&A, proofs version, 18 pages, 10 figures and 3 tables

Journal ref: Astronomy & Astrophysics, volume 682, page A157, 2024

Statistically bias-minimized peculiar velocity catalogs from Gibbs point processes and Bayesian inference

Authors: Jenny G. Sorce, Radu S. Stoica, Elmo Tempel

Abstract: Galaxy peculiar velocities are excellent cosmological probes provided that biases inherent to their measurements are contained before any study. This paper proposes a new algorithm based on an object point process model whose probability density is built to statistically reduce the effects of Malmquist biases and uncertainties due to lognormal errors in radial peculiar velocity catalogs. More prec… ▽ More Galaxy peculiar velocities are excellent cosmological probes provided that biases inherent to their measurements are contained before any study. This paper proposes a new algorithm based on an object point process model whose probability density is built to statistically reduce the effects of Malmquist biases and uncertainties due to lognormal errors in radial peculiar velocity catalogs. More precisely, a simulated annealing algorithm permits maximizing the probability density describing the point process model. The resulting configurations are bias-minimized catalogs. Tests are conducted on synthetic catalogs mimicking the second and third distance modulus catalogs of the Cosmicflows project from which peculiar velocity catalogs are derived. By reducing the local peculiar velocity variance in catalogs by an order of magnitude, the algorithm permits recovering the expected one while preserving the small-scale velocity correlation. It also permits retrieving the expected clustering. The algorithm is then applied to the observational catalogs. The large-scale structure reconstructed with the Wiener-filter technique applied to the bias-minimized observational catalogs matches with great success the local cosmic web as depicted by redshift surveys of local galaxies. These new bias-minimized versions of peculiar velocity catalogs can be used as a starting point for several studies from possibly estimating the most probable Hubble constant, H0, value to the production of simulations constrained to reproduce the local Universe. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A, 26 pages, 22 figures, 3 tables

Journal ref: A&A 679, A1 (2023)

The cosmic web around the Coma cluster from constrained cosmological simulations: I. Filaments connected to Coma at $z = 0$

Authors: Nicola Malavasi, Jenny G. Sorce, Klaus Dolag, Nabila Aghanim

Abstract: Galaxy clusters in the Universe occupy the important position of nodes of the cosmic web. They are connected among them by filaments, elongated structures composed of dark matter, galaxies, and gas. The connection of galaxy clusters to filaments is important, as it is related to the process of matter accretion onto the former. For this reason, investigating the connections to the cosmic web of mas… ▽ More Galaxy clusters in the Universe occupy the important position of nodes of the cosmic web. They are connected among them by filaments, elongated structures composed of dark matter, galaxies, and gas. The connection of galaxy clusters to filaments is important, as it is related to the process of matter accretion onto the former. For this reason, investigating the connections to the cosmic web of massive clusters, especially well known ones for which a lot of information is available, is a hot topic in astrophysics. In a previous work we performed an analysis of the filament connections of the Coma cluster of galaxies, as detected from the observed galaxy distribution. In this work we resort to a numerical simulation whose initial conditions are constrained to reproduce the Local Universe, including the region of the Coma cluster to interpret our observations in an evolutionary context. We detect the filaments connected to the simulated Coma cluster and perform an accurate comparison with the cosmic web configuration we detected in observations. We perform an analysis of the halos' spatial and velocity distributions close to the filaments in the cluster outskirts. We conclude that, although not significantly larger than the average, the flux of accreting matter on the simulated Coma cluster is significantly more collimated close to the filaments with respect to the general isotropic accretion flux. This paper is the first example of such a result and the first installment in a series of publications which will explore the build-up of the Coma cluster system in connection to the filaments of the cosmic web as a function of redshift. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 20 pages, 13 figures, accepted for publication at A&A

Journal ref: A&A 675, A76 (2023)

Simulating the LOcal Web (SLOW): I. Anomalies in the local density field

Authors: Klaus Dolag, Jenny G. Sorce, Sergey Pilipenko, Elena Hernández-Martínez, Milena Valentini, Stefan Gottlöber, Nabila Aghanim, Ildar Khabibullin

Abstract: Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distributio… ▽ More Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Methods (abridged): We used a 500 Mpc/h large constrained simulation of the LU with initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalogue and follow galaxy formation physics directly in the hydro-dynamical simulations to base the comparison on stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 Mpc/h large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that haloes and galaxies in our constrained simulation trace the local DM density field very differently. Thereby, this simulation reproduces the observed 50% under-density of galaxy clusters and groups within the sphere of ~100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ~1.5$σ$ feature. At the same time, the simulation reproduces the observed over-density of massive galaxy clusters within the same sphere, which on its own also corresponds to a ~1.5$σ$ feature. Interestingly, we find that only 44 out of 15635 random realizations (i.e. 0.28%) are matching both anomalies, making the LU to be a ~3$σ$ environment. We finally compared a mock galaxy catalogue with the observed distribution of galaxies in the LU, finding also a match to the observed factor of two over-density at ~16 Mpc as well as the observed 15% under-density at ~40 Mpc distance. △ Less

Submitted 23 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 16 pages, 12 figure, accepted for publication in A&A

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

Cosmic Variance and the Inhomogeneous UV Luminosity Function of Galaxies During Reionization

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Joseph S. W. Lewis, Dominique Aubert, Jenny G. Sorce, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Romain Teyssier, Gustavo Yepes

Abstract: When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- a… ▽ More When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from $z\simeq6$-15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same $z$. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by $z=6$ (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, `rolling-over' at $M_\text{UV}\gtrsim-17$. CoDa II's UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at $z=6$-8, including the faint end, except for the faintest data point at $z=6$, based on one galaxy at $M_\text{UV}=-12.5$. According to CoDa II, the probability of observing the latter is $\sim5\%$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. We find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2-4, primarily by underestimating the variance in halo abundance when the sample volume is small. △ Less

Submitted 8 August, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 16 pages, 16 figures, accepted by MNRAS 07/20/23, comments welcome

Distortions of the Hubble diagram: Line-of-sight signatures of local galaxy clusters

Authors: Jenny G. Sorce, Roya Mohayaee, Nabila Aghanim, Klaus Dolag, Nicola Malavasi

Abstract: The Universe expansion rate is modulated around local inhomogeneities due to their gravitational potential. Velocity waves are then observed around galaxy clusters in the Hubble diagram. This paper studies them in a ~738 Mpc wide, with 2048^3 particles, cosmological simulation of our cosmic environment (a.k.a. CLONE: Constrained LOcal \& Nesting Environment simulation). For the first time, the sim… ▽ More The Universe expansion rate is modulated around local inhomogeneities due to their gravitational potential. Velocity waves are then observed around galaxy clusters in the Hubble diagram. This paper studies them in a ~738 Mpc wide, with 2048^3 particles, cosmological simulation of our cosmic environment (a.k.a. CLONE: Constrained LOcal \& Nesting Environment simulation). For the first time, the simulation shows that velocity waves that arise in the lines-of-sight of the most massive dark matter halos agree with those observed in local galaxy velocity catalogs in the lines-of-sight of Coma and several other local (Abell) clusters. For the best-constrained clusters such as Virgo and Centaurus, i.e. those closest to us, secondary waves caused by galaxy groups, further into the non-linear regime, also stand out. This match is not utterly expected given that before being evolved into a fully non-linear z=0 state, assuming ΛCDM, CLONE initial conditions are constrained with solely linear theory, power spectrum and highly uncertain and sparse local peculiar velocities. Additionally, Gaussian fits to velocity wave envelopes show that wave properties are tightly tangled with cluster masses. This link is complex though and involves the environment and formation history of the clusters. A proposed metric, measuring the distance between the observed and several re-centered simulated lines-of-sight, waves included, is shown to be capable of providing a tight mass range estimate for massive local clusters. Using machine learning techniques to grasp more thoroughly the complex wave-mass relation, velocity waves could in the near future be used to provide additional and independent mass estimates from galaxy dynamics within large cluster radii. △ Less

Submitted 29 April, 2024; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: Accepted for publications in A&A, 17 pages, 12 figures and 5 tables

Euclid legacy science prospects

Authors: Jenny G. Sorce, Antonino Troja, Isaac Tutusaus

Abstract: With the immense number of images, data, and sources that Euclid will deliver, the consortium will be in a unique position to create/provide/construct legacy catalogues. The latter will have exquisite imaging quality and good near-infrared spectroscopy, with impact on many areas of galaxy science. These proceedings review the prospects and scientific output that Euclid will be able to achieve in a… ▽ More With the immense number of images, data, and sources that Euclid will deliver, the consortium will be in a unique position to create/provide/construct legacy catalogues. The latter will have exquisite imaging quality and good near-infrared spectroscopy, with impact on many areas of galaxy science. These proceedings review the prospects and scientific output that Euclid will be able to achieve in areas of galaxy and active galactic nucleus (AGN) evolution, the local and primeval Universe, studies of the Milky Way and stellar populations, supernovae (SN) and transients, Solar System objects, exoplanets, strong lensing and galaxy clusters. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 6 pages, contribution to the ICHEP2022 conference proceedings accompanying the ''Euclid in a nutshell'' and ''Euclid: performance on main cosmological parameter science'' contributions

Euclid in a nutshell

Authors: Antonino Troja, Isaac Tutusaus, Jenny G. Sorce

Abstract: Euclid is a European Space Agency (ESA) mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey in visible and near-infrared bands, covering approximately 15 000 deg2 of the extragalactic sky in six years. Euclid will be equipped with a 1.2 m diameter Silicon Carbide… ▽ More Euclid is a European Space Agency (ESA) mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey in visible and near-infrared bands, covering approximately 15 000 deg2 of the extragalactic sky in six years. Euclid will be equipped with a 1.2 m diameter Silicon Carbide (SiC) mirror telescope feeding two instruments built by the Euclid Consortium: a high-quality panoramic visible imager and a near-infrared photometer and spectrograph. These proceedings briefly describe the satellite and its instruments, which are optimised for pristine point spread function and reduced stray light, producing very crisp images. Furthermore, we summarise the survey strategy, the global scheduling, and the preparations for the satellite commissioning and the Science Data Centers to produce scientific data. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 6 pages, 3 figures, contribution to the ICHEP 2022 conference proceedings, accompanying the "Euclid: performance on main cosmological parameter science" and "Euclid legacy science prospects" contributions

Euclid: performance on main cosmological parameter science

Authors: Isaac Tutusaus, Jenny G. Sorce, Antonino Troja

Abstract: Euclid will observe 15 000 deg$^2$ of the darkest sky, in regions free of contamination by light from our Galaxy and our Solar System. Three "Euclid Deep Fields" surveys covering around 40 deg$^2$ in total will extend the scientific scope of the mission to the high-redshift Universe. The complete survey will be constituted by hundreds of thousands of images and several tens of petabytes of data. A… ▽ More Euclid will observe 15 000 deg$^2$ of the darkest sky, in regions free of contamination by light from our Galaxy and our Solar System. Three "Euclid Deep Fields" surveys covering around 40 deg$^2$ in total will extend the scientific scope of the mission to the high-redshift Universe. The complete survey will be constituted by hundreds of thousands of images and several tens of petabytes of data. About 10 billion sources will be observed. With these images Euclid will probe the expansion history of the Universe and the evolution of cosmic structures. This will be achieved by measuring the effect on galaxy shapes due to dark matter gravitational lensing, and by reconstructing the three-dimensional distribution of cosmic structures from the measured spectroscopic redshifts of galaxies and clusters of galaxies. These proceedings present the implications for cosmology and cosmological constraints of this unprecedented data set. Of particular interest are the expected constraints on the nature of dark energy. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 6 pages, contribution to the ICHEP 2022 conference proceedings, accompanying the "Euclid in a nutshell" and "Euclid legacy science prospects" contributions

Kinematics of the Local Group gas and galaxies in the Hestia simulations

Authors: Luis Biaus, Sebastián E. Nuza, Philipp Richter, Martin Sparre, Cecilia Scannapieco, Mitali Damle, Jenny G. Sorce, Robert J. J. Grand, Elmo Tempel, Noam I. Libeskind, Maan H. Hani

Abstract: We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of $\sim 10^{12}$~M$_\odot$, their satellites and minor… ▽ More We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of $\sim 10^{12}$~M$_\odot$, their satellites and minor isolated galaxies, all sharing the same large-scale motion within a volume of a few Mpc. We characterise the gas and galaxy kinematics within the simulated LGs, from the perspective of the Sun, to compare with observed trends from recent HST/COS absorption-line observations and LG galaxy data. To analyse the velocity pattern of LG gas and galaxies seen in the observational data, we build sky maps from the local standard of rest, and the galactic and local group barycentre frames. Our findings show that the establishment of a radial velocity dipole at low/high latitudes, near the preferred barycentre direction, is a natural outcome of simulation kinematics for material {\it outside} the Milky Way virial radius after removing galaxy rotation when the two main LG galaxies are approaching. Our results favour a scenario where gas and galaxies stream towards the LG barycentre producing a velocity dipole resembling observations. While our study shows in a qualitative way the global matter kinematics in the LG as part of its on-going assembly, quantitative estimates of gas-flow rates and physical conditions of the LG gas have to await a more detailed modeling of the ionization conditions, which will be presented in a follow-up paper. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

Hermeian haloes in cosmological volumes

Authors: Anastasiia Osipova, Sergey Pilipenko, Stefan Gottlöber, Noam I. Libeskind, Oliver Newton, Jenny G. Sorce, Gustavo Yepes

Abstract: Recent studies based on numerical models of the Local Group predict the existence of field haloes and galaxies that have visited two distinct galaxies in the past, called Hermeian haloes. This work presents an analysis of the Hermeian haloes population in two high-resolution dark matter-only N-body simulations from the MultiDark suit (the ESMDPL, VSMDPL). Hermeian haloes make up from 0.5 to 2.5 pe… ▽ More Recent studies based on numerical models of the Local Group predict the existence of field haloes and galaxies that have visited two distinct galaxies in the past, called Hermeian haloes. This work presents an analysis of the Hermeian haloes population in two high-resolution dark matter-only N-body simulations from the MultiDark suit (the ESMDPL, VSMDPL). Hermeian haloes make up from 0.5 to 2.5 per cent of the total number of field haloes depending on their mass. Furthermore, the results of our study suggest that at a sufficiently high resolution simulation, Hermeian haloes may be found around almost every halo, making them interesting for studies of matter exchange between galaxies. We find that about half (22 out of 49) of the selected Local Group analogues contain Hermeian haloes that passed through the haloes of both the Milky Way and M31 if the distance between the two main haloes is below 1 $h^{-1} \;\mathrm{Mpc}$; this fraction drops to one-fifth (24 out of 108) for distances of up to 1.5 $h^{-1} \;\mathrm{Mpc}$. We confirm earlier findings that unlike other field haloes, Hermeians are grouped along the line connecting the primary hosts of Local Group-like systems, which should facilitate their identification in observations. The vast majority of the Hermeian haloes, whose second target is a Milky Way analogue, are currently moving away from it with increased velocity compared to remaining field halo populations. Interestingly, the obtained data admits that NGC 3109 could have passed through the Andromeda Galaxy and the Milky Way earlier. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 13 pages, 13 figures, submitted to MNRAS

Journal ref: Physics of the Dark Universe 42 (2023) 101328

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

Reionization time of the Local Group and Local-Group-like halo pairs

Authors: Jenny G. Sorce, Pierre Ocvirk, Dominique Aubert, Stefan Gottloeber, Paul R. Shapiro, Taha Dawoodbhoy, Gustavo Yepes, Kyungjin Ahn, Ilian T. Iliev, Joseph S. W. Lewis

Abstract: Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper… ▽ More Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS, 11 pages, 5 figures

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)

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

Authors: Hyunbae Park, Hyo Jeong Kim, Kyungjin Ahn, Hyunmi Song, Intae Jung, Pierre Ocvirk, Paul R. Shapiro, Taha Dawoodbhoy, Jenny G. Sorce, Ilian T. Iliev

Abstract: During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previ… ▽ More During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previous literature. Then, we run our code on the simulated IGM of the CoDaII simulation. We find that the scattered light can leave an observable imprint on the emergent spectrum if collected over several square arcminutes. Scattering can redden the emission by increasing the path lengths of photons, but it can also make the photons bluer by upscattering them according to the peculiar motion of the scatterer. The photons emitted on the far blue side of the resonance appear more extended in both frequency and space compared to those emitted near the resonance. This provides a discriminating feature for the blueward emission, which cannot be constrained from the unscattered light coming directly from the source. The ionization state of the IGM also affects the scattered light spectrum. When the source is in a small HII region, the emission goes through more scatterings in the surrounding HI region regardless of the initial frequency and ends up more redshifted and spatially extended. This can result in a weakening of the scattered light toward high $z$ during reionization. Our results provide a framework for interpreting the scattered light to be measured by high-$z$ integral-field-unit surveys. △ Less

Submitted 16 May, 2022; v1 submitted 13 February, 2022; originally announced February 2022.

Comments: Accepted for publication in the ApJ

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Jenny G. Sorce, Yohan Dubois, Dominique Aubert, Luke Conaboy, Paul R. Shapiro, Taha Dawoodbhoy, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Yann Rasera, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Émilie Thélie

Abstract: Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) a… ▽ More Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$). △ Less

Submitted 22 August, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Submitted to MNRAS Accepted 2022 August 08. Received 2022 August 05; in original form 2022 February 11

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

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

I- A hydrodynamical CLONE of the Virgo cluster of galaxies to confirm observationally-driven formation scenarios

Authors: Jenny G. Sorce, Yohan Dubois, Jeremy Blaizot, Sean L. McGee, Gustavo Yepes, Alexander Knebe

Abstract: At ~16-17Mpc from us, the Virgo cluster is a formidable source of information to study cluster formation and galaxy evolution in rich environments. Several observationally-driven formation scenarios arose within the past decade to explain the properties of galaxies that entered the cluster recently and the nature of the last significant merger that the cluster underwent. Confirming these scenarios… ▽ More At ~16-17Mpc from us, the Virgo cluster is a formidable source of information to study cluster formation and galaxy evolution in rich environments. Several observationally-driven formation scenarios arose within the past decade to explain the properties of galaxies that entered the cluster recently and the nature of the last significant merger that the cluster underwent. Confirming these scenarios requires extremely faithful numerical counterparts of the cluster. This paper presents the first CLONE, Constrained LOcal and Nesting Environment, simulation of the Virgo cluster within a ~15Mpc radius sphere. This cosmological hydrodynamical simulation, with feedback from supernovae and active galactic nuclei, with a ~3x10^7Msun dark matter particle mass and a minimum cell size of 350pc in the zoom region, reproduces Virgo within its large scale environment unlike a random cluster simulation. Overall the distribution of the simulated galaxy population matches the observed one including M87. The simulated cluster formation reveals exquisite agreements with observationally-driven scenarios: within the last Gigayear, about 300 small galaxies (M*>10^7Msun) entered the cluster, most of them within the last 500Myr. The last significant merger event occurred about 2 Gigayears ago: a group with a tenth of the mass of today's cluster entered from the far side as viewed from the Milky Way. This excellent numerical replica of Virgo will permit studying different galaxy type evolution (jellyfish, backsplash, etc.) as well as feedback phenomena in the cluster core via unbiased comparisons between simulated and observed galaxies and hot gas phase profiles to understand this great physics laboratory. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: Accepted for publication in MNRAS, 15 pages, 14 figures, 4 tables

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 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

The Past History of Galaxy Clusters told by their present neighbors

Authors: Jenny G. Sorce, Stefan Gottlöber, Gustavo Yepes

Abstract: Galaxy clusters can play a key role in modern cosmology provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of c… ▽ More Galaxy clusters can play a key role in modern cosmology provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass - environment relations. This paper presents a mass-free correlation between the present neighbor distribution of cluster-size halos and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random halos with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size halos to 200 Virgo-like halos extracted from simulations that resemble the local Universe. The number of neighbors at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like halos, clusters with numerous neighbors within a distance of about 2 times their virial radius experience a transition at z~1 between an active period of mass accretion, relative to the mean, and a quiet history. On the contrary, clusters with few neighbors share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about 4 times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbors provides invaluable insights into the past history of these objects. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS, 10 pages, 5 figures, 4 tables

An optimised tiling pattern for multi-object spectroscopic surveys: application to the 4MOST survey

Authors: E. Tempel, T. Tuvikene, M. M. Muru, R. S. Stoica, T. Bensby, C. Chiappini, N. Christlieb, M. -R. L. Cioni, J. Comparat, S. Feltzing, I. Hook, A. Koch, G. Kordopatis, M. Krumpe, J. Loveday, I. Minchev, P. Norberg, B. F. Roukema, J. G. Sorce, J. Storm, E. Swann, E. N. Taylor, G. Traven, C. J. Walcher, R. S. de Jong

Abstract: Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and ext… ▽ More Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: 20 pages, 11 figures, submitted to MNRAS

Efficiently estimating mean, uncertainty and unconstrained large scale fraction of local Universe simulations with paired fixed fields

Authors: Jenny G. Sorce

Abstract: Provided a random realization of the cosmological model, observations of our cosmic neighborhood now allow us to build simulations of the latter down to the non-linear threshold. The resulting local Universe models are thus accurate up to a given residual cosmic variance. Namely some regions and scales are apparently not constrained by the data and seem purely random. Drawing conclusions together… ▽ More Provided a random realization of the cosmological model, observations of our cosmic neighborhood now allow us to build simulations of the latter down to the non-linear threshold. The resulting local Universe models are thus accurate up to a given residual cosmic variance. Namely some regions and scales are apparently not constrained by the data and seem purely random. Drawing conclusions together with their uncertainties involves then statistics implying a considerable amount of computing time. By applying the constraining algorithm to paired fixed fields, this paper diverts the original techniques from their first use to efficiently disentangle and estimate uncertainties on local Universe simulations obtained with random fields. Paired fixed fields differ from random realizations in the sense that their Fourier mode amplitudes are fixed and they are exactly out of phase. Constrained paired fixed fields show that only 20% of the power spectrum on large scales (> tens of megaparsecs) is purely random. Namely 80% of it is partly constrained by the large scale / small scale data correlations. Additionally, two realizations of our local environment obtained with paired fixed fields of the same pair constitute an excellent non-biased average or quasi-linear realization of the latter, namely the equivalent of hundreds of constrained simulations. The variance between these two realizations gives the uncertainty on the achievable local Universe simulations. These two simulations will permit enhancing faster our local cosmic web understanding thanks to a drastically reduced required computational time to appreciate its modeling limits and uncertainties. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS, 11 pages, 9 figures, 2 tables

Lyman-alpha transmission properties of the intergalactic medium in the CoDaII simulation

Authors: Max Gronke, Pierre Ocvirk, Charlotte Mason, Jorryt Matthee, Sarah E. I. Bosman, Jenny G. Sorce, Joseph Lewis, Kyungjin Ahn, Dominique Aubert, Taha Dawoodbhoy, Ilian T. Iliev, Paul R. Shapiro, Gustavo Yepes

Abstract: The decline in abundance of Lyman-$α$ (Ly$α$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$.… ▽ More The decline in abundance of Lyman-$α$ (Ly$α$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$α$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$α$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 1\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$α$ peak. We discuss our results in the context of the interpretation of observations. △ Less

Submitted 21 October, 2021; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: 13 pages, 12 figures. Published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3697-3709

Probabilistic fibre-to-target assignment algorithm for multi-object spectroscopic surveys

Authors: E. Tempel, P. Norberg, T. Tuvikene, T. Bensby, C. Chiappini, N. Christlieb, M. -R. L. Cioni, J. Comparat, L. J. M. Davies, G. Guiglion, A. Koch, G. Kordopatis, M. Krumpe, J. Loveday, A. Merloni, G. Micheva, I. Minchev, B. F. Roukema, J. G. Sorce, E. Starkenburg, J. Storm, E. Swann, W. F. Thi, G. Traven, R. S. de Jong

Abstract: Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We presen… ▽ More Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys. △ Less

Submitted 4 February, 2020; v1 submitted 25 January, 2020; originally announced January 2020.

Comments: 19 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 635, A101 (2020)

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

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

Abstract: Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmissio… ▽ More Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8. △ Less

Submitted 18 June, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Submitted to MNRAS on 09/01/2020. Accepted on 12/06/2020

Virgo: an unlikely cluster of galaxies because of its environment

Authors: Jenny G. Sorce, Jérémy Blaizot, Yohan Dubois

Abstract: Galaxy clusters constitute powerful cosmological probes thanks to comparisons between observed and simulated clusters. As such Virgo constitutes a formidable source for detailed observations facilitated by its proximity. However, the diversity of clusters complicates the comparisons on a one-to-one basis. Simulated clusters must be carefully selected, a daunting task since most properties are unkn… ▽ More Galaxy clusters constitute powerful cosmological probes thanks to comparisons between observed and simulated clusters. As such Virgo constitutes a formidable source for detailed observations facilitated by its proximity. However, the diversity of clusters complicates the comparisons on a one-to-one basis. Simulated clusters must be carefully selected, a daunting task since most properties are unknown. Alternatively, lookalikes produced in the proper large scale environment can be used. Additionally, their statistical study give access to the mean properties of the observed cluster including its most probable history as well as its deviation from an average cluster. This paper presents such a statistical study with 200 Virgo-like and 400+ cluster-size random dark matter halos. Only 18%(0.5%) of these random halos comply within 3(2)-sigma with the mean values (radius, velocity dispersion, number of substructures, spin, velocity, concentration, center of mass offset with respect to the spherical center) of Virgo halos at z=0 and abide by a similar merging history up to redshift 4. None are within 1-sigma because of environmentally induced properties (number of substructures and velocity). For further comparisons, random halos are selected to reproduce the mass distribution of the lookalikes to cancel mass bias effects. Redshift 1 appears then as a turning point: random to Virgo-like property ratios are alternatively smaller/larger than 1. This highlights the importance of studying clusters within their proper large scale environment: simulated galaxy population, grandly affected by the cluster history, can then be compared with the observed one in details. Direct lookalikes simplify grandly the challenge. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: Accepted for publication in MNRAS, 12 pages, 7 figures, 5 tables

4MOST Survey Strategy Plan

Authors: G. Guiglion, C. Battistini, C. P. M. Bell, T. Bensby, T. Boller, C. Chiappini, J. Comparat, N. Christlieb, R. Church, M. -R. L. Cioni, L. Davies, T. Dwelly, R. S. de Jong, S. Feltzing, A. Gueguen, L. Howes, M. Irwin, I. Kushniruk, M. I Lam, J. Liske, R. McMahon, A. Merloni, P. Norberg, A. S. G. Robotham, O. Schnurr , et al. (8 additional authors not shown)

Abstract: The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surv… ▽ More The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 17

Propagation of UHECRs in the local Universe and origin of cosmic magnetic fields

Authors: Stefan Hackstein, Franco Vazza, Marcus Brüggen, Jenny G. Sorce, Stefan Gottlöber

Abstract: We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We investigate the impact of different magneto-genesis scenarios, both, primordial and astrophysical, on the propagation of cosmic rays. Our study shows that different scenarios of magneto-genesis do not have a large i… ▽ More We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We investigate the impact of different magneto-genesis scenarios, both, primordial and astrophysical, on the propagation of cosmic rays. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements. The distribution of nearby sources causes anisotropy at very high energies, independent of the magnetic field model. We compare our results to the dipole signal measured by the Pierre Auger Observatory. All our models could reproduce the observed dipole amplitude with a pure iron injection composition. This is due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of 4 - 8 EeV. △ Less

Submitted 2 February, 2019; originally announced February 2019.

Comments: part of proceedings of XXXth IAU in Vienna, 2018, Focus Meeting 8

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

The Three Hundred project: a large catalogue of theoretically modelled galaxy clusters for cosmological and astrophysical applications

Authors: Weiguang Cui, Alexander Knebe, Gustavo Yepes, Frazer Pearce, Chris Power, Romeel Dave, Alexander Arth, Stefano Borgani, Klaus Dolag, Pascal Elahi, Robert Mostoghiu, Giuseppe Murante, Elena Rasia, Doris Stoppacher, Jesus Vega-Ferrero, Yang Wang, Xiaohu Yang, Andrew Benson, Sofía A. Cora, Darren J. Croton, Manodeep Sinha, Adam R. H. Stevens, Cristian A. Vega-Martínez, Jake Arthur, Anna S. Baldi , et al. (12 additional authors not shown)

Abstract: We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic… ▽ More We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 20 pages, 8 figures, 7 tables. MNRAS published version

Journal ref: Cui, W., Knebe, A., Yepes, G., et al.\ 2018, \mnras, 480, 2898

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

Galaxy clusters in local Universe simulations without density constraints: a long uphill struggle

Authors: Jenny G. Sorce

Abstract: Galaxy clusters are excellent cosmological probes provided that their formation and evolution within the large scale environment are precisely understood. Therefore studies with simulated galaxy clusters have flourished. However detailed comparisons between simulated and observed clusters and their population - the galaxies - are complicated by the diversity of clusters and their surrounding envir… ▽ More Galaxy clusters are excellent cosmological probes provided that their formation and evolution within the large scale environment are precisely understood. Therefore studies with simulated galaxy clusters have flourished. However detailed comparisons between simulated and observed clusters and their population - the galaxies - are complicated by the diversity of clusters and their surrounding environment. An original way initiated by Bertschinger as early as 1987, to legitimize the one-to-one comparison exercise down to the details, is to produce simulations constrained to resemble the cluster under study within its large scale environment. Subsequently several methods have emerged to produce simulations that look like the local Universe. This paper highlights one of these methods and its essential steps to get simulations that not only resemble the local Large Scale Structure but also that host the local clusters. It includes a new modeling of the radial peculiar velocity uncertainties to remove the observed correlation between the decreases of the simulated cluster masses and of the amount of data used as constraints with the distance from us. This method has the particularity to use solely radial peculiar velocities as constraints: no additional density constraints are required to get local cluster simulacra. The new resulting simulations host dark matter halos that match the most prominent local clusters such as Coma. Zoom-in simulations of the latter and of a volume larger than the 30 Mpc/h radius inner sphere become now possible to study local clusters and their effects. Mapping the local Sunyaev-Zel'dovich and Sachs-Wolfe effects can follow. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: Accepted for publication in MNRAS, 9 pages, 6 figures, 2 tables

Bayesian group finder based on marked point processes. Method and feasibility study using the 2MRS data set

Authors: Elmo Tempel, Maarja Kruuse, Rain Kipper, Taavi Tuvikene, Jenny G. Sorce, Radu S. Stoica

Abstract: Galaxy groups and clusters are formidable cosmological probes. They permit the studying of the environmental effects on galaxy formation. A reliable detection of galaxy groups is an open problem and is important for ongoing and future cosmological surveys. We propose a probabilistic galaxy group detection algorithm based on marked point processes with interactions. The pattern of galaxy groups… ▽ More Galaxy groups and clusters are formidable cosmological probes. They permit the studying of the environmental effects on galaxy formation. A reliable detection of galaxy groups is an open problem and is important for ongoing and future cosmological surveys. We propose a probabilistic galaxy group detection algorithm based on marked point processes with interactions. The pattern of galaxy groups in a catalogue is seen as a random set of interacting objects. The positions and the interactions of these objects are governed by a probability density. The estimator of the unknown cluster pattern is given by the configuration of objects maximising the proposed probability density. Adopting the Bayesian framework, the proposed probability density is maximised using a simulated annealing algorithm. The method provides "for free" additional information such as the probabilities that a point or two points in the observation domain belong to the cluster pattern, respectively. These supplementary tools allow the construction of tests and techniques to validate and to refine the detection result. To test the feasibility of the proposed methodology, we applied it to the well-studied 2MRS data set. Compared to previously published Friends-of-Friends (FoF) group finders, the proposed Bayesian group finder gives overall similar results. The proposed Bayesian group finder is tested on a galaxy redshift survey, but more detailed analyses are needed to understand the actual capabilities of the algorithm regarding upcoming cosmological surveys. The presented mathematical framework permits adapting it easily for other data sets (in astronomy and in other fields of sciences). In cosmology, one promising application is the detection of galaxy groups in photometric galaxy redshift surveys, while taking into account the full photometric redshift posteriors. (abridged) △ Less

Submitted 12 July, 2018; v1 submitted 12 June, 2018; originally announced June 2018.

Comments: 19 pages, 14 figures, accepted by A&A

Journal ref: A&A 618, A81 (2018)

Merger types forming the Virgo cluster in recent gigayears

Authors: Mark Olchanski, Jenny G. Sorce

Abstract: As our closest cluster-neighbor, the Virgo cluster of galaxies is intensely studied by observers to unravel the mysteries of galaxy evolution within clusters. At this stage, cosmological numerical simulations of the cluster are useful to efficiently test theories and calibrate model. However, it is not trivial to select the perfect simulacrum of the Virgo cluster to fairly compare in detail its ob… ▽ More As our closest cluster-neighbor, the Virgo cluster of galaxies is intensely studied by observers to unravel the mysteries of galaxy evolution within clusters. At this stage, cosmological numerical simulations of the cluster are useful to efficiently test theories and calibrate model. However, it is not trivial to select the perfect simulacrum of the Virgo cluster to fairly compare in detail its observed and simulated galaxy populations that are affected by the type and history of the cluster. Determining precisely the properties of Virgo for a later selection of simulated clusters becomes essential. It is still not clear how to access some of these properties such as the past history of the Virgo cluster from current observations. Therefore, directly producing effective simulacra of the Virgo cluster is inevitable. Efficient simulacra of the Virgo cluster can be obtained via simulations that resemble the local Universe down to the cluster scale. In such simulations, Virgo-like halos form in the proper local environment and permit assessing the most probable formation history of the cluster. Studies based on these simulations have already revealed that the Virgo cluster has had a quiet merging history over the last seven gigayears and that the cluster accretes matter along a preferential direction. This paper reveals that in addition such Virgo halos have had on average only one merger larger than about a tenth of their mass at redshift zero within the last four gigayears. This second branch (by opposition to main branch) formed in a given sub-region and merged recently (within the last gigayear). These properties are not shared with a set of random halos within the same mass range. This study extends [truncated] △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: Accepted for publication in A&A, 10 pages, 8 figures

Galaxy clusters in simulations of the local Universe: a matter of constraints

Authors: Jenny G. Sorce, Elmo Tempel

Abstract: To study the full formation and evolution history of galaxy clusters and their population, high resolution simulations of the latter are flourishing. However comparing observed clusters to the simulated ones on a one-to-one basis to refine the models and theories down to the details is non trivial. The large variety of clusters limits the comparisons between observed and numerical clusters. Simula… ▽ More To study the full formation and evolution history of galaxy clusters and their population, high resolution simulations of the latter are flourishing. However comparing observed clusters to the simulated ones on a one-to-one basis to refine the models and theories down to the details is non trivial. The large variety of clusters limits the comparisons between observed and numerical clusters. Simulations resembling the local Universe down to the cluster scales permit pushing the limit. Simulated and observed clusters can be matched on a one-to-one basis for direct comparisons provided that clusters are well reproduced besides being in the proper large scale environment. Comparing random and local-Universe like simulations obtained with differently grouped observational catalogs of peculiar velocities, this paper shows that the grouping scheme used to remove non-linear motions in the catalogs that constrain the simulations affects the quality of the numerical clusters. With a less aggressive grouping scheme - galaxies still falling onto clusters are preserved - combined with a bias minimization scheme, the mass of the dark matter halos, simulacra for 5 local clusters - Virgo, Centaurus, Coma, Hydra and Perseus - is increased by 39% closing the gap with observational mass estimates. Simulacra are found on average in 89% of the simulations, an increase of 5% with respect to the previous grouping scheme. The only exception is Perseus. Since the Perseus-Pisces region is not well covered by the used peculiar velocity catalog, the latest release let us foresee a better simulacrum for Perseus in a near future. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS, 9 pages, 5 figures, 1 table

Simulations of ultra-high Energy Cosmic Rays in the local Universe and the origin of Cosmic Magnetic Fields

Authors: Stefan Hackstein, Franco Vazza, Marcus Brüggen, Jenny G. Sorce, Stefan Gottlöber

Abstract: We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We use constrained initial conditions with the cosmological magnetohydrodynamics code {\sc ENZO}. The resulting models of the distribution of magnetic fields in the local Universe are used in the \crpropa code to simul… ▽ More We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We use constrained initial conditions with the cosmological magnetohydrodynamics code {\sc ENZO}. The resulting models of the distribution of magnetic fields in the local Universe are used in the \crpropa code to simulate the propagation of ultra-high energy cosmic rays. We investigate the impact of six different magneto-genesis scenarios, both primordial and astrophysical, on the propagation of cosmic rays over cosmological distances. Moreover, we study the influence of different source distributions around the Milky Way. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements of ultra-high energy cosmic rays. However, at high energies above the GZK-limit, there is anisotropy caused by the distribution of nearby sources, independent of the magnetic field model. This provides a chance to identify cosmic ray sources with future full-sky measurements and high number statistics at the highest energies. Finally, we compare our results to the dipole signal measured by the Pierre Auger Observatory. All our source models and magnetic field models could reproduce the observed dipole amplitude with a pure iron injection composition. Our results indicate that the dipole is observed due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of $4-8 \rm\ EeV$. △ Less

Submitted 5 January, 2018; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 12 pages, 9 figures