Source author record

I. McDonald

I. McDonald 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

33works
6topics
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

33 published item(s)

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid's unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10's stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10's orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2026arXiv

The Gaia All-Sky Stellar Parameters Service (GASPS)

Temperature and luminosity are the two key diagnostics of a star, yet these cannot come directly from survey data, but must be imputed by comparing those data to models. SED fitting offers a high-precision method to obtain both parameters for stars where both their distance and extinction are well known. The recent publication of many all-sky or large-area surveys coincides the publication of parallaxes and 3D extinction cubes from the Gaia satellite, making it possible to perform SED fitting of truly large ($>10^8$) numbers of Galactic stars for the first time. The analysis of this data requires a high level of automation. Here, we describe the ongoing Gaia All-Sky Stellar Parameters Service (GASPS): the fitting of 240 million SEDs from Gaia DR3 and the extraction of temperatures and luminosities for the corresponding stars using the PySSED code. We demonstrate the quality of the initial results, and the promise that these data show, from wavelength-specific information such as the ultraviolet and infrared excess of each star, to stellar classification, to expansion of the project beyond our own Galaxy, and mineralogical mapping of the Milky Way's interstellar medium.

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.

preprint2025arXiv

The Nearby Evolved Stars Survey III: First data release of JCMT CO-line observations

Low- to intermediate-mass ($\sim$0.8$-$8 M$_\odot$) evolved stars contribute significantly to the chemical enrichment of the interstellar medium in the local Universe, making accurate mass-return estimates in their final stages crucial. The Nearby Evolved Stars Survey (NESS) is a large multi-telescope project targeting a volume-limited sample of $\sim$850 stars within 3 kpc in order to derive the dust and gas return rates in the Solar Neighbourhood, and to constrain the physics underlying these processes. We present an initial analysis of the CO-line observations, including detection statistics, carbon isotopic ratios, initial mass-loss rates, and gas-to-dust ratios. We describe a new data reduction pipeline to analyse the available NESS CO data from the JCMT, measuring line parameters and calculating empirical gas mass-loss rates. We present the first release of the available data on 485 sources, one of the largest homogeneous samples of CO data to date. Comparison with a large literature sample finds that high mass-loss rate and especially carbon-rich sources are over-represented in literature, while NESS is probing significantly more sources at low mass-loss rates, detecting 59 sources in CO for the first time and providing useful upper limits. CO line detection rates are 81% for the CO (2--1) line and 75% for CO (3--2). The majority (82%) of detected lines conform to the expected soft parabola shape, while eleven sources show a double wind. Calculated mass-loss rates show power-law relations with both the dust-production rates and expansion velocities up to $\sim 5 \times 10^{-6}$~\msunyr. Median gas-to-dust ratios of 250 and 680 are found for oxygen-rich and carbon-rich sources, respectively. Our analysis of CO observations in this first data release highlights the importance of our volume-limited approach in characterizing the local AGB population as a whole.

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.

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

Optimizing exoplanet atmosphere retrieval using unsupervised machine-learning classification

One of the principal bottlenecks to atmosphere characterisation in the era of all-sky surveys is the availability of fast, autonomous and robust atmospheric retrieval methods. We present a new approach using unsupervised machine learning to generate informed priors for retrieval of exoplanetary atmosphere parameters from transmission spectra. We use principal component analysis (PCA) to efficiently compress the information content of a library of transmission spectra forward models generated using the PLATON package. We then apply a $k$-means clustering algorithm in PCA space to segregate the library into discrete classes. We show that our classifier is almost always able to instantaneously place a previously unseen spectrum into the correct class, for low-to-moderate spectral resolutions, $R$, in the range $R~=~30-300$ and noise levels up to $10$~per~cent of the peak-to-trough spectrum amplitude. The distribution of physical parameters for all members of the class therefore provides an informed prior for standard retrieval methods such as nested sampling. We benchmark our informed-prior approach against a standard uniform-prior nested sampler, finding that our approach is up to a factor two faster, with negligible reduction in accuracy. We demonstrate the application of this method to existing and near-future observatories, and show that it is suitable for real-world application. Our general approach is not specific to transmission spectroscopy and should be more widely applicable to cases that involve repetitive fitting of trusted high-dimensional models to large data catalogues, including beyond exoplanetary science.

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

