Showing 1–38 of 38 results for author: Wurster, J

The Role of Drag and Gravity on Dust Concentration in a Gravitationally Unstable Disc

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru, James Wurster, Hossam Aly, Richard Alexander, Ken Rice, Richard A. Booth

Abstract: We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If w… ▽ More We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If we include both, grains that would otherwise be partially decoupled will become well coupled and trace the spirals. For the dust grains most influenced by drag (with Stokes numbers near unity), the dust disc quickly becomes gravitationally unstable and rapidly forms clumps with masses between 0.15 - 6 Earth masses. A large fraction of clumps are below the threshold where runaway gas accretion can occur. However, if dust self-gravity is neglected, the dust is unable to form clumps, despite still becoming trapped in the gas spirals. When large dust grains are unable to feel either gas gravity or drag, the dust is unable to trace the gas spirals. Hence, full physics is needed to properly simulate dust in gravitationally unstable discs. Dust trapping of large grains in spiral arms of discs stable to gas fragmentation could explain planet formation in very young discs by a population of planetesimals formed due to the combined roles of drag and gravity in the earliest stages of a disc's evolution. Furthermore, it highlights that gravitationally unstable discs are not just important for forming gas giants quickly, it can also rapidly form Earth mass bodies. △ Less

Submitted 24 January, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. Updated references

Star-forming environments in smoothed particle magnetohydrodynamics simulations II: Re-simulating isolated clumps to determine equivalence of extracted clumps and parent simulations

Authors: James Wurster, Connar Rowan

Abstract: What is the numerical reproducibility of a stellar system (including its discs) when evolving only a sub-set of (partially-evolved) smoothed particle hydrodynamics (SPH) particles? To investigate this, we modelled the evolution of 29 star forming clumps that were extracted from our previous simulations that investigated the formation and early evolution of low-mass star clusters. These clumps were… ▽ More What is the numerical reproducibility of a stellar system (including its discs) when evolving only a sub-set of (partially-evolved) smoothed particle hydrodynamics (SPH) particles? To investigate this, we modelled the evolution of 29 star forming clumps that were extracted from our previous simulations that investigated the formation and early evolution of low-mass star clusters. These clumps were evolved using a three-dimensional smoothed particle radiation magnetohydrodynamics code, where we included or excluded non-ideal magnetohydrodynamics to match the cluster simulation. While star formation proceeded as expected, we were unable to identically reproduce any of the systems present at the end of the cluster simulations. However, the final distributions of stellar mass, stellar system mass, disc mass, and disc radii were reproduced statistically; unfortunately, the distribution of average magnetic field strengths in the discs was not reproduced statistically, but this may be a result of our updated algorithms governing the evolution of the magnetic field. Therefore, given that our clumps yield stellar masses that are statistically similar to those in the original low-mass star clusters, we have demonstrated that we can statistically reproduce systems (aside from their magnetic field strength) by evolving a subset of SPH particles. Therefore, clumps such as these can be used as initial conditions to investigate the formation of isolated stars from less-contrived initial environments. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: accepted for publication in MNRAS

Star-forming environments in smoothed particle magnetohydrodynamics simulations I: Clump extraction and properties

Authors: James Wurster, Connar Rowan

Abstract: What is the nature of a star forming clump? Observations reveal these to be chaotic environments being modified and influenced by many physical processes. However, numerical simulations often define these initial star forming clumps to be idealised objects. In this paper, we define and analyse 109 star forming clumps extracted from our previous low-mass star cluster simulations. To define a clump,… ▽ More What is the nature of a star forming clump? Observations reveal these to be chaotic environments being modified and influenced by many physical processes. However, numerical simulations often define these initial star forming clumps to be idealised objects. In this paper, we define and analyse 109 star forming clumps extracted from our previous low-mass star cluster simulations. To define a clump, we identify all the gas in a simulation that ever becomes bound to or accreted onto a star, then follow the gas backwards in time until it decreases to a critical density. This gas, and its neighbouring gas, is defined as our star forming clump. Our clumps span a mass range of $0.15 \lesssim M/$M$_\odot \lesssim 10.2$, while the density range within each clump spans 2--4 orders of magnitude. The gas density distribution is not smooth, indicating that it is highly structured. The clumps are turbulent, with no coherent rotation. Independent of the initial magnetic field strength of the parent cloud, all clumps yield a similar range of field strengths. The clump magnetic field is ordered, but not reflective of the initial field geometry of the parent cloud. In general, most clump properties have a slight trend with clump mass but are independent of (or only very weakly dependent on) the properties of the parent cloud. We conclude that stars are born from a wide variety of environments and there is not a single universal star forming clump. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS

Gas and star kinematics in cloud-cloud collisions

Authors: James Wurster, Ian A. Bonnell

Abstract: We model the collision of molecular clouds to investigate the role of the initial properties on the remnants. Our clouds collide and evolve in a background medium that is approximately ten times less dense than the clouds, and we show that this relatively dense background is dynamically important for the evolution of the collision remnants. Given the motion of the clouds and the remnants through t… ▽ More We model the collision of molecular clouds to investigate the role of the initial properties on the remnants. Our clouds collide and evolve in a background medium that is approximately ten times less dense than the clouds, and we show that this relatively dense background is dynamically important for the evolution of the collision remnants. Given the motion of the clouds and the remnants through the background, we develop, implement, and introduce dynamic boundary conditions. We investigate the effect of the initial cloud mass, velocity, internal turbulence, and impact angle. The initial velocity and its velocity components have the largest affect on the remnant. This affects the spatial extent of the remnant, which affects the number of resulting star clusters and the distribution of their masses. The less extended remnants tend to have fewer, but more massive, clusters. Unlike the clusters, the gas distributions are relatively insensitive to the initial conditions, both the distribution of the bulk gas properties and the gas clumps. In general, cloud collisions are relatively insensitive to their initial conditions when modelled hydrodynamically in a dynamically important background medium. △ Less

Submitted 3 April, 2023; originally announced April 2023.

Comments: 22 pages, 23 figures. Accepted for publication in MNRAS

Rogue planets and brown dwarfs: Predicting the populations of free-floating planetary mass objects observable with JWST

Authors: Aleks Scholz, Koraljka Muzic, Ray Jayawardhana, Lyra Quinlan, James Wurster

