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.

More and more observations indicate that young star clusters could retain imprints of their formation process. In particular, the degree of substructuring and rotation are possibly the direct result of the collapse of the parent molecular cloud from which these systems form. Such properties can, in principle, be washed-out, but they are also expected to have an impact on the relaxation of these systems. We ran and analyzed a set of ten hydrodynamical simulations of the formation of embedded star clusters through the collapse of turbulent massive molecular clouds. We systematically studied the fractality of our star clusters, showing that they are all extremely substructured (fractal dimension $D=1.0-1.8$). We also found that fractality is slowly reduced, with time, on small scales, while it persists on large scales on longer timescales. Signatures of rotation are found in different simulations at every time of the evolution, even for slightly supervirial substructures, proving that the parent molecular gas transfers part of its angular momentum to the new stellar systems.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

Next generation ground-based gravitational-wave detectors will observe binary black hole (BBH) mergers up to redshift $\gtrsim{}10$, probing the evolution of compact binary (CB) mergers across cosmic time. Here, we present a new data-driven model to estimate the cosmic merger rate density (MRD) evolution of CBs, by coupling catalogs of CB mergers with observational constraints on the cosmic star formation rate density and on the metallicity evolution of the Universe. We adopt catalogs of CB mergers derived from recent $N-$body and population-synthesis simulations, to describe the MRD of CBs formed in young star clusters (hereafter, dynamical CBs) and in the field (hereafter, isolated CBs). The local MRD of dynamical BBHs is $\mathcal{R}_{\rm BBH}=67^{+29}_{-23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with the 90% credible interval from the first and second observing run (O1 and O2) of the LIGO-Virgo collaboration, and with the local MRD of isolated BBHs ($\mathcal{R}_{\rm BBH}=49^{+79}_{-37}$ Gpc$^{-3}$ yr$^{-1}$). The local MRD of dynamical and isolated black hole-neutron star binaries (BHNSs) is $\mathcal{R}_{\rm BHNS}=38^{+32}_{-24}$ and $45^{+45}_{-32}$ Gpc$^{-3}$ yr$^{-1}$, respectively. Both values are consistent with the upper limit inferred from O1 and O2. Finally, the local MRD of dynamical binary neutron stars (BNSs, $\mathcal{R}_{\rm BNS}=150^{+56}_{-40}$ Gpc$^{-3}$ yr$^{-1}$) is a factor of two lower than the local MRD of BNSs formed in isolated binaries ($\mathcal{R}_{\rm BNS}=281^{+109}_{-75}$ Gpc$^{-3}$ yr$^{-1}$). The MRD for all CB classes grows with redshift, reaching its maximum at $z \in [1.5,2.5]$, and then decreases. This trend springs from the interplay between cosmic star formation rate, metallicity evolution and delay time of binary compact objects.

The origin of the molecular clouds orbiting SgrA* is one of the most debated questions about our Galactic Centre. Here, we present the hydrodynamic simulation of a molecular cloud infalling towards SgrA*, performed with the adaptive-mesh-refinement code RAMSES. Through such simulation, we propose that the circumnuclear disc and the +20 km/s cloud originated from the same tidal disruption episode, occurred less than 1 Myr ago. We also show that recent star formation is to be expected in the +20 km/s cloud, as also suggested by recent observations.

The high velocity gradient observed in the compact cloud CO-0.40-0.22, at a projected distance of 60 pc from the centre of the Milky Way, has led its discoverers to identify the closeby mm continuum emitter, CO-0.40-0.22*, with an intermediate mass black hole (IMBH) candidate. We describe the interaction between CO-0.40-0.22 and the IMBH, by means of a simple analytical model and of hydrodynamical simulations. Through such calculation, we obtain a lower limit to the mass of CO-0.40-0.22* of few $10^4 \times \; M_{\odot}$. This result tends to exclude the formation of such massive black hole in the proximity of the Galactic Centre. On the other hand, CO-0.40-0.22* might have been brought to such distances in cosmological timescales, if it was born in a dark matter halo or globular cluster around the Milky Way.

The so-called S2 star reached its closest approach to the massive black hole (BH) at around 1500 $R_\mathrm{s}$ in May 2018. It has been proposed that the interaction of its stellar wind with the high-density accretion flow at this distance from Sgr A* will lead to a detectable, month-long X-ray flare. Our goal is to verify whether or not the S2 star wind can be used as a diagnostic tool to infer the properties of the accretion flow towards Sgr A* at its pericentre (an unprobed distance regime), putting important constraints on BH accretion flow models. We run a series of three-dimensional adaptive mesh refinement simulations with the help of the Ramses code which include the realistic treatment of the interaction of S2's stellar wind with the accretion flow along its orbit and - apart from hydrodynamical and thermodynamical effects - include the tidal interaction with the massive BH. These are post-processed to derive the X-ray emission in the observable 2-10 keV window. No significant excess of X-ray emission from Sgr A* is found for typical accretion flow models. A measurable excess is produced for a significantly increased density of the accretion flow. This can, however, be ruled out for standard power-law accretion flow models as in this case the thermal X-ray emission without the S2 wind interaction would already exceed the observed quiescent luminosity. Only a significant change of the wind parameters (increased mass loss rate and decreased wind velocity) might lead to an (marginally) observable X-ray flaring event. Even the detection of an (month-long) X-ray flare during the pericentre passage of the S2 star would not allow for strict constraints to be put on the accretion flow around Sgr A* due to the degeneracy caused by the dependence on multiple parameters (of the accretion flow model as well as the stellar wind).

We investigate the formation of circumnuclear gas structures from the tidal disruption of molecular clouds in galactic nuclei, by means of smoothed particle hydrodynamics simulations. We model galactic nuclei as composed of a supermassive black hole (SMBH) and a nuclear star cluster (NSC) and consider different mass ratios between the two components. We find that the relative masses of the SMBH and the NSC have a deep impact on the morphology of the circumnuclear gas. Extended disks form only inside the sphere of influence of the SMBH. In contrast, compact rings naturally form outside the SMBH's sphere of influence, where the gravity is dominated by the NSC. This result is in agreement with the properties of the Milky Way's circumnuclear ring, which orbits outside the SMBH sphere of influence. Our results indicate that compact circumnuclear rings can naturally form outside the SMBH sphere of influence.

The nature of the gaseous and dusty cloud G2 in the Galactic Centre is still under debate. We present three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of G2, modeled as an outflow from a "compact source" moving on the observed orbit. The construction of mock position-velocity (PV) diagrams enables a direct comparison with observations and allow us to conclude that the observational properties of the gaseous component of G2 could be matched by a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \mathrm{yr^{-1}}$) and slow ($50 \;\mathrm{km \;s^{-1}}$) outflow, as observed for T Tauri stars. In order for this to be true, only the material at larger ($>100 \;\mathrm{AU}$) distances from the source must be actually emitting, otherwise G2 would appear too compact compared to the observed PV diagrams. On the other hand, the presence of a central dusty source might be able to explain the compactness of G2's dust component. In the present scenario, 5-10 years after pericentre the compact source should decouple from the previously ejected material, due to the hydrodynamic interaction of the latter with the surrounding hot and dense atmosphere. In this case, a new outflow should form, ahead of the previous one, which would be the smoking gun evidence for an outflow scenario.