The ALMA detection of CO rotational line emission in AGB stars in the Large Magellanic Cloud

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

preprint2016arXiv

The infrared spectral properties of Magellanic carbon stars

The Infrared Spectrograph on the Spitzer Space Telescope observed 184 carbon stars in the Magellanic Clouds. This sample reveals that the dust-production rate (DPR) from carbon stars generally increases with the pulsation period of the star. The composition of the dust grains follows two condensation sequences, with more SiC condensing before amorphous carbon in metal-rich stars, and the order reversed in metal-poor stars. MgS dust condenses in optically thicker dust shells, and its condensation is delayed in more metal-poor stars. Metal-poor carbon stars also tend to have stronger absorption from C2H2 at 7.5 um. The relation between DPR and pulsation period shows significant apparent scatter, which results from the initial mass of the star, with more massive stars occupying a sequence parallel to lower-mass stars, but shifted to longer periods. Accounting for differences in the mass distribution between the carbon stars observed in the Small and Large Magellanic Clouds reveals a hint of a subtle decrease in the DPR at lower metallicities, but it is not statistically significant. The most deeply embedded carbon stars have lower variability amplitudes and show SiC in absorption. In some cases they have bluer colors at shorter wavelengths, suggesting that the central star is becoming visible. These deeply embedded stars may be evolving off of the AGB and/or they may have non-spherical dust geometries.

preprint2015arXiv

Bulge RR Lyrae stars in the VVV tile $\textit{b201}$

The VISTA Variables in the Vía Láctea (VVV) Survey is one of the six ESO public surveys currently ongoing at the VISTA telescope on Cerro Paranal, Chile. VVV uses near-IR ($ZYJHK_{\rm s}$) filters that at present provide photometry to a depth of $K_{\rm s} \sim 17.0$ mag in up to 36 epochs spanning over four years, and aim at discovering more than 10$^6$ variable sources as well as trace the structure of the Galactic bulge and part of the southern disk. A variability search was performed to find RR Lyrae variable stars. The low stellar density of the VVV tile $\textit{b201}$, which is centered at ($\ell, b$) $\sim$ ($-9^\circ, -9^\circ$), makes it suitable to search for variable stars. Previous studies have identified some RR Lyrae stars using optical bands that served to test our search procedure. The main goal is to measure the reddening, interstellar extinction, and distances of the RR Lyrae stars and to study their distribution on the Milky Way bulge. A total of 1.5 sq deg were analyzed, and we found 39 RR Lyrae stars, 27 of which belong to the ab-type and 12 to the c-type. Our analysis recovers all the previously identified RR Lyrae variables in the field and discovers 29 new RR Lyrae stars. The reddening and extinction toward all the RRab stars in this tile were derived, and distance estimations were obtained through the period--luminosity relation. Despite the limited amount of RR Lyrae stars studied, our results are consistent with a spheroidal or central distribution around $\sim 8.1$ and $\sim 8.5$ kpc. for either the Cardelli or Nishiyama extinction law.

preprint2015arXiv

EU Del: exploring the onset of pulsation-driven winds in giant stars

We explore the wind-driving mechanism of giant stars through the nearby (117 pc), intermediate-luminosity ($L \approx 1600$ L$_\odot$) star EU Del (HIP 101810, HD 196610). APEX observations of the CO (3--2) and (2--1) transitions are used to derive a wind velocity of 9.51 $\pm$ 0.02 km s$^{-1}$, a $^{12}$C/$^{13}$C ratio of 14 $^{+9}_{-4}$, and a mass-loss rate of a few $\times$ 10$^{-8}$ M$_\odot$ yr$^{-1}$. From published spectra, we estimate that the star has a metallicity of [Fe/H] = --0.27 $\pm$ $\sim$0.30 dex. The star's dusty envelope lacks a clear 10-$μ$m silicate feature, despite the star's oxygen-rich nature. Radiative transfer modelling cannot fit a wind acceleration model which relies solely on radiation pressure on condensing dust. We compare our results to VY Leo (HIP 53449), a star with similar temperature and luminosity, but different pulsation properties. We suggest the much stronger mass loss from EU Del may be driven by long-period stellar pulsations, due to its potentially lower mass. We explore the implications for the mass-loss rate and wind velocities of other stars.

