Source author record

Daniel J. Price

Daniel J. Price appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.EP astro-ph.IM astro-ph.GA astro-ph.HE astro-ph.CO physics.comp-ph physics.flu-dyn gr-qc physics.plasm-ph

Catalog footprint

What is connected

65works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Converged simulations of the nozzle shock in tidal disruption events

When debris from a star that experienced a tidal disruption events (TDE) after passing too close to a massive black hole returns to pericenter on the second passage, it is compressed, leading to the formation of nozzle shocks (in the orbital plane) and pancake shocks (perpendicular to the orbital plane). Resolving these shocks is a long-standing problem in the hydrodynamic simulations of parabolic TDEs. Excessive numerical energy dissipation or heating unrealistically expands the stream. In this Letter, we apply adaptive particle refinement to our 3D general relativistic smoothed particle simulations to locally increase the resolution near the pericenter. We achieve resolutions equivalent to $6.55\times10^{11}$ particles, allowing us to converge on the true energy dissipation. We conclude that only $4\times10^{-5}$ of the orbital energy is dissipated in nozzle shocks for a Sun-like star tidally disrupted by a $10^6$ solar-mass black hole, therefore the nozzle shocks are unlikely to be important in the evolution of TDEs.

preprint2024arXiv

Probing initial distributions of orbital eccentricity and disc misalignment via polar discs

In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the distributions of initial eccentricity and mutual inclination in multiple stellar systems at birth.

preprint2022arXiv

Accretion rates in hierarchical triple systems with discs

Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical configurations. Despite its astrophysical relevance, differential accretion in hierarchical systems remains to be understood. In this work, we investigate whether the accretion trends expected in binaries are valid for higher order multiples. We performed a set of 3D Smoothed Particle Hydrodynamics simulations of binaries and of hierarchical triples (HTs) embedded in an accretion disc, with the code Phantom. We identify for the first time accretion trends in HTs and their deviations compared to binaries. These deviations, due to the interaction of the small binary with the infalling material from the circum-triple disc, can be described with a semi-analytical prescription. Generally, the smaller binary of a HT accretes more mass than a single star of the same mass as the smaller binary. We found that in a HT, if the small binary is heavier than the third body, the standard differential accretion scenario (whereby the secondary accretes more of the mass) is hampered. Reciprocally, if the small binary is lighter than the third body, the standard differential accretion scenario is enhanced. The peculiar differential accretion mechanism we find in HTs is expected to affect their mass ratio distribution.

preprint2022arXiv

Common envelopes in massive stars II: The distinct roles of hydrogen and helium recombination

The role of recombination during a common-envelope event has been long debated. Many studies have argued that much of hydrogen recombination energy, which is radiated in relatively cool and optically-thin layers, might not thermalise in the envelope. On the other hand, helium recombination contains 30% of the total recombination energy, and occurs much deeper in the stellar envelope. We investigate the distinct roles played by hydrogen and helium recombination in a common-envelope interaction experienced by a 12 solar mass red supergiant donor. We perform adiabatic, 3D hydrodynamical simulations that (i) include hydrogen, helium, and molecular hydrogen recombination, (ii) include hydrogen and helium recombination, (iii) include only helium recombination, and (iv) do not include recombination energy. By comparing these simulations, we find that the addition of helium recombination energy alone ejects 30% more envelope mass, and leads to a 16% larger post-plunge-in separation. Under the adiabatic assumption, adding hydrogen recombination energy increases the amount of ejected mass by a further 40%, possibly unbinding the entire envelope, but does not affect the post-plunge separation. Most of the ejecta becomes unbound at relatively high (>70%) degrees of hydrogen ionisation, where the hydrogen recombination energy is likely to expand the envelope instead of being radiated away.

preprint2022arXiv

Common envelopes in massive stars: Towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes

We perform 3D hydrodynamical simulations of a common-envelope event involving a 12 solar mass red supergiant donor. Massive stars are expected to be qualitatively different from low-mass stars as their envelopes have significant support from radiation pressure, which increases both the final separation and amount of mass ejected through the common-envelope interaction. We perform adiabatic simulations that include radiation energy through the equation of state, which results in ejecting 60 per cent more mass (up to two thirds of the total envelope mass becoming unbound, or more) and yield a 10 per cent larger final separation compared to simulations that assume an ideal gas. When also including recombination energy, we find that at least three quarters of the envelope, and possibly the entire envelope, may be unbound. The final separation further increases by almost 20 per cent. The additional amount of ejected material is mainly due to energy injected from helium recombination. Hydrogen recombination plays a comparatively small role, as it mainly occurs in gas that has already become unbound. We conclude that the internal energy of the envelope can be a significant energy source for ejecting the common envelope, but ultimately radiation transport and convection need to be included.

preprint2022arXiv

External or internal companion exciting the spiral arms in CQ Tau?

We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$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$σ$ 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.

preprint2022arXiv

Kinematic evidence for an embedded planet in the IM Lupi disc

