Hostname: page-component-7c8c6479df-8mjnm Total loading time: 0 Render date: 2024-03-29T14:59:19.475Z Has data issue: false hasContentIssue false

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

Published online by Cambridge University Press:  22 January 2016

O. Lomax*
Affiliation:
School of Physics and Astronomy, Cardiff University, Cardiff CF24 3AA, UK
A. P. Whitworth
Affiliation:
School of Physics and Astronomy, Cardiff University, Cardiff CF24 3AA, UK
D. A. Hubber
Affiliation:
University Observatory, Ludwig-Maximilians-University Munich, Scheinerstr.1, D-81679 Munich, Germany Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching, Germany
Rights & Permissions [Opens in a new window]

Abstract

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 ~ 104 yr, 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.

Type
Review Article
Copyright
Copyright © Astronomical Society of Australia 2016 

1 INTRODUCTION

Two of the main goals of star formation theory are (i) to understand the origin of the stellar initial mass function (IMF; e.g. Kroupa Reference Kroupa2001; Chabrier Reference Chabrier2003, Reference Chabrier, Corbelli, Palla and Zinnecker2005) and (ii) to explain the properties of stellar multiple systems (e.g. Raghavan et al. Reference Raghavan2010; Janson et al. Reference Janson2012). One possible solution to this problem is the turbulent fragmentation of giant molecular clouds. Here, turbulent flows within molecular clouds produce dense cores of gas (e.g. Padoan & Nordlund Reference Padoan and Nordlund2002; Hennebelle & Chabrier Reference Hennebelle and Chabrier2008, Reference Hennebelle and Chabrier2009) of order 0.01 to 0.1 pc across. These may be Jeans unstable, in which case they collapse to form stars (e.g. Andre, Ward-Thompson, & Barsony Reference Andre, Ward-Thompson and Barsony1993, Reference Andre, Ward-Thompson and Barsony2000). This has been demonstrated in numerical simulations by Bate (Reference Bate1998, Reference Bate2000), Horton, Bate, & Bonnell (Reference Horton, Bate and Bonnell2001), Matsumoto & Hanawa (Reference Matsumoto and Hanawa2003), Goodwin & Whitworth (Reference Goodwin and Whitworth2004), Delgado-Donate, Clarke, & Bate (Reference Delgado-Donate, Clarke and Bate2004a), Delgado-Donate et al. (Reference Delgado-Donate, Clarke, Bate and Hodgkin2004b), Goodwin, Whitworth, & Ward-Thompson (Reference Goodwin, Whitworth and Ward-Thompson2004, Reference Goodwin, Whitworth and Ward-Thompson2006), Walch et al. (Reference Walch, Burkert, Whitworth, Naab and Gritschneder2009, Reference Walch, Naab, Whitworth, Burkert and Gritschneder2010), Walch, Whitworth, & Girichidis (Reference Walch, Whitworth and Girichidis2012), Lomax et al. (Reference Lomax, Whitworth, Hubber, Stamatellos and Walch2014, Reference Lomax, Whitworth, Hubber, Stamatellos and Walch2015b), and Lomax, Whitworth, & Hubber (Reference Lomax, Whitworth and Hubber2015a).

Observations of prestellar cores (e.g. Motte, Andre, & Neri Reference Motte, Andre and Neri1998; Testi & Sargent Reference Testi and Sargent1998; Johnstone et al. Reference Johnstone, Wilson, Moriarty-Schieven, Joncas, Smith, Gregersen and Fich2000; Motte et al. Reference Motte, André, Ward-Thompson and Bontemps2001; Johnstone et al. Reference Johnstone, Fich, Mitchell and Moriarty-Schieven2001; Stanke et al. Reference Stanke, Smith, Gredel and Khanzadyan2006; Enoch et al. Reference Enoch2006; Johnstone & Bally Reference Johnstone and Bally2006; Nutter & Ward-Thompson Reference Nutter and Ward-Thompson2007; Alves, Lombardi, & Lada Reference Alves, Lombardi and Lada2007; Enoch et al. Reference Enoch, Evans, Sargent, Glenn, Rosolowsky and Myers2008; Simpson, Nutter, & Ward-Thompson Reference Simpson, Nutter and Ward-Thompson2008; Rathborne et al. Reference Rathborne, Lada, Muench, Alves, Kainulainen and Lombardi2009; Könyves et al., Reference Könyves2010; Pattle et al. Reference Pattle2015) show that the core mass function (CMF) is very similar in shape to the IMF (i.e. a lognormal distribution with a power-law tail at high mass), albeit shifted upwards in mass by a factor of 3 to 5. This had led to the suggestion that there is a self-similar mapping of the CMF onto the IMF. Statistical analysis by Holman et al. (Reference Holman, Walch, Goodwin and Whitworth2013) suggests that a core should spawn on average four to five stars in order to explain the observed abundance of multiple systems.

A core may fragment into multiple objects via either turbulent fragmentation—similar to the molecular cloud—or disc fragmentation. Observed N2H+ line widths in cores indicate that the internal velocity dispersion in most cores is sub to trans-sonic. This suggests that a typical core is unlikely to collapse and fragment into more than one or two objects through turbulence alone. However, the first protostars to form are usually attended by accretion discs (Kenyon & Hartmann Reference Kenyon and Hartmann1995). These discs may fragment if two criteria are fulfilled. First, the disc must have a sufficiently large surface density Σ(R) so that fragments can overcome thermal and centrifugal support,

(1) $$\begin{eqnarray} \Sigma (R)\gtrsim \frac{c(R)\,\kappa (R)}{\pi \,G}, \end{eqnarray}$$

where κ(R) is the epicyclic frequency and c(R) is the sound speed (Toomre Reference Toomre1964). Second, the cooling time of a fragment must be shorter than the orbital period if it is to avoid being sheared part (Gammie Reference Gammie2001).

The above criteria apply well to low-mass discs. However, disc dynamics are more complex when the disc and the central protostar have similar masses. In cases where the disc is marginally unstable (i.e. $t_\rm {COOL}\sim t_\rm {ORBIT}$ and π G Σ ~ c κ), instabilities develop which bolster the accretion rate onto the central protostar. This lowers the mass (and hence surface density) of the disc, restabilising it. This suggests that an otherwise unstable disc may be able to remain stable by undergoing episodic accretion events (Lodato & Rice Reference Lodato and Rice2005).

When non-local effects such as radiative transfer dominate the disc temperature, local cooling timescales become very difficult—if not impossible—to calculate. Here, again, a seemingly unstable disc may be able to resist fragmentation (e.g. Tsukamoto et al. Reference Tsukamoto, Takahashi, Machida and Inutsuka2015). Forgan & Rice (Reference Forgan and Rice2013) also show that irradiation (from a central protostar and/or interstellar radiation field) increases the jeans mass within disc spiral structures. They deduce that disc fragmentation is far more likely to result in low-mass stars and brown dwarfs than gas giant planets.

In this paper, we analyse simulations which follow the evolution of an ensemble of synthetic cores, based on the properties of those in the Ophiuchus star-forming region (See Lomax et al. Reference Lomax, Whitworth, Hubber, Stamatellos and Walch2014, Reference Lomax, Whitworth and Hubber2015a, Reference Lomax, Whitworth, Hubber, Stamatellos and Walch2015b, for more details). We examine how protostellar feedback affects the fragmentation of discs, and how the results compare with observed stars. We also examine how the ratio of solenoidal to compressive modes in the turbulent velocity field affects the formation of discs and filaments within a prestellar core. In Section 2, we describe the numerical method used to evolve the cores. In Section 3, we describe (i) how to generate realistic core initial conditions and (ii) how they evolve under different prescriptions of radiative feedback. In Section 4, we show how changing the ratio of solenoidal to compressive modes in the velocity field affects core fragmentation. We summarise and conclude in Section 5.

