The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

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 report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

LDN 1157, is one of the several clouds situated in the cloud complex, LDN 1147/1158, represents a coma-shaped morphology with a well-collimated bipolar outflow emanating from a Class 0 protostar, LDN 1157-mm. The main goals of this work are (a) to map the inter-cloud magnetic field (ICMF) geometry of the region surrounding LDN 1157 to investigate its relationship with the cloud morphology, with the outflow direction and with the core magnetic field (CMF) geometry inferred from the mm- and sub-mm polarization results from the literature, and (b) to investigate the kinematic structure of the cloud. We carried out R-band polarization observations of the stars projected on the cloud to map the pc-scale magnetic field geometry and made spectroscopic observations of the entire cloud in 12CO, C18O and N2H+ (J=1-0) lines to investigate its kinematic structure. We obtained a distance of 340$\pm$3 pc to the LDN 1147/1158, complex based on the Gaia DR2 parallaxes and proper motion values of the three YSOs associated with the complex. A single filament of $\sim1.2$ pc in length and $\sim0.09$ pc in width is found to run all along the coma-shaped cloud. Based on the relationships between the ICMF, CMF, filament orientations, outflow direction, and the presence of an hour-glass morphology of the magnetic field, it is likely that the magnetic field had played an important role in the star formation process in LDN 1157. Combining the proper motions of the YSOs and the radial velocity of LDN 1147/1158 and another complex LDN 1172/1174 which is situated $\sim2$\dgr~east of it, we found that both the complexes are moving collectively toward the Galactic plane. The filamentary morphology of the east-west segment of LDN 1157 may have formed as a result of mass lost by ablation due to the interaction of the moving cloud with the ambient interstellar medium.

Radio observations of Asymptotic Giant Branch (AGB) star envelopes frequently show some excess emission at frequencies below 100 GHz which cannot be explained by thermal dust emission (hereafter anomalous microwave emission-AME). Moreover, AGB envelopes are a common place where gas molecules condense to form nanoparticles (e.g., polycyclic aromatic hydrocarbons) and large grains. In this paper, we will study whether electric dipole emission from rapidly spinning nanoparticles can reproduce the AME observed toward AGB stars. To properly model the size distribution of nanoparticles in the AGB envelope, we take into account both the increase of nanoparticles due to rotational disruption of large grains spun-up by radiative torques and the decrease of smallest nanoparticles due to rotational disruption driven by stochastic gas collisions. We then perform detailed modeling of microwave emission from rapidly spinning nanoparticles from both C-rich and O-rich AGB envelopes using the grain size distribution constrained by rotational disruption. We find that spinning dust emission is dominant over thermal dust emission at frequencies below 100 GHz. We attempt to fit the observational data of AME using our spinning dust model and demonstrate that spinning dust can reproduce the observed AME in six AGB stars. Finally, we discuss that microwave emission from spinning dust in AGB envelopes could be observed with high-resolution upcoming radio telescopes such as ngVLA and ALMA Band 1. This would be a major leap for understanding AGB envelopes, formation, evolution, and internal structures of dust. Observations would help to distinguish the carrier of AME from comparing C-rich to O-rich stars, because PAHs are formed in C-rich AGB stars while silicates are formed in O-rich stars.

The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S $(J=2-1)$ and N$_2$H$^+$ $(J=1-0)$ with the NRO 45-m telescope, doubling the number of cores where the CS depletion was probed using C$^{34}$S. In most of our targets, the distribution of C$^{34}$S emission shows features that suggest that the CS molecule is generally depleted in the center of the prestellar cores. The radial profile of the CS abundance with respect to H$_2$ directly measured from the CS emission and the Herschel dust emission indicates that the CS molecule is depleted by a factor of $\sim$3 toward the central regions of the cores with respect to their outer regions. The degree of the depletion is found to be even more enhanced by an order of magnitude when the contaminating effect introduced by the presence of CS molecules in the surrounding envelope that lie along the line-of-sight is removed. Except for L1197 which is classified as relatively the least evolved core in our targets based on its observed physical parameters, we found that the remaining four prestellar cores are suffering from significant CS depletion at their central region regardless of the relative difference in their evolutionary status.

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

