Search

Context. Current implementations of mass loss for hot, massive stars in stellar evolution models include a sharp increase in mass loss when blue supergiants become cooler than Teff 20-22kK. This drastic mass-loss jump has been motivated by the potential presence of a so-called bistability ionisation effect, which may occur for line-driven winds in this temperature region due to recombination of important line-driving ions. Aims. We perform quantitative spectroscopy using UV (ULLYSES program) and optical (XShootU collaboration) data for 17 OB-supergiant stars in the LMC (covering the range Teff 14-32kK), deriving absolute constraints on global stellar, wind, and clumping parameters. We examine whether there are any empirical signs of a mass-loss jump in the investigated region, and we study the clumped nature of the wind. Methods. We use a combination of the model atmosphere code fastwind and the genetic algorithm code Kiwi-GA to fit synthetic spectra of a multitude of diagnostic spectral lines in the optical and UV. Results. We find no signs of any upward mass loss jump anywhere in the examined region. Standard theoretical comparison models, which include a strong bistability jump thus severely over predict the empirical mass-loss rates on the cool side of the predicted jump. Additionally, we find that on average about 40% of the total wind mass seems to reside in the diluted medium in between dense clumps. Conclusions. Our derived mass-loss rates suggest that for applications like stellar evolution one should not include a drastic bistability jump in mass loss for stars in the temperature and luminosity region investigated here. The derived high values of interclump density further suggest that the common assumption of an effectively void interclump medium (applied in the vast majority of spectroscopic studies of hot star winds) is not generally valid in this parameter regime.

Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the SMC - the lowest metallicity conditions in which multiplicity is probed to date (Z = 0.2 Zsun). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods P < 3 yr, (iii) dormant OB+BH binaries, and (iv) a legacy database of physical parameters of massive stars at low metallicity. The stars are observed with the LR02 setup of the giraffe instrument of the Very Large Telescope (3960-4570A, resolving power R=6200; typical signal-to-noise ratio S/N=70-100). This paper utilises the first 9 epochs obtained over a three-month time. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. The sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5

Context. B supergiants (BSGs) represent an important connection between the main sequence and more extreme evolutionary stages of massive stars. Additionally, lying toward the cool end of the hot star regime, determining their wind properties is crucial to constrain the evolution and feedback of massive stars as, for instance, they might manifest the bi-stability jump phenomenon. Aims. We undertake a detailed analysis of a representative sample of 18 Small Magellanic Cloud (SMC) BSGs within the ULLYSES and XShootU datasets. Our UV and optical analysis spans BSGs from B0 to B8 - covering the bi-stability jump region. We aim to evaluate their evolutionary status and verify what their wind properties say about the bi-stability jump in a low-metallicity environment. Methods. We used the CMFGEN to model the spectra and photometry (from UV to infrared) of our sample. We compare our results with different evolutionary models, with previous determinations in the literature of OB stars, and with diverging mass-loss recipes at the bi-stability jump. Additionally, we provide the first BSG models in the SMC including X-rays. Results. (i) Within a single-stellar evolution framework, the evolutionary status of early BSGs seem less clear than that of late BSGs, which agree with H-shell burning models. (ii) UV analysis shows evidence that BSGs contain X-rays in their atmospheres, for which we provide constraints. In general, we find higher X-ray luminosity (close to the standard log(L_X/L) ~ -7) for early BSGs. For cooler BSGs, lower values are preferred, log(L_X/L) ~ -8.5. (iii) The obtained mass-loss rates suggest neither a jump nor a monotonic decrease with temperature. Instead, a rather constant trend is observed, which is at odds with the increase found for Galactic BSGs. (iv) The wind velocity behavior with temperature shows a sharp drop at ~19 kK, similar to what is observed for Galactic BSGs., Comment: 33 pages (23+10)+(22 at Zenodo), 34 figures (21+13)+(21 at Zenodo), 7 tables (3+4)+(1 at Zenodo), accepted for publication

Context: There have been many attempts to identify families of diffuse interstellar bands (DIBs) with perfectly correlating band strengths. Although major efforts have been made to classify broadly based DIB families and important insights have been gained, no family has been identified with sufficient accuracy or statistical significance to prove that a series of selected DIBs originates from the same carrier. This can be attributed in part to the exclusive use of equivalent widths to establish DIB families. Aims: In a change of strategy, we search for DIBs that are highly correlated in both band strength and profile shape. This approach increases the chance of correlating DIBs being members of one family and originating from the same carrier molecule. We also search for correlations between DIB profile families and atomic interstellar lines, with the goal of further chemically constraining possible DIB carriers. Methods: We adapted the well-known method of time-series alignment to perform a spectral alignment; that is, DIB alignment. In a second step, we analysed the alignment results using a clustering analysis. This method required a statistically significant data set of DIB sight lines. The ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES) data were perfectly suited for this application. Results: We report eight DIB families with correlating strengths and profiles, as well as four previously unreported DIBs in the visual range, found using DIB alignment. All profile family members show Pearson correlation coefficients in band strength higher than 0.9. In particular, we report the 6614 - 6521 AA DIB pair, in which both DIBs show the same triple-peak substructure and an unprecedented band strength Pearson correlation coefficient of 0.9935. The presented approach opens up new perspectives that can guide the laboratory search for DIB carriers., Comment: 53 pages, 53 figures, Accepted for publication in Astronomy & Astrophysics

We present new radio continuum images and a source catalogue from the MeerKAT survey in the direction of the Small Magellanic Cloud (SMC). The observations, at a central frequency of 1.3 GHz across a bandwidth of 0.8 GHz, encompass a field of view ~7 x 7 degrees and result in images with resolution of 8 arcsec. The median broad-band Stokes I image Root Mean Squared noise value is ~11 microJy/beam. The catalogue produced from these images contains 108,330 point sources and 517 compact extended sources. We also describe a UHF (544-1088 MHz) single pointing observation. We report the detection of a new confirmed Supernova Remnant (SNR) (MCSNR J0100-7211) with an X-ray magnetar at its centre and 10 new SNR candidates. This is in addition to the detection of 21 previously confirmed SNRs and two previously noted SNR candidates. Our new SNR candidates have typical surface brightness an order of magnitude below those previously known, and on the whole they are larger. The high sensitivity of the MeerKAT survey also enabled us to detect the bright end of the SMC Planetary Nebulae (PNe) sample - point-like radio emission is associated with 38 of 102 optically known PNe, of which 19 are new detections. Lastly, we present the detection of three foreground radio stars amidst 11 circularly polarised sources, and a few examples of morphologically interesting background radio galaxies from which the radio ring galaxy ESO 029-G034 may represent a new type of radio object., Comment: 31 pages, 27 figures

