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

A. P. Whitworth

A. P. Whitworth 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-phastro-ph.IMastro-ph.COastro-ph.EP

Catalog footprint

What is connected

35works
6topics
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

35 published item(s)

preprint2022arXiv

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Protostellar Outflows: a window to the past

During the early phases of low-mass star formation, episodic accretion causes the ejection of high-velocity outflow bullets, which carry a fossil record of the driving protostar's accretion history. We present 44 SPH simulations of $1\,\mathrm{M}_{\odot}$ cores, covering a wide range of initial conditions, and follow the cores for five free-fall times. Individual protostars are represented by sink particles, and the sink particles launch episodic outflows using a subgrid model. The Optics algorithm is used to identify individual episodic bullets within the outflows. The parameters of the overall outflow and the individual bullets are then used to estimate the age and energetics of the outflow, and the accretion events that triggered it; and to evaluate how reliable these estimates are, if observational uncertainties and selection effects (like inclination) are neglected. Of the commonly used methods for estimating outflow ages, it appears that those based on the length and speed of advance of the lobe are the most reliable in the early phases of evolution, and those based on the width of the outflow cavity and the speed of advance are most reliable during the later phases. We describe a new method that is almost as accurate as these methods, and reliable throughout the evolution. In addition we show how the accretion history of the protostar can be accurately reconstructed from the dynamics of the bullets if each lobe contains at least two bullets. The outflows entrain about ten times more mass than originally ejected by the protostar.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The widths of magnetised filaments in molecular clouds

Filaments are an ubiquitous feature of molecular clouds, and appear to play a critical role in assembling the material to form stars. The dominant filaments are observed to have a rather narrow range of widths around $\sim 0.1$ pc, and to be preferentially aligned perpendicularly to the direction of the local magnetic field. We have previously argued that the observed filament widths can be explained if filaments are formed by converging, mildly supersonic flows, resulting from large-scale turbulent motions in the parent molecular cloud. Here we demonstrate that the introduction of a magnetic field perpendicular to the filament long axis does not greatly alter this conclusion, as long as the mass-to-flux ratio is supercritical. The distribution of widths for supercritical magnetised filaments formed via this mechanism is peaked at slightly higher values, and is slightly broader, than for non-magnetised filaments, but still reproduces the basic properties of the width distributions derived from far-infrared observations of molecular clouds. In contrast, subcritical filaments have width distributions with a fundamentally different shape, and typically have much larger widths than those observed. Both subcritical and supercritical filaments are consistent with the observed lack of correlation between filament widths and filament surface densities.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterising lognormal fractional-Brownian-motion density fields with a Convolutional Neural Network

In attempting to quantify statistically the density structure of the interstellar medium, astronomers have considered a variety of fractal models. Here we argue that, to properly characterise a fractal model, one needs to define precisely the algorithm used to generate the density field, and to specify -- at least -- three parameters: one parameter constrains the spatial structure of the field; one parameter constrains the density contrast between structures on different scales; and one parameter constrains the dynamic range of spatial scales over which self-similarity is expected (either due to physical considerations, or due to the limitations of the observational or numerical technique generating the input data). A realistic fractal field must also be noisy and non-periodic. We illustrate this with the exponentiated fractional Brownian motion (xfBm) algorithm, which is popular because it delivers an approximately lognormal density field, and for which the three parameters are, respectively, the power spectrum exponent, $β$, the exponentiating factor, ${\cal S}$, and the dynamic range, ${\cal R}$. We then explore and compare two approaches that might be used to estimate these parameters: Machine Learning and the established $Δ$-Variance procedure. We show that for $2\leqβ\leq 4$ and $0\leq{\cal S}\leq 3$, a suitably trained Convolutional Neural Network is able to estimate objectively both $β$ (with root-mean-square error $ε_{_β}\sim 0.12$) and ${\cal S}$ (with $ε_{_{\cal S}}\sim 0.29$). $\;Δ$-variance is also able to estimate $β$, albeit with a somewhat larger error ($ε_{_β}\sim 0.17$) and with some human intervention, but is not able to estimate ${\cal S}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The apparent anti-correlation between the mass opacity of interstellar dust and the surface-density of interstellar gas

Recent analyses of ${\it Herschel}$ observations suggest that in nearby disc galaxies the dust mass opacity at $500 \, {\rm μm}$, $κ_{500}$, decreases with increasing gas surface density, $Σ_{\rm ISM}$ (Clark et al. 2019). This apparent anti-correlation between $κ_{500}$ and $Σ_{\rm ISM}$ is opposite to the behaviour expected from theoretical dust evolution models; in such models, dust in denser, cooler regions (i.e. regions of increased $Σ_{\rm ISM}$) tends to grow and therefore to have increased $κ_{500}$. We show, using a toy model, that the presence of a range of dust temperatures along the line of sight can lead to spuriously low estimated values of $κ_{500}$. If in regions of higher $Σ_{\rm ISM}$ the range of dust temperatures extends to lower values (as seems likely), the magnitude of this effect may be sufficient to explain the apparent anti-correlation between $κ_{500}$ and $Σ_{\rm ISM}$. Therefore there may not be any need for spatial variation in the intrinsic dust properties that run counter to theoretical expectations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Hi-GAL catalogue of dusty filamentary structures in the Galactic Plane

