Source author record

Ian A. Bonnell

Ian A. Bonnell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

22 published item(s)

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.

preprint2020arXiv

Supernova feedback and the energy deposition in molecular clouds

Feedback from supernovae is often invoked as an important process in limiting star formation, removing gas from galaxies and hence as a determining process in galaxy formation. Here we report on numerical simulations investigating the interaction between supernova explosions and the natal molecular cloud. We also consider the cases with and without previous feedback from the high-mass star in the form of ionising radiation and stellar winds. The supernova is able to find weak points in the cloud and create channels through which it can escape, leaving much of the well shielded cloud largely unaffected. This effect is increased when the channels are pre-existing due to the effects of previous stellar feedback. The expanding supernova deposits its energy in the gas that is in these exposed channels, and hence sweeps up less mass when feedback has already occurred, resulting in faster outflows with less radiative losses. The full impact of the supernova explosion is then able to impact the larger scale of the galaxy in which it abides. We conclude that supernova explosions only have moderate effects on their dense natal environments but that with pre-existing feedback, the energetic effects of the supernova are able to escape and affect the wider scale medium of the galaxy.

preprint2013arXiv

Line Profiles of Cores within Clusters: II Signatures of Dynamical Collapse during High Mass Star Formation

Observations of atomic or molecular lines can provide important information about the physical state of star forming regions. In order to investigate the line profiles from dynamical collapsing massive star forming regions (MSFRs), we model the emission from hydrodynamic simulations of a collapsing cloud in the absence of outflows. By performing radiative transfer calculations, we compute the optically thick HCO+ and optically thin N2H+ line profiles from two collapsing regions at different epochs. Due to large-scale collapse, the MSFRs have large velocity gradients, reaching up to 20 km/s/pc across the central core. The optically thin lines typically contain multiple velocity components resulting from the superposition of numerous density peaks along the line-of-sight. The optically thick lines are only marginally shifted to the blue side of the optically thin line profiles, and frequently do not have a central depression in their profiles due to self-absorption. As the regions evolve the lines become brighter and the optically thick lines become broader. The lower order HCO+ (1-0) transitions are better indicators of collapse than the higher order (4-3) transitions. We also investigate how the beam sizes affect profile shapes. Smaller beams lead to brighter and narrower lines that are more skewed to the blue in HCO+ relative to the true core velocity, but show multiple components in N2H+. High resolution observations (e.g. with ALMA) can test these predictions and provide insights into the nature of MSFRs.

preprint2013arXiv

Misaligned streamers around a galactic centre black hole from a single cloud's infall

We follow the near radial infall of a prolate cloud onto a 4 x 10^6 Msun supermassive black hole in the Galactic Centre using smoothed particle hydrodynamics (SPH). We show that a prolate cloud oriented perpendicular to its orbital plane naturally produces a spread in angular momenta in the gas which can translate into misaligned discs as is seen in the young stars orbiting Sagittarius A*. A turbulent or otherwise highly structured cloud is necessary to avoid cancelling too much angular momentum through shocks at closest approach. Our standard model of a 2 x 10^4 Msun gas cloud brought about the formation of a disc within 0.3 pc from the black hole and a larger, misaligned streamer at 0.5 pc. A total of 1.5 x 10^4 Msun of gas formed these structures. Our exploration of the simulation parameter space showed that when star formation occurred, it resulted in top-heavy IMFs with stars on eccentric orbits with semi-major axes 0.02 to 0.3 pc and inclinations following the gas discs and streamers. We suggest that the single event of an infalling prolate cloud can explain the occurrence of multiple misaligned discs of young stars.

preprint2013arXiv

Primordial triples and collisions of massive stars

Massive stars are known to have a high multiplicity, with examples of higher order multiples among the nearest and best studied objects. In this paper we study hierarchical multiple systems (an inner binary as a component of a wider binary) of massive stars in a clustered environment, in which a system with a size of 100--1000 au will undergo many close encounters during the short lifetime of a massive star. Using two types of N-body experiment we determine the post-formation collision probabilities of these massive hierarchies. We find that, depending on the specifics of the environment, the hierarchy, and the amount of time that is allowed to pass, tens of percent of hierarchies will experience a collision, typically between the two stars of the inner binary. In addition to collisions, clusters hosting a hierarchical massive system produce high velocity runaways at an enhanced rate. The primordial multiplicity specifics of massive stars appear to play a key role in the generation of these relatively small number events in cluster simulations, complicating their use as diagnostics of a cluster's history.