We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disc can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2-3 times the mass of Jupiter can reproduce non-Keplerian velocity perturbations, or `kinks', in the 12CO J=2-1 channel maps. Such a planet can also explain the spiral arms seen in 1.25mm dust continuum emission and 1.6 micron scattered light images. We show that the wake of the planet can be traced in the observed peak velocity map, which appears to closely follow the morphology expected from our simulations and from analytic models of planet-disc interaction.

preprint2022arXiv

Mapping the Planetary Wake in HD 163296 with Kinematics

We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD~163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk major axis confirm that the wake induces predominantly radial motion, as predicted by models of planet-disk interaction. These results provide the first direct confirmation of planet wakes generated by Lindblad resonances. Mapping the wake provides a constraint on the disk aspect ratio, which is required to measure the mass of the planet.

preprint2022arXiv

On the origin of magnetic fields in stars II: The effect of numerical resolution

Are the kG-strength magnetic fields observed in young stars a fossil field left over from their formation or are they generated by a dynamo? Our previous numerical study concluded that magnetic fields must originate by a dynamo process. Here, we continue that investigation by performing even higher numerical resolution calculations of the gravitational collapse of a 1~M$_\odot$ rotating, magnetised molecular cloud core through the first and second collapse phases until stellar densities are reached. Each model includes Ohmic resistivity, ambipolar diffusion, and the Hall effect. We test six numerical resolutions, using between $10^5$ and $3\times10^7$ particles to model the cloud. At all but the lowest resolutions, magnetic walls form in the outer parts of the first hydrostatic core, with the maximum magnetic field strength located within the wall rather than at the centre of the core. At high resolution, this magnetic wall is disrupted by the Hall effect, producing a magnetic field with a spiral-shaped distribution of intensity. As the second collapse occurs, this field is dragged inward and grows in strength, with the maximum field strength increasing with resolution. As the second core forms, the maximum field strength exceeds 1~kG in our highest resolution simulations, and the stellar core field strength exceeds this threshold at the highest resolution. Our resolution study suggests that kG-strength magnetic fields may be implanted in low-mass stars during their formation, and may persist over long timescales given that the diffusion timescale for the magnetic field exceeds the age of the Universe.

preprint2022arXiv

The Origin of the Doppler-flip in HD 100546: a large scale spiral arm generated by an inner binary companion

Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner spiral and the previously reported CO emission Doppler-flip, which has been interpreted as the signature of an embedded protoplanet. Comparisons with hydrodynamical models indicate that this Doppler-flip is instead the kinematic counterpart of the spiral, which is likely generated by an inner companion inside the disk cavity.

preprint2020arXiv

A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics

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.

preprint2020arXiv

Are the spiral arms in the MWC 758 protoplanetary disc driven by a companion inside the cavity?

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.

preprint2020arXiv

Binary-induced spiral arms inside the disc cavity of AB Aurigae

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° 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° and 90°. 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.

preprint2020arXiv

Discovery of a low-mass companion embedded in the disk of the young massive star MWC 297 with VLT/SPHERE

We report the discovery of a low-mass stellar companion around the young Herbig Be star MWC 297. We performed multi-epoch high-contrast imaging in the near infrared (NIR) with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. The companion is found at projected separation of 244.7$\pm$13.2 au and a position angle of 176.4$\pm$0.1 deg. The large separation supports formation via gravitational instability. From the spectrum, we estimate a mass of 0.1-0.5 M$_{\odot}$, the range conveying uncertainties in the extinction of the companion and in evolutionary models at young ages. The orbit coincides with a gap in the dust disk inferred from the Spectral Energy Distribution (SED). The young age ($\lesssim$ 1 Myr) and mass ratio with the central star ($\sim 0.01$) makes the companion comparable to PDS 70~b, suggesting a relation between formation scenarios and disk dynamics.

preprint2020arXiv

Is the gap in the DS Tau disc hiding a planet?

Recent mm-wavelength surveys performed with the Atacama Large Millimeter Array (ALMA) have revealed protoplanetary discs characterized by rings and gaps. A possible explanation for the origin of such rings is the tidal interaction with an unseen planetary companion. The protoplanetary disc around DS Tau shows a wide gap in the ALMA observation at 1.3 mm. We construct a hydrodynamical model for the dust continuum observed by ALMA assuming the observed gap is carved by a planet between one and five Jupiter masses. We fit the shape of the radial intensity profile along the disc major axis varying the planet mass, the dust disc mass, and the evolution time of the system. The best fitting model is obtained for a planet with $M_{\rm p}=3.5\,M_{\rm Jup}$ and a disc with $M_{\rm dust}= 9.6\cdot10^{-5}\,M_{\odot}$. Starting from this result, we also compute the expected signature of the planet in the gas kinematics, as traced by CO emission. We find that such a signature (in the form of a "kink" in the channel maps) could be observed by ALMA with a velocity resolution between $0.2-0.5\,\rm{kms}^{-1}$ and a beam size between 30 and 50 mas.

preprint2020arXiv

Visualizing the Kinematics of Planet Formation

A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the 'Visualizing the Kinematics of Planet Formation' workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.

preprint2016arXiv

A comparison between grid and particle methods on the small-scale dynamo in magnetised supersonic turbulence

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.

preprint2016arXiv

Can non-ideal magnetohydrodynamics solve the magnetic braking catastrophe?

We investigate whether or not the low ionisation fractions in molecular cloud cores can solve the `magnetic braking catastrophe', where magnetic fields prevent the formation of circumstellar discs around young stars. We perform three-dimensional smoothed particle non-ideal magnetohydrodynamics (MHD) simulations of the gravitational collapse of one solar mass molecular cloud cores, incorporating the effects of ambipolar diffusion, Ohmic resistivity and the Hall effect alongside a self-consistent calculation of the ionisation chemistry assuming 0.1 micron grains. When including only ambipolar diffusion or Ohmic resistivity, discs do not form in the presence of strong magnetic fields, similar to the cases using ideal MHD. With the Hall effect included, disc formation depends on the direction of the magnetic field with respect to the rotation vector of the gas cloud. When the vectors are aligned, strong magnetic braking occurs and no disc is formed. When the vectors are anti-aligned, a disc with radius of 13AU can form even in strong magnetic when all three non-ideal terms are present, and a disc of 38 AU can form when only the Hall effect is present; in both cases, a counter-rotating envelope forms around the first hydrostatic core. For weaker, anti-aligned fields, the Hall effect produces massive discs comparable to those produced in the absence of magnetic fields, suggesting that planet formation via gravitational instability may depend on the sign of the magnetic field in the precursor molecular cloud core.

preprint2016arXiv

Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds

We present an updated constrained hyperbolic/parabolic divergence cleaning algorithm for smoothed particle magnetohydrodynamics (SPMHD) that remains conservative with wave cleaning speeds which vary in space and time. This is accomplished by evolving the quantity $ψ/ c_h$ instead of $ψ$. Doing so allows each particle to carry an individual wave cleaning speed, $c_h$, that can evolve in time without needing an explicit prescription for how it should evolve, preventing circumstances which we demonstrate could lead to runaway energy growth related to variable wave cleaning speeds. This modification requires only a minor adjustment to the cleaning equations and is trivial to adopt in existing codes. Finally, we demonstrate that our constrained hyperbolic/parabolic divergence cleaning algorithm, run for a large number of iterations, can reduce the divergence of the field to an arbitrarily small value, achieving $\nabla \cdot B=0$ to machine precision.

preprint2016arXiv

Gas squeezing during the merger of a supermassive black hole binary

