We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium.

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.

The dependence of polarization fraction $p$ on total intensity $I$ in polarized submillimeter emission measurements is typically parameterized as $p\propto I^{-\alpha}$ $(\alpha \leq 1)$, and used to infer dust grain alignment efficiency in star-forming regions, with an index $\alpha=1$ indicating near-total lack of alignment of grains with the magnetic field. In this work we demonstrate that the non-Gaussian noise characteristics of polarization fraction may produce apparent measurements of $\alpha \sim 1$ even in data with significant signal-to-noise in Stokes $Q$, $U$ and $I$ emission, and so with robust measurements of polarization angle. We present a simple model demonstrating this behavior, and propose a criterion by which well-characterized measurements of polarization fraction may be identified. We demonstrate that where our model is applicable, $\alpha$ can be recovered by fitting the $p-I$ relationship with the mean of the Rice distribution, without statistical debiasing of polarization fraction. We apply our model to JCMT BISTRO Survey POL-2 850$\mu$m observations of three clumps in the Ophiuchus Molecular Cloud, finding that in the externally-illuminated Oph A region, $\alpha\approx 0.34$, while in the more isolated Oph B and C, despite their differing star formation histories, $\alpha \sim 0.6-0.7$. Our results thus suggest that dust grain alignment in dense gas is more strongly influenced by incident interstellar radiation field than by star formation history. We further find that grains may remain aligned with the magnetic field at significantly higher gas densities than has previously been believed, thus allowing investigation of magnetic field properties within star-forming clumps and cores.

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.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

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.

The dust properties inferred from the analysis of Planck observations in total and polarized emission challenge current dust models. We propose new dust models compatible with polarized and unpolarized data in extinction and emission for translucent lines of sight ($0.5 < A_V < 2.5$). We amended the DustEM tool to model polarized extinction and emission. We fit the spectral dependence of the mean extinction, polarized extinction, SED and polarized SED with PAHs, astrosilicates and amorphous carbon (a-C). The astrosilicate population is aligned along the magnetic field lines, while the a-C population may be aligned or not. With their current optical properties, oblate astrosilicate grains are not emissive enough to reproduce the emission to extinction polarization ratio $P_{353}/p_V$ derived with Planck data. Models using prolate astrosilicate grains with an elongation $a/b=3$ and an inclusion of 20% of porosity succeed. The spectral dependence of the polarized SED is steeper in our models than in the data. Models perform slightly better when a-C grains are aligned. A small (6%) volume inclusion of a-C in the astrosilicate matrix removes the need for porosity and perfect grain alignment, and improves the fit to the polarized SED. Dust models based on astrosilicates can be reconciled with Planck data by adapting the shape of grains and adding inclusions of porosity or a-C in the astrosilicate matrix.

Polarized extinction and emission from dust in the interstellar medium (ISM) are hard to interpret, as they have a complex dependence on dust optical properties, grain alignment and magnetic field orientation. This is particularly true in molecular clouds. The data available today are not yet used to their full potential. The combination of emission and extinction, in particular, provides information not available from either of them alone. We combine data from the scientific literature on polarized dust extinction with Planck data on polarized emission and we use them to constrain the possible variations in dust and environmental conditions inside molecular clouds, and especially translucent lines of sight, taking into account magnetic field orientation. We focus on the dependence between \lambda_max -- the wavelength of maximum polarization in extinction -- and other observables such as the extinction polarization, the emission polarization and the ratio of the two. We set out to reproduce these correlations using Monte-Carlo simulations where the relevant quantities in a dust model -- grain alignment, size distribution and magnetic field orientation -- vary to mimic the diverse conditions expected inside molecular clouds. None of the quantities chosen can explain the observational data on its own: the best results are obtained when all quantities vary significantly across and within clouds. However, some of the data -- most notably the stars with low emission-to-extinction polarization ratio -- are not reproduced by our simulation. Our results suggest not only that dust evolution is necessary to explain polarization in molecular clouds, but that a simple change in size distribution is not sufficient to explain the data, and point the way for future and more sophisticated models.

We compare the performance of several dust models in reproducing the dust spectral energy distribution (SED) per unit extinction in the diffuse interstellar medium (ISM). We use our results to constrain the variability of the optical properties of big grains in the diffuse ISM, as published by the Planck collaboration. We use two different techniques to compare the predictions of dust models to data from the Planck HFI, IRAS and SDSS surveys. First, we fit the far-infrared emission spectrum to recover the dust extinction and the intensity of the interstellar radiation field (ISRF). Second, we infer the ISRF intensity from the total power emitted by dust per unit extinction, and then predict the emission spectrum. In both cases, we test the ability of the models to reproduce dust emission and extinction at the same time. We identify two issues. Not all models can reproduce the average dust emission per unit extinction: there are differences of up to a factor $\sim2$ between models, and the best accord between model and observation is obtained with the more emissive grains derived from recent laboratory data on silicates and amorphous carbons. All models fail to reproduce the variations in the emission per unit extinction if the only variable parameter is the ISRF intensity: this confirms that the optical properties of dust are indeed variable in the diffuse ISM. Diffuse ISM observations are consistent with a scenario where both ISRF intensity and dust optical properties vary. The ratio of the far-infrared opacity to the $V$ band extinction cross-section presents variations of the order of $\sim20\%$ ($40-50\%$ in extreme cases), while ISRF intensity varies by $\sim30\%$ ($\sim60\%$ in extreme cases). This must be accounted for in future modelling.

