BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthew R. Bate

Matthew R. Bate appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.EPastro-phastro-ph.IMastro-ph.COastro-ph.HEphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

53works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

53 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dust coagulation during the early stages of star formation: molecular cloud collapse and first hydrostatic core evolution

Planet formation in protoplanetary discs requires dust grains to coagulate from the sub-micron sizes that are found in the interstellar medium into much larger objects. For the first time, we study the growth of dust grains during the earliest phases of star formation using three-dimensional hydrodynamical simulations. We begin with a typical interstellar dust grain size distribution and study dust growth during the collapse of a molecular cloud core and the evolution of the first hydrostatic core, prior to the formation of the stellar core. We examine how the dust size distribution evolves both spatially and temporarily. We find that the envelope maintains its initial population of small dust grains with little growth during these phases, except that in the inner few hundreds of au the smallest grains are depleted. However, once the first hydrostatic core forms rapid dust growth to sizes in excess of $100~μ$m occurs within the core (before stellar core formation). Progressively larger grains are produced at smaller distances from the centre of the core. In rapidly-rotating molecular cloud cores, the `first hydrostatic core' that forms is better described as a pre-stellar disc that may be gravitationally unstable. In such cases, grain growth is more rapid in the spiral density waves leading to the larger grains being preferentially found in the spiral waves even though there is no migration of grains relative to the gas. Thus, the grain size distribution can vary substantially in the first core/pre-stellar disc even at these very early times.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Supernovae and photoionizing feedback in spiral arm molecular clouds

We explore the interplay between supernovae and the ionizing radiation of their progenitors in star forming regions. The relative contributions of these stellar feedback processes are not well understood, particularly on scales greater than a single star forming cloud. We focus predominantly on how they affect the interstellar medium. We re-simulate a 500 pc^2 region from previous work that included photoionization and add supernovae. Over the course of 10 Myr more than 500 supernovae occur in the region. The supernovae remnants cool very quickly in the absence of earlier photoionization, but form much larger and more spherical hot bubbles when photoionization is present. Overall, the photoionization has a significantly greater effect on gas morphology and the sites of star formation. However, the two processes are comparable when looking at their effect on velocity dispersion. When combined, the two feedback processes increase the velocity dispersions by more than the sum of their parts, particularly on scales above 5 pc.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ionization and dissociation induced fragmentation of a tidally disrupted star into planets around a supermassive black hole

We show results from the radiation hydrodynamics (RHD) simulations of tidal disruption of a star on a parabolic orbit by a supermassive black hole (SMBH) based on a three-dimensional smoothed particle hydrodynamics code with radiative transfer. We find that such a tidally disrupted star fragment and form clumps soon after its tidal disruption. The fragmentation results from the endothermic processes of ionization and dissociation that reduce the gas pressure, leading to local gravitational collapse. Radiative cooling is less effective because the stellar debris is still highly optically thick in such an early time. Our simulations reveal that a solar-type star with a stellar density profile of n=3 disrupted by a 10^6 solar mass black hole produces $\sim20$ clumps of masses in the range of 0.1 to 12 Jupiter masses. The mass fallback rate decays with time, with pronounced spikes from early to late time. The spikes provide evidence for the clumps of the returning debris, while the clumps on the unbound debris can be potentially freely-floating planets and brown dwarfs. This ionization and dissociation induced fragmentation on a tidally disrupted star are a promising candidate mechanism to form low-mass stars to planets around an SMBH.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoionising feedback in spiral arm molecular clouds

We present simulations of a 500 pc$^2$ region, containing gas of mass 4 $\times$ 10$^6$ M$_\odot$, extracted from an entire spiral galaxy simulation, scaled up in resolution, including photoionising feedback from stars of mass > 18 M$_\odot$. Our region is evolved for 10 Myr and shows clustered star formation along the arm generating $\approx$ 5000 cluster sink particles $\approx$ 5% of which contain at least one of the $\approx$ 4000 stars of mass > 18 M$_\odot$. Photoionisation has a noticeable effect on the gas in the region, producing ionised cavities and leading to dense features at the edge of the HII regions. Compared to the no-feedback case, photoionisation produces a larger total mass of clouds and clumps, with around twice as many such objects, which are individually smaller and more broken up. After this we see a rapid decrease in the total mass in clouds and the number of clouds. Unlike studies of isolated clouds, our simulations follow the long range effects of ionisation, with some already-dense gas becoming compressed from multiple sides by neighbouring HII regions. This causes star formation that is both accelerated and partially displaced throughout the spiral arm with up to 30% of our cluster sink particle mass forming at distances > 5 pc from sites of sink formation in the absence of feedback. At later times, the star formation rate decreases to below that of the no-feedback case.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the dynamics of dust during protostellar collapse

The dynamics of dust and gas can be quite different from each other when the dust is poorly coupled to the gas. In protoplanetary discs, it is well known that this decoupling of the dust and gas can lead to diverse spatial structures and dust-to-gas ratios. In this paper, we study the dynamics of dust and gas during the earlier phase of protostellar collapse, before a protoplanetary disc is formed. We find that for dust grains with sizes < 10 micron, the dust is well coupled during the collapse of a rotating, pre-stellar core and there is little variation of the dust-to-gas ratio during the collapse. However, if larger grains are present, they may have trajectories that are very different from the gas during the collapse, leading to mid-plane settling and/or oscillations of the dust grains through the mid-plane. This may produce variations in the dust-to-gas ratio and very different distributions of large and small dust grains at the very earliest stages of star formation, if large grains are present in pre-stellar cores.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Smoothed Particle Magnetohydrodynamics: A State of the Union

Obtaining a stable magnetohydrodynamical (MHD) formalism in SPH - i.e. smoothed particle magnetohydrodynamics (SPMHD) - has proven remarkably difficult. To implement MHD requires two steps: a modification to the momentum equation and an induction equation, and both present challenges. We first provide an overview of how SPMHD is implemented, and then discuss how this implementation fails and the limitation of various corrective methods - with particular reference to the effects of particle disorder. Although there are many problems for which, with careful choice of corrective measures, good results can be obtained, we then show that, at the very limits of the state of the art, the ability to perform stable MHD calculations in SPH is curtailed by numerical issues.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Combining radiative transfer and diffuse interstellar medium physics to model star formation

We present a method for modelling star-forming clouds that combines two different models of the thermal evolution of the interstellar medium (ISM). In the combined model, where the densities are low enough that at least some part of the spectrum is optically thin, a model of the thermodynamics of the diffuse ISM is more significant in setting the temperatures. Where the densities are high enough to be optically thick across the spectrum, a model of flux limited diffusion is more appropriate. Previous methods either model the low-density interstellar medium and ignore the thermal behaviour at high densities (e.g. inside collapsing molecular cloud cores), or model the thermal behaviour near protostars but assume a fixed background temperature (e.g. approximately 10 K) on large-scales. Our new method treats both regimes. It also captures the different thermal evolution of the gas, dust, and radiation separately. We compare our results with those from the literature, and investigate the dependence of the thermal behaviour of the gas on the various model parameters. This new method should allow us to model the ISM across a wide range of densities and, thus, develop a more complete and consistent understanding of the role of thermodynamics in the star formation process.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The morphology of the Milky Way - II. Reconstructing CO maps from disc galaxies with live stellar distributions

The arm structure of the Milky Way remains somewhat of an unknown, with observational studies hindered by our location within the Galactic disc. In the work presented here we use smoothed particle hydrodynamics (SPH) and radiative transfer to create synthetic longitude-velocity observations. Our aim is to reverse-engineer a top down map of the Galaxy by comparing synthetic longitude-velocity maps to those observed. We set up a system of N-body particles to represent the disc and bulge, allowing for dynamic creation of spiral features. Interstellar gas, and the molecular content, is evolved alongside the stellar system. A 3D-radiative transfer code is then used to compare the models to observational data. The resulting models display arm features that are a good reproduction of many of the observed emission structures of the Milky Way. These arms however are dynamic and transient, allowing for a wide range of morphologies not possible with standard density wave theory. The best fitting models are a much better match than previous work using fixed potentials. They favour a 4-armed model with a pitch angle of approximately 20 degrees, though with a pattern speed that decreases with increasing Galactic radius. Inner bars are lacking however, which appear required to fully reproduce the central molecular zone.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Spatial Structure of Young Stellar Clusters. III. Physical Properties and Evolutionary States

We analyze the physical properties of stellar clusters that are detected in massive star-forming regions in the MYStIX project--a comparative, multiwavelength study of young stellar clusters within 3.6 kpc that contain at least one O-type star. Tabulated properties of subclusters in these regions include physical sizes and shapes, intrinsic numbers of stars, absorptions by the molecular clouds, and median subcluster ages. Physical signs of dynamical evolution are present in the relations of these properties, including statistically significant correlations between subcluster size, central density, and age, which are likely the result of cluster expansion after gas removal. We argue that many of the subclusters identified in Paper I are gravitationally bound because their radii are significantly less than what would be expected from freely expanding clumps of stars with a typical initial stellar velocity dispersion of ~3 km/s for star-forming regions. We explore a model for cluster formation in which structurally simpler clusters are built up hierarchically through the mergers of subclusters--subcluster mergers are indicated by an inverse relation between the numbers of stars in a subcluster and their central densities (also seen as a density vs. radius relation that is less steep than would be expected from pure expansion). We discuss implications of these effects for the dynamical relaxation of young stellar clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Toroidal vortices and the conglomeration of dust into rings in protoplanetary discs

We identify a new hydrodynamical instability in protoplanetary discs that may arise due to variations in the dust-to-gas ratio and may lead to concentration of dust grains within a disc. The instability can arise due to dust settling, which produces a vertical compositional entropy gradient. The entropy gradient drives a baroclinic instability that is capable of creating toroidal gas vortices that gather dust into rings. Such dust rings are potentially observable via continuum emission of the dust or scattered light. Indeed, this instability may offer an explanation for the rings recently observed in the discs around the young stars HL Tau and TW Hya that does not rely on clearing by protoplanets. The instability may also have wider ramifications, potentially aiding dust agglomeration, altering the radial migration of larger planetesimals, and modifying angular momentum transport within a disc.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Toroidal vortices as a solution to the dust migration problem

In an earlier letter, we reported that dust settling in protoplanetary discs may lead to a dynamical dust-gas instability that produces global toroidal vortices. In this letter, we investigate the evolution of a dusty protoplanetary disc with two different dust species (1 mm and 50 cm dust grains), under the presence of the instability. We show how toroidal vortices, triggered by the interaction of mm grains with the gas, stop the radial migration of metre-sized dust, potentially offering a natural and efficient solution to the dust migration problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-fluid dust and gas mixtures in smoothed particle hydrodynamics II: an improved semi-implicit approach

We present an improved version of the Loren-Aguilar & Bate (2014) method to integrate the two-fluid dust/gas equations that correctly captures the limiting velocity of small grains in the presence of net differences (excluding the drag force) between the accelerations of the dust and the gas. A series of accelerated DUSTYBOX tests and a simulation of dust-settling in a protoplanetary disc are performed comparing the performance of the new and old methods. The modified method can accurately capture the correct limiting velocity while preserving all the conservation properties of the original method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Star Cluster Formation and Feedback

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Origin and Universality of the Stellar Initial Mass Function

We review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Spatial Structure of Young Stellar Clusters. I. Subclusters

The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The statistical properties of stars and their dependence on metallicity: the effects of opacity

We report the statistical properties of stars and brown dwarfs obtained from four radiation hydrodynamical simulations of star cluster formation that resolve masses down to the opacity limit for fragmentation. The calculations are identical except for their dust and gas opacities. Assuming dust opacity is proportional to metallicity, the calculations span a range of metallicities from 1/100 to 3 times solar, although we emphasise that changing the metallicity has other thermodynamic effects that the calculations do not capture (e.g. on the thermal coupling between gas and dust). All four calculations produce stellar populations whose statistical properties are difficult to distinguish from observed stellar systems, and we find no significant dependence of stellar properties on opacity. The mass functions and properties of multiple stellar systems are consistent with each other. However, we find protostellar mergers are more common with lower opacities. Combining the results from the three calculations with the highest opacities, we obtain a stellar population consisting of more than 500 stars and brown dwarfs. Many of the statistical properties of this population are in good agreement with those observed in our Galaxy, implying that gravity, hydrodynamics, and radiative feedback may be the primary ingredients for determining the statistical properties of low-mass stars. However, we do find indications that the calculations may be slightly too dissipative. Although further calculations will be required to understand all of the effects of metallicity on stellar properties, we conclude that stellar properties are surprisingly resilient to variations of the dust and gas opacities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Two fluid dust and gas mixtures in SPH: A semi-implicit approach

A method to avoid the explicit time integration of small dust grains in the two fluid gas/dust smoothed particle hydrodynamics (SPH) approach is proposed. By assuming a very simple exponential decay model for the relative velocity between the gas and dust components, all the effective characteristics of the drag force can be reproduced. A series of tests has been performed to compare the accuracy of the method with analytical and explicit integration results. We find that the method performs well on a wide range of tests, and can provide large speed ups over explicit integration when the dust stopping time is small. We have also found that the method is much less dissipative than conventional explicit or implicit two-fluid SPH approaches when modelling dusty shocks.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the convergence of the critical cooling timescale for the fragmentation of self-gravitating discs

We carry out simulations of gravitationally unstable discs using a Smoothed Particle Hydrodynamics (SPH) code and a grid-based hydrodynamics code, FARGO, to understand the previous non-convergent results reported by Meru & Bate (2011a). We obtain evidence that convergence with increasing resolution occurs with both SPH and FARGO and in both cases we find that the critical cooling timescale is larger than previously thought. We show that SPH has a first-order convergence rate while FARGO converges with a second-order rate. We show that the convergence of the critical cooling timescale for fragmentation depends largely on the numerical viscosity employed in both SPH and FARGO. With SPH, particle velocity dispersion may also play a role. We show that reducing the dissipation from the numerical viscosity leads to larger values of the critical cooling time at a given resolution. For SPH, we find that the effect of the dissipation due to the numerical viscosity is somewhat larger than had previously been appreciated. In particular, we show that using a quadratic term in the SPH artificial viscosity (beta_{SPH}) that is too low appears to lead to excess dissipation in gravitationally unstable discs, which may affect any results that sensitively depend on the thermodynamics, such as disc fragmentation. We show that the two codes converge to values of the critical cooling timescale, beta_{crit} > 20 (for a ratio of specific heats of gamma=5/3), and perhaps even as large as beta_{crit} \approx 30. These are approximately 3-5 times larger than has been found by most previous studies. This is equivalent to a maximum gravitational stress that a disc can withstand without fragmenting of alpha_{GI,crit} \approx 0.013-0.02, which is much smaller than the values typically used in the literature. It is therefore easier for self-gravitating discs to fragment than has been concluded from most past studies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The growth and hydrodynamic collapse of a protoplanet envelope

We have conducted three-dimensional self-gravitating radiation hydrodynamical models of gas accretion onto high mass cores (15-33 Earth masses) over hundreds of orbits. Of these models, one case accretes more than a third of a Jupiter mass of gas, before eventually undergoing a hydrodynamic collapse. This collapse causes the density near the core to increase by more than an order of magnitude, and the outer envelope to evolve into a circumplanetary disc. A small reduction in the mass within the Hill radius (R_H) accompanies this collapse as a shock propagates outwards. This collapse leads to a new hydrostatic equilibrium for the protoplanetary envelope, at which point 97 per cent of the mass contained within the Hill radius is within the inner 0.03 R_H which had previously contained less than 40 per cent. Following this collapse the protoplanet resumes accretion at its prior rate. The net flow of mass towards this dense protoplanet is predominantly from high latitudes, whilst at the outer edge of the circumplanetary disc there is net outflow of gas along the midplane. We also find a turnover of gas deep within the bound envelope that may be caused by the establishment of convection cells.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Collapse of a molecular cloud core to stellar densities: the formation and evolution of pre-stellar discs

We report results from radiation hydrodynamical simulations of the collapse of molecular cloud cores to form protostars. The calculations follow the formation and evolution of the first hydrostatic core/disc, the collapse to form a stellar core, and effect of stellar core formation on the surrounding disc and envelope. Past barotropic calculations have shown that rapidly-rotating first cores evolve into `pre-stellar discs' with radii up to ~100 AU that may last thousands of years before a stellar core forms. We investigate how the inclusion of a realistic equation of state and radiative transfer alters this behaviour, finding that the qualitative behaviour is similar, but that the pre-stellar discs may last 1.5-3 times longer in the more realistic calculations. The masses, radii, and lifetimes of the discs increase for initial molecular cloud cores with faster rotation rates. In the most extreme case we model, a pre-stellar disc with a mass of 0.22 Msun and a radius of ~100 AU can form in a solar-mass cloud and last several thousand years before a stellar core is formed. Such large, massive objects may be imaged using ALMA. Fragmentation of these massive discs may also provide an effective route to binary and multiple star formation, before radiative feedback from accretion onto the stellar core can inhibit fragmentation. Once collapse to form a stellar core occurs within the pre-stellar disc, the radiation hydrodynamical simulations produce qualitatively different behaviour from the barotropic calculations due to the accretion energy released. This drives a shock wave through the circumstellar disc and launches a bipolar outflow even in the absence of magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Migration of protoplanets with surfaces through discs with steep temperature gradients

