Source author record

A. M. S. Richards

A. M. S. Richards appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2025arXiv

ATOMIUM: Continuum emission and evidence of dust enhancement from binary motion

Low- and intermediate-mass stars on the asymptotic giant branch (AGB) account for a significant portion of the dust and chemical enrichment in their host galaxy. Here we present ALMA observations of the continuum emission at 1.24 mm around a sample of 17 stars from the ATOMIUM survey. From our analysis of the stellar contributions to the continuum flux, we find that the semi-regular variables all have smaller physical radii and fainter monochromatic luminosities than the Mira variables. Comparing these properties with pulsation periods, we find a positive trend between stellar radius and period only for the Mira variables with periods above 300 days and a positive trend between the period and the monochromatic luminosity only for the red supergiants and the most extreme AGB stars with periods above 500 days. We find that the continuum emission at 1.2 mm can be classified into four groups. "Featureless" continuum emission is confined to the (unresolved) regions close to the star for five stars in our sample, relatively uniform extended flux is seen for four stars, tentative elongated features are seen for three stars, and the remaining five stars have unique or unusual morphological features in their continuum maps. These features can be explained by binary companions to 10 out of the 14 AGB stars in our sample. Based on our results we conclude that there are two modes of dust formation: well established pulsation-enhanced dust formation and our newly proposed companion-enhanced dust formation. If the companion is located close to the AGB star, in the wind acceleration region, then additional dust formed in the wake of the companion can increase the mass lost through the dust driven wind. This explains the different dust morphologies seen around our stars and partly accounts for a large scatter in literature mass-loss rates, especially among semiregular stars with small pulsation periods.

preprint2023arXiv

The Co-Ordinated Radio and Infrared Survey for High-Mass Star Formation. V. The CORNISH-South Survey and Catalogue

We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-South survey achieved a root mean square noise level of $\sim$ 0.11 mJy beam$^{-1}$, using the 6A configuration of the Australia Telescope Compact Array (ATCA). In this paper, we discuss the observations, data processing and measurements of the source properties. Above a 7$σ$ detection limit, 4701 sources were detected, and their ensemble properties show similar distributions with their northern counterparts. The catalogue is highly reliable and is complete to 90 per cent at a flux density level of 1.1 mJy. We developed a new way of measuring the integrated flux densities and angular sizes of non-Gaussian sources. The catalogue primarily provides positions, flux density measurements and angular sizes. All sources with IR counterparts at 8$μm$ have been visually classified, utilizing additional imaging data from optical, near-IR, mid-IR, far-IR and sub-millimetre galactic plane surveys. This has resulted in the detection of 524 H II regions of which 255 are ultra-compact H II regions, 287 planetary nebulae, 79 radio stars and 6 massive young stellar objects. The rest of the sources are likely to be extra-galactic. These data are particularly important in the characterization and population studies of compact ionized sources such as UCHII regions and PNe towards the Galactic mid-plane.

preprint2022arXiv

Self-calibration and improving image fidelity for ALMA and other radio interferometers

This manual is intended to help ALMA and other interferometer users improve images by recognising limitations and how to overcome them and deciding when and how to use self-calibration. The images provided by the ALMA Science Archive are calibrated using standard observing and data processing routines, including a quality assurance process to make sure that the observations meet the proposer's science requirements. This may not represent the full potential of the data, since any interferometry observation can be imaged with a range of resolutions and surface brightness sensitivity. The separation between phase calibration source and target usually limits the target dynamic range to a few hundred (or 50--100 for challenging conditions) but if the noise in the target field has not reached the thermal limit, improvements may be possible using self-calibration. This often requires judgements based on the target properties and is not yet automated for all situations. This manual provides background on the instrumental and atmospheric causes of visibility phase and amplitude errors, their effects on imaging and how to improve the signal to noise ratio and image fidelity by self-calibration. We introduce the conditions for self-calibration to be useful and how to estimate calibration parameter values for a range of observing modes (continuum, spectral line etc.). We also summarise more general error recognition and other techniques to tackle imaging problems. The examples are drawn from ALMA interferometric data processed using CASA, but the principles are generally applicable to most similar cm to sub-mm imaging.

