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.

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 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 $Herschel$ Space telescope carried out an unprecedented survey of nearby stars for debris disks. The dust present in these debris disks scatters and polarizes stellar light in the visible part of the spectrum. We explore what can be learned with aperture polarimetry and detailed radiative transfer modelling about stellar systems with debris disks. We present a polarimetric survey, with measurements from the literature, of candidate stars observed by DEBRIS and DUNES $Herschel$ surveys. We perform a statistical analysis of the polarimetric data with the detection of far-infrared excess by $Herschel$ and $Spitzer$ with a sample of 223 stars. Monte Carlo simulations were performed to determine the effects of various model parameters on the polarization level and find the mass required for detection with current instruments. Eighteen stars were detected with a polarization $0.01 \le P \lesssim 0.1$ per cent and $\ge3\sigma_P$, but only two of them have a debris disk. No statistically significant difference is found between the different groups of stars, with, without, and unknown status for far-infrared excess, and presence of polarization. The simulations show that the integrated polarization is rather small, usually $< 0.01$ per cent for typical masses detected by their far-infrared excess for hot and most warm disks. Masses observed in cold disks can produce polarization levels above $0.01$ per cent since there is usually more dust in them than in closer disks. We list five factors which can explain the observed low-polarization detection rate. Observations with high-precision polarimeters should lead to additional constraints on models of unresolved debris disks.

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.

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.

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.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

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

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.

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.

Thermal emission from cold dust grains in giant molecular clouds can be used to probe the physical properties, such as density, temperature and emissivity in star-forming regions. We present the SCUBA-2 shared-risk observations at 450 $\mu$m and 850 $\mu$m of the Orion A molecular cloud complex taken at the James Clerk Maxwell Telescope (JCMT). Previous studies showed that molecular emission lines can contribute significantly to the measured fluxes in those continuum bands. We use the HARP $^{12}$CO J=3-2 integrated intensity map for Orion A in order to evaluate the molecular line contamination and its effects on the SCUBA-2 maps. With the corrected fluxes, we have obtained a new spectral index $\alpha$ map for the thermal emission of dust in the well-known integral-shaped filament. Furthermore, we compare a sample of 33 sources, selected over the Orion A molecular cloud complex for their high $^{12}$CO J=3-2 line contamination, to 27 previously identified clumps in OMC-4. This allows us to quantify the effect of line contamination on the ratio of 850 $\mu$m to 450 $\mu$m flux densities and how it modifies the deduced spectral index of emissivity $\beta$ for the dust grains. We also show that at least one Spitzer-identified protostellar core in OMC-5 has a $^{12}$CO J=3-2 contamination level of 16 %. Furthermore, we find the strongest contamination level (44 %) towards a young star with disk near OMC-2. This work is part of the JCMT Gould Belt Legacy Survey.

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.

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.

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.

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.

We present near infrared polarimetry data of background stars shining through a selection of starless cores taken in the $K$ band, probing visual extinctions up to $A_{V} \sim 48$. We find that $P_K/{\tau _K}$ continues to decline with increasing $A_{V}$ with a power law slope of roughly -0.5. Examination of published submillimeter (submm) polarimetry of starless cores suggests that by $A_{V} \gtrsim 20$ the slope for $P$ vs. $\tau$ becomes $\sim -1$, indicating no grain alignment at greater optical depths. Combining these two data sets, we find good evidence that, in the absence of a central illuminating source, the dust grains in dense molecular cloud cores with no internal radiation source cease to become aligned with the local magnetic field at optical depths greater than $A_V \sim 20$. A simple model relating the alignment efficiency to the optical depth into the cloud reproduces the observations well.