The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structure of the Galaxy. In order to study this potential link, a statistically significant sample of filaments spread throughout the Galaxy is required. In this work we present the first catalogue of $32,059$ candidate filaments automatically identified in the Hi-GAL survey of the entire Galactic Plane. For these objects we determined morphological (length, $l^{a}$, and geometrical shape) and physical (average column density, $N_{\rm H_{2}}$, and average temperature, $T$) properties. We identified filaments with a wide range of properties: 2$'$\,$\leq l^{a}\leq$\, 100$'$, $10^{20} \leq N_{\rm H_{2}} \leq 10^{23}$\,cm$^{-2}$ and $10 \leq T\leq$ 35\,K. We discuss their association with the Hi-GAL compact sources, finding that the most tenuous (and stable) structures do not host any major condensation and we also assign a distance to $\sim 18,400$ filaments for which we determine mass, physical size, stability conditions and Galactic distribution. When compared to the spiral arms structure, we find no significant difference between the physical properties of on-arm and inter-arm filaments. We compared our sample with previous studies, finding that our Hi-GAL filament catalogue represents a significant extension in terms of Galactic coverage and sensitivity. This catalogue represents an unique and important tool for future studies devoted to understanding the filament life-cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A ram-pressure threshold for star formation

In turbulent fragmentation, star formation occurs in condensations created by converging flows. The condensations must be sufficiently massive, dense and cool to be gravitationally unstable, so that they start to contract; {\it and} they must then radiate away thermal energy fast enough for self-gravity to remain dominant, so that they continue to contract. For the metallicities and temperatures in local star forming clouds, this second requirement is only met robustly when the gas couples thermally to the dust, because this delivers the capacity to radiate across the full bandwidth of the continuum, rather than just in a few discrete spectral lines. This translates into a threshold for vigorous star formation, which can be written as a minimum ram-pressure Pcrit ~ 4 10^-11 dyne. Pcrit is independent of temperature, and corresponds to flows with molecular hydrogen number-density nH2 and velocity v satisfying nH2 v^2 > 800 cm^-3 (km/s)^2. This in turn corresponds to a minimum molecular hydrogen column-density for vigorous star formation, NH2crit ~ 4 10^21 cm^-2 (SIGMAcrit ~ 100 MSun pc^-2), and a minimum visual extinction AVcrit ~ 9. The characteristic diameter and line-density for a star-forming filament when this threshold is just exceeded -- a sweet spot for local star formation regions -- are 2Rfil ~ 0.1 pc and mufil ~ 13 MSun pc^-2. The characteristic diameter and mass for a prestellar core condensing out of such a filament are 2Rcore ~ 0.1 pc, and Mcore ~ MSun. We also show that fragmentation of a shock-compressed layer is likely to commence while the convergent flows creating the layer are still ongoing, and we stress that, under this circumstance, the phenomenology and characteristic scales for fragmentation of the layer are fundamentally different from those derived traditionally for pre-existing layers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A theoretical perspective on the formation and fragmentation of protostellar discs

We discuss the factors influencing the formation and gravitational fragmentation of protostellar discs. We start with a review of how observations of prestellar cores can be analysed statistically to yield plausible initial conditions for simulations of their subsequent collapse. Simulations based on these initial conditions show that, despite the low levels of turbulence in prestellar cores, they deliver primary protostars and associated discs which are routinely subject to stochastic impulsive perturbations; consequently misalignment of the spins and orbits of protostars are common. Also, the simulations produce protostars that collectively have a mass function and binary statistics matching those observed in nearby star formation regions, but only if a significant fraction of the turbulent energy in the core is solenoidal, and accretion onto the primary protostar is episodic with a duty cycle > 3000 yr. Under this circumstance a core typically spawns between 4 and 5 protostars, with high efficiency, and the lower-mass protostars are mainly formed by disc fragmentation. The requirement that a proto-fragment in a disc lose thermal energy on a dynamical timescale dictates that {\bf there is a sweet spot for} disc fragmentation at radii 70 AU < R < 100 AU and temperatures 10 K < T < 20 K, and this might explain the Brown Dwarf Desert.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence that widespread star formation may be underway in G0.253+016, "The Brick"

Image cubes of differential column density as a function of dust temperature are constructed for Galactic Centre molecular cloud G0.253+0.016 ("The Brick") using the recently described PPMAP procedure. The input data consist of continuum images from the Herschel Space Telescope in the wavelength range 70-500 $μ$m, supplemented by previously published interferometric data at 1.3 mm wavelength. While the bulk of the dust in the molecular cloud is consistent with being heated externally by the local interstellar radiation field, our image cube shows the presence, near one edge of the cloud, of a filamentary structure whose temperature profile suggests internal heating. The structure appears as a cool ($\sim 14$ K) tadpole-like feature, $\sim 6$ pc in length, in which is embedded a thin spine of much hotter ($\sim$ 40-50 K) material. We interpret these findings in terms of a cool filament whose hot central region is undergoing gravitational collapse and fragmentation to form a line of protostars. If confirmed, this would represent the first evidence of widespread star formation having started within this cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Forming isolated brown dwarfs by turbulent fragmentation

