Showing 1–50 of 50 results for author: Iorio, G

Tidal mass loss in the Fornax dwarf spheroidal galaxy through N-body simulations with Gaia EDR3-based orbits

Authors: Pierfrancesco Di Cintio, Giuliano Iorio, Francesco Calura, Carlo Nipoti, Marcello Cantari

Abstract: The Fornax dwarf spheroidal galaxy (dSph) represents a challenge for some globular cluster (GC) formation models, because an exceptionally high fraction of its stellar mass is locked in its GC system. In order to shed light on our understanding of GC formation, we aim to constrain the amount of stellar mass that Fornax has lost via tidal interaction with the Milky Way (MW). Exploiting the flexibil… ▽ More The Fornax dwarf spheroidal galaxy (dSph) represents a challenge for some globular cluster (GC) formation models, because an exceptionally high fraction of its stellar mass is locked in its GC system. In order to shed light on our understanding of GC formation, we aim to constrain the amount of stellar mass that Fornax has lost via tidal interaction with the Milky Way (MW). Exploiting the flexibility of effective multi-component $N$-body simulations and relying on state-of-the-art estimates of Fornax's orbital parameters, we study the evolution of the mass distribution of the Fornax dSph in observationally justified orbits in the gravitational potential of the MW over 12 Gyr. We find that, though the dark-matter mass loss can be substantial, the fraction of stellar mass lost by Fornax to the MW is always negligible, even in the most eccentric orbit considered. We conclude that stellar-mass loss due to tidal stripping is not a plausible explanation for the unexpectedly high stellar mass of the GC system of the Fornax dSph and we discuss quantitatively the implications for GC formation scenarios. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 13 pages, 11 figures, Submitted. Comments welcome

Binary black hole mergers from Population III star clusters

Authors: Benedetta Mestichelli, Michela Mapelli, Stefano Torniamenti, Manuel Arca Sedda, Marica Branchesi, Guglielmo Costa, Giuliano Iorio, Filippo Santoliquido