preprint2013arXiv

Shocks, cooling and the origin of star formation rates in spiral galaxies

Understanding star formation is problematic as it originates in the large scale dynamics of a galaxy but occurs on the small scale of an individual star forming event. This paper presents the first numerical simulations to resolve the star formation process on sub-parsec scales, whilst also following the dynamics of the interstellar medium (ISM) on galactic scales. In these models, the warm low density ISM gas flows into the spiral arms where orbit crowding produces the shock formation of dense clouds, held together temporarily by their external pressure. Cooling allows the gas to be compressed to sufficiently high densities that local regions collapse under their own gravity and form stars. The star formation rates follow a Schmidt-Kennicutt Σ_{SFR} ~ Σ_{gas}^{1.4} type relation with the local surface density of gas while following a linear relation with the cold and dense gas. Cooling is the primary driver of star formation and the star formation rates as it determines the amount of cold gas available for gravitational collapse. The star formation rates found in the simulations are offset to higher values relative to the extragalactic values, implying a constant reduction, such as from feedback or magnetic fields, is likely to be required. Intriguingly, it appears that a spiral or other convergent shock and the accompanying thermal instability can explain how star formation is triggered, generate the physical conditions of molecular clouds and explain why star formation rates are tightly correlated to the gas properties of galaxies.

preprint2012arXiv

How long does it take to form a molecular cloud?

We present the first numerical simulations that self-consistently follow the formation of dense molecular clouds in colliding flows. Our calculations include a time-dependent model for the H2 and CO chemistry that runs alongside a detailed treatment of the dominant heating and cooling processes in the ISM. We adopt initial conditions characteristic of the warm neutral medium and study two different flow velocities - a slow flow with v = 6.8 km/s and a fast flow with v = 13.6 km/s. The clouds formed by the collision of these flows form stars, with star formation beginning after 16 Myr in the case of the slower flow, but after only 4.4 Myr in the case of the faster flow. In both flows, the formation of CO-dominated regions occurs only around 2 Myr before the onset of star formation. Prior to this, the clouds produce very little emission in the J = 1 -> 0 transition line of CO, and would probably not be identified as molecular clouds in observational surveys. In contrast, our models show that H2-dominated regions can form much earlier, with the timing depending on the details of the flow. In the case of the slow flow, small pockets of gas become fully molecular around 10 Myr before star formation begins, while in the fast flow, the first H2-dominated regions occur around 3 Myr before the first prestellar cores form. Our results are consistent with models of molecular cloud formation in which the clouds are dominated by "dark" molecular gas for a considerable proportion of their assembly history.

preprint2012arXiv

The rapid dispersal of low-mass virialised clusters

Infant mortality brought about by the expulsion of a star cluster's natal gas is widely invoked to explain cluster statistics at different ages. While a well studied problem, most recent studies of gas expulsion's effect on a cluster have focused on massive clusters, with stellar counts of order $10^4$. Here we argue that the evolutionary timescales associated with the compact low-mass clusters typical of the median cluster in the Solar neighborhood are short enough that significant dynamical evolution can take place over the ages usually associated with gas expulsion. To test this we perform {\it N}-body simulations of the dynamics of a very young star forming region, with initial conditions drawn from a large-scale hydrodynamic simulation of gravitational collapse and fragmentation. The subclusters we analyse, with populations of a few hundred stars, have high local star formation efficiencies and are roughly virialised even after the gas is removed. Over 10 Myr they expand to a similar degree as would be expected from gas expulsion if they were initially gas-rich, but the expansion is purely due to the internal stellar dynamics of the young clusters. The expansion is such that the stellar densities at 2 Myr match those of YSOs in the Solar neighborhood. We argue that at the low-mass end of the cluster mass spectrum, a deficit of clusters at 10s of Myr does not necessarily imply gas expulsion as a disruption mechanism.

preprint2011arXiv

The dynamical state of stellar structure in star-forming regions

