Betelgeuse is the nearest Red Supergiant star and it underwent an unusually deep minimum at optical wavelengths during its most recent pulsation cycle. We present submillimetre observations taken by the James Clerk Maxwell Telescope and Atacama Pathfinder Experiment over a time span of 13 years including the optical dimming. We find that Betelgeuse has also dimmed by \sim20\% at these longer wavelengths during this optical minimum. Using radiative-transfer models, we show that this is likely due to changes in the photosphere (luminosity) of the star as opposed to the surrounding dust as was previously suggested in the literature.

The search for extrasolar planets has driven rapid advances in instrumentation, resulting in cameras such as SPHERE at the VLT, GPI at Gemini South and SCExAO at Subaru, capable of achieving very high contrast ($\sim10^{6}$) around bright stars with small inner working angles ($\sim 0\farcs{1}$). The optimal exploitation of data from these instruments depends on the availability of easy-to-use software to process and analyse their data products. We present a pure-python pipeline, {\sc precision}, which provides fast, memory-efficient reduction of data from the SPHERE/IRDIS near-infrared imager, and can be readily extended to other instruments. We apply {\sc precision} to observations of the extreme red supergiant VX~Sgr, the inner outflow of which is revealed to host complex, asymmetric structure in the near-IR. In addition, optical polarimetric imaging reveals clear extended polarised emission on $\sim0.5^{\prime\prime}$ scales which varies significantly with azimuth, confirming the asymmetry. While not conclusive, this could suggest that the ejecta are confined to a disc or torus, which we are viewing nearly face on, although other non-spherical or clumpy configurations remain possible. VX~Sgr has no known companions, making such a geometry difficult to explain, as there is no obvious source of angular momentum in the system.

The timescales on which astronomical dust grows remain poorly understood, with important consequences for our understanding of processes like circumstellar disk evolution and planet formation.A number of post-asymptotic giant branch stars are found to host optically thick, dust- and gas-rich circumstellar discs in Keplerian orbits. These discs exhibit evidence of dust evolution, similar to protoplanetary discs; however since post-AGB discs have substantially shorter lifetimes than protoplanetary discs they may provide new insights on the grain-growth process. We examine a sample of post-AGB stars with discs to determine the FIR and sub-mm spectral index by homogeneously fitting a sample of data from \textit{Herschel}, the SMA and the literature. We find that grain growth to at least hundreds of micrometres is ubiquitous in these systems, and that the distribution of spectral indices is more similar to that of protoplanetary discs than debris discs. No correlation is found with the mid-infrared colours of the discs, implying that grain growth occurs independently of the disc structure in post-AGB discs. We infer that grain growth to $\sim$mm sizes must occur on timescales $<<10^{5}$ yr, perhaps by orders of magnitude, as the lifetimes of these discs are expected to be $\lesssim10^{5}$~yr and all objects have converged to the same state. This growth timescale is short compared to the results of models for protoplanetary discs including fragmentation, and may provide new constraints on the physics of grain growth.

This white paper discusses recent progress in the field of evolved stars, primarily highlighting the contributions of the James Clerk Maxwell Telescope. It discusses the ongoing project, the \emph{Nearby Evolved Stars Survey} (NESS), and the potential to build upon NESS in the next decade. It then outlines a number of science cases which may become feasible with the proposed 850\,$\mu$m camera which is due to become available at the JCMT in late 2022. These include mapping the extended envelopes of evolved stars, including in polarisation, and time-domain monitoring of their variations. The improved sensitivity of the proposed instrument will facilitate statistical studies that put the morphology, polarisation properties and sub-mm variability in perspective with a relatively modest commitment of time that would be impossible with current instrumentation. We also consider the role that could be played by other continuum wavelengths, heterodyne instruments or other facilities in contributing towards these objectives.