Abstract: Binary black holes (BBHs) born from the evolution of Population III (Pop. III) stars are one of the main high-redshift targets for next-generation ground-based gravitational-wave (GW) detectors. Their predicted initial mass function and lack of metals make them the ideal progenitors of black holes above the upper edge of the pair-instability mass gap, i.e. with a mass higher than $\approx{}134$ (2… ▽ More Binary black holes (BBHs) born from the evolution of Population III (Pop. III) stars are one of the main high-redshift targets for next-generation ground-based gravitational-wave (GW) detectors. Their predicted initial mass function and lack of metals make them the ideal progenitors of black holes above the upper edge of the pair-instability mass gap, i.e. with a mass higher than $\approx{}134$ (241) M$_\odot$ for stars that become (do not become) chemically homogeneous during their evolution. Here, we investigate the effects of cluster dynamics on the mass function of BBHs born from Pop. III stars, by considering the main uncertainties on Pop. III star mass function, orbital properties of binary systems, star cluster's mass and disruption time. In our dynamical models, at least $\sim$5% and up to 100% BBH mergers in Pop. III star clusters have primary mass $m_1$ above the upper edge of the pair-instability mass gap. In contrast, only $\lesssim {} 3$% isolated BBH mergers have primary mass above the gap, unless their progenitors evolved as chemically homogeneous stars. The lack of systems with primary and/or secondary mass inside the gap defines a zone of avoidance with sharp boundaries in the primary mass - mass ratio plane. Finally, we estimate the merger rate density of BBHs and, in the most optimistic case, we find a maximum of $\mathcal{R}\approx200\,{\rm Gpc^{-3}\,yr^{-1}}$ at $z\sim15$ for BBHs formed via dynamical capture. For comparison, the merger rate density of isolated Pop. III BBHs is $\mathcal{R}\leq{}10\,{\rm Gpc^{-3}\,yr^{-1}}$, for the same model of Pop. III star formation history. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 17 pages, 12 figures, 4 tables. Comments welcome

A 3D view on the local gravitational instability of cold gas discs in star-forming galaxies at $0 \lesssim \mathrm{z} \lesssim 5$

Authors: C. Bacchini, C. Nipoti, G. Iorio, F. Roman-Oliveira, F. Rizzo, P. E. Mancera Piña, A. Marasco, A. Zanella, F. Lelli

Abstract: Local gravitational instability (LGI) is considered crucial for regulating star formation and gas turbulence in galaxy discs, especially at high redshift. Instability criteria usually assume infinitesimally thin discs or rely on approximations to include the stabilising effect of the gas disc thickness. We test a new 3D instability criterion for rotating gas discs that are vertically stratified in… ▽ More Local gravitational instability (LGI) is considered crucial for regulating star formation and gas turbulence in galaxy discs, especially at high redshift. Instability criteria usually assume infinitesimally thin discs or rely on approximations to include the stabilising effect of the gas disc thickness. We test a new 3D instability criterion for rotating gas discs that are vertically stratified in an external potential. This criterion reads $Q_{\rm3D}<1$, where $Q_{\rm3D}$ is the 3D analogue of the Toomre parameter $Q$. The advantage of $Q_{\rm3D}$ is that it allows us to study LGI in and above the galaxy midplane in a rigorous and self-consistent way. We apply the criterion to a sample of 44 star-forming galaxies at $0\lesssim\mathrm{z}\lesssim5$ hosting rotating discs of cold gas. The sample is representative of galaxies on the main sequence at $\mathrm{z}\approx 0$ and includes massive star-forming and starburst galaxies at $1\lesssim\mathrm{z}\lesssim5$. For each galaxy, we first apply the Toomre criterion for infinitesimally thin discs, finding 10 unstable systems. We then obtain maps of $Q_{\rm 3D}$ from a 3D model of the gas disc derived in the combined potential of dark matter, stars and the gas itself. According to the 3D criterion, two galaxies with $Q<1$ show no evidence of instability and the unstable regions that are 20% smaller than those where $Q<1$. No unstable disc is found at $0\lesssim\mathrm{z}\lesssim 1$, while $\approx 60$% of the systems at $2\lesssim\mathrm{z}\lesssim5$ are locally unstable. In these latter, a relatively small fraction of the total gas ($\approx 30$%) is potentially affected by the instability. Our results disfavour LGI as the main regulator of star formation and turbulence in moderately star-forming galaxies in the present-day Universe. LGI likely becomes important at high redshift, but the input by other mechanisms seems required [abridged] △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 19 pages (6 figures, 2 tables) + appendixes, accepted by Astronomy & Astrophysics

The boring history of Gaia BH3 from isolated binary evolution

Authors: Giuliano Iorio, Stefano Torniamenti, Michela Mapelli, Marco Dall'Amico, Alessandro A. Trani, Sara Rastello, Cecilia Sgalletta, Stefano Rinaldi, Guglielmo Costa, Bera A. Dhal-Lahtinen, Gaston J. Escobar, Erika Korb, M. Paola Vaccaro, Elena Lacchin, Benedetta Mestichelli, Ugo Niccolò di Carlo, Mario Spera, Manuel Arca Sedda

Abstract: Gaia BH3 is the first observed dormant black hole (BH) with a mass of $\approx{30}$ M$_\odot$ and represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3 focusing on the old and metal-poor stellar population of the Milky Way halo. We use our open-source population synthesis code SEVN to evolve… ▽ More Gaia BH3 is the first observed dormant black hole (BH) with a mass of $\approx{30}$ M$_\odot$ and represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3 focusing on the old and metal-poor stellar population of the Milky Way halo. We use our open-source population synthesis code SEVN to evolve $3.2 \times 10^8$ binaries exploring 16 sets of parameters. We find that systems like Gaia BH3 form preferentially from binaries initially composed of a massive star ($40-60$ M$_\odot$) and a low mass companion ($<1$ M$_\odot$) in a wide ($P>10^3$ days) and eccentric orbit ($e>0.6$). Such progenitors do not undergo any Roche-lobe overflow episode during their entire evolution, so that the final orbital properties of the BH-star system are essentially determined at the core collapse of the primary star. Low natal kicks ($\approx$ 10~km/s) significantly favour the formation of Gaia BH3-like systems, but high velocity kicks up to $\approx 220$ km/s are also allowed. We estimate the formation efficiency for Gaia BH3-like systems in old ($t>10$ Gyr) and metal-poor ($Z<0.01$) populations to be $4 \times 10^{-8}$ M$_\odot^{-1}$, representing $\approx 3\%$ of the whole simulated BH-star population. We expect up to $\sim 3000$ BH-star systems in the Galactic halo formed through isolated evolution, of which $\sim 100$ are compatible with Gaia BH3 alike. Gaia BH3-like systems represent a common product of isolated binary evolution at low metallicity ($Z<0.01$), but given the steep density profile of the Galactic halo we do not expect more than one in the halo at the observed distance of Gaia BH3. Considering the estimated formation efficiency for both the isolated and dynamical formation channel, we conclude that they are almost equally likely to explain the origin of Gaia BH3. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Submitted to A&A, comments welcome

Classifying binary black holes from Population III stars with the Einstein Telescope: a machine-learning approach

Authors: Filippo Santoliquido, Ulyana Dupletsa, Jacopo Tissino, Marica Branchesi, Francesco Iacovelli, Giuliano Iorio, Michela Mapelli, Davide Gerosa, Jan Harms, Mario Pasquato

Abstract: Third-generation (3G) gravitational-wave (GW) detectors like the Einstein Telescope (ET) will observe binary black hole (BBH) mergers at redshifts up to $z\sim 100$. However, unequivocal determination of the origin of high-redshift sources will remain uncertain, due to the low signal-to-noise ratio (SNR) and poor estimate of their luminosity distance. This study proposes a machine learning approac… ▽ More Third-generation (3G) gravitational-wave (GW) detectors like the Einstein Telescope (ET) will observe binary black hole (BBH) mergers at redshifts up to $z\sim 100$. However, unequivocal determination of the origin of high-redshift sources will remain uncertain, due to the low signal-to-noise ratio (SNR) and poor estimate of their luminosity distance. This study proposes a machine learning approach to infer the origins of high-redshift BBHs, specifically differentiating those arising from Population III (Pop. III) stars - likely the first progenitors of stellar-born BBH mergers in the Universe - and those originated from Population I-II (Pop. I-II) stars. We have considered a wide range of state-of-the-art models encompassing current uncertainties on Pop. III BBH mergers. We then estimate parameter errors of detected sources with ET using the Fisher-information-matrix formalism, followed by classification using XGBoost, a machine learning algorithm based on decision trees. For a set of mock observed BBHs, we provide the probability that they belong to the Pop. III class while considering the parameter errors of each source. In our fiducial model, we accurately identify ~10% of detected BBHs originating from Pop. III stars with > 90% precision. Our study demonstrates how machine learning enables to achieve some pivotal aspects of ET science case by exploring the origin of individual high-redshift GW observations. We set the basis for further studies, which will integrate additional simulated populations and account for population modeling uncertainties △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Submitted to Astronomy & Astrophysics. 13 pages, 9 Figures and 5 tables. Comments are welcome

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

From gas to stars: MUSEings on the internal evolution of IC 1613

Authors: S. Taibi, G. Battaglia, M. M. Roth, S. Kamann, G. Iorio, C. Gallart, R. Leaman, E. D. Skillman, N. Kacharov, M. A. Beasley, P. E. Mancera Piña, G. van de Ven

Abstract: The kinematics and chemical composition of stellar populations of different ages provide crucial information about the evolution of a galaxy. We aim to provide such information for IC 1613, an isolated, gas-rich, star-forming dwarf galaxy in the Local Group. We present here the results of a new spectroscopic study performed with MUSE, an integral-field spectrograph on the Very Large Telescope. We… ▽ More The kinematics and chemical composition of stellar populations of different ages provide crucial information about the evolution of a galaxy. We aim to provide such information for IC 1613, an isolated, gas-rich, star-forming dwarf galaxy in the Local Group. We present here the results of a new spectroscopic study performed with MUSE, an integral-field spectrograph on the Very Large Telescope. We extracted from the data cubes more than 2000 sources from which we separated stellar objects for further spectroscopic analysis. The quality of the data set allowed us to obtain accurate classifications and line-of-sight velocities for about 800 stars. Our sample includes not only Red Giant Branch (RGB) and Main Sequence (MS) stars, but also a number of probable Be and C stars. We also obtained reliable metallicities for about 300 RGB stars. The kinematic analysis revealed for the first time the presence of stellar rotation with high significance. We found general agreement with the rotation velocity of the neutral gas component. Examining the kinematics of stars as a function of broad age ranges, we find that the velocity dispersion increases as a function of age, with the behaviour being very clear in the outermost pointings, while the rotation-to-velocity dispersion support decreases. On timescales shorter than a Gyr, the stellar kinematics still follow very closely that of the neutral gas, while the two components decouple on longer timescales. The chemical analysis of the RGB stars revealed average properties comparable to other Local Group dwarf galaxies. We also provide a new estimation of the inclination angle using only independent stellar tracers. Our work provides the largest spectroscopic sample of an isolated LG dwarf galaxy. The results obtained seem to support the scenario in which the stars of a dwarf galaxy are born from a less turbulent gas over time. △ Less

Submitted 7 June, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A. Changes from previous version include: a more structured text, updated figures and the addition of a discussion section (incl. a new figure). Previous results remain unchanged, but are now contextualised with respect to similar dwarf galaxies. A new estimate of the inclination angle has been added. References updated. 18 pages (10 figures, 5 tables) + appendices

Nebular emission from young stellar populations including binary stars

Authors: M. Lecroq, S. Charlot, A. Bressan, G. Bruzual, G. Costa, G. Iorio, M. Spera, M. Mapelli, Y. Chen, J. Chevallard, M. Dall'Amico

Abstract: We investigate the nebular emission produced by young stellar populations using the new GALSEVN model based on the combination of the SEVN population-synthesis code including binary-star processes and the GALAXEV code for the spectral evolution of stellar populations. Photoionization calculations performed with the CLOUDY code confirm that accounting for binary-star processes strongly influences t… ▽ More We investigate the nebular emission produced by young stellar populations using the new GALSEVN model based on the combination of the SEVN population-synthesis code including binary-star processes and the GALAXEV code for the spectral evolution of stellar populations. Photoionization calculations performed with the CLOUDY code confirm that accounting for binary-star processes strongly influences the predicted emission-line properties of young galaxies. In particular, we find that our model naturally reproduces the strong HeII/Hb ratios commonly observed at high Hb equivalent widths in metal-poor, actively star-forming galaxies, which have proven challenging to reproduce using previous models. Including bursty star formation histories broadens the agreement with observations, while the most extreme HeII equivalent widths can be reproduced by models dominated by massive stars. GALSEVN also enables us to compute, for the first time in a way physically consistent with stellar emission, the emission from accretion discs of X-ray binaries (XRBs) and radiative shocks driven by stellar winds and supernova explosions. We find that these contributions are unlikely to prominently affect the predicted HeII/Hb ratio, and that previous claims of a significant contribution by XRBs to the luminosities of high-ionization lines are based on models predicting improbably high ratios of X-ray luminosity to star formation rate, inconsistent with the observed average luminosity function of XRBs in nearby galaxies. The results presented here provide a solid basis for a more comprehensive investigation of the physical properties of observed galaxies with GALSEVN using Bayesian inference. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 26 pages, 18 figures. Accepted for publication in MNRAS

Impact of gas hardening on the population properties of hierarchical black hole mergers in AGN disks

Authors: M. Paola Vaccaro, Michela Mapelli, Carole Périgois, Dario Barone, M. Celeste Artale, Marco Dall'Amico, Giuliano Iorio, Stefano Torniamenti

Abstract: Hierarchical black hole (BH) mergers in active galactic nuclei (AGNs) are unique among formation channels of binary black holes (BBHs) because they are likely associated with electromagnetic counterparts and can efficiently lead to the mass growth of BHs. Here, we explore the impact of gas accretion and migration traps on the evolution of BBHs in AGNs. We have developed a new fast semi-analytic mo… ▽ More Hierarchical black hole (BH) mergers in active galactic nuclei (AGNs) are unique among formation channels of binary black holes (BBHs) because they are likely associated with electromagnetic counterparts and can efficiently lead to the mass growth of BHs. Here, we explore the impact of gas accretion and migration traps on the evolution of BBHs in AGNs. We have developed a new fast semi-analytic model, which allows us to explore the parameter space while capturing the main physical processes involved. We find that effective exchange of energy and angular momentum between the BBH and the surrounding gas (hereafter, gas hardening) during inspiral greatly enhances the efficiency of hierarchical mergers, leading to the formation of intermediate-mass BHs (up to 10.000 solar masses) and triggering spin alignment. Moreover, our models with efficient gas hardening show both an anti-correlation between BBH mass ratio and effective spin, and a correlation between primary BH mass and effective spin. In contrast, if gas hardening is inefficient, the hierarchical merger chain is already truncated after the first two or three generations. We compare the BBH population in AGNs with other dynamical channels as well as isolated binary evolution. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 21 pages, 15 figures, submitted to A&A, comments welcome

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

The Dragon-II simulations -- I. Evolution of single and binary compact objects in star clusters with up to 1 million stars

Authors: Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Mirek Giersz, Peter Berczik, Sara Rastello, Giuliano Iorio, Michela Mapelli, Massimiliano Gatto, Eva K. Grebel

Abstract: We present the first results of the \textsc{Dragon-II} simulations, a suite of 19 $N$-body simulations of star clusters with up to $10^6$ stars, with up to $33\%$ of them initially paired in binaries. In this work, we describe the main evolution of the clusters and their compact objects (COs). All \textsc{Dragon-II} clusters form in their centre a black hole (BH) subsystem with a density $10-100$… ▽ More We present the first results of the \textsc{Dragon-II} simulations, a suite of 19 $N$-body simulations of star clusters with up to $10^6$ stars, with up to $33\%$ of them initially paired in binaries. In this work, we describe the main evolution of the clusters and their compact objects (COs). All \textsc{Dragon-II} clusters form in their centre a black hole (BH) subsystem with a density $10-100$ times larger than the stellar density, with the cluster core containing $50-80\%$ of the whole BH population. In all models, the BH average mass steeply decreases as a consequence of BH burning, reaching values $\langle m_{\rm BH}\rangle < 15$ M$_\odot$ within $10-30$ relaxation times. Generally, our clusters retain only BHs lighter than $30$ M$_\odot$ over $30$ relaxation times. Looser clusters retain a higher binary fraction, because in such environments binaries are less likely disrupted by dynamical encounters. We find that BH-main sequence star binaries have properties similar to recently observed systems. Double CO binaries (DCOBs) ejected from the cluster exhibit larger mass ratios and heavier primary masses than ejected binaries hosting a single CO (SCOBs). Ejected SCOBs have BH masses $m_{\rm BH} = 3-20$ M$_\odot$, definitely lower than those in DCOBs ($m_{\rm BH} = 10-100$ M$_\odot$). △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 22 pages, 21 figures, 4 tables. Comments welcome. Submitted to MNRAS

Dynamical formation of $Gaia$ BH1 in a young star cluster

Authors: Sara Rastello, Giuliano Iorio, Michela Mapelli, Manuel Arca-Sedda, Ugo N. Di Carlo, Gastón J. Escobar, Stefano Torniamenti, Tomer Shenar

Abstract: $Gaia$ BH1, the first quiescent black hole (BH) detected from $Gaia$ data, poses a challenge to most binary evolution models: its current mass ratio is $\approx{0.1}$, and its orbital period seems to be too long for a post-common envelope system and too short for a non-interacting binary system. Here, we explore the hypothesis that $Gaia… ▽ More $Gaia$ BH1, the first quiescent black hole (BH) detected from $Gaia$ data, poses a challenge to most binary evolution models: its current mass ratio is $\approx{0.1}$, and its orbital period seems to be too long for a post-common envelope system and too short for a non-interacting binary system. Here, we explore the hypothesis that $Gaia$ BH1 formed through dynamical interactions in a young star cluster (YSC). We study the properties of BH-main sequence (MS) binaries formed in YSCs with initial mass $3\times{}10^2-3\times{}10^4$ M$_\odot$ at solar metallicity, by means of $3.5\times{}10^4$ direct $N$-body simulations coupled with binary population synthesis. For comparison, we also run a sample of isolated binary stars with the same binary population synthesis code used in the dynamical models. We find that BH-MS systems that form via dynamical exchanges populate the region corresponding to the main orbital properties of $Gaia$ BH1 (period, eccentricity, and masses). In contrast, none of our isolated binary systems matches the orbital period and MS mass of $Gaia$ BH1. Our best matching $Gaia$ BH1--like system forms via repeated dynamical exchanges and collisions involving the BH progenitor star, before it undergoes core collapse. YSCs are at least two orders of magnitude more efficient in forming $Gaia$ BH1--like systems than isolated binary evolution. △ Less

Submitted 25 July, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 10 pages, 7 figures, comments welcome (V2: Fig. 4 and references updated)

Binary neutron star populations in the Milky Way

Authors: Cecilia Sgalletta, Giuliano Iorio, Michela Mapelli, M. Celeste Artale, Lumen Boco, Debatri Chattopadhyay, Andrea Lapi, Andrea Possenti, Stefano Rinaldi, Mario Spera

Abstract: Galactic binary neutron stars (BNSs) are a unique laboratory to probe the evolution of BNSs and their progenitors. Here, we use a new version of the population synthesis code SEVN to evolve the population of Galactic BNSs, by modeling the spin up and down of pulsars self-consistently. We analyze the merger rate $\mathcal{R}_{\rm MW}$, orbital period $P_{\rm orb}$, eccentricity $e$, spin period… ▽ More Galactic binary neutron stars (BNSs) are a unique laboratory to probe the evolution of BNSs and their progenitors. Here, we use a new version of the population synthesis code SEVN to evolve the population of Galactic BNSs, by modeling the spin up and down of pulsars self-consistently. We analyze the merger rate $\mathcal{R}_{\rm MW}$, orbital period $P_{\rm orb}$, eccentricity $e$, spin period $P$, and spin period derivative $\dot{P}$ of the BNS population. Values of the common envelope parameter $α=1 - 3$ and an accurate model of the Milky Way star formation history best reproduce the BNS merger rate in our Galaxy ($\mathcal{R}_{\rm MW}\approx{}30$ Myr$^{-1}$). We apply radio-selection effects to our simulated BNSs and compare them to the observed population. Using a Dirichlet process Gaussian mixture method, we evaluate the four-dimensional likelihood in the $(P_{\rm orb}, e, P, \dot{P})$ space, by comparing our radio-selected simulated pulsars against Galactic BNSs. Our analysis favours an uniform initial distribution for both the magnetic field ($10^{10-13}$ G) and the spin period ($10-100$ ms). The implementation of radio selection effects is critical to match not only the spin period and period derivative, but also the orbital period and eccentricity of Galactic BNSs. According to our fiducial model, the Square Kilometre Array will detect $\sim 20$ new BNSs in the Milky Way. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 20 pages, 8 figures, 8 tables, comments welcome

Confronting fuzzy dark matter with the rotation curves of nearby dwarf irregular galaxies

Authors: Andrés Bañares-Hernández, Andrés Castillo, Jorge Martin Camalich, Giuliano Iorio

Abstract: We investigate phenomenologically the viability of fuzzy dark matter (FDM). We do this by confronting the predictions of the model, in particular, the formation of a solitonic core at the center of dark matter halos, with a homogeneous and robust sample of high-resolution rotation curves from the "LITTLE THINGS in 3D" catalog. This comprises a collection of isolated, dark matter dominated dwarf-ir… ▽ More We investigate phenomenologically the viability of fuzzy dark matter (FDM). We do this by confronting the predictions of the model, in particular, the formation of a solitonic core at the center of dark matter halos, with a homogeneous and robust sample of high-resolution rotation curves from the "LITTLE THINGS in 3D" catalog. This comprises a collection of isolated, dark matter dominated dwarf-irregular galaxies that provides an optimal benchmark for cosmological studies. We use a statistical framework based on Markov chain Monte Carlo techniques that allows us to extract relevant parameters such as the axion mass, the mass of the solitonic core, the mass of the dark matter halo and its concentration parameter with a rather loose set of priors except for the implementation of a core-halo relation that is predicted by simulations. The results of the fits are used to perform various diagnostics on the predictions of the model. FDM provides an excellent fit to the rotation curves of the "LITTLE THINGS in 3D" catalog, with axion masses determined from different galaxies clustering around $m_a\approx2\times10^{-23}$ eV. However, we find two major problems in our analysis. First, the data follow scaling relations of the properties of the core which are not consistent with the predictions of the soliton. This problem is particularly acute in the core radius - mass relation with a tension that, at face value, has a significance $\gtrsim5σ$. The second problem is related to the strong suppression of the linear power spectrum that is predicted by FDM for the axion mass preferred by the data. This can be constrained very conservatively by the galaxy counts in our sample, which leads to a tension exceeding again $5σ$. We estimate the effects of baryons in our analysis and discuss whether they could alleviate the tensions of the model with observations. △ Less

Submitted 14 March, 2024; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 22 pages, 9 figures, 3 tables

Journal ref: A&A 676, A63 (2023)

Binary black hole mergers from Population III stars: uncertainties from star formation and binary star properties

Authors: Filippo Santoliquido, Michela Mapelli, Giuliano Iorio, Guglielmo Costa, Simon C. O. Glover, Tilman Hartwig, Ralf S. Klessen, Lorenzo Merli

Abstract: Population III (Pop. III) binary stars likely produced the first stellar-born binary black hole (BBH) mergers in the Universe. Here, we quantify the main sources of uncertainty for the merger rate density evolution and mass spectrum of Pop. III BBHs by considering four different formation histories and 11 models of the initial orbital properties of Pop. III binary stars. The uncertainty on the orb… ▽ More Population III (Pop. III) binary stars likely produced the first stellar-born binary black hole (BBH) mergers in the Universe. Here, we quantify the main sources of uncertainty for the merger rate density evolution and mass spectrum of Pop. III BBHs by considering four different formation histories and 11 models of the initial orbital properties of Pop. III binary stars. The uncertainty on the orbital properties affects the BBH merger rate density by up to two orders of magnitude, models with shorter orbital periods leading to higher BBH merger rates. The uncertainty on the star formation history has a substantial impact on both the shape and the normalisation of the BBH merger rate density: the peak of the merger rate density shifts from $z\sim{8}$ up to $z\sim{16}$ depending on the assumed star formation rate, while the maximum BBH merger rate density for our fiducial binary population model spans from $\sim{2}$ to $\sim{30}$ Gpc$^{-3}$ yr$^{-1}$. The typical BBH masses are not affected by the star formation rate model and only mildly influenced by the binary population parameters. The primary black holes born from Pop. III stars tend to be rather massive ($30-40$ M$_\odot$) with respect to those born from metal-rich stars ($8-10$ M$_\odot$). We estimate that the Einstein Telescope will detect $10-10^4$ Pop. III BBH mergers per year, depending on the star formation history and binary star properties. △ Less

Submitted 29 January, 2024; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 19 pages, 20 figures, 1 table. We fixed a bug. Results are not significantly affected. See also correction published in MNRAS

Massive binary black holes from Population II and III stars

Authors: Guglielmo Costa, Michela Mapelli, Giuliano Iorio, Filippo Santoliquido, Gastón J. Escobar, Ralf S. Klessen, Alessandro Bressan

Abstract: Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity $Z=10^{-11}$) with o… ▽ More Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity $Z=10^{-11}$) with our population-synthesis code SEVN, and compare them with Population II binary stars ($Z=10^{-4}$). In our models, the lower edge of the pair-instability mass gap corresponds to a BH mass of $\approx{86}$ ($\approx{91}$) M$_\odot$ for single Population III (II) stars. Overall, we find only mild differences between the properties of binary BHs (BBHs) born from Population III and II stars, especially if we adopt the same initial mass function and initial orbital properties. Most BBH mergers born from Population III and II stars have primary BH mass below the pair-instability gap, and the maximum secondary BH mass is $ < 50$ M$_\odot$. Only up to $\approx{3.3}$% ($\approx{0.09}$%) BBH mergers from Population III (II) progenitors have primary mass above the gap. Unlike metal-rich binary stars, the main formation channel of BBH mergers from Population III and II stars involves only stable mass transfer episodes in our fiducial model. △ Less

