Sub-mm Data Research Articles

Census of Young Stars in Galactic Star Forming Region Sh2-88

We present here the census and characterization of Young Stellar Objects (YSOs) associated with the active star-forming H II region Sh2-88 with multi-wavelength data situated in the Vulpecula OB1 association. To understand the extinction properties of Sh2-88 region, we generate an extinction map over 20×20 arcminutes using the neighbourhood method of Gutermuth et al. (2005) from the deep near-infrared HK band photometric archival data of UKIDSS. We find a sizable variable extinction over the region up to AV ≃11.7 mag, with a median value of AV ≃4.3 mag. The high extinction regions trace the dense filamentary structure, which was found earlier in the far-infrared and sub-mm data. Following the YSO classification scheme of Gutermuth et al. (2009), we have identified 74 YSOs in that region using Spitzer-IRAC data from four mid-infrared bands (3.6, 4.5, 5.8 and 8 μm) and 2MASS near-infrared K-band (2.16 μm) data. Among 74 YSOs, six are class I and 68 class II objects. Using IPHAS Hα, r, and i-band photometry, we have identified 41 additional YSOs. Among these 41, four sources are common with IR excess YSOs. From the spatial distribution of YSOs, it is apparent that YSOs are primarily distributed over the relatively higher extinction regions of the cloud, and star-formation activity might be ongoing in those dense parts of the cloud.

The Connection between Different Tracers of the Diffuse Interstellar Medium: Kinematics

Abstract Using visible, radio, microwave, and submillimeter data, we study several lines of sight toward stars generally closer than 1 kpc on a component-by-component basis. We derive the component structure seen in absorption at visible wavelengths from Ca ii, Ca i, K i, CH, CH+, and CN and compare it to emission from H i, CO and its isotopologues, and C+ from the GOT C+ survey. The correspondence between components in emission and absorption helps create a more unified picture of diffuse atomic and molecular gas in the interstellar medium. We also discuss how these tracers are related to the CO-dark H2 gas probed by C+ emission and discuss the kinematic connections among the species observed.

ALMA + VLT observations of a damped Lyman-α absorbing galaxy: massive, wide CO emission, gas-rich but with very low SFR

We are undertaking an ALMA survey of molecular gas in galaxies selected for their strong HI absorption, so-called DLA/sub-DLA galaxies. Here we report CO(2-1) detection from a DLA galaxy at z = 0.716. We also present optical and near-infrared spectra of the galaxy revealing [OII], H{\alpha} and [NII] emission lines shifted by ~170 km/s relative to the DLA, and providing an oxygen abundance 3.2 times solar, similar to the absorption metallicity. We report low unobscured SFR ~1 Msun/yr given the large reservoir of molecular gas, and also modest obscured SFR=4.5(+4.4,-2.6) Msun/yr based on far-IR and sub-mm data. We determine mass components of the galaxy: log[M*/Msun] = 10.80(+0.07,-0.14), log[M mol-gas/Msun] = 10.37 +/-0.04, and log[M dust/Msun] = 8.45(+0.10,-0.30). Surprisingly, this HI absorption-selected galaxy has no equivalent objects in CO surveys of flux-selected samples. The galaxy falls off current scaling relations for the SFR to molecular gas mass and CO Tully-Fisher relation. Detailed comparison of kinematical components of the absorbing, ionized and molecular gas, combined with their spatial distribution, suggests that part of the CO gas is both kinematically and spatially de-coupled from the main galaxy. It is thus possible that a major star burst in the past could explain the wide CO profile as well as the low SFR. Support for this also comes from the SED favouring an instantaneous burst of age ~0.5 Gyr. Our survey will establish whether flux-selected surveys of molecular gas are missing a key stage in the evolution of galaxies and their conversion of gas to stars.

Submillimeter and radio variability of Sagittarius A*

