The formation of brown dwarfs in discs: Physics, numerics, and observations

A large fraction of brown dwarfs and low-mass stars may form by gravitational fragmentation of relatively massive (a few 0.1 Msun), extended (a few hundred AU) discs around Sun-like stars. We present an ensemble of radiative hydrodynamic simulations that examine the conditions for disc fragmentation. We demonstrate that this model can explain the low-mass IMF, the brown dwarf desert, and the binary properties of low-mass stars and brown dwarfs. Observing discs that are undergoing fragmentation is possible but very improbable, as the process of disc fragmentation is short lived (discs fragment within a few thousand years).Comment: 4 pages, for the proceedings of IAU Symposium 270: Computational Star Formation, Barcelona, 201

Can giant planets form by gravitational fragmentation of discs?

Gravitational fragmentation has been proposed as a mechanism for the formation of giant planets in close orbits around solar-type stars. However, it is debatable whether this mechanism can function in the inner regions (R<40 AU) of real discs. We use a newly developed method for treating the energy equation and the equation of state, which accounts for radiative transfer effects in SPH simulations of circumstellar discs. The different chemical and internal states of hydrogen and the properties of dust at different densities and temperatures (ice coated dust grains at low temperatures, ice melting, dust sublimation) are all taken into account by the new method.We present radiative hydrodynamic simulations of the inner regions of massive circumstellar discs and examine two cases: (i) a disc irradiated by a cool background radiation field (T_bgr=10K)and (ii) a disc heated by radiation from its central star (T_bgr~1/R). In neither case does the disc fragment: in the former because it cannot cool fast enough and in the latter because it is not gravitationally unstable. Our results (a) corroborate previous numerical results using different treatments for the hydrodynamics and the radiative transfer, and (b) confirm our own earlier analytic predictions. We conclude that disc fragmentation is unlikely to be able to produce giant planets around solar-type stars at radii <40 AU.Comment: Accepted by A&A, 10 pages, high-resolution available at http://www.astro.cf.ac.uk/pub/Dimitrios.Stamatellos/publications

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

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&

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

Brown dwarf formation by gravitational fragmentation of massive, extended protostellar discs

We suggest that low-mass hydrogen-burning stars like the Sun should sometimes form with massive extended discs; and we show, by means of radiation hydrodynamic simulations, that the outer parts of such discs (R>100 AU) are likely to fragment on a dynamical timescale (10^3 to $10^4 yr), forming low-mass companions: principally brown dwarfs (BDs), but also very low-mass hydrogen-burning stars and planetary-mass objects. A few of the BDs formed in this way remain attached to the primary star, orbiting at large radii. The majority are released into the field, by interactions amongst themselves; in so doing they acquire only a low velocity dispersion (<2 km/s), and therefore they usually retain small discs, capable of registering an infrared excess and sustaining accretion. Some BDs form close BD/BD binaries, and these binaries can survive ejection into the field. This BD formation mechanism appears to avoid some of the problems associated with the `embryo ejection' scenario, and to answer some of the questions not yet answered by the `turbulent fragmentation' scenario.Comment: 5 pages, accepted for publication in MNRAS Letter