New visible and K-band polarization measurements on stars surrounding molecular clouds in Orion A and stars in the BN vicinity are presented. Our results confirm that magnetic fields located inside the Orion A molecular clouds and in their close neighborhood are spatially connected. On and around the BN object, we measured the angular offsets between the K-band polarization data and available submm data. We find high values of the polarization degree, P_{K}, and of the optical depth, \tau_{K}, close to an angular offset position of 90^{\circ} whereas lower values of P_{K} and \tau_{K} are observed for smaller angular offsets. We interpret these results as evidence for the presence of various magnetic field components toward lines of sight in the vicinity of BN. On a larger scale, we measured the distribution of angular offsets between available H-band polarization data and the same submm data set. Here we find an increase of <P_{H}> with angular offset which we interpret as a rotation of the magnetic field by \lesssim 60^{\circ}. This trend generalizes previous results on small scale toward and around lines of sight to BN and is consistent with a twist of the magnetic field on a larger scale towards OMC-1. A comparison of our results with several other studies suggests that a two-component magnetic field, maybe helical, could be wrapping the OMC-1 filament.

Submillimetre continuum radiation allows us to probe cold objects, particularly the earliest, dusty phases of star formation, high-redshift galaxies and circumstellar disks. The submillimetre window gives a unique view of the physical and dynamical conditions in the neutral and molecular interstellar medium. In the next decade a combination of wide-field surveys with single-dish telescopes and targeted follow-up with ALMA and other facilities should enable rapid progress in answering questions about the origins of planetary systems, stars and galaxies.

We present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

The SCUBA polarized 850 microns thermal emission data of the region OMC-2 in Orion A are added to and homogeneously reduced with data already available in the region OMC-3. The data set shows that OMC-2 is a region generally less polarized than OMC-3. Where coincident, most of the 850 microns polarization pattern is similar to that measured in 350 microns polarization data. Only 850 microns polarimetry data have been obtained in and around MMS7, FIR1 & FIR2, and in the region south of FIR6. A realignment of the polarization vectors with the filament can be seen near FIR1 in the region south of OMC-3. An analysis shows that the energy injected by CO outflows and H2 jets associated to OMC-2 and OMC-3 does not appear to alter the polarization patterns at a scale of the 14'' resolution beam. A second order structure function analysis of the polarization position angles shows that OMC-2 is a more turbulent region than OMC-3. OMC-3 appears to be a clear case of a magnetically dominated region with respect to the turbulence. However for OMC-2 it is not clear that this is the case. A more in-depth analysis of five regions displayed along OMC-2/3 indicates a decrease of the mean polarization degree and an increase of the turbulent angular dispersion from north to south. A statistical analysis suggests the presence of two depolarization regimes in our maps. One regime including the effects of the cores, the other one excluding it.

The properties of polarization in scattered light by aligned ellipsoidal grains are investigated with the Fredholm integral equation method (FIM) and the T-matrix method (Tmat), and the results are applied to the observed circular polarization in OMC1. We assume that the grains are composed of silicates and ellipsoidal (oblate, prolate, or tri-axial ellipsoid) in shape with a typical axial ratio of 2:1. The angular dependence of circular polarization p_c on directions of incident and scattered light is investigated with spherical harmonics and associated Legendre polynomials. The degree of circular polarization p_c also depends on the Rayleigh reduction factor R which is a measure of imperfect alignment. We find that p_c is approximately proportional to R for grains with |m|x_{eq} < 3 - 5, where x_{eq} is the dimensionless size parameter and m is the refractive index of the grain. Models that include those grains can explain the observed large circular polarization in the near infrared, ~15%, in the south-east region of the BN object (SEBN) in OMC1, if the directions of incidence and scattering of light is optimal, and if grain alignment is strong, i.e. R > 0.5. Such a strong alignment cannot be explained by the Davis-Greenstein mechanism; we prefer instead an alternative mechanism driven by radiative torques. If the grains are mixed with silicates and ice, the degree of circular polarization p_c decreases in the 3 micron ice feature, while that of linear polarization increases. This wavelength dependence is different from that predicted in a process of dichroic extinction.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