We report on a statistical analysis of the 345 GHz submillimeter (submm) and 100 GHz radio flux density distribution of Sagittarius A* (Sgr A*). The submm data set consists of 345 GHz data obtained from different Large Apex Bolometer Camera (LABOCA) campaigns between 2008 and 2014, and additional literature data from 2004 to 2009 at comparable wavelengths. The radio observations were carried out with the Australia Telescope Compact Array (ATCA) between 2010 and 2014. We used a combined maximum likelihood estimator (MLE) and Kolmogorov-Smirnov (KS) statistics method to test for a possible power-law distribution in the high flux density excursions (flares) at both wavebands. We find that both flux density distributions can be described by a shifted power-law of the form p (x ) ∝ (x − s )− α with α ~ 4 (submm: α = 4.0 ± 1.7; radio: α = 4.7 ± 0.8). The same power-law index was previously found for the near-infrared (NIR) flux density distribution. These results may strengthen our preferred flare emission model: a combined synchrotron self-Compton (SSC) and adiabatically expanding self-absorbed synchrotron blob model where the flaring activity across all wavebands stem from the same source components and the variable emission can be described by a single state red noise process. Within the framework of the expanding blob model the similarity of the radio and the submm flux density distribution may also narrow down possible initial synchrotron turnover frequencies ν 0 to be mainly around 350 GHz and possible expansion velocities v exp to be predominantly around 0.01 c.

Constraints on AGN feedback from its Sunyaev–Zel'dovich imprint on the cosmic background radiation

We derive constraints on feedback by active galactic nuclei (AGN) by setting limits on their thermal Sunyaev-Zel'dovich (SZ) imprint on the cosmic microwave background (CMB). The amplitude of any SZ signature is small and degenerate with the poorly known sub-mm spectral energy distribution of the AGN host galaxy and other unresolved dusty sources along the line of sight. Here we break this degeneracy by combining microwave and sub-mm data from Planck with all-sky far-infrared maps from the AKARI satellite. We first test our measurement pipeline using the Sloan Digital Sky Survey (SDSS) redMaPPer catalogue of galaxy clusters, finding a highly significant detection ($>$$20\sigma$) of the SZ effect together with correlated dust emission. We then constrain the SZ signal associated with spectroscopically confirmed quasi-stellar objects (QSOs) from SDSS data release 7 (DR7) and the Baryon Oscillation Spectroscopic Survey (BOSS) DR12. We obtain a low-significance ($1.6\sigma$) hint of an SZ signal, pointing towards a mean thermal energy of $\simeq 5 \times 10^{60}$ erg, lower than reported in some previous studies. A comparison of our results with high-resolution hydrodynamical simulations including AGN feedback suggests QSO host masses of $M_{200c} \sim 4 \times 10^{12}~h^{-1}M_\odot$, but with a large uncertainty. Our analysis provides no conclusive evidence for an SZ signal specifically associated with AGN feedback.

Dipper discs not inclined towards edge-on orbits

Abstract The so-called dipper stars host circumstellar discs and have optical and infrared light curves that exhibit quasi-periodic or aperiodic dimming events consistent with extinction by transiting dusty structures orbiting in the inner disc. Most of the proposed mechanisms explaining the dips – i.e. occulting disc warps, vortices, and forming planetesimals – assume nearly edge-on viewing geometries. However, our analysis of the three known dippers with publicly available resolved sub-mm data reveals discs with a range of inclinations, most notably the face-on transition disc J1604-2130 (EPIC 204638512). This suggests that nearly edge-on viewing geometries are not a defining characteristic of the dippers and that additional models should be explored. If confirmed by further observations of more dippers, this would point to inner disc processes that regularly produce dusty structures far above the outer disc mid-plane in regions relevant to planet formation.

H-ATLAS: a candidate high redshift cluster/protocluster of star-forming galaxies