We study accretion rates during the gravitational wave-driven merger of a binary supermassive black hole embedded in an accretion disc, formed by gas driven to the centre of the galaxy. We use 3D simulations performed with PHANTOM, a Smoothed Particle Hydrodynamics code. Contrary to previous investigations, we show that there is evidence of a "squeezing phenomenon", caused by the compression of the inner disc gas when the secondary black hole spirals towards the primary. This causes an increase in the accretion rates that always exceed the Eddington rate. We have studied the main features of the phenomenon for a mass ratio $q = 10^{-3}$ between the black holes, including the effects of numerical resolution, the secondary accretion radius and the disc thickness. With our disc model with a low aspect ratio, we show that the mass expelled from the orbit of the secondary is negligible ($< 5\%$ of the initial disc mass), different to the findings of previous 2D simulations with thicker discs. The increase in the accretion rates in the last stages of the merger leads to an increase in luminosity, making it possible to detect an electromagnetic precursor of the gravitational wave signal emitted by the coalescence.

preprint2016arXiv

Grand challenges in protoplanetary disc modelling

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

preprint2016arXiv

On dust entrainment in photoevaporative winds

We investigate dust entrainment by photoevaporative winds in protoplanetary discs using dusty smoothed particle hydrodrodynamics (SPH). We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in EUV driven winds. The maximum grain size is set by dust settling in the disc rather than aerodynamic drag in the wind. More generally, there is a linear relationship between the base flow density and the maximum entrainable grain size in the wind. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between $10$--$20\,$AU. The peak dust density for each grain size occurs at the inner edge of its own entrainment region.

preprint2016arXiv

On the origin of horseshoes in transitional discs

We investigate whether the rings, lopsided features and horseshoes observed at millimetre wavelengths in transitional discs can be explained by the dynamics of gas and dust at the edge of the cavity in circumbinary discs. We use 3D dusty smoothed particle hydrodynamics calculations to show that binaries with mass ratio $q \gtrsim 0.04$ drive eccentricity in the central cavity, naturally leading to a crescent-like feature in the gas density, which is accentuated in the mm dust grain population with intensity contrasts in mm-continuum emission of 10 or higher. We perform mock observations to demonstrate that these features closely match those observed by ALMA, suggesting that the origin of rings, dust horseshoes and other non-axisymmetric structures in transition discs can be explained by the presence of massive companions.

preprint2016arXiv

Suppression of the accretion rate in thin discs around binary black holes

We present three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations investigating the dependence of the accretion rate on the disc thickness around an equal-mass, circular black hole binary system. We find that for thick/hot discs, with $H/R\gtrsim 0.1$, the binary torque does not prevent the gas from penetrating the cavity formed in the disc by the binary (in line with previous investigations). The situation drastically changes for thinner discs, in this case the mass accretion rate is suppressed, such that only a fraction (linearly dependent on $H/R$) of the available gas is able to flow within the cavity and accrete on to the binary. Extrapolating this result to the cold and thin accretion discs expected around supermassive black hole binary systems implies that this kind of systems accretes less material than predicted so far, with consequences not only for the electromagnetic and gravitational waves emissions during the late inspiral phase but also for the recoil speed of the black hole formed after binary coalescence, thus influencing also the evolutionary path both of the binary and of the host galaxy. Our results, being scale-free, are also applicable to equal mass, circular binaries of stellar mass black holes, such as the progenitor of the recently discovered gravitational wave source GW150914.

preprint2016arXiv

The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics

Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve $\nabla \cdot {\bf B}=0$ to machine precision, albeit at significant computational expense.

preprint2016arXiv

Two mechanisms for dust gap opening in protoplanetary discs

We identify two distinct physical mechanisms for dust gap opening by embedded planets in protoplanetary discs based on the symmetry of the drag-induced motion around the planet: I) A mechanism where low mass planets, that do not disturb the gas, open gaps in dust by tidal torques assisted by drag in the inner disc, but resisted by drag in the outer disc; and II) The usual, drag assisted, mechanism where higher mass planets create pressure maxima in the gas disc which the drag torque then acts to evacuate further in the dust. The first mechanism produces gaps in dust but not gas, while the second produces partial or total gas gaps which are deeper in the dust phase. Dust gaps do not necessarily indicate gas gaps.

preprint2015arXiv

A fast and explicit algorithm for simulating the dynamics of small dust grains with smoothed particle hydrodynamics

We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This "diffusion approximation for dust" is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.

preprint2015arXiv

Apsidal precession, disc breaking and viscosity in warped discs

We demonstrate the importance of general relativistic apsidal precession in warped black hole accretion discs by comparing three - dimensional smoothed particle hydrodynamic simulations in which this effect is first neglected, and then included. If apsidal precession is neglected, we confirm the results of an earlier magnetohydrodynamic simulation which made this assumption, showing that at least in this case the $α$ viscosity model produces very similar results to those of simulations where angular momentum transport is due to the magnetorotational instability. Including apsidal precession significantly changes the predicted disc evolution. For moderately inclined discs thick enough that tilt is transported by bending waves, we find a disc tilt which is nonzero at the inner disc edge and oscillates with radius, consistent with published analytic results. For larger inclinations we find disc breaking.

preprint2015arXiv

Disc formation from tidal disruptions of stars on eccentric orbits by Schwarzschild black holes

The potential of tidal disruption of stars to probe otherwise quiescent supermassive black holes cannot be exploited, if their dynamics is not fully understood. So far, the observational appearance of these events has been derived from analytical extrapolations of the debris dynamical properties just after disruption. By means of hydrodynamical simulations, we investigate the subsequent fallback of the stream of debris towards the black hole for stars already bound to the black hole on eccentric orbits. We demonstrate that the debris circularize due to relativistic apsidal precession which causes the stream to self-cross. The circularization timescale varies between 1 and 10 times the period of the star, being shorter for more eccentric and/or deeper encounters. This self-crossing leads to the formation of shocks that increase the thermal energy of the debris. If this thermal energy is efficiently radiated away, the debris settle in a narrow ring at the circularization radius with shock-induced luminosities of $\sim 10-10^3 \, L_{\rm Edd}$. If instead cooling is impeded, the debris form an extended torus located between the circularization radius and the semi-major axis of the star with heating rates $\sim 1-10^2 \, L_{\rm Edd}$. Extrapolating our results to parabolic orbits, we infer that circularization would occur via the same mechanism in $\sim 1$ period of the most bound debris for deeply penetrating encounters to $\sim 10$ for grazing ones. We also anticipate the same effect of the cooling efficiency on the structure of the disc with associated luminosities of $\sim 1-10 \, L_{\rm Edd}$ and heating rates of $\sim 0.1-1 \, L_{\rm Edd}$. In the latter case of inefficient cooling, we deduce a viscous timescale generally shorter than the circularization timescale. This suggests an accretion rate through the disc tracing the fallback rate, if viscosity starts acting promptly.

