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&

Controlling Artificial Viscosity in SPH simulations of accretion disks

We test the operation of two methods for selective application of Artificial Viscosity (AV) in SPH simulations of Keplerian Accretion Disks, using a ring spreading test to quantify effective viscosity, and a correlation coefficient technique to measure the formation of unwanted prograde alignments of particles. Neither the Balsara Switch nor Time Dependent Viscosity work effectively, as they leave AV active in areas of smooth shearing flow, and do not eliminate the accumulation of alignments of particles in the prograde direction. The effect of both switches is periodic, the periodicity dependent on radius and unaffected by the density of particles. We demonstrate that a very simple algorithm activates AV only when truly convergent flow is detected and reduces the unwanted formation of prograde alignments. The new switch works by testing whether all the neighbours of a particle are in Keplerian orbit around the same point, rather than calculating the divergence of the velocity field, which is very strongly affected by Poisson noise in the positions of the SPH particles.Comment: 8 pages, 5 figure

A natural formation scenario for misaligned and short-period eccentric extrasolar planets

Recent discoveries of strongly misaligned transiting exoplanets pose a challenge to the established planet formation theory which assumes planetary systems to form and evolve in isolation. However, the fact that the majority of stars actually do form in star clusters raises the question how isolated forming planetary systems really are. Besides radiative and tidal forces the presence of dense gas aggregates in star-forming regions are potential sources for perturbations to protoplanetary discs or systems. Here we show that subsequent capture of gas from large extended accretion envelopes onto a passing star with a typical circumstellar disc can tilt the disc plane to retrograde orientation, naturally explaining the formation of strongly inclined planetary systems. Furthermore, the inner disc regions may become denser, and thus more prone to speedy coagulation and planet formation. Pre-existing planetary systems are compressed by gas inflows leading to a natural occurrence of close-in misaligned hot Jupiters and short-period eccentric planets. The likelihood of such events mainly depends on the gas content of the cluster and is thus expected to be highest in the youngest star clusters.Comment: 7 pages, 4 figures. Accepted for publication in MNRAS. Updated to match published versio

Simulating star formation in molecular cloud cores IV. The role of turbulence and thermodynamics

We perform SPH simulations of the collapse and fragmentation of low-mass cores having different initial levels of turbulence (alpha_turb=0.05,0.10,0.25). We use a new treatment of the energy equation which captures the transport of cooling radiation against opacity due to both dust and gas (including the effects of dust sublimation, molecules, and H^- ions). We also perform comparison simulations using a standard barotropic equation of state. We find that -- when compared with the barotropic equation of state -- our more realistic treatment of the energy equation results in more protostellar objects being formed, and a higher proportion of brown dwarfs; the multiplicity frequency is essentially unchanged, but the multiple systems tend to have shorter periods (by a factor ~3), higher eccentricities, and higher mass ratios. The reason for this is that small fragments are able to cool more effectively with the new treatment, as compared with the barotropic equation of state. We find that the process of fragmentation is often bimodal. The first protostar to form is usually, at the end, the most massive, i.e. the primary. However, frequently a disc-like structure subsequently forms round this primary, and then, once it has accumulated sufficient mass, quickly fragments to produce several secondaries. We believe that this delayed fragmentation of a disc-like structure is likely to be an important source of very low-mass hydrogen-burning stars and brown dwarfs.Comment: 14 pages, 8 figures. Accepted for publication by A&

Revised analysis of SPIRE observations for 2M1207

We have revised our analysis of the SPIRE observations of 2MASSW J1207334-393254 (2M1207). Recent PACS observations show a bright source located ~25" east of 2M1207. There are issues in terms of the detection/non-detection of the bright source when comparing the Spitzer, WISE, and PACS observations. It is apparently inconsistent, perhaps due to variability or low signal-to-noise of the data. We have looked into the possible misidentification of the target, and have revised the measured SPIRE fluxes and the disc parameters for 2M1207. We have also reviewed which among the various formation mechanisms of this system would still be valid.Comment: Revised SPIRE fluxe

Modelling Herschel observations of infrared-dark clouds in the Hi-GAL survey