We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution.

Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper.

We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$\mu$m polarized dust emission observed with the POL-2 instrument at JCMT. We examine the magnetic field geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the central clumps MM1 and MM2, whereas MM3 in the north has field orientations aligned with its major axis. The overall field orientations are uniform at large (POL-2 at 14$\arcsec$ and SHARP at 10$\arcsec$) to small scales (TADPOL at 2.5$\arcsec$ and SMA at 1.5$\arcsec$) in the MM1 and MM2 regions. SHARP/CSO observations in MM3 at 350\,$\mu$m from Tang et al. show a similar trend as seen in our POL-2 observations. TADPOL observations demonstrate a well-defined field geometry in MM1/MM2 consistent with MHD simulations of accreting filaments. We obtained a plane-of-sky magnetic field strength of 470$\pm$190\,$\mu$G, 100$\pm$40\,$\mu$G, and 60$\pm$34\,$\mu$G in the central, northern and southern regions of G34, respectively, using the updated Davis-Chandrasekhar-Fermi relation. The estimated value of field strength, combined with column density and velocity dispersion values available in the literature, suggests G34 to be marginally critical with criticality parameter $\rm \lambda$ values 0.8$\pm$0.4, 1.1$\pm$0.8, and 0.9$\pm$0.5 in the central, northern, and southern regions, respectively. The turbulent motions in G34 are sub-Alfv\'{e}nic with Alfv\'{e}nic Mach numbers of 0.34$\pm$0.13, 0.53$\pm$0.30, and 0.49$\pm$0.26 in the three regions. The observed aligned B-fields in G34.43+0.24 are consistent with theoretical models suggesting that B-fields play an important role in guiding the contraction of the cloud driven by gravity.

L1521F is found to be forming multiple cores and it is cited as an example of the densest core with an embedded VeLLO in a highly dynamical environment. We present the core-scale magnetic fields (B-fields) in the near vicinity of the VeLLO L1521F-IRS using submm polarization measurements at 850$~\mu$m using JCMT POL-2. This is the first attempt to use high-sensitivity observations to map the sub-parsec scale B-fields in a core with a VeLLO. The B-fields are ordered and very well connected to the parsec-scale field geometry seen in our earlier optical polarization observations and the large-scale structure seen in Planck dust polarization. The core scale B-field strength estimated using Davis-Chandrasekhar-Fermi relation is $\rm 330\pm100~\mu$G which is more than ten times of the value we obtained in the envelope (envelope in this paper is "core envelope"). This indicates that B-fields are getting stronger on smaller scales. The magnetic energies are found to be 1 to 2 orders of magnitude higher than non-thermal kinetic energies in the envelope and core. This suggests that magnetic fields are more important than turbulence in the energy budget of L1521F. The mass-to-flux ratio of 2.3$\pm$0.7 suggests that the core is magnetically-supercritical. The degree of polarization is steadily decreasing towards the denser part of the core with a power law slope of -0.86.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

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.

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

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 present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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.

We present the results of a systematic search for molecular outflows in 68 Very Low Luminosity Objects (VeLLOs) from single-dish observations in CO isotopologues which find 16 VeLLOs showing clear outflow signatures in the CO maps. With additional three VeLLOs from the literature, we analyzed the outflow properties for 19 VeLLOs, identifying 15 VeLLOs as proto-Brown Dwarf (BD) candidates and four VeLLOs as likely faint protostar candidates. The proto-BD candidates are found to have a mass accretion rate ($\sim 10^{-8} - 10^{-7}$ $\rm M_{\odot}$ yr$^{-1}$) lower than that of the protostar candidates ($\gtrsim 10^{-6}$ $\rm M_{\odot}$ yr$^{-1}$). Their accretion luminosities are similar to or smaller than their internal luminosities, implying that many proto-BD candidates might have had either small accretion activity in a quiescent manner throughout their lifetime, or be currently exhibiting a relatively higher (or episodic) mass accretion than the past. There are strong trends that outflows of many proto-BDs are less active if they are fainter or have less massive envelopes. The outflow forces and internal luminosities for more than half of the proto-BD candidates seem to follow an evolutionary track of a protostar with its initial envelope mass of $\sim$0.08 $\rm M_{\odot}$, indicating that some BDs may form in less massive dense cores in a way similar to normal stars. But, because there also exists a significant fraction (about 40%) of proto-BDs with much weaker outflow force than expected by the relations for protostars, we should not rule out the possibility of other formation mechanism for the BDs.