preprint2015arXiv

Evidence for enhanced persistent emission during sub-Eddington thermonuclear bursts

The standard approach for time-resolved X-ray spectral analysis of thermonuclear bursts involves subtraction of the pre-burst emission as background. This approach implicitly assumes that the persistent flux remains constant throughout the burst. We reanalyzed 332 photospheric radius expansion bursts observed from 40 sources by the Rossi X-ray Timing Explorer, introducing a multiplicative factor $f_a$ to the persistent emission contribution in our spectral fits. We found that for the majority of spectra the best-fit value of $f_a$ is significantly greater than 1, suggesting that the persistent emission typically increases during a burst. Elevated $f_a$ values were not found solely during the radius expansion interval of the burst, but were also measured in the cooling tail. The modified model results in a lower average value of the $χ^2$ fit statistic, indicating superior spectral fits, but not yet to the level of formal statistical consistency for all the spectra. We interpret the elevated $f_a$ values as an increase of the mass accretion rate onto the neutron star during the burst, likely arising from the effects of Poynting-Robertson drag on the disk material. We measured an inverse correlation of $f_a$ with the persistent flux, consistent with theoretical models of the disc response. We suggest that this modified approach may provide more accurate burst spectral parameters, as well as offering a probe of the accretion disk structure.

preprint2015arXiv

Post-periapsis pancakes: sustenance for self-gravity in tidal disruption events

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

preprint2015arXiv

Stable smoothed particle magnetohydrodynamics in very steep density gradients

