Search

Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity., Accepted for publication in MNRAS, 20 pages

We perform 3D hydrodynamics simulations of disc-disc stellar flybys with on-the-fly Monte Carlo radiative transfer. We show that pre-existing circumstellar discs around both stars result in fast rising ($\sim$yrs) outbursts lasting 2-5 times longer than for a star-disc flyby. The perturber always goes into outburst ($\dot{M}>10^{-5}~{\rm M_{\odot}~ yr^{-1}}$). Whereas we find that the primary goes into a decades long outburst only when the flyby is retrograde to the circumprimary disc rotation. High accretion rates during the outburst are triggered by angular momentum cancellation in misaligned material generated by the encounter. A large fraction of accreted material is alien., 11 pages, 9 figures, accepted for publication in MNRAS

Planet-disc interactions build up local pressure maxima that may halt the radial drift of protoplanetary dust, and pile it up in rings and crescents. ALMA observations of the HD135344B disc revealed two rings in the thermal continuum stemming from ~mm-sized dust. At higher frequencies the inner ring is brighter relative to the outer ring, which is also shaped as a crescent rather than a full ring. In near-IR scattered light images, the disc is modulated by a 2-armed grand-design spiral originating inside the ALMA inner ring. Such structures may be induced by a massive companion evacuating the central cavity, and by a giant planet in the gap separating both rings, that channels the accretion of small dust and gas through its filamentary wakes while stopping the larger dust from crossing the gap. Here we present ALMA observations in the J=(2-1)CO isotopologue lines and in the adjacent continuum, with up to 12km baselines. Angular resolutions of 0.03" reveal the tentative detection of a filament connecting both rings, and which coincides with a local discontinuity in the pitch angle of the IR spiral, proposed previously as the location of the protoplanet driving this spiral. Line diagnostics suggest that turbulence, or superposed velocity components, is particularly strong in the spirals. The 12CO(2-1) 3-D rotation curve points at stellocentric accretion at radii within the inner dust ring, with a radial velocity of up to ~6%+-0.5% Keplerian, which corresponds to an excessively large accretion rate of ~2E-6M_sun/yr if all of the CO layer follows the 12CO(2-1) kinematics. This suggests that only the surface layers of the disc are undergoing accretion, and that the line broadening is due to superposed laminar flows., accepted to MNRAS

We present ALMA observations of the $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1 transitions and the 1.3\,mm continuum emission for the circumbinary disc around HD 142527, at an angular resolution of $\approx$\,0\farcs3. We observe multiple spiral structures in intensity, velocity and velocity dispersion for the $^{12}$CO and $^{13}$CO gas tracers. A newly detected $^{12}$CO spiral originates from the dust horseshoe, and is rotating at super-Keplerian velocity or vertically ascending, whilst the inter-spiral gas is rotating at sub-Keplerian velocities. This new spiral possibly connects to a previously identified spiral, thus spanning > 360$^\circ$. A spatial offset of ~30 au is observed between the $^{12}$CO and $^{13}$CO spirals, to which we hypothesize that the gas layers are propagating at different speeds (``surfing'') due to a non-zero vertical temperature gradient. Leveraging the varying optical depths between the CO isotopologues, we reconstruct temperature and column density maps of the outer disc. Gas surface density peaks at r\,$\approx$\,180\,au, coincident with the peak of continuum emission. Here the dust grains have a Stokes number of $\approx$\,1, confirming radial and azimuthal trapping in the horseshoe. We measure a cavity radius at half-maximum surface density of $\approx$\,100\,au, and a cavity eccentricity between 0.3 and 0.45.

We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $$20 au). We find that the 1.3 mm radial profiles of these objects are in good agreement with those produced by numerical simulations of dust evolution and planet-disc interactions, which predict the accumulation of mm-sized grains at the edges of planet-induced cavities. Our long-baseline observations resulted in the largest sample of discs observed at $\sim$3-5 au resolution in any given star-forming region (15 objects when combined with Ophiuchus objects in the DSHARP Large Program) and allow for a demographic study of the brightest $\sim5\%$ of the discs in Ophiuchus (i.e. the most likely formation sites of giant planets in the cloud). We use this unique sample to propose an evolutionary sequence and discuss a scenario in which the substructures observed in massive protoplanetary discs are mainly the result of planet formation and dust evolution. If this scenario is correct, the detailed study of disc substructures might provide a window to investigate a population of planets that remains mostly undetectable by other techniques., Comment: 21 pages, 10 figures. Appendix with 3 additional figures. Version 2, identical to previous one, but now accepted for publication (in MNRAS)

How does the cavity size in circumbinary discs depend on disc and binary properties? We investigate by simulating disc cavities carved by binary companions using smoothed particle hydrodynamics (SPH). We find that a cavity is quickly opened on the dynamical time, while the cavity size is set on the viscous time. In agreement with previous findings, we find long term cavity sizes of 2-5 times the binary semi-major axis, increasing with eccentricity and decreasing with disc aspect ratio. When considering binaries inclined with respect to the disc we find three regimes: i) discs that evolve towards a coplanar orbit have a large cavity, slightly smaller than that of an initially coplanar disc; ii) discs that evolve towards a polar orbit by breaking have a small cavity, equal in size to that of an initially polar disc; iii) discs that evolve towards a polar orbit via warping have an intermediate-sized cavity. We find typical gas depletions inside the cavity of $\gtrsim 2$ orders of magnitude in surface density., Comment: 12 pages, 13 figures, accepted to MNRAS

