Sub-mm Observations Research Articles

Submillimetre flux as a probe of molecular ISM mass in high-z galaxies

Abstract Recent long-wavelength observations on the thermal dust continuum suggest that the Rayleigh–Jeans tail can be used as a time-efficient quantitative probe of the dust and interstellar medium (ISM) mass in high-z galaxies. We use high-resolution cosmological simulations from the Feedback in Realistic Environment (FIRE) project to analyse the dust emission of M* ≳ 1010 M⊙ galaxies at z= 2–4. Our simulations (MassiveFIRE) explicitly include various forms of stellar feedback, and they produce the stellar masses and star formation rates of high-z galaxies in agreement with observations. Using radiative transfer modelling, we show that sub-millimetre (sub-mm) luminosity and molecular ISM mass are tightly correlated and that the overall normalization is in quantitative agreement with observations. Notably, sub-mm luminosity traces molecular ISM mass even during starburst episodes as dust mass and mass-weighted temperature evolve only moderately between z = 4 and z = 2, including during starbursts. Our finding supports the empirical approach of using broadband sub-mm flux as a proxy for molecular gas content in high-z galaxies. We thus expect single-band sub-mm observations with ALMA to dramatically increase the sample size of high-z galaxies with reliable ISM masses in the near future.

Luminous, pc-scale CO 6−5 emission in the obscured nucleus of NGC 1377

High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles.

Extreme jet ejections from the black hole X-ray binary V404 Cygni

We present simultaneous radio through sub-mm observations of the black hole X-ray binary (BHXB) V404 Cygni during the most active phase of its June 2015 outburst. Our $4$ hour long set of overlapping observations with the Very Large Array, the Sub-millimeter Array, and the James Clerk Maxwell Telescope (SCUBA-2), covers 8 different frequency bands (including the first detection of a BHXB jet at $666 \,{\rm GHz}/450\mu m$), providing an unprecedented multi-frequency view of the extraordinary flaring activity seen during this period of the outburst. In particular, we detect multiple rapidly evolving flares, which reach Jy-level fluxes across all of our frequency bands. With this rich data set we performed detailed MCMC modeling of the repeated flaring events. Our custom model adapts the van der Laan synchrotron bubble model to include twin bi-polar ejections, propagating away from the black hole at bulk relativistic velocities, along a jet axis that is inclined to the line of sight. The emission predicted by our model accounts for projection effects, relativistic beaming, and the geometric time delay between the approaching and receding ejecta in each ejection event. We find that a total of 8 bi-polar, discrete jet ejection events can reproduce the emission that we observe in all of our frequency bands remarkably well. With our best fit model, we provide detailed probes of jet speed, structure, energetics, and geometry. Our analysis demonstrates the paramount importance of the mm/sub-mm bands, which offer a unique, more detailed view of the jet than can be provided by radio frequencies alone.

A modular set of synthetic spectral energy distributions for young stellar objects

In this paper, I present a new set of synthetic spectral energy distributions (SEDs) for young stellar objects (YSOs) spanning a wide range of evolutionary stages, from the youngest deeply embedded protostars to pre-main-sequence stars with few or no disks. These models include significant improvements on the previous generation of published models: in particular, the new models cover a much wider and more uniform region of parameter space, do not include highly model-dependent parameters, and include a number of improvements that make them more suited to modeling far-infrared and sub-mm observations of forming stars. Rather than all being part of a single monolithic set of models, the new models are split up into sets of varying complexity. The aim of the new set of models is not to provide the most physically realistic models for young stars, but rather to provide deliberately simplified models for initial modeling, which allows a wide range of parameter space to be explored. I present the design of the model set, and show examples of fitting these models to real observations to show how the new grid design can help us better understand what can be determined from limited unresolved observations. The models, as well as a Python-based fitting tool are publicly available to the community.

Magnetic field geometry of the large globule CB 34