The equations of smoothed particle magnetohydrodynamics (SPMHD), even with the various corrections to instabilities so far proposed, have been observed to be unstable when a very steep density gradient is necessarily combined with a variable smoothing length formalism. Here we consider in more detail the modifications made to the SPMHD equations in LBP2015 that resolve this instability by replacing the smoothing length in the induction and anisotropic force equations with an average smoothing length term. We then explore the choice of average used and compare the effects on a test `cylinder-in-a-box' problem and the collapse of a magnetised molecular cloud core. We find that, aside from some benign numerical effects at low resolutions for the quadratic mean, the formalism is robust as to the choice of average but that in complicated models it is essential to apply the average to both equations; in particular, all four averages considered exhibit similar conservation properties. This improved formalism allows for arbitrarily small sink particles and field geometries to be explored, vastly expanding the range of astronomical problems that can be modeled using SPMHD.

preprint2015arXiv

Tearing up a misaligned accretion disc with a binary companion

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.

preprint2015arXiv

Two phase mixtures in SPH - A new approach

We present a new approach to simulating mixtures of gas and dust in smoothed particle hydrodynamics (SPH). We show how the two-fluid equations can be rewritten to describe a single-fluid 'mixture' moving with the barycentric velocity, with each particle carrying a dust fraction. We show how this formulation can be implemented in SPH while preserving the conservation properties (i.e. conservation of mass of each phase, momentum and energy). We also show that the method solves two key issues with the two fluid approach: it avoids over-damping of the mixture when the drag is strong and prevents a problem with dust particles becoming trapped below the resolution of the gas. We also show how the general one-fluid formulation can be simplified in the limit of strong drag (i.e. small grains) to the usual SPH equations plus a diffusion equation for the evolution of the dust fraction that can be evolved explicitly and does not require any implicit timestepping. We present tests of the simplified formulation showing that it is accurate in the small grain/strong drag limit. We discuss some of the issues we have had to solve while developing this method and finally present a preliminary application to dust settling in protoplanetary discs.

preprint2014arXiv

Ambipolar diffusion in smoothed particle magnetohydrodynamics

In partially ionised plasmas, the magnetic field can become decoupled from the neutral gas and diffuse through it in a process known as ambipolar diffusion. Although ambipolar diffusion has been implemented in several grid codes, we here provide an implementation in smoothed particle magnetohydrodynamics (SPMHD). We use the strong coupling approximation in which the ion density is negligible, allowing a single fluid approach. The equations are derived to conserve energy, and to provide a positive definite contribution to the entropy. We test the implementation in both a simple 1D SPMHD code and the fully 3D code PHANTOM. The wave damping test yields agreement within 0.03-2 per cent of the analytical result, depending on the value of the collisional coupling constant. The oblique C-shocks test yields results that typically agree within 4 per cent of the semi-analytical result. Our algorithm is therefore suitable for exploring the effect ambipolar diffusion has on physical processes, such as the formation of stars from molecular clouds.

preprint2014arXiv

Dust and gas mixtures with multiple grain species - a one-fluid approach

We derive the single-fluid evolution equations describing a mixture made of a gas phase and an arbitrary number of dust phases, generalising the approach developed in Laibe & Price (2014a). A generalisation for continuous dust distributions as well as analytic approximations for strong drag regimes are also provided. This formalism lays the foundation for numeri- cal simulations of dust populations in a wide range of astrophysical systems while avoiding limitations associated with a multiple-fluid treatment. The usefulness of the formalism is illustrated on a series of analytical problems, namely the dustybox, dustyshock and dustywave problems as well as the radial drift of grains and the streaming instability in protoplanetary discs. We find physical effects specific to the presence of several dust phases and multiple drag timescales, including non-monotonic evolution of the differential velocity between phases and increased efficiency of the linear growth of the streaming instability. Interestingly, it is found that under certain conditions, large grains can migrate outwards in protoplanetary discs. This may explain the presence of small pebbles at several hundreds of astronomical units from their central star.

preprint2014arXiv

Dusty gas with one fluid

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.

preprint2014arXiv

Dusty gas with one fluid in smoothed particle hydrodynamics

In a companion paper we have shown how the equations describing gas and dust as two fluids coupled by a drag term can be reformulated to describe the system as a single fluid mixture. Here we present a numerical implementation of the one-fluid dusty gas algorithm using Smoothed Particle Hydrodynamics (SPH). The algorithm preserves the conservation properties of the SPH formalism. In particular, the total gas and dust mass, momentum, angular momentum and energy are all exactly conserved. Shock viscosity and conductivity terms are generalised to handle the two-phase mixture accordingly. The algorithm is benchmarked against a comprehensive suit of problems: dustybox, dustywave, dustyshock and dustyoscill, each of them addressing different properties of the method. We compare the performance of the one-fluid algorithm to the standard two-fluid approach. The one-fluid algorithm is found to solve both of the fundamental limitations of the two- fluid algorithm: it is no longer possible to concentrate dust below the resolution of the gas (they have the same resolution by definition), and the spatial resolution criterion h < csts, required in two-fluid codes to avoid over-damping of kinetic energy, is unnecessary. Implicit time stepping is straightforward. As a result, the algorithm is up to ten billion times more efficient for 3D simulations of small grains. Additional benefits include the use of half as many particles, a single kernel and fewer SPH interpolations. The only limitation is that it does not capture multi-streaming of dust in the limit of zero coupling, suggesting that in this case a hybrid approach may be required.

preprint2014arXiv

Giant Outbursts in Be/X-ray Binaries

Be/X-ray binary systems exhibit both periodic (Type I) X-ray outbursts and giant (Type II) outbursts, whose origin has remained elusive. We suggest that Type II X-ray outbursts occur when a highly misaligned decretion disk around the Be star becomes eccentric, allowing the compact object companion to capture a large amount of material at periastron. Using 3D smoothed particle hydrodynamics simulations we model the long term evolution of a representative Be/X-ray binary system. We find that periodic (Type I) X-ray outbursts occur when the neutron star is close to periastron for all disk inclinations. Type II outbursts occur for large misalignment angles and are associated with eccentricity growth that occurs on a timescale of about 10 orbital periods. Mass capture from the eccentric decretion disk results in an accretion disk around the neutron star whose estimated viscous time is long enough to explain the extended duration of Type II outbursts. Previous studies suggested that the outbursts are caused by a warped disk but our results suggest that this is not sufficient, the disk must be both highly misaligned and eccentric to initiate a Type II accretion event.

preprint2014arXiv

The Kozai-Lidov Mechanism in Hydrodynamical Disks

We use three dimensional hydrodynamical simulations to show that a highly misaligned accretion disk around one component of a binary system can exhibit global Kozai-Lidov cycles, where the inclination and eccentricity of the disk are interchanged periodically. This has important implications for accreting systems on all scales, for example, the formation of planets and satellites in circumstellar and circumplanetary disks, outbursts in X-ray binary systems and accretion on to supermassive black holes.

preprint2014arXiv

The morphology of the Milky Way - I. Reconstructing CO maps from simulations in fixed potentials

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.

preprint2013arXiv

A switch to reduce resistivity in smoothed particle magnetohydrodynamics

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 α_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én 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.

preprint2013arXiv

Collapse of a molecular cloud core to stellar densities: stellar core and outflow formation in radiation magnetohydrodynamics simulations

We have performed smoothed particle radiation magnetohydrodynamics (SPRMHD) simulations of the collapse of rotating, magnetised molecular cloud cores to form protostars. The calculations follow the formation and evolution of the first hydrostatic core, the collapse to form a stellar core, the launching of outflows from both the first hydrostatic core and stellar cores, and the breakout of the stellar outflow from the remnant of the first core. We investigate the roles of magnetic fields and thermal feedback on the outflow launching process, finding that both magnetic and thermal forces contribute to the launching of the stellar outflow. We also follow the stellar cores until they grow to masses of up to 20 Jupiter-masses, and determine their properties. We find that at this early stage, before fusion begins, the stellar cores have radii of $\approx 3$ R$_\odot$ with radial entropy profiles that increase outward (i.e. are convectively stable) and minimum entropies per baryon of $s/k_{\rm B} \approx 14$ in their interiors. The structure of the stellar cores is found to be insensitive to variations in the initial magnetic field strength. With reasonably strong initial magnetic fields, accretion on to the stellar cores occurs through inspiralling magnetised pseudo-discs with negligible radiative losses, as opposed to first cores which effectively radiate away the energy liberated in the accretion shocks at their surfaces. We find that magnetic field strengths of >10 kG can be implanted in stellar cores at birth.

preprint2013arXiv

Evidence for accretion rate change during type I X-ray bursts

preprint2013arXiv

Smoothed Particle Magnetohydrodynamics Simulations of Protostellar Jets and Turbulent Dynamos

We presents results from Smoothed Particle Magnetohydrodynamics simulations of collapsing molecular cloud cores, and dynamo amplification of the magnetic field in the presence of Mach 10 magnetised turbulence. Our star formation simulations have produced, for the first time ever, highly collimated magnetised protostellar jets from the first hydrostatic core phase. Up to 40% of the initial core mass may be ejected through this outflow. The primary difficulty in performing these simulations is maintaining the divergence free constraint of the magnetic field, and to address this issue, we have developed a new divergence cleaning method which has allowed us to stably follow the evolution of these protostellar jets for long periods. The simulations performed of supersonic MHD turbulence are able to exponentially amplify magnetic energy by up to 10 orders of magnitude via turbulent dynamo. To reduce numerical dissipation, a new shock detection algorithm is utilised which is able to track magnetic shocks throughout a large range of magnetic field strengths.

preprint2013arXiv

Wave-like warp propagation in circumbinary discs I. Analytic theory and numerical simulations

In this paper we analyse the propagation of warps in protostellar circumbinary discs. We use these systems as a test environment in which to study warp propagation in the bending-wave regime, with the addition of an external torque due to the binary gravitational potential. In particular, we want to test the linear regime, for which an analytic theory has been developed. In order to do so, we first compute analytically the steady state shape of an inviscid disc subject to the binary torques. The steady state tilt is a monotonically increasing function of radius. In the absence of viscosity, the disc does not present any twist. Then, we compare the time-dependent evolution of the warped disc calculated via the known linearised equations both with the analytic solutions and with full 3D numerical simulations, which have been performed with the PHANTOM SPH code using 2 million particles. We find a good agreement both in the tilt and in the phase evolution for small inclinations, even at very low viscosities. Moreover, we have verified that the linearised equations are able to reproduce the diffusive behaviour when α > H/R, where α is the disc viscosity parameter. Finally, we have used the 3D simulations to explore the non-linear regime. We observe a strongly non-linear behaviour, which leads to the breaking of the disc. Then, the inner disc starts precessing with its own precessional frequency. This behaviour has already been observed with numerical simulations in accretion discs around spinning black holes. The evolution of circumstellar accretion discs strongly depends on the warp evolution. Therefore the issue explored in this paper could be of fundamental importance in order to understand the evolution of accretion discs in crowded environments, when the gravitational interaction with other stars is highly likely, and in multiple systems.

preprint2012arXiv

Collimated jets from the first core

We have performed Smoothed Particle Magnetohydrodynamics (SPMHD) simulations demonstrating the production of collimated jets during collapse of 1 solar mass molecular cloud cores to form the `first hydrostatic core' in low mass star formation. Recently a number of candidate first core objects have been observed, including L1448 IRS2E, L1451-mm and Per Bolo 58, although it is not yet clear that these are first hydrostatic cores. Recent observations of Per Bolo 58 in particular appear to show collimated, bipolar outflows which are inconsistent with previous theoretical expectations. We show that low mass first cores can indeed produce tightly collimated jets (opening angles <~ 10 degrees) with speeds of ~2-7 km/s, consistent with some of the observed candidates. We have also demonstrated, for the first time, that such phenomena can be successfully captured in SPMHD simulations.