2 NUMERICAL METHOD

2.1. Smoothed particle hydrodynamics

Core evolution is simulated using the serenh-SPH code (Hubber et al. Reference Hubber, Batty, McLeod and Whitworth2011). Gravitational forces are computed using a tree and artificial viscosity is controlled by the Morris & Monaghan (Reference Morris and Monaghan1997) prescription. In all simulations, the smoothed particle hydrodynamics (SPH) particles have mass $m_\rm {SPH}=10^{-5}\,{\rm M}_{\odot }$ , so that the opacity limit (~ 3 × 10−3 M) is resolved with ~ 300 particles. Gravitationally bound regions with density higher than $\rho _\rm {SINK}=10^{-9}\,{\rm g}\,{\rm cm}^{-3}$ are replaced with sink particles (Hubber, Walch, & Whitworth Reference Hubber, Walch and Whitworth2013). Sink particles have radius $r_\rm {SINK}\simeq 0.2\,{\rm au}$ , corresponding to the smoothing length of an SPH particle with density equal to $\rho _\rm {SINK}$ . The equation of state and the energy equation are treated with the algorithm described in Stamatellos et al. (Reference Stamatellos, Whitworth, Bisbas and Goodwin2007b). Magenetic fields and mechanical feedback (e.g. stellar winds) is not included in these simulations.

2.2. Accretion feedback

Radiative feedback from the protostars (i.e. sink particles formed in the simulations) is included. The dominant contribution to the luminosity of a protostar is usually from accretion,

(2) $$\begin{eqnarray} L_\star \simeq \frac{f\,G\,M_\star \,\dot{M}_\star }{R_\star }, \end{eqnarray}$$

where, f = 0.75 is the fraction of the accreted material’s gravitational energy that is radiated from the surface of the protostar (the rest is presumed to be removed by bipolar jets and outflows; Offner et al. Reference Offner, Klein, McKee and Krumholz2009), M is the mass of the protostar, $\dot{M}_\star$ is the rate of accretion onto the protostar, and R = 3 R is the approximate radius of a protostar (Palla & Stahler Reference Palla and Stahler1993).

We adopt the phenomenological model of episodic accretion, presented by Stamatellos, Whitworth, & Hubber (Reference Stamatellos, Whitworth and Hubber2011, Reference Stamatellos, Whitworth and Hubber2012). This is based on calculations of the magneto-rotational instability (MRI; Zhu, Hartmann, & Gammie Reference Zhu, Hartmann and Gammie2009, Reference Zhu, Hartmann and Gammie2010a; Zhu et al. Reference Zhu, Hartmann, Gammie, Book, Simon and Engelhard2010b). In the outer disc of a protostar (outside the sink radius, and therefore resolved by the simulation), angular momentum is redistributed by gravitational torques and material spirals inwards towards the sink. At distances within the sink radius (unresolved by the simulation), the inner disc is so hot that it is gravitationally stable, and unable to fragment. Material continues to accrete onto this inner disc (i.e. the sink particle) until the gas is hot enough to thermally ionise and couple with the local magnetic field. The MRI cuts in and magnetic torques allow material in the inner disc to rapidly accrete onto the protostar. This results in extended periods of very low accretion luminosity, punctuated by intense, episodic outbursts. The length of the downtime between outbursts—during which disc fragmentation may occur—is given by

(3) $$\begin{eqnarray} && t_\rm {DOWN}\sim 1.3\times 10^4\text{\,yr}\nonumber\\ &&\quad\times \left(\frac{M_\star }{0.2\,{\rm M}_{\odot }}\right)^{2/3}\left(\frac{\dot{M}_\rm {SINK}}{10^{-5}\,{\rm M}_{\odot }{\rm yr}^{-1}}\right)^{-8/9}, \end{eqnarray}$$

where $\dot{M}_\rm {SINK}$ is the rate at which material flows into the sink.

There is also observational motivation for adopting an episodic model. The luminosities of young stars are about an order of magnitude lower than expected from continuous accretion (this is the luminosity problem, first noted by Kenyon et al. Reference Kenyon, Hartmann, Strom and Strom1990). Furthermore, FU Ori-type stars (e.g. Herbig Reference Herbig1977; Greene, Aspin, & Reipurth Reference Greene, Aspin and Reipurth2008; Peneva et al. Reference Peneva, Semkov, Munari and Birkle2010; Green et al. Reference Green2011; Principe et al. Reference Principe, Kastner, Grosso, Hamaguchi, Richmond, Teets and Weintraub2013) can exhibit large increases in luminosity which last ≲ 102 yr. Statistical arguments by Scholz, Froebrich, & Wood (Reference Scholz, Froebrich and Wood2013) suggest that the downtime between outbursts should be of order 104 yr, similar to the timescale given in Equation (3).

3 PRESTELLAR CORES IN OPHIUCHUS

3.1. Initial conditions

Using the observed properties of cores as a basis for numerical simulations presents a difficult inverse problem. The mass, temperature, projected area, and projected aspect ratio of a core can be reasonably inferred from bolometric measurements of the dust in prestellar cores. In addition, the line-of-sight velocity dispersion within the core can be inferred from the width of molecular lines. However, the initial boundary conditions of a core simulation must represent the full spatial and velocity structure of the system. This occupies six dimensions, whereas observational data can only provide information on three (i.e. two spatial and one velocity).

Rather than trying to emulate a specific core—which is arguably an impossibility—we can relatively easily define a distribution of cores which have the same, or at least very similar, statistical properties to those in a given region. We based the synthetic core initial conditions on Ophiuchus. This is a well-studied region, for which many of the aforementioned core properties have been measured.

3.1.1. Mass, size, and velocity dispersion

Only some of the measured core masses in Ophiuchus have both an associated size and velocity dispersion. In order to make the most of the data, we define the following lognormal probability distribution of ${\bf{x}}\equiv (\log (M),\log (R),\log (\sigma _\rm {NT}))$ :

(4) $$\begin{equation} P(\bf{ x})=\frac{1}{(2\pi )^{3/2}|\bf{ \Sigma }|}\exp \left(-\frac{1}{2}(\bf{ x}-\bf{ \mu })^{\rm T}\bf{ \Sigma }^{-1}(\bf{ x}-\bf{ \mu }) \right), \end{equation}$$

where

(5) $$\begin{equation} \bf{ \mu }\equiv \left(\begin{array}{c} \mu _{M} \\ \mu _{R} \\ \mu _{{\sigma _\rm {NT}}} \end{array}\right), \end{equation}$$

and

(6) $$\begin{equation} \bf{ \Sigma }\equiv \left( \begin{array}{c@{\hspace*{8pt}}c@{\hspace*{8pt}}c} \sigma _{M}^2 & \rho _{{M,R}}\,\sigma _{M}\sigma _{R} & \rho _{{M,\sigma _\rm {NT}}}\,\sigma _{M}\sigma _{{\sigma _\rm {NT}}} \\ \rho _{{M,R}}\,\sigma _{M}\sigma _{R} & \sigma _{R}^2 & \rho _{{R,\sigma _\rm {NT}}}\,\sigma _{R}\sigma _{{\sigma _\rm {NT}}} \\ \rho _{{M,\sigma _\rm {NT}}}\,\sigma _{M}\sigma _{{\sigma _\rm {NT}}} & \rho _{{R,\sigma _\rm {NT}}}\,\sigma _{R}\sigma _{{\sigma _\rm {NT}}} & \sigma _{{\sigma _\rm {NT}}}^2 \end{array}\right). \end{equation}$$