During the common envelope binary interaction, the expanding layers of the gaseous common envelope recombine and the resulting recombination energy has been suggested as a contributing factor to the ejection of the envelope. In this paper we perform a comparative study between simulations with and without the inclusion of recombination energy. We use two distinct setups, comprising 0.88-M$_{\odot}$ and 1.8-M$_{\odot}$ giants, that have been studied before and can serve as benchmarks. In so doing we conclude that (i) the final orbital separation is not affected by the choice of equation of state. In other words, simulations that unbind but a small fraction of the envelope result in similar final separations to those that, thanks to recombination energy, unbind a far larger fractions. (ii) The adoption of a tabulated equation of state results in a much greater fraction of unbound envelope and we demonstrate the cause of this to be the release of recombination energy. (iii) The fraction of hydrogen recombination energy that is allowed to do work should be about half of that which our adiabatic simulations use. (iv) However, for the heavier star simulation we conclude that it is helium and not hydrogen recombination energy that unbinds the gas and we determine that all helium recombination energy is thermalised in the envelope and does work. (v)~The outer regions of the expanding common envelope are likely to see the formation of dust. This dust would promote additional unbinding and shaping of the ejected envelope into axisymmetric morphologies., Comment: 18 pages, 11 figures. We welcome any comments or questions

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa)., 11 pages, 5 figures, accepted for publication in MNRAS