Abstract: Free-floating (or rogue) planets are planets that are liberated (or ejected) from their host systems. Although simulations predict their existence in substantial numbers, direct observational evidence for free-floating planets with masses below ~5 MJup is still lacking. Several cycle-1 observing programs with JWST aim to hunt for them in four different star-forming clusters. These surveys are desi… ▽ More Free-floating (or rogue) planets are planets that are liberated (or ejected) from their host systems. Although simulations predict their existence in substantial numbers, direct observational evidence for free-floating planets with masses below ~5 MJup is still lacking. Several cycle-1 observing programs with JWST aim to hunt for them in four different star-forming clusters. These surveys are designed to be sensitive to masses of 1-15 MJup (assuming a hot-start formation), which corresponds to spectral types of early L to late T for the ages of these clusters. If the existing simulations are not wide off the mark, we show here that the planned programs are likely to find up to 10-20 giant rogue planets in moderate density clusters like NGC1333 or IC348, and several dozen to ~100 in high-density regions like NGC2024 and the Orion Nebula Cluster. These numbers correspond to 1-5% of the total cluster population; they could be substantially higher if stars form multiple giant planets at birth. In contrast, the number of free-floating brown dwarfs, formed from core collapse ("like stars") is expected to be significantly lower, only about 0.25% of the number of stars, or 1-7 for the clusters considered here. Below 10 MJup that number drops further by an order of magnitude. We also show that the planned surveys are not at risk of being significantly contaminated by field brown dwarfs in the foreground or background, after spectroscopic confirmation. Taken together, our results imply that if a population of L and T dwarfs were to be found in these JWST surveys, it is expected to be predominantly made up of rogue planets. △ Less

Submitted 23 September, 2022; v1 submitted 19 August, 2022; originally announced August 2022.

Comments: 16 pages, 3 figures, 2 tables, accepted by PASP

The initial magnetic criticality of prestellar cores

Authors: Felix D. Priestley, Charles Yin, James Wurster

Abstract: Direct observational measurements of the magnetic field strength in prestellar cores typically find supercritical mass-to-flux ratios, suggesting that the magnetic field is insufficient to prevent gravitational collapse. These measurements suffer from significant uncertainties; an alternative approach is to utilise the sensitivity of prestellar chemistry to the evolutionary history, and indirectly… ▽ More Direct observational measurements of the magnetic field strength in prestellar cores typically find supercritical mass-to-flux ratios, suggesting that the magnetic field is insufficient to prevent gravitational collapse. These measurements suffer from significant uncertainties; an alternative approach is to utilise the sensitivity of prestellar chemistry to the evolutionary history, and indirectly constrain the degree of magnetic support. We combine non-ideal magnetohydrodynamic simulations of prestellar cores with time-dependent chemistry and radiative transfer modelling, producing synthetic observations of the model cores in several commonly-observed molecular lines. We find that molecules strongly affected by freeze-out, such as CS and HCN, typically have much lower line intensities in magnetically subcritical models compared to supercritical ones, due to the longer collapse timescales. Subcritical models also produce much narrower lines for all species investigated. Accounting for a range of core properties, ages and viewing angles, we find that supercritical models are unable to reproduce the distribution of CS and N$_2$H$^+$ line strengths and widths seen in an observational sample, whereas subcritical models are in good agreement with the available data. This suggests that despite presently having supercritical mass-to-flux ratios, prestellar cores form as magnetically subcritical objects. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 10 pages, 10 figures. MNRAS accepted

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

Authors: James Wurster, Matthew R. Bate, Daniel J. Price, Ian A. Bonnell

Abstract: 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, magnetise… ▽ More 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. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: 19 pages, 17 figures, accepted for publication in MNRAS

The formation and early evolution of embedded star clusters in spiral galaxies

Authors: Steven Rieder, Clare Dobbs, Thomas Bending, Kong You Liow, James Wurster

Abstract: We present Ekster, a new method for simulating star clusters from birth in a live galaxy simulation that combines the smoothed-particle hydrodynamics (SPH) method Phantom with the $N$-body method PeTar. With Ekster, it becomes possible to simulate individual stars in a simulation with only moderately high resolution for the gas, allowing us to study whole sections of a galaxy rather than be restri… ▽ More We present Ekster, a new method for simulating star clusters from birth in a live galaxy simulation that combines the smoothed-particle hydrodynamics (SPH) method Phantom with the $N$-body method PeTar. With Ekster, it becomes possible to simulate individual stars in a simulation with only moderately high resolution for the gas, allowing us to study whole sections of a galaxy rather than be restricted to individual clouds. We use this method to simulate star and star cluster formation in spiral arms, investigating massive GMCs and spiral arm regions with lower mass clouds, from two galaxy models with different spiral potentials. After selecting these regions from pre-run galaxy simulations, we re-sample the particles to obtain a higher resolution. We then re-simulate these regions for 3 Myr to study where and how star clusters form. We analyse the early evolution of the embedded star clusters in these regions. We find that the massive GMC regions, which are more common with stronger spiral arms, form more massive clusters than the sections of spiral arms containing lower mass clouds. Clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

The impact of non-ideal magnetohydrodynamic processes on discs, outflows, counter-rotation and magnetic walls during the early stages of star formation

Authors: James Wurster, Matthew R. Bate, Ian A. Bonnell

Abstract: Non-ideal magnetohydrodynamic (MHD) processes -- namely Ohmic resistivity, ambipolar diffusion and the Hall effect -- modify the early stages of the star formation process and the surrounding environment. Collectively, they have been shown to promote disc formation and promote or hinder outflows. But which non-ideal process has the greatest impact? Using three-dimensional smoothed particle radiati… ▽ More Non-ideal magnetohydrodynamic (MHD) processes -- namely Ohmic resistivity, ambipolar diffusion and the Hall effect -- modify the early stages of the star formation process and the surrounding environment. Collectively, they have been shown to promote disc formation and promote or hinder outflows. But which non-ideal process has the greatest impact? Using three-dimensional smoothed particle radiation non-ideal MHD simulations, we model the gravitational collapse of a rotating, magnetised cloud through the first hydrostatic core phase to shortly after the formation of the stellar core. We investigate the impact of each process individually and collectively. Including any non-ideal process decreases the maximum magnetic field strength by at least an order of magnitude during the first core phase compared to using ideal MHD, and promotes the formation of a magnetic wall. When the magnetic field and rotation vectors are anti-aligned and the Hall effect is included, rotationally supported discs of $r \gtrsim 20$~au form; when only the Hall effect is included and the vectors are aligned, a counter-rotating pseudo-disc forms that is not rotationally supported. Rotationally supported discs of $r \lesssim 4$~au form if only Ohmic resistivity or ambipolar diffusion are included. The Hall effect suppresses first core outflows when the vectors are anti-aligned and suppresses stellar core outflows independent of alignment. Ohmic resistivity and ambipolar diffusion each promote first core outflows and delay the launching of stellar core outflows. Although each non-ideal process influences star formation, these results suggest that the Hall effect has the greatest influence. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: 19 pages, 13 figures. Accepted for publication in MNRAS