preprint2022arXiv

Stringent limits on $^{28}$SiO maser emission from the recurrent nova T Coronae Borealis

There are indications that the third known eruption of the recurrent nova T CrB is imminent, and multi-wavelength observations prior to the eruption are important to characterise the system before it erupts. T CrB is known to display the SiO fundamental vibrational feature at 8$\,μ$m. When the anticipated eruption occurs, it is possible that the shock produced when the ejected material runs into the wind of the red giant in the system may be traced using SiO maser emission. We have used the 100m Effelsberg Radio Telescope to search for $^{28}$SiO emission in the $\upsilon=1$, $\upsilon=2$, $J=1\rightarrow0$ transitions, at 43.122 GHz and 42.820~GHz respectively, while the system is in quiescence. We find no evidence for such emission.

preprint2021arXiv

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst.

preprint2021arXiv

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.

preprint2021arXiv

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust.

preprint2020arXiv

ALMA and VLA reveal the lukewarm chromospheres of the nearby red supergiants Antares and Betelgeuse

We first present spatially resolved ALMA and VLA continuum observations of the early-M red supergiant Antares to search for the presence of a chromosphere at radio wavelengths. We resolve the free-free emission of the Antares atmosphere at 11 unique wavelengths between 0.7 mm (ALMA band 8) and 10 cm (VLA S band). The projected angular diameter is found to continually increase with increasing wavelength, from a low of 50.7 mas at 0.7 mm up to a diameter of 431 mas at 10 cm, which corresponds to 1.35 and 11.6 times the photospheric angular diameter, respectively. All four ALMA measurements show that the shape of the atmosphere is elongated, with a flattening of 15% at a similar position angle. The disk-averaged gas temperature of the atmosphere initially rises from a value of 2700 K at 1.35 $R_{\star}$ (i.e., 0.35 $R_{\star}$ above the photosphere) to a peak value of 3800 K at $\sim$2.5 $R_{\star}$, after which it then more gradually decreases to 1650 K at 11.6 $R_{\star}$. The rise in gas temperature between 1.35 $R_{\star}$ and $\sim$2.5 $R_{\star}$ is evidence for a chromospheric temperature rise above the photosphere of a red supergiant. We detect a clear change in the spectral index across the sampled wavelength range, with the flux density ${S_ν} \propto ν^{1.42}$ between 0.7 mm and 1.4 cm, which we associate with chromosphere-dominated emission, while the flux density ${S_ν} \propto ν^{0.8}$ between 4.3 cm and 10 cm, which we associate with wind-dominated emission. We then perform nonlocal thermal equilibrium modeling of the far-ultraviolet radiation field of another early-M red supergiant, Betelgeuse, and find that an additional hot (i.e., $>7000$ K) chromospheric photoionization component with a much smaller filling factor must also exist throughout the chromospheres of these stars.

preprint2020arXiv

An ALMA view of SO and SO$_2$ around oxygen-rich AGB stars

We present and analyse SO and SO$_2$, recently observed with high angular resolution and sensitivity in a spectral line survey with ALMA, for two oxygen-rich AGB stars: the low mass-loss rate R Dor and high mass-loss rate IK Tau. We analyse 8 lines of SO detected towards both stars, 78 lines of SO$_2$detected towards R Dor and 52 lines of SO$_2$ detected towards IK Tau. We detect several lines of $^{34}$SO, $^{33}$SO and $^{34}$SO$_2$ towards both stars, and tentatively S$^{18}$O towards R Dor, and hence derive isotopic ratios for these species. The spatially resolved observations show us that the two sulphur oxides are co-located towards R Dor and trace out the same wind structures in the circumstellar envelope (CSE). Much of the emission is well reproduced with a Gaussian abundance distribution spatially centred on the star. Emission from the higher energy levels of SO and SO$_2$ towards R Dor provide evidence in support of a rotating inner region of gas identified in earlier work. The new observations allow us to refine the abundance distribution of SO in IK Tau derived from prior observations with single antennas, and confirm the distribution is shell-like with the peak in the fractional abundance not centred on the star. The confirmation of different types of SO abundance distributions will help fine-tune chemical models and allows for an additional method to discriminate between low and high mass-loss rates for oxygen-rich AGB stars.