preprint2012arXiv

Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics

We present a constrained formulation of Dedner et al's hyperbolic/parabolic divergence cleaning scheme for enforcing the \nabla\dot B = 0 constraint in Smoothed Particle Magnetohydrodynamics (SPMHD) simulations. The constraint we impose is that energy removed must either be conserved or dissipated, such that the scheme is guaranteed to decrease the overall magnetic energy. This is shown to require use of conjugate numerical operators for evaluating \nabla\dot B and \nablaψ in the SPMHD cleaning equations. The resulting scheme is shown to be stable at density jumps and free boundaries, in contrast to an earlier implementation by Price & Monaghan (2005). Optimal values of the damping parameter are found to be σ = 0.2-0.3 in 2D and σ = 0.8-1.2 in 3D. With these parameters, our constrained Hamiltonian formulation is found to provide an effective means of enforcing the divergence constraint in SPMHD, typically maintaining average values of h |\nabla\dot B| / |B| to 0.1-1%, up to an order of magnitude better than artificial resistivity without the associated dissipation in the physical field. Furthermore, when applied to realistic, 3D simulations we find an improvement of up to two orders of magnitude in momentum conservation with a corresponding improvement in numerical stability at essentially zero additional computational expense.

preprint2012arXiv

Hyperbolic Divergence Cleaning for SPH

We present SPH formulations of Dedner et al's hyperbolic/parabolic divergence cleaning scheme for magnetic and velocity fields. Our implementation preserves the conservation properties of SPH which is important for stability. This is achieved by deriving an energy term for the Psi field, and imposing energy conservation on the cleaning subsystem of equations. This necessitates use of conjugate operators for divB and gradPsi in the numerical equations. For both the magnetic and velocity fields, the average divergence error in the system is reduced by an order of magnitude with our cleaning algorithm. Divergence errors in SPMHD are maintained to < 1%, even for realistic 3D applications with a corresponding gain in numerical stability. Density errors for an oscillating elliptic water drop using weakly compressible SPH are reduced by a factor of two.

preprint2012arXiv

On the accumulation of planetesimals near disc gaps created by protoplanets