Submitted 24 August, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 15 pages, 17 figures, Comments are welcome. Published in MNRAS

Evidence for supernova feedback sustaining gas turbulence in nearby star-forming galaxies

Authors: C. Bacchini, F. Fraternali, G. Pezzulli, G. Iorio, A. Marasco, C. Nipoti

Abstract: HI and CO observations indicate that the cold gas in galaxies is very turbulent. However, the turbulent energy is expected to be quickly dissipated, implying that some energy source is needed to explain the observations. The nature of such turbulence was long unclear, as even the main candidate, supernova (SN) feedback, seemed insufficient. Other mechanisms have been proposed, but without reaching… ▽ More HI and CO observations indicate that the cold gas in galaxies is very turbulent. However, the turbulent energy is expected to be quickly dissipated, implying that some energy source is needed to explain the observations. The nature of such turbulence was long unclear, as even the main candidate, supernova (SN) feedback, seemed insufficient. Other mechanisms have been proposed, but without reaching a general consensus. The key novelty of our work is considering that the gas disc thickness and flaring increase the dissipation timescale of turbulence, thus reducing the energy injection rate required to sustain it. In excellent agreement with the theoretical expectations, we found that the fraction of the SN energy (a.k.a. SN coupling efficiency) needed to maintain the cold gas turbulence is $\sim 1$%, solving a long-standing conundrum. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 3 pages, 1 figure, Proceedings IAU Symposium 373, accepted

The volumetric star formation law in nearby galaxies

Authors: C. Bacchini, F. Fraternali, G. Pezzulli, G. Iorio, A. Marasco, C. Nipoti

Abstract: Star formation laws are empirical relations between the cold gas (HI+H$_2$) content of a galaxy and its star formation rate (SFR), being crucial for any model of galaxy formation and evolution. A well known example of such laws is the Schmidt-Kennicutt law, which is based on the projected surface densities. However, it has been long unclear whether a more fundamental relation exists between the in… ▽ More Star formation laws are empirical relations between the cold gas (HI+H$_2$) content of a galaxy and its star formation rate (SFR), being crucial for any model of galaxy formation and evolution. A well known example of such laws is the Schmidt-Kennicutt law, which is based on the projected surface densities. However, it has been long unclear whether a more fundamental relation exists between the intrinsic volume densities. By assuming the vertical hydrostatic equilibrium, we infer radial profiles for the thickness of gaseous discs in a sample of 23 local galaxies, and use these measurements to convert the observed surface densities of the gas and the SFR into the de-projected volume densities. We find a tight correlation linking these quantities, that we call the volumetric star formation law. This relation and its properties have crucial implications for our understanding of the physics of star formation. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 3 pages, 1 figure, Proceedings IAU Symposium 373, accepted

Compact object mergers: exploring uncertainties from stellar and binary evolution with SEVN

Authors: Giuliano Iorio, Michela Mapelli, Guglielmo Costa, Mario Spera, Gastón J. Escobar, Cecilia Sgalletta, Alessandro A. Trani, Erika Korb, Filippo Santoliquido, Marco Dall'Amico, Nicola Gaspari, Alessandro Bressan

Abstract: Population-synthesis codes are an unique tool to explore the parameter space of massive binary star evolution and binary compact object (BCO) formation. Most population-synthesis codes are based on the same stellar evolution model, limiting our ability to explore the main uncertainties. Here, we present the new version of the code SEVN, which overcomes this issue by interpolating the main stellar… ▽ More Population-synthesis codes are an unique tool to explore the parameter space of massive binary star evolution and binary compact object (BCO) formation. Most population-synthesis codes are based on the same stellar evolution model, limiting our ability to explore the main uncertainties. Here, we present the new version of the code SEVN, which overcomes this issue by interpolating the main stellar properties from a set of pre-computed evolutionary tracks. We describe the new interpolation and adaptive time-step algorithms of SEVN, and the main upgrades on single and binary evolution. With SEVN, we evolved $1.2\times10^9$ binaries in the metallicity range $0.0001\leq Z \leq 0.03$, exploring a number of models for electron-capture, core-collapse and pair-instability supernovae, different assumptions for common envelope, stability of mass transfer, quasi-homogeneous evolution and stellar tides. We find that stellar evolution has a dramatic impact on the formation of single and binary compact objects. Just by slightly changing the overshooting parameter ($λ_{\rm ov}=0.4,0.5$) and the pair-instability model, the maximum mass of a black hole can vary from $\approx{60}$ to $\approx{100}\ \mathrm{M}_\odot$. Furthermore, the formation channels of BCOs and the merger efficiency we obtain with SEVN show significant differences with respect to the results of other population-synthesis codes, even when the same binary-evolution parameters are used. For example, the main traditional formation channel of BCOs is strongly suppressed in our models: at high metallicity ($Z\gtrsim{0.01}$) only $<20$% of the merging binary black holes and binary neutron stars form via this channel, while other authors found fractions $>70$%. The local BCO merger rate density of our fiducial models is consistent with the most recent estimates by the LIGO--Virgo--KAGRA collaboration. △ Less

Submitted 29 June, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. The SEVN code is available at https://gitlab.com/sevncodes/sevn.git. All the data underlying this article are available in Zenodo at the link https://doi.org/10.5281/zenodo.7794546. All the Jupyter notebooks used to produce the plots in the paper are available in the gitlab repository https://gitlab.com/iogiul/iorio22_plot.git

RR Lyrae From Binary Evolution: Abundant, Young and Metal-Rich

Authors: Alexey Bobrick, Giuliano Iorio, Vasily Belokurov, Joris Vos, Maja Vuckovic, Nicola Giacobbo