We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the large velocity gradient approximation at a resolution of 3.5$''$ (~0.8 pc at the LMC distance). The physical properties of the clouds are divided into two categories: three of them near the shock front show the highest temperatures of ~50 K with densities of ~500-700 cm$^{-3}$, while other clouds slightly distant from the SNR have moderate temperatures of ~20 K with densities of ~800-1300 cm$^{-3}$. The former clouds were heated by supernova shocks, but the latter were dominantly affected by the cosmic-ray heating. These findings are consistent with the efficient production of X-ray recombining plasma in N49 due to thermal conduction between the cold clouds and hot plasma. We also find that the gas pressure is roughly constant except for the three shock-engulfed clouds inside or on the SNR shell, suggesting that almost no clouds have evaporated within the short SNR age of ~4800 yr. This result is compatible with the shock-interaction model with dense and clumpy clouds inside a low-density wind bubble., Comment: 12 pages, 5 figures, 1 table, accepted for publication in The Astrophysical Journal (ApJ)

We present new and archival Australia Telescope Compact Array and Atacama Large Millimeter/submillimeter Array data of the Small Magellanic Cloud supernova remnant 1E 0102.2-7219 at 2100, 5500, 9000, and 108000 MHz; as well as Hi data provided by the Australian Square Kilometre Array Pathfinder. The remnant shows a ring-like morphology with a mean radius of 6.2 pc. The 5500 MHz image reveals a bridge-like structure, seen for the first time in a radio image. This structure is also visible in both optical and X-ray images. In the 9000 MHz image we detect a central feature that has a flux density of 4.3 mJy but rule out a pulsar wind nebula origin, due to the lack of significant polarisation towards the central feature with an upper limit of 4 per cent. The mean fractional polarisation for 1E 0102.2-7219 is 7 +- 1 and 12 +- 2 per cent for 5500 and 9000 MHz, respectively. The spectral index for the entire remnant is -0.61 +- 0.01. We estimate the line-of-sight magnetic field strength in the direction of 1E 0102.2-7219 of ~44 microG with an equipartition field of 65 +- 5 microG. This latter model, uses the minimum energy of the sum of the magnetic field and cosmic ray electrons only. We detect an Hi cloud towards this remnant at the velocity range of ~160-180 km s-1 and a cavity-like structure at the velocity of 163.7-167.6 km s-1. We do not detect CO emission towards 1E 0102.2-7219.

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

High-spatial-resolution HI observations have led to the realisation that the nearby (within few hundreds of parsecs) Galactic atomic filamentary structures are aligned with the ambient magnetic field. Enabled by the high quality data from the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope for the Galactic ASKAP HI (GASKAP-HI) survey, we investigate the potential magnetic alignment of the $\gtrsim 10\,{\rm pc}$-scale HI filaments in the Small Magellanic Cloud (SMC). Using the Rolling Hough Transform (RHT) technique that automatically identifies filamentary structures, combined with our newly devised ray-tracing algorithm that compares the HI and starlight polarisation data, we find that the HI filaments in the northeastern end of the SMC main body ("Bar" region) and the transition area between the main body and the tidal feature ("Wing" region) appear preferentially aligned with the magnetic field traced by starlight polarisation. Meanwhile, the remaining SMC volume lacks starlight polarisation data of sufficient quality to draw any conclusions. This suggests for the first time that filamentary HI structures can be magnetically aligned across a large spatial volume ($\gtrsim\,{\rm kpc}$) outside of the Milky Way. In addition, we generate maps of the preferred orientation of HI filaments throughout the entire SMC, revealing the highly complex gaseous structures of the galaxy likely shaped by a combination of the intrinsic internal gas dynamics, tidal interactions, and star formation feedback processes. These maps can further be compared with future measurements of the magnetic structures in other regions of the SMC., Comment: 24 pages, MNRAS accepted

We present new Hubble Space Telescope WFC3/IR medium-band photometry of the compact elliptical galaxy M32, chemically resolving its thermally pulsating asymptotic giant branch stars. We find 2829 M-type stars and 57 C stars. The carbon stars are likely contaminants from M31. If carbon stars are present in M32 they are so in very low numbers. The uncorrected C/M ratio is 0.020 $\pm$ 0.003; this drops to less than 0.007 after taking into account contamination from M31. As the mean metallicity of M32 is just below solar, this low ratio of C to M stars is unlikely due to a metallicity ceiling for the formation of carbon stars. Instead, the age of the AGB population is likely to be the primary factor. The ratio of AGB to RGB stars in M32 is similar to that of the inner disc of M31 which contain stars that formed 1.5-4 Gyr ago. If the M32 population is at the older end of this age then its lack of C-stars may be consistent with a narrow mass range for carbon star formation predicted by some stellar evolution models. Applying our chemical classifications to the dusty variable stars identified with {\em Spitzer}, we find that the x-AGB candidates identified with Spitzer are predominately M-type stars. This substantially increases the lower limit to the cumulative dust-production rate in M32 to $>$ 1.20 $\times 10^{-5}$ ${\rm M}_{\odot} \, {\rm yr}^{-1}$., Comment: 10 pages, 7 figures, accepted MNRAS

Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus., Comment: Accepted for publication in A&A Letters; 9 pages, 5 figures