We report the results of optical polarimetric observations of a Bok globule CB34 to study magnetic field structure on large scales ($10^5-10^6$ AU), which is combined with archival sub-mm observations to characterize the magnetic field structure of CB34 on small scales ($10^4 - 10^5$ AU). The optical polarization measurements indicate that the magnetic field in the globule is constrained to a maximum radius of $10^5$AU around the core, out to densities not smaller than $10^4$cm$^{-3}$. Our study is mainly concentrated on two submillimeter cores C1 and C2 of CB34. The direction of magnetic field of core C2 is found to be nearly perpendicular to the CO outflow direction of the globule. The magnetic field of core C1 is almost aligned with the minor axis of the core which is typical for magnetically dominated star formation models. The mean value of offset between the minor axis of core C2 and the outflow direction is found to be $14^\circ$ which suggests that the direction of the outflow is almost aligned with the minor axis of core C2. The magnetic field strength in the plane-of-sky for cores C1 and C2 is estimated to be $\approx 34\mu$G and $\approx 70\mu$G.

SOFIA MID-INFRARED IMAGING AND CSO SUBMILLIMETER POLARIMETRY OBSERVATIONS OF G034.43+00.24 MM1*

ABSTRACT We present 11.1 to 37.1 μm imaging observations of the very dense molecular cloud core MM1 in G034.43+00.24 using FORCAST on SOFIA and submillimeter (submm) polarimetry using SHARP on the Caltech Submillimeter Observatory. We find that at the spatial resolution of SOFIA, the point-spread function (PSF) of MM1 is consistent with being a single source, as expected based on millimeter (mm) and submm observations. The spectral energy distributions (SEDs) of MM1 and MM2 have a warm component at the shorter wavelengths not seen in mm and submm SEDs. Examination of H(1.65 μm) stellar polarimetry from the Galactic Plane Infrared Polarization Survey shows that G034 is embedded in an external magnetic field aligned with the Galactic Plane. The SHARP polarimetry at 450 μm shows a magnetic field geometry in the vicinity of MM1 that does not line up with either the Galactic Plane or the mean field direction inferred from the CARMA interferometric polarization map of the central cloud core, but is perpendicular to the long filament in which G034 is embedded. The CARMA polarimetry does show evidence for grain alignment in the central region of the cloud core, and thus does trace the magnetic field geometry near the embedded Class 0 YSO.

A study of the wiggle morphology of HH 211 through numerical simulations

Recent high-resolution high-sensitivity observations of protostellar jets have shown many to possess deviations to their trajectories. HH 211 is one such example where sub-mm observations with the SMA have revealed a clear reflection-symmetric wiggle. The most likely explanation is that the HH 211 jet source could be moving as part of a protobinary system. Here we test this assumption by simulating HH 211 through 3D hydrodynamic jet propagation simulations using the PLUTO code with a molecular chemistry and cooling module, and initial conditions based on an analytical model derived from SMA observations. Our results show the reflection-symmetric wiggle can be recreated through the assumption of a jet source perturbed by binary motion at its base, and that a regular sinusoidal velocity variation in the jet beam can be close to matching the observed knot pattern. However, a more complex model with either additional heating from the protostar, or a shorter period velocity pulsation may be required to account for enhanced emission near the source, and weaker knot emission downstream. Position velocity diagrams along the pulsed jet beam show a complex structure with detectable signatures of knots and show caution must be exercised when interpreting radial velocity profiles through observations. Finally, we make predictions for future HH 211 observations with ALMA.

Warm ISM in the Sagittarius A Complex

We present large-scale submm observations towards the Sgr A Complex covering ~300 arcmin2. These data were obtained in the frame of the HEXGAL Program with the Herschel-HIFI satellite and are complemented with submm observations obtained with the NANTEN2/SMART telescope as part of the NANTEN2/SMART Central Nuclear Zone Survey. The observed species are CO(4-3) observed with the NANTEN2/SMART telescope, and [CI](1-0), [CI](2-1), [NII](1-0), and [CII](3/2-1/2) observed with the Herschel-HIFI satellite. The observations are presented in a 1 km/s spectral resolution and a spatial resolution ranging from 46" to 28". The spectral coverage of the three lower frequency lines is +-200 km/s, while in the two high frequency lines, the upper LSR velocity limit is +94 km/s and +145 km/s for the [NII] and [CII] lines, respectively. The spatial distribution of the emission in all lines is very widespread. The bulk of the CO emission is found towards Galactic latitudes below the Galactic plane, and all the known molecular clouds are identified. Both neutral atomic carbon lines have their brightest emission associated with the +50 km/s cloud. Their spatial distribution at this LSR velocity describes a crescent-shape structure, which is probably the result of interaction with the energetic event (one or several supernovae explosions) that gave origin to the non-thermal Sgr A-East source. The [CII] and [NII] emissions have most of their flux associated with the thermal Arched-Filaments and the H Region and bright spots in [CII] emission towards the Central Nuclear Disk (CND) are detected. Warm Gas at very high (|Vlsr| > 100 km/s) LSR velocities is also detected towards the line of sight to the Sgr A Complex, and it is most probably located outside the region, in the X1 orbits.