Abstract: RR Lyrae are a well-known class of pulsating horizontal branch stars widely used as tracers of old, metal-poor stellar populations. However, mounting observational evidence shows that a significant fraction of these stars may be young and metal-rich. Here, through detailed binary stellar evolution modelling, we show that all such metal-rich RR Lyrae can be naturally produced through binary interac… ▽ More RR Lyrae are a well-known class of pulsating horizontal branch stars widely used as tracers of old, metal-poor stellar populations. However, mounting observational evidence shows that a significant fraction of these stars may be young and metal-rich. Here, through detailed binary stellar evolution modelling, we show that all such metal-rich RR Lyrae can be naturally produced through binary interactions. Binary companions of these RR Lyrae stars formed through binary interactions partly strip their progenitor's envelopes during a preceding red giant phase. As a result, stripped horizontal branch stars become bluer than their isolated stellar evolution counterparts and thus end up in the instability strip. In contrast, in the single evolution scenario, the stars can attain such colours only at large age and low metallicity. While binary-made RR Lyrae can possess any ages and metallicities, their Galactic population is relatively young (1 to 9 Gyr) and dominated by the Thin Disc and the Bulge. We show that Galactic RR Lyrae from binary evolution are produced at rates compatible with the observed metal-rich population and have consistent G-band magnitudes, Galactic kinematics and pulsation properties. Furthermore, these systems dominate the RR Lyrae population in the Solar Neighbourhood. We predict that all metal-rich RR Lyrae have an A, F, G or K-type companion with a long orbital period (P > 1000 d). Observationally characterising the orbital periods and masses of such stellar companions will provide valuable new constraints on mass and angular momentum-loss efficiency for Sun-like accretors and the nature of RR Lyrae populations. △ Less

Submitted 13 January, 2024; v1 submitted 8 August, 2022; originally announced August 2022.

Comments: Accepted version, comments welcome

Revisiting Common Envelope Evolution -- A New Semi-Analytic Model for N-body and Population Synthesis Codes

Authors: Alessandro Alberto Trani, Steven Rieder, Ataru Tanikawa, Giuliano Iorio, Riccardo Martini, Georgii Karelin, Hila Glanz, Simon Portegies Zwart

Abstract: We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory.… ▽ More We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory. Our formalism can model the eccentricity during the common envelope inspiral, and it gives results consistent with smoothed particles hydrodynamical simulations. We apply our formalism to common envelope events from binary population synthesis models and find that the final eccentricity distribution resembles the observed distribution of post-common-envelope binaries. Our model can be used for time-resolved common-envelope evolution in population synthesis calculations or as part of binary interactions in direct N-body simulations of star clusters. △ Less

Submitted 23 August, 2022; v1 submitted 26 May, 2022; originally announced May 2022.

Comments: 14 Pages, 8 Figures, accepted for publication in PRD

The impact of gas disc flaring on rotation curve decomposition and revisiting baryonic and dark-matter relations for nearby galaxies

Authors: Pavel E. Mancera Piña, Filippo Fraternali, Tom Oosterloo, Elizabeth A. K. Adams, Enrico di Teodoro, Cecilia Bacchini, Giuliano Iorio

Abstract: Gas discs of late-type galaxies are flared, with scale heights increasing with the distance from the galaxy centres and often reaching kpc scales. We study the effects of gas disc flaring on the recovered dark matter halo parameters from rotation curve decomposition. For this, we carefully select a sample of 32 dwarf and spiral galaxies with high-quality neutral gas, molecular gas, and stellar mas… ▽ More Gas discs of late-type galaxies are flared, with scale heights increasing with the distance from the galaxy centres and often reaching kpc scales. We study the effects of gas disc flaring on the recovered dark matter halo parameters from rotation curve decomposition. For this, we carefully select a sample of 32 dwarf and spiral galaxies with high-quality neutral gas, molecular gas, and stellar mass profiles, robust H\,{\sc i} rotation curves obtained via 3D kinematic modelling, and reliable bulge-disc decomposition. By assuming vertical hydrostatic equilibrium, we derive the scale heights of the atomic and molecular gas discs and fit dark matter haloes to the rotation curves self-consistently. We find that the effect of the gas flaring in the rotation curve decomposition can play an important role only for the smallest, gas-dominated dwarfs, while for most of the galaxies the effect is minor and can be ignored. We revisit the stellar- and baryon-to-halo mass relations ($M_\ast-M_{200}$ and $M_{\rm bar}-M_{200}$). Both relations increase smoothly up to $M_{200} \approx 10^{12}~\rm{ M_\odot}$, with galaxies at this end having high $M_\ast/M_{200}$ and $M_{\rm bar}/M_{200}$ ratios approaching the cosmological baryon fraction. At higher $M_{200}$ the relations show a larger scatter. Most haloes of our galaxy sample closely follow the concentration-mass ($c_{200}-M_{\rm 200}$) relation resulting from N-body cosmological simulations. Interestingly, the galaxies deviating above and below the relation have the highest and lowest stellar and baryon factions, respectively, which suggests that the departures from the $c_{200}-M_{\rm 200}$ law are regulated by adiabatic contraction and an increasing importance of feedback. △ Less

Submitted 25 May, 2022; originally announced May 2022.

Comments: accepted for publication in MNRAS. Supplementary material will be available from MNRAS and meantime at this link: https://www.dropbox.com/sh/768phag6ttgym0j/AABZwerBgFofkkZ5npHW17YBa

The Smallest Scale of Hierarchy Survey (SSH). II. Extended star formation and bar-like features in the dwarf galaxy NGC 3741: recent merger or ongoing gas accretion?

Authors: F. Annibali, C. Bacchini, G. Iorio, M. Bellazzini, R. Pascale, G. Beccari, M. Cignoni, L. Ciotti, C. Nipoti, E. Sacchi, M. Tosi, F. Cusano, S. Bisogni, A. Gargiulo, D. Paris

Abstract: Using Large Binocular Telescope deep imaging data from the Smallest Scale of Hierarchy Survey (SSH) and archival Hubble Space Telescope data, we reveal the presence of two elongated stellar features contiguous to a bar-like stellar structure in the inner regions of the dwarf irregular galaxy NGC 3741. These structures are dominated by stars younger than a few hundred Myr and collectively are about… ▽ More Using Large Binocular Telescope deep imaging data from the Smallest Scale of Hierarchy Survey (SSH) and archival Hubble Space Telescope data, we reveal the presence of two elongated stellar features contiguous to a bar-like stellar structure in the inner regions of the dwarf irregular galaxy NGC 3741. These structures are dominated by stars younger than a few hundred Myr and collectively are about twice as extended as the old stellar component. These properties are very unusual for dwarf galaxies in the nearby Universe and difficult to explain by hydro-dynamical simulations. From the analysis of archival 21-cm observations, we find that the young stellar "bar" coincides with an HI high-density region proposed by previous studies to be a purely gaseous bar; we furthermore confirm radial motions of a few km/s, compatible with an inflow/outflow, and derive a steeply-rising rotation curve and high HI surface density at the center, indicating a very concentrated mass distribution. We propose that the peculiar properties of the stellar and gaseous components of NGC 3741 may be explained by a recent merger or ongoing gas accretion from the intergalactic medium, which caused gas inflows towards the galaxy center and triggered star formation a few hundred Myr ago. This event may explain the young and extended stellar features, the bar-like structure, the very extended HI disc and the central HI spiral arms. The high central HI density and the steeply rising rotation curve suggest that NGC 3741 may be the progenitor or the descendant of a starburst dwarf. △ Less

Submitted 25 February, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS. Corrected typo in the list of authors

Hierarchical generative models for star clusters from hydro-dynamical simulations

Authors: Stefano Torniamenti, Mario Pasquato, Pierfrancesco Di Cintio, Alessandro Ballone, Giuliano Iorio, M. Celeste Artale, Michela Mapelli

Abstract: Star formation in molecular clouds is clumpy, hierarchically subclustered. Fractal structure also emerges in hydro-dynamical simulations of star-forming clouds. Simulating the formation of realistic star clusters with hydro-dynamical simulations is a computational challenge, considering that only the statistically averaged results of large batches of simulations are reliable, due to the chaotic na… ▽ More Star formation in molecular clouds is clumpy, hierarchically subclustered. Fractal structure also emerges in hydro-dynamical simulations of star-forming clouds. Simulating the formation of realistic star clusters with hydro-dynamical simulations is a computational challenge, considering that only the statistically averaged results of large batches of simulations are reliable, due to the chaotic nature of the gravitational N-body problem. While large sets of initial conditions for N-body runs can be produced by hydro-dynamical simulations of star formation, this is prohibitively expensive in terms of computational time. Here we address this issue by introducing a new technique for generating many sets of new initial conditions from a given set of star masses, positions and velocities from a hydro-dynamical simulation. We use hierarchical clustering in phase space to learn a tree representation of the spatial and kinematic relations between stars. This constitutes the basis for the random generation of new sets of stars which share the same clustering structure of the original ones but have individually different masses, positions, and velocities. We apply this method to the output of a number of hydro-dynamical star-formation simulations, comparing the generated initial conditions to the original ones through a series of quantitative tests, including comparing mass and velocity distributions and fractal dimension. Finally, we evolve both the original and the generated star clusters using a direct N-body code, obtaining a qualitatively similar evolution. △ Less

Submitted 24 December, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: Updated version of the manuscript: "Seeing the forest for the trees: hierarchical generative models for star clusters from hydro-dynamical simulations". 15 pages, 15 figures, 2 tables. Comments welcome

Dynamics of binary black holes in low-mass young star clusters

Authors: Sara Rastello, Michela Mapelli, Ugo N. di Carlo, Giuliano Iorio, Alessandro Ballone, Nicola Giacobbo, Filippo Santoliquido, Stefano Torniamenti

Abstract: Young star clusters are dynamically active stellar systems and are a common birthplace for massive stars. Low-mass star clusters ($\sim{}300-10^3$ M$_\odot$) are more numerous than massive systems and are characterized by a two-body relaxation time scale of a few Myr: the most massive stars sink to the cluster core and dynamically interact with each other even before they give birth to compact obj… ▽ More Young star clusters are dynamically active stellar systems and are a common birthplace for massive stars. Low-mass star clusters ($\sim{}300-10^3$ M$_\odot$) are more numerous than massive systems and are characterized by a two-body relaxation time scale of a few Myr: the most massive stars sink to the cluster core and dynamically interact with each other even before they give birth to compact objects. Here, we explore the properties of black holes (BHs) and binary black holes (BBHs) formed in low-mass young star clusters, by means of a suite of $10^5$ direct $N$-body simulations with a high original binary fraction (100 % for stars with mass $>5$ M$_\odot$). Most BHs are ejected in the first $\sim{}20$ Myr by dynamical interactions. Dynamical exchanges are the main formation channel of BBHs, accounting for $\sim{}40-80$ % of all the systems. Most BBH mergers in low-mass young star clusters involve primary BHs with mass $<40$ M$_\odot$ and low mass ratios are extremely more common than in the field. Comparing our data with those of more massive star clusters ($10^3-3\times{}10^4$ M$_\odot$), we find a strong dependence of the percentage of exchanged BBHs on the mass of the host star cluster. In contrast, our results show just a mild correlation between the mass of the host star cluster and the efficiency of BBH mergers. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 14 pages, 13 figures, comments welcome

Intermediate mass black holes from stellar mergers in young star clusters

Authors: Ugo N. Di Carlo, Michela Mapelli, Mario Pasquato, Sara Rastello, Alessandro Ballone, Marco Dall'Amico, Nicola Giacobbo, Giuliano Iorio, Mario Spera, Stefano Torniamenti, Francesco Haardt