We use Smoothed Particle Hydrodynamics to explore the circumstances under which an isolated very-low-mass prestellar core can be formed by colliding turbulent flows and collapse to form a brown-dwarf. Our simulations suggest that the flows need not be very fast, but do need to be very strongly convergent, i.e. the gas must flow in at comparable speeds from all sides, which seems rather unlikely. We therefore revisit the object Oph-B11, which Andre, Ward-Thompson and Greaves (2012) have identified as a prestellar core with mass between $\sim 0.020\,\mathrm{M_\odot}$ and $\sim 0.030\,\mathrm{M_\odot}$. We reanalyse the observations using a Markov-chain Monte Carlo method that allows us (i) to include the uncertainties on the distance, temperature and dust mass opacity, and (ii) to consider different Bayesian prior distributions of the mass. We estimate that the posterior probability that Oph-B11 has a mass below the hydrogen burning limit at $\sim 0.075\,\mathrm{M_\odot}$, is between 0.66 and 0.86. We conclude that, if Oph-B11 is destined to collapse, it probably will form a brown dwarf. However, the flows required to trigger this appear to be so contrived that it is difficult to envisage this being the only way, or even a major way, of forming isolated brown dwarfs. Moreover, Oph-B11 could easily be a transient, bouncing, prolate core, seen end-on; there could, indeed should, be many such objects masquerading as very low-mass prestellar cores.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SPAMCART: a code for smoothed particle Monte Carlo radiative transfer

We present a code for generating synthetic SEDs and intensity maps from Smoothed Particle Hydrodynamics simulation snapshots. The code is based on the Lucy (1999) Monte Carlo Radiative Transfer method, i.e. it follows discrete luminosity packets, emitted from external and/or embedded sources, as they propagate through a density field, and then uses their trajectories to compute the radiative equilibrium temperature of the ambient dust. The density is not mapped onto a grid, and therefore the calculation is performed at exactly the same resolution as the hydrodynamics. We present two example calculations using this method. First, we demonstrate that the code strictly adheres to Kirchhoff's law of radiation. Second, we present synthetic intensity maps and spectra of an embedded protostellar multiple system. The algorithm uses data structures that are already constructed for other purposes in modern particle codes. It is therefore relatively simple to implement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Are the majority of Sun-like stars single?

It has recently been suggested that, in the field, $\sim\!\!56\%$ of Sun-like stars ($0.8\,{\rm M}_{_\odot}\lesssim M_\star\lesssim 1.2\,{\rm M}_{_\odot}$) are single. We argue here that this suggestion may be incorrect, since it appears to be based on the multiplicity frequency of systems with Sun-like primaries, and therefore takes no account of Sun-like stars that are secondary (or higher-order) components in multiple systems. When these components are included in the reckoning, it seems likely that only $\sim\!46\%$ of Sun-like stars are single. This estimate is based on a model in which the system mass function has the form proposed by Chabrier, with a power-law Salpeter extension to high masses; there is a flat distribution of mass ratios; and the probability that a system of mass $M$ is a binary is $\,0.50 + 0.46\log_{_{10}}\!\left(M/{\rm M}_{_\odot}\right)\,$ for $\,0.08\,{\rm M}_{_\odot}\leq M\leq 12.5\,{\rm M}_{_\odot}$, $\,0\,$ for $\,M<0.08\,{\rm M}_{_\odot}$, and $\,1\,$ for $\,M>12.5\,{\rm M}_{_\odot}$. The constants in this last relation are chosen so that the model also reproduces the observed variation of multiplicity frequency with primary mass. However, the more qualitative conclusion, that a minority of Sun-like stars are single, holds up for virtually all reasonable values of the model parameters. Parenthetically, it is still likely that the majority of {\it all} stars in the field are single, but that is because most M Dwarfs probably are single.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the effects of solenoidal and compressive turbulence in prestellar cores

We present the results of an ensemble of SPH simulations that follow the evolution of prestellar cores for $0.2\,{\rm Myr}$. All the cores have the same mass, and start with the same radius, density profile, thermal and turbulent energy. Our purpose is to explore the consequences of varying the fraction of turbulent energy, $δ_\mathrm{sol}$, that is solenoidal, as opposed to compressive; specifically we consider $δ_\mathrm{sol}=1,\,2/3,\,1/3,\,1/9\;{\rm and}\;0$. For each value of $δ_\mathrm{sol}$, we follow ten different realisations of the turbulent velocity field, in order also to have a measure of the stochastic variance blurring any systematic trends. With low $δ_\mathrm{sol}(<\!1/3)$ filament fragmentation dominates and delivers relatively high mass stars. Conversely, with high values of $δ_\mathrm{sol}(>\!1/3)$ disc fragmentation dominates and delivers relatively low mass stars. There are no discernible systematic trends in the multiplicity statistics obtained with different $δ_\mathrm{sol}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Star Formation triggered by cloud-cloud collisions