We study the evolution of G2 in a \textit{Compact Source Scenario}, where G2 is the outflow from a low-mass central star moving on the observed orbit. This is done through 3D AMR simulations of the hydrodynamic interaction of G2 with the surrounding hot accretion flow. A comparison with observations is done by means of mock position-velocity (PV) diagrams. We found that a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \; \mathrm{yr^{-1}}$) and slow ($v_\mathrm{w}=50 \;\mathrm{km\; s^{-1}}$) outflow can reproduce G2's properties. A faster outflow ($v_\mathrm{w}=400 \;\mathrm{km\; s^{-1}}$) might also be able to explain the material that seems to follow G2 on the same orbit.

With the help of 3D AMR hydrodynamical simulations we aim at understanding G2's nature, recent evolution and fate in the coming years. By exploring the possible parameter space of the diffuse cloud scenario, we find that a starting point within the disc of young stars is favoured by the observations, which may hint at G2 being the result of stellar wind interactions.

Observations of the gas component of the cloud G2 in the Galactic Center have revealed its connection to a tail (G2t) lying on the same orbit. More recent studies indicate a connection between G2 and G1, another cloud detected on the blueshifted side of G2's orbit, suggesting a scenario in which G2 is a denser clump in a stream of gas. In this Letter we show that a simulation of an outflow by a central source (possibly a T Tauri star) moving on G2's orbit and interacting with a hot atmosphere surrounding SgrA* can have G2 and G2t as a byproduct. G2 would be the bow-shock formed in the head of the source, while G2t might be the result of the stripping of the rest of the shocked material by the ram pressure of the surrounding hot gas and of its successive accumulation in the trailing region. Mock position-velocity diagrams for the Br$\gamma$ emission for this simulation can indeed reproduce the correct position and velocity of G2t, as well as the more tenuous material in between. Though some tension between the observations and the simulated model remains, we argue that this might be due to issues in the construction of observed position-velocity (PV) diagrams and/or to a poor treatment of some physical processes - like hydrodynamic mixing - in our simulation.