The fraction of star formation that results in bound star clusters is influenced by the density spectrum in which stars are formed and by the response of the stellar structure to gas expulsion. We analyse hydrodynamical simulations of turbulent fragmentation in star-forming regions to assess the dynamical properties of the resulting population of stars and (sub)clusters. Stellar subclusters are identified using a minimum spanning tree algorithm. When considering only the gravitational potential of the stars and ignoring the gas, we find that the identified subclusters are close to virial equilibrium (the typical virial ratio Q_vir~0.59, where virial equilibrium would be Q_vir~0.5). This virial state is a consequence of the low gas fractions within the subclusters, caused by the accretion of gas onto the stars and the accretion-induced shrinkage of the subclusters. Because the subclusters are gas-poor, up to a length scale of 0.1-0.2 pc at the end of the simulation, they are only weakly affected by gas expulsion. The fraction of subclusters that reaches the high density required to evolve to a gas-poor state increases with the density of the star-forming region. We extend this argument to star cluster scales, and suggest that the absence of gas indicates that the early disruption of star clusters due to gas expulsion (infant mortality) plays a smaller role than anticipated, and is potentially restricted to star-forming regions with low ambient gas densities. We propose that in dense star-forming regions, the tidal shocking of young star clusters by the surrounding gas clouds could be responsible for the early disruption. This `cruel cradle effect' would work in addition to disruption by gas expulsion. We suggest possible methods to quantify the relative contributions of both mechanisms.

preprint2010arXiv

A quantification of the non-spherical geometry and accretion of collapsing cores

We present the first detailed classification of the structures of Class 0 cores in a high resolution simulation of a giant molecular cloud. The simulated cloud contains 10^4 solar masses and produces over 350 cores which allows for meaningful statistics. Cores are classified into three types according to how much they depart from spherical symmetry. We find that three quarters of the cores are better described as irregular filaments than as spheres. Recent Herschel results have shown that cores are formed within a network of filaments, which we find has had a significant impact on the resulting core geometries. We show that the column densities and ram pressure seen by the protostar are not uniform and generally peak along the axes of the filament. The angular momentum vector of the material in the cores varies both in magnitude and direction, which will cause the rotation vector of the central source to fluctuate during the collapse of the core. In the case of the more massive stars, accretion from the environment outside the original core volume is even more important than that from the core itself. This additional gas is primarily accreted onto the cores along the dense filaments in which the cores are embedded, and the sections of the surfaces of the cores which do not coincide with a filament have very little additional material passing through them. The assumption of spherical symmetry cannot be applied to the majority of collapsing cores, and is never a good description of how stars accrete gas from outside the original core radius. This has ramifications for our understanding of collapsing cores, in particular their line profiles, the effect of radiation upon them and their ability to fragment.

preprint2010arXiv

Escaping stars from young low-N clusters

With the use of N-body calculations the amount and properties of escaping stars from low-N (N = 100 and 1000) young embedded star clusters prior to gas expulsion are studied over the first 5 Myr of their existence. Besides the number of stars also different initial radii and binary populations are examined as well as virialised and collapsing clusters. It is found that these clusters can loose substantial amounts (up to 20%) of stars within 5 Myr with considerable velocities up to more than 100 km/s. Even with their mean velocities between 2 and 8 km/s these stars will still be travelling between 2 and 30 pc during the 5 Myr. Therefore can large amounts of distributed stars in star-forming regions not necessarily be counted as evidence for the isolated formation of stars.

preprint2010arXiv

Super-star clusters versus OB associations

Super Star Clusters (Mecl > 10^5 Msol) are the largest stellar nurseries in our local Universe, containing hundreds of thousands to millions of young stars within a few light years. Many of these systems are found in external galaxies, especially in pairs of interacting galaxies, and in some dwarf galaxies, but relatively few in disk galaxies like our own Milky Way. We show that a possible explanation for this difference is the presence of shear in normal spiral galaxies which impedes the formation of the very large and dense super star clusters but prefers the formation of loose OB associations possibly with a less massive cluster at the center. In contrast, in interacting galaxies and in dwarf galaxies, regions can collapse without having a large-scale sense of rotation. This lack of rotational support allows the giant clouds of gas and stars to concentrate into a single, dense and gravitationally bound system.

preprint2010arXiv

The Effect of Environment on Massive Star Formation

In this contribution we review our recent numerical work discussing the essential role of the local cluster environment in assembling massive stars. First we show that massive stars are formed from low mass pre-stellar cores and become massive due to accretion. Proto-stars that benefit from this accretion are those situated at the centre of a cluster's potential well, which is the focal point of the contraction of the cluster gas. Given that most of the mass which makes up a massive star in this model comes from the cluster environment rather than the core, it is important to model the molecular cloud environment accurately. Preliminary results of a simulation which accurately treats the chemistry and time-dependent thermodynamics of a molecular cloud show quantitatively similar star formation to previous models, but allow a true comparison to be made between simulation and observations. This method can also be applied to cases with varying metallicities allowing star formation in primordial gas to be studied. In general, these numerical studies of clustered star formation yield IMFs which are compatible with the Salpeter mass function. The only possible exception to this is in low density unbound regions of molecular clouds which lack very low and high mass stars.

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.

preprint2009arXiv

Does sub-cluster merging accelerate mass segregation in local star formation?

The nearest site of massive star formation in Orion is dominated by the Trapezium subsystem, with its four OB stars and numerous companions. The question of how these stars came to be in such close proximity has implications for our understanding of massive star formation and early cluster evolution. A promising route toward rapid mass segregation was proposed by McMillan et al. (2007), who showed that the merger product of faster-evolving sub clusters can inherit their apparent dynamical age from their progenitors. In this paper we briefly consider this process at a size and time scale more suited for local and perhaps more typical star formation, with stellar numbers from the hundreds to thousands. We find that for reasonable ages and cluster sizes, the merger of sub-clusters can indeed lead to compact configurations of the most massive stars, a signal seen both in Nature and in large-scale hydrodynamic simulations of star formation from collapsing molecular clouds, and that sub-virial initial conditions can make an un-merged cluster display a similar type of mass segregation. Additionally, we discuss a variation of the minimum spanning tree mass-segregation technique introduced by Allison et al. (2009).

preprint2009arXiv

Fragmentation in Molecular Clouds and its connection to the IMF

We present an analysis of star-forming gas cores in an SPH simulation of a Giant Molecular Cloud. We identify cores using their deep potential wells. This yields a smoother distribution with clearer boundaries than density. Additionally, this gives an indication of future collapse, as bound potential cores (p-cores) represent the earliest stages of fragmentation in molecular clouds. We find that the mass function of the p-cores resembles the stellar IMF and the observed clump mass function, although p-core masses (~0.7 Msol) are smaller than typical density clumps. The bound p-cores are generally subsonic, have internal substructure, and are only quasi-spherical. We see no evidence of massive bound cores supported by turbulence. We trace the evolution of the p-cores forward in time, and investigate the connection between the original p-core mass and the stellar mass that formed from it. We find that there is a poor correlation, with considerable scatter suggesting accretion onto the core is dependent on more factors than just the initial core mass. During the accretion process the p-cores accrete from beyond the region first bound, highlighting the importance of the core environment to its subsequent evolution.

preprint2009arXiv

Limits on initial mass segregation in young clusters

Mass segregation is observed in many star clusters, including several that are less than a few Myr old. Timescale arguments are frequently used to argue that these clusters must be displaying primordial segregation, because they are too young to be dynamically relaxed. Looking at this argument from the other side, the youth of these clusters and the limited time available to mix spatially distinct populations of stars can provide constraints on the amount of initial segregation that is consistent with current observations. We present n-body experiments testing this idea, and discuss the implications of our results for theories of star formation. For system ages less than a few crossing times, we show that star formation scenarios predicting general primordial mass segregation are inconsistent with observed segregation levels.

preprint2005arXiv

Competitive Accretion and the IMF

Competitive accretion occurs when stars in a cluster accrete from a shared reservoir of gas. The competition arises due to the relative attraction of stars as a function of their mass and location in the cluster. The low relative motions of the stars and gas in young, gas dominated clusters results in a tidal limit to the accretion whereas in the stellar dominated cluster cores, the high relative velocities results in Bondi-Hoyle accretion. The combination of these two accretion processes produces a two power-law IMF with $γ\approx -1.5$, for low-mass stars which accrue their mass in the gas dominated regime, and a steeper, $γ\approx -2.5$, IMF for higher-mass stars that form in the core of a cluster. Simulations of the fragmentation and formation of a stellar cluster show that the final stellar masses, and IMF, are due to competitive accretion. Competitive accretion also naturally results in a mass segregated cluster and in a direct correlation between the richness of a cluster and the mass of the most massive star therein. The {\sl knee} where the IMF slope changes occurs near the Jeans mass of the system.

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.

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.

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.

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.