Investigating the role of magnetic fields in star formation using molecular line profiles

Authors: Charles Yin, Felix D. Priestley, James Wurster

Abstract: Determining the importance of magnetic fields in star forming environments is hampered by the difficulty of accurately measuring both field strength and gas properties in molecular clouds. We post-process three-dimensional non-ideal magnetohydrodynamic simulations of prestellar cores with a time-dependent chemical network, and use radiative transfer modelling to calculate self-consistent molecular… ▽ More Determining the importance of magnetic fields in star forming environments is hampered by the difficulty of accurately measuring both field strength and gas properties in molecular clouds. We post-process three-dimensional non-ideal magnetohydrodynamic simulations of prestellar cores with a time-dependent chemical network, and use radiative transfer modelling to calculate self-consistent molecular line profiles. Varying the initial mass-to-flux ratio from sub- to super-critical results in significant changes to both the intensity and shape of several observationally important molecular lines. We identify the peak intensity ratio of N$_2$H$^+$ to CS lines, and the CS $J=2-1$ blue-to-red peak intensity ratio, as promising diagnostics of the initial mass-to-flux ratio, with N$_2$H$^+$/CS values of $>0.6$ ($<0.2$) and CS blue/red values of $<3$ ($>5$) indicating subcritical (supercritial) collapse. These criteria suggest that, despite presently being magnetically supercritical, L1498 formed from subcritical initial conditions. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 8 pages of text, 10 pages in total. Accepted in MNRAS

The properties of clusters, and the orientation of magnetic fields relative to filaments, in magnetohydrodynamic simulations of colliding clouds

Authors: C. L. Dobbs, J. Wurster

Abstract: We have performed Smoothed Particle Magneto-Hydrodynamics (SPMHD) calculations of colliding clouds to investigate the formation of massive stellar clusters, adopting a timestep criterion to prevent large divergence errors. We find that magnetic fields do not impede the formation of young massive clusters (YMCs), and the development of high star formation rates, although we do see a strong dependen… ▽ More We have performed Smoothed Particle Magneto-Hydrodynamics (SPMHD) calculations of colliding clouds to investigate the formation of massive stellar clusters, adopting a timestep criterion to prevent large divergence errors. We find that magnetic fields do not impede the formation of young massive clusters (YMCs), and the development of high star formation rates, although we do see a strong dependence of our results on the direction of the magnetic field. If the field is initially perpendicular to the collision, and sufficiently strong, we find that star formation is delayed, and the morphology of the resulting clusters is significantly altered. We relate this to the large amplification of the field with this initial orientation. We also see that filaments formed with this configuration are less dense. When the field is parallel to the collision, there is much less amplification of the field, dense filaments form, and the formation of clusters is similar to the purely hydrodynamical case. Our simulations reproduce the observed tendency for magnetic fields to be aligned perpendicularly to dense filaments, and parallel to low density filaments. Overall our results are in broad agreement with past work in this area using grid codes. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 13 pages, 9 figures

Do we need non-ideal magnetohydrodynamics to model protostellar discs?

Authors: James Wurster

Abstract: We investigate and discuss protostellar discs in terms of where the various non-ideal magnetohydrodynamics (MHD) processes are important. We find that the traditional picture of a magnetised disc (where Ohmic resistivity is dominant near the mid-plane, surrounded by a region dominated by the Hall effect, with the remainder of the disc dominated by ambipolar diffusion) is a great oversimplification… ▽ More We investigate and discuss protostellar discs in terms of where the various non-ideal magnetohydrodynamics (MHD) processes are important. We find that the traditional picture of a magnetised disc (where Ohmic resistivity is dominant near the mid-plane, surrounded by a region dominated by the Hall effect, with the remainder of the disc dominated by ambipolar diffusion) is a great oversimplification. In simple parameterised discs, we find that the Hall effect is typically the dominant term throughout the majority of the disc. More importantly, we find that in much of our parameterised discs, at least two non-ideal processes have coefficients within a factor of 10 of one another, indicating that both are important and that naming a dominant term underplays the importance of the other terms. Discs that were self-consistently formed in our previous studies are also dominated by the Hall effect, and the ratio of ambipolar diffusion and Hall coefficients is typically less than 10, suggesting that both terms are equally important and listing a dominant term is misleading. These conclusions become more robust once the magnetic field geometry is taken into account. In agreement with the literature we review, we conclude that non-ideal MHD processes are important for the formation and evolution of protostellar discs. Ignoring any of the non-ideal processes, especially ambipolar diffusion and the Hall effect, yields an incorrect description of disc evolution. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS; 20 pages, 15 figures

Non-ideal magnetohydrodynamics vs turbulence II: Which is the dominant process in stellar core formation?

Authors: James Wurster, Benjamin T. Lewis

Abstract: Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of protostars by following the gravitational collapse of 1~M$_\odot$ gas clouds through the first hydrostatic core to stellar densities. The clouds are imposed with both rotational and turbulent velocities, and are th… ▽ More Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of protostars by following the gravitational collapse of 1~M$_\odot$ gas clouds through the first hydrostatic core to stellar densities. The clouds are imposed with both rotational and turbulent velocities, and are threaded with a magnetic field that is parallel/anti-parallel or perpendicular to the rotation axis; we investigate two rotation rates and four Mach numbers. The initial radius and mass of the stellar core are only weakly dependent on the initial parameters. In the models that include ideal MHD, the magnetic field strength implanted in the protostar at birth is much higher than observed, independent of the initial level of turbulence; only non-ideal MHD can reduce this strength to near or below the observed levels. This suggests that not only is ideal MHD an incomplete picture of star formation, but that the magnetic fields in low mass stars are implanted later in life by a dynamo process. Non-ideal MHD suppresses magnetically launched stellar core outflows, but turbulence permits thermally launched outflows to form a few years after stellar core formation. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS. 12 pages, 9 figures

Non-ideal magnetohydrodynamics vs turbulence I: Which is the dominant process in protostellar disc formation?

Authors: James Wurster, Benjamin T. Lewis