Stellar winds and supernova (SN) explosions of massive stars ("stellar feedback") create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work we estimate the amount of feedback energy retained as kinetic energy when the bubble walls have decelerated to the sound speed of the ambient medium. We show that the feedback of the most massive star outweighs the feedback from less massive stars. For a giant molecular cloud (GMC) mass of 1e5 solar masses (as e.g. found in the Orion GMCs) and a star formation efficiency of 8% the initial mass function predicts a most massive star of approximately 60 solar masses. For this stellar evolution model we test the dependence of the retained kinetic energy of the cold GMC gas on the inclusion of stellar winds. In our model winds insert 2.34 times the energy of a SN and create stellar wind bubbles serving as pressure reservoirs. We find that during the pressure driven phases of the bubble evolution radiative losses peak near the contact discontinuity (CD), and thus, the retained energy depends critically on the scales of the mixing processes across the CD. Taking into account the winds of massive stars increases the amount of kinetic energy deposited in the cold ISM from 0.1% to a few percent of the feedback energy.

The gas cloud G2 is currently being tidally disrupted by the Galactic Centre super-massive black hole, Sgr A*. The region around the black hole is populated by $\sim 30$ Wolf-Rayet stars, which produce strong outflows. We explore the possibility that gas clumps, such as G2, originate from the collision of stellar winds via the non-linear thin shell instability. Following an analytical approach, we study the thermal evolution of slabs formed in the symmetric collision of winds, evaluating whether instabilities occur, and estimating possible clump masses. We find that the collision of relatively slow ($< 750$ km s$^{-1}$) and strong ($\sim 10^{-5}$ Msun yr$^{-1}$) stellar winds from stars at short separations ($<10$ mpc) is a process that indeed could produce clumps of G2's mass and above. Such short separation encounters of single stars along their known orbits are not common in the Galactic Centre, making this process a possible but unlikely origin for G2. We also discuss clump formation in close binaries such as IRS 16SW and in asymmetric encounters as promising alternatives that deserve further numerical study.

The dusty, ionized gas cloud G2 is currently passing the massive black hole in the Galactic Center at a distance of roughly 2400 Schwarzschild radii. We explore the possibility of a starting point of the cloud within the disks of young stars. We make use of the large amount of new observations in order to put constraints on G2's origin. Interpreting the observations as a diffuse cloud of gas, we employ three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations with the PLUTO code and do a detailed comparison with observational data. The simulations presented in this work update our previously obtained results in multiple ways: (1) high resolution three-dimensional hydrodynamical AMR simulations are used, (2) the cloud follows the updated orbit based on the Brackett-$\gamma$ data, (3) a detailed comparison to the observed high-quality position-velocity diagrams and the evolution of the total Brackett-$\gamma$ luminosity is done. We concentrate on two unsolved problems of the diffuse cloud scenario: the unphysical formation epoch only shortly before the first detection and the too steep Brackett-$\gamma$ light curve obtained in simulations, whereas the observations indicate a constant Brackett-$\gamma$ luminosity between 2004 and 2013. For a given atmosphere and cloud mass, we find a consistent model that can explain both, the observed Brackett-$\gamma$ light curve and the position-velocity diagrams of all epochs. Assuming initial pressure equilibrium with the atmosphere, this can be reached for a starting date earlier than roughly 1900, which is close to apo-center and well within the disks of young stars.

We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.

In 2011, we discovered a compact gas cloud ("G2") with roughly three Earth masses that is falling on a near-radial orbit toward the massive black hole in the Galactic Center. The orbit is well constrained and pericenter passage is predicted for early 2014. Our data beautifully show that G2 gets tidally sheared apart due to the massive black hole's force. During the next months, we expect that in addition to the tidal effects, hydrodynamics get important, when G2 collides with the hot ambient gas around Sgr A*. Simulations show that ultimately, the cloud's material might fall into the massive black hole. Predictions for the accretion rate and luminosity evolution, however, are very difficult due to the many unknowns. Nevertheless, this might be a unique opportunity in the next years to observe how gas feeds a massive black hole in a galactic nucleus.

