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.

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 $\eta$, 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 $\eta$ 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 \eta_\mathrm{atom} \rangle=0.015_{-0.008}^{+0.018}$ for the atomic gas and $\langle \eta_\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.

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.

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.

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.

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.

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 $\rho_\mathrm{SFR} \propto \rho_\mathrm{gas}^\alpha$ with $\alpha \approx 2$ is valid in the Milky Way as well as in nearby disc galaxies.

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

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 $\omega$~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.

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.

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.

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.

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.

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 $\alpha=-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.

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 $\Lambda$ 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 $\Lambda$ 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.

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 ($\sigma_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 $\sigma_v$) that we make publicly available. These can be used to perform dynamical studies and improve our understanding of these low-mass galaxies.

The shape and diversity of dwarf galaxy rotation curves is at apparent odds with dark matter halos in a $\Lambda$ Cold Dark Matter ($\Lambda$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 $\Lambda$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 $\Lambda$CDM appears to give an excellent match to dwarf galaxy rotation curves.