The Galactic centre (GC) is a crowded environment: observations have revealed the presence of (molecular, atomic and ionized) gas, of a cusp of late-type stars, and of ~100 early-type stars, about half of which lying in one or possibly two discs. In this paper, we study the perturbations exerted on a thin stellar disc (with outer radius ~0.4 pc) by a molecular cloud that falls towards the GC and is disrupted by the supermassive black hole (SMBH). The initial conditions for the stellar disc were drawn from the results of previous simulations of molecular cloud infall and disruption in the SMBH potential. We find that most of the gas from the disrupted molecular cloud settles into a dense and irregular disc surrounding the SMBH. If the gas disc and the stellar disc are slightly misaligned (~5-20 deg), the precession of the stellar orbits induced by the gas disc significantly increases the inclinations of the stellar orbits (by a factor of ~3-5 in 1.5 Myr) with respect to the normal vector to the disc. Furthermore, the distribution of orbit inclinations becomes significantly broader. These results might be the clue to explain the broad distribution of observed inclinations of the early-type stars with respect to the normal vector of the main disc. We discuss the implications for the possibility that fresh gas was accreted by the GC after the formation of the disc(s) of early-type stars.

We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries, and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outward, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year timescale, which is smaller than timescales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse timescales can be shorter than inward migration timescales, suggesting that clumps could survive tidal disruption and deliver a proto-gas giant to distances of even a few AU from the central star.

The formation of the massive young stars surrounding SgrA* is still an open question. In this paper, we simulate the infall of a turbulent molecular cloud towards the Galactic Center (GC). We adopt two different cloud masses (4.3x10^4 and 1.3x10^5 solar masses). We run five simulations: the gas is assumed to be isothermal in four runs, whereas radiative cooling is included in the fifth run. In all the simulations, the molecular cloud is tidally disrupted, spirals towards the GC, and forms a small, dense and eccentric disk around SgrA*. With high resolution simulations, we follow the fragmentation of the gaseous disk. Star candidates form in a ring at ~0.1-0.4 pc from the super-massive black hole (SMBH) and have moderately eccentric orbits (~0.2-0.4), in good agreement with the observations. The mass function of star candidates is top-heavy only if the local gas temperature is high (>~100 K) during the star formation and if the parent cloud is sufficiently massive (>~10^5 solar masses). Thus, this study indicates that the infall of a massive molecular cloud is a viable scenario for the formation of massive stars around SgrA*, provided that the gas temperature is kept sufficiently high (>~100 K).

We present a novel implementation of Smoothed Particle Hydrodynamics (SPHS) that uses the spatial derivative of the velocity divergence as a higher order dissipation switch. Our switch -- which is second order accurate -- detects flow convergence before it occurs. If particle trajectories are going to cross, we switch on the usual SPH artificial viscosity, as well as conservative dissipation in all advected fluid quantities (for example, the entropy). The viscosity and dissipation terms (that are numerical errors) are designed to ensure that all fluid quantities remain single-valued as particles approach one another, to respect conservation laws, and to vanish on a given physical scale as the resolution is increased. SPHS alleviates a number of known problems with `classic' SPH, successfully resolving mixing, and recovering numerical convergence with increasing resolution. An additional key advantage is that -- treating the particle mass similarly to the entropy -- we are able to use multimass particles, giving significantly improved control over the refinement strategy. We present a wide range of code tests including the Sod shock tube, Sedov-Taylor blast wave, Kelvin-Helmholtz Instability, the `blob test', and some convergence tests. Our method performs well on all tests, giving good agreement with analytic expectations.