preprint2015arXiv

Photometric properties of carbon stars in the Small Magellanic Cloud

The Optical Gravitational Lensing Experiment identified over 1,800 carbon-rich Mira and semi-regular variables in the Small Magellanic Cloud. Multi-epoch infrared photometry reveals that the semi-regulars and Miras follow different sequences in color-color space when using colors sensitive to molecular absorption bands. The dustiest Miras have the strongest pulsation amplitudes and longest periods. Efforts to determine bolometric magnitudes reveal possible systematic errors with published bolometric corrections.

preprint2015arXiv

Spitzer Infrared Spectrographic point source classification in the Small Magellanic Cloud

The Magellanic clouds are uniquely placed to study the stellar contribution to dust emission. Individual stars can be resolved in these systems even in the mid-infrared, and they are close enough to allow detection of infrared excess caused by dust.We have searched the Spitzer Space Telescope data archive for all Infrared Spectrograph (IRS) staring-mode observations of the Small Magellanic Cloud (SMC) and found that 209 Infrared Array Camera (IRAC) point sources within the footprint of the Surveying the Agents of Galaxy Evolution in the Small Magellanic Cloud (SAGE-SMC) Spitzer Legacy programme were targeted, within a total of 311 staring mode observations. We classify these point sources using a decision tree method of object classification, based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership and variability information. We find 58 asymptotic giant branch (AGB) stars, 51 young stellar objects (YSOs), 4 post-AGB objects, 22 Red Supergiants (RSGs), 27 stars (of which 23 are dusty OB stars), 24 planetary nebulae (PNe), 10Wolf-Rayet (WR) stars, 3 Hii regions, 3 R Coronae Borealis (R CrB) stars, 1 Blue Supergiant and 6 other objects, including 2 foreground AGB stars. We use these classifications to evaluate the success of photometric classification methods reported in the literature.

preprint2014arXiv

A Spitzer Space Telescope survey of extreme Asymptotic Giant Branch stars in M32

We investigate the population of cool, evolved stars in the Local Group dwarf elliptical galaxy M32, using Infrared Array Camera observations from the Spitzer Space Telescope. We construct deep mid-infrared colour-magnitude diagrams for the resolved stellar populations within 3.5 arcmin of M32's centre, and identify those stars that exhibit infrared excess. Our data is dominated by a population of luminous, dust-producing stars on the asymptotic giant branch (AGB) and extend to approximately 3 mag below the AGB tip. We detect for the first time a sizeable population of `extreme' AGB stars, highly enshrouded by circumstellar dust and likely completely obscured at optical wavelengths. The total dust-injection rate from the extreme AGB candidates is measured to be $7.5 \times 10^{-7}$ ${\rm M}_{\odot} \, {\rm yr}^{-1}$, corresponding to a gas mass-loss rate of $1.5 \times 10^{-4}$ ${\rm M}_{\odot} \, {\rm yr}^{-1}$. These extreme stars may be indicative of an extended star-formation epoch between 0.2 and 5 Gyr ago.

preprint2014arXiv

A study of rotating globular clusters - the case of the old, metal-poor globular cluster NGC 4372

