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

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

Revisiting the large-scale CMB anomalies: The impact of the SZ signal from the Local Universe

Authors: Gabriel Jung, Nabila Aghanim, Jenny Sorce, Benjamin Seidel, Klaus Dolag, Marian Douspis

Abstract: The full sky measurements of the Cosmic Microwave Background (CMB) temperature anisotropies by $\textit{WMAP}$ and $\textit{Planck}$ have highlighted the presence of several unexpected isotropy-breaking features on the largest angular scales. In this work, we investigate the impact of the local large-scale structure on these anomalies through the thermal and kinetic Sunyaev-Zeldovich effects. We… ▽ More The full sky measurements of the Cosmic Microwave Background (CMB) temperature anisotropies by $\textit{WMAP}$ and $\textit{Planck}$ have highlighted the presence of several unexpected isotropy-breaking features on the largest angular scales. In this work, we investigate the impact of the local large-scale structure on these anomalies through the thermal and kinetic Sunyaev-Zeldovich effects. We use a constrained hydrodynamical simulation that reproduces the local Universe in a box of $500\,h^{-1}\,$Mpc to construct full sky maps of the temperature anisotropies produced by these two CMB secondary effects and discuss their statistical properties on large angular scales. We show the significant role played by the Virgo cluster on these scales, and compare it to theoretical predictions and random patches of the universe obtained from the hydrodynamical simulation $\textit{Magneticum}$. We explore three of the main CMB large-scale anomalies -- i.e., lack of correlation, quadrupole-octopole alignment and hemispherical asymmetry -- , both in the latest $\textit{Planck}$ data (PR4), where they are detected at a similar level to the previous releases, and using the simulated secondaries from the local Universe, verifying their negligible impact. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 15 pages, 20 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 2 September, 2024; v1 submitted 27 March, 2024; originally announced March 2024.

Comments: Published in A&A, proofs version, 14 pages, 10 figures and 1 table

Journal ref: A&A 689, A19 (2024)

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

The undiscovered ultra-diffuse galaxies of the Local Group

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

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

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

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

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

Cosmic filaments delay quenching inside clusters

Authors: S. Kotecha, C. Welker, Z. Zhou, J. Wadsley, K. Kraljic, J. Sorce, E. Rasia, I. Roberts, M. Gray, G. Yepes, W. Cui

Abstract: We investigate how large-scale cosmic filaments impact the quenching of galaxies within one virial radius of 324 simulated clusters from The Three Hundred project. We track cosmic filaments with the versatile, observation-friendly program DisPerSE and identify halos hosting galaxies with VELOCIRaptor. We find that cluster galaxies close to filaments tend to be more star-forming, bluer, and contain… ▽ More We investigate how large-scale cosmic filaments impact the quenching of galaxies within one virial radius of 324 simulated clusters from The Three Hundred project. We track cosmic filaments with the versatile, observation-friendly program DisPerSE and identify halos hosting galaxies with VELOCIRaptor. We find that cluster galaxies close to filaments tend to be more star-forming, bluer, and contain more cold gas than their counterparts further away from filaments. This effect is recovered at all stellar masses. This is in stark contrast with galaxies residing outside of clusters, where galaxies close to filaments show clear signs of density related pre-processing. We first show that the density contrast of filaments is reduced inside the intra-cluster medium. Moreover, examination of flows around and into cluster galaxies shows that the gas flows in intra-cluster filaments are colder and tend to stream along with galaxies in their midst, partially shielding them from strangulation by the hot, dense intra-cluster medium. This also preserves accretion onto satellites and limit ram pressure. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 20 pages, 21 figures, submitted to MNRAS

NIHAO-LG: The uniqueness of Local Group dwarf galaxies

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

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

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

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

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

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

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

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

Submitted 17 September, 2019; originally announced September 2019.

Comments: 9 pages, 3 figures, 1 table

The Origin of Lopsided Satellite Galaxy Distribution in Galaxy Pairs

Authors: Chen Chris Gong, Noam I Libeskind, Elmo Tempel, Quan Guo, Stefan Gottloeber, Gustavo Yepes, Peng Wang, Jenny Sorce, Marcel Pawlowski

Abstract: It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $Λ$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromed… ▽ More It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $Λ$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $Λ$CDM simulations indicating that such a set up is not inconsistent with $Λ$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events. △ Less

Submitted 18 June, 2019; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 10 pages, 11 figures

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 Consortium Survey 8: Cosmology Redshift Survey (CRS)

Authors: J. Richard, J. -P. Kneib, C. Blake, A. Raichoor, J. Comparat, T. Shanks, J. Sorce, M. Sahlén, C. Howlett, E. Tempel, R. McMahon, M. Bilicki, B. Roukema, J. Loveday, D. Pryer, T. Buchert, C. Zhao

Abstract: The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects ov… ▽ More The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission. △ 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 the 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 50

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

4MOST Scientific Operations

Authors: C. J. Walcher, M. Banerji, C. Battistini, C. P. M. Bell, O. Bellido-Tirado, T. Bensby, J. M. Bestenlehner, T. Boller, J. Brynnel, A. Casey, C. Chiappini, N. Christlieb, R. Church, M. -R. L. Cioni, S. Croom, J. Comparat, L. J. M. Davies, R. S. de Jong, T. Dwelly, H. Enke, S. Feltzing, D. Feuillet, M. Fouesneau, D. Ford, S. Frey , et al. (43 additional authors not shown)

Abstract: The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature… ▽ More The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent. △ 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, 12

4MOST: Project overview and information for the First Call for Proposals

Authors: R. S. de Jong, O. Agertz, A. Agudo Berbel, J. Aird, D. A. Alexander, A. Amarsi, F. Anders, R. Andrae, B. Ansarinejad, W. Ansorge, P. Antilogus, H. Anwand-Heerwart, A. Arentsen, A. Arnadottir, M. Asplund, M. Auger, N. Azais, D. Baade, G. Baker, S. Baker, E. Balbinot, I. K. Baldry, M. Banerji, S. Barden, P. Barklem , et al. (313 additional authors not shown)

Abstract: We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolut… ▽ More We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = λ/Δλ\sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations. △ 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, 3

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