We perform three-dimensional self-gravitating radiative transfer simulations of protoplanet migration in circumstellar discs to explore the impact upon migration of the radial temperature profiles in these discs. We model protoplanets with masses ranging between 10-100 M\bigoplus, in discs with surface density profiles of Σ \varpropto r^-1/2, and temperature profiles of the form T \varpropto r^-β, where β ranges 0-2. We find that steep (β > 1) temperature profiles lead to outward migration of low mass protoplanets in interstellar grain opacity discs, but in more optically thin discs the migration is always inwards. The trend in migration rates with changing β obtained from our models shows good agreement with those obtained using recent analytic descriptions which include consideration of the co-orbital torques and their saturation. We find that switching between two models of the protoplanet, one in which accretion acts by evacuating gas and one in which gas piles up on a surface to form an atmosphere, leads to a small shift in the migration rates. If comparing these models in discs with conditions which lead to a marginally inward migration, the small shift can lead to outward migration. However, the direction and speed of migration is dominated by disc conditions rather than by the specific prescription used to model the flow near the protoplanet.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stellar, brown dwarf and multiple star properties from a radiation hydrodynamical simulation of star cluster formation

We report the statistical properties of stars, brown dwarfs and multiple systems obtained from the largest radiation hydrodynamical simulation of star cluster formation to date that resolves masses down to the opacity limit for fragmentation (a few Jupiter masses). The initial conditions are identical to those of previous barotropic calculations published by Bate, but this time the calculation is performed using a realistic equation of state and radiation hydrodynamics. The calculation uses sink particles to model 183 stars and brown dwarfs, including 28 binaries and 12 higher-order multiple systems, the properties of which are compared the results from observational surveys. We find that the radiation hydrodynamical/sink particle simulation reproduces many observed stellar properties very well. In particular, whereas using a barotropic equation of state produces more brown dwarfs than stars, the inclusion of radiative feedback results in a stellar mass function and a ratio of brown dwarfs to stars in good agreement with observations of Galactic star-forming regions. In addition, many of the other statistical properties of the stars and brown dwarfs are in reasonable agreement with observations, including multiplicity as a function of primary mass, the frequency of very-low-mass binaries, and general trends for the mass ratio and separation distributions of binaries. We also examine the velocity dispersion of the stars, the distributions of disc truncation radii due to dynamical interactions, and coplanarity of orbits and sink particle spins in multiple systems. Overall, the calculation produces a cluster of stars whose statistical properties are difficult to distinguish from observed systems, implying that gravity, hydrodynamics, and radiative feedback are the primary ingredients for determining the origin of the statistical properties of low-mass stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Collapse of a molecular cloud core to stellar densities: the radiative impact of stellar core formation on the circumstellar disc