How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $\mu$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex.

This paper presents our observational attempts to precisely measure the central mass of a proto-brown dwarf candidate, L328-IRS, in order to investigate whether L328-IRS is in the substellar mass regime. Observations were made for the central region of L328-IRS with the dust continuum and CO isotopologue line emission at ALMA band 6, discovering the detailed outflow activities and a deconvolved disk structure of a size of $\sim$87 AU $\times$ 37 AU. We investigated the rotational velocities as a function of the disk radius, finding that its motions between 130 AU and 60 AU are partially fitted with a Keplerian orbit by a stellar object of $\sim0.30 ~\rm\,M_\odot$, while the motions within 60 AU do not follow any Keplerian orbit at all. This makes it difficult to lead a reliable estimation of the mass of L328-IRS. Nonetheless, our ALMA observations were useful enough to well constrain the inclination angle of the outflow cavity of L328-IRS as $\sim 66^{\circ}$, enabling us to better determine the mass accretion rate of $\sim 8.9\times 10^{-7} \rm\,M_\odot~yr^{-1}$. From assumptions that the internal luminosity of L328-IRS is mostly due to this mass accretion process in the disk, or that L328-IRS has mostly accumulated the mass through this constant accretion rate during its outflow activity, its mass was estimated to be $\sim 0.012 - 0.023\rm\,M_\odot$, suggesting L328-IRS to be a substellar object. However, we leave our identification of L328-IRS as a proto-brown dwarf to be tentative because of various uncertainties especially regarding the mass accretion rate.

The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

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.

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

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 made R-band polarization measurements of 234 stars towards the direction of MBM 33-39 cloud complex. The distance of MBM 33-39 complex was determined as $120\pm10$ pc using polarization results and near-infrared photometry from the 2MASS survey. The magnetic field geometry of the individual cloud inferred from our polarimetric results reveals that the field lines are in general consistent with the global magnetic field geometry of the region obtained from the previous studies. This implies that the clouds in the complex are permeated by the interstellar magnetic field. Multi-wavelength polarization measurements of a few stars projected on the complex suggest that the size of the dust grains in these clouds is similar to those found in the normal interstellar medium of the Milky Way. We studied a possible formation scenario of MBM 33-39 complex by combining the polarization results from our study and from the literature and by identifying the distribution of ionized, atomic and molecular (dust) components of material in the region.

Bright-rimmed clouds form on the edges of H II regions affected by the high energy radiation from a central ionizing source. The UV radiation from the ionizing source results in compression and ionization causing either cloud disruption or further star formation. In this work, we present R-band polarization measurements towards four bright-rimmed clouds, IC1396A, BRC 37, BRC 38, and BRC 39, located in the different directions of the H II region, Sh2-131, in order to map magnetic fields (B-fields) in the plane of the sky. These BRCs are illuminated by the O star HD206267 and present a range of projected on sky geometries. This provides an opportunity to understand the magnetized evolution of BRCs. The B-field geometries of the clouds deduced from the polarization data, after correction for foreground ISM contamination, are seen to be connected to the ambient B-fields on the large scale. They seem to play an important role in shaping the cloud IC1396A and BRC 37. BRCs 38 and 39 show a broader and snubber head morphology possibly due to the B-fields being aligned with incoming radiation as explained in the simulations. A good general agreement is noted on comparing our observational results with the simulations supporting the importance of B-fields in BRC evolution. This work is the first step towards systematic mapping the B-fields morphology in multiple BRCs in an expanding H II region, extending the work presented by Soam et al. (2017b).

Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of $T_{\rm d}$ = 11-14 K, and H$_{2}$ column densities of $N_{\rm H_{2}}$ = 0.36-2.5$\times10^{22}$ cm$^{-2}$. We found that not all PGCCs contain prestellar objects. In general, the dense cores in PGCCs are usually at their earliest evolutionary stages. All the dense cores have non-thermal velocity dispersions larger than the thermal velocity dispersions from molecular line data, suggesting that the dense cores may be turbulence-dominated. We have calculated the virial parameter $\alpha$ and found that 14 of the dense cores have $\alpha$ $<$ 2, while 16 of the dense cores have $\alpha$ $>$ 2. This suggests that some of the dense cores are not bound in the absence of external pressure and magnetic fields. The column density profiles of dense cores were fitted. The sizes of the flat regions and core radii decrease with the evolution of dense cores. CO depletion was found to occur in all the dense cores, but is more significant in prestellar core candidates than in protostellar or starless cores. The protostellar cores inside the PGCCs are still at a very early evolutionary stage, sharing similar physical and chemical properties with the prestellar core candidates.

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

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.