We present the implementation of a dust growth and fragmentation module in the public Smoothed Particle Hydrodynamics (SPH) code PHANTOM. This module is made available for public use with this paper. The coagulation model considers locally monodisperse dust size distributions around single values that are carried by the SPH particles. Along with the presentation of the model, implementation and tests, we showcase growth and fragmentation in a few typical circumstellar disc simulations and revisit previous results. The module is also interfaced with the radiative transfer code MCFOST, which facilitates the comparison between simulations and ALMA observations by generating synthetic maps. Circumstellar disc simulations with growth and fragmentation reproduce the `self-induced dust trap' mechanism first proposed by Gonzalez et al., which supports its existence. Synthetic images of discs featuring this mechanism suggest it would be detectable by ALMA as a bright axisymmetric ring at several tens of au from the star. With this paper, our aim is to provide a public tool to be able to study and explore dust growth in a variety of applications related to planet formation., 22 pages, 27 figures, accepted for publication in MNRAS

IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$\alpha$ imaging and Keck astrometric measures of the binary to develop a coherent dynamical model of this system. The system features an azimuthal asymmetry detected at the western edge of the cavity in Band~7 observations and a wiggling outflow. Dust emission in ALMA Band 4 from the proximity of the individual stars suggests the presence of marginally resolved circumstellar discs. We estimate the binary orbital parameters from the measured arc of the orbit from Keck and ALMA astrometry. We further constrain these estimates using considerations from binary-disc interaction theory. We finally perform three SPH gas + dust simulations based on the theoretical constraints; we post-process the hydrodynamic output using radiative transfer Monte Carlo methods and directly compare the models with observations. Our results suggest that a highly eccentric $e\sim 0.5\textrm{--}0.7$ equal mass binary, with a semi-major axis of $\sim 55$ au, and small/moderate orbital plane vs. circumbinary disc inclination $\theta\lesssim 30^\circ$ provides a good match with observations. A dust mass of $\sim 1.5\times 10^{-4} {\rm M_\odot}$ best reproduces the flux in Band 7 continuum observations. Synthetic CO line emission maps qualitatively capture both the emission from the central region and the non-Keplerian nature of the gas motion in the binary proximity., Comment: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

Spiral arms in protoplanetary discs are thought to be linked to the presence of companions. We test the hypothesis that the double spiral arm morphology observed in the transition disc MWC 758 can be generated by an $\approx 10$ M$_{\rm Jup}$ companion on an eccentric orbit internal to the spiral arms. Previous studies on MWC 758 have assumed an external companion. We compare simulated observations from three dimensional hydrodynamics simulations of disc-companion interaction to scattered light, infrared and CO molecular line observations, taking into account observational biases. The inner companion hypothesis is found to explain the double spiral arms, as well as several additional features seen in MWC 758 -- the arc in the northwest, substructures inside the spiral arms, the cavity in CO isotopologues, and the twist in the kinematics. Testable predictions include detection of fainter spiral structure, detection of a point source south-southeast of the primary, and proper motion of the spiral arms., Comment: 13 pages, 6 figures, resubmitted to MNRAS

What cosmic ray ionisation rate is required such that a non-ideal magnetohydrodynamics (MHD) simulation of a collapsing molecular cloud will follow the same evolutionary path as an ideal MHD simulation or as a purely hydrodynamics simulation? To investigate this question, we perform three-dimensional smoothed particle non-ideal magnetohydrodynamics simulations of the gravitational collapse of rotating, one solar mass, magnetised molecular cloud cores, that include Ohmic resistivity, ambipolar diffusion, and the Hall effect. We assume a uniform grain size of $a_\text{g} = 0.1\mu$m, and our free parameter is the cosmic ray ionisation rate, $\zeta_\text{cr}$. We evolve our models, where possible, until they have produced a first hydrostatic core. Models with $\zeta_\text{cr}\gtrsim10^{-13}$ s$^{-1}$ are indistinguishable from ideal MHD models and the evolution of the model with $\zeta_\text{cr}=10^{-14}$ s$^{-1}$ matches the evolution of the ideal MHD model within one per cent when considering maximum density, magnetic energy, and maximum magnetic field strength as a function of time; these results are independent of $a_\text{g}$. Models with very low ionisation rates ($\zeta_\text{cr}\lesssim10^{-24}$ s$^{-1}$) are required to approach hydrodynamical collapse, and even lower ionisation rates may be required for larger $a_\text{g}$. Thus, it is possible to reproduce ideal MHD and purely hydrodynamical collapses using non-ideal MHD given an appropriate cosmic ray ionisation rate. However, realistic cosmic ray ionisation rates approach neither limit, thus non-ideal MHD cannot be neglected in star formation simulations., Comment: 13 pages, 15 figures, accepted for publication in MNRAS

We present results from radiation non-ideal magnetohydrodynamics (MHD) calculations that follow the collapse of rotating, magnetised, molecular cloud cores to stellar densities. These are the first such calculations to include all three non-ideal effects: ambipolar diffusion, Ohmic resistivity and the Hall effect. We employ an ionisation model in which cosmic ray ionisation dominates at low temperatures and thermal ionisation takes over at high temperatures. We explore the effects of varying the cosmic ray ionisation rate from $\zeta_\text{cr}= 10^{-10}$ to $10^{-16}$ s$^{-1}$. Models with ionisation rates $\gtrsim 10^{-12}$ s$^{-1}$ produce results that are indistinguishable from ideal MHD. Decreasing the cosmic ray ionisation rate extends the lifetime of the first hydrostatic core up to a factor of two, but the lifetimes are still substantially shorter than those obtained without magnetic fields. Outflows from the first hydrostatic core phase are launched in all models, but the outflows become broader and slower as the ionisation rate is reduced. The outflow morphology following stellar core formation is complex and strongly dependent on the cosmic ray ionisation rate. Calculations with high ionisation rates quickly produce a fast (~14km s$^{-1}$) bipolar outflow that is distinct from the first core outflow, but with the lowest ionisation rate a slower (~3-4 km s$^{-1}$) conical outflow develops gradually and seamlessly merges into the first core outflow., Comment: 22 pages, 24 figures. Accepted for publication in MNRAS

We investigate the effect of non-ideal magnetohydrodynamics (MHD) on the formation of binary stars using a suite of three-dimensional smoothed particle magnetohydrodynamics simulations of the gravitational collapse of one solar mass, rotating, perturbed molecular cloud cores. Alongside the role of Ohmic resistivity, ambipolar diffusion and the Hall effect, we also examine the effects of magnetic field strength, orientation and amplitude of the density perturbation. When modelling sub-critical cores, ideal MHD models do not collapse whereas non-ideal MHD models collapse to form single protostars. In super-critical ideal MHD models, increasing the magnetic field strength or decreasing the initial density perturbation amplitude decreases the initial binary separation. Strong magnetic fields initially perpendicular to the rotation axis suppress the formation of binaries and yield discs with magnetic fields ~10 times stronger than if the magnetic field was initially aligned with the rotation axis. When non-ideal MHD is included, the resulting discs are larger and more massive, and the binary forms on a wider orbit. Small differences in the super-critical cores caused by non-ideal MHD effects are amplified by the binary interaction near periastron. Overall, the non-ideal effects have only a small impact on binary formation and early evolution, with the initial conditions playing the dominant role., Accepted to MNRAS. 18 pages, 20 figures, 1 appendix

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., Comment: 6 pages, 7 figures, accepted to MNRAS

We study the evolution of galactic magnetic fields using 3D smoothed particle magnetohydrodynamics (SPMHD) simulations of galaxies with an imposed spiral potential. We consider the appearance of reversals of the field, and amplification of the field. We find magnetic field reversals occur when the velocity jump across the spiral shock is above $\approx$20km s$^{-1}$, occurring where the velocity change is highest, typically at the inner Lindblad resonance (ILR) in our models. Reversals also occur at corotation, where the direction of the velocity field reverses in the co-rotating frame of a spiral arm. They occur earlier with a stronger amplitude spiral potential, and later or not at all with weaker or no spiral arms. The presence of a reversal at a radii of around 4--6 kpc in our fiducial model is consistent with a reversal identified in the Milky Way, though we caution that alternative Galaxy models could give a similar reversal. We find that relatively high resolution, a few million particles in SPMHD, is required to produce consistent behaviour of the magnetic field. Amplification of the magnetic field occurs in the models, and while some may be genuinely attributable to differential rotation or spiral arms, some may be a numerical artefact. We check our results using Athena, finding reversals but less amplification of the field, suggesting that some of the amplification of the field with SPMHD is numerical., 15 pages, 13 figures, accepted for publication in MNRAS

We perform a comparison between the smoothed particle magnetohydrodynamics (SPMHD) code, Phantom, and the Eulerian grid-based code, Flash, on the small-scale turbulent dynamo in driven, Mach 10 turbulence. We show, for the first time, that the exponential growth and saturation of an initially weak magnetic field via the small-scale dynamo can be successfully reproduced with SPMHD. The two codes agree on the behaviour of the magnetic energy spectra, the saturation level of magnetic energy, and the distribution of magnetic field strengths during the growth and saturation phases. The main difference is that the dynamo growth rate, and its dependence on resolution, differs between the codes, caused by differences in the numerical dissipation and shock capturing schemes leading to differences in the effective Prandtl number in Phantom and Flash., Accepted for publication in MNRAS, 17 pages, 16 figures

A tidal disruption event, which occurs when a star is destroyed by the gravitational field of a supermassive black hole, produces a stream of debris, the evolution of which ultimately determines the observational properties of the event. Here we show that a post-periapsis caustic -- a location where the locus of gas parcels comprising the stream would collapse into a two-dimensional surface if they evolved solely in the gravitational field of the hole -- occurs when the pericenter distance of the star is on the order of the tidal radius of the hole. It is demonstrated that this "pancake" induces significant density perturbations in the debris stream, and, for stiffer equations of state (adiabatic index $\gamma \gtrsim 5/3$), these fluctuations are sufficient to gravitationally destabilize the stream, resulting in its fragmentation into bound clumps. The results of our findings are discussed in the context of the observational properties of tidal disruption events., Comment: 18 pages, 14 figures. Accepted for publication in MNRAS

We have developed a modified form of the equations of smoothed particle magnetohydrodynamics which are stable in the presence of very steep density gradients. Using this formalism, we have performed simulations of the collapse of magnetised molecular cloud cores to form protostars and drive outflows. Our stable formalism allows for smaller sink particles (< 5 AU) than used previously and the investigation of the effect of varying the angle, {\theta}, between the initial field axis and the rotation axis. The nature of the outflows depends strongly on this angle: jet-like outflows are not produced at all when {\theta} > 30{\deg}, and a collimated outflow is not sustained when {\theta} > 10{\deg}. No substantial outflows of any kind are produced when {\theta} > 60{\deg}. This may place constraints on the geometry of the magnetic field in molecular clouds where bipolar outflows are seen., Comment: Accepted for publication in MNRAS, 13 pages, 14 figures. Animations can be found at http://www.astro.ex.ac.uk/people/blewis/research/outflows_misaligned_fields.html

Accretion discs are common in binary systems, and they are often found to be misaligned with respect to the binary orbit. The gravitational torque from a companion induces nodal precession in misaligned disc orbits. We calculate whether this precession is strong enough to overcome the internal disc torques communicating angular momentum. For typical parameters precession wins: the disc breaks into distinct planes that precess effectively independently. We run hydrodynamical simulations to check these results, and confirm that disc breaking is widespread and generally enhances accretion on to the central object. This applies in many cases of astrophysical accretion, e.g. supermassive black hole binaries and X--ray binaries., 8 pages, 6 figures, accepted for publication in MNRAS

We investigate the effect of black hole spin on warped or misaligned accretion discs - in particular i) whether or not the inner disc edge aligns with the black hole spin and ii) whether the disc can maintain a smooth transition between an aligned inner disc and a misaligned outer disc, known as the Bardeen-Petterson effect. We employ high resolution 3D smoothed particle hydrodynamics simulations of $\alpha$-discs subject to Lense-Thirring precession, focussing on the bending wave regime where the disc viscosity is smaller than the aspect ratio $\alpha \lesssim H/R$. We first address the controversy in the literature regarding possible steady-state oscillations of the tilt close to the black hole. We successfully recover such oscillations in 3D at both small and moderate inclinations ($\lesssim 15^{\circ}$), provided both Lense-Thirring and Einstein precession are present, sufficient resolution is employed, and provided the disc is not so thick so as to simply accrete misaligned. Second, we find that discs inclined by more than a few degrees in general steepen and break rather than maintain a smooth transition, again in contrast to previous findings, but only once the disc scale height is adequately resolved. Finally, we find that when the disc plane is misaligned to the black hole spin by a large angle, the disc 'tears' into discrete rings which precess effectively independently and cause rapid accretion, consistent with previous findings in the diffusive regime ($\alpha \gtrsim H/R$). Thus misalignment between the disc and the spin axis of the black hole provides a robust mechanism for growing black holes quickly, regardless of whether the disc is thick or thin., Comment: 15 pages, 18 figures, movies available at http://users.monash.edu.au/~rnealon/ or YouTube

We present an investigation into the morphological features of the Milky Way. We use smoothed particle hydrodynamics (SPH) to simulate the interstellar medium (ISM) in the Milky Way under the effect of a number of different gravitational potentials representing spiral arms and bars, assuming the Milky Way is grand design in nature. The gas is subject to ISM cooling and chemistry, enabling us to track the evolution of molecular gas. We use a 3D radiative transfer code to simulate the emission from the SPH output, allowing for the construction of synthetic longitude-velocity (l-v) emission maps as viewed from the Earth. By comparing these maps with the observed emission in CO from the Milky Way, we infer the arm/bar geometry that provides a best fit to our Galaxy. We find that it is possible to reproduce nearly all features of the l-v diagram in CO emission. There is no model, however, that satisfactorily reproduces all of the features simultaneously. Models with 2 arms cannot reproduce all the observed arm features, while 4 armed models produce too bright local emission in the inner Galaxy. Our best fitting models favour a bar pattern speed within 50-60km/s/kpc and an arm pattern speed of approximately 20km/s/kpc, with a bar orientation of approximately 45 degrees and arm pitch angle between 10-15 degrees., Comment: 25 pages, 27 Figures. Accepted for publication in MNRAS. Animations of a 2-armed and 4-armed simulation can be found at http://youtube.com/watch?v=9dAZIvo9FH8 and http://youtube.com/watch?v=mmYOMgEoaE0

Artificial resistivity is included in Smoothed Particle Magnetohydrodynamics simulations to capture shocks and discontinuities in the magnetic field. Here we present a new method for adapting the strength of the applied resistivity so that shocks are captured but the dissipation of the magnetic field away from shocks is minimised. Our scheme utilises the gradient of the magnetic field as a shock indicator, setting {\alpha}_B = h|gradB|/|B|, such that resistivity is switched on only where strong discontinuities are present. The advantage to this approach is that the resistivity parameter does not depend on the absolute field strength. The new switch is benchmarked on a series of shock tube tests demonstrating its ability to capture shocks correctly. It is compared against a previous switch proposed by Price & Monaghan (2005), showing that it leads to lower dissipation of the field, and in particular, that it succeeds at capturing shocks in the regime where the Alfv\'en speed is much less than the sound speed (i.e., when the magnetic field is very weak). It is also simpler. We also demonstrate that our recent constrained divergence cleaning algorithm has no difficulty with shock tube tests, in contrast to other implementations., Comment: 8 pages, 7 figures, accepted for publication in MNRAS

We investigate the evolution of the surface density of a circumbinary accretion disc after the mass loss induced by the merger of two supermassive black holes. We first introduce an analytical model, under the assumption of a disc composed of test particles, to derive the surface density evolution of the disc following the mass loss. The model predicts the formation of sharp density peaks in the disc; the model also allows us to compute the typical timescale for the formation of these peaks. To test and validate the model, we run numerical simulations of the process using the Smoothed Particle Hydrodynamics (SPH) code PHANTOM, taking fluid effects into account. We find good agreement in the shape and position of the peaks between the model and the simulations. In a fluid disc, however, the epicyclic oscillations induced by the mass loss can dissipate, and only some of the predicted peaks form in the simulation. To quantify how fast this dissipation proceeds, we introduce an appropriate parameter, and we show that it is effective in explaining the differences between the analytical, collisionless model and a real fluid disc.

We have performed three-dimensional two-fluid (gas–dust) hydrodynamical models of circumstellar discs with embedded protoplanets (3–333 M⊕) and small solid bodies (radii 10 cm to 10 m). We find that high-mass planets (≳ Saturn mass) open sufficiently deep gaps in the gas disc such that the density maximum at the outer edge of the gap can very efficiently trap metre-sized solid bodies. This allows the accumulation of solids at the outer edge of the gap as solids from large radii spiral inwards to the trapping region. This process of accumulation occurs fastest for those bodies that spiral inwards most rapidly, typically metre-sized boulders, whilst smaller and larger objects will not migrate sufficiently rapidly in the discs lifetime to benefit from the process. Around a Jupiter mass planet we find that bound clumps of solid material, as large as several Earth masses, may form, potentially collapsing under self-gravity to form planets or planetesimals. These results are in agreement with Lyra et al., supporting their finding that the formation of a second generation of planetesimals or of terrestrial mass planets may be triggered by the presence of a high-mass planet.

Accretion in the nuclei of active galaxies may occur chaotically. This can produce accretion discs which are counter‐rotating or strongly misaligned with respect to the spin of the central supermassive black hole (SMBH), or the axis of a close SMBH binary. Accordingly we consider the cancellation of angular momentum in accretion discs with a significant change of plane (tilt) between inner and outer parts. We estimate analytically the maximum accretion rate through such discs and compare this with the results of Smoothed Particle Hydrodynamics (SPH) simulations. These suggest that accretion rates on to supermassive black holes may be larger by factors& 100 if the disc is internally tilted in this way rather than pl anar. This offers a natural way of driving the rapid growth of supermassive black holes, and the coalescence of SMBH binaries.

In a companion paper, we have presented an algorithm for simulating two-fluid gas and dust mixtures in smoothed particle hydrodynamics (SPH). In this paper, we develop an implicit timestepping method that preserves the exact conservation of the both linear and angular momenta in the underlying SPH algorithm, but unlike previous schemes, allows the iterations to converge to arbitrary accuracy and is suited to the treatment of non-linear drag regimes. The algorithm presented in Paper I is also extended to deal with realistic astrophysical drag regimes, including both linear and non-linear Epstein and Stokes drag. The scheme is benchmarked against the test suite presented in Paper I, including (i) the analytic solutions of the dustybox problem and (ii) solutions of the dustywave, dustyshock, dustysedov and dustydisc obtained with explicit timestepping. We find that the implicit method is 1–10 times faster than the explicit temporal integration when the ratio r between the timestep and the drag stopping time is 1 ≲r≲ 1000.

We present a new algorithm for simulating two-fluid gas and dust mixtures in smoothed particle hydrodynamics (SPH), systematically addressing a number of key issues including the generalized SPH density estimate in multifluid systems, the consistent treatment of variable smoothing-length terms, finite particle size, timestep stability, thermal coupling terms and the choice of kernel and smoothing length used in the drag operator. We find that using double-hump-shaped kernels improves the accuracy of the drag interpolation by a factor of several hundred compared to the use of standard SPH bell-shaped kernels, at no additional computational expense. In order to benchmark our algorithm, we have developed a comprehensive suite of standardized, simple test problems for gas and dust mixtures: dustybox, dustywave, dustyshock, dustysedov and dustydisc, the first three of which have known analytic solutions. We present the validation of our algorithm against all of these tests. In doing so, we show that the spatial resolution criterion Δ≲csts is a necessary condition in all gas+dust codes that becomes critical at high drag (i.e. small stopping time ts) in order to correctly predict the dynamics. Implicit timestepping and the implementation of realistic astrophysical drag regimes are addressed in a companion paper.

In this paper we present the analytic solutions for two test problems involving two-fluid mixtures of dust and gas in an astrophysical context. The solutions provide a means of benchmarking numerical codes designed to simulate the non-linear dynamics of dusty gas. The first problem, dustybox, consists of two interpenetrating homogeneous fluids moving with relative velocity difference. We provide exact solutions to the full non-linear problem for a range of drag formulations appropriate to astrophysical fluids (i.e. various prescriptions for Epstein and Stokes drag in different regimes). The second problem, dustywave, consists of the propagation of linear acoustic waves in a two-fluid gas–dust mixture. We provide the analytic solution for the case when the two fluids are interacting via a linear drag term. Both test problems are simple to set up in any numerical code and can be run with periodic boundary conditions. The solutions we derive are completely general with respect to both the dust-to-gas ratio and the amplitude of the drag coefficient. A stability analysis of waves in a gas–dust system is also presented, showing that sound waves in an astrophysical dust–gas mixture are linearly stable.

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative smoothed particle hydrodynamics (SPH) formulation developed by Price and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10: 1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.

We introduce and test an expression for calculating the variance of a physical field in three dimensions using only information contained in the two-dimensional projection of the field. The method is general but assumes statistical isotropy. To test the method we apply it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds, and demonstrate that it can recover the 3D normalised density variance with ~10% accuracy if the assumption of isotropy is valid. We show that the assumption of isotropy breaks down at low sonic Mach number if the turbulence is sub-Alfvenic. Theoretical predictions suggest that the 3D density variance should increase proportionally to the square of the Mach number of the turbulence. Application of our method will allow this prediction to be tested observationally and therefore constrain a large body of analytic models of star formation that rely on it.

We examine the effect of magnetic fields on star cluster formation by performing simulations following the self-gravitating collapse of a turbulent molecular cloud to form stars in ideal MHD. The collapse of the cloud is computed for global mass-to-flux ratios of infinity, 20, 10, 5 and 3, that is using both weak and strong magnetic fields. Whilst even at very low strengths the magnetic field is able to significantly influence the star formation process, for magnetic fields with plasma beta < 1 the results are substantially different to the hydrodynamic case. In these cases we find large-scale magnetically-supported voids imprinted in the cloud structure; anisotropic turbulent motions and column density structure aligned with the magnetic field lines, both of which have recently been observed in the Taurus molecular cloud. We also find strongly suppressed accretion in the magnetised runs, leading to up to a 75% reduction in the amount of mass converted into stars over the course of the calculations and a more quiescent mode of star formation. There is also some indication that the relative formation efficiency of brown dwarfs is lower in the strongly magnetised runs due to the reduction in the importance of protostellar ejections., 16 pages, 9 figures, 8 very pretty movies, MNRAS, accepted. Version with high-res figures + movies available from http://www.astro.ex.ac.uk/people/dprice/pubs/mcluster/index.html

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D MHD simulations of galactic discs subject to a spiral potential. Recent hydrodynamic simulations have demonstrated the formation of inter-arm spurs as well as spiral arm molecular clouds provided the ISM model includes a cold HI phase. We find that the main effect of adding a magnetic field to these calculations is to inhibit the formation of structure in the disc. However, provided a cold phase is included, spurs and spiral arm clumps are still present if $\beta \gtrsim 0.1$ in the cold gas. A caveat to two phase calculations though is that by assuming a uniform initial distribution, $\beta \gtrsim 10$ in the warm gas, emphasizing that models with more consistent initial conditions and thermodynamics are required. Our simulations with only warm gas do not show such structure, irrespective of the magnetic field strength. Furthermore, we find that the introduction of a cold HI phase naturally produces the observed degree of disorder in the magnetic field, which is again absent from simulations using only warm gas. Whilst the global magnetic field follows the large scale gas flow, the magnetic field also contains a substantial random component that is produced by the velocity dispersion induced in the cold gas during the passage through a spiral shock. Without any cold gas, the magnetic field in the warm phase remains relatively well ordered apart from becoming compressed in the spiral shocks. Our results provide a natural explanation for the observed high proportions of disordered magnetic field in spiral galaxies and we thus predict that the relative strengths of the random and ordered components of the magnetic field observed in spiral galaxies will depend on the dynamics of spiral shocks.

We present a new, completely Lagrangian magnetohydrodynamics code that is based on the SPH method. The equations of self-gravitating hydrodynamics are derived self-consistently from a Lagrangian and account for variable smoothing length (``grad-h''-) terms in both the hydrodynamic and the gravitational acceleration equations. The evolution of the magnetic field is formulated in terms of so-called Euler potentials which are advected with the fluid and thus guarantee the MHD flux-freezing condition. This formulation is equivalent to a vector potential approach and therefore fulfills the $\vec{\nabla}\cdot\vec{B}=0$-constraint by construction. Extensive tests in one, two and three dimensions are presented. The tests demonstrate the excellent conservation properties of the code and show the clear superiority of the Euler potentials over earlier magnetic SPH formulations.

In this paper we describe an adaptive softening length formalism for collisionless N-body and self-gravitating Smoothed Particle Hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially variable softening lengths can be used in N-body calculations without secular increases in energy. The formalism requires the calculation of a small additional term to the gravitational force related to the gradient of the softening length. The extra term is similar in form to the usual SPH pressure force (although opposite in direction) and is therefore straightforward to implement in any SPH code at almost no extra cost. For N-body codes some additional cost is involved as the formalism requires the computation of the density via a summation over neighbouring particles using the smoothing kernel. The results of numerical tests demonstrate that, for homogeneous mass distributions, the use of adaptive softening lengths gives a softening which is always close to the `optimal' choice of fixed softening parameter, removing the need for fine-tuning. For a heterogeneous mass distribution (as may be found in any large scale N-body simulation) we find that the errors on the least-dense component are lowered by an order of magnitude compared to the use of a fixed softening length tuned to the densest component. For SPH codes our method presents a natural and elegant choice of softening formalism which makes a small improvement to both the force resolution and the total energy conservation at almost zero additional cost.

