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.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

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.

The JCMT Gould Belt Survey was one of the first Legacy Surveys with the James Clerk Maxwell Telescope in Hawaii, mapping 47 square degrees of nearby (< 500 pc) molecular clouds in both dust continuum emission at 850 $\mu$m and 450 $\mu$m, as well as a more-limited area in lines of various CO isotopologues. While molecular clouds and the material that forms stars have structures on many size scales, their larger-scale structures are difficult to observe reliably in the submillimetre regime using ground-based facilities. In this paper, we quantify the extent to which three subsequent data-reduction methods employed by the JCMT GBS accurately recover emission structures of various size scales, in particular, dense cores which are the focus of many GBS science goals. With our current best data-reduction procedure, we expect to recover 100% of structures with Gaussian sigma sizes of $\le$30" and intensity peaks of at least five times the local noise for isolated peaks of emission. The measured sizes and peak fluxes of these compact structures are reliable (within 15% of the input values), but source recovery and reliability both decrease significantly for larger emission structures and for fainter peaks. Additional factors such as source crowding have not been tested in our analysis. The most recent JCMT GBS data release includes pointing corrections, and we demonstrate that these tend to decrease the sizes and increase the peak intensities of compact sources in our dataset, mostly at a low level (several percent), but occasionally with notable improvement.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present the first high-resolution, submillimeter-wavelength polarimetric observations of -- and thus direct observations of the magnetic field morphology within -- the dense gas of the Pillars of Creation in M16. These 850$\,\mu$m observations, taken as part of the BISTRO (B-Fields in Star-forming Region Observations) Survey using the POL-2 polarimeter on the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT), show that the magnetic field runs along the length of the pillars, perpendicular to, and decoupled from, the field in the surrounding photoionized cloud. Using the Chandrasekhar-Fermi method we estimate a plane-of-sky magnetic field strength of $170-320\,\mu$G in the Pillars, consistent with their having been formed through compression of gas with initially weak magnetization. The observed magnetic field strength and morphology suggests that the magnetic field may be slowing the pillars' evolution into cometary globules. We thus hypothesize that the evolution and lifetime of the Pillars may be strongly influenced by the strength of the coupling of their magnetic field to that of their parent photoionized cloud -- i.e. that the Pillars' longevity results from magnetic support

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

The crust of accreting neutron stars plays a central role in many different observational phenomena. In these stars, heavy elements produced by H-He burning in the rapid proton capture (rp-) process continually freeze to form new crust. In this paper, we explore the expected composition of the solid phase. We first demonstrate using molecular dynamics that two distinct types of chemical separation occur, depending on the composition of the rp-process ashes. We then calculate phase diagrams for three-component mixtures and use them to determine the allowed crust compositions. We show that, for the large range of atomic numbers produced in the rp-process ($Z\sim 10$--$50$), the solid that forms has only a small number of available compositions. We conclude that accreting neutron star crusts should be polycrystalline, with domains of distinct composition. Our results motivate further work on the size of the compositional domains, and have implications for crust physics and accreting neutron star phenomenology.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

Nucleosynthesis in the oceans of accreting neutron stars can produce novel mixtures of nuclides, whose composition is dependent on the exact astrophysical conditions. Many simulations have now been done to determine the nucleosynthesis products in the ocean, but the phase separation at the base of the ocean, which determines the composition of the crust, has not been as well studied. In this work, we simulate the phase separation of a composition, which was predicted to produce a crust enriched in light nuclei, in contrast with past work which predicts that crust is enriched in heavy nuclei. We perform molecular dynamics simulations of the phase separation of this mixture using the methods of Horowitz $\textit{et. al.}$ (2007). We find good agreement with the predictions of Mckinven $\textit{et al.}$ (2016) for the phase separation of this mixture. Moreover, this supports their method as a computationally efficient alternative to molecular dynamics for calculating phase separation for a wider regime of astrophysical conditions.

We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases the neutrino signal at late times of 10 or more seconds after stellar core collapse. This signal, for a galactic SN, should be clearly visible in large detectors such as Super-Kamiokande.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

