Triple trouble for XZ Tau : deep imaging with the Jansky Very Large Array

DF gratefully acknowledges support from STFC grant ST/J001422/1. RJI acknowledges support in the form of ERC Advanced Investigator programme, cosmicism. EI acknowledges funding from CONICYT/FONDECYT postdoctoral project no.: 3130504.We present new observations of the XZ Tau system made at high angular resolution (55 mas) with the Karl G. Jansky Very Large Array (VLA) at a wavelength of 7 mm. Observations of XZ Tau made with the VLA in 2004 appeared to show a triple-star system, with XZ Tau A resolved into two sources, XZ Tau A and XZ Tau C. The angular separation of XZ Tau A and C (0.09 arcsec) suggested a projected orbital separation of around 13 au with a possible orbital period of around 40 yr. Our follow-up observations were obtained approximately 8 yr later, a fifth of this putative orbital period, and should therefore allow us to constrain the orbital parameters of XZ Tau C, and evaluate the possibility that a recent periastron passage of C coincided with the launch of extended optical outflows from XZ Tau A. Despite improved sensitivity and resolution, as compared with the 2004 observations, we find no evidence of XZ Tau C in our data. Components A and B are detected with a signal-to-noise ratio greater than 10; their orbital motions are consistent with previous studies of the system, although the emission from XZ Tau A appears to be weaker. Three possible interpretations are offered: either XZ Tau C is transiting XZ Tau A, which is broadly consistent with the periastron passage hypothesis, or the emission seen in 2004 was that of a transient, or XZ Tau C does not exist. A fourth interpretation, that XZ Tau C was ejected from the system, is dismissed due to the lack of angular momentum redistribution in the orbits of XZ Tau A and XZ Tau B that would result from such an event. Transients are rare but cannot be ruled out in a T Tauri system known to exhibit variable behaviour. Our observations are insufficient to distinguish between the remaining possibilities, at least not until we obtain further VLA observations at a sufficiently later time. A further non-detection would allow us to reject the transit hypothesis, and the periastron passage of XZ Tau C as agent of XZ Tau A's outflows.Publisher PDFPeer reviewe

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

Numerical Simulations of Supernova Dust Destruction. II. Metal-Enriched Ejecta Knots

Following our previous work, we investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We vary the relative velocity between the reverse shock and ejecta clump to explore the effects of shock-heating and cloud compression. Because supernova ejecta will be metal-enriched, we consider gas metallicities from Z/Zsun = 1 to 100 and their influence on cooling properties of the cloud and the thermal sputtering rates of embedded dust grains. We post-process the simulation output to calculate grain sputtering for a variety of species and size distributions. In the metallicity regime considered in this paper, the balance between increased radiative cooling and increased grain erosion depends on the impact velocity of the reverse shock. For slow shocks (velocity less than or equal to 3000 km/s), the amount of dust destruction is comparable across metallicities, or in some cases is decreased with increased metallicity. For higher shock velocities (velocity greater than or equal to 5000 km/s), an increase in metallicity from Z/Zsun = 10 to 100 can lead to an additional 24% destruction of the initial dust mass. While the total dust destruction varies widely across grain species and simulation parameters, our most extreme cases result in complete destruction for some grain species and only 44% dust mass survival for the most robust species. These survival rates are important in understanding how early supernovae contribute to the observed dust masses in high-redshift galaxies.Comment: 10 pages, 6 figures, 2 tables, changes made to the text and figures as suggested by the anonymous referee, accepted by the Astrophysical Journa

Evolution from protoplanetary to debris discs: The transition disc around HD 166191

HD 166191 has been identified by several studies as hosting a rare and extremely bright warm debris disc with an additional outer cool disc component. However, an alternative interpretation is that the star hosts a disc that is currently in transition between a full gas disc and a largely gas-free debris disc. With the help of new optical to mid-IR spectra and Herschel imaging, we argue that the latter interpretation is supported in several ways: i) we show that HD 166191 is co-moving with the ~4 Myr-old Herbig Ae star HD 163296, suggesting that the two have the same age, ii) the disc spectrum of HD 166191 is well matched by a standard radiative transfer model of a gaseous protoplanetary disc with an inner hole, and iii) the HD 166191 mid-IR silicate feature is more consistent with similarly primordial objects. We note some potential issues with the debris disc interpretation that should be considered for such extreme objects, whose lifetime at the current brightness is mush shorter than the stellar age, or in the case of the outer component requires a mass comparable to the solid component of the Solar nebula. These aspects individually and collectively argue that HD 166191 is a 4-5 Myr old star that hosts a gaseous transition disc. Though it does not argue in favour of either scenario, we find strong evidence for 3-5 um disc variability. We place HD 166191 in context with discs at different evolutionary stages, showing that it is a potentially important object for understanding the protoplanetary to debris disc transition.Comment: accepted to MNRAS, fixed typos in abstract and axis labe