We have performed three-dimensional two-fluid (gas-dust) hydrodynamical models of circumstellar discs with embedded protoplanets (3 - 333 M\oplu) and small solid bodies (radii 10cm to 10m). We find that high mass planets (\gtrsim 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. (2009), 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.

preprint2012arXiv

Response of a circumbinary accretion disc to black hole mass loss

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.

preprint2011arXiv

Dusty gas with SPH - I. Algorithm and test suite

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 generalised SPH density estimate in multi-fluid systems, the consistent treatment of variable smoothing length terms, finite particle size, time step 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 standardised, 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 Δ< cs ts 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.

preprint2011arXiv

Dusty gas with SPH - II. Implicit timestepping and astrophysical drag regimes

In a companion paper (Laibe & Price 2011b), 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 momentum 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 the timestep and the drag stopping time is 1 < r < 1000.

preprint2011arXiv

DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures

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.

preprint2011arXiv

Resolving high Reynolds numbers in SPH simulations of subsonic turbulence

Accounting for the Reynolds number is critical in numerical simulations of turbulence, particularly for subsonic flow. For Smoothed Particle Hydrodynamics (SPH) with constant artificial viscosity coefficient alpha, it is shown that the effective Reynolds number in the absence of explicit physical viscosity terms scales linearly with the Mach number - compared to mesh schemes, where the effective Reynolds number is largely independent of the flow velocity. As a result, SPH simulations with alpha=1 will have low Reynolds numbers in the subsonic regime compared to mesh codes, which may be insufficient to resolve turbulent flow. This explains the failure of Bauer and Springel (2011, arXiv:1109.4413v1) to find agreement between the moving-mesh code AREPO and the GADGET SPH code on simulations of driven, subsonic (v ~ 0.3 c_s) turbulence appropriate to the intergalactic/intracluster medium, where it was alleged that SPH is somehow fundamentally incapable of producing a Kolmogorov-like turbulent cascade. We show that turbulent flow with a Kolmogorov spectrum can be easily recovered by employing standard methods for reducing alpha away from shocks.

preprint2011arXiv

Smoothed Particle Hydrodynamics: Things I wish my mother taught me

I discuss the key features of Smoothed Particle Hydrodynamics (SPH) as a numerical method - in particular the key differences between SPH and more standard grid based approaches - that are important to the practitioner. These include the exact treatment of advection, the absence of intrinsic dissipation, exact conservation and more subtle properties that arise from its Hamiltonian formulation such as the existence of a minimum energy state for the particles. The implications of each of these are discussed, showing how they can be both advantages and disadvantages.

preprint2010arXiv

A comparison between grid and particle methods on the statistics of driven, supersonic, isothermal turbulence

We compare the statistics of driven, supersonic turbulence at high Mach number using FLASH a widely used Eulerian grid-based code and PHANTOM, a Lagrangian 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 rho^1/3 v as suggested by Kritsuk et al. and a log-normal PDF with a width that scales with 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 subgrid 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 max. densities at a resolution of 128^3 particles that were similar those 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 rho^1/3 v. For SPH simulations of high Mach number turbulence we find it important to use sufficient non-linear beta-viscosity to prevent particle interpenetration in shocks (we require beta=4 instead of the default beta=2).

preprint2010arXiv

A method for reconstructing the PDF of a 3D turbulent density field from 2D observations

We introduce a method for calculating the probability density function (PDF) of a turbulent density field in three dimensions using only information contained in the projected two-dimensional column density field. We test the method by applying it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds. To a good approximation, the PDF of log(normalised column density) is a compressed, shifted version of the PDF of log(normalised density). The degree of compression can be determined observationally from the column density power spectrum, under the assumption of statistical isotropy of the turbulence.

preprint2010arXiv

Magnetic fields and radiative feedback in the star formation process

Star formation is a complex process involving the interplay of many physical effects, including gravity, turbulent gas dynamics, magnetic fields and radiation. Our understanding of the process has improved substantially in recent years, primarily as a result of our increased ability to incorporate the relevant physics in numerical calculations of the star formation process. In this contribution we present an overview of our recent studies of star cluster formation in turbulent, magnetised clouds using self-gravitating radiation-magnetohydrodynamics calculations (Price and Bate 2008, 2009). Our incorporation of magnetic fields and radiative transfer into the Smoothed Particle Hydrodynamics method are discussed. We highlight how magnetic fields and radiative heating of the gas around newborn stars can solve several of the key puzzles in star formation, including an explanation for why star formation is such a slow and inefficient process. However, the presence of magnetic fields at observed strengths in collapsing protostellar cores also leads to problems on smaller scales, including a difficulty in forming protostellar discs and binary stars (Price and Bate 2007, Hennebelle and Teyssier 2008), which suggests that our understanding of the role of magnetic fields in star formation is not yet complete.

preprint2010arXiv

Magnetic fields and Turbulence in Star Formation using Smoothed Particle Hydrodynamics

Firstly, we give a historical overview of attempts to incorporate magnetic fields into the Smoothed Particle Hydrodynamics method by solving the equations of Magnetohydrodynamics (MHD), leading an honest assessment of the current state-of-the-art in terms of the limitations to performing realistic calculations of the star formation process. Secondly, we discuss the results of a recent comparison we have performed on simulations of driven, supersonic turbulence with SPH and Eulerian techniques. Finally we present some new results on the relationship between the density variance and the Mach number in supersonic turbulent flows, finding sigma^2_{ln rho} = ln (1 + b^2 M^2) with b=0.33 up to Mach~20, consistent with other numerical results at lower Mach number (Lemaster and Stone 2008) but inconsistent with observational constraints on sigma_rho and M in Taurus and IC5146.

preprint2010arXiv

Smoothed Particle Hydrodynamics and Magnetohydrodynamics

This paper presents an overview and introduction to Smoothed Particle Hydrodynamics and Magnetohydrodynamics in theory and in practice. Firstly, we give a basic grounding in the fundamentals of SPH, showing how the equations of motion and energy can be self-consistently derived from the density estimate. We then show how to interpret these equations using the basic SPH interpolation formulae and highlight the subtle difference in approach between SPH and other particle methods. In doing so, we also critique several `urban myths' regarding SPH, in particular the idea that one can simply increase the `neighbour number' more slowly than the total number of particles in order to obtain convergence. We also discuss the origin of numerical instabilities such as the pairing and tensile instabilities. Finally, we give practical advice on how to resolve three of the main issues with SPMHD: removing the tensile instability, formulating dissipative terms for MHD shocks and enforcing the divergence constraint on the particles, and we give the current status of developments in this area. Accompanying the paper is the first public release of the NDSPMHD SPH code, a 1, 2 and 3 dimensional code designed as a testbed for SPH/SPMHD algorithms that can be used to test many of the ideas and used to run all of the numerical examples contained in the paper.

preprint2010arXiv

Smoothed Particle Hydrodynamics: Turbulence and MHD

In this paper we discuss recent applications of the Smoothed Particle Hydrodynamics (SPH) method to the simulation of supersonic turbulence in the interstellar medium, as well as giving an update on recent algorithmic developments in solving the equations of magnetohydrodynamics (MHD) in SPH. Using high resolution calculations (up to 134 million particles), we find excellent agreement with grid-based results on a range of measures including the power spectrum slope in both the velocity field and the density-weighted velocity rho^(1/3) v, the latter showing a Kolmogorov-like k^-5/3 scaling as proposed by Kritsuk et al. (2007). We also find good agreement on the statistics of the Probability Distribution Function (PDF) and structure functions, independently confirming the scaling found by Schmidt, Federrath & Klessen (2008). On Smoothed Particle Magnetohydrodynamics (SPMHD) we have recently wasted a great deal of time and effort investigating the vector potential as an alternative to the Euler potentials formulation, in the end concluding that using the vector potential has even more severe problems than the standard (B-field based) SPMHD approach.

preprint2010arXiv

The density variance -- Mach number relation in supersonic, isothermal turbulence

We examine the relation between the density variance and the mean-square Mach number in supersonic, isothermal turbulence, assumed in several recent analytic models of the star formation process. From a series of calculations of supersonic, hydrodynamic turbulence driven using purely solenoidal Fourier modes, we find that the `standard' relationship between the variance in the log of density and the Mach number squared, i.e., sigma^2_(ln rho/rhobar)=ln (1+b^2 M^2), with b = 1/3 is a good fit to the numerical results in the supersonic regime up to at least Mach 20, similar to previous determinations at lower Mach numbers. While direct measurements of the variance in linear density are found to be severely underestimated by finite resolution effects, it is possible to infer the linear density variance via the assumption of log-normality in the Probability Distribution Function. The inferred relationship with Mach number, consistent with sigma_(rho/rhobar) ~ b M with b=1/3, is, however, significantly shallower than observational determinations of the relationship in the Taurus Molecular Cloud and IC5146 (both consistent with b~ 0.5), implying that additional physics such as gravity is important in these clouds and/or that turbulent driving in the ISM contains a significant compressive component. Magnetic fields are not found to change this picture significantly, in general reducing the measured variances and thus worsening the discrepancy with observations.

preprint2009arXiv

Inefficient star formation: The combined effects of magnetic fields and radiative feedback

We investigate the effects of magnetic fields and radiative protostellar feedback on the star formation process using self-gravitating radiation magnetohydrodynamical calculations. We present results from a series of calculations of the collapse of 50 solar mass molecular clouds with various magnetic field strengths and with and without radiative transfer. We find that both magnetic fields and radiation have a dramatic impact on star formation, though the two effects are in many ways complementary. Magnetic fields primarily provide support on large scales to low density gas, whereas radiation is found to strongly suppress small-scale fragmentation by increasing the temperature in the high-density material near the protostars. With strong magnetic fields and radiative feedback the net result is an inefficient star formation process with a star formation rate of ~< 10% per free-fall time that approaches the observed rate, although we have only been able to follow the calculations for ~1/3 of a free-fall time beyond the onset of star formation.

preprint2009arXiv

Smoothed Particle Magnetohydrodynamics IV - Using the Vector Potential

In this paper we investigate the use of the vector potential as a means of maintaining the divergence constraint in the numerical solution of the equations of Magnetohydrodynamics (MHD) using the Smoothed Particle Hydrodynamics (SPH) method. We derive a self-consistent formulation of the equations of motion using a variational principle that is constrained by the numerical formulation of both the induction equation and the curl operator used to obtain the magnetic field, which guarantees exact and simultaneous conservation of momentum, energy and entropy in the numerical scheme. This leads to a novel formulation of the MHD force term, unique to the vector potential, which differs from previous formulations. We also demonstrate how dissipative terms can be correctly formulated for the vector potential such that the contribution to the entropy is positive definite and the total energy is conserved. On a standard suite of numerical tests in one, two and three dimensions we find firstly that the consistent formulation of the vector potential equations is unstable to the well-known SPH tensile instability, even more so than in the standard Smoothed Particle Magnetohydrodynamics (SPMHD) formulation where the magnetic field is evolved directly. Furthermore we find that, whilst a hybrid approach based on the vector potential evolution equation coupled with a standard force term gives good results for one and two dimensional problems (where dAz/dt = 0), such an approach suffers from numerical instability in three dimensions related to the unconstrained evolution of vector potential components. We conclude that use of the vector potential is not a viable approach for Smoothed Particle Magnetohydrodynamics.

Daniel J. Price

What is connected

Connect this record

See the researcher in context

Building this map preview

65 published item(s)

Converged simulations of the nozzle shock in tidal disruption events

Probing initial distributions of orbital eccentricity and disc misalignment via polar discs

Accretion rates in hierarchical triple systems with discs

Common envelopes in massive stars II: The distinct roles of hydrogen and helium recombination

Common envelopes in massive stars: Towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes

External or internal companion exciting the spiral arms in CQ Tau?

Kinematic evidence for an embedded planet in the IM Lupi disc

Mapping the Planetary Wake in HD 163296 with Kinematics

On the origin of magnetic fields in stars II: The effect of numerical resolution

The Origin of the Doppler-flip in HD 100546: a large scale spiral arm generated by an inner binary companion

A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics

Are the spiral arms in the MWC 758 protoplanetary disc driven by a companion inside the cavity?

Binary-induced spiral arms inside the disc cavity of AB Aurigae

Discovery of a low-mass companion embedded in the disk of the young massive star MWC 297 with VLT/SPHERE

Is the gap in the DS Tau disc hiding a planet?

Visualizing the Kinematics of Planet Formation

A comparison between grid and particle methods on the small-scale dynamo in magnetised supersonic turbulence

Can non-ideal magnetohydrodynamics solve the magnetic braking catastrophe?

Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds

Gas squeezing during the merger of a supermassive black hole binary

Grand challenges in protoplanetary disc modelling

On dust entrainment in photoevaporative winds

On the origin of horseshoes in transitional discs

Suppression of the accretion rate in thin discs around binary black holes

The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics

Two mechanisms for dust gap opening in protoplanetary discs

A fast and explicit algorithm for simulating the dynamics of small dust grains with smoothed particle hydrodynamics

Apsidal precession, disc breaking and viscosity in warped discs

Disc formation from tidal disruptions of stars on eccentric orbits by Schwarzschild black holes

Evidence for enhanced persistent emission during sub-Eddington thermonuclear bursts

Post-periapsis pancakes: sustenance for self-gravity in tidal disruption events

Stable smoothed particle magnetohydrodynamics in very steep density gradients

Tearing up a misaligned accretion disc with a binary companion

Two phase mixtures in SPH - A new approach

Ambipolar diffusion in smoothed particle magnetohydrodynamics

Dust and gas mixtures with multiple grain species - a one-fluid approach

Dusty gas with one fluid

Dusty gas with one fluid in smoothed particle hydrodynamics

Giant Outbursts in Be/X-ray Binaries

The Kozai-Lidov Mechanism in Hydrodynamical Disks

The morphology of the Milky Way - I. Reconstructing CO maps from simulations in fixed potentials

A switch to reduce resistivity in smoothed particle magnetohydrodynamics

Collapse of a molecular cloud core to stellar densities: stellar core and outflow formation in radiation magnetohydrodynamics simulations

Evidence for accretion rate change during type I X-ray bursts

Smoothed Particle Magnetohydrodynamics Simulations of Protostellar Jets and Turbulent Dynamos

Wave-like warp propagation in circumbinary discs I. Analytic theory and numerical simulations

Collimated jets from the first core

Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics

Hyperbolic Divergence Cleaning for SPH

On the accumulation of planetesimals near disc gaps created by protoplanets

Response of a circumbinary accretion disc to black hole mass loss

Dusty gas with SPH - I. Algorithm and test suite

Dusty gas with SPH - II. Implicit timestepping and astrophysical drag regimes

DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures

Resolving high Reynolds numbers in SPH simulations of subsonic turbulence

Smoothed Particle Hydrodynamics: Things I wish my mother taught me

A comparison between grid and particle methods on the statistics of driven, supersonic, isothermal turbulence

A method for reconstructing the PDF of a 3D turbulent density field from 2D observations

Magnetic fields and radiative feedback in the star formation process

Magnetic fields and Turbulence in Star Formation using Smoothed Particle Hydrodynamics

Smoothed Particle Hydrodynamics and Magnetohydrodynamics

Smoothed Particle Hydrodynamics: Turbulence and MHD

The density variance -- Mach number relation in supersonic, isothermal turbulence

Inefficient star formation: The combined effects of magnetic fields and radiative feedback

Smoothed Particle Magnetohydrodynamics IV - Using the Vector Potential