We present a new Australian Square Kilometre Array Pathfinder (ASKAP) sample of 14 radio Supernova Remnant (SNR) candidates in the Large Magellanic Cloud (LMC). This new sample is a significant increase to the known number of older, larger and low surface brightness LMC SNRs. We employ a multi-frequency search for each object and found possible traces of optical and occasionally X-ray emission in several of these 14 SNR candidates. One of these 14 SNR candidates (MCSNR J0522-6543) has multi-frequency properties that strongly indicate a bona fide SNR. We also investigate a sample of 20 previously suggested LMC SNR candidates and confirm the SNR nature of MCSNR J0506-6815. We detect lower surface brightness SNR candidates which were likely formed by a combination of shock waves and strong stellar winds from massive progenitors (and possibly surrounding OB stars). Some of our new SNR candidates are also found in a lower density environments in which SNe type Ia explode inside a previously excavated interstellar medium (ISM)., Comment: accepted for publication in MNRAS, this version corrects arXiv metadata only

We present proper motion (PM) measurements within the central region of the Large Magellanic Cloud (LMC) using near-infrared data from the VISTA survey of the Magellanic Cloud system (VMC). This work encompasses 18 VMC tiles covering a total sky area of $\sim$28~deg$^2$. We computed absolute stellar PMs from multi-epoch observations in the $K_s$ filter over time baselines between $\sim$12 and 47 months. Our final catalogue contains $\sim$6,322,000 likely LMC member stars with derived PMs. We employed a simple flat-rotating disc model to analyse and interpret the PM data. We found a stellar centre of rotation ($\alpha_0$ = 79.95 deg +0.22 -0.23, $\delta_0$ = -69.31 deg +0.12 -0.11) that is in agreement with that resulting from Hubble Space Telescope data. The inferred viewing angles of the LMC disc (i = 33.5 deg +1.2 -1.3, $\Theta$ = 129.8 deg +1.9 -1.9) are in good agreement with values from the literature but suggest a higher inclination of the central parts of the LMC. Our data confirm a higher rotation amplitude for the young ($\lesssim$0.5~Gyr) stars compared to the intermediate-age/old ($\gtrsim$1~Gyr) population, which can be explained by asymmetric drift. We constructed spatially resolved velocity maps of the intermediate-age/old and young populations. Intermediate-age/old stars follow elongated orbits parallel to the bar's major axis, providing first observational evidence for $x_1$ orbits within the LMC bar. In the innermost regions, the motions show more chaotic structures. Young stars show motions along a central filamentary bar structure., Comment: 18 pages, 17 figures, accepted for publication in MNRAS

We employed the {\it VISTA near-infrared $YJK_\mathrm{s}$ survey of the Magellanic System} (VMC), to analyse the $Y,\,J,\,K_\mathrm{s}$ light curves of $\delta$ Cepheid stars (DCEPs) in the Large Magellanic Cloud (LMC). Our sample consists of 4408 objects accounting for 97 per cent of the combined list of OGLE\,IV and {\it Gaia}\,DR2 DCEPs. We determined a variety of period-luminosity ($PL$) and period-Wesenheit $PW$ relationships for Fundamental (F) and First Overtone (1O) pulsators. We discovered for the first time a break in these relationships for 1O DCEPs at $P$=0.58 d. We derived relative individual distances for DCEPs in the LMC with a precision of $\sim$1 kpc, calculating the position angle of the line of nodes and inclination of the galaxy: $\theta$=145.6$\pm$1.0 deg and $i$=25.7$\pm$0.4 deg. The bar and the disc are seen under different viewing angles. We calculated the ages of the pulsators, finding two main episodes of DCEP formation lasting $\sim$40 Myr which happened 93 and 159 Myr ago. Likely as a result of its past interactions with the SMC, the LMC shows a non-planar distribution, with considerable structuring: the bar is divided into two distinct portions, the eastern and the western displaced by more than 1 kpc from each other. Similar behaviour is shown by the spiral arms. The LMC disc appears "flared" and thick, with a disc scale height of $h\sim 0.97$ kpc. This feature can be explained by strong tidal interactions with the Milky Way and/or the Small Magellanic Cloud or past merging events with now disrupted LMC satellites., Comment: 22 pages, 21 figures, accepted for publication on MNRAS. Tables 1 and 2 are available in advance of publication upon request from the first author

We report the discovery of J0624-6948, a low-surface brightness radio ring, lying between the Galactic Plane and the Large Magellanic Cloud (LMC). It was first detected at 888 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP), and with a diameter of ~196 arcsec. This source has phenomenological similarities to Odd Radio Circles (ORCs). Significant differences to the known ORCs - a flatter radio spectral index, the lack of a prominent central galaxy as a possible host, and larger apparent size - suggest that J0624-6948 may be a different type of object. We argue that the most plausible explanation for J0624-6948 is an intergalactic supernova remnant due to a star that resided in the LMC outskirts that had undergone a single-degenerate type Ia supernova, and we are seeing its remnant expand into a rarefied, intergalactic environment. We also examine if a massive star or a white dwarf binary ejected from either galaxy could be the supernova progenitor. Finally, we consider several other hypotheses for the nature of the object, including the jets of an active galactic nucleus (AGN) or the remnant of a nearby stellar super-flare., Comment: 20 pages accepted to MNRAS

We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal HI in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K (1.6 mJy/beam) per 0.98 km s$^{-1}$ spectral channel with an angular resolution of 30$''$ ($\sim$10 pc). We discuss the calibration scheme and the custom imaging pipeline that utilizes a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire $\sim$25 deg$^2$ field-of-view. We provide an overview of the data products and characterize several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power-law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high velocity cloud with previous ASKAP+Parkes HI test observations., Comment: Accepted for publication in PASA, 34 pages, 18 figures, 5 tables

Using the Australian Square Kilometre Array Pathfinder to measure 21-cm absorption spectra toward continuum background sources, we study the cool phase of the neutral atomic gas in the far outer disk, and in the inner Galaxy near the end of the Galactic bar at longitude 340 degrees. In the inner Galaxy the cool atomic gas has a smaller scale height than in the solar neighborhood, similar to the molecular gas and the superthin stellar population in the bar. In the outer Galaxy the cool atomic gas is mixed with the warm, neutral medium, with the cool fraction staying roughly constant with Galactic radius. The mean spin temperature, i.e. the ratio of the emission brightness temperature to the absorption, is roughly constant for velocities corresponding to Galactic radius greater than about twice the solar circle radius. The ratio has a value of about 300 K, but this does not correspond to a physical temperature in the gas. If the gas causing the absorption has kinetic temperature of about 100 K, as in the solar neighborhood, then the value 300 K indicates that the fraction of the gas mass in this phase is one-third of the total HI mass., Comment: 29 pages, 12 figures, submitted to Astrophysical Journal

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 $\mu$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., Comment: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