We present the results of SPH simulations in which two clouds, each having mass $M_{_{\rm{o}}}\!=\!500\,{\rm M}_{_\odot}$ and radius $R_{_{\rm{o}}}\!=\!2\,{\rm pc}$, collide head-on at relative velocities of $Δv_{_{\rm{o}}} =2.4,\;2.8,\;3.2,\;3.6\;{\rm and}\;4.0\,{\rm km}\,{\rm s}^{-1}$. There is a clear trend with increasing $Δv_{_{\rm{o}}}$. At low $Δv_{_{\rm{o}}}$, star formation starts later, and the shock-compressed layer breaks up into an array of predominantly radial filaments; stars condense out of these filaments and fall, together with residual gas, towards the centre of the layer, to form a single large-$N$ cluster, which then evolves by competitive accretion, producing one or two very massive protostars and a diaspora of ejected (mainly low-mass) protostars; the pattern of filaments is reminiscent of the hub and spokes systems identified recently by observers. At high $Δv_{_{\rm{o}}}$, star formation occurs sooner and the shock-compressed layer breaks up into a network of filaments; the pattern of filaments here is more like a spider's web, with several small-$N$ clusters forming independently of one another, in cores at the intersections of filaments, and since each core only spawns a small number of protostars, there are fewer ejections of protostars. As the relative velocity is increased, the {\it mean} protostellar mass increases, but the {\it maximum} protostellar mass and the width of the mass function both decrease. We use a Minimal Spanning Tree to analyse the spatial distributions of protostars formed at different relative velocities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Temperature as a third dimension in column-density mapping of dusty astrophysical structures associated with star formation

We present PPMAP, a Bayesian procedure that uses images of dust continuum emission at multiple wavelengths to produce resolution-enhanced image cubes of differential column-density as a function of dust temperature and position. PPMAP is based on the generic 'point process' formalism, whereby the system of interest (in this case, a dusty astrophysical structure such as a filament or prestellar core) is represented by a collection of points in a suitably defined state space. It can be applied to a variety of observational data, such as Herschel images, provided only that the image intensity is delivered by optically thin dust in thermal equilibrium. PPMAP takes full account of the instrumental point spread functions and does not require all images to be degraded to the same resolution. We present the results of testing using simulated data for a prestellar core and a fractal turbulent cloud, and demonstrate its performance with real data from the Hi-GAL survey. Specifically, we analyse observations of a large filamentary structure in the CMa OB1 giant molecular cloud. Histograms of differential column-density indicate that the warm material (T > 13 K) is distributed log-normally, consistent with turbulence, but the column-densities of the cooler material are distributed as a high density tail, consistent with the effects of self-gravity. The results illustrate the potential of PPMAP to aid in distinguishing between different physical components along the line of sight in star-forming clouds, and aid the interpretation of the associated PDFs of column density.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Role of Discs in the Collapse and Fragmentation of Prestellar Cores

Disc fragmentation provides an important mechanism for producing low mass stars in prestellar cores. Here, we describe Smoothed Particle Hydrodynamics simulations which show how populations of prestellar cores evolve into stars. We find the observed masses and multiplicities of stars can be recovered under certain conditions. First, protostellar feedback from a star must be episodic. The continuous accretion of disc material on to a central protostar results in local temperatures which are too high for disc fragmentation. If, however, the accretion occurs in intense outbursts, separated by a downtime of $\sim10^4\,\mathrm{years}$, gravitational instabilities can develop and the disc can fragment. Second, a significant amount of the cores' internal kinetic energy should be in solenoidal turbulent modes. Cores with less than a third of their kinetic energy in solenoidal modes have insufficient angular momentum to form fragmenting discs. In the absence of discs, cores can fragment but results in a top heavy distribution of masses with very few low mass objects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating star formation in Ophiuchus

We have simulated star formation in prestellar cores, using SPH and initial conditions informed by observations of the cores in Ophiuchus. Because the observations are limited to two spatial dimensions plus radial velocity, we cannot infer initial conditions for the collapse of a particular core. However, with a minimum of assumptions (isotropic turbulence with a power-law spectrum, a thermal mix of compressive and solenoidal modes, a critical Bonnor-Ebert density profile) we can generate initial conditions that match, in a statistical sense, the distributions of mass, projected size and aspect ratio, thermal and non-thermal one-dimensional velocity dispersion, observed in Ophiuchus. The time between core-core collisions in Ophiuchus is sufficiently long, that we can simulate single cores evolving is isolation, and therefore we are able to resolve masses well below the opacity limit. We generate an ensemble of 100 cores, and evolve them with no radiative feedback from the stars formed, then with continuous radiative feedback, and finally with episodic radiative feedback. With no feedback the simulations produce too many brown dwarfs, and with continuous feedback too few. With episodic radiative feedback, both the peak of the protostellar mass function (at ~ 0.2 M_sun) and the ratio of H-burning stars to brown dwarfs are consistent with observations. The mass of a star is not strongly related to the mass of the core in which it forms. Low-mass cores (M ~ 0.1 M_sun) tend to collapse into single objects, whereas high-mass cores (1 > M_sun) usually fragment into several objects with a broad mass range.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulations of star formation in Ophiuchus, II: Multiplicity