New visible polarization data combined with existing IR and FIR polarization data are used to study how the magnetic field threading the filamentary molecular cloud GF 9 connects to larger structures in its general environment. We find that when both visible and NIR polarization data are plotted as a function of extinction, there is no evidence for a plateau or a saturation effect in the polarization at Av ~ 1.3 as seen in dark clouds in Taurus. This lack of saturation effect suggests that even in the denser parts of GF 9 we are still probing the magnetic field. The visible polarization is smooth and has a well-defined orientation. The IR data are also well defined but with a different direction, and the FIR data in the core region are well defined and with yet another direction, but are randomly distributed in the filament region. On the scale of a few times the mean radial dimension of the molecular cloud, it is as if the magnetic field were `blind' to the spatial distribution of the filaments while on smaller scales within the cloud, in the core region near the IRAS point source PSC 20503+6006, polarimetry shows a rotation of the magnetic field lines in these denser phases. Hence, in spite of the fact that the spatial resolution is not the same in the visible/NIR and in the FIR data, all the data put together indicate that the field direction changes with the spatial scale. Finally, the Chandrasekhar and Fermi method is used to evaluate the magnetic field strength, indicating that the core region is approximately magnetically critical. A global interpretation of the results is that in the core region an original poloidal field could have been twisted by a rotating elongated (core+envelope) structure. There is no evidence for turbulence and ambipolar diffusion does not seem to be effective at the present time.

We present a new continuum 3D radiative transfer code, MCFOST, based on a Monte-Carlo method. MCFOST can be used to calculate (i) monochromatic images in scattered light and/or thermal emission, (ii) polarisation maps, (iii) interferometric visibilities, (iv) spectral energy distributions and (v) dust temperature distributions of protoplanetary disks. Several improvements to the standard Monte Carlo method are implemented in MCFOST to increase efficiency and reduce convergence time, including wavelength distribution adjustments, mean intensity calculations and an adaptive sampling of the radiation field. The reliability and efficiency of the code are tested against a previously defined benchmark, using a 2D disk configuration. No significant difference (no more than 10%, and generally much less) is found between the temperatures and SEDs calculated by MCFOST and by other codes included in the benchmark. A study of the lowest disk mass detectable by Spitzer, around young stars, is presented and the colours of ``representative'' parametric disks are compared to recent IRAC and MIPS Spitzer colours of solar-like young stars located in nearby star forming regions.

An implicit method for radiative transfer in SPH is described. The diffusion approximation is used, and the hydrodynamic calculations are performed by a fully three--dimensional SPH code. Instead of the energy equation of state for an ideal gas, various energy states and the dissociation of hydrogen molecules are considered in the energy calculation for a more realistic temperature and pressure determination. In order to test the implicit code, we have performed non--isothermal collapse simulations of a centrally condensed cloud, and have compared our results with those of finite difference calculations performed by MB93. The results produced by the two completely different numerical methods agree well with each other.

SCUBA-2, which replaces SCUBA (the Submillimeter Common User Bolometer Array) on the James Clerk Maxwell Telescope (JCMT) in 2006, is a large-format bolometer array for submillimeter astronomy. Unlike previous detectors which have used discrete bolometers, SCUBA-2 has two dc-coupled, monolithic, filled arrays with a total of ~10,000 bolometers. It will offer simultaneous imaging of a 50 sq-arcmin field of view at wavelengths of 850 and 450 microns. SCUBA-2 is expected to have a huge impact on the study of galaxy formation and evolution in the early Universe as well as star and planet formation in our own Galaxy. Mapping the sky to the same S/N up to 1000 times faster than SCUBA, it will also act as a pathfinder for the new submillimeter interferometers such as ALMA. SCUBA-2's absorber-coupled pixels use superconducting transition edge sensors operating at 120 mK for performance limited by the sky background photon noise. The monolithic silicon detector arrays are deep-etched by the Bosch process to isolate the pixels on silicon nitride membranes. Electrical connections are made through indium bump bonds to a SQUID time-domain multiplexer (MUX). We give an overview of the SCUBA-2 system and an update on its status, and describe some of the technological innovations that make this unique instrument possible.

We present simultaneous UBV polarimetric and photometric observations of the pre-main-sequence binary AK Sco, obtained over 12 nights, slightly less than the orbital period of 13.6 days. The polarization is a sum of interstellar and intrinsic polarization, with a significant intrinsic polarization of 1% at 5250A, indicating the presence of circumstellar matter distributed in an asymmetric geometry. The polarization and its position angle are clearly variable on time scales of hours and nights, in all 3 wavelengths, with a behavior related to the orbital motion. The variations have the highest amplitudes seen so far for pre-main-sequence binaries (~1%, ~30deg) and are sinusoidal with periods similar to the orbital period and half of it. The polarization variations are generally correlated with the photometric ones: when the star gets fainter, it also gets redder and its polarization increases. The color-magnitude diagram B-V, V exhibits a ratio of total to selective absorption R=4.3 higher than in normal interstellar clouds (R=3.1). The interpretation of the simultaneous photometric and polarimetric observations is that a cloud of circumstellar matter passes in front of the star, decreasing the amount of direct, unpolarized light, and hence increasing the contribution of scattered (blue) light. We show that the large amplitude of the polarization variations can not be reproduced with a single scattering model and axially symmetric circumbinary or circumstellar disks.

We use extensive 350 um polarimetry and continuum maps obtained with Hertz and SHARC II along with HCN and HCO+ spectroscopic data to trace the orientation of the magnetic field in the Orion A star-forming region. Using the polarimetry data, we find that the direction of the projection of the magnetic field in the plane of the sky relative to the orientation of the integral-shaped filament varies considerably as one moves from north to south. While in IRAS 05327-0457 and OMC-3 MMS 1-6 the projection of the field is primarily perpendicular to the filament it becomes better aligned with it at OMC-3 MMS 8-9 and well aligned with it at OMC-2 FIR 6. The OMC-2 FIR 4 cloud, located between the last two, is a peculiar object where we find almost no polarization. The projected angle of the field is more complicated in OMC-1 where it exhibits smooth variations in its orientation across the face of this massive complex. By combining the polarimetry and spectroscopic data we were able to measure a set of average values for the inclination angle of the magnetic field relative to the line of sight. We find that the field is oriented quite close to the plane of the sky in most places. More precisely, the inclination of the magnetic field is ~73 deg. around OMC-3 MMS 6, ~74 deg. at OMC-3 MMS 8-9, ~80 deg. at OMC-2 FIR 4, ~65 deg. in the northeastern part of OMC-1, and ~49 deg. in the Bar. We also present polarimetry data for the OMC-4 region located some 13 arcminutes south of OMC-1.

We demonstrate that the combination of Zeeman, polarimetry and ion-to-neutral molecular line width ratio measurements permits the determination of the magnitude and orientation of the magnetic field in the weakly ionized parts of molecular clouds. Zeeman measurements provide the strength of the magnetic field along the line of sight, polarimetry measurements give the field orientation in the plane of the sky and the ion-to-neutral molecular line width ratio determines the angle between the magnetic field and the line of sight. We apply the technique to the M17 star-forming region using a HERTZ 350 um polarimetry map and HCO+-to-HCN molecular line width ratios to provide the first three-dimensional view of the magnetic field in M17.

The magnetic field is believed to play a crucial role in the process of star formation. From the support it provides during the initial collapse of molecular clouds to the creation of strong collimated jets responsible for large mass losses, current theories predict its importance in many different stages during the formation of stars. Here we report on observational evidence which tests one aspect that can be inferred from these theories: the alignment between the local magnetic field and collimated bipolar outflows in such environments. There is good evidence of an alignment in the case of NGC 2071.

We present further observational evidence in support of our earlier proposal (Houde et al. 2000) for detecting the presence of the magnetic field in molecular clouds by comparing spectra of molecular ions with those of neutral molecules. The ion lines tend to be narrower and do not show the wings due to flows, when the magnetic field is sufficiently strong. We obtained spectra for the optically thin lines of the H13CN and H13CO+ species in a sample of ten molecular clouds and found the results to be in agreement with our previous observations of the main isotopic species, HCN and HCO+, made in OMC-1, OMC-2, OMC-3 and DR21OH, thus eliminating the possibility of optical depth effects playing a role in the ion line narrowing. HCS+ was also detected in four of these star forming regions. We also discuss previously published results by (Benson et al. 1998) of N2H+ detections in a large sample of dark clouds. We show that the similarity in line widths between ion and neutral species in their sample is consistent with the relatively small amount of turbulence and other flows observed in these clouds.

Observations of the effect of the magnetic field on its environment are usually achieved with techniques which rely on the interaction with the spin of the particles under study. Because of the relative weakness of this effect, extraction of the field characteristics proves to be a most challenging task. We take a totally different approach to the problem and show that the manifestation of the magnetic field can be directly observed by means of a comparison of the spectra of molecular ions with those of neutral molecules. This takes advantage of the strong cyclotron interaction between the ions and the field, but requires the presence of flows or turbulent motion in the gas. We compare our theory to data obtained on the OMC-1, OMC-2, OMC-3 and DR21(OH) molecular clouds.

We present polarimetric observations of 7 pre-main-sequence (PMS) spectroscopic binaries located in the rho Oph and Upper Sco star forming regions (SFRs). The average observed polarizations at 7660A are between 0.5% and 3.5%. After estimates of the interstellar polarization are removed, all binaries have an intrinsic polarization above 0.4%. Two binaries, NTTS162814-2427 and NTTS162819-2423S, present high levels of intrinsic polarization between 1.5% and 2.1%. All 7 PMS binaries have a statistically variable or possibly variable polarization. Combining these results with our previous sample of binaries located in the Tau, Aur and Ori SFRs, 68% of the binaries have an intrinsic polarization above 0.5%, and 90% of the binaries are polarimetrically variable or possibly variable. NTTS160814-1857, NTTS162814-2427, and NTTS162819-2423S are clearly polarimetrically variable. The first two also exhibit phase-locked variations over ~10 and ~40 orbits respectively. NTTS160905-1859 shows periodic variations that are not phased-locked and only present for short intervals of time. The amplitudes of the variations reach a few tenths of a percent. The high-eccentricity system NTTS162814-2427 shows single-periodic variations, in agreement with our previous numerical simulations. Non-periodic events introduce stochastic noise that partially masks the periodic variations and prevents the Brown, McLean, & Emslie (1978) formalism from finding a reasonable estimate of the inclination.

We present polarimetric observations of 14 pre-main-sequence (PMS) binaries located in the Taurus, Auriga, and Orion star forming regions. The majority of the average observed polarizations are below 0.5%, and none are above 0.9%. After removal of estimates of the interstellar polarization, about half the binaries have an intrinsic polarization above 0.5%, even though most of them do not present other evidences for the presence of circumstellar dust. Various tests reveal that 77% of the PMS binaries have or possibly have a variable polarization. LkCa3, Par1540, and Par2494 present detectable periodic and phase-locked variations. The periodic polarimetric variations are noisier and of a lesser amplitude (~0.1%) than for other types of binaries, such as hot stars. This could be due to stochastic events that produce deviations in the average polarization, a non-favorable geometry (circumbinary envelope), or the nature of the scatterers (dust grains are less efficient polarizers than electrons). Par1540 is a Weak-line TTauri Star, but nonetheless has enough dust in its environment to produce detectable levels of polarization and variations. A fourth interesting case is W134, which displays rapid changes in polarization that could be due to eclipses. We compare the observations with some of our numerical simulations, and also show that an analysis of the periodic polarimetric variations with the Brown, McLean, & Emslie (BME) formalism to find the orbital inclination is for the moment premature: non-periodic events introduce stochastic noise that partially masks the periodic low-amplitude variations and prevents the BME formalism from finding a reasonable estimate of the orbital inclination.

We present polarimetric observations of a massive pre-main sequence short-period binary star of the Herbig Ae/Be type, MWC 1080. The mean polarization at 7660 A is 1.60% at 81.6 deg, or 0.6% at 139 deg if an estimate of the interstellar polarization is subtracted. The intrinsic polarization points to an asymmetric geometry of the circumstellar or circumbinary environment while the 139 deg intrinsic position angle traces the axis of symmetry of the system and is perpendicular to the position angle of the outflow cavity. The polarization and its position angle are clearly variable, at all wavelengths, and on time scales of hours, days, months, and years. Stochastic variability is accompanied by periodic variations caused by the orbital motion of the stars in their dusty environment. These periodic polarimetric variations are the first phased-locked ones detected for a pre-main sequence binary. The variations are not simply double-periodic (seen twice per orbit) but include single-periodic (seen once per orbit) and higher-order variations. The presence of single-periodic variations could be due to non equal mass stars, the presence of dust grains, an asymmetric configuration of the circumstellar or circumbinary material, or the eccentricity of the orbit. MWC 1080 is an eclipsing binary with primary and secondary eclipses occurring at phases 0.0 and 0.55. The signatures of the eclipses are seen in the polarimetric observations.

We present numerical simulations of the periodic polarimetric variations produced by a binary star placed at the center of an empty spherical cavity inside a circumbinary ellipsoidal and optically thin envelope made of dust grains. Mie single-scattering is considered along with pre- and post-scattering extinction factors which produce a time-varying optical depth and affect the morphology of the periodic variations. We are interested in the effects that various parameters will have on the average polarization, the amplitude of the polarimetric variations, and the morphology of the variability. We show that the absolute amplitudes of the variations are smaller for Mie scattering than for Thomson scattering. Among the four grain types that we have studied, the highest polarizations are produced by grains with sizes in the range 0.1-0.2 micron. In general, the variations are seen twice per orbit. In some cases, because spherical dust grains have an asymmetric scattering function, the polarimetric curves produced also show variations seen once per orbit. Circumstellar disks produce polarimetric variations of greater amplitude than circumbinary envelopes. Another goal of these simulations is to see if the 1978 BME (Brown, McLean, & Emslie, ApJ, 68, 415) formalism, which uses a Fourier analysis of the polarimetric variations to find the orbital inclination for Thomson-scattering envelopes, can still be used for Mie scattering. We find that this is the case, if the amplitude of the variations is sufficient and the true inclinations is i_true > 45 deg. For eccentric orbits, the first-order coefficients of the Fourier fit, instead of second-order ones, can be used to find almost all inclinations.

We present a study of the scaling properties of bipolar outflows, generic to the circumstellar environment of young stellar objects, and which are generally believed to be magnetohydrodynamically driven. Our work consists of a multifractal analysis of near-infrared light scattered by dust grains in the lobes of seven bipolar nebulae, namely GGD 18, V380 Orionis, LkH\alpha 101/NGC 1579, LkH\alpha 233, PV Cephei, V645 Cygni (GL 2789), and V633 Cassiopeiae (LkH\alpha 198). We argue that these objects pertain to a multifractal universality class, with the corresponding universal parameters averaged over the ensemble yielding alpha= 1.96 +/- 0.02, C_1=0.04 +/- 0.02, and H=0.7 +/- 0.2; these results might suggest that the dynamics of the observed objects are similar. We then proceed to investigate the existence of anisotropy in the scaling of GGD 18 using the Generalized Scale Invariance (GSI) formalism. It is found that differential rotation and stratification are involved in the outflow mechanism. The stratification is believed to be dynamical in nature i.e., gravity is apparently not the dominant stratifying force.

We present a study of the polarizing power of the dust in cold dense regions (dark clouds) compared to that of dust in the general interstellar medium (ISM). Our study uses new polarimetric, optical, and spectral classification data for 36 stars to carefully study the relation between polarization percentage (p) and extinction (A_V) in the Taurus dark cloud complex. We find two trends in our p-A_V study: (1) stars background to the warm ISM show an increase in p with A_V; and (2) the percentage of polarization of stars background to cold dark clouds does not increase with extinction. We detect a break in the p-A_V relation at an extinction 1.3 +/- 0.2 mag, which we expect corresponds to a set of conditions where the polarizing power of the dust associated with the Taurus dark clouds drops precipitously. This breakpoint places important restrictions on the use of polarimetry in studying interstellar magnetic fields.