H-ATLAS/GAMA: the nature and characteristics of optically red galaxies detected at submillimetre wavelengths

We combine Herschel/SPIRE sub-millimeter (submm) observations with existing multi-wavelength data to investigate the characteristics of low redshift, optically red galaxies detected in submm bands. We select a sample of galaxies in the redshift range 0.01$\leq$z$\leq$0.2, having >5$\sigma$ detections in the SPIRE 250 micron submm waveband. Sources are then divided into two sub-samples of $red$ and $blue$ galaxies, based on their UV-optical colours. Galaxies in the $red$ sample account for $\approx$4.2 per cent of the total number of sources with stellar masses M$_{*}\gtrsim$10$^{10}$ Solar-mass. Following visual classification of the $red$ galaxies, we find that $\gtrsim$30 per cent of them are early-type galaxies and $\gtrsim$40 per cent are spirals. The colour of the $red$-spiral galaxies could be the result of their highly inclined orientation and/or a strong contribution of the old stellar population. It is found that irrespective of their morphological types, $red$ and $blue$ sources occupy environments with more or less similar densities (i.e., the $\Sigma_5$ parameter). From the analysis of the spectral energy distributions (SEDs) of galaxies in our samples based on MAGPHYS, we find that galaxies in the $red$ sample (of any morphological type) have dust masses similar to those in the $blue$ sample (i.e. normal spiral/star-forming systems). However, in comparison to the $red$-spirals and in particular $blue$ systems, $red$-ellipticals have lower mean dust-to-stellar mass ratios. Besides galaxies in the $red$-elliptical sample have much lower mean star-formation/specific-star-formation rates in contrast to their counterparts in the $blue$ sample. Our results support a scenario where dust in early-type systems is likely to be of an external origin.

Magnetic and gaseous spiral arms in M83