We present a comprehensive multi-frequency catalogue of radio sources behind the Large Magellanic Cloud between 0.2 and 20 GHz, gathered from a combination of new and legacy radio continuum surveys. This catalogue covers an area of $\sim$144~deg$^2$ at angular resolutions from 45 arcsec to $\sim$3 arcmin. We find 6434 discrete radio sources in total, of which 3789 are detected at two or more radio frequencies. We estimate the median spectral index ($\alpha$; where $S_{v}\sim\nu^\alpha$) of $\alpha = -0.89 $ and mean of $-0.88 \pm 0.48$ for 3636 sources detected exclusively at two frequencies (0.843 and 1.384 GHz) with similar resolution (FWHM $\sim$40-45 arcsec). The large frequency range of the surveys makes it an effective tool to investigate Gigahertz Peak Spectrum (GPS), Compact Steep Spectrum (CSS) and Infrared Faint Radio sources populations within our sample. We find 10 GPS candidates with peak frequencies near 5 GHz, from which we estimate their linear size. 1866 sources from our catalogue are (CSS) candidates with $\alpha <-0.8$. We found six candidates for High Frequency Peaker (HFP) sources, whose radio fluxes peak above 5 GHz and no sources with unconstrained peaks and $\alpha~>0.5$. We found optical counterparts for 343 of the radio continuum sources, of which 128have a redshift measurement. Finally, we investigate the population of 123 Infrared Faint Radio Sources (IFRSs) found in this study., Comment: Accepted for publication in MNRAS

We present results from an analysis of $\sim$ 29,000 RR Lyrae stars located in the Large Magellanic Cloud (LMC). For these objects, near-infrared time-series photometry from the VISTA survey of the Magellanic Clouds system (VMC) and optical data from the OGLE (Optical Gravitational Lensing Experiment) IV survey and the Gaia Data Release 2 catalogue of confirmed RR Lyrae stars were exploited. Using VMC and OGLE IV magnitudes we derived period-luminosity (PL), period-luminosity-metallicity (PLZ), period-Wesenheit (PW) and period-Wesenheit-metallicity (PWZ) relations in all available bands. More that ~7,000 RR Lyrae were discarded from the analysis because they appear to be overluminous with respect to the PL relations. The $PL_{K_{\mathrm{s}}}$ relation was used to derive individual distance to $\sim 22,000$ RR Lyrae stars, and study the three-dimensional structure of the LMC. The distribution of the LMC RR Lyrae stars is ellipsoidal with the three axis $S_1$=6.5 kpc, $S_2$=4.6 kpc and $S_3$=3.7 kpc, inclination i=$22\pm4^{\circ }$ relative to the plane of the sky and position angle of the line of nodes $\theta=167\pm7^{\circ }$ (measured from north to east). The north-eastern part of the ellipsoid is closer to us and no particular associated substructures are detected as well as any metallicity gradient., Comment: accepted for publication in MNRAS

We used data from the near-infrared VISTA survey of the Magellanic Cloud system (VMC) to measure proper motions (PMs) of stars within the Small Magellanic Cloud (SMC). The data analysed in this study comprise 26 VMC tiles, covering a total contiguous area on the sky of ~40 deg$^2$. Using multi-epoch observations in the Ks band over time baselines between 13 and 38 months, we calculated absolute PMs with respect to ~130,000 background galaxies. We selected a sample of ~2,160,000 likely SMC member stars to model the centre-of-mass motion of the galaxy. The results found for three different choices of the SMC centre are in good agreement with recent space-based measurements. Using the systemic motion of the SMC, we constructed spatially resolved residual PM maps and analysed for the first time the internal kinematics of the intermediate-age/old and young stellar populations separately. We found outward motions that point either towards a stretching of the galaxy or stripping of its outer regions. Stellar motions towards the North might be related to the "Counter Bridge" behind the SMC. The young populations show larger PMs in the region of the SMC Wing, towards the young Magellanic Bridge. In the older populations, we further detected a coordinated motion of stars away from the SMC in the direction of the Old Bridge as well as a stream towards the SMC., Comment: 20 pages, 20 figures, accepted for publication in MNRAS

We present observations of a region of the Galactic plane taken during the Early Science Program of the Australian Square Kilometre Array Pathfinder (ASKAP). In this context, we observed the SCORPIO field at 912 MHz with an uncompleted array consisting of 15 commissioned antennas. The resulting map covers a square region of ~40 deg^2, centred on (l, b)=(343.5{\deg}, 0.75{\deg}), with a synthesized beam of 24"x21" and a background rms noise of 150-200 {\mu}Jy/beam, increasing to 500-600 {\mu}Jy/beam close to the Galactic plane. A total of 3963 radio sources were detected and characterized in the field using the CAESAR source finder. We obtained differential source counts in agreement with previously published data after correction for source extraction and characterization uncertainties, estimated from simulated data. The ASKAP positional and flux density scale accuracy were also investigated through comparison with previous surveys (MGPS, NVSS) and additional observations of the SCORPIO field, carried out with ATCA at 2.1 GHz and 10" spatial resolution. These allowed us to obtain a measurement of the spectral index for a subset of the catalogued sources and an estimated fraction of (at least) 8% of resolved sources in the reported catalogue. We cross-matched our catalogued sources with different astronomical databases to search for possible counterparts, finding ~150 associations to known Galactic objects. Finally, we explored a multiparametric approach for classifying previously unreported Galactic sources based on their radio-infrared colors., Comment: 31 pages, 15 figures, 15 tables