The origin of the dense gas cloud G2 discovered in the Galactic Center (Gillessen et al. 2012) is still a debated puzzle. G2 might be a diffuse cloud or the result of an outflow from an invisible star embedded in it. We present here detailed simulations of the evolution of winds on G2's orbit. We find that the hydrodynamic interaction with the hot atmosphere present in the Galactic Center and the extreme gravitational field of the supermassive black hole must be taken in account when modeling such a source scenario. We find that the hydrodynamic interaction with the hot atmosphere present in the Galactic Center and the extreme gravitational field of the supermassive black hole must be taken in account when modeling such a source scenario. We also find that in this scenario most of the Br\gamma\ luminosity is expected to come from the highly filamentary densest shocked wind material. G2's observational properties can be used to constrain the properties of the outflow and our best model has a mass outflow rate of Mdot,w=8.8 x 10^{-8} Msun/yr and a wind velocity of vw = 50 km/s. These values are compatible with those of a young TTauri star wind, as already suggested by Scoville & Burkert (2013).

The origin of the dense gas cloud G2 discovered in the Galactic Center (Gillessen et al. 2012) is still a debated puzzle. G2 might be a diffuse cloud or the result of an outflow from an invisible star embedded in it. We present hydrodynamical simulations of the evolution of different spherically symmetric winds of a stellar object embedded in G2. We find that the interaction with the ambient medium and with the extreme gravitational field of the supermassive black hole in the Galactic Center must be taken into account for such a source scenario. The thermal pressure of the hot and dense atmosphere confines the wind, while its ram pressure shapes it via stripping along the orbit, with the details depending on the wind parameters. Tidal forces squeeze the wind near pericenter, reducing it to a thin and elongated filament. We also find that in this scenario most of the Br\gamma\ luminosity is expected to come from the densest part of the wind, which has a highly filamentary structure with low filling factor. The observations can be best matched by a mass outflow rate of Mdot,w=8.8 x 10^{-8} Msun/yr and a wind velocity of vw = 50 km/s. These values are compatible with those of a young T Tauri star wind, as already suggested by (Scoville & Burkert 2013).

We have further followed the evolution of the orbital and physical properties of G2, the object currently falling toward the massive black hole in the Galactic Center on a near-radial orbit. New, very sensitive data were taken in April 2013 with NACO and SINFONI at the ESO VLT . The 'head' of G2 continues to be stretched ever further along the orbit in position-velocity space. A fraction of its emission appears to be already emerging on the blue-shifted side of the orbit, past pericenter approach. Ionized gas in the head is now stretched over more than 15,000 Schwarzschild radii RS around the pericenter of the orbit, at ~ 2000 RS ~ 20 light hours from the black hole. The pericenter passage of G2 will be a process stretching over a period of at least one year. The Brackett-{\gamma} luminosity of the head has been constant over the past 9 years, to within +- 25%, as have the line ratios Brackett-{\gamma} / Paschen-{\alpha} and Brackett-{\gamma} / Helium-I. We do not see any significant evidence for deviations of G2's dynamical evolution, due to hydrodynamical interactions with the hot gas around the black hole, from a ballistic orbit of an initially compact cloud with moderate velocity dispersion. The constant luminosity and the increasingly stretched appearance of the head of G2 in the position-velocity plane, without a central peak, is not consistent with several proposed models with continuous gas release from an initially bound zone around a faint star on the same orbit as G2.

We investigate the origin and fate of the recently discovered gas cloud G2 close to the Galactic Center. Our hydrodynamical simulations focussing on the dynamical evolution of the cloud in combination with currently available observations favour two scenarios: a Compact Cloud which started around the year 1995 and an extended Spherical Shell of gas, with an apocenter distance within the disk(s) of young stars. The former is able to explain the detected signal of G2 in the position-velocity-diagram of the year 2008.5 and 2011.5 data. The latter can account for both, G2's signal as well as the fainter extended tail-like structure G2t seen at larger distances to the black hole and smaller velocities. From these first idealised simulations we expect a rise of the current activity of Sgr A* shortly after the closest approach and a constant feeding through a nozzle-like structure over a long period. The near future evolution of the cloud will be a sensitive probe of the conditions of the gas distribution in the milli-parsec environment of the massive black hole in the Galactic Center and will also give us invaluable information of the feeding of black holes and the activation of the central source.