Magnetic fields are ubiquitous in our Universe, but remain poorly understood in many branches of astrophysics. A key tool for inferring astrophysical magnetic field properties is dust emission polarimetry. The James Clerk Maxwell Telescope (JCMT) is planning a new 850$\mu$m camera consisting of an array of 7272 paired Microwave Kinetic Inductance Detectors (MKIDs), which will inherently acquire linear polarization information. The camera will allow wide-area polarization mapping of dust emission at 14$^{\prime\prime}$-resolution, allowing magnetic field properties to be studied in a wide range of environments, including all stages of the star formation process, Asymptotic Giant Branch stellar envelopes and planetary nebula, external galaxies including starburst galaxies and analogues for the Milky Way, and the environments of active galactic nuclei (AGN). Time domain studies of AGN and protostellar polarization variability will also become practicable. Studies of the polarization properties of the interstellar medium will also allow detailed investigation of dust grain properties and physics. These investigations would benefit from a potential future upgrade adding 450$\mu$m capability to the camera, which would allow inference of spectral indices for polarized dust emission in a range of environments. The enhanced mapping speed and polarization capabilities of the new camera will transform the JCMT into a true submillimetre polarization survey instrument, offering the potential to revolutionize our understanding of magnetic fields in the cold Universe.

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

This white paper gives a brief summary of the time domain science that has been performed with the JCMT in recent years and highlights the opportunities for continuing work in this field over the next decade. The main focus of this document is the JCMT Transient Survey, a large program initiated in 2015 to measure the frequency and amplitude of variability events associated with protostars in nearby star-forming regions. After summarising the major accomplishments so far, an outline is given for extensions to the current survey, featuring a discussion on what will be possible with the new 850 micron camera that is expected to be installed in late 2022. We also discuss possible applications of submillimetre monitoring to active galactic nuclei, X-ray binaries, asymptotic giant branch stars, and flare stars.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago.

We present the sub-mm variability of two of the most well studied AGB stars, IRC+10216 and $o$ Ceti. The data are obtained at $450~\micron$ and $850~\micron$ as part of pointing calibration observations for the James Clerk Maxwell Telescope's SCUBA-2 instrument over a span of 7 years. The periods are derived using non-parametric methods, \texttt{Gatspy Supersmoother} and \texttt{P4J} in order not to assume an underlying shape to the periodicity. These were compared to two Lomb-Scargle parametric methods. We find that for both sources and wavelengths the periods derived from all methods are consistent within $1\sigma$. The $850~\micron$ phase folded light curves of IRC+10216 show a time lag of $\sim 540$ days compared to its optical counterpart. We explore the origins of the sub-mm variability and the phase lag using radiative transfer models. Combining the modelling with findings in the literature, we find that the sub-mm emission and phase lag can be partially attributed to the dust formation/destruction cycle. A second, unknown mechanism must be invoked; we defer an investigation of the origin and nature of this mechanism to a future work.

We present the SCUBA-2 850 ${\mu}m$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 ${\mu}m$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data. We measure total integrated fluxes for the entire JINGLE sample in 10 infrared/submillimetre bands, including all WISE, Herschel-PACS, Herschel-SPIRE and SCUBA-2 850 ${\mu}m$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 ${\mu}m$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 ${\mu}m$ with a $\geq$ 5${\sigma}$ significance. In the SCUBA-2 850 ${\mu}m$ band we detect 126 galaxies with $\geq$ 3${\sigma}$ significance. The distribution of the JINGLE galaxies in far-infrared/sub-millimetre colour-colour plots reveals that the sample is not well fit by single modified-blackbody models that assume a single dust-emissivity index $(\beta)$. Instead, our new 850 ${\mu}m$ data suggest either that a large fraction of our objects require $\beta < 1.5$, or that a model allowing for an excess of sub-mm emission (e.g., a broken dust emissivity law, or a very cold dust component 10 K) is required. We provide relations to convert far-infrared colours to dust temperature and $\beta$ for JINGLE-like galaxies. For JINGLE the FIR colours correlate more strongly with star-formation rate surface-density rather than the stellar surface-density, suggesting heating of dust is greater due to younger rather than older stellar-populations, consistent with the low proportion of early-type galaxies in the sample.

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity.