Bright-rimmed clouds (BRCs) are formed at the periphery of H$~$II regions as the radiation from the central star interacts with dense gas. The ionization and resulting compression of the clouds may lead to cloud disruption causing secondary star formation depending on the stellar and gas parameters. Here we use R-band polarimetry to probe the plane-of-the sky magnetic field in the two near-by BRCs IC\,59 and IC\,63. Both nebulae are illuminated by $\gamma$ Cas with the direction of ionizing radiation being orientated parallel or perpendicular to the local magnetic field, allowing us to probe the importance of magnetic field pressure in the evolution of BRCs. Because of the proximity of the system ($\sim$200pc) we have acquired a substantial sample of over 500 polarization measurements for stars background to the nebulae. On large scales, the magnetic field geometries of both clouds are anchored to the ambient magnetic field. For IC 63, the magnetic field is aligned parallel to the head-tail morphology of the main condensation, with convex morphology relative to the direction of the ionizing radiation. We estimate the plane of the sky magnetic field strength in IC\,63 to be $\sim90\mu$G. In IC\,59, the projected magnetic field follows the M shape morphology of the cloud. Here, field lines present a concave shape with respect to the direction of the ionizing radiation from $\gamma$ Cas. Comparing our observations to published theoretical models we find good general agreement, supporting the importance of magnetic fields in BRC evolution.

We present the R-band polarimetric results towards two nebulae L1415 and L1389 containing low luminosity stars. Aim of this study is to understand the role played by magnetic fields in formation of low luminosity objects. Linear polarization arise due to dichroism of the background starlight projected on the cloud providing the plane-of-the sky magnetic field orientation. The offsets between mean magnetic field directions obtained towards L1415 and L1389 and the projected outflow axes are found to be 35$^{\circ}$ and 12$^{\circ}$, respectively. The offset between cloud minor axes and mean envelope magnetic field direction in L1415 and L1389 are 50$^{\circ}$ and 87$^{\circ}$, respectively. To estimate the magnetic field strength by using the updated Chandrasekhar-Fermi relation, we obtained the $^{12}$CO(J=1-0) line velocity dispersion value towards L1415 cloud using the TRAO single dish observations. The values of B$_{pos}$ in L1415 and L1389 are found to be 28$~\mu$G and 149$~\mu$G using CF technique and 23$~\mu$G and 140$~\mu$G using structure function analysis, respectively. The values of B$_{pos}$ in these clouds are found to be consistent using both the techniques. By combining the present results with those obtained from our previous study of magnetic fields in cores with VeLLOs, we attempt to improve the sample of cores with low luminosity protostars and bridge the gap between the understanding of importance of magnetic fields in cores with VeLLOs and low luminosity protostars. The Results of this work and that of our previous work show that the outflow directions are aligned with envelope magnetic fields of the clouds.