The coefficients of μ and Σ are calculated from the observed Ophiuchus data (Motte et al. Reference Motte, Andre and Neri1998; André et al. Reference André, Belloche, Motte and Peretto2007) and are given in Table 1. From P( x ), we are able to draw any number of masses, sizes, and velocity dispersions, all of which were statistically similar to those in Ophiuchus. The distribution of P(x) is shown in Figure 1.

Table 1. Arithmetic means, standard deviations, and correlation coefficients of log (M), log (R), and $\log (\sigma _\rm {NT })$ for cores in Ophiuchus.

Figure 1. The multivariate lognormal distribution, P( x ) where $\bf{ x}=(\log (M),\log (R),\log (\sigma _\rm {NT}))$ . The figure shows the projections through $\log (\sigma _\rm {NT})$ , through log (M), and through log (R). The concentric ellipses show the 1σ, 2σ, and 3σ regions of the distribution. The green circles are randomly drawn points from P( x ). The red squares are the observational data from Motte et al. (Reference Motte, Andre and Neri1998) and André et al. (Reference André, Belloche, Motte and Peretto2007). See Lomax et al. (Reference Lomax, Whitworth, Hubber, Stamatellos and Walch2014) for the original version of this figure.

3.1.2. Shape

Molecular cloud cores often have elongated, irregular shapes. We include this in the simulations by assuming that each intrinsic core shape can be drawn from a family of triaxial ellipsoids. Each ellipsoid had axes:

(7) $$\begin{eqnarray} A&=&1,\nonumber\\ B&=&\exp (\tau \mathcal {G}_\rm {B}),\\ C&=&\exp (\tau \mathcal {G}_\rm {C}), \nonumber \end{eqnarray}$$

where $\mathcal {G}_\rm {B}$ and $\mathcal {G}_\rm {C}$ are random numbers drawn from a Gaussian distribution with zero mean and unit standard deviation. The scale-parameter τ ≈ 0.6 is a fit to the distribution of projected aspect ratios in Ophiuchus (Lomax, Whitworth, & Cartwright Reference Lomax, Whitworth and Cartwright2013). The axes are normalised to a given R, giving the dimensions of the core,

(8) $$\begin{eqnarray} A_\rm {CORE}&=&\frac{R}{(BC)^{{1}/{3}}},\nonumber\\ B_\rm {CORE}&=&BA_\rm {CORE},\\ C_\rm {CORE}&=&CA_\rm {CORE}. \nonumber \end{eqnarray}$$

3.1.3. Density profile

Density profiles of cores are often well fitted by those of critical Bonnor–Ebert spheres (e.g. Bonnor Reference Bonnor1956; Alves, Lada, & Lada Reference Alves, Lada and Lada2001; Harvey et al. Reference Harvey, Wilner, Lada, Myers, Alves and Chen2001; Lada et al. Reference Lada, Muench, Rathborne, Alves and Lombardi2008). We use such a profile for the ellipsoidal cores. Here, ρ = ρCe−ψ(ξ), where ρC is the central density, ψ is the Isothermal Function, ξ is the dimensionless radius, ξB = 6.451 is the boundary of the sphere. The density at any given point (x, y, z), where (0, 0, 0) is the centre of the core, is given by

(9) $$\begin{eqnarray} \xi &=&\xi _{{\rm FWHM}}\left(\frac{x^2}{A_\rm {CORE}^2}+\frac{y^2}{B_\rm {CORE}^2}+\frac{z^2}{C_\rm {CORE}^2}\right)^{1/2}, \end{eqnarray}$$
(10) $$\begin{eqnarray} \rho (\xi )&=&\frac{M_\rm {CORE}\xi _{{\rm B}}{\rm e}^{-\psi (\xi )}}{4\pi A_\rm {CORE}B_\rm {CORE}C_\rm {CORE}\psi ^{\prime }(\xi _{{\rm B}})},\hspace{14.22636pt}\xi <\xi _{{\rm B}}, \end{eqnarray}$$

where ψ′ is the first derivative of ψ and ξFWHM = 2.424 is the full width at half maximum of the column density through a critical Bonnor–Ebert sphere.

3.1.4. Velocity field

Each core is given a turbulent velocity field with power spectrum Pk −4 in three dimensions. We include bulk rotation and radial excursion by modifying the amplitudes a and phases φ of the k = 1 modes:

(11) $$\begin{eqnarray} \left[\begin{array}{c}\bf{ a}(1,0,0) \\ \bf{ a}(0,1,0) \\ \bf{ a}(0,0,1) \end{array}\right]=\left[\begin{array}{c@{\hspace*{8pt}}c@{\hspace*{8pt}}c}r_x & \omega _z & -\omega _y \\ -\omega _z & r_y & \omega _x \\ \omega _y & -\omega _x & r_z \end{array}\right], \end{eqnarray}$$
(12) $$\begin{eqnarray} \left[ \begin{array}{c} \bf{ \varphi }(1,0,0) \\ \bf{ \varphi }(0,1,0) \\ \bf{ \varphi }(0,0,1) \end{array}\right]= \left[\begin{array}{c@{\hspace*{8pt}}c@{\hspace*{8pt}}c}\pi /2 & \pi /2 & \pi /2 \\ \pi /2 & \pi /2 & \pi /2 \\ \pi /2 & \pi /2 & \pi /2 \end{array}\right]. \end{eqnarray}$$

The amplitude components rx, ry, rz , ω x , ω y , and ω z are drawn independently from a Gaussian distribution with zero mean and unit variance. The r terms define the amount of excursion along a given axis. The ω terms define the amount rotation about a given axis. The fields are generated on a 1283 grid and interpolated onto the SPH particles. The velocity dispersion of the particles is normalised to a given value of $\sigma _\rm {NT}$ .

3.2. Results

One hundred synthetic cores have been evolved for 2 × 105 yr. This is roughly an order of magnitude greater than the average core free-fall time and of the same order as the estimated core–core collisional timescale (André et al. Reference André, Belloche, Motte and Peretto2007). Of the one hundred cores, sixty are prestellar. Each simulation has been performed three times: once with no radiative feedback from accretion (NRF); again with episodic feedback from accretion (ERF); and finally with continuous feedback from accretion (CRF). With CRF, the protostellar luminosity was calculated using Equation (2), with $\dot{M}_\star$ set to $\dot{M}_\rm {SINK}$ .

3.2.1. Stellar masses

We find two very noticeable trends when progressing along the series NRF → ERF → CRF. First, the median number stars formed per core decreases: $\mathcal {N}_\rm {S/C}^\rm {NRF}=6.0^{+4.0}_{-2.0}$ , $\mathcal {N}_\rm {S/C}^\rm {ERF}=3.5^{+3.5}_{-2.5}$ , $\mathcal {N}_\rm {S/C}^\rm {CRF}=1.0^{+0.0}_{-0.0}$ . Second, the median protostellar mass (see Figure 2) shifts upwards: $\log _{10} M_\star ^\rm {NRF}=-1.1^{+0.4}_{-0.6}$ , $\log _{10} M_\star ^\rm {ERF}=-0.8^{+0.2}_{-0.4}$ , $\log _{10} M_\star ^\rm {CRF}=-0.5^{+0.2}_{-0.2}$ . Footnote 1