Abstract: Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of circumstellar discs that form nearly simultaneously with the formation of the protostar. This is done by modelling the gravitational collapse of a 1~M$_\odot$ gas cloud that is threaded with a magnetic field and im… ▽ More Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of circumstellar discs that form nearly simultaneously with the formation of the protostar. This is done by modelling the gravitational collapse of a 1~M$_\odot$ gas cloud that is threaded with a magnetic field and imposed with both rotational and turbulent velocities. We investigate magnetic fields that are parallel/anti-parallel and perpendicular to the rotation axis, two rotation rates and four Mach numbers. Disc formation occurs preferentially in the models that include non-ideal MHD where the magnetic field is anti-parallel or perpendicular to the rotation axis. This is independent of the initial rotation rate and level of turbulence, suggesting that subsonic turbulence plays a minimal role in influencing the formation of discs. Aside from first core outflows which are influenced by the initial level of turbulence, non-ideal MHD processes are more important than turbulent processes during the formation of discs around low-mass stars. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS. 12 pages, 8 figures

There is no magnetic braking catastrophe: Low-mass star cluster and protostellar disc formation with non-ideal magnetohydrodynamics

Authors: James Wurster, Matthew R. Bate, Daniel J. Price

Abstract: We present results from the first radiation non-ideal magnetohydrodynamics (MHD) simulations of low-mass star cluster formation that resolve the fragmentation process down to the opacity limit. We model 50~M$_\odot$ turbulent clouds initially threaded by a uniform magnetic field with strengths of 3, 5 10 and 20 times the critical mass-to-magnetic flux ratio, and at each strength, we model both an… ▽ More We present results from the first radiation non-ideal magnetohydrodynamics (MHD) simulations of low-mass star cluster formation that resolve the fragmentation process down to the opacity limit. We model 50~M$_\odot$ turbulent clouds initially threaded by a uniform magnetic field with strengths of 3, 5 10 and 20 times the critical mass-to-magnetic flux ratio, and at each strength, we model both an ideal and non-ideal (including Ohmic resistivity, ambipolar diffusion and the Hall effect) MHD cloud. Turbulence and magnetic fields shape the large-scale structure of the cloud, and similar structures form regardless of whether ideal or non-ideal MHD is employed. At high densities ($10^6 \lesssim n_{\rm H} \lesssim 10^{11}$~cm$^{-3}$), all models have a similar magnetic field strength versus density relation, suggesting that the field strength in dense cores is independent of the large-scale environment. Albeit with limited statistics, we find no evidence for the dependence of the initial mass function on the initial magnetic field strength, however, the star formation rate decreases for models with increasing initial field strengths; the exception is the strongest field case where collapse occurs primarily along field lines. Protostellar discs with radii $\gtrsim 20$~au form in all models, suggesting that disc formation is dependent on the gas turbulence rather than on magnetic field strength. We find no evidence for the magnetic braking catastrophe, and find that magnetic fields do not hinder the formation of protostellar discs. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 24 pages, 19 figures. Accepted for publication in MNRAS. Videos available at https://www.youtube.com/watch?v=CLhmaOhj5RU&list=PLwI7am9c6sBjW69kfhfRDMR3Gz6nzOPOi

Ambipolar diffusion and the molecular abundances in prestellar cores

Authors: Felix D. Priestley, James Wurster, Serena Viti

Abstract: We investigate differences in the molecular abundances between magnetically super- and sub-critical prestellar cores, performing three-dimensional non-ideal magnetohydrodynamical (MHD) simulations with varying densities and magnetic field strengths, and post-processing the results with a time-dependent gas-grain chemical code. Most molecular species show significantly more central depletion in sub… ▽ More We investigate differences in the molecular abundances between magnetically super- and sub-critical prestellar cores, performing three-dimensional non-ideal magnetohydrodynamical (MHD) simulations with varying densities and magnetic field strengths, and post-processing the results with a time-dependent gas-grain chemical code. Most molecular species show significantly more central depletion in subcritical models, due to the longer duration of collapse. However, the directly observable quantities - the molecule to hydrogen column density ratios - are generally too similar for observational data to discriminate between models. The profiles of N$_2$H$^+$ and HCO$^+$ show qualitative differences between supercritical and subcritical models on scales of $0.01 \, {\rm pc}$, which may allow the two cases to be distinguished. However, this requires knowledge of the hydrogen column density, which is not directly measureable, and predicted line intensity profiles from radiative transfer modelling are similar for these molecules. Other commonly observed species, such as HCN and CH$_3$OH, have line intensity profiles which differ more strongly between models, and so are more promising as tracers of the mechanism of cloud collapse. △ Less

Submitted 5 July, 2019; originally announced July 2019.

Comments: 8 pages, 11 figures. MNRAS accepted 05/07/19

Resolving numerical star formation: A cautionary tale

Authors: James Wurster, Matthew R. Bate

Abstract: Resolution studies of test problems set baselines and help define minimum resolution requirements, however, resolution studies must also be performed on scientific simulations to determine the effect of resolution on the specific scientific results. We perform a resolution study on the formation of a protostar by modelling the collapse of gas through 14 orders of magnitude in density. This is done… ▽ More Resolution studies of test problems set baselines and help define minimum resolution requirements, however, resolution studies must also be performed on scientific simulations to determine the effect of resolution on the specific scientific results. We perform a resolution study on the formation of a protostar by modelling the collapse of gas through 14 orders of magnitude in density. This is done using compressible radiative non-ideal magnetohydrodynamics. Our suite consists of an ideal magnetohydrodynamics (MHD) model and two non-ideal MHD models, and we test three resolutions for each model. The resulting structure of the ideal MHD model is approximately independent of resolution, although higher magnetic field strengths are realised in higher resolution models. The non-ideal MHD models are more dependent on resolution, specifically the magnetic field strength and structure. Stronger magnetic fields are realised in higher resolution models, and the evolution of detailed structures such as magnetic walls are only resolved in our highest resolution simulation. In several of the non-ideal MHD models, there is an off-set between the location of the maximum magnetic field strength and the maximum density, which is often obscured or lost at lower resolutions. Thus, understanding the effects of resolution on numerical star formation is imperative for understanding the formation of a star. △ Less

Submitted 28 June, 2019; originally announced June 2019.

Comments: 8 pages, 5 figures. Proceedings of the "14th international SPHERIC workshop", Exeter, United Kingdom, 25-27 June 2019

Disc formation and fragmentation using radiative non-ideal magnetohydrodynamics

