Monte Carlo Radiative Transfer in Embedded Prestellar Cores

We implement a Monte Carlo radiative transfer method, that uses a large number of monochromatic luminosity packets to represent the radiation transported through a system. These packets are injected into the system and interact stochastically with it. We test our code against various benchmark calculations and determine the number of packets required to obtain accurate results under different circumstances. We then use this method to study cores that are directly exposed to the interstellar radiation field (non-embedded cores) and find similar results with previous studies. We also explore a large number of models of cores that are embedded in the centre of a molecular cloud. Our study indicates that the temperature profiles in embedded cores are less steep than those in non-embedded cores. Deeply embedded cores (ambient cloud with visual extinction larger than 15-25) are almost isothermal at around 7-8 K. The temperature inside cores surrounded by an ambient cloud of even moderate thickness (Av~5) is less than 12 K, which is lower than previous studies have assumed. Thus, previous mass calculations of embedded cores (for example in the rho Ophiuchi protocluster), based on mm continuum observations, may underestimate core masses by up to a factor of 2. 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. The upcoming HERSCHEL mission (ESA, 2007) will, in principle, be able to detect these features and test our models.Comment: 15 pages, 18 figures, accepted by A&

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.Comment: 6 pages. Accepted by MNRA

High-resolution simulations of clump-clump collisions using SPH with Particle Splitting

We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using an SPH code which satisfies the Jeans condition by invoking On-the-Fly Particle Splitting. Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor-Ebert, spheres. Collisions are characterised by M_0 (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius), and M (Mach Number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling, and (iii) adiabatic heating in optically thick protostellar fragments. The efficiency of star formation is found to vary between 10% and 30% in the different collisions studied and it appears to increase with decreasing M_0, and/or decreasing b, and/or increasing M. For b<0.5 collisions produce shock compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high, 1 to 5 x 10^{-5} M_sun yr^{-1}, for the first 1 to 3 x 10^4 yrs. The densities in the filaments, n >~ 5 x 10^5 cm^{-3}, are sufficient that they could be mapped in NH_3 or CS line radiation, in nearby star formation regions.Comment: Accepted for publication in MNRAS; 21 pages; 25 figures. Four figures are provided separately in reduced jpg format due to their large original ps size: click on "PostScript" to have direct access to the 4 jpg figures; full size ps files for these 4 figures can be found at http://www.aip.de/People/skitsionas/papers

An empirical model for protostellar collapse

We propose a new analytic model for the initial conditions of protostellar collapse in relatively isolated regions of star formation. The model is non-magnetic, and is based on a Plummer-like radial density profile as its initial condition. It fits: the observed density profiles of pre-stellar cores and Class 0 protostars; recent observations in pre-stellar cores of roughly constant contraction velocities over a wide range of radii; and the lifetimes and accretion rates derived for Class 0 and Class I protostars. However, the model is very simple, having in effect only 2 free parameters, and so should provide a useful framework for interpreting observations of pre-stellar cores and protostars, and for calculations of radiation transport and time-dependent chemistry. As an example, we model the pre-stellar core L1544.Comment: To appear in Astrophysical Journal, Jan 20th, 2001. 18 pages incl. 3 fig

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, $\delta_\mathrm{sol}$, that is solenoidal, as opposed to compressive; specifically we consider $\delta_\mathrm{sol}=1,\,2/3,\,1/3,\,1/9\;{\rm and}\;0$. For each value of $\delta_\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 $\delta_\mathrm{sol}(<\!1/3)$ filament fragmentation dominates and delivers relatively high mass stars. Conversely, with high values of $\delta_\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 $\delta_\mathrm{sol}$.Comment: 9 pages. Accepted by MNRA

Simulating star formation in molecular cloud cores I. The influence of low levels of turbulence on fragmentation and multiplicity

We present the results of an ensemble of simulations of the collapse and fragmentation of dense star-forming cores. We show that even with very low levels of turbulence the outcome is usually a binary, or higher-order multiple, system. We take as the initial conditions for these simulations a typical low-mass core, based on the average properties of a large sample of observed cores. All the simulated cores start with a mass of $M = 5.4 M_{\odot}$, a flattened central density profile, a ratio of thermal to gravitational energy $\alpha_{\rm therm} = 0.45$ and a ratio of turbulent to gravitational energy $\alpha_{\rm turb} = 0.05$. Even this low level of turbulence is sufficient to produce multiple star formation in 80% of the cores; the mean number of stars and brown dwarfs formed from a single core is 4.55, and the maximum is 10. At the outset, the cores have no large-scale rotation. The only difference between each individual simulation is the detailed structure of the turbulent velocity field. The multiple systems formed in the simulations have properties consistent with observed multiple systems. Dynamical evolution tends preferentially to eject lower mass stars and brown dwarves whilst hardening the remaining binaries so that the median semi-major axis of binaries formed is $\sim 30$ au. Ejected objects are usually single low-mass stars and brown dwarfs, yielding a strong correlation between mass and multiplicity. Our simulations suggest a natural mechanism for forming binary stars that does not require large-scale rotation, capture, or large amounts of turbulence.Comment: 20 pages, 12 figures submitted to A&