Figure 2. The black histograms show stellar mass functions for (a) NRF, (b) ERF, and (c) CRF. The blue-dotted straight lines, and the red-dashed lognormal curve, show, respectively, the Chabrier (Reference Chabrier, Corbelli, Palla and Zinnecker2005) and Kroupa (Reference Kroupa2001) fits to the observed IMF. The vertical-dashed line shows the hydrogen burning limit at M = 0.08 M. See Lomax et al. (Reference Lomax, Whitworth, Hubber, Stamatellos and Walch2015b) for the original version of this figure.

These trends occur because disc fragmentation is strongly affected by protostellar feedback. In the case with NRF, discs are relatively cold and easily fragment. Recall from Equation (1) that fragmentation can occur if the disc has a sufficiently large column density or the sound speed is sufficiently low. With ERF, fragmentation is occasionally interrupted by the episodic outbursts. With CRF, the discs are constantly heated, and fragmentation becomes difficult. As a consequence, the central protostar or binary usually accretes the entire mass of its disc. Protostars with NRF, and to a slightly lesser extent ERF, are often attended by multiple low-mass companions which partially starve the primary protostar of the remaining gas. Of the three sets of simulations, those with ERF best reproduce the Chabrier (Reference Chabrier, Corbelli, Palla and Zinnecker2005) IMF.

In star-forming regions, the observed ratio of stars to brown dwarfs is

(13) $$\begin{equation} \mathcal {A}=\frac{N(0.08\,{\rm M}_{\odot }<M\le 1.0\,{\rm M}_{\odot })}{N(0.03\,{\rm M}_{\odot }<M\le 0.08\,{\rm M}_{\odot })}=4.3\pm 1.6, \end{equation}$$

(See Andersen et al. Reference Andersen, Meyer, Greissl and Aversa2008). This figure is best reproduced with ERF: $\mathcal {A}_\rm {ERF}=3.9\pm 0.6$ .Footnote 2 Simulations with NRF and CRF yield $\mathcal {A}_\rm {NRF}=2.2\pm 0.3$ and $\mathcal {A}_\rm {CRF}=17\pm 8$ . While the figure for NRF is not completely incompatible with observations, the value with CRF is far too high. This is because brown dwarfs are unable to form via disc fragmentation (e.g. Stamatellos, Hubber, & Whitworth Reference Stamatellos, Hubber and Whitworth2007a; Stamatellos & Whitworth Reference Stamatellos and Whitworth2009).

3.2.2. Stellar multiplicities

The core simulations also produce a wide variety of multiple systems. These are either simple binary systems or high order (N ⩾ 3) hierarchical multiples. A high-order system can be thought of as a binary where one or both components is another binary system. Multiple systems of protostars are extracted from the end state simulations if they have been tidally stable for at least one orbital period.

There are many ways to statistically describe the multiplicity of a population of stellar systems. Here, following Reipurth & Zinnecker (Reference Reipurth and Zinnecker1993), we use the multiplicity frequency and the pairing factor. The multiplicity frequency is fraction of systems which is multiple

(14) $$\begin{equation} {\rm mf}=\frac{B+T+Q\dots }{S+B+T+Q\dots }, \end{equation}$$

where S is the number of single systems, B is the number of binaries, T is the number of triples, etc. The pairing factor is the average number of orbits per system

(15) $$\begin{equation} {\rm pf}=\frac{B+2T+3Q\dots }{S+B+T+Q\dots }. \end{equation}$$

These two quantities are particularly useful in conjunction, as their ratio gives the average number of objects per multiple system

(16) $$\begin{equation} \mathcal {O}_\rm {SYS}=1+\frac{{\rm pf}}{{\rm mf}}. \end{equation}$$

By definition, $\mathcal {O}_\rm {SYS}\ge 2$ .

Figure 3 shows mf, pf, and $\mathcal {O}_\rm {SYS}$ for systems with very low mass or brown dwarf primaries, M-dwarf primaries, and solar-type primariesFootnote 3 . For comparison, we have also included the same figure for stars in the field (See Duchêne & Kraus Reference Duchêne and Kraus2013, and references therein). Observations of pre-main-sequence stars in Ophiuchus and Taurus suggest that multiplicity is high when stars are young (e.g. Leinert, Zinnecker, Weitzel, Christou, Ridgway, Jameson, Haas Reference Leinert, Zinnecker, Weitzel, Christou, Ridgway, Jameson, Haas and Lenzen1993; Ratzka, Köhler, & Leinert Reference Ratzka, Köhler and Leinert2005; Kraus et al. Reference Kraus, Ireland, Martinache and Hillenbrand2011). As stars age, particularly if they are in a clustered environment, dynamical interactions erode these systems and the multiplicity frequency drops to that observed in the field (e.g. Kroupa Reference Kroupa1995; Parker et al. Reference Parker, Goodwin, Kroupa and Kouwenhoven2009; Parker & Goodwin Reference Parker and Goodwin2011, Reference Parker and Goodwin2012). The multiplicity of protostars (i.e. the objects formed in these simulations) should therefore be higher than the multiplicity of field stars. This requirement is met for the simulations with NRF and ERF. However, with CRF, the multiplicity frequency is too low for solar-type primaries.

Figure 3. Multiplicity frequency (a), pairing factor (b), and mean system order (c) for systems with very low mass, M-dwarf and solar-type primaries. The red boxes give the values for the NRF simulations, blue for the ERF, and green for the CRF. The black points give the values observed in field main-sequence stars. In all cases, the width of a box shows the extent of the mass bin, and the height shows the uncertainty.

Very high order systems, e.g. sextuples, are also form in these simulations. Figure 4 shows a sextuple system formed from a solar mass core with ERF. Initially the core fragments into two objects, which form a binary system with a circumbinary disc. This disc fragments into a further five objects. It finally settles into a sextuple system (a binary system orbiting a quadruple) with separations ranging from ~ 600 AU for the outer orbit to ~ 0.1 AU for the innermost orbits. Sextuple systems similar to this are found in star-forming regions (Kraus et al. Reference Kraus, Ireland, Martinache and Hillenbrand2011) and the field Eggleton & Tokovinin (Reference Eggleton and Tokovinin2008); Tokovinin (Reference Tokovinin2008).

Figure 4. A sequence of column density maps of a core during disc fragmentation. The initial core has M = 1.3 M, R = 3000 AU, $\sigma _\rm {NT}=0.3\,{\rm km\,s^{-1}}$ and is evolved with ERF. The colour-bar gives shows the column density in units of g cm−2. The black dots show the positions of sink particles, i.e. protostars. Figure 4(a) shows gravitational instabilities developing in a circumbinary disc. Figure 4(b) shows seven protostars in an unstable configuration. In Figure 4(c), the protostars are configured in a quadruple system (right) and a binary (left); a single protostar is being ejected (centre). Figure 4(d) shows a stable sextuple system which lasts until the end of the simulation. (a) t = 2.2 × 104 yr (b) t = 2.7 × 104 yr (c) t = 3.2 × 104 yr (d) t = 3.7 × 104 yr

All three sets of simulations produce systems with semi-major axes ranging between roughly 0.1 AU (the resolution limit of the simulation) and 1000 AU. This upper limit corresponds to that observed in star-forming regions (e.g. King et al. Reference King, Parker, Patience and Goodwin2012a, Reference King, Goodwin, Parker and Patience2012b). Systems with wider orbits ( $\gtrsim 10^3\,{\rm \text{AU}}$ ) are probably assembled later through dynamical interactions in clustered environments (e.g. Kouwenhoven et al. Reference Kouwenhoven, Goodwin, Parker, Davies, Malmberg and Kroupa2010).