We present results from the first three-dimensional radiation hydrodynamical calculations to follow the collapse of a molecular cloud core beyond the formation of the stellar core. We find the energy released by the formation of the stellar core, within the optically-thick first hydrostatic core, is comparable to the binding energy of the disc-like first core. This heats the inner regions of the disc, drives a shock wave through the disc, dramatically decreases the accretion rate on to the stellar core, and launches a temporary bipolar outflow perpendicular to the rotation axis that travels in excess of 50 AU into the infalling envelope. This outburst may assist the young protostar in launching a conventional magnetic jet. Furthermore, if these events are cyclic, they may provide a mechanism for intense bursts of accretion separated by long periods of relatively quiescent accretion which can potentially solve both the protostellar luminosity problem and the apparent age spread of stars in young clusters. Such outbursts may also provide a formation mechanism for the chondrules found in meteorites, with the outflow transporting them to large distances in the circumstellar disc.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring the conditions required to form giant planets via gravitational instability in massive protoplanetary discs

We carry out global three-dimensional radiation hydrodynamical simulations of self-gravitating accretion discs to determine if, and under what conditions, a disc may fragment to form giant planets. We explore the parameter space (in terms of the disc opacity, temperature and size) and include the effect of stellar irradiation. We find that the disc opacity plays a vital role in determining whether a disc fragments. Specifically, opacities that are smaller than interstellar Rosseland mean values promote fragmentation (even at small radii, R < 25AU) since low opacities allow a disc to cool quickly. This may occur if a disc has a low metallicity or if grain growth has occurred. With specific reference to the HR 8799 planetary system, given its star is metal-poor, our results suggest that the formation of its imaged planetary system could potentially have occurred by gravitational instability. We also find that the presence of stellar irradiation generally acts to inhibit fragmentation (since the discs can only cool to the temperature defined by stellar irradiation). However, fragmentation may occur if the irradiation is sufficiently weak that it allows the disc to attain a low Toomre stability parameter.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-convergence of the critical cooling timescale for fragmentation of self-gravitating discs

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results certainly suggest that beta_{crit} is larger than previously thought. However, the absence of convergence also questions whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-convergence of the critical cooling timescale for fragmentation of self-gravitating discs

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results suggest that at the very least, the critical value of the cooling timescale is longer than previously thought. However, the absence of convergence also raises the question of whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the evolution of a star cluster and its multiple stellar systems following gas dispersal