In two previous papers (Price & Monaghan 2004a,b) (papers I,II) we have described an algorithm for solving the equations of Magnetohydrodynamics (MHD) using the Smoothed Particle Hydrodynamics (SPH) method. The algorithm uses dissipative terms in order to capture shocks and has been tested on a wide range of one dimensional problems in both adiabatic and isothermal MHD. In this paper we investigate multidimensional aspects of the algorithm, refining many of the aspects considered in papers I and II and paying particular attention to the code's ability to maintain the div B = 0 constraint associated with the magnetic field. In particular we implement a hyperbolic divergence cleaning method recently proposed by Dedner et al. (2002) in combination with the consistent formulation of the MHD equations in the presence of non-zero magnetic divergence derived in papers I and II. Various projection methods for maintaining the divergence-free condition are also examined. Finally the algorithm is tested against a wide range of multidimensional problems used to test recent grid-based MHD codes. A particular finding of these tests is that in SPMHD the magnitude of the divergence error is dependent on the number of neighbours used to calculate a particle's properties and only weakly dependent on the total number of particles. Whilst many improvements could still be made to the algorithm, our results suggest that the method is ripe for application to problems of current theoretical interest, such as that of star formation., Here is the latest offering in my quest for a decent SPMHD algorithm. 26 pages, 15 figures, accepted for publication in MNRAS. Version with high res figures available from http://www.astro.ex.ac.uk/people/dprice/pubs/spmhd/spmhdpaper3.pdf