We demonstrate the use of the 3D Monte Carlo radiative transfer code PHAETHON to model infrared-dark clouds (IRDCs) that are externally illuminated by the interstellar radiation field (ISRF). These clouds are believed to be the earliest observed phase of high-mass star formation, and may be the high-mass equivalent of lower-mass prestellar cores. We model three different cases as examples of the use of the code, in which we vary the mass, density, radius, morphology and internal velocity field of the IRDC. We show the predicted output of the models at different wavelengths chosen to match the observing wavebands of Herschel and Spitzer. For the wavebands of the long- wavelength SPIRE photometer on Herschel, we also pass the model output through the SPIRE simulator to generate output images that are as close as possible to the ones that would be seen using SPIRE. We then analyse the images as if they were real observations, and compare the results of this analysis with the results of the radiative transfer models. We find that detailed radiative transfer modelling is necessary to accurately determine the physical parameters of IRDCs (e.g. dust temperature, density profile). This method is applied to study G29.55+00.18, an IRDC observed by the Herschel Infrared Galactic Plane survey (Hi-GAL), and in the future it will be used to model a larger sample of IRDCs from the same survey.Comment: MNRAS accepted, High resolution paper available at http://www.astro.cardiff.ac.uk/pub/Dimitrios.Stamatellos/Publications.htm

How to identify the youngest protostars

We study the transition from a prestellar core to a Class 0 protostar, using SPH to simulate the dynamical evolution, and a Monte Carlo radiative transfer code to generate the SED and isophotal maps. For a prestellar core illuminated by the standard interstellar radiation field, the luminosity is low and the SED peaks at ~190 micron. Once a protostar has formed, the luminosity rises (due to a growing contribution from accretion onto the protostar) and the peak of the SED shifts to shorter wavelengths (~80-100 micron). However, by the end of the Class 0 phase, the accretion rate is falling, the luminosity has decreased, and the peak of the SED shifts back towards longer wavelengths (90-150 micron). In our simulations, the density of material around the protostar remains sufficiently high well into the Class 0 phase that the protostar only becomes visible in the NIR if it is displaced from the centre dynamically. Raw submm/mm maps of Class 0 protostars tend to be much more centrally condensed than those of prestellar cores. However, when convolved with a typical telescope beam, the difference in central concentration is less marked, although the Class 0 protostars appear more circular. Our results suggest that, if a core is deemed to be prestellar on the basis of having no associated IRAS source, no cm radio emission, and no outflow, but it has a circular appearance and an SED which peaks at wavelengths below ~170 micron, it may well contain a very young Class 0 protostar.Comment: Accepted by A&A (avaliable with high-res images at http://carina.astro.cf.ac.uk/pub/Dimitrios.Stamatellos/publications

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.Comment: 10 pages, 8 figures, Accepted for publication in Astrophysics and Space Scienc

A VLA search for young protostars embedded in dense cores

Four dense cores, L1582A, L1689A, B133 and B68, classified as prestellar in terms of the absence of detectable NIR emission, are observed at radio wavelengths to investigate whether they nurture very young protostars. No definite young protostars were discovered in any of the four cores observed. A few radio sources were discovered close to the observed cores, but these are most likely extragalactic sources or YSOs unrelated to the cores observed. In L1582A we discovered a weak radio source near the centre of the core with radio characteristics and offset from the peak of the submillimeter emission similar to that of the newly discovered protostar in the core L1014, indicating a possible protostellar nature for this source. This needs to be confirmed with near- and/or mid-infrared observations (e.g. with Spitzer). Hence based on the current observations we are unable to confirm unequivocally that L1582A is starless. In L1689A a possible 4.5-sigma radio source was discovered at the centre of the core, but needs to be confirmed with future observations. In B133 a weak radio source, possibly a protostar, was discovered at the edge of the core on a local peak of the core submm emission, but no source was detected at the centre of the core. Thus, B133 is probably starless, but may have a protostar at its edge. In B68 no radio sources were discovered inside or at the edge of the core, and thus B68 is indeed starless. Four more radio sources with spectral indices characteristic of young protostars were discovered outside the cores but within the extended clouds in which these cores reside. Conclusions: We conclude that the number of cores misclassified as prestellar is probably very small and does not significantly alter the estimated lifetime of the prestellar phase.Comment: Accepted by A&