Lomax et al. have constructed an ensemble of 60 prestellar cores having masses, sizes, projected shapes, temperatures and non-thermal radial velocity dispersions that match, statistically, the cores in Ophiuchus; and have simulated the evolution of these cores using SPH. Each core has been evolved once with no radiative feedback from stars, once with continuous radiative feedback, and once with episodic radiative feedback. Here we analyse the multiplicity statistics from these simulations. With episodic radiative feedback, (i) the multiplicity frequency is ~60% higher than in the field; (ii) the multiplicity frequency and the mean semi-major axis both increase with primary mass; (iii) one third of multiple systems are hierarchical systems with more than two components; (iv) in these hierarchical systems the inner pairings typically have separations of a few au and mass ratios concentrated towards unity, whereas the outer pairings have separations of order 100 au and a flatter distribution of mass ratios. The binary statistics are compatible with observations of young embedded populations, and -- if wider orbits are disrupted preferentially by external perturbations -- with observations of mature field populations. With no radiative feedback, the results are similar to those from simulations with episodic feedback. With continuous radiative feedback, brown dwarfs are under-produced, the number of multiple systems is too low, and the statistical properties of multiple systems are at variance with observation. This suggests that star formation in Ophiuchus may only be representative of global star formation if accretion onto protostars, and hence radiative feedback, is episodic.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An improved sink particle algorithm for SPH simulations

Numerical simulations of star formation frequently rely on the implementation of sink particles, (a) to avoid expending computational resource on the detailed internal physics of individual collapsing protostars, (b) to derive mass functions, binary statistics and clustering kinematics (and hence to make comparisons with observation), and (c) to model radiative and mechanical feedback; sink particles are also used in other contexts, for example to represent accreting black holes in galactic nuclei. We present a new algorithm for creating and evolving sink particles in SPH simulations, which appears to represent a significant improvement over existing algorithms {\refrpt -- particularly in situations where sinks are introduced after the gas has become optically thick to its own cooling radiation and started to heat up by adiabatic compression}. (i) It avoids spurious creation of sinks. (ii) It regulates the accretion of matter onto a sink so as to mitigate non-physical perturbations in the vicinity of the sink. (iii) Sinks accrete matter, but the associated angular momentum is transferred back to the surrounding medium. With the new algorithm -- and modulo the need to invoke sufficient resolution to capture the physics preceding sink formation -- the properties of sinks formed in simulations are essentially independent of the user-defined parameters of sink creation, or the number of SPH particles used.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Clumps and triggered star formation in ionised molecular clouds

Infrared shells and bubbles are ubiquitous in the Galaxy and can generally be associated with HII regions formed around young, massive stars. In this paper, we use high-resolution 3D SPH simulations to explore the effect of a single O7 star emitting photons at 10^49 1/s and located at the centre of a molecular cloud with mass 10^4 M_sun and radius 6.4 pc; the internal structure of the cloud is characterised by its fractal dimension, D (with 2.0 <= D <= 2.8), and the variance of its (log-normal) density distribution, sigma_0^2 (with 0.36 <= sigma_0^2 <= 1.42). Our study focuses on the morphology of the swept-up cold gas and the distribution and statistics of the resulting star formation. If the fractal dimension is low, the border of the HII region is dominated by extended shell-like structures, and these break up into a small number of massive high-density clumps which then spawn star clusters; star formation occurs relatively quickly, and delivers somewhat higher stellar masses. Conversely, if the fractal dimension is high, the border of the HII region is dominated by a large number of pillars and cometary globules, which contain compact dense clumps and tend to spawn single stars or individual multiple systems; star formation occurs later, the stellar masses are somewhat lower, and the stars are more widely distributed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mapping the core mass function onto the stellar IMF: multiplicity matters

Observations indicate that the central portions of the Present-Day Prestellar Core Mass Function (CMF) and the Stellar Initial Mass Function (IMF) both have approximately log-normal shapes, but that the CMF is displaced to higher mass than the IMF by a factor F = 4+/-1. This has lead to suggestions that the shape of the IMF is directly inherited from the shape of the CMF - and therefore, by implication, that there is a self-similar mapping from the CMF onto the IMF. If we assume a self-similar mapping, it follows (i) that F = N0/eta, where eta is the mean fraction of a core's mass that ends up in stars, and N0 is the mean number of stars spawned by a single core; and (ii) that the stars spawned by a single core must have an approximately log-normal distribution of relative masses, with universal standard deviation sigma0. Observations can be expected to deliver ever more accurate estimates of F, but this still leaves a degeneracy between eta and N0; and sigma0 is also unconstrained by observation. Here we show that these parameters can be estimated by invoking binary statistics. Specifically, if (a) each core spawns one long-lived binary system, and (b) the probability that a star of mass M is part of this long-lived binary is proportional to M^alpha, current observations of the binary frequency as a function of primary mass and the distribution of mass ratios strongly favour eta = 1.0 +/- 0.3, N0 = 4.3 +/- 0.4, sigma0 = 0.3 +/- 0.03 and alpha = 0.9 +/- 0.6. The mapping from CMF to IMF is not necessarily self-similar - there are many possible motivations for a non self-similar mapping - but if it is not, then the shape of the IMF cannot be inherited from the CMF. Given the limited observational constraints currently available and the ability of a self-similar mapping to satisfy them, the possibility that the shape of the IMF is inherited from the CMF cannot be ruled out at this juncture.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of the non-linear thin shell instability