The Orion A molecular cloud is one of the most well-studied nearby star-forming regions, and includes regions of both highly clustered and more dispersed star formation across its full extent. Here, we analyze dense, star-forming cores identified in the 850 {\mu}m and 450 {\mu}m SCUBA-2 maps from the JCMT Gould Belt Legacy Survey. We identify dense cores in a uniform manner across the Orion A cloud and analyze their clustering properties. Using two independent lines of analysis, we find evidence that clusters of dense cores tend to be mass segregated, suggesting that stellar clusters may have some amount of primordial mass segregation already imprinted in them at an early stage. We also demonstrate that the dense core clusters have a tendency to be elongated, perhaps indicating a formation mechanism linked to the filamentary structure within molecular clouds.

We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

In view of increased interest in object-oriented systems for describing coordinate information, we present a description of the data model used by the Starlink AST library. AST provides a comprehensive range of facilities for attaching world co-ordinate systems to astronomical data, and for retrieving and interpreting that information in a variety of formats, including FITS-WCS. AST is a mature system that has been in use for more than 17 years, and may consequently be useful as a means of informing development of similar systems in the future.

[Background] The "pasta" phase of nuclear matter may play an important role in the structure and evolution of neutron stars. Recent works suggest nuclear pasta has a high resistivity which could be explained by the presence of long lived topological defects. The defects act as impurities that decrease thermal and electrical conductivity of the pasta. [Purpose] To quantify how topological defects affect transport properties of nuclear pasta and estimate this effect using an impurity parameter $Q_{\text{imp}}$. [Methods] Contrast molecular dynamics simulations of up to 409\,600 nucleons arranged in parallel nuclear pasta slabs (perfect pasta) with simulations of pasta slabs connected by topological defects (impure pasta). From these simulations compare the viscosity and heat conductivity of perfect and impure pasta to obtain an effective impurity parameter $Q_{\text{imp}}$ due to the presence of defects. [Results] Both the viscosity and thermal conductivity calculated for both perfect and impure pasta are anisotropic, peaking along directions perpendicular to the slabs and reaching a minimum close to zero parallel to them. In our 409\,600 nucleon simulation topological defects connecting slabs of pasta reduce both the thermal conductivity and viscosity on average by about 37\%. We estimate an effective impurity parameter due to the defects of order $Q_{\text{imp}}\sim30$. [Conclusions] Topological defects in the pasta phase of nuclear matter have an effect similar to impurities in a crystal lattice. The irregularities introduced by the defects reduce the thermal and electrical conductivities and the viscosity of the system. This effect can be parameterized by a large impurity parameter $Q_{\text{imp}}\sim30$.

The JCMT Gould Belt Legacy Survey obtained SCUBA-2 observations of dense cores within three sub-regions of Orion B: LDN 1622, NGC 2023/2024, and NGC 2068/2071, all of which contain clusters of cores. We present an analysis of the clustering properties of these cores, including the two-point correlation function and Cartwright's Q parameter. We identify individual clusters of dense cores across all three regions using a minimal spanning tree technique, and find that in each cluster, the most massive cores tend to be centrally located. We also apply the independent M-Sigma technique and find a strong correlation between core mass and the local surface density of cores. These two lines of evidence jointly suggest that some amount of mass segregation in clusters has happened already at the dense core stage.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artifacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software, Starlink, but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only 3 times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

A striking shape was recently observed for the cellular organelle endoplasmic reticulum consisting of stacked sheets connected by helical ramps. This shape is interesting both for its biological function, to synthesize proteins using an increased surface area for ribosome factories, and its geometric properties that may be insensitive to details of the microscopic interactions. In the present work, we find very similar shapes in our molecular dynamics simulations of the nuclear pasta phases of dense nuclear matter that are expected deep in the crust of neutron stars. There are dramatic differences between nuclear pasta and terrestrial cell biology. Nuclear pasta is 14 orders of magnitude denser than the aqueous environs of the cell nucleus and involves strong interactions between protons and neutrons, while cellular scale biology is dominated by the entropy of water and complex assemblies of biomolecules. Nonetheless the very similar geometry suggests both systems may have similar coarse-grained dynamics and that the shapes are indeed determined by geometrical considerations, independent of microscopic details. Many of our simulations self-assemble into flat sheets connected by helical ramps. These ramps may impact the thermal and electrical conductivities, viscosity, shear modulus, and breaking strain of neutron star crust. The interaction we use, with Coulomb frustration, may provide a simple model system that reproduces many biologically important shapes.