Authors: James Wurster, Matthew R. Bate

Abstract: We investigate the formation and fragmentation of discs using a suite of three-dimensional smoothed particle radiative magnetohydrodynamics simulations. Our models are initialised as 1M$_\odot$ rotating Bonnor-Ebert spheres that are threaded with a uniform magnetic field. We examine the effect of including ideal and non-ideal magnetic fields, the orientation and strength of the magnetic field, and… ▽ More We investigate the formation and fragmentation of discs using a suite of three-dimensional smoothed particle radiative magnetohydrodynamics simulations. Our models are initialised as 1M$_\odot$ rotating Bonnor-Ebert spheres that are threaded with a uniform magnetic field. We examine the effect of including ideal and non-ideal magnetic fields, the orientation and strength of the magnetic field, and the initial rotational rate. We follow the gravitational collapse and early evolution of each system until the final classification of the protostellar disc can be determined. Of our 105 models, 41 fragment, 21 form a spiral structure but do not fragment, and another 12 form smooth discs. Fragmentation is more likely to occur for faster initial rotation rates and weaker magnetic fields. For stronger magnetic field strengths, the inclusion of non-ideal MHD promotes disc formation, and several of these models fragment, whereas their ideal MHD counterparts do not. For the models that fragment, there is no correlation between our parameters and where or when the fragmentation occurs. Bipolar outflows are launched in only 17 models, and these models have strong magnetic fields that are initially parallel to the rotation axis. Counter-rotating envelopes form in four slowly-rotating, strong-field models -- including one ideal MHD model -- indicating they form only in a small fraction of the parameter space investigated. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 17 pages, 10 figures, accepted for publication in MNRAS

The role of magnetic fields in the formation of protostellar discs

Authors: James Wurster, Zhi-Yun Li

Abstract: Truncated abstract: The formation of a protostellar disc is a natural outcome during the star formation process. As gas in a molecular cloud core collapses under self-gravity, the angular momentum of the gas will slow its collapse on small scales and promote the formation of a protostellar disc. Although the angular momenta of dense star-forming cores remain to be fully characterized observational… ▽ More Truncated abstract: The formation of a protostellar disc is a natural outcome during the star formation process. As gas in a molecular cloud core collapses under self-gravity, the angular momentum of the gas will slow its collapse on small scales and promote the formation of a protostellar disc. Although the angular momenta of dense star-forming cores remain to be fully characterized observationally, existing data indicates that typical cores have enough angular momenta to form relatively large, rotationally supported discs. However, molecular clouds are observed to be permeated by magnetic fields, which can strongly affect the evolution of angular momentum through magnetic braking. Indeed, in the ideal MHD limit, magnetic braking has been shown to be so efficient as to remove essentially all of the angular momentum of the material close to the forming star such that disc formation is suppressed. This is known as the magnetic braking catastrophe. The catastrophe must be averted in order for the all-important rotationally supported discs to appear, but when and how this happens remains debated. We review the resolutions proposed to date, with emphasis on misalignment, turbulence and especially non-ideal effects. The dissipative non-ideal effects weaken the magnetic field, and the dispersive term redirects it to promote or hinder disc formation. When self-consistently applying non-ideal processes, rotationally supported discs of at least tens of au form, thus preventing the magnetic braking catastrophe. The non-ideal processes are sensitive to the magnetic field strength, cosmic ray ionization rate, and gas and dust grain properties, thus a complete understanding of the host molecular cloud is required. Therefore, the properties of the host molecular cloud -- and especially its magnetic field -- cannot be ignored when numerically modelling the formation and evolution of protostellar discs. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: Review article. Published version available at https://www.frontiersin.org/articles/10.3389/fspas.2018.00039/full

Journal ref: Frontiers in Astronomy and Space Science, 5:39, December 2018

On the origin of magnetic fields in stars

Authors: James Wurster, Matthew R. Bate, Daniel J. Price

Abstract: 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? We use radiation non-ideal magnetohydrodynamics simulations of the gravitational collapse of a rotating, magnetized molecular cloud core over 17 orders of magnitude in density, past the first hydrostatic core to the formation of the second, stellar core, to e… ▽ More 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? We use radiation non-ideal magnetohydrodynamics simulations of the gravitational collapse of a rotating, magnetized molecular cloud core over 17 orders of magnitude in density, past the first hydrostatic core to the formation of the second, stellar core, to examine the fossil field hypothesis. Whereas in previous work we found that magnetic fields in excess of 10 kG can be implanted in stars at birth, this assumed ideal magnetohydrodynamics (MHD), i.e. that the gas is coupled to the magnetic field. Here we present non-ideal MHD calculations which include Ohmic resistivity, ambipolar diffusion and the Hall effect. For realistic cosmic ray ionization rates, we find that magnetic field strengths of $\lesssim$ kG are implanted in the stellar core at birth, ruling out a strong fossil field. While these results remain sensitive to resolution, they cautiously provide evidence against a fossil field origin for stellar magnetic fields, suggesting instead that magnetic fields in stars originate in a dynamo process. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 8 pages, 8 figures. Accepted for publication in MNRAS

Hall effect-driven formation of gravitationally unstable discs in magnetized molecular cloud cores

Authors: James Wurster, Matthew R. Bate, Daniel J. Price

Abstract: We demonstrate the formation of gravitationally unstable discs in magnetized molecular cloud cores with initial mass-to-flux ratios of 5 times the critical value, effectively solving the magnetic braking catastrophe. We model the gravitational collapse through to the formation of the stellar core, using Ohmic resistivity, ambipolar diffusion and the Hall effect and using the canonical cosmic ray i… ▽ More We demonstrate the formation of gravitationally unstable discs in magnetized molecular cloud cores with initial mass-to-flux ratios of 5 times the critical value, effectively solving the magnetic braking catastrophe. We model the gravitational collapse through to the formation of the stellar core, using Ohmic resistivity, ambipolar diffusion and the Hall effect and using the canonical cosmic ray ionization rate of $ζ_\text{cr} = 10^{-17}$ s$^{-1}$. When the magnetic field and rotation axis are initially aligned, a $\lesssim1$~au disc forms after the first core phase, whereas when they are anti-aligned, a gravitationally-unstable 25~au disc forms during the first core phase. The aligned model launches a 3~km~s$^{-1}$ first core outflow, while the anti-aligned model launches only a weak $\lesssim 0.3$~km~s$^{-1}$ first core outflow. Qualitatively, we find that models with $ζ_\text{cr} = 10^{-17}$ s$^{-1}$ are similar to purely hydrodynamical models if the rotation axis and magnetic field are initially anti-aligned, whereas they are qualitatively similar to ideal magnetohydrodynamical models if initially aligned. △ Less