Aims: We present the first in-depth study of the kinematic properties and derive the structural parameters of NGC 4372 based on the fit of a Plummer profile and a rotating, physical model. We explore the link between internal rotation to different cluster properties and together with similar studies of more GCs, we put these in the context of globular cluster formation and evolution. Methods: We present radial velocities for 131 cluster member stars measured from high-resolution FLAMES/GIRAFFE observations. Their membership to the GC is additionally confirmed from precise metallicity estimates. Using this kinematic data set we build a velocity dispersion profile and a systemic rotation curve. Additionally, we obtain an elliptical number density profile of NGC 4372 based on optical images using a MCMC fitting algorithm. From this we derive the cluster's half-light radius and ellipticity as r_h=3.4'+/-0.04' and e=0.08+/-0.01. Finally, we give a physical interpretation of the observed morphological and kinematic properties of this GC by fitting an axisymmetric, differentially rotating, dynamical model. Results: Our results show that NGC 4372 has an unusually high ratio of rotation amplitude to velocity dispersion (1.2 vs. 4.5 km/s) for its metallicity. This, however, puts it in line with two other exceptional, very metal-poor GCs - M 15 and NGC 4590. We also find a mild flattening of NGC 4372 in the direction of its rotation. Given its old age, this suggests that the flattening is indeed caused by the systemic rotation rather than tidal interactions with the Galaxy. Additionally, we estimate the dynamical mass of the GC M_dyn=2.0+/-0.5 x 10^5 M_Sun based on the dynamical model, which constrains the mass-to-light ratio of NGC 4372 between 1.4 and 2.3 M_Sun/L_Sun, representative of an old, purely stellar population.

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

An Infrared Census of DUST in Nearby Galaxies with Spitzer (DUSTiNGS), II. Discovery of Metal-poor Dusty AGB Stars

The DUSTiNGS survey (DUST in Nearby Galaxies with Spitzer) is a 3.6 and 4.5 micron imaging survey of 50 nearby dwarf galaxies designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. Using two epochs, spaced approximately six months apart, we identify a total of 526 dusty variable AGB stars (sometimes called "extreme" or x-AGB stars; [3.6]-[4.5]>0.1 mag). Of these, 111 are in galaxies with [Fe/H] < -1.5 and 12 are in galaxies with [Fe/H] < -2.0, making them the most metal-poor dust-producing AGB stars known. We compare these identifications to those in the literature and find that most are newly discovered large-amplitude variables, with the exception of approximately 30 stars in NGC 185 and NGC 147, one star in IC 1613, and one star in Phoenix. The chemical abundances of the x-AGB variables are unknown, but the low metallicities suggest that they are more likely to be carbon-rich than oxygen-rich and comparisons with existing optical and near-IR photometry confirms that 70 of the x-AGB variables are confirmed or likely carbon stars. We see an increase in the pulsation amplitude with increased dust production, supporting previous studies suggesting that dust production and pulsation are linked. We find no strong evidence linking dust production with metallicity, indicating that dust can form in very metal-poor environments.

preprint2014arXiv

ExELS: an exoplanet legacy science proposal for the ESA Euclid mission. II. Hot exoplanets and sub-stellar systems

The Exoplanet Euclid Legacy Survey (ExELS) proposes to determine the frequency of cold exoplanets down to Earth mass from host separations of ~1 AU out to the free-floating regime by detecting microlensing events in Galactic Bulge. We show that ExELS can also detect large numbers of hot, transiting exoplanets in the same population. The combined microlensing+transit survey would allow the first self-consistent estimate of the relative frequencies of hot and cold sub-stellar companions, reducing biases in comparing "near-field" radial velocity and transiting exoplanets with "far-field" microlensing exoplanets. The age of the Bulge and its spread in metallicity further allows ExELS to better constrain both the variation of companion frequency with metallicity and statistically explore the strength of star-planet tides. We conservatively estimate that ExELS will detect ~4100 sub-stellar objects, with sensitivity typically reaching down to Neptune-mass planets. Of these, ~600 will be detectable in both Euclid's VIS (optical) channel and NISP H-band imager, with ~90% of detections being hot Jupiters. Likely scenarios predict a range of 2900-7000 for VIS and 400-1600 for H-band. Twice as many can be expected in VIS if the cadence can be increased to match the 20-minute H-band cadence. The separation of planets from brown dwarfs via Doppler boosting or ellipsoidal variability will be possible in a handful of cases. Radial velocity confirmation should be possible in some cases, using 30-metre-class telescopes. We expect secondary eclipses, and reflection and emission from planets to be detectable in up to ~100 systems in both VIS and NISP-H. Transits of ~500 planetary-radius companions will be characterised with two-colour photometry and ~40 with four-colour photometry (VIS,YJH), and the albedo of (and emission from) a large sample of hot Jupiters in the H-band can be explored statistically.