Galaxy interactions distort the distribution of baryonic matter and can affect star formation. The nearby Magellanic Clouds are a prime example of an ongoing galaxy interaction process. Here we use the intermediate-age ($\sim1$-$10$ Gyr) red clump stars to map the three-dimensional structure of the Small Magellanic Cloud (SMC) and interpret it within the context of its history of interaction with the Large Magellanic Cloud (LMC) and the Milky Way. Red clump stars are selected from near-infrared colour-magnitude diagrams based on data from the VISTA survey of the Magellanic Clouds. Interstellar reddening is measured and removed, and the corrected brightness is converted to a distance, on a star-by-star basis. A flat plane fitted to the spatial distribution of red clump stars has an inclination $i=35\deg$-$48\deg$ and position angle PA$=170\deg$-$186\deg$. However, significant deviations from this plane are seen, especially in the periphery and on the eastern side of the SMC. In the latter part, two distinct populations are present, separated in distance by as much as 10 kpc. Distant red clump stars are seen in the North of the SMC, and possibly also in the far West; these might be associated with the predicted `Counter-Bridge'. We also present a dust reddening map, which shows that dust generally traces stellar mass. The structure of the intermediate-age stellar component of the SMC bears the imprints of strong interaction with the LMC a few Gyr ago, which cannot be purely tidal but must have involved ram pressure stripping., Comment: 16 pages, 11 figures; accepted for publication in MNRAS

(abridged)Variability is a key property of stars on the asymptotic giant branch (AGB). Their pulsation period is related to the luminosity and mass-loss rate of the star. The long-period variables (LPVs) and Mira variables are the most prominent of all types of variability of evolved stars. The reddest, most obscured AGB stars are too faint in the optical and have eluded large variability surveys. Selection criteria are derived based on colour-colour and colour-magnitude diagrams from the combination of VISTA Magellanic Cloud (VMC) survey, Spitzer IRAC and AllWISE data. After eliminating LPVs with known periods shorter than 450 days, a sample of 1299 candidate obscured AGB stars is selected. K-band LCs are constructed combining the epoch photometry available in the VMC survey with literature data, analysed for variability and fitted with a single period sine curve to derive mean magnitudes, amplitudes and periods. A subset of 254 stars are either new variables, or known variables where the period we find is better determined than the literature value, or longer than 1000 days. The spectral energy distributions (SEDs) of these stars are fitted to a large number of templates. For this purpose the SEDs and Spitzer IRS spectra of some non-AGB stars are also fitted to have templates of the most likely contaminants in the sample. A sample of 217 likely LPVs is found. Thirty-four stars have periods longer than 1000 days although some of them have alternative shorter periods. Mass-loss rates and luminosities are estimated from the template fitting. Period-luminosity relations are presented for C- and O-rich Miras that appear to be extensions of relations derived in the literature for shorter periods. The fit for the C-stars is particularly well defined (with 182 objects) and reads Mbol = (-2.27 \pm 0.20) log P + (1.45 \pm 0.54)mag with an rms of 0.41 mag., Comment: A&A accepted

We present Herschel Space Observatory Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE FTS) spectroscopy of a sample of twenty massive Young Stellar Objects (YSOs) in the Large and Small Magellanic Clouds (LMC and SMC). We analyse the brightest far infrared (far-IR) emission lines, that diagnose the conditions of the heated gas in the YSO envelope and pinpoint their physical origin.We compare the properties of massive Magellanic and Galactic YSOs.We find that [OI] and [CII] emission, that originates from the photodissociation region associated with the YSOs, is enhanced with respect to the dust continuum in the Magellanic sample. Furthermore the photoelectric heating efficiency is systematically higher for Magellanic YSOs, consistent with reduced grain charge in low metallicity environments. The observed CO emission is likely due to multiple shock components. The gas temperatures, derived from the analysis of CO rotational diagrams, are similar to Galactic estimates. This suggests a common origin to the observed CO excitation, from low-luminosity to massive YSOs, both in the Galaxy and the Magellanic Clouds. Bright far-IR line emission provides a mechanism to cool the YSO environment. We find that, even though [OI], CO and [CII] are the main line coolants, there is an indication that CO becomes less important at low metallicity, especially for the SMC sources. This is consistent with a reduction in CO abundance in environments where the dust is warmer due to reduced ultraviolet-shielding. Weak H$_2$O and OH emission is detected, consistent with a modest role in the energy balance of wider massive YSO environments., Comment: 38 pages; accepted for publication by MNRAS; full integrated version; journal version will include appendices on-line only

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies., Comment: 18 pages, 23 figues, 9 tables, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time., Comment: A white paper submitted for the Voyage 2050 long-term plan in the ESA Science Programme. 21 pages, 1 table, 6 figures. arXiv admin note: substantial text overlap with arXiv:1903.05235

The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of $\log (M_{\rm dyn}/M_{\odot})=$3.8$\pm$0.3 for Hodge 301, and an upper limit of $\log (M_{\rm dyn}/M_{\odot})$<3.1$\pm$0.8 for SL 639, assuming Virial equilibrium. Analysis of the multi-epoch data reveals one RV-variable, potential binary candidate (VFTS744), which is likely a semi-regular variable asymptotic giant branch star. We estimate an upper limit on the observed binary fraction for our sample of 0.3, where we are sensitive to maximum periods for individual objects in the range of 1 to 10 000 days and mass-ratios above 0.3 depending on the data quality. From simulations of the RV measurements from binary systems given the current data we conclude that systems within the parameter range q>0.3, $\log$P[days]<3.5, would be detected by our variability criteria, at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters where orbital periods are uniformly distributed in the 3.3<$\log$P[days]<4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. [Abridged], Comment: accepted in A&A

We estimate physical parameters for the late-type massive stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud (LMC). The observational sample comprises 20 candidate red supergiants (RSGs) which are the reddest (($B-V$) $>$ 1 mag) and brightest ($V$ $<$ 16 mag) objects in the VFTS. We use optical and near-IR photometry to estimate their temperatures and luminosities, and introduce the luminosity-age diagram to estimate their ages. We derive physical parameters for our targets, including temperatures from a new calibration of $(J-K_{\rm s})_{0}$ colour for luminous cool stars in the LMC, luminosities from their $J$-band magnitudes (thence radii), and ages from comparisons with state-of-the-art evolutionary models. We show that interstellar extinction is a significant factor for our targets, highlighting the need to take it into account in analysis of the physical parameters of RSGs. We find that some of the candidate RSGs could be massive AGB stars. The apparent ages of the RSGs in the Hodge 301 and SL 639 clusters show a significant spread (12-24 Myr). We also apply our approach to the RSG population of the relatively nearby NGC 2100 cluster, finding a similarly large spread. We argue that the effects of mass-transfer in binaries may lead to more massive and luminous RSGs (which we call `red stragglers') than expected from single-star evolution, and that the true cluster ages correspond to the upper limit of the estimated RSG ages. In this way, the RSGs can serve as a new and potentially reliable age tracer in young star clusters. The corresponding analysis yields ages of 24$^{+5}_{-3}$ Myr for Hodge 301, 22$^{+6}_{-5}$ Myr for SL 639, and 23$^{+4}_{-2}$ Myr for NGC 2100., Comment: 13 pages, 8 figures, accepted for publication in Astronomy and Astrophysics journal