With the advent of modern multi-detector heterodyne instruments that can result in observations generating thousands of spectra per minute it is no longer feasible to reduce these data as individual spectra. We describe the automated data reduction procedure used to generate baselined data cubes from heterodyne data obtained at the James Clerk Maxwell Telescope. The system can automatically detect baseline regions in spectra and automatically determine regridding parameters, all without input from a user. Additionally it can detect and remove spectra suffering from transient interference effects or anomalous baselines. The pipeline is written as a set of recipes using the ORAC-DR pipeline environment with the algorithmic code using Starlink software packages and infrastructure. The algorithms presented here can be applied to other heterodyne array instruments and have been applied to data from historical JCMT heterodyne instrumentation.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard is, however, showing its age. Developed in the late 1970s the FITS authors made a number of implementation choices for the format that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not appreciate the challenges which we would be facing today in astronomical computing. Difficulties we now face include, but are not limited to, having to address the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets and the need to capture significantly more complex metadata and data relationships. There are members of the community today who find some (or all) of these limitations unworkable, and have decided to move ahead with storing data in other formats. This reaction should be taken as a wakeup call to the FITS community to make changes in the FITS standard, or to see its usage fall. In this paper we detail some selected important problems which exist within the FITS standard today. It is not our intention to prescribe specific remedies to these issues; rather, we hope to call attention of the FITS and greater astronomical computing communities to these issues in the hopes that it will spur action to address them.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them.

Background: Exotic non-spherical nuclear pasta shapes are expected in nuclear matter at just below saturation density because of competition between short range nuclear attraction and long range Coulomb repulsion. Purpose: We explore the impact of nuclear pasta on nucleosynthesis, during neutron star mergers, as cold dense nuclear matter is ejected and decompressed. Methods: We perform classical molecular dynamics simulations with 51200 and 409600 nucleons, that are run on GPUs. We expand our simulation region to decompress systems from an initial density of 0.080 fm^{-3} down to 0.00125 fm^{-3}. We study proton fractions of Y_P=0.05, 0.10, 0.20, 0.30, and 0.40 at T =0.5, 0.75, and 1.0 MeV. We calculate the composition of the resulting systems using a cluster algorithm. Results: We find final compositions that are in good agreement with nuclear statistical equilibrium models for temperatures of 0.75 and 1 MeV. However, for proton fractions greater than Y_P=0.2 at a temperature of T = 0.5 MeV, the MD simulations produce non-equilibrium results with large rod-like nuclei. Conclusions: Our MD model is valid at higher densities than simple nuclear statistical equilibrium models and may help determine the initial temperatures and proton fractions of matter ejected in mergers.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

This paper describes the FellWalker algorithm, a watershed algorithm that segments a 1-, 2- or 3-dimensional array of data values into a set of disjoint clumps of emission, each containing a single significant peak. Pixels below a nominated constant data level are assumed to be background pixels and are not assigned to any clump. FellWalker is thus equivalent in purpose to the CLUMPFIND algorithm. However, unlike CLUMPFIND, which segments the array on the basis of a set of evenly-spaced contours and thus uses only a small fraction of the available data values, the FellWalker algorithm is based on a gradient-tracing scheme which uses all available data values. Comparisons of CLUMPFIND and FellWalker using a crowded field of artificial Gaussian clumps, all of equal peak value and width, suggest that the results produced by FellWalker are less dependent on specific parameter settings than are those of CLUMPFIND.