This paper describes a simple model of a star where compressibility is retained but the gravitational force is replaced by a force between pairs which is directed along their line of centres and proportional to their separation. This force was analysed by Newton. It may be considered the simplest many body force because the system reduces to a set of independent particles moving in a harmonic oscillator potential. We call such models ‘toy stars’. In this paper we simplify further and focus on the one-dimensional problem. We show that a non-linear solution exists where the velocity is a linear function of the coordinate and the density a quadratic function of the coordinates. We generalize this result to include a magnetic field. Our results provide a very useful benchmark for algorithms designed to simulate gravitating gas and we show that Smoothed Particle Hydrodynamics (SPH) simulations of toy stars with or without magnetic fields give results in good agreement with theory.

In this paper we show how a Lagrangian variational principle can be used to derive the SPMHD (smoothed particle magnetohydrodynamics) equations for ideal MHD. We also consider the effect of a variable smoothing length in the SPH kernels after which we demonstrate by numerical tests that the consistent treatment of terms relating to the gradient of the smoothing length in the SPMHD equations significantly improves the accuracy of the algorithm. Our results complement those obtained in a companion paper (Price and Monaghan 2003a, paper I) for non ideal MHD where artificial dissipative terms were included to handle shocks., Comment: 14 pages, 4 figures, accepted to MNRAS