We investigate the evolution, following gas dispersal, of a star cluster produced from a hydrodynamical calculation. We find that when the gas, initially comprising 60% of the mass, is removed, the system settles into a bound cluster containing 30-40% of the stellar mass surrounding by an expanding halo of ejected stars. The bound cluster expands from an initial radius of <0.05 pc to 1-2 pc over 4-10 Myr, depending on how quickly the gas is removed, implying that stellar clusters may begin with far higher stellar densities than usually assumed. With rapid gas dispersal the most massive stars are found to be mass segregated for the first ~1 Myr of evolution, but classical mass segregation only develops for cases with long gas removal timescales. Eventually, many of the most massive stars are expelled from the bound cluster. Despite the high initial stellar density and the extensive dynamical evolution of the system, we find that the stellar multiplicity is almost constant during the 10 Myr of evolution. This is because the primordial multiple systems are formed in a clustered environment and, thus, by their nature are already resistant to further evolution. The majority of multiple system evolution is confined to the decay of high-order systems and the formation of a significant population of very wide (10^4-10^5 AU) multiple systems in the expanding halo. This formation mechanism for wide binaries potentially solves the problem of how most stars apparently form in clusters and yet a substantial population of wide binaries exist in the field. Many of these wide binaries and the binaries produced by the decay of high-order multiple systems have unequal mass components, potentially solving the problem that hydrodynamical simulations of star formation are found to under-produce unequal-mass solar-type binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the fragmentation criteria of self-gravitating protoplanetary discs

