We perform 3D general relativistic smoothed particle hydrodynamics (GRSPH) simulations of tidal disruption events involving 1 $M_\odot$ stars and $10^6 M_\odot$ rotating supermassive black holes. We consider stars on initially elliptical orbits both in, and inclined to, the black hole equatorial plane. We confirm that stream-stream collisions caused by relativistic apsidal precession rapidly circularise the disrupted material into a disc. For inclined trajectories we find that nodal precession induced by the black hole spin (i.e. Lense-Thirring precession) inhibits stream-stream collisions only in the first orbit, merely causing a short delay in forming a disc, which is inclined to the black hole equatorial plane. We also investigate the effect of radiative cooling on the remnant disc structure. We find that with no cooling a thick, extended, slowly precessing torus is formed, with a radial extent of 5 au (for orbits with a high penetration factor). Radiatively efficient cooling produces a narrow, rapidly precessing ring close to pericentre. We plot the energy dissipation rate, which tracks the pancake shock, stream-stream collisions and viscosity. We compare this to the effective luminosity due to accretion onto the black hole. We find energy dissipation rates of $\sim10^{45}$ erg s$^{-1}$ for stars disrupted at the tidal radius, and up to $\sim10^{47}$ erg s$^{-1}$ for deep encounters.

We present a method for general relativistic smoothed particle hydrodynamics (GRSPH), based on an entropy-conservative form of the general relativistic hydrodynamic equations for a perfect fluid. We aim to replace approximate treatments of general relativity in current SPH simulations of tidal disruption events and accretion discs. We develop an improved shock capturing formulation that distinguishes between shock viscosity and conductivity in relativity. We also describe a new Hamiltonian time integration algorithm for relativistic orbital dynamics and GRSPH. Our method correctly captures both Einstein and spin-induced precession around black holes. We benchmark our scheme in 1D and 3D against mildly and ultra relativistic shock tubes, exact solutions for epicyclic and vertical oscillation frequencies, and Bondi accretion. We assume fixed background metrics (Minkowski, Schwarzschild and Kerr in Cartesian Boyer-Lindquist coordinates) but the method lays the foundation for future direct coupling with numerical relativity.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

We test the hypothesis that the initial mass function (IMF) is determined by the density probability distribution function (PDF) produced by supersonic turbulence. We compare 14 simulations of star cluster formation in 50 solar mass molecular cloud cores where the initial turbulence contains either purely solenoidal or purely compressive modes, in each case resolving fragmentation to the opacity limit to determine the resultant IMF. We find statistically indistinguishable IMFs between the two sets of calculations, despite a factor of two difference in the star formation rate and in the standard deviation of $\log(\rho)$. This suggests that the density PDF, while determining the star formation rate, is not the primary driver of the IMF.

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parameterised approaches and find good agreement with semi-empirical measures for the "explodability" of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.