In this paper JK_s data from the VISTA Magellanic Cloud (VMC) survey are used to investigate the tip of the red giant branch (TRGB) as a distance indicator. A linear fit to recent theoretical models is used which reads M_{K_s} = -4.196 -2.013 (J-K_s), valid in the colour range 0.75 < (J-K_s) < 1.3 mag and in the 2MASS system. The observed TRGB is found based on a classical first-order and a second-order derivative filter applied to the binned luminosity function using the "sharpened" magnitude that takes the colour term into account. Extensive simulations are carried out to investigate any biases and errors in the derived distance modulus (DM). Based on these simulations criteria are established related to the number of stars per bin in the 0.5 magnitude range below the TRGB and related to the significance with which the peak in the filter response curve is determined such that the derived distances are unbiased. The DMs based on the second-order derivative filter are found to be more stable and are therefore adopted, although this requires twice as many stars per bin. The TRGB method is applied to specific lines-of-sight where independent distance estimates exist, based on detached eclipsing binaries in the LMC and SMC, classical Cepheids in the LMC, RR Lyrae stars in the SMC, and fields in the SMC where the star formation history (together with reddening and distance) has been derived from deep VMC data. The analysis shows that the theoretical calibration is consistent with the data, that the systematic error on the DM is approximately 0.045 mag, and that random errors of 0.015 mag are achievable. Reddening is an important element in deriving the distance: we find mean DMs ranging from 18.92 (for a typical E(B-V) of 0.15 mag) to 19.07 mag (E(B-V) about 0.04) for the SMC, and ranging from 18.48 (E(B-V) about 0.12 mag) to 18.57 mag (E(B-V) about 0.05) for the LMC., Comment: A&A accepted

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus., Comment: Authors' version of a response to a technical comment published in Science; 8 pages, 1 figure

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged], Comment: 20 pages (incl. appendix), 15 figures, 2 tables; accepted for publication in A&A

N49 (LHA 120-N49) is a bright X-ray supernova remnant (SNR) in the Large Magellanic Cloud. We present new $^{12}$CO($J$ = 1-0, 3-2), HI, and 1.4 GHz radio-continuum observations of the SNR N49 using Mopra, ASTE, ALMA, and ATCA. We have newly identified three HI clouds using ATCA with an angular resolution of ~20": one associated with the SNR and the others located in front of the SNR. Both the CO and HI clouds in the velocity range from 280-291 km s$^{-1}$ are spatially correlated with both the soft X-rays (0.2-1.2 keV) and the hard X-rays (2.0-7.0 keV) of N49 on a ~10 pc scale. CO 3-2/1-0 intensity ratios indicate higher values of the CO cloud toward the SNR shell with an angular resolution of ~45", and thus a strong interaction was suggested. Using the ALMA, we have spatially resolved CO clumps embedded within or along the southeastern rim of N49 with an angular resolution of ~3''. Three of the CO clumps are rim-brightened on a 0.7-2 pc scale in both hard X-rays and the radio continuum$:$ this provides further evidence for dynamical interactions between the CO clumps and the SNR shock wave. The enhancement of the radio synchrotron radiation can be understood in terms of magnetic-field amplification around the CO clumps via a shock-cloud interaction. We also present a possible scenario in which the recombining plasma that dominates the hard X-rays from N49 was formed via thermal conduction between the SNR shock waves and the cold$/$dense molecular clumps., Comment: 28 pages, 12 figures, 1 table, accepted for publication in The Astrophysical Journal (ApJ)

We present the first spatially resolved map of stellar proper motions within the central ($\sim$3.1 $\times$ 2.4 kpc) regions of the Small Magellanic Cloud (SMC). The data used for this study encompasses four tiles from the ongoing near-infrared VISTA survey of the Magellanic Clouds system and covers a total contiguous area on the sky of $\sim$6.81 deg$^2$. Proper motions have been calculated independently in two dimensions from the spatial offsets in the $K_s$ filter over time baselines between 22 and 27 months. The reflex motions of $\sim$33~000 background galaxies are used to calibrate the stellar motions to an absolute scale. The resulting catalog is composed of more than 690 000 stars which have been selected based on their position in the $(J-K_s, K_s)$ color-magnitude diagram. For the median absolute proper motion of the SMC, we find ($\mu_{\alpha}\mathrm{cos}(\delta)$, $\mu_{\delta}$) = (1.087 $\pm$ 0.192 (sys.) $\pm$ 0.003 (stat.), $-$1.187 $\pm$ 0.008 (sys.) $\pm$ 0.003 (stat.)) mas yr$^{-1}$, consistent with previous studies. Mapping the proper motions as a function of position within the SMC reveals a non uniform velocity pattern indicative of a tidal feature behind the main body of the SMC and a flow of stars in the South-East moving predominantly along the line-of-sight., Comment: 4 pages, 2 figures, accepted for publication in A&A

