Researcher profile

D. Stamatellos

D. Stamatellos contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2015arXiv

The dynamical evolution of low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects formed through disc fragmentation

Theory and simulations suggest that it is possible to form low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects via disc fragmentation. As disc fragmentation results in the formation of several bodies at comparable distances to the host star, their orbits are generally unstable. Here, we study the dynamical evolution of these objects. We set up the initial conditions based on the outcomes of the SPH simulations of Stamatellos & Whitworth, and for comparison we also study the evolution of systems resulting from lower-mass fragmenting discs. We refer to these two sets of simulations as set 1 and set 2. At 10 Myr, approximately half of the host stars have one companion left, and approximately 22% (set 1) to 9.8% (set 2) of the host stars are single. Systems with multiple secondaries in relatively stable configurations are common (about 30% and 44%, respectively). The majority of the companions are ejected within 1 Myr with velocities mostly below 5 km/s, with some runaway escapers with velocities over 30 km/s. About 6% (set 1) and 2% (set 2) of the companions pair up into very low-mass binary systems. The majority of these pairs escape as very low-mass binaries, while others remain bound to the host star in hierarchical configurations (often with retrograde inner orbits). Physical collisions with the host star (0.43 and 0.18 events per host star for set 1 and set 2) and between companions (0.08 and 0.04 events per host star for set 1 and set 2) are relatively common and their frequency increases with increasing disc mass. Our study predicts observable properties of very low-mass binaries, low-mass hierarchical systems, the brown dwarf desert, and free-floating brown dwarfs and planetary-mass objects in and near young stellar groupings, which can be used to distinguish between different formation scenarios of very low-mass stars, brown dwarfs and planetary-mass objects.

preprint2014arXiv

Properties of Starless and Prestellar Cores in Taurus Revealed by Herschel SPIRE/PACS Imaging

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $μ$m, 160 $μ$m, 250 $μ$m, 350 $μ$m and 500 $μ$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

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.

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.

preprint2013arXiv

Reconstructing the density and temperature structure of prestellar cores from $Herschel$ data: A case study for B68 and L1689B

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $κ_λ =0.1(λ/300 μm)^{-2}$, agrees to within 50% with the dust extinction constraints

preprint2012arXiv

A search for pre-substellar cores and proto-brown dwarf candidates in Taurus: multiwavelength analysis in the B213-L1495 clouds

In an attempt to study whether the formation of brown dwarfs (BDs) takes place as a scaled-down version of low-mass stars, we conducted IRAM30m/MAMBO-II observations at 1.2 mm in a sample of 12 proto-BD candidates selected from Spitzer/IRAC data in the B213-L1495 clouds in Taurus. Subsequent observations with the CSO at 350 micron, VLA at 3.6 and 6 cm, and IRAM30m/EMIR in the 12CO(1-0), 13CO(1-0), and N2H+(1-0) transitions were carried out toward the two most promising Spitzer/IRAC source(s), J042118 and J041757. J042118 is associated with a compact (<10 arcsec or <1400 AU) and faint source at 350 micron, while J041757 is associated with a partially resolved (~16 arcsec or ~2000 AU) and stronger source emitting at centimetre wavelengths with a flat spectral index. The corresponding masses of the dust condensations are ~1 and ~5 Mjup for J042118 and J041757, respectively. In addition, about 40 arcsec to the northeast of J041757 we detect a strong and extended submillimetre source, J041757-NE, which is not associated with NIR/FIR emission down to our detection limits, but is clearly detected in 13CO and N2H+ at ~7 km/s, and for which we estimated a total mass of ~100 Mjup, close to the mass required to be gravitationally bound. In summary, our observational strategy has allowed us to find in B213-L1495 two proto-BD candidates and one pre-substellar core candidate, whose properties seem to be consistent with a scaled-down version of low-mass stars.

preprint2012arXiv

Cores in Infra-Red Dark Clouds (IRDCs) seen in the Hi-GAL survey between l = 300° and l = 330°