Abstract: Intermediate mass black holes (IMBHs) in the mass range $10^2-10^5\,\mathrm{M_{\odot}}$ bridge the gap between stellar black holes (BHs) and supermassive BHs. Here, we investigate the possibility that IMBHs form in young star clusters via runaway collisions and BH mergers. We analyze $10^4$ simulations of dense young star clusters, featuring up-to-date stellar wind models and prescriptions for cor… ▽ More Intermediate mass black holes (IMBHs) in the mass range $10^2-10^5\,\mathrm{M_{\odot}}$ bridge the gap between stellar black holes (BHs) and supermassive BHs. Here, we investigate the possibility that IMBHs form in young star clusters via runaway collisions and BH mergers. We analyze $10^4$ simulations of dense young star clusters, featuring up-to-date stellar wind models and prescriptions for core collapse and (pulsational) pair instability. In our simulations, only 9 IMBHs out of 218 form via binary BH mergers, with a mass $\sim{}100-140$ M$_\odot$. This channel is strongly suppressed by the low escape velocity of our star clusters. In contrast, IMBHs with masses up to $\sim{}438$ M$_{\odot}$ efficiently form via runaway stellar collisions, especially at low metallicity. Up to $\sim{}0.2$~% of all the simulated BHs are IMBHs, depending on progenitor's metallicity. The runaway formation channel is strongly suppressed in metal-rich ($Z=0.02$) star clusters, because of stellar winds. IMBHs are extremely efficient in pairing with other BHs: $\sim{}70$% of them are members of a binary BH at the end of the simulations. However, we do not find any IMBH-BH merger. More massive star clusters are more efficient in forming IMBHs: $\sim{}8$% ($\sim{}1$%) of the simulated clusters with initial mass $10^4-3\times{}10^4$ M$_\odot$ ($10^3-5\times{}10^3$ M$_\odot$) host at least one IMBH. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 12 pages, 16 figures, 2 tables. Comments welcome

Hierarchical black hole mergers in young, globular and nuclear star clusters: the effect of metallicity, spin and cluster properties

Authors: Michela Mapelli, Marco Dall'Amico, Yann Bouffanais, Nicola Giacobbo, Manuel Arca Sedda, M. Celeste Artale, Alessandro Ballone, Ugo N. Di Carlo, Giuliano Iorio, Filippo Santoliquido, Stefano Torniamenti

Abstract: We explore hierarchical black hole (BH) mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs), accounting for both original and dynamically assembled binary BHs (BBHs). We find that the median mass of both first- and nth-generation dynamical mergers is larger in GCs and YSCs with respect to NSCs, because the lighter BHs are ejected by supernova kicks from… ▽ More We explore hierarchical black hole (BH) mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs), accounting for both original and dynamically assembled binary BHs (BBHs). We find that the median mass of both first- and nth-generation dynamical mergers is larger in GCs and YSCs with respect to NSCs, because the lighter BHs are ejected by supernova kicks from the lower-mass clusters. Also, first- and nth-generation BH masses are strongly affected by the metallicity of the progenitor stars: the median mass of the primary BH of a nth-generation merger is $\sim{}24-38$ M$_\odot$ ($\sim{}9-15$ M$_\odot$) in metal-poor (metal-rich) NSCs. The maximum BH mass mainly depends on the escape velocity: BHs with mass up to several thousand M$_\odot$ form in NSCs, while YSCs and GCs host BHs with mass up to several hundred M$_\odot$. Furthermore, we calculate the fraction of mergers with at least one component in the pair-instability mass gap ($f_{\rm PI}$) and in the intermediate-mass BH regime ($f_{\rm IMBH}$). In the fiducial model for dynamical BBHs with metallicity $Z=0.002$, we find $f_{\rm PI}\approx{}0.05$, $0.02$ and $0.007$ ($f_{\rm IMBH}\approx{}0.01$, $0.002$ and $0.001$) in NSCs, GCs and YSCs, respectively. Both $f_{\rm PI}$ and $f_{\rm IMBH}$ drop by at least one order of magnitude at solar metallicity. Finally, we investigate the formation of GW190521 by assuming that it is either a nearly equal-mass BBH or an intermediate-mass ratio inspiral. △ Less

Submitted 6 May, 2021; v1 submitted 8 March, 2021; originally announced March 2021.

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

New insights on binary black hole formation channels after GWTC-2: young star clusters versus isolated binaries

Authors: Yann Bouffanais, Michela Mapelli, Filippo Santoliquido, Nicola Giacobbo, Ugo N. Di Carlo, Sara Rastello, M. Celeste Artale, Giuliano Iorio

Abstract: With the recent release of the second gravitational-wave transient catalogue (GWTC-2), which introduced dozens of new detections, we are at a turning point of gravitational wave astronomy, as we are now able to directly infer constraints on the astrophysical population of compact objects. Here, we tackle the burning issue of understanding the origin of binary black hole (BBH) mergers. To this effe… ▽ More With the recent release of the second gravitational-wave transient catalogue (GWTC-2), which introduced dozens of new detections, we are at a turning point of gravitational wave astronomy, as we are now able to directly infer constraints on the astrophysical population of compact objects. Here, we tackle the burning issue of understanding the origin of binary black hole (BBH) mergers. To this effect, we make use of state-of-the-art population synthesis and N-body simulations, to represent two distinct formation channels: BBHs formed in the field (isolated channel) and in young star clusters (dynamical channel). We then use a Bayesian hierarchical approach to infer the distribution of the mixing fraction $f$, with $f=0$ ($f=1$) in the pure dynamical (isolated) channel. %that controls the proportion of isolated and dynamical BBHs. We explore the effects of additional hyper-parameters of the model, such as the spread in metallicity $σ_{\text{Z}}$ and the parameter $σ_{\text{sp}}$, describing the distribution of spin magnitudes. We find that the dynamical model is slightly favoured with a median value of $f=0.26$, when $σ_{\text{sp}}=0.1$ and $σ_{\text{Z}}=0.4$. Models with higher spin magnitudes tend to strongly favour dynamically formed BBHs ($f\le{}0.1$ if $σ_{\text{sp}}=0.3$). Furthermore, we show that hyper-parameters controlling the rates of the model, such as $σ_{\rm Z}$, have a large impact on the inference of the mixing fraction, which rises from $0.18$ to $0.43$ when we increase $σ_{\text{Z}}$ from 0.2 to 0.6, for a fixed value of $σ_{\text{sp}}=0.1$. Finally, our current set of observations is better described by a combination of both formation channels, as a pure dynamical scenario is excluded at the $99\%$ credible interval, except when the spin magnitude is high. △ Less

Submitted 5 December, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 13 pages, 10 figures, 2 tables, published in MNRAS

Effective N-body models of composite collisionless stellar systems

Authors: Carlo Nipoti, Giacomo Cherchi, Giuliano Iorio, Francesco Calura

Abstract: Gas-poor galaxies can be modelled as composite collisionless stellar systems, with a dark matter halo and one or more stellar components, representing different stellar populations. The dynamical evolution of such composite systems is often studied with numerical N-body simulations, whose initial conditions typically require realizations with particles of stationary galaxy models. We present a nov… ▽ More Gas-poor galaxies can be modelled as composite collisionless stellar systems, with a dark matter halo and one or more stellar components, representing different stellar populations. The dynamical evolution of such composite systems is often studied with numerical N-body simulations, whose initial conditions typically require realizations with particles of stationary galaxy models. We present a novel method to conceive these N-body realizations, which allows one to exploit at best a collisionless N-body simulation that follows their evolution. The method is based on the use of an effective N-body model of a composite system, which is in fact realized as a one-component system of particles that is interpreted a posteriori as a multi-component system, by assigning in post-processing fractions of each particle's mass to different components. Examples of astrophysical applications are N-body simulations that aim to reproduce the observed properties of interacting galaxies, satellite galaxies and stellar streams. As a case study we apply our method to an N-body simulation of tidal stripping of a two-component (dark matter and stars) satellite dwarf galaxy orbiting in the gravitational potential of the Milky Way. △ Less

Submitted 11 March, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 10 pages, 9 figures, accepted for publication in MNRAS. Figures 3, 7 and 9 illustrate the method applied to tidal streams

MSC Class: 85A05 ACM Class: J.2

From hydrodynamics to N-body simulations of star clusters: mergers and rotation

Authors: Alessandro Ballone, Stefano Torniamenti, Michela Mapelli, Ugo N. Di Carlo, Mario Spera, Sara Rastello, Nicola Gaspari, Giuliano Iorio

Abstract: We present a new method to obtain more realistic initial conditions for N-body simulations of young star clusters. We start from the outputs of hydrodynamical simulations of molecular cloud collapse, in which star formation is modelled with sink particles. In our approach, we instantaneously remove gas from these hydrodynamical simulation outputs to mock the end of the gas-embedded phase, induced… ▽ More We present a new method to obtain more realistic initial conditions for N-body simulations of young star clusters. We start from the outputs of hydrodynamical simulations of molecular cloud collapse, in which star formation is modelled with sink particles. In our approach, we instantaneously remove gas from these hydrodynamical simulation outputs to mock the end of the gas-embedded phase, induced by stellar feedback. We then enforce a realistic initial mass function by splitting or joining the sink particles based on their mass and position. Such initial conditions contain more consistent information on the spatial distribution and the kinematical and dynamical states of young star clusters, which are fundamental to properly study these systems. For example, by applying our method to a set of previously run hydrodynamical simulations, we found that the early evolution of young star clusters is affected by gas removal and by the early dry merging of sub-structures. This early evolution can either quickly erase the rotation acquired by our (sub-)clusters in their embedded phase or "fuel" it by feeding of angular momentum by sub-structure mergers, before two-body relaxation acts on longer timescales. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 15 pages, 15 figures, accepted on MNRAS

Constraining accretion efficiency in massive binary stars with LIGO-Virgo black holes

Authors: Yann Bouffanais, Michela Mapelli, Filippo Santoliquido, Nicola Giacobbo, Giuliano Iorio, Guglielmo Costa

Abstract: The growing sample of LIGO--Virgo black holes (BHs) opens new perspectives for the study of massive binary evolution. Here, we study the impact of mass accretion efficiency and common envelope on the properties of binary BH (BBH) mergers, by means of population synthesis simulations. We model mass accretion efficiency with the parameter $f_{\rm MT}\in[0.05,1]$, which represents the fraction of mas… ▽ More The growing sample of LIGO--Virgo black holes (BHs) opens new perspectives for the study of massive binary evolution. Here, we study the impact of mass accretion efficiency and common envelope on the properties of binary BH (BBH) mergers, by means of population synthesis simulations. We model mass accretion efficiency with the parameter $f_{\rm MT}\in[0.05,1]$, which represents the fraction of mass lost from the donor which is effectively accreted by the companion. Lower values of $f_{\rm MT}$ result in lower BBH merger rate densities and produce mass spectra skewed towards lower BH masses. Our hierarchical Bayesian analysis, applied to BBH mergers in the first and second gravitational wave transient catalogue, yields zero support for values of $f_{\rm MT}\lesssim{}0.6$, with a lower boundary of the 99\% credible intervals equal to $f_{\rm MT}= 0.59$. This result holds for all the values of the common-envelope efficiency parameter we considered in this study $α_{\rm CE} \in [1,10]$. This confirms that gravitational-wave data can be used to put constraints on several uncertain binary evolution processes. △ Less

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

Comments: 9 pages, 5 figures, published in MNRAS

Formation of GW190521 from stellar evolution: the impact of the hydrogen-rich envelope, dredge-up and $^{12}$C($α$, $γ$)$^{16}$O rate on the pair-instability black hole mass gap

Authors: Guglielmo Costa, Alessandro Bressan, Michela Mapelli, Paola Marigo, Giuliano Iorio, Mario Spera