We investigate the fragmentation criterion in massive self-gravitating discs. We present new analysis of the fragmentation conditions which we test by carrying out global three-dimensional numerical simulations. Whilst previous work has placed emphasis on the cooling timescale in units of the orbital timescale, β, we find that at a given radius the surface mass density (i.e. disc mass and profile) and star mass also play a crucial role in determining whether a disc fragments or not as well as where in the disc fragments form. We find that for shallow surface mass density profiles (p<2, where Σ\propto R^{-p}), fragments form in the outer regions of the disc. However for steep surface mass density profiles (p \gtrsim 2), fragments form in the inner regions of a disc. In addition, we also find that the critical value of the cooling timescale in units of the orbital timescale found in previous simulations is only applicable to certain disc surface mass density profiles and for particular disc radii and is not a general rule for all discs. We find an empirical fragmentation criteria between the cooling timescale in units of the orbital timescale, β, the surface mass density, the star mass and the radius.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Planet migration: self-gravitating radiation hydrodynamical models of protoplanets with surfaces

We calculate radial migration rates of protoplanets in laminar minimum mass solar nebula discs using three-dimensional self-gravitating radiation hydrodynamical (RHD) models. The protoplanets are free to migrate, whereupon their migration rates are measured. For low mass protoplanets (10-50 M_\oplus) we find increases in the migration timescales of up to an order of magnitude between locally-isothermal and RHD models. In the high-mass regime the migration rates are changed very little. These results are arrived at by calculating migration rates in locally-isothermal models, before sequentially introducing self-gravity, and radiative transfer, allowing us to isolate the effects of the additional physics. We find that using a locally-isothermal equation of state, without self-gravity, we reproduce the migration rates obtained by previous analytic and numerical models. We explore the impact of different protoplanet models, and changes to their assumed radii, upon migration. The introduction of self-gravity gives a slight reduction of the migration rates, whilst the inertial mass problem, which has been proposed for high mass protoplanets with circumplanetary discs, is reproduced. Upon introducing radiative transfer to models of low mass protoplanets (\approx 10 M_\oplus), modelled as small radius accreting point masses, we find outward migration with a rate of approximately twice the analytic inward rate. However, when modelling such a protoplanet in a more realistic manner, with a surface which enables the formation of a deep envelope, this outward migration is not seen.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar and brown dwarf properties from numerical simulations