preprint2020arXiv

COBRaS: The e-MERLIN 21 cm Legacy survey of Cygnus OB2

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

preprint2020arXiv

The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5μ$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

preprint2020arXiv

The Mass-Loss History of the Red Hypergiant VY CMa

Imaging and spectroscopy of the knots, clumps, and extended arcs in the complex ejecta of VY CMa confirm a record of high mass loss events over the past few hundred years. HST/STIS spectroscopy of numerous small knots close to the star allow us to measure their radial velocities from the strong K I emission and determine their separate motions, spatial orientations, and time since ejecta. Their ages concentrate around 70, 120, 200 and 250 years ago. A K I emission knot only 50 mas from the star ejected as recently as 1985 -- 1995 may coincide with an H2O maser. Comparison with VY CMa's historic light curve from 1800 to the present, shows several knots with ejection times that correspond with extended periods of variability and deep minima. The similarity of this correspondence in VY CMa with the remarkable recent dimming of Betelgeuse and an outflow of gas is apparent. The evidence for similar outflows from the surface of a more typical red supergiant suggests that discrete ejections are more common and surface or convective activity is a major source of mass loss for red supergiants.

preprint2016arXiv

ALMA observations of the nearby AGB star L2 Puppis - I. Mass of the central star and detection of a candidate planet

Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System.

preprint2016arXiv

ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY CMa

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

preprint2015arXiv

ALMA Long Baseline Observations of the Strongly Lensed Submillimeter Galaxy HATLAS J090311.6+003906 at z=3.042

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

preprint2015arXiv

ALMA Observations of Asteroid 3 Juno at 60 Kilometer Resolution

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

preprint2015arXiv

ALMA observations of TiO$_2$ around VY Canis Majoris

Titanium dioxide, TiO$_2$, is a refractory species that could play a crucial role in the dust-condensation sequence around oxygen-rich evolved stars. To date, gas phase TiO$_2$ has been detected only in the complex environment of the red supergiant VY CMa. We aim to constrain the distribution and excitation of TiO$_2$ around VY CMa in order to clarify its role in dust formation. We analyse spectra and channel maps for TiO$_2$ extracted from ALMA science verification data. We detect 15 transitions of TiO$_2$, and spatially resolve the emission for the first time. The maps demonstrate a highly clumpy, anisotropic outflow in which the TiO$_2$ emission likely traces gas exposed to the stellar radiation field. A roughly east-west oriented, accelerating bipolar-like structure is found, of which the blue component runs into and breaks up around a solid continuum component. A distinct tail to the south-west is seen for some transitions, consistent with features seen in the optical and near-infrared. We find that a significant fraction of TiO$_2$ remains in the gas phase outside the dust-formation zone and suggest that this species might play only a minor role in the dust-condensation process around extreme oxygen-rich evolved stars like VY CMa.

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

preprint2015arXiv

First Results from High Angular Resolution ALMA Observations Toward the HL Tau Region

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($α$), which ranges from $α\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

preprint2015arXiv

Measuring Magnetic Fields Near and Far with the SKA via the Zeeman Effect

The measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.

preprint2015arXiv

The Physics of Water Masers observable with ALMA and SOFIA: Model Predictions for Evolved Stars