Abstract: Pair-instability (PI) is expected to open a gap in the mass spectrum of black holes (BHs) between $\approx{}40-65$ M$_\odot$ and $\approx{}120$ M$_\odot$. The existence of the mass gap is currently being challenged by the detection of GW190521, with a primary component mass of $85^{+21}_{-14}$ M$_{\odot}$. Here, we investigate the main uncertainties on the PI mass gap: the $^{12}$C($α$, $γ$)… ▽ More Pair-instability (PI) is expected to open a gap in the mass spectrum of black holes (BHs) between $\approx{}40-65$ M$_\odot$ and $\approx{}120$ M$_\odot$. The existence of the mass gap is currently being challenged by the detection of GW190521, with a primary component mass of $85^{+21}_{-14}$ M$_{\odot}$. Here, we investigate the main uncertainties on the PI mass gap: the $^{12}$C($α$, $γ$)$^{16}$O reaction rate and the H-rich envelope collapse. With the standard $^{12}$C($α$, $γ$)$^{16}$O rate, the lower edge of the mass gap can be 70 M$_\odot$ if we allow for the collapse of the residual H-rich envelope at metallicity $Z\leq{}0.0003$. Adopting the uncertainties given by the STARLIB database, for models computed with the $^{12}$C($α$, $γ$)$^{16}$O rate $-1\, σ$, we find that the PI mass gap ranges between $\approx{}80$ M$_\odot$ and $\approx{}150$ M$_\odot$. Stars with $M_{\rm ZAMS}>110$ M$_\odot$ may experience a deep dredge-up episode during the core helium-burning phase, that extracts matter from the core enriching the envelope. As a consequence of the He-core mass reduction, a star with $M_{\rm ZAMS} =160$ M$_\odot$ may avoid the PI and produce a BH of 150 M$_\odot$. In the $-2\,{}σ$ case, the PI mass gap ranges from 92 M$_\odot$ to 110 M$_\odot$. Finally, in models computed with $^{12}$C($α$, $γ$)$^{16}$O $-3\,{}σ$, the mass gap is completely removed by the dredge-up effect. The onset of this dredge-up is particularly sensitive to the assumed model for convection and mixing. The combined effect of H-rich envelope collapse and low $^{12}$C($α$, $γ$)$^{16}$O rate can lead to the formation of BHs with masses consistent with the primary component of GW190521. △ Less

Submitted 22 December, 2020; v1 submitted 5 October, 2020; originally announced October 2020.

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

Chemo-kinematics of the $Gaia$ RR Lyrae: the halo and the disc

Authors: G. Iorio, V. Belokurov

Abstract: We present the results of a multi-component kinematic model of a large sample of RR Lyrae detected by $Gaia$. By imposing a four-fold symmetry and employing $Gaia$ proper motions, we are able to infer the behaviour of the velocity ellipsoid between $\approx3$ and $\approx30$ kpc from the centre of the Galaxy. We detect the presence of two distinct components: a dominant non-rotating halo-like popu… ▽ More We present the results of a multi-component kinematic model of a large sample of RR Lyrae detected by $Gaia$. By imposing a four-fold symmetry and employing $Gaia$ proper motions, we are able to infer the behaviour of the velocity ellipsoid between $\approx3$ and $\approx30$ kpc from the centre of the Galaxy. We detect the presence of two distinct components: a dominant non-rotating halo-like population and a much smaller rotating disc-like population. We demonstrate that the halo RR Lyrae can be described as a superposition of an isotropic and radially-biased parts. The radially-biased portion of the halo is characterised by a high orbital anisotropy $β\approx0.9$ and contributes between 50\% and 80\% of the halo RR Lyrae at $5<R$(kpc)$<25$. In line with previous studies, we interpret this high-$β$ component as the debris cloud of the ancient massive merger also known as the $Gaia$ Sausage (GS) whose orbital extrema we constrain. The lightcurve properties of the RR Lyrae support the kinematic decomposition: the GS stars are more metal-rich and boast higher fractions of Oosterhoff Type 1 and high amplitude short period (HASP) variables compared to the isotropic halo component. The metallicity/HASP maps reveal that the inner 10 kpc of the halo is likely inhabited by the RR Lyrae born in-situ. The mean azimuthal speed and the velocity dispersion of the disc RR Lyrae out to $R\approx30$ kpc are consistent with the behaviour of a young and metal-rich thin disc stellar population. △ Less

Submitted 2 January, 2021; v1 submitted 5 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS. The data underlying this article are available in Zenodo at http://doi.org/10.5281/zenodo.3972287

Evidence for supernova feedback sustaining gas turbulence in nearby star-forming galaxies

Authors: Cecilia Bacchini, Filippo Fraternali, Giuliano Iorio, Gabriele Pezzulli, Antonino Marasco, Carlo Nipoti

Abstract: It is well known that gas in galaxy discs is highly turbulent, but there is much debate on which mechanism can energetically maintain this turbulence. Among the possible candidates, supernova (SN) explosions are likely the primary drivers but doubts remain on whether they can be sufficient in regions of moderate star formation activity, in particular in the outer parts of discs. In this paper, we… ▽ More It is well known that gas in galaxy discs is highly turbulent, but there is much debate on which mechanism can energetically maintain this turbulence. Among the possible candidates, supernova (SN) explosions are likely the primary drivers but doubts remain on whether they can be sufficient in regions of moderate star formation activity, in particular in the outer parts of discs. In this paper, we measure the SN efficiency $η$, namely the fraction of the total SN energy needed to sustain turbulence in galaxies, and verify that SNe can indeed be the sole driving mechanism. The key novelty of our approach is that we take into account the increased turbulence dissipation timescale associated to the flaring in outer regions of gaseous discs. We analyse the distribution and kinematics of HI and CO in 10 nearby star-forming galaxies to obtain the radial profiles of the kinetic energy per unit area, for both the atomic gas and the molecular gas. We use a theoretical model to reproduce the observed energy with the sum of turbulent energy from SNe, as inferred from the observed star formation rate (SFR) surface density, and the gas thermal energy. We find that the observed kinetic energy is remarkably well reproduced by our model across the whole extent of the galactic discs, assuming $η$ constant with the galactocentric radius. Taking into account the uncertainties on the SFR surface density and on the atomic gas phase, we obtain that the median SN efficiencies for our sample of galaxies are $\langle η_\mathrm{atom} \rangle=0.015_{-0.008}^{+0.018}$ for the atomic gas and $\langle η_\mathrm{mol} \rangle = 0.003_{-0.002}^{+0.006}$ for the molecular gas. We conclude that SNe alone can sustain gas turbulence in nearby galaxies with only few percent of their energy and that there is essentially no need for any further source of energy. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: 18 pages, 6 figures, and 2 appendices; accepted by Astronomy & Astrophysics; abridged abstract

Journal ref: A&A 641, A70 (2020)

Unresolved stellar companions with Gaia DR2 astrometry

Authors: Vasily Belokurov, Zephyr Penoyre, Semyeong Oh, Giuliano Iorio, Simon Hodgkin, N. Wyn Evans, Andrew Everall, Sergey E. Koposov, Christopher A. Tout, Robert Izzard, Cathie J. Clarke, Anthony G. A. Brown

Abstract: For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to… ▽ More For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: submitted to MNRAS, comments welcome

Dynamics of black hole - neutron star binaries in young star clusters

Authors: Sara Rastello, Michela Mapelli, Ugo N. Di Carlo, Nicola Giacobbo, Filippo Santoliquido, Mario Spera, Alessandro Ballone, Giuliano Iorio

Abstract: Young star clusters are likely the most common birthplace of massive stars across cosmic time and influence the formation of compact binaries in several ways. Here, we simulate the formation of black hole -- neutron star binaries (BHNSs) in young star clusters, by means of the binary population synthesis code \texttt{MOBSE} interfaced with the $N$-body code \texttt{NBODY6++GPU}. BHNSs formed in yo… ▽ More Young star clusters are likely the most common birthplace of massive stars across cosmic time and influence the formation of compact binaries in several ways. Here, we simulate the formation of black hole -- neutron star binaries (BHNSs) in young star clusters, by means of the binary population synthesis code \texttt{MOBSE} interfaced with the $N$-body code \texttt{NBODY6++GPU}. BHNSs formed in young star clusters (dynamical BHNSs) are significantly more massive than BHNSs formed from isolated binaries (isolated BHNSs): $\sim{}40$~\% of the dynamical BHNS mergers have total mass $>15$ M$_\odot$, while only $\sim{}0.01$~\% of the isolated BHNS mergers have mass in excess of this value. Hence, our models strongly support a dynamical formation scenario for GW190814, given its total mass $\sim{}26$ M$_\odot$, if this event is a BHNS merger. All our dynamical BHNSs are ejected from their parent star cluster before they reach coalescence. Thus, a significant fraction of BHNS mergers occurring in the field might have originated in a young star cluster. The mass spectrum of BHNS mergers from gravitational-wave detections will provide a clue to differentiate between dynamical and isolated formation of BHNSs. △ Less

Submitted 30 June, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Comments: 9 pages, 4 figures, 2 tables, resubmitted after addressing referee's comments

The Tucana dwarf spheroidal galaxy: not such a massive failure after all

Authors: S. Taibi, G. Battaglia, M. Rejkuba, R. Leaman, N. Kacharov, G. Iorio, P. Jablonka, M. Zoccali

Abstract: Isolated Local Group (LG) dwarf galaxies have evolved most or all of their life unaffected by interactions with the large LG spirals and therefore offer the opportunity to learn about the intrinsic characteristics of this class of objects. Here we explore the internal kinematic and metallicity properties of one of the three isolated LG dwarf spheroidal galaxies, i.e. the Tucana dSph. This is an in… ▽ More Isolated Local Group (LG) dwarf galaxies have evolved most or all of their life unaffected by interactions with the large LG spirals and therefore offer the opportunity to learn about the intrinsic characteristics of this class of objects. Here we explore the internal kinematic and metallicity properties of one of the three isolated LG dwarf spheroidal galaxies, i.e. the Tucana dSph. This is an intriguing system, as it has been found in the literature to have an internal rotation of up to 16 km/s, a much higher velocity dispersion than other dSphs of similar luminosity, and a possible exception to the too-big-too-fail problem. We present results for a new VLT/FORS2 spectroscopic dataset in the CaII triplet region for 50 candidate red giant branch stars in the direction of Tucana, which yielded line-of-sight velocity and metallicity ([Fe/H]) measurements of 39 effective members. This doubles the number of Tucana's stars with such measurements. In addition, we re-reduce and include in our analysis the other two spectroscopic datasets presented in the literature, the VLT/FORS2 sample by Fraternali et al. (2009) and the VLT/FLAMES one by Gregory et al. (2019). We measure a systemic velocity of $180.0\pm1.3$ km/s, consistently across the various datasets analyzed, and find that a dispersion-only model is moderately favored over models accounting also for internal rotation. Our best estimate of the internal velocity dispersion is $6.2_{-1.3}^{+1.6}$ km/s, much smaller than the values reported in the literature and in line with similarly luminous dSphs; this is consistent with Tucana not being an exception to the too-big-to-fail problem, nor living in a dark matter halo much more massive than those of its siblings. As for the metallicity properties, we do not find anything unusual; there are hints of the presence of a [Fe/H] gradient but more data are needed to pin its presence down. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 26 pages, 9 figures, accepted to A&A

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

Kinematic and metallicity properties od the Aquarius dwarf galaxy from FORS2 MXU spectroscopy