We review the statistical properties of stars and brown dwarfs obtained from the first hydrodynamical simulation of star cluster formation to produce more than a thousand stars and brown dwarfs while simultaneously resolving the lowest mass brown dwarfs (those with masses set by the opacity limit for fragmentation), binaries with separations down to 1 AU, and discs with radii greater than 10 AU. In particular, we present the eccentricity distribution of the calculation's very-low-mass and brown dwarf binaries which has not been previously published.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The efficiency of star formation in clustered and distributed regions

We investigate the formation of both clustered and distributed populations of young stars in a single molecular cloud. We present a numerical simulation of a 10,000 solar mass elongated, turbulent, molecular cloud and the formation of over 2500 stars. The stars form both in stellar clusters and in a distributed mode which is determined by the local gravitational binding of the cloud. A density gradient along the major axis of the cloud produces bound regions that form stellar clusters and unbound regions that form a more distributed population. The initial mass function also depends on the local gravitational binding of the cloud with bound regions forming full IMFs whereas in the unbound, distributed regions the stellar masses cluster around the local Jeans mass and lack both the high-mass and the low-mass stars. The overall efficiency of star formation is ~ 15 % in the cloud when the calculation is terminated, but varies from less than 1 % in the the regions of distributed star formation to ~ 40 % in regions containing large stellar clusters. Considering that large scale surveys are likely to catch clouds at all evolutionary stages, estimates of the (time-averaged) star formation efficiency for the giant molecular cloud reported here is only ~ 4 %. This would lead to the erroneous conclusion of 'slow' star formation when in fact it is occurring on a dynamical timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three-dimensional molecular line transfer: A simulated star-forming region