3.3. Summary

We find that the mass distribution and multiplicity statistics of young stars can be recovered from simulations if radiative feedback from protostellar accretion is episodic. The periods of low luminosity provide a window of opportunity, during which protostellar discs can fragment. When feedback is continuous, the disc is too warm to permit disc fragmentation. In this instance, the simulations result in (i) a top-heavy IMF, and (ii) too few stars per core to satisfactorily reproduce observed multiplicity statistics. When there is no radiative feedback, good multiplicity statistics are recovered, but the resultant IMF has too many brown dwarfs.

4 TURBULENCE: DISCS AND FILAMENTS

To examine how the structure of a core’s velocity field affects the star formation process, we take a single core set-up and vary the structure of the velocity field. The core is spherical, with M = 3 M, R = 3000 AU, and $\sigma _\rm {NT}=0.44\,{\rm km\,s^{-1}}$ . These dimensions are similar to those of SM1 in Oph-A region (Motte et al. Reference Motte, Andre and Neri1998; André et al. Reference André, Belloche, Motte and Peretto2007). We alter the partition of kinetic energy in solenoidal modes (i.e. shear and rotation) and compressive modes (i.e. compression and rarefaction).

4.1. Initial conditions

The amplitude of a turbulent mode $\bf {a}(\bf{ k})$ can be split into its longitudinal (compressive) component $\bf {a_\rm {L}}(\bf{ k})$ and transverse (solenoidal) component $\bf {a_\rm {T}}(\bf{ k})$ using Helmholtz decomposition:

(17) $$\begin{eqnarray} \bf{ a_\rm {L}}(\bf{ k})&=&\hat{\bf{k}}(\bf{ a}(\bf{ k})\cdot \hat{\bf{k}}),\nonumber\\[-6pt] &&\\[-6pt] \bf{ a_\rm {T}}(\bf{ k})&=&\bf{ a}(\bf{ k})-\hat{\bf{k}}(\bf{ a}(\bf{ k})\cdot \hat{\bf{k}}).\nonumber \end{eqnarray}$$

When the magnitude and direction of the amplitude is random, there is on average twice as much energy in transverse modes as there is in solenoidal modes.

We define the parameter $\delta _\rm {SOL}$ as the average fraction of solenoidal kinetic energy in a velocity field. We modify the field to have given $\delta _\rm {SOL}$ by performing the transformation

(18) $$\begin{equation} \bf{ a}(\bf{ k})\rightarrow \sqrt{3\,(1-\delta _\rm {SOL})}\,\bf{ a_\rm {L}}(\bf{ k})+\sqrt{\frac{3}{2}\,\delta _\rm {SOL}}\,\bf{ a_\rm {T}}(\bf{ k}). \end{equation}$$

We generate ten initial cores, each with a unique random velocity seed. For each core, we apply the transformation in Equation (18) with values $\delta _\rm {SOL}=0,\,\frac{1}{9},\,\frac{1}{3},\,\frac{2}{3},\,1$ , yielding a total fifty core set-ups.

4.2. Results

Figure 5 shows a montage of simulation snapshots where the random seed is fixed and $\delta _\rm {SOL}$ is varied from 0 to 1. When the field is purely compressive ( $\delta _\rm {SOL}=0$ ), protostars form within a network of filaments. Due to the low angular momentum of the system, the resultant protostars are only attended by small discs. In contrast, when the field is purely solenoidal ( $\delta _\rm {SOL}=1)$ , a single protostar forms and is attended by a coherent disc structure. The disc proceeds to fragment into multiple objects. A smooth transition between filament fragmentation and disc fragmentation is seen in the snapshots with intermediate values of $\delta _\rm {SOL}$ .

Figure 5. Column density maps of the central 820 AU × 820 AU of the (x, y)-plane, from the simulations with fixed random seed and different values of $\delta _\rm {SOL}$ , at times $t=1.00,1.25,1.50,\text{ and }1.75\times 10^4\,{\rm yr}$ . The colour scale gives the logarithmic column density in units of g cm−2. Sink particles are represented by black dots. See Lomax et al. (Reference Lomax, Whitworth and Hubber2015a) for the original version of this figure. (a) $\delta _\rm {SOL}=0$ . (b) $\delta _\rm {SOL}=1/9$ . (c) $\delta _\rm {SOL}=1/3$ . (d) $\delta _\rm {SOL}=2/3$ . (e) $\delta _\rm {SOL}=1$ .

Figure 6 shows the fraction of protostars formed by filament fragmentation (i.e. in relative isolation) and disc fragmentation (i.e. in discs around more massive protostars) as a function of $\delta _\rm {SOL}$ . These values are averaged over all random seeds. Here, we see that the occurrence of disc fragmentation increases monotonically with $\delta _\rm {SOL}$ . Filament fragmentation, therefore, decreases monotonically with $\delta _\rm {SOL}$ .

Figure 6. The fraction of stars formed by filament fragmentation (red crosses) and disc fragmentation (green boxes) for different values of $\delta _\rm {SOL}$ , averaged over all random seeds. The error bars show the Poison counting uncertainties. See Lomax et al. (Reference Lomax, Whitworth and Hubber2015a) for the original version of this figure.

On average, the number of protostars spawned per core ranges from 5.4 ± 7 when $\delta _\rm {SOL}=0$ to 8.1 ± 9 when $\delta _\rm {SOL}=1$ . The median and interquartile range of mass is shown in Figure 7. Here, we see that increasing $\delta _\rm {SOL}$ pushes the median sink mass down from roughly 0.6 M when $\delta _\rm {SOL}=0$ to 0.3 M when $\delta _\rm {SOL}=1$ . We also find that purely compressive fields form very few brown dwarfs. We note that—even for this limited set of initial conditions—the interquartile range of protostellar masses when $\delta _\rm {SOL}\gtrsim 2/3$ is very similar to that of the Chabrier (Reference Chabrier, Corbelli, Palla and Zinnecker2005) IMF.

Figure 7. The black points show the median stellar mass, and the vertical black bars show interquartile range of mass, for different values of $\delta _\rm {SOL}$ , averaged over all random seeds. The solid and dashed horizontal red lines show the median and interquartile range for the Chabrier (Reference Chabrier, Corbelli, Palla and Zinnecker2005) IMF. See Lomax et al. (Reference Lomax, Whitworth and Hubber2015a) for the original version of this figure.

4.3. Summary

We show that the collapse and fragmentation of prestellar cores is strongly influenced by the structure of the velocity field. Disc formation and fragmentation dominates the star formation process when $\delta _\rm {SOL}\gtrsim 1/3$ . At values below this, stars form mostly through the fragmentation of filamentary structures. The value of $\delta _\rm {SOL}$ also affects distribution of protostellar masses. The distribution most resembles the observed IMF when $\delta _\rm {SOL}\gtrsim 2/3$ . Reducing $\delta _\rm {SOL}$ reduces the level of disc fragmentation, resulting in smaller number of objects with greater average mass. In extreme cases (i.e. $\delta _\rm {SOL}<1/9$ ), the formation of low-mass stars and brown dwarfs is heavily suppressed. This suggests that disc fragmentation may be a requirement for forming these objects.

5 SUMMARY AND CONCLUSIONS