preprint2014arXiv

Herschel Planetary Nebula Survey (HerPlaNS) - First Detection of OH+ in Planetary Nebulae

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$μ$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$μ$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

preprint2014arXiv

Modelling the alumina abundance of oxygen-rich evolved stars in the Large Magellanic Cloud

In order to determine the composition of the dust in the circumstellar envelopes of oxygen-rich asymptotic giant branch (AGB) stars we have computed a grid of modust radiative-transfer models for a range of dust compositions, mass-loss rates, dust shell inner radii and stellar parameters. We compare the resulting colours with the observed oxygen-rich AGB stars from the SAGE-Spec Large Magellanic Cloud (LMC) sample, finding good overall agreement for stars with a mid-infrared excess. We use these models to fit a sample of 37 O-rich AGB stars in the LMC with optically thin circumstellar envelopes, for which 5$-$35-$μ$m Spitzer infrared spectrograph (IRS) spectra and broadband photometry from the optical to the mid-infrared are available. From the modelling, we find mass-loss rates in the range $\sim 8\times10^{-8}$ to $5\times10^{-6}$ M$_{\odot}\ \mathrm{yr}^{-1}$, and we show that a grain mixture consisting primarily of amorphous silicates, with contributions from amorphous alumina and metallic iron provides a good fit to the observed spectra. Furthermore, we show from dust models that the AKARI [11]$-$[15] versus [3.2]$-$[7] colour-colour diagram, is able to determine the fractional abundance of alumina in O-rich AGB stars.

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

ExELS: an exoplanet legacy science proposal for the ESA Euclid mission I. Cold exoplanets

Euclid is the second M-class mission of the ESA Cosmic Vision programme, with the principal science goal of studying dark energy. Euclid is also expected to undertake additional Legacy Science programmes. One proposal is the Exoplanet Euclid Legacy Survey (ExELS) which will be the first survey able to measure the abundance of exoplanets down to Earth mass for host separations from ~1AU out to the free-floating (unbound) regime. The cold and free-floating exoplanet regimes represent a crucial discovery space for testing planet formation theories. ExELS will use the gravitational microlensing technique and will detect 1000 microlensing events per month over 1.6 deg^2 of the Galactic bulge. We assess how many of these events will have detectable planetary signatures using a detailed multi-wavelength microlensing simulator (MABuLS) which incorporates the Besancon Galactic model with 3D extinction. MABuLS is the first theoretical simulation of microlensing to treat the effects of point spread function (PSF) blending self-consistently with the underlying Galactic model. We use MABuLS, together with current numerical models for the Euclid PSFs, to explore a number of designs and de-scope options for ExELS, including the exoplanet yield as a function of filter choice and slewing time, and the effect of systematic photometry errors. Using conservative extrapolations of current empirical exoplanet mass functions determined from ground-based microlensing and radial velocity surveys, ExELS can expect to detect a few hundred cold exoplanets around mainly G, K and M-type stellar hosts, including ~45 Earth-mass planets and ~6 Mars-mass planets for an observing programme totalling 10 months. ExELS will be capable of measuring the cold exoplanet mass function down to Earth mass or below, with orbital separations from ~1AU to the free-floating regime. (Abridged)

preprint2012arXiv

Carbon-rich dust production in metal-poor galaxies in the Local Group

