Submillimeter Spectral Energy Distributions Research Articles

ArTéMiS imaging of the filamentary infrared dark clouds G1.75-0.08 and G11.36+0.80: Dust-based physical properties of the clouds and their clumps

Context. Filamentary infrared dark clouds (IRDCs) are a useful class of interstellar clouds for studying the cloud fragmentation mechanisms on different spatial scales. Determination of the physical properties of the substructures in IRDCs can also provide useful constraints on the initial conditions and early stages of star formation, including those of high-mass stars. Aims. We aim to determine the physical characteristics of two filamentary IRDCs, G1.75-0.08 and G11.36+0.80, and their clumps. We also attempt to understand how the IRDCs are fragmented into clumps. Methods. We imaged the target IRDCs at 350 and 450 µm using the bolometer called Architectures de bolomètres pour des Télescopes à grand champ de vue dans le domaine sub-Millimétrique au Sol (ArTéMiS). These data were used in conjunction with our previous 870 µm observations with the Large APEX BOlometer CAmera (LABOCA) and archival Spitzer and Berschel data. The LABOCA clump positions in G11.36+0.80 were also observed in the N2H+(1–0) transition with the Institut de Radioastronomie Millimétrique (IRAM) 30-metre telescope. Results. On the basis of their far-IR to submillimetre spectral energy distributions (SEDs), G1.75-0.08 was found to be composed of two cold (~14.5 K), massive (several ~103 M⊙) clumps that are projectively separated by ~3.7 pc. Both clumps are 70 µm dark, but they do not appear to be bounded by self-gravity. The G1.75-0.08 filament was found to be subcritical by a factor of ~14 with respect to its critical line mass, but the result is subject to uncertain gas velocity dispersion. The IRDC G11.36+0.80 was found to be moderately (by a factor of ~2) supercritical and composed of four clumps that are detected at all wavelengths observed with the ground-based bolometers. The SED-based dust temperatures of the clumps are ~13–15 K, and their masses are in the range ~232–633 M⊙. All the clumps are gravitationally bound and they appear to be in somewhat different stages of evolution on the basis of their luminosity-to-mass ratio. The projected, average separation of the clumps is ~1 pc. At least three clumps in our sample show hints of fragmentation into smaller objects in the ArTéMiS images. Conclusions. A configuration that is observed in G1.75-0.08, namely two clumps at the ends of the filament, could be the result of gravitational focussing acting along the cloud. The two clumps fulfil the mass-radius threshold for high-mass star formation, but if their single-dish-based high velocity dispersion is confirmed, their gravitational potential energy would be strongly overcome by the internal kinetic energy, and the clumps would have to be confined by external pressure to survive. Owing to the location of G1.75-0.08 near the Galactic centre (~270 pc), environmental effects such as a high level of turbulence, tidal forces, and shearing motions could affect the cloud dynamics. The observed clump separation in G11.36+0.80 can be understood in terms of a sausage instability, which conforms to the findings in some other IRDC filaments. The G11.36+0.80 clumps do not lie above the mass-radius threshold where high-mass star formation is expected to be possible, and hence lower-mass star formation seems more likely. The substructure observed in one of the clumps in G11.36+0.80 suggests that the IRDC has fragmented in a hierarchical fashion with a scale-dependent physical mechanism. This conforms to the filamentary paradigm for Galactic star formation.

A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics

ABSTRACT We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori, and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our Bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbour high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of ∼10 for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003−2238 and Mrk 1014 show evidence of dual AGNs in their SED fits, suggesting a minimum dual AGN fraction in the sample of 5 per cent.

Dense cores in the Seahorse infrared dark cloud: physical properties from modified blackbody fits to the far-infrared–submillimetre spectral energy distributions