Previous numerical work (e.g. Stamatellos et al. Reference Stamatellos, Hubber and Whitworth2007a; Stamatellos & Whitworth Reference Stamatellos and Whitworth2009) shows that disc fragmentation is important mechanism for reproducing the properties of low-mass stars. We demonstrate that disc fragmentation also plays an important role in the conversion of prestellar cores into stars. Importantly, the observed masses and multiplicities of stars can be recovered from simulations if the following criteria are satisfied:

  • Radiative feedback from accretion onto protostars is episodic. Simulations with episodic radiative feedback produce both an IMF and multiplicity statistics in good agreement with those observed. Simulations with continuous radiative feedback fail to produce both the observed number of brown dwarfs and the multiplicity statistics associated with young objects. Furthermore, continuous radiative feedback produces protostellar luminosities greater than those observed in young stars. Simulations with no radiative feedback can produce good multiplicity statistics and a reasonable—albeit bottom-heavy—IMF, but are unrealistic.

  • A significant proportion of the core’s internal kinetic energy is in solenoidal turbulent modes. Cores with more than a third of their kinetic energy in solenoidal modes are able to easily produce stars via disc fragmentation. Decreasing this fraction results in more stars forming in fragmenting filaments. These objects tend to be of greater mass, resulting in a top-heavy mass distribution relative to the observed IMF. Furthermore, filament fragmentation struggles to produce brown dwarfs and very low mass stars.

ACKNOWLEDGEMENTS

OL and APW gratefully acknowledge the support of a consolidated grant (ST/K00926/1) from the UK STFC. This work was performed using the computational facilities of the Advanced Research Computing @ Cardiff (ARCCA) Division, Cardiff University. All false-colour images have been rendered with splash (Price Reference Price2007).

Footnotes

1 The uncertainties give the interquartile range of the distribution.

2 Here, the uncertainty is calculated from the Poison counting error.

3 Here, we define very low mass stars and brown dwarfs as stars with 0.06 MM < 0.1 M, M-dwarfs as stars with 0.1 MM < 0.5 M and solar types as stars with 0.7 MM < 1.3 M.

References

REFERENCES

