Anomalous Microwave Emission: Theory, Modeling, and ObservationsView this Special Issue
Research Article | Open Access
Investigating the Source of Planck-Detected AME: High-Resolution Observations at 15 GHz
The Planck 28.5 GHz maps were searched for potential Anomalous Microwave Emission (AME) regions on the scale of ~3° or smaller, and several new regions of interest were selected. Ancillary data at both lower and higher frequencies were used to construct spectral energy distributions (SEDs), which seem to confirm an excess consistent with spinning dust models. Here we present higher resolution observations of two of these new regions with the Arcminute Microkelvin Imager Small Array (AMI SA) between 14 and 18 GHz to test for the presence of a compact (~10 arcmin or smaller) component. For AME-G107.1+5.2, dominated by the Hii region S140, we find evidence for the characteristic rising spectrum associated with either the spinning dust mechanism for AME or an ultra- /hypercompact Hii region across the AMI frequency band; however, for AME-G173.6+208 we find no evidence for AME on scales of ~2–10 arcmin.
The Planck satellite [1, 2] observes the sky in nine frequency bands, covering a range from 30 to 857 GHz. Its wide frequency range potentially allows the detection of AME since the high-frequency data above 100 GHz can be used to constrain the thermal emission, while the lower frequency data is close to the theoretical peak of the spinning dust emission. When combined with ancillary data at lower frequencies, the spectra of AME regions can be accurately determined and used to probe their properties on large (1 degree) scales. The Planck maps were used to detect several new potential AME regions  by subtracting a spatial model of known emission mechanisms (synchrotron, free-free, and thermal dust) extrapolated from observational or theoretical predictions. Two of these regions, AME-G173.6+2.8 and AME-G107.1+5.2, were selected, and ancillary data were used at both higher and lower frequencies to construct SEDs, which contain suggestions of AME consistent with spinning dust emission.
The Arcminute Microkelvin Imager Small Array (AMI SA) is a radio interferometer situated near Cambridge, UK. Primarily an SZ survey instrument, the AMI-SA, is specifically designed to have high sensitivity to low-surface-brightness emission on scales of 2–10 arcmin. It operates between 14 and 18 GHz, close in frequency to the theoretical peak of the spinning dust emission, the position of which varies between 10–50 GHz depending on grain size, composition, and ambient conditions. The AMI SA has previously been used both to identify and to characterize spinning dust regions in multiple galactic (e.g. [4, 5]) and extragalactic  sources. The higher angular resolution and lower frequency coverage of the AMI SA make it a highly complementary instrument to the Planck satellite for studies of AME. The synthesised beam of the SA, which is an effective measure of the resolution, is ≃2 arcmin at FWHM, while the Planck maps and ancillary data were smoothed to a common resolution of ≃1° for AME detection.
2. Observations and Data Reduction
The AMI SA is situated at the Mullard Radio Astronomy Observatory, Cambridge, UK . It consists of ten 3.7 m diameter dishes with a baseline range of ≃5–20 m and observes in the band 12–18 GHz with eight 0.75-GHz bandwidth channels. In practice, the lowest two frequency channels are unused due to a low response in this frequency range and interference from geostationary satellites.
The two sources listed in Table 1 were observed from 2011 Jan 21–Feb 09, for ≃5 hours each.
Data reduction was performed using the local software tool REDUCE, which flags interference, shadowing, and hardware errors, and applies phase and amplitude calibrations. The correlator data is then Fourier transformed to synthesize the frequency channels, before output to disk in uv FITS format. Flux calibration was performed using short observations of 3C48, 3C286, or 3C147 near the beginning and the end of each run. The assumed flux densities for 3C286 were converted from very large array total-intensity measurements provided by Perley (private communication) and are consistent with the Rudy et al.  model of Mars transferred on to absolute scale, using results from the Wilkinson Microwave Anisotropy Probe. The assumed flux densities for 3C48 and 3C147 are based on long-term monitoring with the AMI SA using 3C286 for flux calibration (see Table 2). A correction for changing airmass is also applied using a noise-injection system, the “rain gauge.”
Bright, nearby point sources selected from the Very Long Baseline Array Calibrator Survey  were observed during each observation at hourly intervals for phase calibration purposes (see Table 1 for phase calibrators used for the AMI SA observations). The reduced visibility data were imaged using AIPS (http://aips.nrao.edu/), from the individual channel datasets (for channels 3 to 8 inclusive), as well as from the combined channels at a central frequency of 15.75 GHz. Gaussians were fitted to the sources detected at 5σ in the maps using the AIPS task JMFIT. Errors on AMI SA flux density values were estimated by adding in quadrature the error output from JMFIT, which folds in an estimate of the r.m.s. map noise and the error associated with the Gaussian fit, and the error on flux calibration (including rain-gauge correction) of ≃5 percent of the integrated flux density.
2.1. Matching Spatial Scales
Since the sources detected were extended to the AMI SA beam, flux loss corrections had to be applied. These were calculated by sampling Canadian Galactic Plane Survey (CGPS, ) 1.4 GHz total power maps with the uv coverage of the AMI observations, followed by mapping and fitting Gaussians to the sampled maps in the same manner as the AMI observations. Figure 2 shows the sampled visibilities compared to the AMI observed visibilities and the ratios derived from fitting elliptical Gaussians to the maps. As a consistency check, the flux loss corrections were also calculated by modeling the sources as Gaussians based on the size parameters derived from the AMI continuum maps. The channel 3 flux density was taken as the reference, and channels 4–8 and the CGPS 1.42 GHz flux density were corrected to the channel 3 scale using the fitted flux loss percentages.
Lower frequency flux densities at 0.408, 2.7, and 4.85 GHz are estimated from the lower resolution CGPS, Effelsberg, and GB6  maps by convolving the CGPS 1.42 GHz maps to the appropriate resolution, measuring a flux density on the convolved CGPS and lower resolution maps, and calculating a spectral index. The lower resolution flux density is then normalised to the flux density estimated from the CGPS map sampled with AMI channel 3 uv coverage using the calculated spectral index. These fluxes are, however, considered as upper limits since the amount of extra background flux density measured by the lower resolution maps is uncertain and is not accounted for in the flux loss correction.
AMI SA maps of both sources are shown in Figures 3 and 7. The maps displayed are not corrected for attenuation due to the primary beam; the flux densities reported have been so corrected, where spectral indices are quoted, the convention is used, is flux density, and is frequency. Errors quoted are .
This region is centred on the HII region S235 , part of the FVW 172.8+1.5 complex within Auriga. FVW 172.8+1.5 is an HII complex composed of Sharpless HII regions and OB associations in the Perseus arm with individual objects organized along a large (~7–4 deg) filamentary structure, the morphology of which has been compared to a bow tie . Several HII regions within this complex have been the subject of previous study (e.g. ). The north-eastern regions: S231, S232, S233, and S235, are known to be associated with a giant molecular cloud at an adopted distance of 1.8 kpc (; individual distances determined by spectrophotometry to each of the exciting stars within the HII regions range from 1.0 to 2.3 kpc) and have systematic CO radial velocities of −18.1 to −23.0 km.
It has recently been proposed that these four HII regions lie on the shell of a SNR . This premise, rather than that of an expanding HII region, is supported by the fast expansion velocity of the shell and high kinetic energy associated with it. In addition the filamentary structures joining individual HII regions were shown to have nonthermal spectral contributions, inconsistent with an HII bubble.
Figure 3 shows the AMI SA channel-averaged map and the corresponding CGPS 1.4 GHz map. The smoothing resolution of the Planck analysis is shown for comparison. It can be seen that there are many components within the smoothing radius, both diffuse and compact, which could potentially contribute to the excess emission seen by Planck.
A number of discrete sources are detected within the AMI SA primary beam towards this region. The field is dominated by the north-south string of HII regions: S235, S235A (BFS46), and BFS47. BFS46 (S235A) was previously part of an AME study of HII regions  where it was found to exhibit no excess, having a spectrum consistent with optically thin free-free emission with . This provides a useful check for any systematic offsets in individual channel calibration that may affect the spectrum of S235. Although it lies away from the centre of the pointing where phase errors are expected to have greater effect, the flux densities derived from the AME-G173.6+2.8 map for BFS46 are consistent with those from the previous study.
Here we also confirm that the spectrum of S235 is consistent with optically thin free-free emission. After flux loss correction using the CGPS sampled maps, a spectral index of is fitted to the AMI channel flux densities. Using a simulated elliptical Gaussian with the deconvolved source size and position angle from the AMI continuum map as a consistency check gives a spectral index of , and it is likely that the difference indicates that this source is not well modelled by a Gaussian. Figure 1(a) shows the CGPS sampled visibilities compared to the AMI visibilities and shows the correspondence derived from fitting to the maps. It can be seen that the correlation between the visibilities is very good, indicating a good morphological correlation between 1.42 and 16 GHz.
Although a spectrum with an index of is steeper than the canonical free-free index of , such indices are expected for bremsstrahlung from cooler (heavy) plasma . However, it is more likely in this instance that the steeper spectrum across the AMI band arises as a consequence of inexact flux loss corrections due to differences in the structure of the source between 1.4 GHz and the frequency range covered by AMI. Fitting to the CGPS 1.42 GHz catalogue flux for S235 and the AMI channel fluxes gives a spectral index of , consistent with the radio recombination line measurements of Lockman  who found kmfor this region, indicative of optically thin plasma. Figure 2(a) shows the SED indicating the corrections for flux loss; it can be clearly seen that AMI sees no AME. Flux densities from the Effelsberg 11 cm and GB6 4.85 GHz surveys should be considered upper limits as the flux loss corrections calculated from the CGPS data fitted at a resolution of 1 arcmin cannot account for any contribution from proximal low-surface-brightness emission which has contributed to the flux density of sources fitted at resolutions of 4.3 and 3 arcmin, respectively. Indeed, when considered in conjunction with lower frequency data, Figure 2(b), we find that the AMI data is well fitted by a free-free spectrum (; see Scaife 2012, this volume). Fixing the electron temperature at K, we find an emission measure of pc, assuming a source with deconvolved size as determined from the AMI data ( P.A. = ).
This region is centred on the HII region S140 . At a distance of 910 pc, S140 is a bright-rimmed cloud that forms the interface between an HII region and the molecular cloud L1204 . S140 contains an IR cluster of at least three sources . IRS1 is the progenitor of a large molecular outflow in the SE-NW direction with a reflection nebula associated with the blue-shifted SE lobe (; see also ). The activity of this outflow has caused significant disturbance in the surrounding area where a number of sweeping feather-like structures are visible in high angular resolution K-band images . Perpendicular to the outflow a disk has been detected at cm wavelengths, confirmed by K-band polarimetric imaging ([26, 27]).
The AMI observation of this region is shown in Figure 7. After flux loss correction, a spectral index of is fitted to the AMI channel fluxes. The spectral index calculated using a simulated Gaussian source with size derived from the AMI channel-averaged map for flux loss correction is ; however, it can be seen that the source is highly non-Gaussian which is likely to account for the difference. Figure 1 shows the CGPS sampled visibilities compared to the AMI visibilities and the flux density ratio derived from fitting to the maps. It can be seen that the correlation between the visibilities is very good, indicating a good morphological correlation between 1.42 and 15 GHz and supporting the rising spectrum derived by using the CGPS data to correct the AMI flux densities for flux loss.
Fitting to the sampled CGPS 1.42 GHz flux density and AMI channel 3 uv coverage gives a spectral index of , consistent with optically thin free-free emission. Fluxes at lower resolution were not estimated due to the more complex background of this source; the CGPS 408 MHz catalogue flux density adjusted by the same flux loss percentage as the 1.42 GHz value is included as a reference; however, the errors are extremely large. Figure 4(a) shows the SED; it seems likely that the AMI fluxes indicate the presence of AME or an ultra-/hypercompact HII region. Figure 7 shows the AMI SA channel-averaged map and the corresponding CGPS 1.4 GHz map.
Two scenarios are possible that may explain the rising spectrum across the AMI band: scenario 1 is the presence of a hypercompact HII (HCHII) region and scenario 2 is the presence of a spinning dust component. Scenario 1 is illustrated in Figure 4(b) where the lower frequency data is dominated by free-free emission from plasma with an emission measure of pc assuming an electron temperature of K. This would indicate a region of classical HII; however, the relatively high-emission measure implied by the turnover frequency lying between 408 MHz and 1.4 GHz requires that the source have a much smaller angular size (0.1 arcmin) than that fitted to the AMI map (P.A. = ) to be consistent with the magnitude of the flux density. In addition the best fitting emission measure for the high-frequency component is pc, which is extremely high.
Fixing the source size to that measured from the AMI data produces a poorer fit to the low-frequency data, with an emission measure of pc, see Figure 5. The rising spectrum across the AMI band requires a second free-free component with pc, consistent with an ultra- or hypercompact HII region. In this case we find a best fitting emission measure for the high-frequency component of pc with a size of ~0.001 pc consistent with an object on the ultra- /hypercompact border. In spite of the poorer fit to the low-frequency data we consider this case more likely than the first as it is consistent with the measured source extent and predicts emission measures which are more physically feasible.
Scenario 2 is illustrated in Figure 6(a) where the rising spectrum in the AMI data is accounted for by the introduction of a spinning dust model, here taken as the Draine and Lazarian Warm Ionized Medium (WIM) model . For the small angular size of the emitting region the column density of the spinning dust emitting region needs to be reasonably high at . A further possibility is that the 408 MHz datum defines an optically thin free-free component, whilst the 1.4 GHz point represents a contribution from a more compact region of HII, and the rising spectrum across the AMI band indicates a further component arising from even more compact HII or spinning dust. This possibility is shown in Figure 6(b); interestingly the optically thin free-free component shown has a similar emission measure to that fitted to the low-frequency data on larger scales by , although the uncertainty on the fitted emission measure due to the large error on the 408 MHz point makes the comparison weak.
With only the upturn in the SED visible over the AMI band it is not possible to distinguish between these two scenarios using the AMI data alone. We do, however, confirm that a rising spectrum is present at frequencies above 10 GHz. However, we also note that the magnitude of the emission required to explain the rising AMI spectrum is approximately an order of magnitude, too small to explain the excess observed with Planck at ≈1° resolution.
Detections of AME from Planck data are based on flux densities integrated within an aperture of 1-degree radius. Such measurements will contain contributions from a variety of discrete sources, in addition to the large-scale diffuse background. The data presented here allow us to examine the population of discrete sources that contribute to these detections and investigate the physics behind their emission over a frequency range sensitive to the proposed spinning dust mechanism, which is generally considered to provide the physical basis for AME.
The data presented here have shown that the emission from S235, the dominant source within the G173.6+2.8 region, has a spectrum consistent with optically thin free-free emission and shows no indication of a spinning dust component. For the source S140, dominant within the G107.1+5.2 region, the data presented here have shown that the spectrum of this source rises across the AMI band and is not consistent with an extrapolation of optically thin free-free from lower frequencies.
4.1. The Spectrum of G107.1+5.2
On large scales the spectrum of G107.1+5.2  is ambiguous in nature. Here we examine that ambiguity by comparing different physical scenarios through their spectra (data points from Dickinson, priv. comm.). In Figure 8(a) we show a scenario where the excess emission at microwave frequencies is provided by the presence of an ultracompact HII region, and in Figure 8(b) we show a scenario where the excess emission at microwave frequencies is provided by a spinning dust component. In both cases maximum likelihood (ML) fitted models are shown. These fits were found using an MCMC-based method which calculates both ML parameters with associated uncertainties, as well as the evidence values for the specific model using a simulated annealing approach .
We consider models of the same form as the original Planck analysis with the electron temperature of the free-free emitting gas fixed to K but leave the emission measure free to vary. Similarly we fix the temperature of the dust component to be K  as this is determined largely by the peak of the dust greybody; we retain the opacity index, , as a free parameter as well as a normalization coefficient. To these two initial components we add a second free-free component, parameterized by a second emission measure and by a size which is defined relative to the size of the first free-free emitting region. The electron temperature of this second component is also fixed at K.
In the second scenario we assume an identical first free-free component and similarly a dust greybody with a fixed dust temperature. To these two components we add a spinning dust component parameterized by a column density, , and include a potential CMB contribution in a similar fashion to . The spinning dust component that we use is the Draine and Lazarian WIM model . The model parameters we obtain differ slightly from those of  but are consistent within errors.
In the first scenario we find that the data are well fitted by a two-component free-free model where the second component has an emission measure of pc and a source size which is greater than a factor of 1200, smaller than the size of the first emitting region. Both constraints are consistent with the presence of an ultracompact HII region as the relative source size indicates an object of 0.1 pc at a distance of 910 pc. The lack of data in the frequency range 1.4–20 GHz makes these two parameters highly degenerate as the optical depth of such an object can reach unity anywhere over an order of magnitude in frequency. Indeed a rising spectrum over the AMI band, as observed, would indicate a turn over frequency at GHz and consequently an emission measure of pc, consistent with the fitted constraints.
In the second scenario we find parameter values similar to those determined by , with the exception of the CMB contribution where we find a maximum likelihood value of approximately half of that determined by . Maximum likelihood parameters for both models are listed in Tables 3 and 4. These models are shown to be fitted to the data in Figures 8(a) and 8(b). The logarithmic difference in the evidence between the two models, as determined over a prior volume, is in favour of the spinning dust model.
|†Radius of circular aperture. ‡Listed parameters represent the best fit values, and uncertainties are highly nonsymmetric. |
|†Radius of circular aperture. |
Assuming that the optically thin free-free component which dominates the Planck SED arises on scales which are not visible to AMI, we examine the possibility that the rising spectrum observed across the AMI band and the excess emission detected by Planck can both be explained by a single hypercompact HII region. If this is the case, we predict that the emission seen across the AMI band in a Planck-size aperture would be simply the sum of an extrapolation of the optically thin component to these frequencies plus the measured AMI flux densities. This is subject to a couple of strong assumptions: firstly, the emission from the HCHII region arises on scales which are compact to the AMI beam, and there will consequently be no flux loss relative to the Planck measurement; secondly, the HCHII is the dominant contribution to the AMI flux densities. The first of these assumptions is reliable as HCHII exists on scales which are highly compact to both the AMI and Planck beams. The second assumption is less reliable as it depends on the parameters of the HCHII model, which are not possible to constrain absolutely from the AMI data alone. Considering this to be the case it is immediately obvious that the magnitude of the AMI flux densities is insufficient to be compatible with the excess seen by Planck being contributed by a single HCHII region.
By the same token it is also not possible for the Planck excess to be explained by the presence of a spinning dust region with angular scales entirely recovered by AMI. Since the possibility of an extended region of HCHII has no precedent in the literature, it seems likely then that the excess measured by Planck towards G107.1+5.2 arises from a spinning dust region, the majority of which is extended on scales larger than those measured by AMI. This conclusion is supported by the strong evidence in favour of a spinning dust component found by the model comparison.
In the case of both G173.6+2.8 (S235) and G107.1+5.2 (S140) we conclude that the bulk of the excess emission seen by Planck must arise on scales larger than 10 arcmin; the magnitude of the emission seen by AMI is 2% of that expected at 15 GHz at ≈1° resolution. In the case of G173.6+2.8 we confirm that the dominant source within the Planck aperture, S235, has a spectrum on scales of 2–10 arcmin consistent with optically thin free-free emission. In the case of G107.1+5.2 we demonstrate that the dominant source, S140, has a rising spectrum across the AMI band. This spectrum is consistent with either spinning dust emission or the presence of UC/HCHII. With only the low-frequency data available on these scales it is not possible to distinguish between the two mechanisms; however, we conclude that the magnitude of the contribution causing the rising spectral index across the AMI band is not sufficient to explain the measured Planck excess, which must once again arise on scales larger than those probed by AMI. In this case, it is therefore possible that the rising spectrum over the AMI band is also a consequence of spinning dust emission.
The authors thank the staff of the Mullard Radio Astronomy Observatory for their invaluable assistance in the commissioning and operation of AMI, which is supported by Cambridge University and the STFC. They also thank Clive Dickinson for the data points for the large-scale spectrum of G107.1+5.2 and the referee for useful comments. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. Y. C. Perrott acknowledges the support of a CCT/Cavendish Laboratory studentship.
- J. A. Tauber, N. Mandolesi, J. L. Puget et al., “Planck pre-launch status: the Planck mission,” Astronomy & Astrophysics, vol. 520, p. A1, 2010.
- P. A. R. Ade, N. Aghanim, M. Arnaud et al., “Planck early results. I. The Planck mission,” Astronomy & Astrophysics, vol. 536, article A1, 16 pages, 2011.
- P. A. R. Ade, N. Aghanim, M. Arnaud et al., “Planck early results. XX. New light on anomalous microwave emission from spinning dust grains,” Astronomy & Astrophysics, vol. 536, article A20, 17 pages, 2011.
- A. M. M. Scaife, N. Hurley-Walker, D. A. Green et al., “An excess of emission in the dark cloud LDN1111 with the Arcminute Microkelvin Imager,” Monthly Notices of the Royal Astronomical Society, vol. 394, no. 1, pp. L46–L50, 2009.
- A. M. M. Scaife, D. A. Green, G. G. Pooley et al., “High-resolution AMI Large Array imaging of spinning dust sources: spatially correlated 8 μm emission and evidence of a stellar wind in L675,” Monthly Notices of the Royal Astronomical Society, vol. 403, no. 1, pp. L46–L50, 2010.
- E. J. Murphy, G. Helou, J. J. Condon et al., “The detection of anomalous dust emission in the nearby galaxy NGC 6946,” The Astrophysical Journal, vol. 709, no. 2, p. L108, 2010.
- J. T. L. Zwart, R. W. Barker, P. Biddulph et al., “The arcminute microkelvin imager,” Monthly Notices of the Royal Astronomical Society, vol. 391, no. 4, pp. 1545–1558, 2008.
- D. J. Rudy, D. O. Muhleman, G. L. Berge, B. M. Jakosky, and P. R. Christensen, “Mars: VLA observations of the Northern Hemisphere and the north polar region at wavelengths of 2 and 6 cm,” Icarus, vol. 71, no. 1, pp. 159–177, 1987.
- A. B. Peck and A. J. Beasley, “A VLBA calibrator survey,” in Radio Emission from Galactic and Extragalactic Compact Sources: Proceedings of IAU Colloquium 164, Held in Socorro, New Mexico, USA, 21–26 April 1997, vol. 144, p. 155, 1998.
- P. C. Gregory, W. K. Scott, K. Douglas, and J. J. Condon, “The GB6 catalog of radio sources,” The Astrophysical Journal Supplement Series, vol. 103, pp. 427–432, 1996.
- S. Sharpless, “A catalogue of H II regions,” The Astrophysical Journal Supplement Series, vol. 4, p. 257, 1959.
- C. R. Kerton, J. Murphy, and J. Patterson, “A sharper view of the outer Galaxy at 1420 and 408 MHz from the Canadian Galactic Plane Survey II: the catalogue of extended radio sources,” Monthly Notices of the Royal Astronomical Society, vol. 379, no. 1, pp. 289–296, 2007.
- F. P. Israel and M. Felli, “Aperture synthesis observations of galactic H II regions. VIII—S106 and S235: regions of star formatiom,” Astronomy and Astrophysics,, vol. 63, no. 3, pp. 325–334, 1978.
- N. J. Evans II and G. N. Blair, “The energetics of molecular clouds. III—The S235 molecular cloud,” The Astrophysical Journal, vol. 246, pp. 394–408, 1981.
- L. Blitz, M. Fich, and A. A. Stark, “Catalog of CO radial velocities toward galactic H ii regions,” The Astrophysical Journal Supplement Series, vol. 49, pp. 183–206, 1982.
- J.-H. Kang, B.-C. Koo, and C. Salter, “An old supernova remnant within an H II complex at ℓ≈ 173°: FVW 172.8+1.5,” The Astronomical Journal, vol. 143, no. 3, article 75, 2012.
- A. M. M. Scaife, N. Hurley-Walker, M. L. Davies et al., “AMI limits on 15-GHz excess emission in northern H ii regions,” Monthly Notices of the Royal Astronomical Society, vol. 385, no. 2, pp. 809–822, 2008.
- C. Dickinson, R. D. Davies, and R. J. Davis, “Towards a free-free template for CMB foregrounds,” Monthly Notices of the Royal Astronomical Society, vol. 341, pp. 369–384, 2003.
- F. J. Lockman, “A survey of radio H II regions in the northern sky,” The Astrophysical Journal Supplement Series, vol. 71, pp. 469–479, 1989.
- L. Oster, “Free-free emission in the radio-frequency range,” Astronomical Journal, vol. 70, p. 50, 1961.
- D. Crampton and W. A. Fisher, “Spectroscopic observations of stars in H II regions,” Publications of the Dominion Astrophysical Observatory, Victoria, British Columbia, vol. 14, no. 12, pp. 283–304, 1974.
- C. A. Beichman, E. E. Becklin, and C. G. Wynn-Williams, “New multiple systems in molecular clouds,” The Astrophysical Journal, vol. 232, pp. L47–L51, 1979.
- M. Hayashi, T. Hasegawa, T. Omodaka, S. S. Hayashi, and R. Miyawaki, “The bright-rimmed molecular cloud around S140 IRS. II—bipolar outflow from S140 IRS 1,” The Astrophysical Journal, vol. 312, pp. 327–336, 1987.
- M. G. Hoare and J. Franco, “Massive star formation,” in Diffuse Matter from Star Forming Regions to Active Galaxies—A Volume Honouring John Dyson, T. W. Hartquist, J. M. Pittard, and S. A. E. G. Falle, Eds., Astrophysics and Space Science Proceedings, p. 61, Springer, Dordrecht, The Netherlands, 2007.
- W. J. Forrest and M. A. Shure, “Unipolar bubbles in star-forming regions,” The Astrophysical Journal, vol. 311, pp. L81–L84, 1986.
- M. G. Hoare, “An equatorial wind from the massive young stellar object S140 IRS 1,” The Astrophysical Journal, vol. 649, p. 856, 2006.
- Z. Jiang, M. Tamura, M. G. Hoare et al., “Disks around massive young stellar objects: are they common?” The Astrophysical Journal Letters, vol. 673, no. 2, pp. L175–L179, 2008.
- B. T. Draine and A. Lazarian, “Diffuse galactic emission from spinning dust grains,” The Astrophysical Journal Letters, vol. 494, no. 1, pp. L19–L22, 1998.
- M. P. Hobson and J. E. Baldwin, “Markov-chain monte carlo approach to the design of multilayer thin-film optical coatings,” Applied Optics, vol. 43, no. 13, pp. 2651–2660, 2004.
Copyright © 2013 Yvette C. Perrott et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.