Context. Infrared dark clouds (IRDCs) can be the birth sites of high-mass stars, and hence determining the physical properties of dense cores in IRDCs is useful to constrain the initial conditions and theoretical models of high-mass star formation. Aims. We aim to determine the physical properties of dense cores in the filamentary Seahorse IRDC G304.74+01.32. Methods. We used data from the Wide-field Infrared Survey Explorer (WISE), Infrared Astronomical Satellite (IRAS), and Herschel in conjuction with our previous 350 and 870 μm observations with the Submillimetre APEX Bolometer Camera (SABOCA) and Large APEX BOlometer CAmera, and constructed the far-IR to submillimetre spectral energy distributions (SEDs) of the cores. The SEDs were fitted using single or two-temperature modified blackbody emission curves to derive the dust temperatures, masses, and luminosities of the cores. Results. For the 12 analysed cores, which include two IR dark cores (no WISE counterpart), nine IR bright cores, and one H II region, the mean dust temperature of the cold (warm) component, the mass, luminosity, H2 number density, and surface density were derived to be 13.3 ± 1.4 K (47.0 ± 5.0 K), 113 ± 29 M⊙, 192 ± 94 L⊙, (4.3 ± 1.2) × 105 cm−3, and 0.77 ± 0.19 g cm−3, respectively. The H II region IRAS 13039-6108a was found to be the most luminous source in our sample ((1.1 ± 0.4) × 103 L⊙). All the cores were found to be gravitationally bound (i.e. the virial parameter αvir < 2). Two out of the nine analysed IR bright cores (22%) were found to follow an accretion luminosity track under the assumptions that the mass accretion rate is 10−5 M⊙ yr−1, the stellar mass is 10% of the parent core mass, and the radius of the central star is 5 R⊙. Most of the remaing ten cores were found to lie within 1 dex below this accretion luminosity track. Seven out of 12 of the analysed cores (58%) were found to lie above the mass-radius thresholds of high-mass star formation proposed in the literature. The surface densities of Σ > 0.4 g cm−3 derived for these seven cores also exceed the corresponding threshold for high-mass star formation. Five of the analysed cores (42%) show evidence of fragmentation into two components in the SABOCA 350 μm image. Conclusions. In addition to the H II region source IRAS 13039-6108a, some of the other cores in Seahorse also appear to be capable of giving birth to high-mass stars. The 22 μm dark core SMM 9 is likely to be the youngest source in our sample that has the potential to form a high-mass star (96 ± 23 M⊙ within a radius of ~0.1 pc). The dense core population in the Seahorse IRDC has comparable average properties to the cores in the well-studied Snake IRDC G11.11-0.12 (e.g. Tdust and L agree within a factor of ~1.8); furthermore, the Seahorse, which lies ~60 pc above the Galactic plane, appears to be a smaller (e.g. three times shorter in projection, ~100 times less massive) version of the Snake. The Seahorse core fragmentation mechanisms appear to be heterogenous, including cases of both thermal and non-thermal Jeans instability. High-resolution follow-up studies are required to address the fragmented cores’ genuine potential of forming high-mass stars.

Starbursts and dusty tori in distant 3CR radio galaxies

We present a study of the complete ultraviolet to submillimetre spectral energy distributions (SEDs) of twelve 3CR radio galaxy hosts in the redshift range $1.0 < z < 2.5$, which were all detected in the far-infrared by the Herschel Space Observatory. The study employs the new spectro-chemical evolutionary code P\'EGASE.3, in combination with recently published clumpy AGN torus models. We uncover the properties of the massive host galaxy stellar populations, the AGN torus luminosities, and the properties of the recent starbursts, which had earlier been inferred in these objects from their infrared SEDs. The P\'EGASE.3 fitting yields very luminous (up to 10$^{13}$L$_{\odot}$) young stellar populations with ages of several hundred million years in hosts with masses exceeding 10$^{11}$M$_{\odot}$. Dust masses are seen to increase with redshift, and a surprising correlation -- or better upper envelope behaviour -- is found between the AGN torus luminosity and the starburst luminosity, as revealed by their associated dust components. The latter consistently exceeds the former by a constant factor, over a range of one order of magnitude in both quantities.