Alves, J. F., Lada, C. J., & Lada, E. A. 2001, Nature, 409, 159 2001Natur.409..159A Google Scholar
Alves, J., Lombardi, M., & Lada, C. J. 2007, A&A, 462, L17 10.1051/0004-6361:20066389 2007A&A...462L..17A Google Scholar
Andersen, M., Meyer, M. R., Greissl, J., & Aversa, A. 2008, ApJ, 683, L183 10.1086/591473 2008ApJ...683L.183A CrossRefGoogle Scholar
André, P., Belloche, A., Motte, F., & Peretto, N. 2007, A&A, 472, 519 10.1051/0004-6361:20077422 2007A&A...472..519A Google Scholar
Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 10.1086/172425 1993ApJ...406..122A Google Scholar
Andre, P., Ward-Thompson, D., & Barsony, M. 2000, Protostars and Planets IV (Tucson: Univ. Arizona Press), 59 2000prpl.conf...59A Google Scholar
Bate, M. R. 1998, ApJ, 508, L95 10.1086/311719 1998ApJ...508L..95B Google Scholar
Bate, M. R. 2000, MNRAS, 314, 33 10.1046/j.1365-8711.2000.03333.x 2000MNRAS.314...33B CrossRefGoogle Scholar
Bonnor, W. B. 1956, MNRAS, 116, 351 1956MNRAS.116..351B Google Scholar
Chabrier, G. 2003, ApJ, 586, L133 10.1086/374879 2003ApJ...586L.133C CrossRefGoogle Scholar
Chabrier, G. 2005, in Astrophysics and Space Science Library, Vol. 327, The Initial Mass Function 50 Years Later, eds. Corbelli, E., Palla, F., Zinnecker, H. (Dordrecht: Springrer), p. 41 (arXiv:astro-ph/0409465)astro-ph/0409465 CrossRefGoogle Scholar
Delgado-Donate, E. J., Clarke, C. J., & Bate, M. R. 2004a, MNRAS, 347, 759 10.1111/j.1365-2966.2004.07259.x 2004MNRAS.347..759D Google Scholar
Delgado-Donate, E. J., Clarke, C. J., Bate, M. R., & Hodgkin, S. T. 2004b, MNRAS, 351, 617 10.1111/j.1365-2966.2004.07803.x 2004MNRAS.351..617D CrossRefGoogle Scholar
Duchêne, G., & Kraus, A. 2013, ARA&A, 51, 269 10.1146/annurev-astro-081710-102602 2013ARA&A..51..269D Google Scholar
Eggleton, P. P., & Tokovinin, A. A. 2008, MNRAS, 389, 869 10.1111/j.1365-2966.2008.13596.x 2008MNRAS.389..869E Google Scholar
Enoch, M. L., Evans, II, N. J., Sargent, A. I., Glenn, J., Rosolowsky, E., & Myers, P. 2008, ApJ, 684, 1240 10.1086/589963 2008ApJ...684.1240E Google Scholar
Enoch, M. L., et al. 2006, ApJ, 638, 293 10.1086/498678 2006ApJ...638..293E CrossRefGoogle Scholar
Forgan, D., & Rice, K. 2013, MNRAS, 430, 2082 10.1093/mnras/stt032 2013MNRAS.430.2082F CrossRefGoogle Scholar
Gammie, C. F. 2001, ApJ, 553, 174 10.1086/320631 2001ApJ...553..174G Google Scholar
Goodwin, S. P., & Whitworth, A. P. 2004, A&A, 413, 929 10.1051/0004-6361:20031529 2004A&A...413..929G Google Scholar
Goodwin, S. P., Whitworth, A. P., & Ward-Thompson, D. 2004, A&A, 423, 169 10.1051/0004-6361:20040285 2004A&A...423..169G Google Scholar
Goodwin, S. P., Whitworth, A. P., & Ward-Thompson, D. 2006, A&A, 452, 487 10.1051/0004-6361:20054026 2006A&A...452..487G Google Scholar
Green, J. D., et al. 2011, ApJ, 731, L25 10.1088/2041-8205/731/2/L25 2011ApJ...731L..25G Google Scholar
Greene, T. P., Aspin, C., & Reipurth, B. 2008, ApJ, 135, 1421 10.1088/0004-6256/135/4/1421 2008AJ....135.1421G CrossRefGoogle Scholar
Harvey, D. W. A., Wilner, D. J., Lada, C. J., Myers, P. C., Alves, J. F., & Chen, H. 2001, ApJ, 563, 903 10.1086/324076 2001ApJ...563..903H Google Scholar
Hennebelle, P., & Chabrier, G. 2008, ApJ, 684, 395 10.1086/589916 2008ApJ...684..395H Google Scholar
Hennebelle, P., & Chabrier, G. 2009, ApJ, 702, 1428 10.1088/0004-637X/702/2/1428 2009ApJ...702.1428H Google Scholar
Herbig, G. H. 1977, ApJ, 217, 693 10.1086/155615 1977ApJ...217..693H Google Scholar
Holman, K., Walch, S. K., Goodwin, S. P., & Whitworth, A. P. 2013, MNRAS, 432, 3534 10.1093/mnras/stt705 2013MNRAS.432.3534H Google Scholar
Horton, A. J., Bate, M. R., & Bonnell, I. A. 2001, MNRAS, 321, 585 10.1046/j.1365-8711.2001.04071.x 2001MNRAS.321..585H CrossRefGoogle Scholar
Hubber, D. A., Batty, C. P., McLeod, A., & Whitworth, A. P. 2011, A&A, 529, A27 10.1051/0004-6361/201014949 2011A&A...529A..27H Google Scholar
Hubber, D. A., Walch, S., & Whitworth, A. P. 2013, MNRAS, 430, 3261 10.1093/mnras/stt128 2013MNRAS.430.3261H CrossRefGoogle Scholar
Janson, M., et al. 2012, ApJ, 754, 44 10.1088/0004-637X/754/1/44 2012ApJ...754...44J CrossRefGoogle Scholar
Johnstone, D., & Bally, J. 2006, ApJ, 653, 383 10.1086/508852 2006ApJ...653..383J CrossRefGoogle Scholar
Johnstone, D., Fich, M., Mitchell, G. F., & Moriarty-Schieven, G. 2001, ApJ, 559, 307 10.1086/322323 2001ApJ...559..307J Google Scholar
Johnstone, D., Wilson, C. D., Moriarty-Schieven, G., Joncas, G., Smith, G., Gregersen, E., & Fich, M. 2000, ApJ, 545, 327 10.1086/317790 2000ApJ...545..327J CrossRefGoogle Scholar
Kenyon, S. J., & Hartmann, L. 1995, ApJ, 101, 117 10.1086/192235 1995ApJS..101..117K Google Scholar
Kenyon, S. J., Hartmann, L. W., Strom, K. M., & Strom, S. E. 1990, AJ, 99, 869 10.1086/115380 1990AJ.....99..869K CrossRefGoogle Scholar
King, R. R., Goodwin, S. P., Parker, R. J., & Patience, J. 2012b, MNRAS, 427, 2636 10.1111/j.1365-2966.2012.22108.x 2012MNRAS.427.2636K Google Scholar
King, R. R., Parker, R. J., Patience, J., & Goodwin, S. P. 2012a, MNRAS, 421, 2025 10.1111/j.1365-2966.2012.20437.x 2012MNRAS.421.2025K CrossRefGoogle Scholar
Könyves, V., et al. 2010, A&A, 518, L106 10.1051/0004-6361/201014689 2010A&A...518L.106K Google Scholar
Kouwenhoven, M. B. N., Goodwin, S. P., Parker, R. J., Davies, M. B., Malmberg, D., & Kroupa, P. 2010, MNRAS, 404, 1835 10.1111/j.1365-2966.2010.16399.x 2010MNRAS.404.1835K Google Scholar
Kraus, A. L., Ireland, M. J., Martinache, F., & Hillenbrand, L. A. 2011, ApJ, 731, 8 10.1088/0004-637X/731/1/8 2011ApJ...731....8K Google Scholar
Kroupa, P. 1995, MNRAS, 277, 1491 1995MNRAS.277.1491K Google Scholar
Kroupa, P. 2001, MNRAS, 322, 231 10.1046/j.1365-8711.2001.04022.x 2001MNRAS.322..231K Google Scholar
Lada, C. J., Muench, A. A., Rathborne, J., Alves, J. F., & Lombardi, M. 2008, ApJ, 672, 410 10.1086/523837 2008ApJ...672..410L CrossRefGoogle Scholar
Leinert, C., Zinnecker, H., Weitzel, N., Christou, J., Ridgway, S. T., Jameson, R., Haas, M., & Lenzen, R. 1993, A&A, 278, 129 1993A&A...278..129L Google Scholar
Lodato, G., & Rice, W. K. M. 2005, MNRAS, 358, 1489 10.1111/j.1365-2966.2005.08875.x 2005MNRAS.358.1489L Google Scholar
Lomax, O., Whitworth, A. P., & Cartwright, A. 2013, MNRAS, 436, 2680 10.1093/mnras/stt1764 2013MNRAS.tmp.2484L Google Scholar
Lomax, O., Whitworth, A. P., & Hubber, D. A. 2015a, MNRAS, 449, 662 10.1093/mnras/stv310 2015MNRAS.449..662L CrossRefGoogle Scholar
Lomax, O., Whitworth, A. P., Hubber, D. A., Stamatellos, D., & Walch, S. 2014, MNRAS, 439, 3039 10.1093/mnras/stu177 2014MNRAS.439.3039L CrossRefGoogle Scholar
Lomax, O., Whitworth, A. P., Hubber, D. A., Stamatellos, D., & Walch, S. 2015b, MNRAS, 447, 1550 10.1093/mnras/stu2530 2015MNRAS.447.1550L Google Scholar
Matsumoto, T., & Hanawa, T. 2003, ApJ, 595, 913 10.1086/377367 2003ApJ...595..913M Google Scholar
Morris, J. P., & Monaghan, J. J. 1997, JCoPh, 136, 41 10.1006/jcph.1997.5690 1997JCoPh.136...41M Google Scholar
Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 1998A&A...336..150M Google Scholar
Motte, F., André, P., Ward-Thompson, D., & Bontemps, S. 2001, A&A, 372, L41 10.1051/0004-6361:20010543 2001A&A...372L..41M Google Scholar
Nutter, D., & Ward-Thompson, D. 2007, MNRAS, 374, 1413 10.1111/j.1365-2966.2006.11246.x 2007MNRAS.374.1413N Google Scholar
Offner, S. S. R., Klein, R. I., McKee, C. F., & Krumholz, M. R. 2009, ApJ, 703, 131 10.1088/0004-637X/703/1/131 2009ApJ...703..131O Google Scholar
Padoan, P., & Nordlund, Å. 2002, ApJ, 576, 870 10.1086/341790 2002ApJ...576..870P Google Scholar
Palla, F., & Stahler, S. W. 1993, ApJ, 418, 414 10.1086/173402 1993ApJ...418..414P Google Scholar
Parker, R. J., & Goodwin, S. P. 2011, MNRAS, 411, 891 10.1111/j.1365-2966.2010.17722.x 2011MNRAS.411..891P Google Scholar
Parker, R. J., & Goodwin, S. P. 2012, MNRAS, 424, 272 10.1111/j.1365-2966.2012.21190.x 2012MNRAS.424..272P Google Scholar
Parker, R. J., Goodwin, S. P., Kroupa, P., & Kouwenhoven, M. B. N. 2009, MNRAS, 397, 1577 10.1111/j.1365-2966.2009.15032.x 2009MNRAS.397.1577P Google Scholar
Pattle, K., et al. 2015, MNRAS, 450, 1094 10.1093/mnras/stv376 2015MNRAS.450.1094P Google Scholar
Peneva, S. P., Semkov, E. H., Munari, U., & Birkle, K. 2010, A&A, 515, A24 10.1051/0004-6361/201014092 2010A&A...515A..24P Google Scholar
Price, D. J. 2007, PASA, 24, 159 10.1071/AS07022 2007PASA...24..159P Google Scholar
Principe, D. A., Kastner, J. H., Grosso, N., Hamaguchi, K., Richmond, M., Teets, W. K., & Weintraub, D. A. 2013, preprint (arXiv:1311.5232)2013arXiv1311.5232P 1311.5232 Google Scholar
Raghavan, D., et al. 2010, ApJ, 190, 1 10.1088/0067-0049/190/1/1 2010ApJS..190....1R Google Scholar
Rathborne, J. M., Lada, C. J., Muench, A. A., Alves, J. F., Kainulainen, J., & Lombardi, M. 2009, ApJ, 699, 742 10.1088/0004-637X/699/1/742 2009ApJ...699..742R Google Scholar
Ratzka, T., Köhler, R., & Leinert, C. 2005, A&A, 437, 611 10.1051/0004-6361:20042107 2005A&A...437..611R Google Scholar
Reipurth, B., & Zinnecker, H. 1993, A&A, 278, 81 1993A&A...278...81R Google Scholar
Scholz, A., Froebrich, D., & Wood, K. 2013, MNRAS, 430, 2910 10.1093/mnras/stt091 2013MNRAS.430.2910S Google Scholar
Simpson, R. J., Nutter, D., & Ward-Thompson, D. 2008, MNRAS, 391, 205 10.1111/j.1365-2966.2008.13750.x 2008MNRAS.391..205S Google Scholar
Stamatellos, D., Hubber, D. A., & Whitworth, A. P. 2007a, MNRAS, 382, L30 10.1111/j.1745-3933.2007.00383.x 2007MNRAS.382L..30S Google Scholar
Stamatellos, D., & Whitworth, A. P. 2009, MNRAS, 392, 413 10.1111/j.1365-2966.2008.14069.x 2009MNRAS.392..413S Google Scholar
Stamatellos, D., Whitworth, A. P., Bisbas, T., & Goodwin, S. 2007b, A&A, 475, 37 10.1051/0004-6361:20077373 2007A&A...475...37S Google Scholar
Stamatellos, D., Whitworth, A. P., & Hubber, D. A. 2011, ApJ, 730, 32 10.1088/0004-637X/730/1/32 2011ApJ...730...32S Google Scholar
Stamatellos, D., Whitworth, A. P., & Hubber, D. A. 2012, MNRAS, 427, 1182 10.1111/j.1365-2966.2012.22038.x 2012MNRAS.427.1182S Google Scholar
Stanke, T., Smith, M. D., Gredel, R., & Khanzadyan, T. 2006, A&A, 447, 609 10.1051/0004-6361:20041331 2006A&A...447..609S Google Scholar
Testi, L., & Sargent, A. I. 1998, ApJ, 508, L91 10.1086/311724 1998ApJ...508L..91T Google Scholar
Tokovinin, A. 2008, MNRAS, 389, 925 10.1111/j.1365-2966.2008.13613.x 2008MNRAS.389..925T Google Scholar
Toomre, A. 1964, ApJ, 139, 1217 10.1086/147861 1964ApJ...139.1217T Google Scholar
Tsukamoto, Y., Takahashi, S. Z., Machida, M. N., & Inutsuka, S. 2015, MNRAS, 446, 1175 10.1093/mnras/stu2160 2015MNRAS.446.1175T Google Scholar
Walch, S., Burkert, A., Whitworth, A., Naab, T., & Gritschneder, M. 2009, MNRAS, 400, 13 10.1111/j.1365-2966.2009.15293.x 2009MNRAS.400...13W CrossRefGoogle Scholar
Walch, S., Naab, T., Whitworth, A., Burkert, A., & Gritschneder, M. 2010, MNRAS, 402, 2253 10.1111/j.1365-2966.2009.16058.x 2010MNRAS.402.2253W Google Scholar
Walch, S., Whitworth, A. P., & Girichidis, P. 2012, MNRAS, 419, 760 10.1111/j.1365-2966.2011.19741.x 2012MNRAS.419..760W Google Scholar
Zhu, Z., Hartmann, L., & Gammie, C. 2009, ApJ, 694, 1045 10.1088/0004-637X/694/2/1045 2009ApJ...694.1045Z Google Scholar
Zhu, Z., Hartmann, L., & Gammie, C. 2010a, ApJ, 713, 1143 10.1088/0004-637X/713/2/1143 2010ApJ...713.1143Z Google Scholar
Zhu, Z., Hartmann, L., Gammie, C. F., Book, L. G., Simon, J. B., & Engelhard, E. 2010b, ApJ, 713, 1134 10.1088/0004-637X/713/2/1134 2010ApJ...713.1134Z Google Scholar
Figure 0

