Stellar cluster formation in a Milky Way-sized galaxy at z>4 -- I. The proto-globular cluster population and the imposter amongst us

The formation history of globular clusters (GCs) at redshift $z > 4$ remains an unsolved problem. In this work, we use the cosmological, $N$-body hydrodynamical ``zoom-in'' simulation GigaEris to study the properties and formation of proto-GC candidates in the region surrounding the progenitor of a Milky Way-sized galaxy. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion and allows to resolve systems at the scale of star clusters. We define proto-GC candidate systems as gravitationally bound stellar systems with baryonic mass fraction $F_{\rm b} \geq 0.75$ and stellar velocity dispersion $\sigma_{\star} < 20$ km s$^{-1}$. At $z=4.4$ we identify 9 systems which satisfy our criteria, all of which form between 10 kpc to 30 kpc from the centre of the main host. Their baryonic masses are in the range $10^5$- $10^7$ M$_{\odot}$. By the end of the simulation, they still have a relatively low stellar mass ($M_{\star} \sim 10^4$--$10^5$ M$_{\odot}$) and a metallicity ($-1.8 \lesssim {\rm [Fe/H]} \lesssim -0.8$) similar to the blue Galactic GCs. All of the identified systems except one appear to be associated with gas filaments accreting onto the main galaxy in the circum-galactic region, and formed at $z=5-4$. The exception is the oldest object, which appears to be a stripped compact dwarf galaxy that has interacted with the main halo between $z = 5.8$ and $z=5.2$ and has lost its entire dark matter content due to tidal mass loss.Comment: 11 pages, 8 figures, accepted for publication in MNRA

The Dawn of Disk Formation in a Milky Way-sized Galaxy Halo: Thin Stellar Disks at z > 4

We present results from GigaEris, a cosmological, N-body hydrodynamical "zoom-in" simulation of the formation of a Milky Way-sized galaxy halo with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We focus on the early assembly of the galaxy, down to redshift z = 4.4. The simulated galaxy has properties consistent with extrapolations of the main sequence of star-forming galaxies to higher redshifts and levels off to a star formation rate of ∼60 M⊙ yr−1 at z = 4.4. A compact, thin rotating stellar disk with properties analogous to those of low-redshift systems arises already at z ∼ 8. The galaxy rapidly develops a multi-component structure, and the disk, at least at these early stages, does not grow "upside-down" as often reported in the literature. Rather, at any given time, newly born stars contribute to sustain a thin disk. The kinematics reflect the early, ubiquitous presence of a thin disk, as a stellar disk component with vϕ/σR larger than unity is already present at z ∼ 9–10. Our results suggest that high-resolution spectro-photometric observations of very high-redshift galaxies should find thin rotating disks, consistent with the recent discovery of cold rotating gas disks by ALMA. Finally, we present synthetic images for the James Webb Space Telescope NIRCam camera, showing how the early disk would be easily detectable already at those early times

Astrophysics with the Laser Interferometer Space Antenna

Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy as it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and other space-based instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed: ultra-compact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help make progress in the different areas. New research avenues that LISA itself, or its joint exploitation with studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe

Black holes, gravitational waves and fundamental physics:a roadmap

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress

A Subgrid Model for the Growth of Dust Particles in Hydrodynamical Simulations of Protoplanetary Disks

We present the first 2D hydrodynamical finite-volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a subgrid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust size, with and without an embedded Jupiter-mass planet that triggers gap formation. We use the two-fluid polar Godunov-type code RoSSBi developed by Surville et al. combined with a new local subgrid method for dust evolution based on the model by Birnstiel et al. We find striking differences between simulations with variable and fixed dust sizes. The timescales for dust depletion differ significantly and yield a completely different evolution of the dust surface density. In general, sharp features such as pileups of dust in the inner disk and near gap edges, when a massive planet is present, become much weaker. This has important implications for the interpretation of observed substructure in disks, suggesting that the presence of a massive planet does not necessarily cause sharp gaps and rings in the dust component. Also, particles with different dust sizes show a different distribution, pointing to the importance of multiwavelength synthetic observations in order to compare with observations by ALMA and other instruments. We also find that simulations adopting fixed intermediate particle sizes, in the range of 10−2 to 10−1 cm, best approximate the surface density evolution seen in simulations with dust evolution

The lifetime of binary black holes in Sérsic galaxy models

In the local Universe, black holes of $10^{5-6}\, {\rm M_\odot }$ are hosted in galaxies displaying a variety of stellar profiles and morphologies. These black holes are the anticipated targets of LISA, the Laser Interferometer Space Antenna that will detect the low-frequency gravitational-wave signal emitted by binary black holes in this mass interval. In this paper, we infer upper limits on the lifetime of binary black holes of $10^{5-6}\, {\rm M_\odot }$ and up to $10^8\, {\rm M_\odot }$, forming in galaxy mergers, exploring two underlying stellar density profiles, by Dehnen and by Prugniel &amp; Simien, and by exploiting local scaling relations between the mass of the black holes and several quantities of their hosts. We focus on the phase of the dynamical evolution when the binary is transitioning from the hardening phase ruled by the interaction with single stars to the phase driven by the emission of gravitational waves. We find that different stellar profiles predict very distinct trends with binary mass, with lifetimes ranging between fractions of a Gyr to more than 10 Gyr, and with a spread of about one order of magnitude, given by the uncertainties in the observed correlations, which are larger in the low-mass tail of the observed black hole population

Stellar cluster formation in a Milky Way-sized galaxy at z > 4 – I. The proto-globular cluster population and the imposter amongst us

The formation history of globular clusters (GCs) at redshift z &gt; 4 remains an unsolved problem. In this work, we use the cosmological, N-body hydrodynamical ‘zoom-in’ simulation GigaEris to study the properties and formation environment of proto-GC candidates in the region surrounding the progenitor of a Milky Way-sized galaxy. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We define proto-GC candidate systems as gravitationally bound stellar systems with baryonic mass fraction Fb ≥ 0.75 and stellar velocity dispersion σ⋆ &lt; 20 km s−1. At z = 4.4, we identify nine systems that satisfy our criteria, all of which form between 10 and 30 kpc from the centre of the main host. Their baryonic masses are in the range 105–107 M⊙. By the end of the simulation, they still have a relatively low stellar mass (M⋆ ∼ 104–105 M⊙) and a metallicity (−1.8 ≲ [Fe/H] ≲ −0.8) similar to the blue Galactic GCs. All of the identified systems except one appear to be associated with gas filaments accreting onto the main galaxy in the circum-galactic region and formed at z = 5–4. The exception is the oldest object, which appears to be a stripped compact dwarf galaxy that has interacted with the main halo between z = 5.8 and z = 5.2 and has lost its entire dark matter content due to tidal mass loss