Circumbinary disks are common around post-asymptotic giant branch (post-AGB) stars with a stellar companion on orbital time scales of a few 100 to few 1000 days. The presence of a disk is usually inferred from the system's spectral energy distribution (SED), and confirmed, for a sub-sample, by interferometric observations. We used the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument on the Very Large Telescope to obtain extreme adaptive optics assisted scattered light images of the post-AGB binary system AR Puppis. Data have been obtained in the V, I, and H bands. Our observations have produced the first resolved images of AR Puppis' circumbinary disk and confirm its edge-on orientation. In our high angular-resolution and high dynamic-range images we identify several structural components such as a dark mid-plane, the disk surface, and arc-like features. We discuss the nature of these components and use complementary photometric monitoring to relate them to the orbital phase of the binary system. Because the star is completely obscured by the disk at visible wavelengths, we conclude that the long-term photometric variability of the system must be caused by variable scattering, not extinction, of star light by the disk over the binary orbit. Finally, we discuss how the short disk life times and fast evolution of the host stars compared to the ages at which protoplanetary disks are typically observed make systems like AR Puppis valuable extreme laboratories to study circumstellar disk evolution and constrain the time scale of dust grain growth during the planet formation process.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

We present results of the analysis of a sample of 22 stars of spectral types from O7 to B5 and luminosity classes I-V for which Spitzer/IRS spectra are available. The IRS spectra of these stars are examined for signs of excess infrared (IR) emission by comparison with stellar atmospheric spectra. We find that the spectra of half of the studied stars are dominated by excess emission in the far-IR, including all six super- and bright giants. In order to examine the origin of the far-IR excess, we supplement the Spitzer data with optical high-resolution echelle spectroscopy ($\lambda/\Delta \lambda \sim 10^5$), near-IR high-contrast coronagraphic imaging taken with the SPHERE instrument at VLT with a spatial resolution of 0.05", and WISE and Herschel photometry. In the optical region, we detect various absorption and emission lines (H$\alpha$, CIII, and NIII) irrespective of the far-IR excess. Pfund($\alpha$) and Humphrey($\alpha$) lines are observed at the same time as the far-IR excess. These lines are stronger in stars with far-IR excess than in stars without excess. A scattered-light disk in the central r < 2.5" region of the far-IR excess stars HD149404, HD151804, and HD154368 can be excluded from H band imaging down to a 1$\sigma$ contrast of $F(r)/F_{*} \sim 10^{-6}$. The far-IR excess is fit either by a free-free component from ionized gas as for the winds of hot stars or a large (1pc) circumstellar dust shell. The putative dust envelopes required to explain the excess have a visual extinction as low as a few hundred $\mu$-mag.

We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts.

We present JCMT SCUBA-2 $450\mu$m and $850\mu$m observations of 14 Asymptotic Giant Branch (AGB) stars (9 O--rich, 4 C-rich and 1 S--type) and one Red Supergiant (RSG) in the Solar Neighbourhood. We combine these observations with \emph{Herschel}/PACS observations at $70\mu$m and $160\mu$m and obtain azimuthally-averaged surface-brightness profiles and their PSF subtracted residuals. The extent of the SCUBA-2 850 $\mu$m emission ranges from 0.01 to 0.16 pc with an average of $\sim40\%$ of the total flux being emitted from the extended component. By fitting a modified black-body to the four-point SED at each point along the radial profile we derive the temperature ($T$), spectral index of dust emissivity ($\beta$) and dust column density ($\Sigma$) as a function of radius. For all the sources, the density profile deviates significantly from what is expected for a constant mass-loss rate, showing that all the sources have undergone variations in mass-loss during this evolutionary phase. In combination with results from CO line emission, we determined the dust-to-gas mass ratio for all the sources in our sample. We find that, when sources are grouped according to their chemistry, the resulting average dust-to-gas ratios are consistent with the respective canonical values. However we see a range of values with significant scatter which indicate the importance of including spatial information when deriving these numbers.

We analyze results from the first eighteen months of monthly sub-mm monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC53, shows behavior well above the expected measurement uncertainty. Another four sources, two disks and two protostars, show moderately-enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, fifty of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with sigma ~ 0.023. Most sources are not rapidly brightening or fading in the sub-mm. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change sigma = 0.005 would be unobservable in the present sample, whereas an underlying distribution with sigma = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a larger time sample, with a factor of two improvement expected by the conclusion of our 36-month survey.

Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.