We examine the hypothesis that there exists a simple scaling between the observed velocities of jets found in Young Stellar Objects (YSOs) and jets found in Active Galactic Nuclei (AGN). We employ a very simplified physical model of the jet acceleration process. We use time-dependent, spherically symmetric wind models in Newtonian and relativistic gravitational fields to ask whether the energy input rates required to produce the jet velocities observed in YSOs (of about 2 $\times$ the escape velocity from the central object) can also produce AGN jet velocities (Lorentz factors of $\gamma \sim$ 10). Such a scaling would be expected if there is a common production mechanism for such jets. We demonstrate that such a scaling does exist, provided that the energy input process takes place sufficiently deep in the gravitational potential well, enabling physical use to be made of the speed of light as a limiting velocity, and provided that the energy released in the accretion process is imparted to a small fraction of the available accreting material.

In this paper, we show how the two-fluid equations describing the evolution of a dust and gas mixture can be reformulated to describe a single fluid moving with the barycentric velocity of the mixture. This leads to evolution equations for the total density, momentum, the differential velocity between the dust and the gas phases and either the dust-to-gas ratio or the dust fraction. The equations are similar to the usual equations of gas dynamics, providing a convenient way to extend existing codes to simulate two-fluid mixtures without modifying the code architecture. Our approach avoids the inherent difficulties related to the standard approach where the two phases are separate and coupled via a drag term. In particular, the requirements of infinite spatial and temporal resolution as the stopping time tends to zero are no longer necessary. This means that both small and large grains can be straightforwardly treated with the same method, with no need for complicated implicit schemes. Since there is only one resolution scale the method also avoids the problem of unphysical trapping of one fluid (e.g. dust) below the resolution of the other. We also derive a simplified set of equations applicable to the case of strong drag/small grains, consisting of the standard fluid equations with a modified sound speed, plus an advection-diffusion equation for the dust-to-gas ratio. This provides a simple and fast way to evolve the mixture when the stopping time is smaller than the Courant timestep. We present a Smoothed Particle Hydrodynamics implementation in a companion paper., Comment: Accepted for publication in MNRAS (very minor revisions included)