Submitted 13 August, 2018; originally announced August 2018.

Comments: 9 pages, 9 figures. Accepted for publication in MNRAS

The effect of extreme ionisation rates during the initial collapse of a molecular cloud core

Authors: James Wurster, Matthew R. Bate, Daniel J. Price

Abstract: 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 rota… ▽ More 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μ$m, and our free parameter is the cosmic ray ionisation rate, $ζ_\text{cr}$. We evolve our models, where possible, until they have produced a first hydrostatic core. Models with $ζ_\text{cr}\gtrsim10^{-13}$ s$^{-1}$ are indistinguishable from ideal MHD models and the evolution of the model with $ζ_\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 ($ζ_\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. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: 13 pages, 15 figures, accepted for publication in MNRAS

The collapse of a molecular cloud core to stellar densities using radiation non-ideal magnetohydrodynamics

Authors: James Wurster, Matthew R. Bate, Daniel J. Price

Abstract: 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 temper… ▽ More 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 $ζ_\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. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Comments: 22 pages, 24 figures. Accepted for publication in MNRAS

Investigating prescriptions for artificial resistivity in smoothed particle magnetohydrodynamics

Authors: James Wurster, Matthew R. Bate, Daniel J. Price, Terrence S. Tricco

Abstract: In numerical simulations, artificial terms are applied to the evolution equations for stability. To prove their validity, these terms are thoroughly tested in test problems where the results are well known. However, they are seldom tested in production-quality simulations at high resolution where they interact with a plethora of physical and numerical algorithms. We test three artificial resistivi… ▽ More In numerical simulations, artificial terms are applied to the evolution equations for stability. To prove their validity, these terms are thoroughly tested in test problems where the results are well known. However, they are seldom tested in production-quality simulations at high resolution where they interact with a plethora of physical and numerical algorithms. We test three artificial resistivities in both the Orszag-Tang vortex and in a star formation simulation. From the Orszag-Tang vortex, the Price et. al. (2017) artificial resistivity is the least dissipative thus captures the density and magnetic features; in the star formation algorithm, each artificial resistivity algorithm interacts differently with the sink particle to produce various results, including gas bubbles, dense discs, and migrating sink particles. The star formation simulations suggest that it is important to rely upon physical resistivity rather than artificial resistivity for convergence. △ Less

Submitted 23 June, 2017; originally announced June 2017.

Comments: 8 pages, 7 figures. Proceedings of the "12th international SPHERIC workshop", Ourense, Spain, 13-15 June 2017

Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics

Authors: Daniel J. Price, James Wurster, Terrence S. Tricco, Chris Nixon, Stéven Toupin, Alex Pettitt, Conrad Chan, Daniel Mentiplay, Guillaume Laibe, Simon Glover, Clare Dobbs, Rebecca Nealon, David Liptai, Hauke Worpel, Clément Bonnerot, Giovanni Dipierro, Giulia Ballabio, Enrico Ragusa, Christoph Federrath, Roberto Iaconi, Thomas Reichardt, Duncan Forgan, Mark Hutchison, Thomas Constantino, Ben Ayliffe , et al. (2 additional authors not shown)

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

Submitted 15 June, 2018; v1 submitted 13 February, 2017; originally announced February 2017.

Comments: 88 pages, 60 figures, accepted to PASA. Code available from https://phantomsph.bitbucket.io/

The impact of non-ideal magnetohydrodynamics on binary star formation

Authors: James Wurster, Daniel J. Price, Matthew R. Bate

Abstract: 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 f… ▽ More 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. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: Accepted to MNRAS. 18 pages, 20 figures, 1 appendix

Does turbulence determine the initial mass function?

Authors: David Liptai, Daniel J. Price, James Wurster, Matthew R. Bate

Abstract: 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 l… ▽ More 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(ρ)$. This suggests that the density PDF, while determining the star formation rate, is not the primary driver of the IMF. △ Less

Submitted 26 October, 2016; v1 submitted 24 October, 2016; originally announced October 2016.

Comments: 6 pages, 7 figures, accepted to MNRAS

NICIL: A stand alone library to self-consistently calculate non-ideal magnetohydrodynamic coefficients in molecular cloud cores

Authors: James Wurster

Abstract: In this paper, we introduce Nicil: Non-Ideal magnetohydrodynamics Coefficients and Ionisation Library. Nicil is a stand-alone Fortran90 module that calculates the ionisation values and the coefficients of the non-ideal magnetohydrodynamics terms of Ohmic resistivity, the Hall effect, and ambipolar diffusion. The module is fully parameterised such that the user can decide which processes to include… ▽ More In this paper, we introduce Nicil: Non-Ideal magnetohydrodynamics Coefficients and Ionisation Library. Nicil is a stand-alone Fortran90 module that calculates the ionisation values and the coefficients of the non-ideal magnetohydrodynamics terms of Ohmic resistivity, the Hall effect, and ambipolar diffusion. The module is fully parameterised such that the user can decide which processes to include and decide upon the values of the free parameters, making this a versatile and customisable code. The module includes both cosmic ray and thermal ionisation; the former includes two ion species and three species of dust grains (positively charged, negatively charged and neutral), and the latter includes five elements which can be doubly ionised. We demonstrate tests of the module, and then describe how to implement it into an existing numerical code. △ Less

Submitted 6 September, 2016; v1 submitted 2 August, 2016; originally announced August 2016.

Comments: 13 pgs; accepted to PASA; source code for NICIL: https://bitbucket.org/jameswurster/nicil ; Corrected typos in eqn 47 and in initial conditions of section 4

Comparing Simulations of AGN Feedback

Authors: Mark L. A. Richardson, Evan Scannapieco, Julien Devriendt, Adrianne Slyz, Robert J. Thacker, Yohan Dubois, James Wurster, Joseph Silk

Abstract: We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slight… ▽ More We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core. △ Less

Submitted 11 May, 2016; originally announced May 2016.

Comments: 22 pages, 20 figures, 3 tables, Accepted to ApJ, comments welcome

The effect of a wider initial separation on common envelope binary interaction simulations

Authors: Roberto Iaconi, Thomas Reichardt, Jan Staff, Orsola De Marco, Jean-Claude Passy, Daniel Price, James Wurster, Falk Herwig

Abstract: We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the sim… ▽ More We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the simulation outcome. We conclude that for both codes the post-common envelope separation is somewhat larger and the amount of unbound mass slightly greater when the initial separation is wide enough that the giant does not yet overflow or just overflows its Roche lobe. PHANTOM has been adapted to the common envelope problem here for the first time and a full comparison with ENZO is presented, including an investigation of convergence as well as energy and angular momentum conservation. We also set our simulations in the context of past simulations. This comparison reveals that it is the expansion of the giant before rapid in-spiral and not spinning up of the star that causes a larger final separation. We also suggest that the large range in unbound mass for different simulations is difficult to explain and may have something to do with simulations that are not fully converged. △ Less

Submitted 16 September, 2016; v1 submitted 7 March, 2016; originally announced March 2016.

Comments: Accepted by MNRAS

Can non-ideal magnetohydrodynamics solve the magnetic braking catastrophe?

Authors: James Wurster, Daniel J. Price, Matthew R. Bate

Abstract: 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 th… ▽ More 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. △ Less

Submitted 12 January, 2016; v1 submitted 4 December, 2015; originally announced December 2015.

Comments: Accepted by MNRAS. Many figures have been slightly modified from the previous version. 26 pages, 22 figures

Ambipolar diffusion in smoothed particle magnetohydrodynamics

Authors: James Wurster, Daniel J. Price, Ben A. Ayliffe

Abstract: 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, allo… ▽ More 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. △ Less

Submitted 8 August, 2014; originally announced August 2014.

Comments: Accepted for publication by MNRAS

AGN Feedback models: Correlations with star formation and observational implications of time evolution

Authors: Robert J. Thacker, C. MacMackin, James Wurster, Alexander Hobbs

Abstract: We examine the correlation between the star formation rate (SFR) and black hole accretion rate (BHAR) across a suite of different AGN feedback models, using the time evolution of a merger simulation. By considering three different stages of evolution, and a distinction between the nuclear and outer regions of star formation, we consider 63 different cases. Despite many of the feedback models fitti… ▽ More We examine the correlation between the star formation rate (SFR) and black hole accretion rate (BHAR) across a suite of different AGN feedback models, using the time evolution of a merger simulation. By considering three different stages of evolution, and a distinction between the nuclear and outer regions of star formation, we consider 63 different cases. Despite many of the feedback models fitting the M-σ relationship well, there are often distinct differences in the SFR-BHAR correlations, with close to linear trends only being present after the merger. Some of the models also show evolution in the SFR-BHAR parameter space that is at times directly across the long-term averaged SFR-BHAR correlation. This suggests that the observational SFR-BHAR correlation found for ensembles of galaxies is an approximate statistical trend, as suggested by Hickox et al. Decomposing the SFR into nuclear and outer components also highlights notable differences between models and there is only modest agreement with observational studies examining this in Seyfert galaxies. For the fraction of the black hole mass growth from the merger event relative to the final black hole mass, we find as much as a factor of three variation among models. This also translates into a similar variation in the post-starburst black hole mass growth. Overall, we find that while qualitative features are often similar amongst models, precise quantitative analysis shows there can be quite distinct differences. △ Less

Submitted 2 July, 2014; originally announced July 2014.

Comments: Accepted to MNRAS. Comments welcome

Cloud angular momentum and effective viscosity in global SPH simulations with feedback

Authors: David J. Williamson, Robert J. Thacker, James Wurster, Brad K. Gibson

Abstract: We examine simulations of isolated galaxies to analyse the effects of localised feedback on the formation and evolution of molecular clouds. Feedback contributes to turbulence and the destruction of clouds, leading to a population of clouds that is younger, less massive, and with more retrograde rotation. We investigate the evolution of clouds as they interact with each other and the diffuse ISM,… ▽ More We examine simulations of isolated galaxies to analyse the effects of localised feedback on the formation and evolution of molecular clouds. Feedback contributes to turbulence and the destruction of clouds, leading to a population of clouds that is younger, less massive, and with more retrograde rotation. We investigate the evolution of clouds as they interact with each other and the diffuse ISM, and determine that the role of cloud interactions differs strongly with the presence of feedback: in models without feedback, scattering events dramatically increase the retrograde fraction, but in models with feedback, mergers between clouds may slightly increase the prograde fraction. We also produce an estimate of the viscous time-scale due to cloud-cloud collisions, which increases with increasing strength of feedback (~20 Gyr vs ~10 Gyr), but is still much smaller than previous estimates (~1000 Gyr); although collisions become more frequent with feedback, less energy is lost in each collision than in the models without feedback. △ Less

Submitted 12 June, 2014; v1 submitted 22 April, 2014; originally announced April 2014.

Comments: Accepted by MNRAS

A Comparative Study of AGN Feedback Algorithms

Authors: James Wurster, Robert J. Thacker

Abstract: Modelling AGN feedback in numerical simulations is both technically and theoretically challenging, with numerous approaches having been published in the literature. We present a study of five distinct approaches to modelling AGN feedback within gravitohydrodynamic simulations of major mergers of Milky Way-sized galaxies. To constrain differences to only be between AGN feedback models, all simulati… ▽ More Modelling AGN feedback in numerical simulations is both technically and theoretically challenging, with numerous approaches having been published in the literature. We present a study of five distinct approaches to modelling AGN feedback within gravitohydrodynamic simulations of major mergers of Milky Way-sized galaxies. To constrain differences to only be between AGN feedback models, all simulations start from the same initial conditions and use the same star formation algorithm. Most AGN feedback algorithms have five key aspects: black hole accretion rate, energy feedback rate and method, particle accretion algorithm, black hole advection algorithm, and black hole merger algorithm. All models follow different accretion histories, with accretion rates that differ by up to three orders of magnitude at any given time. We consider models with either thermal or kinetic feedback, with the associated energy deposited locally around the black hole. Each feedback algorithm modifies the gas properties near the black hole to different extents. The particle accretion algorithms usually (but not always) maintain good agreement between the mass accreted by \dot{M} dt and the mass of gas particles removed from the simulation. The black hole advection algorithms dampen inappropriate dragging of the black holes by two-body interactions. Advecting the black hole a limited distance based upon local mass distributions has many desirably properties. The black holes merge when given criteria are met, and we find a range of merger times for different criteria. Using the M_{BH}-σrelation as a diagnostic of the remnants yields three models that lie within the one-sigma scatter of the observed relation and two that fall below it. The wide variation in accretion behaviours of the models reinforces the fact that there remains much to be learnt about the evolution of galactic nuclei. (abridged) △ Less

Submitted 21 February, 2013; originally announced February 2013.

Comments: 24 pages, accepted to MNRAS

Accretion Disc Particle Accretion in Major Merger Simulations

Authors: James Wurster, Rob J. Thacker

Abstract: A recent approach to simulating localized feedback from active galactic nuclei by Power et al. (2011) uses an accretion disc particle to represent both the black hole and its accretion disc. We have extrapolated and adapted this approach to simulations of Milky Way-sized galaxy mergers containing black holes and explored the impact of the various parameters in this model as well as its resolution… ▽ More A recent approach to simulating localized feedback from active galactic nuclei by Power et al. (2011) uses an accretion disc particle to represent both the black hole and its accretion disc. We have extrapolated and adapted this approach to simulations of Milky Way-sized galaxy mergers containing black holes and explored the impact of the various parameters in this model as well as its resolution dependence. The two key parameters in the model are an effective accretion radius, which determines the radius within which gas particles are added to the accretion disc, and a viscous time-scale which determines how long it takes for material in the accretion disc to accrete on to the black hole itself. We find that there is a limited range of permitted accretion radii and viscous time-scales, with unphysical results produced outside this range. For permitted model parameters, the nuclear regions of simulations with the same resolution follow similar evolutionary paths, producing final black hole masses that are consistent within a factor of two. When comparing the resolution dependence of the model, there is a trend towards higher resolution producing slightly lower mass black holes, but values for the two resolutions studied again agree within a factor of two. We also compare these results to two other AGN feedback algorithms found in the literature. While the evolution of the systems vary, most notably the intermediate total black hole mass, the final black hole masses differ by less than a factor of five amongst all of our models, and the remnants exhibit similar structural parameters. The implication of this accretion model is that, unlike most accretion algorithms, a decoupling of the accretion rate on to the black hole and the local gas properties is permitted and obtained; this allows for black hole growth even after feedback has prevented additional accretion events on to the disc. △ Less

Submitted 30 January, 2013; originally announced January 2013.

Comments: 17 pages, accepted to MNRAS

Magnetically-Regulated Fragmentation Induced by Nonlinear Flows and Ambipolar Diffusion

Authors: Shantanu Basu, Glenn E. Ciolek, Wolf Dapp, James Wurster

Abstract: We present a parameter study of simulations of fragmentation regulated by gravity, magnetic fields, ambipolar diffusion, and nonlinear flows. The thin-sheet approximation is employed with periodic lateral boundary conditions, and the nonlinear flow field ("turbulence") is allowed to freely decay. In agreement with previous results in the literature, our results show that the onset of runaway col… ▽ More We present a parameter study of simulations of fragmentation regulated by gravity, magnetic fields, ambipolar diffusion, and nonlinear flows. The thin-sheet approximation is employed with periodic lateral boundary conditions, and the nonlinear flow field ("turbulence") is allowed to freely decay. In agreement with previous results in the literature, our results show that the onset of runaway collapse (formation of the first star) in subcritical clouds is significantly accelerated by nonlinear flows in which a large-scale wave mode dominates the power spectrum. In addition, we find that a power spectrum with equal energy on all scales also accelerates collapse, but by a lesser amount. For a highly super-Alfvenic initial velocity field with most power on the largest scales, the collapse occurs promptly during the initial compression wave. However, for trans-Alfvenic perturbations, a subcritical magnetic field causes a rebound from the initial compression, and the system undergoes several oscillations before runaway collapse occurs. Models that undergo prompt collapse have highly supersonic infall motions at the core boundaries. Cores in magnetically subcritical models with trans-Alfvenic initial perturbations also pick up significant systematic speeds by inheriting motions associated with magnetically-driven oscillations. Core mass distributions are much broader than in models with small-amplitude initial perturbations, although the disturbed structure of cores that form due to nonlinear flows does not guarantee subsequent monolithic collapse. Our simulations also demonstrate that significant power can (if present initially) be maintained with negligible dissipation in large-scale compressive modes of a magnetic thin sheet, in the limit of perfect flux freezing. △ Less

Submitted 13 January, 2009; v1 submitted 4 October, 2008; originally announced October 2008.

Comments: v2, 26 pages, 12 figures, to appear in New Astronomy, animations at http://www.astro.uwo.ca/~basu/pb.htm

Journal ref: New Astronomy 14 (2009) 483-495

Nonlinear Evolution of Gravitational Fragmentation Regulated by Magnetic Fields and Ambipolar Diffusion

Authors: Shantanu Basu, Glenn E. Ciolek, James Wurster

Abstract: We present results from an extensive set of simulations of gravitational fragmentation in the presence of magnetic fields and ambipolar diffusion. The average fragmentation spacing in the nonlinear phase of evolution is in excellent agreement with the prediction of linear perturbation theory. The time scale for nonlinear growth and runaway of the first core is $\approx 10$ times the calculated g… ▽ More We present results from an extensive set of simulations of gravitational fragmentation in the presence of magnetic fields and ambipolar diffusion. The average fragmentation spacing in the nonlinear phase of evolution is in excellent agreement with the prediction of linear perturbation theory. The time scale for nonlinear growth and runaway of the first core is $\approx 10$ times the calculated growth time $\taugm$ of the eigenmode with minimum growth time, when starting from a uniform background state with small-amplitude white-noise perturbations. Subcritical and transcritical models typically evolve on a significantly longer time scale than the supercritical models. Infall motions in the nonlinear fully-developed contracting cores are subsonic on the core scale in subcritical and transcritical clouds, but are somewhat supersonic in supercritical clouds. Core mass distributions are sharply peaked with a steep decline to large masses, consistent with the existence of a preferred mass scale for each unique set of dimensionless free parameters. However, a sum total of results for various initial mass-to-flux ratios yields a broad distribution reminiscent of observed core mass distributions. Based on our results, we conclude that fragmentation spacings, magnitude of infall motions, core shapes, and, especially, the curvature of magnetic field morphology, may serve as indirect observational means of determining a cloud's ambient mass-to-flux ratio. △ Less

Submitted 11 July, 2008; v1 submitted 15 June, 2008; originally announced June 2008.

Comments: v2, 13 figures, New Astronomy, animations can be obtained at http://www.astro.uwo.ca/~basu/pb.htm

Journal ref: NewAstron.14:221-237,2009