The tori around active galactic nuclei (AGN) are potential formation sites for large amounts of dust, and they may help resolve the so-called dust budget crisis at high redshift. We investigate the dust composition in 53 of the 87 Palomar Green (PG) quasars showing the 9.7 micron silicate feature in emission. By simultaneously fitting the mid-infrared spectroscopic features and the underlying continuum, we estimate the mass fraction in various amorphous and crystalline dust species. We find that the dust consists predominantly of alumina and amorphous silicates, with a small fraction in crystalline form. The mean crystallinity is 8 +/- 6%, with more than half of the crystallinities greater than 5%, well above the upper limit determined for the Galaxy. Higher values of crystallinity are found for higher oxide fractions and for more luminous sources.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We performed a detailed spectroscopic analysis of the fullerene C60-containing planetary nebula (PN) Lin49 in the Small Magellanic Cloud using XSHOOTER at the ESO VLT and the Spitzer/IRS instruments. We derived nebular abundances for nine elements. We used TLUSTY to derive photospheric parameters for the central star. Lin49 is C-rich and metal-deficient PN (Z~0.0006). The nebular abundances are in good agreement with Asymptotic Giant Branch nucleosynthesis models for stars with initial mass 1.25 Msun and metallicity Z = 0.001. Using the TLUSTY synthetic spectrum of the central star to define the heating and ionising source, we constructed the photoionisation model with CLOUDY that matches the observed spectral energy distribution (SED) and the line fluxes in the UV to far-IR wavelength ranges simultaneously. We could not fit the ~1-5 um SED using a model with 0.005-0.1 um-sized graphite grains and a constant hydrogen density shell owing to the prominent near-IR excess, while at other wavelengths the model fits the observed values reasonably well. We argue that the near-IR excess might indicate either (1) the presence of very small particles in the form of small carbon clusters, small graphite sheets, or fullerene precursors, or (2) the presence of a high-density structure surrounding the central star. We found that SMC C60 PNe show a near-IR excess component to lesser or greater degree. This suggests that these C60 PNe might maintain a structure nearby their central star.

Many physical parameters in astrophysics are derived using the ratios of two observed quantities. If the relative uncertainties on measurements are small enough, uncertainties can be propagated analytically using simplifying assumptions, but for large normally distributed uncertainties, the probability distribution of the ratio become skewed, with a modal value offset from that expected in Gaussian uncertainty propagation. Furthermore, the most likely value of a ratio A/B is not equal to the reciprocal of the most likely value of B/A. The effect is most pronounced when the uncertainty on the denominator is larger than that on the numerator. We show that this effect is seen in an analysis of 12,126 spectra from the Sloan Digital Sky Survey. The intrinsically fixed ratio of the [O III] lines at 4959 and 5007 ${\AA}$ is conventionally expressed as the ratio of the stronger line to the weaker line. Thus, the uncertainty on the denominator is larger, and non-Gaussian probability distributions result. By taking this effect into account, we derive an improved estimate of the intrinsic 5007/4959 ratio. We obtain a value of 3.012 $\pm$ 0.008, which is slightly but statistically significantly higher than the theoretical value of 2.98. We further investigate the suggestion that fluxes measured from emission lines at low signal to noise are strongly biased upwards. We were unable to detect this effect in the SDSS line flux measurements, and we could not reproduce the results of Rola and Pelat who first described this bias. We suggest that the magnitude of this effect may depend strongly on the specific fitting algorithm used.