We investigate the region around the Planck-detected z=3.26 gravitationally lensed galaxy HATLAS J114637.9-001132 (hereinafter HATLAS12-00) using both archival Herschel data from the H-ATLAS survey and using submm data obtained with both LABOCA and SCUBA2. The lensed source is found to be surrounded by a strong overdensity of both Herschel-SPIRE sources and submm sources. We detect 17 bright (S_870 >~7 mJy) sources at >4sigma closer than 5 arcmin to the lensed object at 850/870 microns. Ten of these sources have good cross-identifications with objects detected by Herschel-SPIRE which have redder colours than other sources in the field, with 350 micron flux > 250 micron flux, suggesting that they lie at high redshift. Submillimeter Array (SMA) observations localise one of these companions to ~1 arcsecond, allowing unambiguous cross identification with a 3.6 and 4.5 micron Spitzer source. The optical/near-IR spectral energy distribution (SED) of this source is measured by further observations and found to be consistent with z>2, but incompatible with lower redshifts. We conclude that this system may be a galaxy cluster/protocluster or larger scale structure that contains a number of galaxies undergoing starbursts at the same time.

Herschel Hi-GAL imaging of massive young stellar objects

We used Herschel Hi-GAL survey data to determine whether massive young stellar objects (MYSOs) are resolved at 70$\mu$m and to study their envelope density distribution. Our analysis of three relatively isolated sources in the l=30{\deg} and l=59{\deg} Galactic fields show that the objects are partially resolved at 70$\mu$m. The Herschel Hi-GAL survey data have a high scan velocity which makes unresolved and partially resolved sources appear elongated in the 70$\mu$m images. We analysed the two scan directions separately and examine the intensity profile perpendicular to the scan direction. Spherically symmetric radiative transfer models with a power law density distribution were used to study the circumstellar matter distribution. Single dish sub-mm data were also included to study how different spatial information affects the fitted density distribution. The density distribution which best fits both the 70$\mu$m intensity profile and SED has an average index of ~0.5. This index is shallower than expected and is probably due to the dust emission from bipolar outflow cavity walls not accounted for in the spherical models. We conclude that 2D axisymmetric models and Herschel images at low scan speeds are needed to better constrain the matter distribution around MYSOs.

Conditions for star formation in nearby AGN and QSO hosts observed with near-infrared integral-field spectroscopy

Integral-field spectroscopy in the near-infrared (NIR) is a powerful tool to analyze the gaseous and stellar distributions and kinematics, as well as the excitation mechanisms in the centers of galaxies. The unique combination of NIR and sub-mm data at comparable high angular resolution, which has just been possible with SINFONI and ALMA, allows to trace warm and cold gas reservoirs. Only the NIR gives an unobscured view to the center and allows to study the conditions and impact of star formation in the centers of galaxies in a spatially resolved way. Here, we present recent studies of nearby Seyferts and low-luminosity QSOs performed by our group.

Black hole accretion preferentially occurs in gas-rich galaxies*

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$\sigma$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation.

THE MYSTERIOUS SICKLE OBJECT IN THE CARINA NEBULA: A STELLAR WIND INDUCED BOW SHOCK GRAZING A CLUMP?

Optical and near-infrared images of the Carina Nebula show a peculiar arc-shaped feature, which we call the "Sickle", next to the B-type star Trumpler 14 MJ 218. We use multi-wavelength observations to explore and constrain the nature and origin of the nebulosity. Using sub-mm data from APEX/LABOCA as well as Herschel far-infrared maps, we discovered a dense, compact clump with a mass of ~ 40 Msun located close to the apex of the Sickle. We investigate how the B-star MJ 218, the Sickle, and the clump are related. Our numerical simulations show that, in principle, a B-type star located near the edge of a clump can produce a crescent-shaped wind shock front, similar to the observed morphology. However, the observed proper motion of MJ 218 suggest that the star moves with high velocity (~ 100 km/s) through the ambient interstellar gas. We argue that the star is just about to graze along the surface of the clump, and the Sickle is a bow shock induced by the stellar wind, as the object moves supersonically through the density gradient in the envelope of the clump.

Triggered/sequential star formation? A multi-phase ISM study around the prominent IRDC G18.93-0.03

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

Constraining clumpy dusty torus models using optimized filter sets