Numerical Simulations of Supernova Dust Destruction. I. Cloud-crushing and Post-processed Grain Sputtering

We investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We then post-process the simulation output to calculate the grain sputtering for a variety of species and size distributions. We explore the parameter space appropriate for this problem by altering the over-density of the ejecta clumps and the speed of the reverse shocks. Since radiative cooling could lower the temperature of the medium in which the dust is embedded and potentially protect the dust by slowing or halting grain sputtering, we study the effects of different cooling methods over the time scale of the simulations. In general, our results indicate that grains with radii less than 0.1 microns are sputtered to much smaller radii and often destroyed completely, while larger grains survive their interaction with the reverse shock. We also find that, for high ejecta densities, the percentage of dust that survives is strongly dependent on the relative velocity between the clump and the reverse shock, causing up to 50% more destruction for the highest velocity shocks. The fraction of dust destroyed varies widely across grain species, ranging from total destruction of Al2O3 grains to minimal destruction of Fe grains (only 20% destruction in the most extreme cases). C and SiO2 grains show moderate to strong sputtering as well, with 38% and 80% mass loss. The survival rate of grains formed by early supernovae is crucial in determining whether or not they can act as the "dust factories" needed to explain high-redshift dust.Comment: 18 pages, 16 figures, 2 tables, changes make to the text and figures as suggested by the anonymous referee, accepted by the Astrophysical Journa

AKARI and BLAST Observations of the Cassiopeia A Supernova Remnant and Surrounding Interstellar Medium

We use new large area far infrared maps ranging from 65 - 500 microns obtained with the AKARI and the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) missions to characterize the dust emission toward the Cassiopeia A supernova remnant (SNR). Using the AKARI high resolution data we find a new "tepid" dust grain population at a temperature of ~35K and with an estimated mass of 0.06 solar masses. This component is confined to the central area of the SNR and may represent newly-formed dust in the unshocked supernova ejecta. While the mass of tepid dust that we measure is insufficient by itself to account for the dust observed at high redshift, it does constitute an additional dust population to contribute to those previously reported. We fit our maps at 65, 90, 140, 250, 350, and 500 microns to obtain maps of the column density and temperature of "cold" dust (near 16 K) distributed throughout the region. The large column density of cold dust associated with clouds seen in molecular emission extends continuously from the surrounding interstellar medium to project on the SNR, where the foreground component of the clouds is also detectable through optical, X-ray, and molecular extinction. At the resolution available here, there is no morphological signature to isolate any cold dust associated only with the SNR from this confusing interstellar emission. Our fit also recovers the previously detected "hot" dust in the remnant, with characteristic temperature 100 K.Comment: Accepted for publication in the Astrophysical Journal. Maps and related data are available at http://blastexperiment.info

Herschel ATLAS : the cosmic star formation history of quasar host galaxies

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0 z 6 and absolute I-band magnitudes -22 > I-AB > -32 We use the science demonstration observations of the first similar to 16 deg(2) from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at > 5 sigma at 250, 350 and 500 mu m. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22 > I-AB > -24 have a comoving star formation rate (derived from 100 mu m rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I-AB -26 have a maximum contribution to the star formation density at z similar to 3. The volume-averaged star formation rate of -22 > IAB > -24 quasars evolves as (1 + z)(2.3 +/- 0.7) at z 2, but the evolution at higher luminosities is much faster reaching (1 + z)(10 +/- 1) at -26 > I-AB > -28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation

Herschel observations of embedded protostellar clusters in the Rosette Molecular Cloud

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.Comment: Astronomy & Astrophysics letter, 6 pages, 4 figures, accepted for publication in the Special Issue for Herschel first result