Given that a majority of stars form in multiple systems, in order to fully understand the star- and planet-formation processes we must seek to understand them in multiple stellar systems. With this in mind, we present an analysis of the enigmatic binary T-Tauri system VV Corona Australis, in which both components host discs, but only one is visible at optical wavelengths. We seek to understand the peculiarities of this system by searching for a model for the binary which explains all the available continuum observations of the system. We present new mid-infrared interferometry and near-infrared spectroscopy along with archival millimetre-wave observations, which resolve the binary at 1.3mm for the first time. We compute a grid of pre-main-sequence radiative transfer models and calculate their posterior probabilities given the observed spectral energy distributions and mid-infrared interferometric visibilities of the binary components, beginning with the assumption that the only differences between the two components are their inclination and position angles. Our best-fitting solution corresponds to a relatively low luminosity T-tauri binary, with each component's disc having a large scale height and viewed at moderate inclination ($\sim 50^\circ$), with the infrared companion inclined by $\sim 5^\circ$ degrees more than the primary. Comparing the results of our model to evolutionary models suggests stellar masses $\sim 1.7$M$_\odot$ and an age for the system of 3.5Myr, towards the upper end of previous estimates. Combining these results with accretion indicators from near-IR spectroscopy, we determine an accretion rate of 4.0$\times 10^{-8}$M$_\odot$ yr$^{-1}$ for the primary. We suggest that future observations of VV~CrA and similar systems should prioritise high angular resolution sub-mm and near-IR imaging of the discs and high resolution optical/NIR spectroscopy of the central stars.

We present millimetre (SMA) and sub-millimetre (SCUBA-2) continuum observations of the peculiar star KIC 8462852 which displayed several deep and aperiodic dips in brightness during the Kepler mission. Our observations are approximately confusion-limited at 850 $\mu$m and are the deepest millimetre and sub-millimetre photometry of the star that has yet been carried out. No significant emission is detected towards KIC 8462852. We determine upper limits for dust between a few 10$^{-6}$ M$_{\oplus}$ and 10$^{-3}$ M$_{\oplus}$ for regions identified as the most likely to host occluding dust clumps and a total overall dust budget of $<$7.7 M$_{\oplus}$ within a radius of 200 AU. Such low limits for the inner system make the catastrophic planetary disruption hypothesis unlikely. Integrating over the Kepler lightcurve we determine that at least 10$^{-9}$ M$_{\oplus}$ of dust is required to cause the observed Q16 dip. This is consistent with the currently most favoured cometary breakup hypothesis, but nevertheless implies the complete breakup of $\sim$ 30 Comet 1/P Halley type objects. Finally, in the wide SCUBA-2 field-of-view we identify another candidate debris disc system that is potentially the largest yet discovered.

Massive stars live short lives, losing large amounts of mass through their stellar wind. Their mass is a key factor determining how and when they explode as supernovae, enriching the interstellar medium with heavy elements and dust. During the red supergiant phase, mass-loss rates increase prodigiously, but the driving mechanism has proven elusive. Here we present high-contrast optical polarimetric-imaging observations of the extreme red supergiant VY Canis Majoris and its clumpy, dusty, mass-loss envelope, using the new extreme-adaptive-optics instrument SPHERE at the VLT. These observations allow us to make the first direct and unambiguous detection of submicron dust grains in the ejecta; we derive an average grain radius $\sim$ 0.5 $\mu$m, 50 times larger than in the diffuse ISM, large enough to receive significant radiation pressure by photon scattering. We find evidence for varying grain sizes throughout the ejecta, highlighting the dynamical nature of the envelope. Grains with 0.5 $\mu$m sizes are likely to reach a safe distance from the eventual explosion of VY Canis Majoris; hence it may inject upwards of 10$^{-2}$ M$_\odot$ of dust into the ISM.