We use deep multi-epoch PSF photometry taken with the Visible and Infrared Survey Telescope for Astronomy (VISTA) to measure and analyze the proper motions of stars within the Galactic globular cluster 47 Tucanae (47 Tuc, NGC 104). The observations are part of the ongoing near-infrared VISTA survey of the Magellanic Cloud system (VMC). The data analyzed here corresponds to one VMC tile and covers a total sky area of 1.77 deg^2. Absolute proper motions with respect to ~9070 background galaxies are calculated from a linear regression model applied to the stellar positions in 11 epochs in the Ks filter. The data extend over a total time baseline of ~17 months. We found a median proper motion of the stars within 47 Tuc of (mu_a cos(d), mu_d) = (+5.89 +/- 0.02(stat) +/- 0.13(sys), -2.14 +/- 0.02(stat) +/- 0.08(sys)) mas/yr, based on the measurements of ~35000 individual sources between 5' and 42' from the cluster center. We compared our result to the proper motions from the newest US Naval Observatory CCD Astrograph Catalog (UCAC5) which includes data from the Gaia data release 1. Selecting cluster members (~2700 stars) we found a median proper motion of (mu_a cos(d), mu_d) = (+5.30 +/- 0.03(stat) +/- 0.70(sys), -2.70 +/- 0.03(stat) +/- 0.70(sys)) mas/yr. The values derived from the VMC data are consistent with the UCAC5 result, and are close to literature measurements obtained using HST. We combined our results with literature radial velocity measurements and reconstructed the orbit of 47 Tuc, finding that the cluster is on a low-ellipticity orbit and is confined within the inner ~7.5 kpc of the Galaxy. We show that the use of an increased time baseline combined with PSF-determined stellar centroids in crowded regions significantly improves the accuracy of the method. In future works, we will apply this method to more VMC tiles to study in detail the kinematics of the Magellanic Clouds., Comment: 16 pages, 15 figures, accepted for publication in A&A

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35., Comment: Authors' version of paper published in Science at http://science.sciencemag.org/content/359/6371/69 ; 15 pages main text (3 figures), 47 pages supplementary materials (10 figures, 3 tables)

The survey for DUST in Nearby Galaxies with Spitzer (DUSTiNGS) identified several candidate Asymptotic Giant Branch (AGB) stars in nearby dwarf galaxies and showed that dust can form even in very metal-poor systems (Z ~ 0.008 $Z_\odot$). Here, we present a follow-up survey with WFC3/IR on the Hubble Space Telescope (HST), using filters that are capable of distinguishing carbon-rich (C-type) stars from oxygen-rich (M-type) stars: F127M, F139M, and F153M. We include six star-forming DUSTiNGS galaxies (NGC 147, IC 10, Pegasus dIrr, Sextans B, Sextans A, and Sag DIG), all more metal-poor than the Magellanic Clouds and spanning 1 dex in metallicity. We double the number of dusty AGB stars known in these galaxies and find that most are carbon rich. We also find 26 dusty M-type stars, mostly in IC 10. Given the large dust excess and tight spatial distribution of these M-type stars, they are most likely on the upper end of the AGB mass range (stars undergoing Hot Bottom Burning). Theoretical models do not predict significant dust production in metal-poor M-type stars, but we see evidence for dust excess around M-type stars even in the most metal-poor galaxies in our sample (12+log(O/H) = 7.26-7.50). The low metallicities and inferred high stellar masses (up to ~10 $M_\odot$) suggest that AGB stars can produce dust very early in the evolution of galaxies (~30 Myr after they form), and may contribute significantly to the dust reservoirs seen in high-redshift galaxies., Comment: 18 pages, 9 figures, accepted for publication in ApJ

We present results from the analysis of 2997 fundamental mode RR Lyrae variables located in the Small Magellanic Cloud (SMC). For these objects near-infrared time-series photometry from the VISTA survey of the Magellanic Clouds system (VMC) and visual light curves from the OGLE IV survey are available. In this study the multi-epoch $K_{\rm s}$-band VMC photometry was used for the first time to derive intensity-averaged magnitudes of the SMC RR Lyrae stars. We determined individual distances to the RR Lyrae stars from the near-infrared period-absolute magnitude-metallicity ($PM_{K_{\rm s}}Z$) relation, which has a number of advantages in comparison with the visual absolute magnitude-metallicity ($M_{V}-{\rm [Fe/H]}$) relation, such as a smaller dependence of the luminosity on interstellar extinction, evolutionary effects and metallicity. The distances we have obtained were used to study the three-dimensional structure of the SMC. The distribution of the SMC RR Lyrae stars is found to be ellipsoidal. The actual line-of-sight depth of the SMC is in the range from 1 to 10 kpc, with an average depth of 4.3 $\pm$ 1.0 kpc. We found that RR Lyrae stars in the eastern part of the SMC are affected by interactions of the Magellanic Clouds. However, we do not see a clear bimodality in the distribution of RR Lyrae stars as observed for red clump (RC) stars., Comment: 18 pages, 15 figures, 3 tables. Accepted for publication in MNRAS

Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained., Comment: ApJL accepted, 9 pages, 5 figures

The Infrared Spectrograph (IRS) on the {\em Spitzer Space Telescope} observed nearly 800 point sources in the Large Magellanic Cloud (LMC), taking over 1,000 spectra. 197 of these targets were observed as part of the Sage-Spec Spitzer Legacy program; the remainder are from a variety of different calibration, guaranteed time and open time projects. We classify these point sources into types according to their infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, using a decision-tree classification method. We then refine the classification using supplementary information from the astrophysical literature. We find that our IRS sample is comprised substantially of YSO and H\,{\sc ii} regions, post-Main Sequence low-mass stars: (post-)AGB stars and planetary nebulae and massive stars including several rare evolutionary types. Two supernova remnants, a nova and several background galaxies were also observed. We use these classifications to improve our understanding of the stellar populations in the Large Magellanic Cloud, study the composition and characteristics of dust species in a variety of LMC objects, and to verify the photometric classification methods used by mid-IR surveys. We discover that some widely-used catalogues of objects contain considerable contamination and others are missing sources in our sample., Comment: Accepted for publication in MNRAS

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae., Comment: Accepted by MNRAS, 16 pages