The dense neutron-rich matter found in supernovae and neutron stars is expected to form complex nonuniform phases referred to as nuclear pasta. The pasta shapes depend on density, temperature and proton fraction and determine many transport properties in supernovae and neutron star crusts. We use two recently developed hybrid CPU/GPU codes to perform large scale molecular dynamics (MD) simulations with $51200$ and $409600$ nucleons of nuclear pasta. From the output of the MD simulations we characterize the topology and compute two observables, the radial distribution function $g(r)$ and the structure factor $S(q)$, for systems with proton fractions $Y_p=0.10, 0.20, 0.30$ and $0.40$ at about one third of nuclear saturation density and temperatures near $1.0$ MeV. We observe that the two lowest proton fraction systems simulated, $Y_p=0.10$ and $0.20$, equilibrate quickly and form liquid-like structures. Meanwhile, the two higher proton fraction systems, $Y_p=0.30$ and $0.40$, take a longer time to equilibrate and organize themselves in solid-like periodic structures. Furthermore, the $Y_p=0.40$ system is made up of slabs, lasagna phase, interconnected by defects while the $Y_p=0.30$ systems consist of a stack of perforated plates, the nuclear waffle phase. The periodic configurations observed in our MD simulations for proton fractions $Y_p\ge0.30$ have important consequences for the structure factors $S(q)$ of protons and neutrons, which relate to many transport properties of supernovae and neutron star crust. A detailed study of the waffle phase and how its structure depends on temperature, size of the simulation and the screening length showed that finite-size effects appear to be under control and, also, that the plates in the waffle phase merge at temperatures slightly above $1.0$ MeV and the holes in the plates form an hexagonal lattice at temperatures slightly lower than $1.0$ MeV.

The extensible N-Dimensional Data Format (NDF) was designed and developed in the late 1980s to provide a data model suitable for use in a variety of astronomy data processing applications supported by the UK Starlink Project. Starlink applications were used extensively, primarily in the UK astronomical community, and form the basis of a number of advanced data reduction pipelines today. This paper provides an overview of the historical drivers for the development of NDF and the lessons learned from using a defined hierarchical data model for many years in data reduction software, data pipelines and in data acquisition systems.

Nuclear pasta, with non-spherical shapes, is expected near the base of the crust in neutron stars. Large scale molecular dynamics simulations of pasta show long lived topological defects that could increase electron scattering and reduce both the thermal and electrical conductivities. We model a possible low conductivity pasta layer by increasing an impurity parameter Q_{imp}. Predictions of light curves for the low mass X-ray binary MXB 1659-29, assuming a large Q_{imp}, find continued late time cooling that is consistent with Chandra observations. The electrical and thermal conductivities are likely related. Therefore observations of late time crust cooling can provide insight on the electrical conductivity and the possible decay of neutron star magnetic fields (assuming these are supported by currents in the crust).

The formation of complex nonuniform phases of nuclear matter, known as nuclear pasta, is studied with molecular dynamics simulations containing 51200 nucleons. A phenomenological nuclear interaction is used that reproduces the saturation binding energy and density of nuclear matter. Systems are prepared at an initial density of 0.10fm$^{-3}$ and then the density is decreased by expanding the simulation volume at different rates to densities of 0.01 fm$^{-3}$ or less. An originally uniform system of nuclear matter is observed to form spherical bubbles ("swiss cheese"), hollow tubes, flat plates ("lasagna"), thin rods ("spaghetti") and, finally, nearly spherical nuclei with decreasing density. We explicitly observe nucleation mechanisms, with decreasing density, for these different pasta phase transitions. Topological quantities known as Minkowski functionals are obtained to characterize the pasta shapes. Different pasta shapes are observed depending on the expansion rate. This indicates non equilibrium effects. We use this to determine the best ways to obtain lower energy states of the pasta system from MD simulations and to place constrains on the equilibration time of the system.

The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is an instrument operating on the 15-m James Clerk Maxwell Telescope, nominally consisting of 5120 bolometers in each of two simultaneous imaging bands centred over 450 and 850 um. The camera is operated by scanning across the sky and recording data at a rate of 200 Hz. As the largest of a new generation of multiplexed kilopixel bolometer cameras operating in the (sub)millimetre, SCUBA-2 data analysis represents a significant challenge. We describe the production of maps using the Sub-Millimetre User Reduction Facility (SMURF) in which we have adopted a fast, iterative approach to map-making that enables data reduction on single, modern, high-end desktop computers, with execution times that are typically shorter than the observing times. SMURF is used in an automated setting, both at the telescope for real-time feedback to observers, as well as for the production of science products for the JCMT Science Archive at the Canadian Astronomy Data Centre. Three detailed case studies are used to: (i) explore convergence properties of the map-maker using simple prior constraints (Uranus -- a point source); (ii) achieve the white-noise limit for faint point-source studies (extragalactic blank-field survey of the Lockman Hole); and (iii) demonstrate that our strategy is capable of recovering angular scales comparable to the size of the array footprint (approximately 5 arcmin) for bright extended sources (star-forming region M17).

SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 microns, the vast increase in pixel count means that SCUBA-2 maps the sky 100-150 times faster than the previous SCUBA instrument. In this paper we present an overview of the instrument, discuss the physical characteristics of the superconducting detector arrays, outline the observing modes and data acquisition, and present the early performance figures on the telescope. We also showcase the capabilities of the instrument via some early examples of the science SCUBA-2 has already undertaken. In February 2012, SCUBA-2 began a series of unique legacy surveys for the JCMT community. These surveys will take 2.5 years and the results are already providing complementary data to the shorter wavelength, shallower, larger-area surveys from Herschel. The SCUBA-2 surveys will also provide a wealth of information for further study with new facilities such as ALMA, and future telescopes such as CCAT and SPICA.

SCUBA-2 is a 10000-bolometer submillimetre camera on the James Clerk Maxwell Telescope (JCMT). The instrument commissioning was completed in September 2011, and full science operations began in October 2011. To harness the full potential of this powerful new astronomical tool, the instrument calibration must be accurate and well understood. To this end, the algorithms for calculating the line-of-sight opacity have been improved, and the derived atmospheric extinction relationships at both wavebands of the SCUBA-2 instrument are presented. The results from over 500 primary and secondary calibrator observations have allowed accurate determination of the flux conversion factors (FCF) for the 850 and 450 micron arrays. Descriptions of the instrument beam-shape and photometry methods are presented. The calibration factors are well determined, with relative calibration accuracy better than 5 per cent at 850 microns and 10 per cent at 450 microns, reflecting the success of the derived opacity relations as well as the stability of the performance of the instrument over several months. The sample-size of the calibration observations and accurate FCFs have allowed the determination of the 850 and 450 micron fluxes of several well-known submillimetre sources, and these results are compared with previous measurements from SCUBA.

The neutron rich isotope 22Ne may be a significant impurity in carbon and oxygen white dwarfs and could impact how the stars freeze. We perform molecular dynamics simulations to determine the influence of 22Ne in carbon-oxygen-neon systems on liquid-solid phase equilibria. Both liquid and solid phases are present simultaneously in our simulation volumes. We identify liquid, solid, and interface regions in our simulations using a bond angle metric. In general we find good agreement for the composition of liquid and solid phases between our MD simulations and the semi analytic model of Medin and Cumming. The trace presence of a third component, neon, does not appear to strongly impact the chemical separation found previously for two component carbon and oxygen systems. This suggests that small amounts of 22Ne may not qualitatively change how the material in white dwarf stars freezes. However, we do find systematically lower melting temperatures (higher Gamma) in our MD simulations compared to the semi analytic model. This difference seems to grow with impurity parameter Q_imp and suggests a problem with simple corrections to the linear mixing rule for the free energy of multicomponent solid mixtures that is used in the semi analytic model.

Recent developments in the AST library are described, including a Python interface, support for the FITS-WCS "-TAB" system for storing tabular co-ordinate information, and extended support for representing distortions in spatial projections, using several schemes in common use (IRAF TNX/ZPX, Spitzer SIP, NOAO TPV and SCAMP).

SCUBA-2 is the largest submillimetre wide-field bolometric camera ever built. This 43 square arc-minute field-of-view instrument operates at two wavelengths (850 and 450 microns) and has been installed on the James Clerk Maxwell Telescope on Mauna Kea, Hawaii. SCUBA-2 has been successfully commissioned and operational for general science since October 2011. This paper presents an overview of the on-sky performance of the instrument during and since commissioning in mid-2011. The on-sky noise characteristics and NEPs of the 450 and 850 micron arrays, with average yields of approximately 3400 bolometers at each wavelength, will be shown. The observing modes of the instrument and the on-sky calibration techniques are described. The culmination of these efforts has resulted in a scientifically powerful mapping camera with sensitivities that allow a square degree of sky to be mapped to 10 mJy/beam rms at 850 micron in 2 hours and 60 mJy/beam rms at 450 micron in 5 hours in the best weather.