SPECTRAL ENERGY DISTRIBUTIONS OF LOCAL LUMINOUS AND ULTRALUMINOUS INFRARED GALAXIES

Luminous and ultraluminous infrared galaxies ((U)LIRGs) are the most extreme star forming galaxies in the universe. The local (U)LIRGs provide a unique opportunity to study their multi-wavelength properties in detail for comparison to their more numerous counterparts at high redshifts. We present common large aperture photometry at radio through X-ray wavelengths, and spectral energy distributions (SEDs) for a sample of 53 nearby LIRGs and 11 ULIRGs spanning log (LIR/Lsun) = 11.14-12.57 from the flux-limited Great Observatories All-sky LIRG Survey (GOALS). The SEDs for all objects are similar in that they show a broad, thermal stellar peak and a dominant FIR thermal dust peak, where nuLnu(60um) / nuLnu(V) increases from ~2-30 with increasing LIR. When normalized at IRAS-60um, the largest range in the luminosity ratio, R(lambda)=log[nuLnu(lambda)/nuLnu(60um)] observed over the full sample is seen in the Hard X-rays (HX=2-10 keV). A small range is found in the Radio (1.4GHz), where the mean ratio is largest. Total infrared luminosities, LIR(8-1000um), dust temperatures, and dust masses were computed from fitting thermal dust emission modified blackbodies to the mid-infrared (MIR) through submillimeter SEDs. The new results reflect an overall ~0.02 dex lower luminosity than the original IRAS values. Total stellar masses were computed by fitting stellar population synthesis models to the observed near-infrared (NIR) through ultraviolet (UV) SEDs. Mean stellar masses are found to be log(M/Msun) = 10.79+/-0.40. Star formation rates have been determined from the infrared (SFR_IR~45Msun/yr) and from the monochromatic UV luminosities (SFR_UV~1.3Msun/yr), respectively. Multiwavelength AGN indicators have be used to select putative AGN: about 60% of the ULIRGs would have been classified as an AGN by at least one of the selection criteria.

STAR AND DUST FORMATION ACTIVITIES IN AzTEC-3, A STARBURST GALAXY ATz= 5.3

Analyses of high-redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observational constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. An important observational constraint neglected in the analysis is the mass of dust giving rise to the IR emission. In this paper we add this constraint to the analysis of AzTEC-3. Adopting an upper limit to the mass of stars and a bolometric luminosity for this object, we construct stellar and chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminosity before the associated stellar mass exceeds the observational limit. We find that the model with a Top Heavy IMF provided the most plausible scenario consistent with the observational constraints. In this scenario the dust formed over a period of ~200 Myr, with a SFR of ~500 Msun/yr. These values for the age and SFR in AzTEC-3 are significantly higher and lower, respectively, from those derived without the dust mass constraint. However, this scenario is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of our models is that all scenarios require most of the radiating dust mass to have been accreted in molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and of the use of dust evolution models in constraining the age and the star formation activity and history in galaxies.

Adaptive Optics Imaging of the AU Microscopii Circumstellar Disk: Evidence for Dynamical Evolution

We present an H-band image of the light scattered from circumstellar dust around the nearby (10 pc) young M star AU Microscopii (AU Mic, GJ 803, HD 197481), obtained with the Keck adaptive optics system. We resolve the disk both vertically and radially, tracing it over 17-60 AU from the star. Our high angular resolution (0.04 or 0.4 AU per resolution element) observations thus offer the opportunity to probe the morphology of the AU Mic debris disk on solar system scales. We identify substructures (dust clumps and gaps) in the disk that may point to the existence of orbiting planets. No planets are seen in our H-band image down to a limiting mass of 1MJ at >20 AU, although the existence of smaller planets cannot be excluded from the current data. Modeling of the disk surface brightness distribution at H band and R band, and of the optical to submillimeter spectral energy distribution, allows us to constrain the disk geometry and the dust grain properties. We confirm the nearly edge-on orientation of the disk inferred from previous observations and deduce that the disk may be clear inward of 1-10 AU. We find evidence for a lack of small grains at <50 AU, either as a result of grain growth or because of destruction by Poynting-Robertson and/or corpuscular drag. A change in the power-law index of the surface brightness profile is observed near 33 AU, similar to a feature known in the profile of the β Pic circumstellar debris disk. By comparing the timescales for interparticle collisions and Poynting-Robertson drag between the two systems, we argue that the breaks are linked to one or both of these processes.