Table 1. Arithmetic means, standard deviations, and correlation coefficients of log (M), log (R), and $\log (\sigma _\rm {NT })$ for cores in Ophiuchus.

Figure 1

Figure 1. The multivariate lognormal distribution, P(x) where $\bf{ x}=(\log (M),\log (R),\log (\sigma _\rm {NT}))$. The figure shows the projections through $\log (\sigma _\rm {NT})$, through log (M), and through log (R). The concentric ellipses show the 1σ, 2σ, and 3σ regions of the distribution. The green circles are randomly drawn points from P(x). The red squares are the observational data from Motte et al. (1998) and André et al. (2007). See Lomax et al. (2014) for the original version of this figure.

Figure 2

Figure 2. The black histograms show stellar mass functions for (a) NRF, (b) ERF, and (c) CRF. The blue-dotted straight lines, and the red-dashed lognormal curve, show, respectively, the Chabrier (2005) and Kroupa (2001) fits to the observed IMF. The vertical-dashed line shows the hydrogen burning limit at M = 0.08 M. See Lomax et al. (2015b) for the original version of this figure.

Figure 3

Figure 3. Multiplicity frequency (a), pairing factor (b), and mean system order (c) for systems with very low mass, M-dwarf and solar-type primaries. The red boxes give the values for the NRF simulations, blue for the ERF, and green for the CRF. The black points give the values observed in field main-sequence stars. In all cases, the width of a box shows the extent of the mass bin, and the height shows the uncertainty.

Figure 4

Figure 4. A sequence of column density maps of a core during disc fragmentation. The initial core has M = 1.3 M, R = 3000 AU, $\sigma _\rm {NT}=0.3\,{\rm km\,s^{-1}}$ and is evolved with ERF. The colour-bar gives shows the column density in units of g cm−2. The black dots show the positions of sink particles, i.e. protostars. Figure 4(a) shows gravitational instabilities developing in a circumbinary disc. Figure 4(b) shows seven protostars in an unstable configuration. In Figure 4(c), the protostars are configured in a quadruple system (right) and a binary (left); a single protostar is being ejected (centre). Figure 4(d) shows a stable sextuple system which lasts until the end of the simulation. (a) t = 2.2 × 104 yr (b) t = 2.7 × 104 yr (c) t = 3.2 × 104 yr (d) t = 3.7 × 104 yr

Figure 5

Figure 5. Column density maps of the central 820 AU × 820 AU of the (x, y)-plane, from the simulations with fixed random seed and different values of $\delta _\rm {SOL}$, at times $t=1.00,1.25,1.50,\text{ and }1.75\times 10^4\,{\rm yr}$. The colour scale gives the logarithmic column density in units of g cm−2. Sink particles are represented by black dots. See Lomax et al. (2015a) for the original version of this figure. (a) $\delta _\rm {SOL}=0$. (b) $\delta _\rm {SOL}=1/9$. (c) $\delta _\rm {SOL}=1/3$. (d) $\delta _\rm {SOL}=2/3$. (e) $\delta _\rm {SOL}=1$.

Figure 6

Figure 6. The fraction of stars formed by filament fragmentation (red crosses) and disc fragmentation (green boxes) for different values of $\delta _\rm {SOL}$, averaged over all random seeds. The error bars show the Poison counting uncertainties. See Lomax et al. (2015a) for the original version of this figure.

Figure 7

Figure 7. The black points show the median stellar mass, and the vertical black bars show interquartile range of mass, for different values of $\delta _\rm {SOL}$, averaged over all random seeds. The solid and dashed horizontal red lines show the median and interquartile range for the Chabrier (2005) IMF. See Lomax et al. (2015a) for the original version of this figure.