We determine the liquid-solid phase diagram for carbon-oxygen and oxygen-selenium plasma mixtures using two-phase MD simulations. We identified liquid, solid, and interface regions using a bond angle metric. To study finite size effects, we perform 27648 and 55296 ion simulations. To help monitor non-equilibrium effects, we calculate diffusion constants $D_i$. For the carbon-oxygen system we find that $D_O$ for oxygen ions in the solid is much smaller than $D_C$ for carbon ions and that both diffusion constants are 80 or more times smaller than diffusion constants in the liquid phase. There is excellent agreement between our carbon-oxygen phase diagram and that predicted by Medin and Cumming. This suggests that errors from finite size and non-equilibrium effects are small and that the carbon-oxygen phase diagram is now accurately known. The oxygen-selenium system is a simple two-component model for more complex rapid proton capture nucleosynthesis ash compositions for an accreting neutron star. Diffusion of oxygen, in a predominately selenium crystal, is remarkably fast, comparable to diffusion in the liquid phase. We find a somewhat lower melting temperature for the oxygen-selenium system than that predicted by Medin and Cumming. This is probably because of electron screening effects.

In this book chapter we review plasma crystals in the laboratory, in the interior of white dwarf stars, and in the crust of neutron stars. We describe a molecular dynamics formalism and show results for many neutron star crust properties including phase separation upon freezing, diffusion, breaking strain, shear viscosity and dynamics response of nuclear pasta. We end with a summary and discuss open questions and challenges for the future.

Diffusion in coulomb crystals can be important for the structure of neutron star crusts. We determine diffusion constants $D$ from molecular dynamics simulations. We find that $D$ for coulomb crystals with relatively soft-core $1/r$ interactions may be larger than $D$ for Lennard-Jones or other solids with harder-core interactions. Diffusion, for simulations of nearly perfect body-centered-cubic lattices, involves the exchange of ions in ring-like configurations. Here ions "hop" in unison without the formation of long lived vacancies. Diffusion, for imperfect crystals, involves the motion of defects. Finally, we find that diffusion, for an amorphous system rapidly quenched from coulomb parameter $\Gamma=175$ to coulomb parameters up to $\Gamma=1750$, is fast enough so that the system starts to crystallize during long simulation runs. These results strongly suggest that coulomb solids in cold white dwarf stars, and the crust of neutron stars, will be crystalline and not amorphous.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

Sedimentation of the neutron rich isotope $^{22}$Ne may be an important source of gravitational energy during the cooling of white dwarf stars. This depends on the diffusion constant for $^{22}$Ne in strongly coupled plasma mixtures. We calculate self-diffusion constants $D_i$ from molecular dynamics simulations of carbon, oxygen, and neon mixtures. We find that $D_i$ in a mixture does not differ greatly from earlier one component plasma results. For strong coupling (coulomb parameter $\Gamma>$ few), $D_i$ has a modest dependence on the charge $Z_i$ of the ion species, $D_i \propto Z_i^{-2/3}$. However $D_i$ depends more strongly on $Z_i$ for weak coupling (smaller $\Gamma$). We conclude that the self-diffusion constant $D_{\rm Ne}$ for $^{22}$Ne in carbon, oxygen, and neon plasma mixtures is accurately known so that uncertainties in $D_{\rm Ne}$ should be unimportant for simulations of white dwarf cooling.

SCUBA-2 is the largest submillimetre array camera in the world and was commissioned on the James Clerk Maxwell Telescope (JCMT) with two arrays towards the end of 2009. A period of shared-risks observing was then completed and the full planned complement of 8 arrays, 4 at 850 microns and 4 at 450 microns, are now installed and ready to be commissioned. SCUBA-2 has 10,240 bolometers, corresponding to a data rate of 8 MB/s when sampled at the nominal rate of 200 Hz. The pipeline produces useful maps in near real time at the telescope and often publication quality maps in the JCMT Science Archive (JSA) hosted at the Canadian Astronomy Data Centre (CADC).

We determine the phase diagram for dense carbon/ oxygen mixtures in White Dwarf (WD) star interiors using molecular dynamics simulations involving liquid and solid phases. Our phase diagram agrees well with predictions from Ogata et al. and Medin and Cumming and gives lower melting temperatures than Segretain et al. Observations of WD crystallization in the globular cluster NGC 6397 by Winget et al. suggest that the melting temperature of WD cores is close to that for pure carbon. If this is true, our phase diagram implies that the central oxygen abundance in these stars is less than about 60%. This constraint, along with assumptions about convection in stellar evolution models, limits the effective S factor for the $^{12}$C($\alpha,\gamma$)$^{16}$O reaction to S_{300} <= 170 keV barns.