We present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space in number density, kinetic temperature and dust temperature. Results, in the form of maser optical depths, have been summarized in a master table, Table 6. Maser transitions identified in these models were grouped according to loci of inverted regions in the density/kinetic temperature plane, a property clearly related to the dominant mode of pumping. A more detailed study of the effect of dust temperature on maser optical depth enabled us to divide the maser transitions into three groups: those with both collisional and radiative pumping schemes (22,96,209,321,325,395,941 and 1486\,GHz), a much larger set that are predominantly radiatively pumped, and another large group with a predominantly collisional pump. The effect of accelerative and decelerative velocity shifts of up to 5km/s was found to be generally modest, with the primary effect of reducing computed maser optical depths. More subtle asymmetric effects, dependent on line overlap, include maximum gains offset from zero shift by >1km/s, but these effects were predominantly found under conditions of weak amplification. These models will allow astronomers to use multi-transition water maser observations to constrain physical conditions down to the size of individual masing clouds (size of a few astronomical units).

preprint2014arXiv

ALMA data suggest the presence of a spiral structure in the inner wind of CW Leo

(abbreviated) We aim to study the inner wind of the well-known AGB star CW Leo. Different diagnostics probing different geometrical scales have pointed toward a non-homogeneous mass-loss process: dust clumps are observed at milli-arcsec scale, a bipolar structure is seen at arcsecond-scale and multi-concentric shells are detected beyond 1". We present the first ALMA Cycle 0 band 9 data around 650 GHz. The full-resolution data have a spatial resolution of 0".42x0".24, allowing us to study the morpho-kinematical structure within ~6". Results: We have detected 25 molecular lines. The emission of all but one line is spatially resolved. The dust and molecular lines are centered around the continuum peak position. The dust emission has an asymmetric distribution with a central peak flux density of ~2 Jy. The molecular emission lines trace different regions in the wind acceleration region and suggest that the wind velocity increases rapidly from about 5 R* almost reaching the terminal velocity at ~11 R*. The channel maps for the brighter lines show a complex structure; specifically for the 13CO J=6-5 line different arcs are detected within the first few arcseconds. The curved structure present in the PV map of the 13CO J=6-5 line can be explained by a spiral structure in the inner wind, probably induced by a binary companion. From modeling the ALMA data, we deduce that the potential orbital axis for the binary system lies at a position angle of ~10-20 deg to the North-East and that the spiral structure is seen almost edge-on. We infer an orbital period of 55 yr and a binary separation of 25 au (or ~8.2 R*). We tentatively estimate that the companion is an unevolved low-mass main-sequence star. The ALMA data hence provide us for the first time with the crucial kinematical link between the dust clumps seen at milli-arcsecond scale and the almost concentric arcs seen at arcsecond scale.

preprint2014arXiv

ALMA Observations of Anisotropic Dust Mass-loss in the Inner Circumstellar Environment of the Red Supergiant VY Canis Majoris

The processes leading to dust formation and the subsequent role it plays in driving mass loss in cool evolved stars is an area of intense study. Here we present high resolution ALMA Science Verification data of the continuum emission around the highly evolved oxygen-rich red supergiant VY CMa. These data enable us to study the dust in its inner circumstellar environment at a spatial resolution of 129 mas at 321 GHz and 59 mas at 658 GHz, thus allowing us to trace dust on spatial scales down to 11 R$_{\star}$ (71 AU). Two prominent dust components are detected and resolved. The brightest dust component, C, is located 334 mas (61 R$_{\star}$) South East of the star and has a dust mass of at least $2.5\times 10^{-4}$ M$_{\odot}$. It has a dust emissivity spectral index of $β=-0.1$ at its peak, implying that it is optically thick at these frequencies with a cool core of $T_{d}\lesssim 100$ K. Interestingly, not a single molecule in the ALMA data has emission close to the peak of this massive dust clump. The other main dust component, VY, is located at the position of the star and contains a total dust mass of $4.0 \times 10^{-5} $M$_{\odot}$. It also contains a weaker dust feature extending over $60$ R$_{\star}$ to the North with the total component having a typical dust emissivity spectral index of $β=0.7$. We find that at least $17\%$ of the dust mass around VY CMa is located in clumps ejected within a more quiescent roughly spherical stellar wind, with a quiescent dust mass loss rate of $5 \times 10^{-6}$ M$_{\odot} $yr$^{-1}$. The anisotropic morphology of the dust indicates a continuous, directed mass loss over a few decades, suggesting that this mass loss cannot be driven by large convection cells alone.