There is evidence to suggest that clathrate hydrates have a significant effect on the surface geology of icy bodies in the Solar System. However the aqueous environments believed to be present on these bodies are likely to be saline rather than pure water. Laboratory work to underpin the properties of clathrates in such environments is lacking. We fill this gap by carrying out a laboratory investigation of the physical properties of CO$_2$ clathrates produced in weak aqueous solutions of MgSO$_4$. We use synchrotron X-ray powder diffraction to investigate clathrates formed at high CO$_2$ pressure in ice that has formed from aqueous solutions of MgSO$_4$. We measure the thermal expansion, density and dissociation properties of the clathrates under temperature conditions similar to those on icy Solar System bodies. We find that the sulphate solution inhibits the formation of clathrates by lowering their dissociation temperatures. Hysteresis is found in the thermal expansion coefficients as clathrates are cooled and heated; we attribute this to the presence of the salt in solution. The density derived from X-ray powder diffraction is temperature and pressure dependent. When comparing the density of CO$_2$ clathrates to that of the solution in which they formed, we conclude that they sink in the oceans in which they form. We also find that the polymorph of ice present at low temperatures is Ih rather than Ic, which we attribute to the presence of the MgSO$_4$. We 1) conclude that the clathrate density has implications for their behaviour in satellite oceans as their sinking and floating capabilities are temperature and pressure dependent, 2) conclude that the presence of MgSO$_4$ inhibits the formation of clathrates and in some cases may even affect their structure and 3) report the dominance of Ih throughout the experimental procedure despite Ic being the stable phase at low temperature., Comment: Accepted for publication by Astronomy and Astrophysics on 25 January 2017. Abstract truncated

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., Comment: Accepted for publication by A&A

We present $K$-band integral field spectroscopic observations towards 17 massive young stellar objects (YSOs) in the low metallicity Small Magellanic Cloud (SMC) and two YSO candidates in the compact H ii regions N81 and N88 A (also in the SMC). These sources, originally identified using $Spitzer$ photometry and/or spectroscopy, have been resolved into 29 $K$-band continuum sources. By comparing Br$\gamma$ emission luminosities with those presented for a Galactic sample of massive YSOs, we find tentative evidence for increased accretion rates in the SMC. Around half of our targets exhibit emission line (Br$\gamma$, He i and H$_2$) morphologies which extend significantly beyond the continuum source and we have mapped both the emission morphologies and the radial velocity fields. This analysis also reveals evidence for the existence of ionized low density regions in the centre outflows from massive YSOs. Additionally we present an analysis of optical spectra towards a similar sample of massive YSOs in the SMC, revealing that the optical emission is photo-excited and originates near the outer edges of molecular clouds, and is therefore consistent with a high mean-free path of UV photons in the interstellar medium (ISM) of the SMC. Finally, we discuss the sample of YSOs in an evolutionary context incorporating the results of previous infrared and radio observations, as well as the near-infrared and optical observations presented in this work. Our spectroscopic analysis in both the $K$-band and the optical regimes, combined with previously obtained infrared and radio data, exposes differences between properties of massive YSOs in our own Galaxy and the SMC, including tracers of accretion, discs and YSO--ISM interactions., Comment: 42 pages, 22 figures, Accepted for publication in MNRAS 2016 September 16

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

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., Comment: Accepted for publication by ApJ, 20 pages, 20 figures. Full versions of the tables are available at: http://isc.astro.cornell.edu/~sloan/library/2016/mcc

We present observations of CK Vul obtained with the Spitzer Space Telescope. The infrared spectrum reveals a warm dust continuum with nebular, molecular hydrogen and HCN lines superimposed, together with the "Unidentified Infrared" (UIR) features. The nebular lines are consistent with emission by a low density gas. We conclude that the Spitzer data, combined with other information, are incompatible with CK Vul being a classical nova remnant in "hibernation" after the event of 1670, a "Very Late Thermal Pulse", a "Luminous Red Variable" such as V838 Mon, or a "Diffusion-induced nova". The true nature of CK Vul remains a mystery., Comment: MNRAS letters, in press. Paper consists of 6 pages, 4 figures, 2 tables

The \textit{Spitzer} SAGE survey has allowed the identification and analysis of significant samples of Young Stellar Object (YSO) candidates in the Large Magellanic Cloud (LMC). However the angular resolution of \textit{Spitzer} is relatively poor meaning that at the distance of the LMC, it is likely that many of the \textit{Spitzer} YSO candidates in fact contain multiple components. We present high resolution \textit{K}-band integral field spectroscopic observations of the three most prominent massive YSO candidates in the N113 H\,{\sc ii} region using VLT/SINFONI. We have identified six \textit{K}-band continuum sources within the three \textit{Spitzer} sources and we have mapped the morphology and velocity fields of extended line emission around these sources. Br$\gamma$, He\,{\sc i} and H$_2$ emission is found at the position of all six \textit{K}-band sources; we discuss whether the emission is associated with the continuum sources or whether it is ambient emission. H$_2$ emission appears to be mostly ambient emission and no evidence of CO emission arising in the discs of YSOs has been found. We have mapped the centroid velocities of extended Br$\gamma$ emission and He {\sc i} emission and found evidence of two expanding compact H\,{\sc ii} regions. One source shows compact and strong H$_2$ emission suggestive of a molecular outflow. The diversity of spectroscopic properties observed is interpreted in the context of a range of evolutionary stages associated with massive star formation., Comment: 20 pages, 12 figures, Accepted by MNRAS 2015 October 19

We present spectral classifications from optical spectroscopy of 263 massive stars in the north-eastern region of the Large Magellanic Cloud. The observed two-degree field includes the massive 30 Doradus star-forming region, the environs of SN1987A, and a number of star-forming complexes to the south of 30 Dor. These are the first classifications for the majority (203) of the stars and include eleven double-lined spectroscopic binaries. The sample also includes the first examples of early OC-type spectra (AAOmega 30 Dor 248 and 280), distinguished by the weakness of their nitrogen spectra and by C IV 4658 emission. We propose that these stars have relatively unprocessed CNO abundances compared to morphologically normal O-type stars, indicative of an earlier evolutionary phase. From analysis of observations obtained on two consecutive nights, we present radial-velocity estimates for 233 stars, finding one apparent single-lined binary and nine (>3sigma) outliers compared to the systemic velocity; the latter objects could be runaway stars or large-amplitude binary systems and further spectroscopy is required to investigate their nature., Comment: Accepted by A&A