In this paper we show how the equations of motion for the smoothed particle hydrodynamics (SPH) method may be derived from a variational principle for both non-relativistic and relativistic motion when there is no dissipation. Because the SPH density is a function of the coordinates the derivation of the equations of motion through variational principles is simpler than in the continuum case where the density is defined through the continuity equation. In particular, the derivation of the general relativistic equations is more direct and simpler than that of Fock. The symmetry properties of the Lagrangian lead immediately to the familiar additive conservation laws of linear and angular momentum and energy. In addition, we show that there is an approximately conserved quantity which, in the continuum limit, is the circulation.

We compare the statistics of driven, supersonic turbulence at a high Mach number using FLASH, a widely used Eulerian grid-based code, and PHANTOM, a Lagrangian smoothed particle hydrodynamics (SPH) code at resolutions of up to 512 3 in both grid cells and SPH particles. We find excellent agreement between codes on the basic statistical properties: a slope of k −1.95 in the velocity power spectrum for hydrodynamic, Mach 10 turbulence, evidence in both codes for a Kolmogorov-like slope of k −5/3 in the variable ρ 1/3 v as suggested by Kritsuk et al. and a lognormal probability distribution function (PDF) with a width that scales with the Mach number and proportionality constant b = 0.33–0.5 in the density variance–Mach number relation. The measured structure function slopes are not converged in either code at 512 3 elements. We find that for measuring volumetric statistics such as the power spectrum slope and structure function scaling, SPH and grid codes give roughly comparable results when the number of SPH particles is approximately equal to the number of grid cells. In particular, to accurately measure the power spectrum slope in the inertial range, in the absence of sub-grid turbulence models, requires at least 512 3 computational elements in either code. On the other hand the SPH code was found to be better at resolving dense structures, giving maximum densities at a resolution of 128 3 particles that were similar to the maximum densities resolved in the grid code at 512 3 cells, reflected also in the high density tail of the PDF. We find SPH to be more dissipative at comparable numbers of computational elements in statistics of the velocity field, but correspondingly less dissipative than the grid code in the statistics of density-weighted quantities such as ρ 1/3 v. For SPH simulations of high Mach number turbulence, we find it important to use sufficient non-linear β-viscosity in order to prevent particle interpenetration in shocks (we require βvisc = 4 instead of the widely used default value, βvisc = 2).