We have observed a sample of 19 carbon stars in the Sculptor, Carina, Fornax, and Leo I dwarf spheroidal galaxies with the Infrared Spectrograph on the Spitzer Space Telescope. The spectra show significant quantities of dust around the carbon stars in Sculptor, Fornax, and Leo I, but little in Carina. Previous comparisons of carbon stars with similar pulsation properties in the Galaxy and the Magellanic Clouds revealed no evidence that metallicity affected the production of dust by carbon stars. However, the more metal-poor stars in the current sample appear to be generating less dust. These data extend two known trends to lower metallicities. In more metal-poor samples, the SiC dust emission weakens, while the acetylene absorption strengthens. The bolometric magnitudes and infrared spectral properties of the carbon stars in Fornax are consistent with metallicities more similar to carbon stars in the Magellanic Clouds than in the other dwarf spheroidals in our sample. A study of the carbon budget in these stars reinforces previous considerations that the dredge-up of sufficient quantities of carbon from the stellar cores may trigger the final superwind phase, ending a star's lifetime on the asymptotic giant branch.

preprint2012arXiv

On the metallicity dependence of crystalline silicates in oxygen-rich asymptotic giant branch stars and red supergiants

We investigate the occurrence of crystalline silicates in oxygen-rich evolved stars across a range of metallicities and mass-loss rates. It has been suggested that the crystalline silicate feature strength increases with increasing mass-loss rate, implying a correlation between lattice structure and wind density. To test this, we analyse Spitzer IRS and Infrared Space Observatory SWS spectra of 217 oxygen-rich asymptotic giant branch stars and 98 red supergiants in the Milky Way, the Large and Small Magellanic Clouds and Galactic globular clusters. These encompass a range of spectral morphologies from the spectrally-rich which exhibit a wealth of crystalline and amorphous silicate features to 'naked' (dust-free) stars. We combine spectroscopic and photometric observations with the GRAMS grid of radiative transfer models to derive (dust) mass-loss rates and temperature. We then measure the strength of the crystalline silicate bands at 23, 28 and 33 microns. We detect crystalline silicates in stars with dust mass-loss rates which span over 3 dex, down to rates of ~10^-9 solar masses/year. Detections of crystalline silicates are more prevalent in higher mass-loss rate objects, though the highest mass-loss rate objects do not show the 23-micron feature, possibly due to the low temperature of the forsterite grains or it may indicate that the 23-micron band is going into absorption due to high column density. Furthermore, we detect a change in the crystalline silicate mineralogy with metallicity, with enstatite seen increasingly at low metallicity.

preprint2012arXiv

The Dust Budget of the SMC: Are AGB Stars the Primary Dust Source at Low Metallicity?

We estimate the total dust input from the cool evolved stars in the Small Magellanic Cloud (SMC), using the 8 micron excess emission as a proxy for the dust-production rate. We find that Asymptotic Giant Branch (AGB) and red supergiant (RSG) stars produce (8.6-9.5) x 10^7 solar masses per year of dust, depending on the fraction of far-infrared sources that belong to the evolved star population (with 10%-50% uncertainty in individual dust-production rates). RSGs contribute the least (<4%), while carbon-rich AGB stars (especially the so-called "extreme" AGB stars) account for 87%-89% of the total dust input from cool evolved stars. We also estimate the dust input from hot stars and supernovae (SNe), and find that if SNe produce 10^-3 solar masses of dust each, then the total SN dust input and AGB input are roughly equivalent. We consider several scenarios of SNe dust production and destruction and find that the interstellar medium (ISM) dust can be accounted for solely by stellar sources if all SNe produce dust in the quantities seen around the dustiest examples and if most SNe explode in dense regions where much of the ISM dust is shielded from the shocks. We find that AGB stars contribute only 2.1% of the ISM dust. Without a net positive contribution from SNe to the dust budget, this suggests that dust must grow in the ISM or be formed by another unknown mechanism.

preprint2010arXiv

Is Dust Forming on the Red Giant Branch in 47 Tuc?