Recent success in explaining several properties of the dusty torus around the central engine of active galactic nuclei has been gathered with the assumption of clumpiness. The properties of such clumpy dusty tori can be inferred by analyzing spectral energy distributions (SEDs), sometimes with scarce sampling given that large aperture telescopes and long integration times are needed to get good spatial resolution and signal. We aim at using the information already present in the data and the assumption of clumpy dusty torus, in particular, the CLUMPY models of Nenkova et al., to evaluate the optimum next observation such that we maximize the constraining power of the new observed photometric point. To this end, we use the existing and barely applied idea of Bayesian adaptive exploration, a mixture of Bayesian inference, prediction and decision theories. The result is that the new photometric filter to use is the one that maximizes the expected utility, which we approximate with the entropy of the predictive distribution. In other words, we have to sample where there is larger variability in the SEDs compatible with the data with what we know of the model parameters. We show that Bayesian adaptive exploration can be used to suggest new observations, and ultimately optimal filter sets, to better constrain the parameters of the clumpy dusty torus models. In general, we find that the region between 10 and 200 um produces the largest increase in the expected utility, although sub-mm data from ALMA also prove to be useful. It is important to note that here we are not considering the angular resolution of the data, which is key when constraining torus parameters. Therefore, the expected utilities derived from this methodology must be weighted with the spatial resolution of the data.

PATIAL DISTRIBUTION OF STAR FORMATION ACTIVITY ON NGC 253 BY FIR AND RADIO EMISSION LINES

The aim of this research is to reveal the spatial distribution of the star formation activity of nearby galaxies by comparing CO molecular emission lines with the large area observation in far-infrared (FIR) lines. We report the imaging observations of NGC 253 by FIR forbidden lines via FIS-FTS and CO molecular lines from low to high excitation levels with ASTE, which are good tracers of star forming regions or photo-dissociation regions, especially spiral galaxies, in order to derive the information of the physical conditions of the ambient interstellar radiation fields. The combination of spatially resolved FIR and sub-mm data leads to the star formation efficiency within galaxy. The ratio between the FIR luminosity and molecular gas mass, <TEX>$L_{FIR}/M_{H_2}$</TEX>, is expected to be proportional to the number of stars formed in the galaxy per unit molecular gas mass and time. Moreover the FIR line ux shows current star formation activity directly. Furthermore these can be systematic and statistical data for star formation history and evolution of spiral galaxies.

ANHSTIMAGING SURVEY OF LOW-MASS STARS IN THE CHAMAELEON I STAR-FORMING REGION

We present new HST/WFPC2 observations of 20 fields centered around T Tauri stars in the Chamaeleon I star forming region. Images have been obtained in the F631N ([OI]6300A), F656N (Ha) and F673N ([SII]6716A+6731A) narrow-band filters, plus the Johnson V-band equivalent F547M filter. We detect 31 T Tauri stars falling within our fields. We discuss the optical morphology of 10 sources showing evidence of either binarity, circumstellar material, or mass loss. We supplement our photometry with a compilation of optical, infrared and sub-millimeter data from the literature, together with new sub-mm data for three objects, to build the Spectral Energy Distributions (SED) of 19 single sources. Using an SED model fitting tool, we self-consistently estimate a number of stellar and disk parameters, while mass accretion rates are directly derived from our Ha photometry. We find that bolometric luminosities derived from dereddened optical data tend to be underestimated in systems with high alpha(2-24} IR spectral index, suggesting that disks seen nearly edge-on may occasionally be interpreted as low luminosity (and therefore more evolved) sources. On the other hand, the same alpha(2-24) spectral index, a tracer of the amount of dust in the warmer layers of the circumstellar disks, and the mass accretion rate appear to decay with the isocronal stellar age, suggesting that the observed age spread (~0.5-5 Myr) within the cluster is real. Our sample contains a few outliers that may have dissipated their circumstellar disks on shorter time-scale.

Confronting a solar irradiance reconstruction with solar and stellar data