We present the first non-LTE, co-moving frame molecular line calculations of a star-forming cluster simulated using smoothed particle hydrodynamics (SPH), from which we derive high-resolution synthetic observations. We have resampled a particle representation onto an adaptive mesh and self-consistently solved the equations of statistical equilibrium in the co-moving frame, using TORUS, a three-dimensional adaptive mesh refined (AMR) radiative transfer (RT) code. We verified the applicability of the code to the conditions of the SPH simulation by testing its output against other codes. We find that the level populations obtained for optically thick and thin scenarios closely match the ensemble average of the other codes. We have used the code to obtain non-LTE level populations of multiple molecular species throughout the cluster and have created three-dimensional velocity-resolved spatial maps of the emergent intensity. Line profiles of cores traced by N2H+ (1-0) are compared to probes of low density gas, 13CO (1-0) and C18O (1-0), surrounding the cores along the line of sight. The relative differences of the line-centre velocities are shown to be small compared to the velocity dispersion, matching recent observations. We conclude that one cannot reject competitive accretion as a viable theory of star formation based on observed velocity profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Circumplanetary disc properties obtained from radiation hydrodynamical simulations of gas accretion by protoplanets

We investigate the properties of circumplanetary discs formed in three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, in solar nebulae of differing grain opacities. We find that the analytical prediction of circumplanetary disc radii in an evacuated gap (R_Hill/3) from Quillen & Trilling (1998) yields a good estimate for discs formed by high mass protoplanets. The radial density profiles of the circumplanetary discs may be described by power-laws between r^-2 and r^-3/2. We find no evidence for the ring-like density enhancements that have been found in some previous models of circumplanetary discs. Temperature profiles follow a ~r^-7/10 power-law regardless of protoplanet mass or nebula grain opacity. The discs invariably have large scaleheights (H/r > 0.2), making them thick in comparison with their encompassing circumstellar discs, and they show no flaring.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The dependence of star formation on initial conditions and molecular cloud structure

We investigate the dependence of stellar properties on the initial kinematic structure of the gas in star-forming molecular clouds. We compare the results from two large-scale hydrodynamical simulations of star cluster formation that resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell and Bromm. The initial conditions of the two calculations are identical, but in the new simulation the power spectrum of the velocity field imposed on the cloud initially and allowed to decay is biased in favour of large-scale motions. Whereas the calculation of Bate et al. began with a power spectrum P(k) ~ k^{-4} to match the Larson scaling relations for the turbulent motions observed in molecular clouds, the new calculation begins with a power spectrum P(k) ~ k^{-6}. Despite this change to the initial motions in the cloud and the resulting density structure of the molecular cloud, the resulting stellar properties resulting from the two calculations are indistinguishable. This demonstrates that the results of such hydrodynamical calculations of star cluster formation are relatively insensitive to the initial conditions. It is also consistent with the fact that the statistical properties of stars and brown dwarfs (e.g. the stellar initial mass function) are observed to be relatively invariant within our Galaxy and do not appear to depend on environment.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Evolution of Giant Planets in Eccentric Disks

We investigate the interaction between a giant planet and a viscous circumstellar disk by means of high-resolution, two-dimensional hydrodynamical simulations. We consider planet masses that range from 1 to 3 Jupiter masses (Mjup) and initial orbital eccentricities that range from 0 to 0.4. We find that a planet can cause eccentricity growth in a disk region adjacent to the planet's orbit, even if the planet's orbit is circular. Disk-planet interactions lead to growth in a planet's orbital eccentricity. The orbital eccentricities of a 2 Mjup and a 3 Mjup planet increase from 0 to 0.11 within about 3000 orbits. Over a similar time period, the orbital eccentricity of a 1 Mjup planet grows from 0 to 0.02. For a case of a 1 Mjup planet with an initial eccentricity of 0.01, the orbital eccentricity grows to 0.09 over 4000 orbits. Radial migration is directed inwards, but slows considerably as a planet's orbit becomes eccentric. If a planet's orbital eccentricity becomes sufficiently large, e > ~0.2, migration can reverse and so be directed outwards. The accretion rate towards a planet depends on both the disk and the planet orbital eccentricity and is pulsed over the orbital period. Planet mass growth rates increase with planet orbital eccentricity. For e~0.2 the mass growth rate of a planet increases by approximately 30% above the value for e=0. For e > ~0.1, most of the accretion within the planet's Roche lobe occurs when the planet is near the apocenter. Similar accretion modulation occurs for flow at the inner disk boundary which represents accretion toward the star.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Star formation in unbound giant molecular clouds: the origin of OB associations?