The Planck-HFI all-sky survey from 353 to 857GHz combined with the 100 microns IRAS show that the dust properties vary in the diffuse ISM at high Galactic latitude (1e19<NH<2.5e20 H/cm2). Our aim is to explain these variations with changes in the ISM properties and grain evolution. Our starting point is the latest core-mantle dust model. It consists of small aromatic-rich carbon grains, larger amorphous carbon grains with aliphatic-rich cores and aromatic-rich mantles, and amorphous silicates with Fe/FeS nano-inclusions covered by aromatic-rich carbon mantles. We explore whether variations in the radiation field or in the gas density distribution in the diffuse ISM could explain the observations. The dust properties are also varied in terms of mantle thickness, Fe/FeS inclusions, carbon abundance, and size distribution. Variations in the radiation field intensity and gas density distribution cannot explain the observed variations but radiation fields harder than the standard ISRF may participate in creating part of them. We further show that variations in the grain mantle thickness coupled with changes in the grain size distribution can reproduce most of the observations. We put a limit on the mantle thickness of the silicates (~10-15nm), and find that aromatic-rich mantles are needed for the carbon grains (at least 5-7.5nm thick). We also find that changes in the carbon abundance in the grains could explain part of the observed variations. Finally, we show that varying the composition of Fe/FeS inclusions in the silicates cannot account for the variations. With small variations in the dust properties, we are able to explain most of the variations in the dust emission observed by Planck-HFI in the diffuse ISM. We also find that the small realistic changes in the dust properties that we consider almost perfectly match the anti-correlation and scatter in the observed beta-T relation.

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $\tau_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/\tau_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required.

Context. The evolution of amorphous hydrocarbon materials, a-C(:H), principally resulting from ultraviolet (UV) photon absorption- induced processing, are likely at the heart of the variations in the observed properties of dust in the interstellar medium. Aims. The consequences of the size-dependent and compositional variations in a-C(:H), from aliphatic-rich a-C:H to aromatic-rich a-C, are studied within the context of the interstellar dust extinction and emission. Methods. Newly-derived optical property data for a-C(:H) materials, combined with that for an amorphous forsterite-type silicate with iron nano-particle inclusions, a-SilFe, are used to explore dust evolution in the interstellar medium. Results. We present a new dust model that consists of a power-law distribution of small a-C grains and log-normal distributions of large a-SilFe and a-C(:H) grains. The model, which is firmly anchored by laboratory-data, is shown to quite naturally explain the variations in the infrared (IR) to far-ultraviolet (FUV) extinction, the 217 nm UV bump, the IR absorption and emission bands and the IR-mm dust emission. Conclusions. The major strengths of the new model are its inherent simplicity and built-in capacity to follow dust evolution in interstellar media. We show that mantle accretion in molecular clouds and UV photo-processing in photo-dominated regions are likely the major drivers of dust evolution.

Observational evidence seems to indicate that the depletion of interstellar carbon into dust shows rather wide variations and that carbon undergoes rather rapid recycling in the interstellar medium (ISM). Small hydrocarbon grains are processed in photo-dissociation regions by UV photons, by ion and electron collisions in interstellar shock waves and by cosmic rays. A significant fraction of hydrocarbon dust must therefore be re-formed by accretion in the dense, molecular ISM. A new dust model (Jones et al., Astron. Astrophys., 2013, 558, A62) shows that variations in the dust observables in the diffuse interstellar medium (nH = 1000 cm^3), can be explained by systematic and environmentally-driven changes in the small hydrocarbon grain population. Here we explore the consequences of gas-phase carbon accretion onto the surfaces of grains in the transition regions between the diffuse ISM and molecular clouds (e.g., Jones, Astron. Astrophys., 2013, 555, A39). We find that significant carbonaceous dust re-processing and/or mantle accretion can occur in the outer regions of molecular clouds and that this dust will have significantly different optical properties from the dust in the adjacent diffuse ISM. We conclude that the (re-)processing and cycling of carbon into and out of dust is perhaps the key to advancing our understanding of dust evolution in the ISM.

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $\psi$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]