(abridged) The dust content of the universe is primarily explored via its interaction with stellar photons, producing interstellar extinction. However, owing to the physical extension of the observing beam, observations may detect scattered photons, resulting in a change in the observed (or effective) extinction, depending on the spatial distribution of the dust and the resolution of the instrument. We investigate the influence of clumpy dust distributions on effective extinction toward embedded sources and those in the diffuse ISM. We use Monte Carlo radiative transfer to examine effective extinction for various geometries. By varying the number, optical depth and volume-filling factor of clumps in models of spherical shells and the diffuse ISM, we explore the evolution of extinction. Depending on the number of scattering events in the beam, the extinction curve steepens in homogeneous media and flattens in clumpy media. As a result, clumpy dust distributions can to reproduce extinction curves with arbitrary R_v, the effective ratio of total-to-selective extinction. The flattening `washes out' the 2175 \AA bump and shifts the peak to shorter wavelengths. The mean R_v of a shell correlates with the optical depth of an individual clump and the wavelength at which a clump becomes optically thick. Similar behaviour is seen for edge-on discs or tori. However, at grazing inclinations the combination of extinction and strong forward scattering results in chaotic behaviour. Caution is therefore advised when measuring extinction in, for example, AGN tori or toward SNIa or GRB afterglows. In face-on discs, the shape of the scattered continuum changes with clumpiness, however, unlike absorption features, individual features in the scattering cross-sections are preserved. Finally, we show that diffuse interstellar extinction is not modified by scattering on scales of a few kpc.

CONTEXT. Little is known about the planetary systems around single white dwarfs although there is strong evidence that they do exist. AIMS. We performed a pilot study with the extreme-AO system on the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) on the Very Large Telescopes (VLT) to look for giant planets around a young white dwarf, GD 50. METHODS. We were awarded science verification time on the new ESO instrument SPHERE. Observations were made with the InfraRed Dual-band Imager and Spectrograph in classical imaging mode in H band. RESULTS. Despite the faintness of the target (14.2 mag in R band), the AO loop was closed and a strehl of 37\% was reached in H band. No objects were detected around GD 50. We achieved a 5-sigma contrast of 6.2, 8.0 and 8.25 mags at 0{\farcs}2, 0{\farcs}4 and 0{\farcs}6 and beyond, respectively. We exclude any substellar objects more massive than 4.0 M$_\textrm{J}$ at 6.2 AU, 2.9 M$_\textrm{J}$ at 12.4 AU and 2.8 M$_\textrm{J}$ at 18.6 AU and beyond. This rivals the previous upper limit set by Spitzer. We further show that SPHERE is the most promising instrument available to search for close-in substellar objects around nearby white dwarfs.

Young stars show evidence of accretion discs which evolve quickly and disperse with an e-folding time of $\sim$ 3Myr. This is in striking contrast with recent observations that suggest evidence for numerous $>30$ Myr old stars with an accretion disc in large star-forming complexes. We consider whether these observations of apparently old accretors could be explained by invoking Bondi-Hoyle accretion to rebuild a new disc around these stars during passage through a clumpy molecular cloud. We combine a simple Monte Carlo model to explore the capture of mass by such systems with a viscous evolution model to infer the levels of accretion that would be observed. We find that a significant fraction of stars may capture enough material via the Bondi-Hoyle mechanism to rebuild a disc of mass $\gtrsim$ 1 minimum-mass solar nebula, and $\lesssim 10\%$ accrete at observable levels at any given time. A significant fraction of the observed old accretors may be explained with our proposed mechanism. Such accretion may provide a chance for a second epoch of planet formation, and have unpredictable consequences for planetary evolution.

Whenever observations are compared to theories, an estimate of the uncertainties associated with the observations is vital if the comparison is to be meaningful. However, many determinations of temperatures, densities and abundances in photoionized nebulae do not quote the associated uncertainty. Those that do typically propagate the uncertainties using analytical techniques which rely on assumptions that generally do not hold. Motivated by this issue, we have developed NEAT (Nebular Empirical Analysis Tool), a new code for calculating chemical abundances in photoionized nebulae. The code carries out an analysis of lists of emission lines using long-established techniques to estimate the amount of interstellar extinction, calculate representative temperatures and densities, compute ionic abundances from both collisionally excited lines and recombination lines, and finally to estimate total elemental abundances using an ionization correction scheme. NEAT uses a Monte Carlo technique to robustly propagate uncertainties from line flux measurements through to the derived abundances. We show that for typical observational data, this approach is superior to analytic estimates of uncertainties. NEAT also accounts for the effect of upward biasing on measurements of lines with low signal to noise, allowing us to accurately quantify the effect of this bias on abundance determinations. We find not only that the effect can result in significant over-estimates of heavy element abundances derived from weak lines, but that taking it into account reduces the uncertainty of these abundance determinations. Finally, we investigate the effect of possible uncertainties in R, the ratio of selective to total extinction, on abundance determinations. We find that the uncertainty due to this parameter is negligible compared to the statistical uncertainties due to typical line flux measurement uncertainties.