We use three-dimensional smoothed particle hydrodynamics simulations to study the non-linear thin shell instability (NTSI) in supersonic colliding flows. We show that for flows with monochromatic perturbations and for flows with white-noise perturbations, growth speeds approximate quite well to the analytic predictions of Vishniac (1994). For flows with subsonic turbulence, growth speeds match Vishniac's predictions only at short wavelengths where the turbulence is weaker. We find that supersonic turbulence, of a lower Mach number than the colliding flows, completely suppresses the NTSI. Our results provide a diagnostic for identifying the presence of the NTSI in colliding flows with turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dispersal of molecular clouds by ionising radiation

The role of feedback from massive stars is believed to be a key element in the evolution of molecular clouds. We use high-resolution 3D SPH simulations to explore the dynamical effects of a single O7 star located at the centre of a molecular cloud with mass 10^4M_sun and radius 6.4pc. The initial internal structure of the cloud is characterised by its fractal dimension, D=2.0 - 2.8, and its log-normal density PDF. (i) As regards star formation, in the short term ionising feedback is positive, in the sense that star formation occurs much more quickly in gas that is compressed by the high pressure of the ionised gas. However, in the long term ionising feedback is negative, in the sense that most of the cloud is dispersed with an outflow rate of up to ~0.01M_sun/yr, on a timescale comparable with the sound-crossing time for the ionised gas (~1-2Myr), and triggered star formation is therefore limited to a few percent of the cloud's mass. (ii) As regards the morphology of the ionisation fronts (IFs) bounding the HII region and the systematics of outflowing gas, we distinguish two regimes. For low D<=2.2, the initial cloud is dominated by large-scale structures, so the neutral gas tends to be swept up into a few extended coherent shells, and the ionised gas blows out through a few large holes between these shells; we term these HII regions "shell-dominated". Conversely, for high D>=2.6, the initial cloud is dominated by small-scale structures, and these are quickly overrun by the advancing IF, thereby producing neutral pillars whilst the ionised gas blows out through a large number of small holes between the pillars; we term these HII regions "pillar-dominated". (iii) As regards the injection of bulk kinetic energy, by ~1Myr, the expansion of the HII region has delivered a rms velocity of ~6km/s; this represents less than 0.1% of the total energy radiated by the O7 star.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Interactions between brown-dwarf binaries and Sun-like stars

Several mechanisms have been proposed for the formation of brown dwarfs, but there is as yet no consensus as to which -- if any -- are operative in nature. Any theory of brown dwarf formation must explain the observed statistics of brown dwarfs. These statistics are limited by selection effects, but they are becoming increasingly discriminating. In particular, it appears (a) that brown dwarfs that are secondaries to Sun-like stars tend to be on wide orbits, $a\ga 100\,{\rm AU}$ (the Brown Dwarf Desert), and (b) that these brown dwarfs have a significantly higher chance of being in a close ($a\la 10\,{\rm AU}$) binary system with another brown dwarf than do brown dwarfs in the field. This then raises the issue of whether these brown dwarfs have formed {\it in situ}, i.e. by fragmentation of a circumstellar disc; or have formed elsewhere and subsequently been captured. We present numerical simulations of the purely gravitational interaction between a close brown-dwarf binary and a Sun-like star. These simulations demonstrate that such interactions have a negligible chance ($<0.001$) of leading to the close brown-dwarf binary being captured by the Sun-like star. Making the interactions dissipative by invoking the hydrodynamic effects of attendant discs might alter this conclusion. However, in order to explain the above statistics, this dissipation would have to favour the capture of brown-dwarf binaries over single brown-dwarfs, and we present arguments why this is unlikely. The simplest inference is that most brown-dwarf binaries -- and therefore possibly also most single brown dwarfs -- form by fragmentation of circumstellar discs around Sun-like protostars, with some of them subsequently being ejected into the field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SEREN - A new SPH code for star and planet formation simulations