We present a fix to the overdamping problem found by Laibe & Price (2012) when simulating strongly coupled dust-gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective., 6 pages, 5 figures, accepted to MNRAS

We present new high-contrast images in near-infrared wavelengths ($\lambda$ = 1.04, 1.24, 1.62, 2.18 and 3.78$\mu$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M$_{\rm Jup}$ immediately outside the spiral arms at $\sim$60 au and above 2-3 M$_{\rm Jup}$ beyond 100 au to 5$\sigma$ confidence assuming hot-start models. We do, however, detect four spiral arms, a double-arc and evidence for shadows in scattered light cast by a misaligned inner disc. Our observations may be explained by an unseen close-in companion on an inclined and eccentric orbit. Such a hypothesis would also account for the disc CO cavity and disturbed kinematics., Comment: 10 figures, 14 pages. Accepted for publication in MNRAS on 25 July 2022

In this work we demonstrate that the inner spiral structure observed in AB Aurigae can be created by a binary star orbiting inside the dust cavity. We find that a companion with a mass-ratio of 0.25, semi-major axis of 40 au, eccentricity of 0.5, and inclination of 90{\deg} produces gaseous spirals closely matching the ones observed in $^{12}$CO (2-1) line emission. Based on dust dynamics in circumbinary discs (Poblete, Cuello, and Cuadra 2019), we constrain the inclination of the binary with respect to the circumbinary disc to range between 60{\deg} and 90{\deg}. We predict that the stellar companion is located roughly 0.18 arcsec from the central star towards the east-southeast, above the plane of the disc. Should this companion be detected in the near future, our model indicates that it should be moving away from the primary star at a rate of 6 mas/yr on the plane of the sky. Since our companion is inclined, we also predict that the spiral structure will appear to change with time, and not simply co-rotate with the companion., Comment: Accepted for publication in MNRAS. 11 Pages, 7 Figures

We present an online, open and comprehensive template library of gravitational waveforms produced during the tidal disruptions of stars by massive black holes, spanning a broad space of parameters. We build this library thanks to a new feature that we implement in the general relativistic version of PHANTOM, a smoothed particle hydrodynamics code for three dimensional simulations in general relativity. We first perform a series of numerical tests to show that the gravitational wave (GW) signal obtained is in excellent agreement with the one expected from theory. This benchmark is done for well studied scenarios (such as binary stellar systems). We then apply our code to calculate the GW signals from tidal disruption events (TDEs), finding that our results are consistent with the theoretical estimates obtained in previous studies for selected parameters. We illustrate interesting results from the catalogue, where we stress how the gravitational signal is affected by variations of some parameters (like black hole spin, stellar orbital eccentricity and inclination). The full catalogue is available online. It is intended to be a living catalogue., Comment: Accepted for publication on MNRAS

The protoplanetary disc HD 100453 exhibits a curious combination of spirals, shadows and a relative misalignment between the observed outer disc and inferred inner disc. This disc is accompanied by a secondary star on a bound orbit exterior to the disc. Recent observations have suggested there may be an additional low-mass companion residing within the disc inner cavity. In our companion paper the orbit of the secondary was shown to be misaligned by 61 degrees to the plane of the outer disc. Here we investigate the properties of the inner companion and the origin of the misalignment between the inner and outer disc. Using numerical simulations and synthetic observations, we show that the disc structure and kinematics are consistent with a $\lesssim$ 5 Jupiter mass planet located at 15-20au. We find that the disc evolution over around 50 binary orbits (about 10$^5$ yrs) is governed by differential precession and to a lesser extent, the Kozai-Lidov effect. In our proposed model the misalignment observed between the outer and inner disc arises naturally as a result of the misaligned outer companion driving the outer disc to precess more rapidly than the inner disc., Comment: Accepted for publication in MNRAS, 12 pages, 6 figures. Movie of simulations available at https://youtu.be/MvH8TO21W_4