Authors: L. Hermosa Muñoz, S. Taibi, G. Battaglia, G. Iorio, M. Rejkuba, R. Leaman, A. A. Cole, M. Irwin, P. Jablonka, N. Kacharov, A. McConnachie, E. Starkenburg, E. Tolstoy

Abstract: Dwarf galaxies found in isolation in the Local Group (LG) are unlikely to have interacted with the large LG spirals, and therefore environmental effects should not be the main drivers of their evolution. We aim to provide insight into the internal mechanisms shaping LG dwarf galaxies by increasing our knowledge of the internal properties of isolated systems. We focus on the evolved stellar compone… ▽ More Dwarf galaxies found in isolation in the Local Group (LG) are unlikely to have interacted with the large LG spirals, and therefore environmental effects should not be the main drivers of their evolution. We aim to provide insight into the internal mechanisms shaping LG dwarf galaxies by increasing our knowledge of the internal properties of isolated systems. We focus on the evolved stellar component of the Aquarius dwarf, whose kinematic and metallicity properties have only recently started to be explored. We have obtained spectroscopic data in the the near-infrared CaII triplet lines region with FORS2 at the Very Large Telescope for 53 red giant branch (RGB) stars, to derive line-of-sight velocities and [Fe/H] of the individual RGB stars. We have derived a systemic velocity of $-142.2^{+1.8}_{-1.8}$ km s$^{-1}$, in agreement with previous measurements from both the HI gas and stars. The internal kinematics of Aquarius appears to be best modelled by a combination of random motions (l.o.s. velocity dispersion of $10.3^{+1.6}_{-1.3}$ km s$^{-1}$) and linear rotation (with a gradient $-5.0^{+1.6}_{-1.9}$ km s$^{-1}$ arcmin$^{-1}$) along a P.A.=$139_{-27}^{+17}$ deg, consistent with the optical projected major axis. This rotation signal is significantly misaligned or even counter-rotating to that derived from the HI gas. We also find the tentative presence of a mild negative metallicity gradient and indications that the metal-rich stars have a colder velocity dispersion than the metal-poor ones. This work represents a significant improvement with respect to previous measurements of the RGB stars of Aquarius, as it doubles the number of member stars studied in the literature. We speculate that the misaligned rotation between the HI gas and evolved stellar component might have been the result of recent accretion of HI gas or re-accretion after gas-loss due to internal stellar feedback. △ Less

Submitted 29 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A. 15 pages, 9 figures

Journal ref: A&A 634, A10 (2020)

The volumetric star formation law in the Milky Way

Authors: C. Bacchini, F. Fraternali, G. Pezzulli, A. Marasco, G. Iorio, C. Nipoti

Abstract: Several open questions on galaxy formation and evolution have their roots in the lack of a universal star formation law, that could univocally link the gas properties, e.g. its density, to the star formation rate (SFR) density. In a recent paper, we used a sample of nearby disc galaxies to infer the volumetric star formation (VSF) law, a tight correlation between the gas and the SFR volume densiti… ▽ More Several open questions on galaxy formation and evolution have their roots in the lack of a universal star formation law, that could univocally link the gas properties, e.g. its density, to the star formation rate (SFR) density. In a recent paper, we used a sample of nearby disc galaxies to infer the volumetric star formation (VSF) law, a tight correlation between the gas and the SFR volume densities derived under the assumption of hydrostatic equilibrium for the gas disc. However, due to the dearth of information about the vertical distribution of the SFR in these galaxies, we could not find a unique slope for the VSF law, but two alternative values. In this paper, we use the scale height of the SFR density distribution in our Galaxy adopting classical Cepheids (age$\lesssim 200$ Myr) as tracers of star formation. We show that this latter is fully compatible with the flaring scale height expected from gas in hydrostatic equilibrium. These scale heights allowed us to convert the observed surface densities of gas and SFR into the corresponding volume densities. Our results indicate that the VSF law $ρ_\mathrm{SFR} \propto ρ_\mathrm{gas}^α$ with $α\approx 2$ is valid in the Milky Way as well as in nearby disc galaxies. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: 9 pages, 4 figures, + 4 pages of Appendix, accepted by A&A

Journal ref: A&A 632, A127 (2019)

Kinematics of the Palomar 5 stellar stream from RR Lyrae stars

Authors: Adrian M. Price-Whelan, Cecilia Mateu, Giuliano Iorio, Sarah Pearson, Ana Bonaca, Vasily Belokurov

Abstract: Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and su… ▽ More Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically-luminous standard candles and their distances can be inferred with small relative precision ($\sim3\%$). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, $d = 20.6 \pm 0.2~\textrm{kpc}$. We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work. △ Less

Submitted 5 February, 2020; v1 submitted 1 October, 2019; originally announced October 2019.

Comments: 13 pages, 4 figures. Published in AJ

The effect of tides on the Sculptor dwarf spheroidal galaxy

Authors: G. Iorio, C. Nipoti, G. Battaglia, A. Sollima

Abstract: Dwarf spheroidal galaxies (dSphs) appear to be some of the most dark matter dominated objects in the Universe. Their dynamical masses are commonly derived using the kinematics of stars under the assumption of equilibrium. However, these objects are satellites of massive galaxies (e.g.\ the Milky Way) and thus can be influenced by their tidal fields. We investigate the implication of the assumption… ▽ More Dwarf spheroidal galaxies (dSphs) appear to be some of the most dark matter dominated objects in the Universe. Their dynamical masses are commonly derived using the kinematics of stars under the assumption of equilibrium. However, these objects are satellites of massive galaxies (e.g.\ the Milky Way) and thus can be influenced by their tidal fields. We investigate the implication of the assumption of equilibrium focusing on the Sculptor dSph by means of ad-hoc $N$-body simulations tuned to reproduce the observed properties of Sculptor following the evolution along some observationally motivated orbits in the Milky Way gravitational field. For this purpose, we used state-of-the-art spectroscopic and photometric samples of Sculptor's stars. We found that the stellar component of the simulated object is not directly influenced by the tidal field, while $\approx 30\%-60\%$ the mass of the more diffuse DM halo is stripped. We conclude that, considering the most recent estimate of the Sculptor proper motion, the system is not affected by the tides and the stellar kinematics represents a robust tracer of the internal dynamics. In the simulations that match the observed properties of Sculptor, the present-day dark-to-luminous mass ratio is $\approx 6$ within the stellar half-light radius ($\approx0.3$ kpc) and $>50$ within the maximum radius of the analysed dataset ($\approx1.5^\circ\approx2$ kpc). △ Less

Submitted 28 June, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 19 pages, 10 figures, accepted for publication in MNRAS. V3: updated after editor comments See our playlist for simulation videos: https://av.tib.eu/series/633/supplemental+videos+of+the+paper+the+effect+of+tides+on+the+sculptor+dwarf+spheroidal+galaxy

Evidence for Two Early Accretion Events That Built the Milky Way Stellar Halo

Authors: G. C. Myeong, E. Vasiliev, G. Iorio, N. W. Evans, V. Belokurov

Abstract: The Gaia Sausage is the major accretion event that built the stellar halo of the Milky Way galaxy. Here, we provide dynamical and chemical evidence for a second substantial accretion episode, distinct from the Gaia Sausage. The Sequoia Event provided the bulk of the high energy retrograde stars in the stellar halo, as well as the recently discovered globular cluster FSR 1758. There are up to 6 fur… ▽ More The Gaia Sausage is the major accretion event that built the stellar halo of the Milky Way galaxy. Here, we provide dynamical and chemical evidence for a second substantial accretion episode, distinct from the Gaia Sausage. The Sequoia Event provided the bulk of the high energy retrograde stars in the stellar halo, as well as the recently discovered globular cluster FSR 1758. There are up to 6 further globular clusters, including $ω$~Centauri, as well as many of the retrograde substructures in Myeong et al. (2018), associated with the progenitor dwarf galaxy, named the Sequoia. The stellar mass in the Sequoia galaxy is $\sim 5 \times 10^{7} M_\odot$, whilst the total mass is $\sim 10^{10} M_\odot$, as judged from abundance matching or from the total sum of the globular cluster mass. Although clearly less massive than the Sausage, the Sequoia has a distinct chemo-dynamical signature. The strongly retrograde Sequoia stars have a typical eccentricity of $\sim0.6$, whereas the Sausage stars have no clear net rotation and move on predominantly radial orbits. On average, the Sequoia stars have lower metallicity by $\sim 0.3$ dex and higher abundance ratios as compared to the Sausage. We conjecture that the Sausage and the Sequoia galaxies may have been associated and accreted at a comparable epoch. △ Less

Submitted 10 June, 2019; v1 submitted 4 April, 2019; originally announced April 2019.

Comments: MNRAS, revised version in light of referee's comments. Minor modifications and clarifications to the argument, as well as an updated list of Sausage Globular Clusters

Volumetric star formation laws of disc galaxies

Authors: Cecilia Bacchini, Filippo Fraternali, Giuliano Iorio, Gabriele Pezzulli

Abstract: Star formation (SF) laws are fundamental relations between the gas content of a galaxy and its star formation rate (SFR) and play key roles in galaxy evolution models. In this paper, we present new empirical SF laws of disc galaxies based on volume densities. Following the assumption of hydrostatic equilibrium, we calculated the radial growth of the thickness of the gaseous discs in the combined g… ▽ More Star formation (SF) laws are fundamental relations between the gas content of a galaxy and its star formation rate (SFR) and play key roles in galaxy evolution models. In this paper, we present new empirical SF laws of disc galaxies based on volume densities. Following the assumption of hydrostatic equilibrium, we calculated the radial growth of the thickness of the gaseous discs in the combined gravitational potential of dark matter, stars, and gas for 12 nearby star-forming galaxies. This allowed us to convert the observed surface densities of gas and SFR into the deprojected volume densities. We found a tight correlation with slope in the range 1.3-1.9 between the volume densities of gas (HI+H$_2$) and the SFR with a significantly smaller scatter than the surface-based (Kennicutt) law and no change in the slope over five orders of magnitude. This indicates that taking into account the radial increase of the thickness of galaxy discs is crucial to reconstruct their three-dimensional density profiles, in particular in their outskirts. Moreover, our result suggests that the break in the slope seen in the Kennicutt law is due to disc flaring rather than to a drop of the SF efficiency at low surface densities. Surprisingly, we discovered an unexpected correlation between the volume densities of HI and SFR, indicating that the atomic gas is a good tracer of the cold star-forming gas, especially in low density HI-dominated environments. △ Less

Submitted 17 December, 2018; v1 submitted 8 October, 2018; originally announced October 2018.

Comments: 24 pages, 12 figures, accepted by A&A

Journal ref: A&A 622, A64 (2019)

The shape of the Galactic halo with $Gaia$ DR2 RR Lyrae. Anatomy of an ancient major merger

Authors: G. Iorio, V. Belokurov