Isotropic and anisotropic wavelet transforms are used to decompose the images of the spiral galaxy M83 in various tracers to quantify structures in a range of scales from 0.2 to 10 kpc. We used radio polarization observations at {\lambda}6 cm and 13 cm obtained with the VLA, Effelsberg and ATCA telescopes and APEX sub-mm observations at 870 {\mu}m, which are first published here, together with maps of the emission of warm dust, ionized gas, molecular gas, and atomic gas. The spatial power spectra are similar for the tracers of dust, gas, and total magnetic field, while the spectra of the ordered magnetic field are significantly different. The wavelet cross-correlation between all material tracers and total magnetic field is high, while the structures of the ordered magnetic field are poorly correlated with those of other tracers. -- The magnetic field configuration in M83 contains pronounced magnetic arms. Some of them are displaced from the corresponding material arms, while others overlap with the material arms. The magnetic field vectors at {\lambda}6 cm are aligned with the outer material arms, while significant deviations occur in the inner arms and in the bar region, possibly due to non-axisymmetric gas flows. Outside the bar region, the typical pitch angles of the material and magnetic spiral arms are very close to each other at about 10{\deg}. The typical pitch angle of the magnetic field vectors is about 20{\deg} larger than that of the material spiral arms. One of the main magnetic arms in M83 is displaced from the gaseous arms, while the other main arm overlaps a gaseous arm. We propose that a regular spiral magnetic field generated by a mean-field dynamo is compressed in material arms and partly aligned with them. The interaction of galactic dynamo action with a transient spiral pattern is a promising mechanism for producing such complicated spiral patterns as in M83.

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($\nu_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $\nu_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($\nu$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

Massive envelopes and filaments in the NGC 3603 star forming region

The formation of massive stars and their arrival on the zero-age main-sequence occurs hidden behind dense clouds of gas and dust. In the giant Hii region NGC 3603, the radiation of a young cluster of OB stars has dispersed dust and gas in its vicinity. At a projected distance of 2:5 pc from the cluster, a bright mid-infrared (mid-IR) source (IRS 9A) had been identified as a massive young stellar object (MYSO), located on the side of a molecular clump (MM2) of gas facing the cluster. We investigated the physical conditions in MM2, based on APEX sub-mm observations using the SABOCA and SHFI instruments, and archival ATCA 3 mm continuum and CS spectral line data. We resolved MM2 into several compact cores, one of them closely associated with IRS 9A. These are likely infrared dark clouds as they do not show the typical hot-core emission lines and are mostly opaque against the mid-IR background. The compact cores have masses of up to several hundred times the solar mass and gas temperatures of about 50 K, without evidence of internal ionizing sources. We speculate that IRS 9A is younger than the cluster stars, but is in an evolutionary state after that of the compact cores.

NUMERICAL SIMULATIONS OF HH 211: A REFLECTION-SYMMETRIC BIPOLAR OUTFLOW

Recent high-resolution, high-sensitivity observations of protostellar jets have shown many to possess an underlying ‘wiggle’ structure. HH 211 is one such example where recent sub-mm observations revealed a clear reflection-symmetric wiggle. An explanation for this is that the HH211 jet source is moving as part of a protobinary system. Here we test this assumption by simulating HH211 through 3D hydrodynamic simulations using the pluto code with a molecular chemistry and cooling module, and initial conditions based on an analytical model derived from SMA observations. Molecular chemistry allows us to accurately plot synthetic molecular emission maps and position-velocity diagrams for direct comparison to observations, enabling us to test the observational assumptions and put constraints on the physical parameters of HH211. Our preliminary results show that the reflection-symmetric wiggle can be recreated through the assumption of a jet source being part of a binary system.

NGC 4370: a case study for testing our ability to infer dust distribution and mass in nearby galaxies

A fraction of the early-type galaxy population hosts a prominent dust lane. Methods to quantify the dust content of these systems based on optical imaging data usually yield dust masses which are an order of magnitude lower than dust masses derived from the observed FIR emission. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/submm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods to determine the dust content in galaxies and explore new ones. We aim to derive the dust mass in NGC 4370 from both optical and FIR data, and investigate the need to invoke a putative diffuse dust component. We create color and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses are also derived from SED fits to FIR/submm observations. Finally, inverse radiative transfer fitting is performed to investigate more complex dust geometries. The empirical methods applied to the optical data yield lower limits of 3.4e5 solar masses, an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses which are slightly lower, but fully consistent with the FIR-derived mass. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, using inverse radiative transfer fitting, we are able to constrain most of the parameters describing these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/submm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution for the interpretation of dust masses and optical depths based on optical data alone, using overly simplistic star-dust geometries and the neglect of scattering effects. [ABRIDGED]

Coordinated Hubble Space Telescope and Venus Express Observations of Venus’ upper cloud deck

Hubble Space Telescope Imaging Spectrograph (HST/STIS) UV observations of Venus’ upper cloud tops were obtained between 20N and 40S latitude on December 28, 2010; January 22, 2011 and January 27, 2011 in coordination with the Venus Express (VEx) mission. The high spectral (0.27nm) and spatial (40–60km/pixel) resolution HST/STIS data provide the first direct and simultaneous record of the latitude and local time distribution of Venus’ 70–80km SO and SO2 (SOx) gas density on Venus’ morning quadrant. These data were obtained simultaneously with (a) VEx/SOIR occultation and/or ground-based James Clerk Maxwell Telescope sub-mm observations that record respectively, Venus’ near-terminator SO2 and dayside SOx vertical profiles between ∼75 and 100km; and (b) 0.36μm VEx/VMC images of Venus’ cloud-tops. Updating the (Marcq, E. et al. [2011]. Icarus 211, 58–69) radiative transfer model SO2 gas column densities of ∼2–10μm-atm and ∼0.4–1.8μm-atm are retrieved from the December 2010 and January 2011 HST observations, respectively on Venus’ dayside (i.e., at solar zenith angles (SZA)<60°); SO gas column densities of 0.1–0.11μm-atm, 0.03–0.31μm-atm and 0.01–0.13μm-atm are also retrieved from the respective December 28, 2010, January 22, 2011 and January 27, 2011 HST observations. A decline in the observed low-latitude 0.24 and 0.36μm cloud top brightness paralleled the declining SOx gas densities. On December 28, 2010 SO2 VMR values ∼280–290ppb are retrieved between 74 and 81km from the HST and SOIR data obtained near Venus’ morning terminator (at SZAs equal to 70° and 90°, respectively); these values are 10× higher than the HST-retrieved January 2011 near terminator values. Thus, the cloud top SO2 gas abundance declined at all local times between the three HST observing dates. On all dates the average dayside SO2/SO ratio inferred from HST between 70 and 80km is higher than that inferred from the sub-mm the JCMT data above 84km confirming that SOx photolysis is more efficient at higher altitudes. The direct correlation of the SOx gases provides the first clear evidence that SOx photolysis is not the only source for Venus’ 70–80km sulfur reservoir. The cloud top SO2 gas density is dependent in part on the vertical transport of the gas from the lower atmosphere; and the 0.24μm cloud top brightness levels are linked to the density of the sub-micron haze. Thus, the new results may suggest a correlation between Venus’ cloud-top sub-micron haze density and the vertical transport rate. These new results must be considered in models designed to simulate and explore the relationship between Venus’ sulfur chemistry cycle, H2SO4 cloud formation rate and climate evolution. Additionally, we present the first photochemical model that uniquely tracks the transition of the SO2 atmosphere from steady to non-steady state with increasing SZA, as function of altitude within Venus’ mesosphere, showing the photochemical and dynamical basis for the factor of ∼2 enhancements in the SOx gas densities observed by HST near the terminator above that observed at smaller SZA. These results must also be considered when modeling the long-term evolution of Venus’ atmospheric chemistry and dynamics.

ALMA RESOLVES THE PROPERTIES OF STAR-FORMING REGIONS IN A DENSE GAS DISK AT z ∼ 3

We exploit long-baseline ALMA sub-mm observations of the lensed star-forming galaxy SDP 81 at z=3.042 to investigate the properties of inter-stellar medium on scales of 50-100pc. The kinematics of the CO gas within this system are well described by a rotationally-supported disk with an inclination-corrected rotation speed, v=320+/-20km/s and a dynamical mass of M=(3.5+/-1.0)x10^10Mo within a radius of 1.5 kpc. The disk is gas rich and unstable, with a Toomre parameter, Q=0.30+/-0.10 and so should collapse in to star-forming regions with Jeans length L_J~130pc. We identify five star-forming regions within the ISM on these scales and show that their scaling relations between luminosity, line-widths and sizes are significantly offset from those typical of molecular clouds in local Galaxies (Larson's relations). These offsets are likely to be caused by the high external hydrostatic pressure for the interstellar medium (ISM), P/kB=(40+/-20)x10^7K/cm3, which is ~10,000x higher than the typical ISM pressure in the Milky Way. The physical conditions of the star-forming ISM and giant molecular clouds appears to be similar to the those found in the densest environments in the local Universe, such as those in the Galactic center.

SUB-mm JET PROPERTIES OF THE X-RAY BINARY SWIFT J1745–26

We present the results of our observations of the early stages of the 2012--2013 outburst of the transient black hole X-ray binary (BHXRB), Swift J1745$-$26, with the VLA, SMA, and JCMT (SCUBA--2). Our data mark the first multiple-band mm & sub-mm observations of a BHXRB. During our observations the system was in the hard accretion state producing a steady, compact jet. The unique combination of radio and mm/sub-mm data allows us to directly measure the spectral indices in and between the radio and mm/sub-mm regimes, including the first mm/sub-mm spectral index measured for a BHXRB. Spectral fitting revealed that both the mm (230 GHz) and sub-mm (350 GHz) measurements are consistent with extrapolations of an inverted power-law from contemporaneous radio data (1--30 GHz). This indicates that, as standard jet models predict, a power-law extending up to mm/sub-mm frequencies can adequately describe the spectrum, and suggests that the mechanism driving spectral inversion could be responsible for the high mm/sub-mm fluxes (compared to radio fluxes) observed in outbursting BHXRBs. While this power-law is also consistent with contemporaneous optical data, the optical data could arise from either jet emission with a jet spectral break frequency of $\nu_{{\rm break}}\gtrsim1\times10^{14}\,{\rm Hz}$ or the combination of jet emission with a lower jet spectral break frequency of $\nu_{{\rm break}}\gtrsim2\times10^{11}\,{\rm Hz}$ and accretion disc emission. Our analysis solidifies the importance of the mm/sub-mm regime in bridging the crucial gap between radio and IR frequencies in the jet spectrum, and justifies the need to explore this regime further.

From flux to dust mass: Does the grain-temperature distribution matter for estimates of cold dust masses in supernova remnants?

The amount of dust estimated from infrared to sub-millimetre (submm) observations strongly depends on assumptions of different grain sizes, compositions and optical properties. Here we use a simple model of thermal emission from cold silicate/carbon dust at a range of dust grain temperatures and fit the spectral energy distribution (SED) of the Crab Nebula as a test. This can lower the derived dust mass for the Crab by ~50% and 30-40% for astronomical silicates and amorphous carbon grains compared to recently published values (0.25M_sun -> 0.12M_sun and 0.12M_sun -> 0.072M_sun, respectively), but the implied dust mass can also increase by as much as almost a factor of six (0.25M_sun -> 1.14M_sun and 0.12M_sun -> 0.71M_sun) depending on assumptions regarding the sizes/temperatures of the coldest grains. The latter values are clearly unrealistic due to the expected metal budget, though. Furthermore, we show by a simple numerical experiment that if a cold-dust component does have a grain-temperature distribution, it is almost unavoidable that a two-temperature fit will yield an incorrect dust mass estimate. But we conclude that grain temperatures is not a greater uncertainty than the often poorly constrained emissivities (i.e., material properties) of cosmic dust, although there is clearly a need for improved dust emission models. The greatest complication associated with deriving dust masses still arises in the uncertainty in the dust composition.

Submillimetre observations of WISE/radio-selected AGN and their environments

We present JCMT SCUBA-2 850microns submillimetre (submm) observations of 30 mid-infrared (mid-IR) luminous AGN, detected jointly by the WISE all-sky IR survey and the NVSS/FIRST radio survey. These rare sources are selected by their extremely red mid-infrared spectral energy distributions (SEDs) and compact radio counterparts. Further investigations show that they are highly obscured, have abundant warm AGN-heated dust and are thought to be experiencing intense AGN feedback. These galaxies appear to be consistent with an AGN-dominated galaxy, and could be a transient phase of merging galaxies. When comparing the number of submm galaxies (SMGs) detected serendipitously in the surrounding 1.5-arcmin to those in blank-field submm surveys, there is a very significant overdensity, of order 5, but no sign of radial clustering centred at our primary objects. The WISE/radio-selected AGN thus reside in 10-Mpc-scale overdense environments, that could be forming in pre-viralised clusters of galaxies. WISE/radio-selected AGNs appear to be the strongest signposts of high-density regions of active, luminous and dusty galaxies. SCUBA-2 850microns observations indicate that their submm fluxes are low compared to many popular AGN SED templates, hence the WISE/radio-selected AGNs have either less cold and/or more warm dust emission than normally assumed for typical AGN. Most of the targets are not detected, only four targets are detected at SCUBA-2 850microns, and have total IR luminosities >= 10^13 L_solar, if their redshifts are consistent with the subset of the 10 SCUBA-2 undetected targets with known redshifts, z ~ 0.44 - 2.86.

Submillimetre observations of WISE-selected high-redshift, luminous, dusty galaxies

We present SCUBA-2 850um submillimetre (submm) observations of the fields of 10 dusty, luminous galaxies at z ~ 1.7 - 4.6, detected at 12um and/or 22um by the WISE all-sky survey, but faint or undetected at 3.4um and 4.6um; dubbed hot, dust-obscured galaxies (Hot DOGs). The six detected targets all have total infrared luminosities greater than 10^13 L_sun, with one greater than 10^14 L_sun. Their spectral energy distributions (SEDs) are very blue from mid-infrared to submm wavelengths and not well fitted by standard AGN SED templates, without adding extra dust extinction to fit the WISE 3.4um and 4.6um data. The SCUBA-2 850um observations confirm that the Hot DOGs have less cold and/or more warm dust emission than standard AGN templates, and limit an underlying extended spiral or ULIRG-type galaxy to contribute less than about 2% or 55% of the typical total Hot DOG IR luminosity, respectively. The two most distant and luminous targets have similar observed submm to mid-infrared ratios to the rest, and thus appear to have even hotter SEDs. The number of serendipitous submm galaxies (SMGs) detected in the 1.5-arcmin-radius SCUBA-2 850um maps indicates there is a significant over-density of serendipitous sources around Hot DOGs. These submm observations confirm that the WISE-selected ultra-luminous galaxies have very blue mid-infrared to submm SEDs, suggesting that they contain very powerful AGN, and are apparently located in unusual arcmin-scale overdensities of very luminous dusty galaxies.