We have used data taken as part of the Herschel infrared Galactic Plane survey (Hi-GAL) to study 3171 infrared-dark cloud (IRDC) candidates that were identified in the mid-infrared (8 μm) by Spitzer (we refer to these as &#39;Spitzer-dark&#39; regions). They all lie in the range l=300 - 330 \circ and |b| 6 1 \circ. Of these, only 1205 were seen in emission in the far-infrared (250-500 μm) by Herschel (we call these &#39;Herschel-bright&#39; clouds). It is predicted that a dense cloud will not only be seen in absorption in the mid-infrared, but will also be seen in emission in the far-infrared at the longest Herschel wavebands (250-500 μm). If a region is dark at all wavelengths throughout the mid-infrared and far-infrared, then it is most likely to be simply a region of lower background infrared emission (a &#39;hole in the sky&#39;). Hence, it appears that previous surveys, based on Spitzer and other mid-infrared data alone, may have over-estimated the total IRDC population by a factor of 2. This has implications for estimates of the star formation rate in IRDCs in the Galaxy.We studied the 1205 Herschel-bright IRDCs at 250 μm, and found that 972 of them had at least one clearly defined 250-μm peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-μm point source somewhere within the cloud, 149 (15 per cent) contained a 24-μm point source but no 8-μm source, and 170 (18 per cent) contained no 24-μm or 8-μm point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

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.

preprint2012arXiv

Isolated starless cores in IRDCs in the Hi-GAL survey

In a previous paper we identified cores within infrared dark clouds (IRDCs). We regarded those without embedded sources as the least evolved, and labelled them starless. Here we identify the most isolated starless cores and model them using a three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive the cores&#39; physical parameters and discuss the relation between the mass, temperature, density, size and the surrounding interstellar radiation field (ISRF) for the cores. The masses of the cores were found not to correlate with their radial size or central density. The temperature at the surface of a core was seen to depend almost entirely on the level of the ISRF surrounding the core. No correlation was found between the temperature at the centre of a core and its local ISRF. This was seen to depend, instead, on the density and mass of the core.

preprint2012arXiv

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&#34; 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.

preprint2011arXiv

Disc Frequencies for Brown Dwarfs in the Upper Scorpius OB Association: Implications for Brown Dwarf Formation Theories

We have investigated the brown dwarf (BD) and stellar disc fractions in the Upper Scorpius OB Association (USco) and compared them with several other young regions. We have compiled the most complete sample of of all spectroscopically confirmed BDs in USco, and have made use of the WISE catalog to identify the disc candidates. We report on the discovery of 12 new BD discs in USco, with spectral type (SpT) between M6 and M8.5. The WISE colors for the new discs are similar to the primordial (transition) discs earlier detected in USco. Combining with previous surveys, we find the lowest inner disc fractions (~20-25%) for a wide range in stellar masses (~0.01-4.0 Msun) in the USco association. The low disc fractions for high-mass stars in USco (and the other clusters) are consistent with an evolutionary decline in inner disc frequency with age. However, BD disc fractions are higher than those for the stars in 1-3 Myr clusters, but very low in the ~5 Myr old USco. Also, primordial BD discs are still visible in the ~10 Myr TW Hydrae association, whereas the higher mass stars have all transitioned to the debris stage by this age. The disc frequencies for BDs do not show any dependence on the stellar density or the BD/star number ratio in a cluster. We suggest that the large differences in the observed BD disc fractions between regions may well be due to different BD formation mechanisms and therefore different initial disc fractions/properties. We also present a WISE SED classification scheme, based on the Ks and WISE bands of 3.4-12micron, to distinguish between the Class I/II and the Class III sequence. Our work includes a comparison of the sensitivities of WISE and Spitzer disc surveys. We estimate that WISE can be incomplete for discs at SpT later than M8 in distant clusters such as SOri. WISE should be able to recover the M8-M9 discs in the nearby young clusters.

preprint2011arXiv

The initial conditions of high-mass star formation: radiative transfer models of IRDCs seen in the Herschel Hi-GAL survey

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi-wavelength, spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G031.03+00.26 and G031.03+00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We found masses ranging from 90 to 290 solar masses with temperatures from 8 to 11K at the centre of each core and 18 to 28K at the surface. The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high mass star formation having yet occurred in three of our cores.

preprint2010arXiv

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call &#39;cores&#39; in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

preprint2010arXiv

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.

preprint2010arXiv

Toward understanding the formation of multiple systems - A pilot IRAM-PdBI survey of Class 0 objects

The formation process of binary stars and multiple systems is poorly understood. Here, we seek to determine the typical outcome of protostellar collapse and to constrain models of binary formation by core fragmentation during collapse, using high-resolution millimeter continuum imaging of very young (Class 0) protostars observed at the beginning of the main accretion phase. We carried out a pilot high-resolution study of 5 Class 0 objects, using the most extended (A) configuration of the IRAM PdBI at 1.3 mm, which allow us to probe the multiplicity of Class 0 protostars down to separations a ~50 AU and circumstellar mass ratios q ~0.07. We show that our PdBI observations revealed only wide (>1500 AU) protobinary systems and/or outflow-generated features. When combined with previous millimeter interferometric observations of Class 0 protostars, our pilot PdBI study tentatively suggests that the binary fraction in the ~ 75-1000 AU range increases from the Class 0 to the Class I stage. It also seems to argue against purely hydrodynamic models of binary star formation. We briefly discuss possible alternative scenarios to reconcile the low multiplicity rate of Class 0 protostars on small scales with the higher binary fraction observed at later evolutionary stages.

preprint2009arXiv

The initial conditions of isolated star formation -- IX. Akari mapping of an externally heated pre-stellar core

We present observations of L1155 and L1148 in the Cepheus molecular cloud, taken using the FIS instrument on the Akari satellite. We compare these data to submillimetre data taken using the SCUBA camera on the JCMT, and far-infrared data taken with the ISOPHOT camera on board the ISO satellite. All of the data show a relation between the position of the peak of emission and the wavelength for the core of L1155. We interpret this as a temperature gradient. We fit modified blackbody curves to the spectral energy distributions at two positions in the core and see that the central core in L1155 (L1155C) is approximately 2 degrees warmer at one edge than it is in the centre. We consider a number of possible heating sources and conclude that the A6V star BD+67 1263 is the most likely candidate. This star is at a distance of 0.7 pc from the front of L1155C in the plane of the sky. We carry out radiative transfer modelling of the L1155C core including the effects from the nearby star. We find that we can generate a good fit to the observed data at all wavelengths, and demonstrate that the different morphologies of the core at different wavelengths can be explained by the observed 2 degree temperature gradient. The L1148 core exhibits a similar morphology to that of L1155C, and the data are also consistent with a temperature gradient across the core. In this case, the most likely heating source is the star BD197053. Our findings illustrate very clearly that the apparent observed morphology of a pre-stellar core can be highly dependent on the wavelength of the observation, and that temperature gradients must be taken into account before converting images into column density distributions. This is important to note when interpreting Akari and Spitzer data and will also be significant for Herschel data.

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.