We investigate the formation of star clusters in an unbound GMC, where the supporting kinetic energy is twice as large as the cloud's self-gravity. This cloud manages to form a series of star clusters and disperse, all within roughly 2 crossing times (10 Myr), supporting recent claims that star formation is a rapid process. Simple assumptions about the nature of the star formation occurring in the clusters allows us to place an estimate for the star formation efficiency at about 5 to 10 %, consistent with observations. We also propose that unbound clouds can act as a mechanism for forming OB associations. The clusters that form in the cloud behave as OB subgroups. These clusters are naturally expanding from one another due to unbound nature of the flows that create them. The properties of the cloud we present here are are consistent with those of classic OB associations.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Origin of the Initial Mass Function and Its Dependence on the Mean Jeans Mass in Molecular Clouds

We investigate the dependence of stellar properties on the mean thermal Jeans mass in molecular clouds. We compare the results from the two largest hydrodynamical simulations of star formation to resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical except for the radii of the clouds, which are chosen so that the mean densities and mean thermal Jeans masses of the clouds differ by factors of nine and three, respectively. We find that the denser cloud, with the lower mean thermal Jeans mass, produces a higher proportion of brown dwarfs and has a lower characteristic (median) mass of the stars and brown dwarfs. This dependence of the initial mass function (IMF) on the density of the cloud may explain the observation that the Taurus star-forming region appears to be deficient in brown dwarfs when compared with the Orion Trapezium cluster. The new calculation also produces wide binaries (separations >20 AU), one of which is a wide binary brown dwarf system. Based on the hydrodynamical calculations, we develop a simple accretion/ejection model for the origin of the IMF. In the model, all stars and brown dwarfs begin with the same mass (set by the opacity limit for fragmentation) and grow in mass until their accretion is terminated stochastically by their ejection from the cloud through dynamically interactions. The model predicts that the main variation of the IMF in different star-forming environments should be in the location of the peak (due to variations in the mean thermal Jeans mass of the cloud) and in the substellar regime. However, the slope of the IMF at high-masses may depend on the dispersion in the accretion rates of protostars.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

The Formation of a Star Cluster: Predicting the Properties of Stars and Brown Dwarfs

We present results from the largest numerical simulation of star formation to resolve the fragmentation process down to the opacity limit. The simulation follows the collapse and fragmentation of a large-scale turbulent molecular cloud to form a stellar cluster and, simultaneously, the formation of circumstellar discs and binary stars. This large range of scales enables us to predict a wide variety of stellar properties for comparison with observations. The calculation clearly demonstrates that star formation is a highly-dynamic and chaotic process. Star-disc encounters form binaries and truncate discs. Stellar encounters disrupt bound multiple systems. The cloud produces roughly equal numbers of stars and brown dwarfs, with masses down to the opacity limit for fragmentation (~5 Jupiter masses). The initial mass function is consistent with a Salpeter slope (Gamma=-1.35) above 0.5 Msun, a roughly flat distribution (Gamma=0) in the range 0.006-0.5 Msun, and a sharp cutoff below ~0.005 Msun. This is consistent with recent observational surveys. The brown dwarfs form by the dynamical ejection of low-mass fragments from dynamically unstable multiple systems before the fragments have been able to accrete to stellar masses. Close binary systems (with separations <10 AU) are not formed by fragmentation in situ. Rather, they are produced by hardening of initially wider multiple systems through a combination of dynamical encounters, gas accretion, and/or the interaction with circumbinary and circumtriple discs. Finally, we find that the majority of circumstellar discs have radii less than 20 AU due to truncation in dynamical encounters. This is consistent with observations of the Orion Trapezium Cluster and implies that most stars and brown dwarfs do not form large planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

The formation of close binary systems

A viable solution to the origin of close binary systems, unaccounted for in recent theories, is presented. Fragmentation, occurring at the end of the secondary collapse phase (during which molecular hydrogen is dissociating), can form binary systems with separations less than 1 au. Two fragmentation modes are found to occur after the collapse is halted. The first consists of the fragmentation of a protostellar disc due to rotational instabilities in a protostellar core, involving both an $m=1$ and an $m=2$ mode. This fragmentation mechanism is found to be insensitive to the initial density distribution: it can occur in both centrally condensed and uniform initial conditions. The second fragmentation mode involves the formation of a rapidly rotating core at the end of the collapse phase which is unstable to the axisymmetric perturbations. This core bounces into a ring which quickly fragments into several components. The binary systems thus formed contain less than 1 per cent of a solar mass and therefore will need to accrete most of their final mass if they are to form a binary star system. Their orbital properties will thus be determined by the properties of the accreted matter.

0 citations0 reviewsNo code linkedOpen access