We present high-resolution 3D smoothed particle hydrodynamics simulations of the formation and evolution of protostellar discs in a turbulent molecular cloud. Using a piecewise polytropic equation of state, we perform two sets of simulations. In both cases we find that isolated systems undergo a fundamentally different evolution than members of binary or multiple systems. When formed, isolated systems must accrete mass and increase their specific angular momentum, leading to the formation of massive, extended discs, which undergo strong gravitational instabilities and are susceptible to disc fragmentation. Fragments with initial masses of 5.5 M_jup, 7.4 M_jup and 12 M_jup are produced in our simulations. In binaries and small clusters, we observe that due to competition for material from the parent core, members do not accrete significant amounts of high specific angular momentum gas relative to isolated systems. We find that discs in multiple systems are strongly self-gravitating but that they are stable against fragmentation due to disc truncation and mass profile steeping by tides, accretion of high specific angular momentum gas by other members, and angular momentum being redirected into members' orbits. In general, we expect disc fragmentation to be less likely in clusters and to be more a feature of isolated systems.

Standard formulations of smoothed particle hydrodynamics (SPH) are unable to resolve mixing at fluid boundaries. We use an error and stability analysis of the generalised SPH equations of motion to prove that this is due to two distinct problems. The first is a leading order error in the momentum equation. This should decrease with increasing neighbour number, but does not because numerical instabilities cause the kernel to be irregularly sampled. We identify two important instabilities: the clumping instability and the banding instability, and we show that both are cured by a suitable choice of kernel. The second problem is the local mixing instability (LMI). This occurs as particles attempt to mix on the kernel scale, but are unable to due to entropy conservation. The result is a pressure discontinuity at boundaries that pushes fluids of different entropy apart. We cure the LMI by using a weighted density estimate that ensures that pressures are single valued throughout the flow. This also gives a better volume estimate for the particles, reducing errors in the continuity and momentum equations. We demonstrate mixing in our new Optimised Smoothed Particle Hydrodynamics (OSPH) scheme using a Kelvin Helmholtz instability (KHI) test with density contrast 1:2, and the 'blob test' - a 1:10 density ratio gas sphere in a wind tunnel - finding excellent agreement between OSPH and Eulerian codes.

We explore the initial conditions for fragments in the extended regions ($r\gtrsim 50$ AU) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H$_2$ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.

The formation of brown dwarfs (BDs) due to the fragmentation of proto-stellar disks undergoing pairwise encounters was investigated. High resolution allowed the use of realistic initial disk models where both the vertical structure and the local Jeans mass were resolved. The results show that objects with masses ranging from giant planets to low mass stars can form during such encounters from initially stable disks. The parameter space of initial spin-orbit orientations and the azimuthal angles for each disk was explored. An upper limit on the initial Toomre Q value of ~2 was found for fragmentation to occur. Depending on the initial configuration, shocks, tidal-tail structures and mass inflows were responsible for the condensation of disk gas. Retrograde disks were generally more likely to fragment. When the interaction timescale was significantly shorter than the disks' dynamical timescales, the proto-stellar disks tended to be truncated without forming objects. The newly-formed objects had masses ranging from 0.9 to 127 Jupiter masses, with the majority in the BD regime. They often resided in star-BD multiples and in some cases also formed hierarchical orbiting systems. Most of them had large angular momenta and highly flattened, disk-like shapes. The objects had radii ranging from 0.1 to 10 AU. The disk gas was assumed to be locally isothermal, appropriate for the short cooling times in extended proto-stellar disks, but not for condensed objects. An additional case with explicit cooling that reduced to zero for optically thick gas was simulated to test the extremes of cooling effectiveness and it was still possible to form objects in this case. Detailed radiative transfer is expected to lengthen the internal evolution timescale for these objects, but not to alter our basic results.

The formation of the massive young stars surrounding SgrA* is still an open question. In this paper, we simulate the infall of an isothermal, turbulent molecular cloud towards the Galactic Centre (GC). As it spirals towards the GC, the molecular cloud forms a small and dense disc around SgrA*. Efficient star formation (SF) is expected to take place in such a dense disc. We model SF by means of sink particles. At ~6x10^5 yr, ~6000 solar masses of stars have formed, and are confined within a thin disc with inner and outer radius of 0.06 and 0.5 pc, respectively. Thus, this preliminary study shows that the infall of a molecular cloud is a viable scenario for the formation of massive stars around SgrA*. Further studies with more realistic radiation physics and SF will be required to better constrain this intriguing scenario.