Submillimeter Observations of the Low-Metallicity Galaxy NGC 4214

Results of submillimeter (450 and 850 μm) observations of a nearby dwarf irregular galaxy NGC 4214 with SCUBA on JCMT are presented. We aimed at examining the far-infrared–to–submillimeter spectral energy distribution (SED) and properties of dust thermal emission in a low-metallicity environment by choosing NGC 4214, in which the gas metallicity (log O/H + 12) is 8.34. We found that the SED is quite similar to those of the IRAS Bright Galaxies Sample (IBGS), which are local bright star-forming galaxies with metallicities comparable to the solar abundance. The dust temperature and emissivity index for NGC 4214 obtained by a fitting to the single temperature graybody model are Td = 35 ± 0.8 K and β = 1.4 ± 0.1, respectively, which are typical values for IBGS galaxies. Compiling the previous studies on similar nearby dwarf irregular galaxies, we found that NGC 1569 shows similar results to those of NGC 4214, while NGC 4449 and IC 10 SE show different SEDs and low emissivity indices. There seems to be a variety of SEDs among metal-poor dwarf irregular galaxies. We examined the dependence on the intensity of interstellar radiation field, as well as a two-temperature model, but the origin of the difference is not clear. Some mechanism(s) other than metallicity and the interstellar radiation field must be responsible for controlling dust emission properties.

Cirrus models for local and high-z SCUBA galaxies

We present a model for the ultraviolet to submillimetre emission from stars embedded in the general interstellar dust in galaxies (the ‘infrared cirrus’ component). Such emission is characterized by relatively low optical depths of dust and by cool (<30 K) dust temperatures. The model incorporates the stellar population synthesis model of Bruzual & Charlot and the dust model of Siebenmorgen & Krügel which includes the effects of small grains/polycyclic aromatic hydrocarbons. We apply the model to fit the optical to submillimetre spectral energy distributions (SEDs) of nearby galaxies, which are dominated by cirrus emission, and we find that our simple model is quite adequate to explain the observed SEDs. We also, more controversially, apply this cirrus model to the SEDs of high redshift sources detected in blank-field submillimetre surveys with the Submillimetre Common User Bolometer Array (SCUBA). Surprisingly, an excellent fit is found for many of these sources, with typical values for the optical depth AV and the surface brightness of the stellar radiation field ψ being only a factor of 2–3 higher than for nearby galaxies. This increase is not unreasonable given the expected evolution of dust optical depth in currently favoured star-formation history models. We conclude that the tendency to interpret the high-z SCUBA galaxies as very highly obscured starbursts may be premature and that these galaxies may be more closely linked to optically selected high redshift galaxies than previously assumed.

Dust Temperatures in theInfrared Space ObservatoryAtlas of Bright Spiral Galaxies

We examine far-infrared and submillimeter spectral energy distributions for galaxies in the Infrared Space Observatory Atlas of Bright Spiral Galaxies. For the 71 galaxies where we had complete 60-180 micron data, we fit blackbodies with lambda^-1 emissivities and average temperatures of 31 K or lambda^-2 emissivities and average temperatures of 22 K. Except for high temperatures determined in some early-type galaxies, the temperatures show no dependence on any galaxy characteristic. For the 60-850 micron range in eight galaxies, we fit blackbodies with lambda^-1, lambda-2, and lambda^-beta (with beta variable) emissivities to the data. The best results were with the lambda^-beta emissivities, where the temperatures were ~30 K and the emissivity coefficient beta ranged from 0.9 to 1.9. These results produced gas to dust ratios that ranged from 150 to 580, which were consistent with the ratio for the Milky Way and which exhibited relatively little dispersion compared to fits with fixed emissivities.