Using Spitzer IRAC observations from the SAGE-SMC Legacy program and archived Spitzer IRAC data, we investigate dust production in 47 Tuc, a nearby massive Galactic globular cluster. A previous study detected infrared excess, indicative of circumstellar dust, in a large population of stars in 47 Tuc, spanning the entire Red Giant Branch (RGB). We show that those results suffered from effects caused by stellar blending and imaging artifacts and that it is likely that no stars below about 1 mag from the tip of the RGB are producing dust. The only stars that appear to harbor dust are variable stars, which are also the coolest and most luminous stars in the cluster.

preprint2010arXiv

Line-profile tomography of exoplanet transits -- II. A gas-giant planet transiting a rapidly-rotating A5 star

Most of our knowledge of extrasolar planets rests on precise radial-velocity measurements, either for direct detection or for confirmation of the planetary origin of photometric transit signals. This has limited our exploration of the parameter space of exoplanet hosts to solar- and later-type, sharp-lined stars. Here we extend the realm of stars with known planetary companions to include hot, fast-rotating stars. Planet-like transits have previously been reported in the lightcurve obtained by the SuperWASP survey of the A5 star HD15082 (WASP-33; V=8.3, v sin i = 86 km/sec). Here we report further photometry and time-series spectroscopy through three separate transits, which we use to confirm the existence of a gas giant planet with an orbital period of 1.22d in orbit around HD15082. From the photometry and the properties of the planet signal travelling through the spectral line profiles during the transit we directly derive the size of the planet, the inclination and obliquity of its orbital plane, and its retrograde orbital motion relative to the spin of the star. This kind of analysis opens the way to studying the formation of planets around a whole new class of young, early-type stars, hence under different physical conditions and generally in an earlier stage of formation than in sharp-lined late-type stars. The reflex orbital motion of the star caused by the transiting planet is small, yielding an upper mass limit of 4.1 Jupiter masses on the planet. We also find evidence of a third body of sub-stellar mass in the system, which may explain the unusual orbit of the transiting planet. In HD 15082, the stellar line profiles also show evidence of non-radial pulsations, clearly distinct from the planetary transit signal. This raises the intriguing possibility that tides raised by the close-in planet may excite or amplify the pulsations in such stars.

preprint2010arXiv

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

preprint2010arXiv

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud. Point source classification I

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

preprint2009arXiv

WASP-12b: The hottest transiting planet yet discovered

We report on the discovery of WASP-12b, a new transiting extrasolar planet with $R_{\rm pl}=1.79 \pm 0.09 R_J$ and $M_{\rm pl}=1.41 \pm 0.1 M_J$. The planet and host star properties were derived from a Monte Carlo Markov Chain analysis of the transit photometry and radial velocity data. Furthermore, by comparing the stellar spectrum with theoretical spectra and stellar evolution models, we determined that the host star is a super-solar metallicity ([M/H]$=0.3^{+0.05}_{-0.15}$), late-F (T$_{\rm eff}=6300^{+200}_{-100}$ K) star which is evolving off the zero age main sequence. The planet has an equilibrium temperature of T$_{\rm eq}$=2516 K caused by its very short period orbit ($P=1.09$ days) around the hot, 12th magnitude host star. WASP-12b has the largest radius of any transiting planet yet detected. It is also the most heavily irradiated and the shortest period planet in the literature.

preprint2008arXiv

The sub-Jupiter mass transiting exoplanet WASP-11b

We report the discovery of a sub-Jupiter mass exoplanet transiting a magnitude V=11.7 host star 1SWASP J030928.54+304024.7. A simultaneous fit to the transit photometry and radial-velocity measurements yield a planet mass M_p=0.53+-0.07M_J, radius R_p=0.91^{+0.06}_{-0.03}R_J and an orbital period of 3.722465^{+0.000006}_{-0.000008} days. The host star is of spectral type K3V, with a spectral analysis yielding an effective temperature of 4800+-100K and log g=4.45+-0.2. It is amongst the smallest, least massive and lowest luminosity stars known to harbour a transiting exoplanet. WASP-11b is the third least strongly irradiated transiting exoplanet discovered to date, experiencing an incident flux F_p=1.9x10^8 erg s^{-1} cm^{-2} and having an equilibrium temperature T_eq=960+-70K.