Context. A recent paper by Shapiro and colleagues (2011, A&A, 529, A67) reconstructs spectral and total irradiance variations of the Sun during the holocene. Aims. In this note, we comment on why their methodology leads to large (0.5%) variations in the solar TSI on century-long time scales, in stark contrast to other reconstructions which have ≲ 0.1% variations.Methods. We examine the amplitude of the irradiance variations from the point of view of both solar and stellar data.Results. Shapiro et al.’s large amplitudes arise from differences between the irradiances computed from models A and C of Fontenla and colleagues, and from their explicit assumption that the radiances of the quiet Sun vary with the cosmic ray modulation potential. We suggest that the upper photosphere, as given by model A, is too cool, and discuss relative contributions of local vs. global dynamos to the magnetism and irradiance of the quiet Sun. We compare the slow (>22 yr) components of the irradiance reconstructions with secular changes in stellar photometric data that span 20 years or less, and find that the Sun, if varying with such large amplitudes, would still lie within the distribution of stellar photometric variations measured over a 10−20 year period. However, the stellar time series are individually too short to see if the reconstructed variations will remain consistent with stellar variations when observed for several decades more.Conclusions. By adopting model A, Shapiro et al. have over-estimated quiet-Sun irradiance variations by about a factor of two, based upon a re-analysis of sub-mm data from the James Clerk Maxwell telescope. But both estimates are within bounds set by current stellar data. It is therefore vital to continue accurate photometry of solar-like stars for at least another decade, to reveal secular and cyclic variations on multi-decadal time scales of direct interest to the Sun.

Herschelobservations of az∼ 2 stellar mass selected galaxy sample drawn from the GOODS NICMOS Survey

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.

Fomalhaut debris disk emission at 7 millimeters: constraints on the collisional models of planetesimals

We present new spatially resolved observations of the dust thermal emission at 7 mm from the Fomalhaut debris disk obtained with the Australia Telescope Compact Array. These observations provide the longest wavelength detection of the Fomalhaut debris disk to date. We combined the new data to literature sub-mm data to investigate the spectral index of the dust thermal emission in the sub-millimeter and constrained the $q$-slope of the power-law grain size distribution. We derived a value for $q = 3.48 \pm 0.14$ for grains with sizes around 1 mm. This is consistent with the classical prediction for a collisional cascade at the steady-state. The same value cannot be explained by more recent collisional models of planetesimals in which either the velocity distribution of the large bodies or their tensile strength is a strong function of the body size.

Radio continuum observations of Class I protostellar discs in Taurus: constraining the greybody tail at centimetre wavelengths★

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in the Taurus molecular cloud. These objects have previously been extensively studied in the sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical parametres.We use this new data to constrain the properties of the long-wavelength tail of the greybody spectrum, which is expected to be dominated by emission from large dust grains in the protostellar disk. We find spectra consistent with the opacity indices expected for such a population, with an average opacity index of beta = 0.26+/-0.22 indicating grain growth within the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave dust contribution. We find that disk masses derived from these flux densities under assumptions consistent with the literature are systematically higher than those calculated from sub-mm data, and meet the criteria for giant planet formation in a number of cases.

Study of the dense molecular gas surrounding the ‘extended green object’ G35.03+0.35

We present the results of a new study of the molecular gas associated with the "extended green object" (EGO) G35.03+0.35. This object, very likely a massive young stellar object, is embedded in a molecular cloud at the border of an HII region. The observations were performed with the Atacama Submillimeter Telescope Experiment (ASTE) in the 12CO and 13CO J=3-2, HCO+ J=4-3, and CS J=7-6 lines with an angular resolution about 22". From the 12CO J=3-2 line we discovered outflowing activity of the massive young stellar object. We obtained a total mass and kinetic energy for the outflows of 30 M_sun and 3000 M_sun (km/s)^2 (6 x 10^{46} ergs), respectively. We discovered a HCO+ and CS clump towards the EGO G35.03+0.35. The detection of these molecular species supports the presence of molecular outflows and a dense molecular envelope with temperatures and densities above 40 K and 6 x 10^{6} cm^{-3}, respectively. Using public near- and mid-IR, and sub-mm data we investigated the spectral energy distribution confirming that EGO G35.03+0.35 is a massive young stellar object at the earliest evolutionary stage (i.e. a class I young stellar object). By anlysing radio continuum archival data we found three radio sources towards the object, suggesting the presence of several young stellar objects in the region. Our radio continuum analysis is consistent with the presence of at least one ultracompact HII region and an hypercompact HII region or a constant-velocity ionized wind source.