This paper describes a new Heterodyne Array Receiver Programme (HARP) and Auto-Correlation Spectral Imaging System (ACSIS) that have recently been installed and commissioned on the James Clerk Maxwell Telescope (JCMT). The 16-element focal-plane array receiver, operating in the submillimetre from 325 to 375 GHz, offers high (three-dimensional) mapping speeds, along with significant improvements over single-detector counterparts in calibration and image quality. Receiver temperatures are $\sim$120 K across the whole band and system temperatures of $\sim$300K are reached routinely under good weather conditions. The system includes a single-sideband filter so these are SSB figures. Used in conjunction with ACSIS, the system can produce large-scale maps rapidly, in one or more frequency settings, at high spatial and spectral resolution. Fully-sampled maps of size 1 square degree can be observed in under 1 hour. The scientific need for array receivers arises from the requirement for programmes to study samples of objects of statistically significant size, in large-scale unbiased surveys of galactic and extra-galactic regions. Along with morphological information, the new spectral imaging system can be used to study the physical and chemical properties of regions of interest. Its three-dimensional imaging capabilities are critical for research into turbulence and dynamics. In addition, HARP/ACSIS will provide highly complementary science programmes to wide-field continuum studies, and produce the essential preparatory work for submillimetre interferometers such as the SMA and ALMA.

Using molecular dynamics simulations, we determine the structure of neutron star crust made of rapid proton capture nucleosynthesis material. We find a regular body centered cubic lattice, even with the large number of impurities that are present. Low charge $Z$ impurities tend to occupy interstitial positions, while high $Z$ impurities tend to occupy substitutional lattice sites. We find strong attractive correlations between low $Z$ impurities that could significantly increase the rate of pycnonuclear (density driven) nuclear reactions. The thermal conductivity is significantly reduced by electron impurity scattering. Our results will be used in future work to study the effects of impurities on mechanical properties such as the shear modulus and breaking strain.

Recently, crust cooling times have been measured for neutron stars after extended outbursts. These observations are very sensitive to the thermal conductivity $\kappa$ of the crust and strongly suggest that $\kappa$ is large. We perform molecular dynamics simulations of the structure of the crust of an accreting neutron star using a complex composition that includes many impurities. The composition comes from simulations of rapid proton capture nucleosynthesys followed by electron captures. We find that the thermal conductivity is reduced by impurity scattering. In addition, we find phase separation. Some impurities with low atomic number $Z$ are concentrated in a subregion of the simulation volume. For our composition, the solid crust must separate into regions of different compositions. This could lead to an asymmetric star with a quadrupole deformation that radiates gravitational waves. Observations of crust cooling can constrain impurity concentrations.

We calculate the shear viscosity $\eta$ and thermal conductivity $\kappa$ of a nuclear pasta phase in neutron star crusts. This involves complex non-spherical shapes. We use semiclassical molecular dynamics simulations involving 40,000 to 100,000 nucleons. The viscosity $\eta$ can be simply expressed in terms of the height $Z^*$ and width $\Delta q$ of the peak in the static structure factor $S_p(q)$. We find that $\eta$ increases somewhat, compared to a lower density phase involving spherical nuclei, because $Z^*$ decreases from form factor and ion screening effects. However, we do not find a dramatic increase in $\eta$ from non-spherical shapes, as may occur in conventional complex fluids.

Fusion reactions in the crust of an accreting neutron star are an important source of heat, and the depth at which these reactions occur is important for determining the temperature profile of the star. Fusion reactions depend strongly on the nuclear charge $Z$. Nuclei with $Z\le 6$ can fuse at low densities in a liquid ocean. However, nuclei with Z=8 or 10 may not burn until higher densities where the crust is solid and electron capture has made the nuclei neutron rich. We calculate the $S$ factor for fusion reactions of neutron rich nuclei including $^{24}$O + $^{24}$O and $^{28}$Ne + $^{28}$Ne. We use a simple barrier penetration model. The $S$ factor could be further enhanced by dynamical effects involving the neutron rich skin. This possible enhancement in $S$ should be studied in the laboratory with neutron rich radioactive beams. We model the structure of the crust with molecular dynamics simulations. We find that the crust of accreting neutron stars may contain micro-crystals or regions of phase separation. Nevertheless, the screening factors that we determine for the enhancement of the rate of thermonuclear reactions are insensitive to these features. Finally, we calculate the rate of thermonuclear $^{24}$O + $^{24}$O fusion and find that $^{24}$O should burn at densities near $10^{11}$ g/cm$^3$. The energy released from this and similar reactions may be important for the temperature profile of the star.