We carried out optical polarimetry of an isolated cloud, Gal 110-13, to map the plane-of-the-sky magnetic field geometry. The main aim of the study is to understand the most plausible mechanism responsible for the unusual cometary shape of the cloud in the context of its magnetic field geometry. When unpolarized starlight passes through the intervening interstellar dust grains that are aligned with their short axes parallel to the local magnetic field, it gets linearly polarized. The plane-of-the-sky magnetic field component can therefore be traced by doing polarization measurements of background stars projected on clouds. Because the light in the optical wavelength range is most efficiently polarized by the dust grains typically found in the outer layers of the molecular clouds, optical polarimetry enables us to trace the magnetic field geometry of the outer layers of the clouds. We made R-band polarization measurements of 207 stars in the direction of Gal 110-13. The distance of Gal 110-13 was determined as $\sim450\pm80$ pc using our polarization and 2MASS near-infrared data. The foreground interstellar contribution was removed from the observed polarization values by observing a number of stars located in the vicinity of Gal 110-13 which has Hipparcos parallax measurements. The plane-of-the-sky magnetic field lines are found to be well ordered and aligned with the elongated structure of Gal 110-13. Using structure function analysis, we estimated the strength of the plane-of-the-sky component of the magnetic field as $\sim25\mu$G. Based on our results and comparing them with those from simulations, we conclude that compression by the ionization fronts from 10 Lac is the most plausible cause of the comet-like morphology of Gal 110-13 and of the initiation of subsequent star formation.

[Abridged] We carried out optical polarimetry of five dense cores, (IRAM 04191, L1521F, L328, L673-7, and L1014) which are found to harbour VeLLO. This study was conducted mainly to understand the role played by the magnetic field in the formation of very low and substellar mass range objects using optical polarisation. The angular offsets between the envelope magnetic field direction (inferred from optical polarisation measurements) and the outflow position angles from the VeLLOs in IRAM 04191, L1521F, L328, L673-7, and L1014 are found to be 84$^\circ$, 53$^\circ$, 24$^\circ$, 08$^\circ$, and 15$^\circ$, respectively. The mean value of the offsets for all the five clouds is $\sim37^\circ$. If we exclude IRAM 04191, the mean value reduces to become $\sim25^\circ$. In IRAM 04191, the offset between the projected envelope and the inner magnetic field (inferred from the submillimetre data from SCUPOL) is found to be $\sim68^\circ$. The inner magnetic field, however, is found to be nearly aligned with the projected position angles of the minor axis, the rotation axis of the cloud, and the outflow from the IRAM 04191-IRS. We discuss a possible explanation for the nearly perpendicular orientation between the envelope and core scale magnetic fields in IRAM04191. The angular offset between the envelope magnetic field direction and the minor axis of IRAM 04191, L1521F, L673-7, and L1014 are 82$^\circ$, 60$^\circ$, 47$^\circ$, and 55$^\circ$, respectively. The mean value of the offsets between the envelope magnetic field and the minor axis position angles for the four cores is found to be $\sim60^\circ$. The results obtained from our study on the limited sample of five cores with VeLLOs show that the outflows in three of them tend to nearly align with the envelope magnetic field.

We present results of our $R-$band polarimetry of a cometary globule, LBN 437 (Gal96-15, $\ell$ $=$ 96$\degree$, \textit{b} $=-15\degree$), to study magnetic field geometry of the cloud. We estimated a distance of $360\pm65$ pc to LBN 437 (also one additional cloud, CB 238) using near-IR photometric method. Foreground contribution to the observed polarisation values was subtracted by making polarimetric observations of stars that are located in the direction of the cloud and with known distances from the Hipparcos parallax measurements. The magnetic field geometry of LBN 437 is found to follow the curved shape of the globule head. This could be due to the drag that the magnetic field lines could have experienced because of the ionisation radiation from the same exciting source that caused the cometary shape of the cloud. The orientation of the outflow from the Herbig A4e star, LkH$\alpha$ 233 (or V375 Lac), located at the head of LBN 437, is found to be parallel to both the initial (prior to the ionising source was turned on) ambient magnetic field (inferred from a star HD 214243 located just in front of the cloud) and the Galactic plane.