We present SEREN, a new hybrid Smoothed Particle Hydrodynamics and N-body code designed to simulate astrophysical processes such as star and planet formation. It is written in Fortran 95/2003 and has been parallelised using OpenMP. SEREN is designed in a flexible, modular style, thereby allowing a large number of options to be selected or disabled easily and without compromising performance. SEREN uses the conservative `grad-h' formulation of SPH, but can easily be configured to use traditional SPH or Godunov SPH. Thermal physics is treated either with a barotropic equation of state, or by solving the energy equation and modelling the transport of cooling radiation. A Barnes-Hut tree is used to obtain neighbour lists and compute gravitational accelerations efficiently, and an hierarchical time-stepping scheme is used to reduce the number of computations per timestep. Dense gravitationally bound objects are replaced by sink particles, to allow the simulation to be evolved longer, and to facilitate the identification of protostars and the compilation of stellar and binary properties. At the termination of a hydrodynamical simulation, SEREN has the option of switching to a pure N-body simulation, using a 4th-order Hermite integrator, and following the ballistic evolution of the sink particles (e.g. to determine the final binary statistics once a star cluster has relaxed). We describe in detail all the algorithms implemented in SEREN and we present the results of a suite of tests designed to demonstrate the fidelity of SEREN and its performance and scalability. Further information and additional tests of SEREN can be found at the web-page http://www.astro.group.shef.ac.uk/seren.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The influence of the turbulent perturbation scale on prestellar core fragmentation and disk formation

The collapse of weakly turbulent prestellar cores is a critical stage in the process of star formation. Being highly non-linear and stochastic, the outcome of collapse can only be explored theoretically by performing large ensembles of numerical simulations. Standard practice is to quantify the initial turbulent velocity field in a core in terms of the amount of turbulent energy (or some equivalent) and the exponent in the power spectrum (n \equiv -d log Pk /d log k). In this paper, we present a numerical study of the influence of the details of the turbulent velocity field on the collapse of an isolated, weakly turbulent, low-mass prestellar core. We show that, as long as n > 3 (as is usually assumed), a more critical parameter than n is the maximum wavelength in the turbulent velocity field, λ_MAX. This is because λ_MAX carries most of the turbulent energy, and thereby influences both the amount and the spatial coherence of the angular momentum in the core. We show that the formation of dense filaments during collapse depends critically on λ_MAX, and we explain this finding using a force balance analysis. We also show that the core only has a high probability of fragmenting if λ_MAX > 0.5 R_CORE (where R_CORE is the core radius); that the dominant mode of fragmentation involves the formation and break-up of filaments; and that, although small protostellar disks (with radius R_DISK <= 20 AU) form routinely, more extended disks are rare. In turbulent, low-mass cores of the type we simulate here, the formation of large, fragmenting protostellar disks is suppressed by early fragmentation in the filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The initial conditions of isolated star formation - X. A suggested evolutionary diagram for prestellar cores

We propose an evolutionary path for prestellar cores on the radius-mass diagram, which is analogous to stellar evolutionary paths on the Hertzsprung-Russell Diagram. Using James Clerk Maxwell Telescope (JCMT) observations of L1688 in the Ophiuchus star-forming complex, we analyse the HCO+ (J=4\rightarrow3) spectral line profiles of prestellar cores. We find that of the 58 cores observed, 14 show signs of infall in the form of a blue-asymmetric double-peaked line profile. These 14 cores all lie beyond the Jeans mass line for the region on a radius-mass plot. Furthermore another 10 cores showing tentative signs of infall, in their spectral line profile shapes, appear on or just over the Jeans mass line. We therefore propose the manner in which a prestellar core evolves across this diagram. We hypothesise that a core is formed in the low-mass, low-radius region of the plot. It then accretes quasistatically, increasing in both mass and radius. When it crosses the limit of gravitational instability it begins to collapse, decreasing in radius, towards the region of the diagram where protostellar cores are seen.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical comparison of clouds and star clusters

The extent to which the projected distribution of stars in a cluster is due to a large-scale radial gradient, and the extent to which it is due to fractal sub-structure, can be quantified -- statistically -- using the measure ${\cal Q} = \bar{m}/\bar{s}$. Here $\bar{m}$ is the normalized mean edge length of its minimum spanning tree (i.e. the shortest network of edges connecting all stars in the cluster) and $\bar{s}$ is the correlation length (i.e. the normalized mean separation between all pairs of stars). We show how ${\cal Q}$ can be indirectly applied to grey-scale images by decomposing the image into a distribution of points from which $\bar{m}$ and $\bar{s}$ can be calculated. This provides a powerful technique for comparing the distribution of dense gas in a molecular cloud with the distribution of the stars that condense out of it. We illustrate the application of this technique by comparing ${\cal Q}$ values from simulated clouds and star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The direction of outflows from filaments : constraints on core formation

Context: It is generally presumed that the outflows from a YSO are directed close to its rotation axis (i.e. along its angular momentum vector and orthogonal to any attendant accretion disc). Many YSOs are formed from dense prestellar cores embedded in filaments, and therefore the relative orientations of outflows and filaments may place useful constraints on the dynamics of core formation. Aims: We explore this possibility, from the viewpoint of what it may tell us about the angular momentum delivered to a core forming in a filament. We stress that we are not here addressing the issue of the relationship of filaments and outflows to the prevailing magnetic field direction, although this is evidently also an interesting issue. Results: Using the limited data available, and neglecting any possible selection effects or correlations between nearby outflows, we infer that the observed outflows have a tendency to be orthogonal to the filaments that contain their driving sources. Specfically, in the cases where the directions of the filaments and outflows are well defined, we infer statistically that 72% of outflows are within $45^{\rm o}$ of being orthogonal to the filament, and only 28% are within $45^{\rm o}$ of being parallel to the filament. Conclusions: This suggests that the prestellar cores which spawned the YSOs driving the observed outflows had angular momenta which were approximately orthogonal to the filaments out of which the cores formed. We briefly discuss the implications of this for two proposed core formation mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Adaptive smoothing lengths in SPH

Context: There is a need to improve the fidelity of SPH simulations of self-gravitating gas dynamics. Aims: We remind users of SPH that, if smoothing lengths are adjusted so as to keep the number of neighbours, ${\cal N}$, in the range ${\cal N}_{_{\rm NEIB}}\pmΔ{\cal N}_{_{\rm NEIB}}$, the tolerance, $Δ{\cal N}_{_{\rm NEIB}}$, should be set to zero, as first noted by Nelson & Papaloizou. We point out that this is a very straightforward and computationally inexpensive constraint to implement. Methods: We demonstrate this by simulating acoustic oscillations of a self-gravitating isentropic monatomic gas-sphere (cf. Lucy), using ${\cal N}_{_{\rm TOT}}\sim6,000$ particles and ${\cal N}_{_{\rm NEIB}}=50$. Results: We show that there is a marked reduction in the rates of numerical dissipation and diffusion as $Δ{\cal N}_{_{\rm NEIB}}$ is reduced from 10 to zero. Moreover this reduction incurs a very small computational overhead. Conclusions: We propose that this should become a standard test for codes used in simulating star formation. It is a highly relevant test, because pressure waves generated by the switch from approximate isothermality to approximate adiabaticity play a critical role in the fragmentation of collapsing prestellar cores. Since many SPH simulations in the literature use ${\cal N}_{_{\rm NEIB}}=50$ and $Δ{\cal N}_{_{\rm NEIB}}\geq10$, their results must be viewed with caution.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Radiative Transfer in Prestellar Cores: A Monte Carlo Approach

We use our Monte Carlo radiative transfer code to study non-embedded prestellar cores and cores that are embedded at the centre of a molecular cloud. Our study indicates that the temperature inside embedded cores is lower than in isolated non-embedded cores, and generally less than 12 K, even when the cores are surrounded by an ambient cloud of small visual extinction (Av~5). Our study shows that the best wavelength region to observe embedded cores is between 400 and 500 microns, where the core is quite distinct from the background. We also predict that very sensitive observations (~1-3 MJy/sr) at 170-200 microns can be used to estimate how deeply a core is embedded in its parent molecular cloud. Finally, we present preliminary results of asymmetric models of non-embedded cores.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Numerical simulations of protostellar encounters I. Star-disc encounters

It appears that most stars are born in clusters, and that at birth most stars have circumstellar discs which are comparable in size to the separations between the stars. Interactions between neighbouring stars and discs are therefore likely to play a key role in determining disc lifetimes, stellar masses, and the separations and eccentricities of binary orbits. Such interactions may also cause fragmentation of the discs, thereby triggering the formation of additional stars. We have carried out a series of simulations of disc-star interactions using an SPH code which treats self-gravity, hydrodynamic and viscous forces. We find that interactions between discs and stars provide a mechanism for removing energy from, or adding energy to, the orbits of the stars, and for truncating the discs. However, capture during such encounters is unlikely to be an important binary formation mechanism. A more significant consequence of such encounters is that they can trigger fragmentation of the disc, via tidally and compressionally induced gravitational instabilities, leading to the formation of additional stars. When the disc-spins and stellar orbits are randomly oriented, encounters lead to the formation of new companions to the original star in 20% of encounters. If most encounters are prograde and coplanar, as suggested by simulations of dynamically-triggered star formation, then new companions are formed in approximately 50% of encounters.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Numerical simulations of protostellar encounters II. Coplanar disc-disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. Such interactions must have a significant impact upon the evolution of the protostars and their discs. We have carried out a series of simulations of coplanar encounters between two stars, each possessing a massive circumstellar disc, using an SPH code that models gravitational, hydrodynamic and viscous forces. We find that during a coplanar encounter, disc material is swept up into a shock layer between the two interacting stars, and the layer then fragments to produce new protostellar condensations. The truncated remains of the discs may subsequently fragment; and the outer regions of the discs may be thrown off to form circumbinary disc-like structures around the stars. Thus coplanar disc-disc encounters lead efficiently to the formation of multiple star systems and small-N clusters.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Numerical simulations of protostellar encounters III. Non-coplanar disc-disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. If protostars are formed individually within a cluster before falling together and interacting, there should be no preferred orientation for such interactions. As star formation within clusters is believed to be coeval, it is probable that during interactions, both protostars possess massive, extended discs. We have used an SPH code to carry out a series of simulations of non-colpanar disc-disc interactions. We find that non-coplanar interactions trigger gravitational instabilities in the discs, which may then fragment to form new companions to the existing stars. (This is different from coplanar interactions, in which most of the new companion stars form after material in the discs has been swept up into a shock layer, and this then fragments.) The original stars may also capture each other, leading to the formation of a small-N cluster. If every star undergoes a randomly oriented disc-disc interaction, then the outcome will be the birth of many new stars. Approximately two-thirds of the stars will end up in multiple systems.

0 citations0 reviewsNo code linkedOpen access