Nucleosynthesis, on the surface of accreting neutron stars, produces a range of chemical elements. We perform molecular dynamics simulations of crystallization to see how this complex composition forms new neutron star crust. We find chemical separation, with the liquid ocean phase greatly enriched in low atomic number elements compared to the solid crust. This phase separation should change many crust properties such as the thermal conductivity and shear modulus. The concentration of carbon, if present, is enriched in the ocean. This may allow unstable thermonuclear burning of the carbon and help explain the ignition of the very energetic explosions known as superbursts.

Neutrinos in core collapse supernovae are likely trapped by neutrino-nucleus elastic scattering. Using molecular dynamics simulations, we calculate neutrino mean free paths and ion-ion correlation functions for heterogeneous plasmas. Mean free paths are systematically shorter in plasmas containing a mixture of ions compared to a plasma composed of a single ion species. This is because neutrinos can scatter from concentration fluctuations. The dynamical response function of a heterogeneous plasma is found to have an extra peak at low energies describing the diffusion of concentration fluctuations. Our exact molecular dynamics results for the static structure factor reduce to the Debye Huckel approximation, but only in the limit of very low momentum transfers.

The nuclear pasta -- a novel state of matter having nucleons arranged in a variety of complex shapes -- is expected to be found in the crust of neutron stars and in core-collapse supernovae at subnuclear densities of about $10^{14}$ g/cm$^3$. Due to frustration, a phenomenon that emerges from the competition between short-range nuclear attraction and long-range Coulomb repulsion, the nuclear pasta displays a preponderance of unique low-energy excitations. These excitations could have a strong impact on many transport properties, such as neutrino propagation through stellar environments. The excitation spectrum of the nuclear pasta is computed via a molecular-dynamics simulation involving up to 100,000 nucleons. The dynamic response of the pasta displays a classical plasma oscillation in the 1-2 MeV region. In addition, substantial strength is found at low energies. Yet this low-energy strength is missing from a simple ion model containing a single-representative heavy nucleus. The low-energy strength observed in the dynamic response of the pasta is likely to be a density wave involving the internal degrees of freedom of the clusters.

Nonuniform neutron-rich matter present in both core-collapse supernovae and neutron-star crusts is described in terms of a semiclassical model that reproduces nuclear-matter properties and includes long-range Coulomb interactions. The neutron-neutron correlation function and the corresponding static structure factor are calculated from molecular dynamics simulations involving 40,000 to 100,000 nucleons. The static structure factor describes coherent neutrino scattering which is expected to dominate the neutrino opacity. At low momentum transfers the static structure factor is found to be small because of ion screening. In contrast, at intermediate momentum transfers the static structure factor displays a large peak due to coherent scattering from all the neutrons in a cluster. This peak moves to higher momentum transfers and decreases in amplitude as the density increases. A large static structure factor at zero momentum transfer, indicative of large density fluctuations during a first-order phase transition, may increase the neutrino opacity. However, no evidence of such an increase has been found. Therefore, it is unlikely that the system undergoes a simple first-order phase transition. It is found that corrections to the commonly used single heavy nucleus approximation first appear at a density of the order of $10^{13}$ g/cm$^3$ and increase rapidly with increasing density. Thus, neutrino opacities are overestimated in the single heavy nucleus approximation relative to the complete molecular dynamics simulations.

A polarimeter has been built for use with the Submillimetre Common-User Bolometer Array (SCUBA), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. SCUBA is the first of a new generation of highly sensitive submillimetre cameras, and the UK/Japan Polarimeter adds a polarimetric imaging/photometry capability in the wavelength range 350 to 2000 microns. Early science results range from measuring the synchrotron polarization of the black hole candidate Sgr A* to mapping magnetic fields inferred from polarized dust emission in Galactic star-forming clouds. We describe the instrument design, performance, observing techniques and data reduction processes, along with an assessment of the current and future scientific capability.