Quiescent Giant Molecular Cloud Cores in the Galactic Center

We have used the Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO) to map the far-infrared continuum emission (45-175 micrometer) toward several massive Giant Molecular Cloud (GMC) cores located near the Galactic center. The observed far-infrared and submillimeter spectral energy distributions imply low temperatures (approx. 15 - 22 K) for the bulk of the dust in all the sources, consistent with external heating by the diffuse ISRF and suggest that these GMCs do not harbor high- mass star-formation sites, in spite of their large molecular mass. Observations of FIR atomic fine structure lines of C(sub II) and O(sub I) indicate an ISRF enhancement of approx. 10(exp 3) in the region. Through continuum radiative transfer modeling we show that this radiation field strength is in agreement with the observed FIR and submillimeter spectral energy distributions, assuming primarily external heating of the dust with only limited internal luminosity (approx. 2 x 10(exp 5) solar luminosity). Spectroscopic observations of millimeter-wave transitions of H2CO, CS, and C-34S carried out with the Caltech Submillimeter Observatory (CSO) and the Institut de Radio Astronomie Millimetrique (IRAM) 30-meter telescope indicate a gas temperature of approx. 80 K, significantly higher than the dust temperatures, and density of approx. 1 x 10(exp 5)/cc in GCM0.25 + 0.01, the brightest submillimeter source in the region. We suggest that shocks caused by cloud collisions in the turbulent interstellar medium in the Galactic center region are responsible for heating the molecular gas. This conclusion is supported by the presence of wide-spread emission from molecules such as SiO, SO, and CH3OH, which are considered good shock tracers. We also suggest that the GMCs studied here are representative of the typical, pre-starforming cloud population in the Galactic center.

A molecular line study of NGC 1333/IRAS 4.

Molecular line surveys and fully sampled spectral line maps at 1.3 and 0.87 mm are used to examine the physical and chemical characteristics of the extreme Class I sources IRAS 4A and 4B in the L1450/NGC 1333 molecular cloud complex. A very well collimated, jetlike molecular outflow emanates from IRAS 4A, with a dynamical age of a few thousand years. Symmetric, clumpy structure along the outflow lobes suggests that there is considerable variability in the mass-loss rate or wind velocity even at this young age. Molecular emission lines toward IRAS 4A and 4B are observed to be weak in the velocity range corresponding to quiescent material surrounding the young stellar objects (YSOs). Depletion factors of 10-20 are observed for all molecules, including CO, even for even for very conservative mass estimates from the measured millimeter and submillimeter dust continuum. However, abundances scaled with respect to CO are similar to other dark molecular cloud cores. Such depletions could be mimicked by high dust optical depths or increased grain emissivities at the observing frequencies of 230 and 345 GHz, but the millimeter and submillimeter spectral energy distributions suggest that this is unlikely over the single-dish size scales of 5000-10,000 AU. Dense, outflowing gas is found to be kinematically, but not spatially, distinct from the quiescent material on these size scales. If CO is used as a chemical standard for the high-velocity gas, we find substantial enhancements in the abundances of several molecules in outflowing material, most notably CS, SiO, and CH3OH. The SiO emission is kinematically well displaced from the bulk cloud velocity and likely arises from directly shocked material. As is the case for CO, however, the outflow features from more volatile species are centered near the cloud velocity and are often characterized by quite low rotational temperatures. We suggest that grain-grain collisions induced by velocity shear zones surrounding the outflow axes transiently desorb the grain mantles, resulting in large abundance enhancements of selected species. Similar results have recently been obtained in several other low-mass YSOs, where the outflowing gas is often both kinematically and spatially distinct, and are illustrative of the ability of accretion and outflow processes to simultaneously modify the composition of the gas and dust surrounding young stars.