Abstract: We use the $Gaia$ DR2 RR Lyrae sample to gain an uninterrupted view of the Galactic stellar halo. We dissect the available volume in slices parallel to the Milky Way's disc to show that within $\sim30$ kpc from the Galactic centre the halo is triaxial, with the longest axis misaligned by $\sim70^{\circ}$ with respect to the Galactic $x$-axis. This anatomical procedure exposes two large diffuse ove… ▽ More We use the $Gaia$ DR2 RR Lyrae sample to gain an uninterrupted view of the Galactic stellar halo. We dissect the available volume in slices parallel to the Milky Way's disc to show that within $\sim30$ kpc from the Galactic centre the halo is triaxial, with the longest axis misaligned by $\sim70^{\circ}$ with respect to the Galactic $x$-axis. This anatomical procedure exposes two large diffuse over-densities aligned with the semi-major axis of the halo: the Hercules-Aquila Cloud and the Virgo Over-density. We reveal the kinematics of the entire inner halo by mapping out the amplitudes and directions of the RR Lyrae proper motions. These are then compared to simple models with different anisotropies to demonstrate that the inner halo is dominated by stars on highly eccentric orbits. We interpret the shape of the density and the kinematics of the $Gaia$ DR2 RR Lyrae as evidence in favour of a scenario in which the bulk of the halo was deposited in a single massive merger event. △ Less

Submitted 15 October, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

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

S0 galaxies are faded spirals: clues from their angular momentum content

Authors: Francesca Rizzo, Filippo Fraternali, Giuliano Iorio

Abstract: The distribution of galaxies in the stellar specific angular momentum versus stellar mass plane ($j_{\star}$-$M_{\star}$) provides key insights into their formation mechanisms. In this paper, we determine the location in this plane of a sample of ten field/group unbarred lenticular (S0) galaxies from the CALIFA survey. We performed a bulge-disc decomposition both photometrically and kinematically… ▽ More The distribution of galaxies in the stellar specific angular momentum versus stellar mass plane ($j_{\star}$-$M_{\star}$) provides key insights into their formation mechanisms. In this paper, we determine the location in this plane of a sample of ten field/group unbarred lenticular (S0) galaxies from the CALIFA survey. We performed a bulge-disc decomposition both photometrically and kinematically to study the stellar specific angular momentum of the disc components alone and understand the evolutionary links between S0s and other Hubble types. We found that eight of our S0 discs have a distribution in the $j_{\star}$-$M_{\star}$ plane that is fully compatible with that of spiral discs, while only two have values of $j_{\star}$ lower than the spirals. These two outliers show signs of recent merging. Our results suggest that merger and interaction processes are not the dominant mechanisms in S0 formation in low-density environments. Instead, S0s appear to be the result of secular processes and the fading of spiral galaxies after the shutdown of star formation. △ Less

Submitted 8 February, 2018; v1 submitted 28 November, 2017; originally announced November 2017.

Comments: 35 pages, 22 figures. Accepted for publication in MNRAS

A Magellanic origin for the Virgo substructure

Authors: D. Boubert, V. Belokurov, D. Erkal, G. Iorio

Abstract: Iorio et al. (2018) mapped out the Milky Way halo using a sample of RR Lyrae stars drawn from a cross-match of Gaia with 2MASS. We investigate the significant residual in their model which we constrain to lie at Galactocentric radii $12<R<27\;\mathrm{kpc}$ and extend over $2600\;\mathrm{deg}^2$ of the sky. A counterpart of this structure exists in both the Catalina Real Time Survey and the sample… ▽ More Iorio et al. (2018) mapped out the Milky Way halo using a sample of RR Lyrae stars drawn from a cross-match of Gaia with 2MASS. We investigate the significant residual in their model which we constrain to lie at Galactocentric radii $12<R<27\;\mathrm{kpc}$ and extend over $2600\;\mathrm{deg}^2$ of the sky. A counterpart of this structure exists in both the Catalina Real Time Survey and the sample of RR Lyrae variables identified in Pan-STARRS by Hernitschek et al. (2016), demonstrating that this structure is not caused by the spatial inhomogeneity of Gaia. The structure is likely the Virgo Stellar Stream and/or Virgo Over-Density. We show the structure is aligned with the Magellanic Stream and suggest that it is either debris from a disrupted dwarf galaxy that was a member of the Vast Polar Structure or that it is SMC debris from a tidal interaction of the SMC and LMC $3\;\mathrm{Gyr}$ ago. If the latter then the sub-structure in Virgo may have a Magellanic origin. △ Less

Submitted 31 October, 2018; v1 submitted 20 November, 2017; originally announced November 2017.

Comments: 9 pages, 7 figures, accepted to MNRAS 31/10/2018

The first all-sky view of the Milky Way stellar halo with Gaia+2MASS RR Lyrae

Authors: G. Iorio, V. Belokurov, D. Erkal, S. E. Koposov, C. Nipoti, F. Fraternali

Abstract: We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo (… ▽ More We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $α=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre. △ Less

Submitted 2 November, 2017; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS

LITTLE THINGS in 3D: robust determination of the circular velocity of dwarf irregular galaxies

Authors: Giuliano Iorio, Filippo Fraternali, Carlo Nipoti, Enrico Di Teodoro, Justin I. Read, Giuseppina Battaglia

Abstract: Dwarf Irregular galaxies (dIrrs) are the smallest stellar systems with extended HI discs. The study of the kinematics of such discs is a powerful tool to estimate the total matter distribution at these very small scales. In this work, we study the HI kinematics of 17 galaxies extracted from the `Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey' (LITTLE THINGS). Our appr… ▽ More Dwarf Irregular galaxies (dIrrs) are the smallest stellar systems with extended HI discs. The study of the kinematics of such discs is a powerful tool to estimate the total matter distribution at these very small scales. In this work, we study the HI kinematics of 17 galaxies extracted from the `Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey' (LITTLE THINGS). Our approach differs significantly from previous studies in that we directly fit 3D models (two spatial dimensions plus one spectral dimension) using the software $^\text{3D}$BAROLO, fully exploiting the information in the HI datacubes. For each galaxy we derive the geometric parameters of the HI disc (inclination and position angle), the radial distribution of the surface density, the velocity-dispersion ($σ_v$) profile and the rotation curve. The circular velocity (V$_{\text{c}}$), which traces directly the galactic potential, is then obtained by correcting the rotation curve for the asymmetric drift. As an initial application, we show that these dIrrs lie on a baryonic Tully-Fisher relation in excellent agreement with that seen on larger scales. The final products of this work are high-quality, ready-to-use kinematic data ($\textrm{V}_\textrm{c}$ and $σ_v$) that we make publicly available. These can be used to perform dynamical studies and improve our understanding of these low-mass galaxies. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 36 pages, 28 figures, 2 tables. Submitted to MNRAS (revised version after the referee report). The final rotation curves can be downloaded from http://www.filippofraternali.com/styled-9/index.html

Journal ref: Mon Not R Astron Soc (2017) 466 (4): 4159-4192

The stellar mass-halo mass relation of isolated field dwarfs: a critical test of $Λ$CDM at the edge of galaxy formation

Authors: J. I. Read, G. Iorio, O. Agertz, F. Fraternali

Abstract: We fit the rotation curves of isolated dwarf galaxies to directly measure the stellar mass-halo mass relation ($M_*-M_{200}$) over the mass range $5 \times 10^5 < M_{*}/{\rm M}_\odot < 10^{8}$. By accounting for cusp-core transformations due to stellar feedback, we find a monotonic relation with little scatter. Such monotonicity implies that abundance matching should yield a similar $M_*-M_{200}$… ▽ More We fit the rotation curves of isolated dwarf galaxies to directly measure the stellar mass-halo mass relation ($M_*-M_{200}$) over the mass range $5 \times 10^5 < M_{*}/{\rm M}_\odot < 10^{8}$. By accounting for cusp-core transformations due to stellar feedback, we find a monotonic relation with little scatter. Such monotonicity implies that abundance matching should yield a similar $M_*-M_{200}$ if the cosmological model is correct. Using the 'field galaxy' stellar mass function from the Sloan Digital Sky Survey (SDSS) and the halo mass function from the $Λ$ Cold Dark Matter Bolshoi simulation, we find remarkable agreement between the two. This holds down to $M_{200} \sim 5 \times 10^9$M$_\odot$, and to $M_{200} \sim 5 \times 10^8$M$_\odot$ if we assume a power law extrapolation of the SDSS stellar mass function below $M_* \sim 10^7$M$_\odot$. However, if instead of SDSS we use the stellar mass function of nearby galaxy groups, then the agreement is poor. This occurs because the group stellar mass function is shallower than that of the field below $M_* \sim 10^9$M$_\odot$, recovering the familiar 'missing satellites' and 'too big to fail' problems. Our result demonstrates that both problems are confined to group environments and must, therefore, owe to 'galaxy formation physics' rather than exotic cosmology. Finally, we repeat our analysis for a $Λ$ Warm Dark Matter cosmology, finding that it fails at 68% confidence for a thermal relic mass of $m_{\rm WDM} < 1.25$keV, and $m_{\rm WDM} < 2$keV if we use the power law extrapolation of SDSS. We conclude by making a number of predictions for future surveys based on these results. △ Less

Submitted 8 February, 2017; v1 submitted 11 July, 2016; originally announced July 2016.

Comments: 22 pages; 2 Tables; 10 Figures. This is the version accepted for publication in MNRAS. Key changes: (i) added substantially more information on the surveys used to measure the stellar mass functions; (ii) added tests of the robustness of our results. Results and conclusions unchanged from previously. Minor typos corrected from previous version

Understanding the shape and diversity of dwarf galaxy rotation curves in LCDM

Authors: J. I. Read, G. Iorio, O. Agertz, F. Fraternali

Abstract: The shape and diversity of dwarf galaxy rotation curves is at apparent odds with dark matter halos in a $Λ$ Cold Dark Matter ($Λ$CDM) cosmology. We use mock data from isolated dwarf galaxy simulations to show that this owes to three main effects. Firstly, stellar feedback heats dark matter, leading to a 'coreNFW' dark matter density profile with a slowly rising rotation curve. Secondly, if close t… ▽ More The shape and diversity of dwarf galaxy rotation curves is at apparent odds with dark matter halos in a $Λ$ Cold Dark Matter ($Λ$CDM) cosmology. We use mock data from isolated dwarf galaxy simulations to show that this owes to three main effects. Firstly, stellar feedback heats dark matter, leading to a 'coreNFW' dark matter density profile with a slowly rising rotation curve. Secondly, if close to a recent starburst, large HI bubbles push the rotation curve out of equilibrium, deforming the rotation curve shape. Thirdly, when galaxies are viewed near face-on, their best fit inclination is biased high. This can lead to a very shallow rotation curve that falsely implies a large dark matter core. All three problems can be avoided, however, by a combination of improved mass models and a careful selection of target galaxies. Fitting our coreNFW model to mock rotation curve data, we show that we can recover the rotation curve shape, dark matter halo mass $M_{200}$ and concentration parameter $c$ within our quoted uncertainties. We fit our coreNFW model to real data for four isolated dwarf irregulars, chosen to span a wide range of rotation curve shapes. We obtain an excellent fit for NGC 6822 and WLM, with tight constraints on $M_{200}$, and $c$ consistent with $Λ$CDM. However, IC 1613 and DDO 101 give a poor fit. For IC 1613, we show that this owes to disequilibria and its uncertain inclination $i$; for DDO 101, it owes to its uncertain distance $D$. If we assume $i_{\rm IC1613} \sim 15^\circ$ and $D_{\rm DDO101} \sim 12$ Mpc, consistent with current uncertainties, we are able to fit both galaxies very well. We conclude that $Λ$CDM appears to give an excellent match to dwarf galaxy rotation curves. △ Less

Submitted 1 July, 2016; v1 submitted 21 January, 2016; originally announced January 2016.

Comments: 21 pages; 1 Table; 11 Figures. New section added discussing how our results relate to the Baryonic Tully Fisher Relation; appendices added exploring our concentration parameter prior, and allowing the dark matter core size to freely vary for WLM; all conclusions unchanged. To appear in MNRAS