preprint2014arXiv

ALMA sub-mm maser and dust distribution of VY Canis Majoris

Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries.

preprint2014arXiv

SKA studies of nearby galaxies: star-formation, accretion processes and molecular gas across all environments

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $μ$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

preprint2013arXiv

e-MERLIN resolves Betelgeuse at wavelength 5 cm

Convection, pulsation and magnetic fields have all been suggested as mechanisms for the transport of mass and energy from the optical photosphere of red supergiants, out to the region where the stellar wind is launched. We imaged the red supergiant Betelgeuse at 0.06-0.18 arcsec resolution, using e-MERLIN at 5.5--6.0 GHz, with a sensitivity of ~0.01 mJy/beam. Most of the radio emission comes from within an ellipse (0.235x0.218) arcsec^2 (~5x the optical radius), with a flux density of 1.62 mJy, giving an average brightness temperature ~1250 K. This radio photosphere contains two hotspots of 0.53 and 0.79 mJy/beam, separated by 90 milli-arcsec, with brightness temperatures 5400+/-600 K and 3800+/-500 K. Similar hotspots, at more than double the distance from the photosphere of those seen in any other regime, were detected by the less-sensitive `old' MERLIN in 1992, 1995 and 1996 and many exceed the photospheric temperature of 3600 K. Such brightness temperatures are high enough to emanate from pockets of chromospheric plasma. Other possibilities include local shock heating, the convective dredge-up of hot material or exceptionally cool, low density regions, transparent down to the hottest layer at ~40 milliarcsec radius. We also detect an arc 0.2--0.3 arcsec to the SW, brightness temperature ~150 K, in a similar direction to extensions seen on both smaller and larger scales in the infra-red and in CO at mm wavelengths. These preliminary results will be followed by further e-MERLIN, VLA and ALMA observations to help resolve the problem of mass elevation from 1 to 10 R* in red supergiants.

preprint2013arXiv

OH (1720 MHz) Masers: A Multiwavelength Study of the Interaction between the W51C Supernova Remnant and the W51B Star Forming Region

We present a comprehensive view of the W51B HII region complex and the W51C supernova remnant (SNR) using new radio observations from the VLA, VLBA, MERLIN, JCMT, and CSO along with archival data from Spitzer, ROSAT, ASCA, and Chandra. Our VLA data include the first 400 cm (74 MHz) continuum image of W51 at high resolution (88 arcsec). The 400 cm image shows non-thermal emission surrounding the G49.2-0.3 HII region, and a compact source of non-thermal emission (W51B_NT) coincident with the previously-identified OH (1720 MHz) maser spots, non-thermal 21 and 90 cm emission, and a hard X-ray source. W51B_NT falls within the region of high likelihood for the position of TeV gamma-ray emission. Using the VLBA three OH (1720 MHz) maser spots are detected in the vicinity of W51B_NT with sizes of 60 to 300 AU and Zeeman effect magnetic field strengths of 1.5 to 2.2 mG. The multiwavelength data demonstrate that the northern end of the W51B HII region complex has been partly enveloped by the advancing W51C SNR and this interaction explains the presence of W51B_NT and the OH masers. This interaction also appears in the thermal molecular gas which partially encircles W51B_NT and exhibits narrow pre-shock (DeltaV 5 km/s) and broad post-shock (DeltaV 20 km/s) velocity components. RADEX radiative transfer modeling of these two components yield physical conditions consistent with the passage of a non-dissociative C-type shock. Confirmation of the W51B/W51C interaction provides additional evidence in favor of this region being one of the best candidates for hadronic particle acceleration known thus far.

preprint2013arXiv

Polarization morphology of SiO masers in the circumstellar envelope of the AGB star R Cassiopeiae

Silicon monoxide maser emission has been detected in the circumstellar envelopes of many evolved stars in various vibrationally-excited rotational transitions. It is considered a good tracer of the wind dynamics close to the photosphere of the star. We have investigated the polarization morphology in the circumstellar envelope of an AGB star, R Cas. We mapped the linear and circular polarization of SiO masers in the v=1, J=1-0 transition. The linear polarization is typically a few tens of percent while the circular polarization is a few percent. The fractional polarization tends to be higher for emission of lower total intensity. We found that, in some isolated features the fractional linear polarization appears to exceed 100%. We found the Faraday rotation is not negligible but is ~15 deg., which could produce small scale structure in polarized emission whilst total intensity is smoother and partly resolved out. The polarization angles vary considerably from feature to feature but there is a tendency to favour the directions parallel or perpendicular to the radial direction with respect to the star. In some features, the polarization angle abruptly flips 90 deg. We found that our data are in the regime where the model of Goldreich et al (1973) can be applied and the polarization angle flip is caused when the magnetic field is at close to 55 deg. to the line of sight. The polarization angle configuration is consistent with a radial magnetic field although other configurations are not excluded.

preprint2013arXiv

Probing the accelerating Universe with radio weak lensing in the JVLA Sky Survey

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

preprint2012arXiv

Evolved star water maser cloud size determined by star size

Cool, evolved stars undergo copious mass loss but the details of how the matter is returned to the ISM are still under debate. We investigated the structure and evolution of the wind at 5 to 50 stellar radii from Asymptotic Giant Branch and Red Supergiant stars. 22-GHz water masers around seven evolved stars were imaged using MERLIN, at sub-AU resolution. Each source was observed at between 2 and 7 epochs (several stellar periods). We compared our results with long-term Pushchino single dish monitoring. The 22-GHz emission is located in ~spherical, thick, unevenly filled shells. The outflow velocity doubles between the inner and outer shell limits. Water maser clumps could be matched at successive epochs separated by <2 years for AGB stars, or at least 5 years for RSG. This is much shorter than the decades taken for the wind to cross the maser shell, and comparison with spectral monitoring shows that some features fade and reappear. In 5 sources, most of the matched features brighten or dim in concert from one epoch to the next. One cloud in W Hya was caught in the act of passing in front of a background cloud leading to 50-fold, transient amplification. The masing clouds are 1-2 orders of magnitude denser than the wind average and contain a substantial fraction of the mass loss in this region, with a filling factor <1%. The RSG clouds are ~10x bigger than those round the AGB stars. Proper motions are dominated by expansion, with no systematic rotation. The maser clouds survive for decades (the shell crossing time) but the masers are not always beamed in our direction. Radiative effects cause changes in flux density throughout the maser shells on short timescales. Cloud size is proportional to parent star size; clouds have a similar radius to the star in the 22-GHz maser shell. Stellar properties such as convection cells must determine the clumping scale.

preprint2011arXiv

Flares and Proper Motions of Ground-State OH Masers in W75N

The star-forming region W75N hosts bright OH masers that are observed to be variable. We present observations taken in 2008 of the ground-state OH maser transitions with the Very Long Baseline Array (VLBA) and the Multi-Element Radio-Linked Interferometer Network (MERLIN) and with the Nancay Radio Telescope in 2011. Several of the masers in W75N were observed to be flaring, with the brightest 1720-MHz maser in excess of 400 Jy. The 1720-MHz masers appear to be associated with the continuum source VLA 1, unlike the bright flaring 1665- and 1667-MHz masers, which are associated with VLA 2. The 1720-MHz masers are located in an outflow traced by water masers and are indicative of very dense molecular material near the H II region. The magnetic field strengths are larger in the 1720-MHz maser region than in most regions hosting only main-line OH masers. The density falls off along the outflow, and the order of appearance of different transitions of OH masers is consistent with theoretical models. The 1665- and 1667-MHz VLBA data are compared against previous epochs over a time baseline of over 7 years. The median maser motion is 3.5 km/s, with a scatter that is comparable to thermal turbulence. The general pattern of maser proper motions observed in the 1665- and 1667-MHz transitions is consistent with previous observations.

preprint2011arXiv

The 43GHz SiO maser in the circumstellar envelope of the AGB star R Cassiopeiae

We present multi-epoch, total intensity, high-resolution images of 43GHz, v=1, J=1-0 SiO maser emission toward the Mira variable R Cas. In total we have 23 epochs of data for R Cas at approximate monthly intervals over an optical pulsation phase range from 0.158 to 1.78. These maps show a ring-like distribution of the maser features in a shell, which is assumed to be centred on the star at a radius of 1.6 to 2.3 times the stellar radii. It is clear from these images that the maser emission is significantly extended around the star. At some epochs a faint outer arc can be seen at 2.2 stellar radii. The intensity of the emission waxes and wanes during the stellar phase. Some maser features are seen infalling as well as outflowing. We have made initial comparisons of our data with models by Gray et. al. (2009).

preprint2010arXiv

Discrete Source Survey of 6 GHz OH emission from PNe & pPNe and first 6 GHz images of K 3-35

The aim of this study is to investigate the physical properties of molecular envelopes of planetary nebulae in their earliest stages of evolution. Using the 100m telescope at Effelsberg, we have undertaken a high sensitivity discrete source survey for the first excited state of OH maser emission (J=5/2, 2PI3/2 at 6GHz) in the direction of planetary and proto-planetary nebulae exhibiting 18cm OH emission (main and/or satellite lines), and we further validate our detections using the Nançay radio telescope at 1.6-1.7GHz and MERLIN interferometer at 1.6-1.7 and 6GHz. Two sources have been detected at 6035MHz (5cm), both of them are young (or very young) planetary nebulae. The first one is a confirmation of the detection of a weak 6035MHz line in Vy 2-2. The second one is a new detection, in K 3-35, which was already known to be an exceptional late type star because it exhibits 1720MHz OH emission. The detection of 6035MHz OH maser emission is confirmed by subsequent observations made with the MERLIN interferometer. These lines are very rarely found in evolved stars. The 1612MHz masers surround but are offset from the 1720 and 6035MHz masers which in turn lie close to a compact 22GHz continuum source embedded in the optical nebula.

preprint2010arXiv

Observational evidence for the shrinking of bright maser spots

The nature of maser emission means that the apparent angular size of an individual maser spot is determined by the amplification process as well as by the instrinsic size of the emitting cloud. Highly sensitive MERLIN radio interferometry images spatially and spectrally resolve water maser clouds around evolved stars. We measured the properties of clouds around the red supergiant S Per and the AGB stars IK Tau, RT Vir, U Her and U Ori, to test maser beaming theory. Spherical clouds are expected to produce an inverse relationship between maser intensity and apparent size, which would not be seen from cylindrical or slab-like regions. We analysed the maser properties, in order to estimate the saturation state, and investigated the variation of observed spot size with intensity and across the spectral line profiles. Circumstellar masers emanate from discrete clouds from about one to 20 AU in diameter depending on the star. Most of the maser features have negative excitation temperatures close to zero and modest optical depths, showing that they are mainly unsaturated. Around S Per and (at most epochs) RT Vir and IK Tau, the maser component size shrinks with increasing intensity. In contrast, the masers around U Ori and U Her tend to increase in size, with a larger scatter. The water masers from S Per, RT Vir and IK Tau are mainly beamed into spots with an observed angular size much smaller than the emitting clouds and smallest of all at the line peaks. This suggests that the masers are amplification-bounded, emanating from approximately spherical clouds. Many of the masers around U Her and U Ori have apparent sizes which are more similar to the emitting clouds and have less or no dependence on intensity, suggesting that these masers are matter-bounded. This is consistent with an origin in flattened clouds and these two stars have shown other behaviour indicating the presence of shocks.

preprint2010arXiv

Quasi-Periodic Formaldehyde Maser Flares in the Massive Protostellar Object IRAS18566+0408

We report results of an extensive observational campaign of the 6 cm formaldehyde maser in the young massive stellar object IRAS18566+0408 (G37.55+0.20) conducted from 2002 to 2009. Using Arecibo, VLA, and GBT, we discovered quasi-periodic formaldehyde flares (P ~ 237 days). Based on Arecibo observations, we also discovered correlated variability between formaldehyde (H2CO) and methanol (CH3OH) masers. The H2CO and CH3OH masers are not spatially coincident, as demonstrated by different line velocities and high angular resolution MERLIN observations. The flares could be caused by variations in the infrared radiation field, possibly modulated by periodic accretion onto a young binary system.

preprint2009arXiv

The diversity of methanol maser morphologies from VLBI observations

We investigate which structures the 6.7 GHz methanol masers trace in the environment of high-mass protostar candidates by observing a homogenous sample of methanol masers selected from Torun surveys. We also probed their origins by looking for associated H II regions and IR emission. We selected 30 methanol sources with improved position accuracies achieved using MERLIN and another 3 from the literature. We imaged 31 of these using the European VLBI Network's expanded array of telescopes with 5-cm (6-GHz) receivers. We used the VLA to search for 8.4 GHz radio continuum counterparts and inspected Spitzer GLIMPSE data at 3.6-8 um from the archive. High angular resolution images allowed us to analyze the morphology and kinematics of the methanol masers in great detail and verify their association with radio continuum and mid-infrared emission. A new class of "ring-like" methanol masers in star--forming regions appeared to be suprisingly common, 29 % of the sample. The new morphology strongly suggests that methanol masers originate in the disc or torus around a proto- or a young massive star. However, the maser kinematics indicate the strong influence of outflow or infall. This suggests that they form at the interface between the disc/torus and a flow. This is also strongly supported by Spitzer results because the majority of the masers coincide with 4.5 um emission to within less than 1 arcsec. Only four masers are associated with the central parts of UC H II regions. This implies that 6.7 GHz methanol maser emission occurs before H II region observable at cm wavelengths is formed.

preprint2006arXiv

High-z radio starbursts host obscured X-ray AGN

We use Virtual Observatory methods to investigate the association between radio and X-ray emission at high redshifts. Fifty-five of the 92 HDF(N) sources resolved by combining MERLIN+VLA data were detected by Chandra, of which 18 are hard enough and bright enough to be obscured AGN. The high-z population of microJy radio sources is dominated by starbursts an order of magnitude more active and more extended than any found at z<1 and at least a quarter of these simultaneously host highly X-ray-luminous obscured AGN.

preprint2002arXiv

Discovery of the host galaxy of HDF850.1, the brightest sub-mm source in the Hubble Deep Field

Despite extensive observational efforts, the brightest sub--mm source in the Hubble Deep Field, HDF850.1, has failed to yield a convincing optical/infrared identification almost 4 years after its discovery. This failure is all the more notable given the availability of supporting multi-frequency data of unparalleled depth, and sub-arcsec positional accuracy for the sub-mm/mm source. Consequently, HDF850.1 has become a test case of the possibility that the most violently star-forming objects in the universe are too red and/or distant to be seen in the deepest optical images. Here we report the discovery of the host galaxy of HDF850.1. This object has been revealed by careful analysis of a new, deep K-prime image of the HDF obtained with the Subaru 8.2-m telescope. Its reality is confirmed by a similar analysis of the HST NICMOS F160W image of the same region. This object is extremely faint (K=23.5), clumpy (on sub-arcsec scales) and very red (I-K > 5.2; H-K = 1.4 +- 0.35). The likelihood that it is the correct identification is strongly reinforced by a reanalysis of the combined MERLIN+VLA 1.4-GHz map of the field which yields a new radio detection of HDF850.1 only 0.1 arcsec from the new near-ir counterpart, and with sufficient positional accuracy to exclude all previously considered alternative optical candidates. We have calculated new confidence limits on the estimated redshift of HDF850.1 and find z = 4.1 +- 0.5. We also calculate that the flux density of HDF850.1 has been boosted by a factor of ~3 through lensing by the intervening elliptical 3-586.0, consistent with predictions that a small but significant fraction of blank-field sub-mm sources are lensed by foreground galaxies. We discuss the wider implications of these results for the sub-mm population and cosmic star-formation history.