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 mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

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.

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 report the first high-resolution (sub-arcminute) large-scale mapping $^{12}$CO and $^{13}$CO observations of the molecular clouds associated with the giant outer Galaxy HII region CTB~102 (KR 1). These observations were made using a newly commissioned receiver system on the 13.7-m radio telescope at the Taeduk Radio Astronomy Observatory. Our observations show that the molecular clouds have a spatial extent of $60 \times 35$ pc and a total mass of $10^{4.8} - 10^{5.0}$ M$_\odot$. Infrared data from WISE and 2MASS were used to identify and classify the YSO population associated with ongoing star formation activity within the molecular clouds. We directly detect 18 class I/class II YSOs and six transition disk objects. Moving away from the HII region, there is an age/class gradient consistent with sequential star formation. The infrared and molecular-line data were combined to estimate the star formation efficiency (SFE) of the entire cloud as well as the SFE for various sub-regions of the cloud. We find that the overall SFE is between $\sim5 - 10$%, consistent with previous observations of giant molecular clouds. One of the sub-regions, region 1a, is a clear outlier, with a SFE of 17 $-$ 35% on a 5 pc spatial scale. This high SFE is more typical for much smaller (sub-pc scale) star-forming cores, and we think region 1a is likely an embedded massive protocluster.

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

The discrepancy between the theoretical prediction of primordial lithium abundances and astronomical observations is called the Lithium Problem. We find that extra contributions from non-thermal hydrogen and helium during Big Bang nucleosynthesis can explain the discrepancy, for both Li-7 and Li-6, and will change the deuterium abundance only little. The allowed parameter space of such an amount of non-thermal particles and the energy range is shown. The hypothesis is stable regardless of the cross-section uncertainty of relevant reactions and the explicit shape of the energy spectrum.

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

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J = 1 - 0 lines, and dense cores in L1641 were observed in the 2.0-0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of 35 km/s and a precession period of 1600 years. The CO outflow length gives a dynamical timescale of 6300 years, and the protostar must be extremely young. The inferred binary separation of 6-70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of beta = 0.3-2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of beta = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices.

During the protostellar phase of stellar evolution, accretion onto the star is expected to be variable, but this suspected variability has been difficult to detect because protostars are deeply embedded. In this paper, we describe a sub-mm luminosity burst of the Class I protostar EC 53 in Serpens Main, the first variable found during our dedicated JCMT/SCUBA-2 monitoring program of eight nearby star-forming regions. EC 53 remained quiescent for the first 6 months of our survey, from February to August 2016. The sub-mm emission began to brighten in September 2016, reached a peak brightness of $1.5$ times the faint state, and has been decaying slowly since February 2017. The change in sub-mm brightness is interpreted as dust heating in the envelope, generated by a luminosity increase of the protostar of a factor of $\ge 4$. The 850~$\mu$m lightcurve resembles the historical $K$-band lightcurve, which varies by a factor of $\sim 6$ with a 543 period and is interpreted as accretion variability excited by interactions between the accretion disk and a close binary system. The predictable detections of accretion variability observed at both near-infrared and sub-mm wavelengths make the system a unique test-bed, enabling us to capture the moment of the accretion burst and to study the consequences of the outburst on the protostellar disk and envelope.

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 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 observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

We present observations of HDCO and H2CO emission toward a sample of 15 Class 0 protostars in the Orion A and B clouds. Of these, eleven protostars are Herschel-identified PACS Bright Red Sources (PBRS) and four are previously identified protostars. Our observations revealed the chemical properties of the PBRS envelope for the first time. The column densities of HDCO and H2CO are derived from single dish observations at an angular resolution of ~20 arcsec (~8400 AU). The degree of deuteration in H2CO ([HDCO]/[H2CO]) was estimated to range from 0.03 to 0.31. The deuterium fractionation of most PBRS (70%) is similar to that of the non-PBRS sources. Three PBRS (30%) exhibit high deuterium fractionation, larger than 0.15. The large variation of the deuterium fractionation of H2CO in the whole PBRS sample may reflect the diversity in the initial conditions of star forming cores. There is no clear correlation between the [HDCO]/[H2CO] ratio and the evolutionary sequence of protostars.

Star formation activities in the NGC 2024 FIR 4 region were studied by imaging centimeter continuum sources and water maser sources using several archival data sets from the Very Large Array. The continuum source VLA 9 is elongated in the northwest-southeast direction, consistent with the FIR 4 bipolar outflow axis, and has a flat spectrum in the 6.2-3.6 cm interval. The three water maser spots associated with FIR 4 are also distributed along the outflow axis. One of the spots is located close to VLA 9, and another one is close to an X-ray source. Examinations of the positions of compact objects in this region suggest that the FIR 4 cloud core contains a single low-mass protostar. VLA 9 is the best indicator of the protostellar position. VLA 9 may be a radio thermal jet driven by this protostar, and it is unlikely that FIR 4 contains a high-mass young stellar object (YSO). A methanol 6.7 GHz maser source is located close to VLA 9, at a distance of about 100 AU. The FIR 4 protostar must be responsible for the methanol maser action, which suggests that methanol class II masers are not necessarily excited by high-mass YSOs. Also discussed are properties of other centimeter continuum sources in the field of view and the water masers associated with FIR 6n. Some of the continuum sources are radio thermal jets, and some are magnetically active young stars.

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

We investigate the spatial density configuration of stars around four metal-poor globular clusters (NGC 6266, NGC 6626, NGC 6642 and NGC 6723) in the Galactic bulge region using wide-field deep J, H, and K imaging data obtained with the WFCAM near-infrared array on United Kingdom Infrared Telescope. Statistical weighted filtering algorithm for the stars on the color-magnitude diagram is applied in order to sort cluster member candidates from the field star contamination. In two-dimensional isodensity contour maps of the clusters, we find that all of the four globular clusters exhibit strong evidence of tidally stripping stellar features beyond tidal radius, in the form of tidal tail or small density lobes or chunk. The orientations of the extended stellar substructures are likely to be associated with the effect of the dynamic interaction with the Galaxy and the cluster space motion. The observed radial density profiles of the four globular clusters also describe the extended substructures; they depart from theoretical King and Wilson models and have an overdensity feature with a break in a slope of profile at the outer region of clusters. The observed results could imply that four globular clusters in the Galactic bulge region have experienced strong environmental effect such as tidal force or bulge or disk shock of the Galaxy in the dynamical evolution of the globular clusters. These observational results provide us further constraints to understand the evolution of clusters in the Galactic bulge region as well as the formation of the Galaxy.

Ten very low luminosity objects were observed multiple times in the 8.5 GHz continuum in search of protostellar magnetic activities. A radio outburst of IRAM 04191+1522 IRS was detected, and the variability timescale was about 20 days or shorter. The results of this survey and archival observations suggest that IRAM 04191+1522 IRS is in active states about half the time. Archival data show that L1014 IRS and L1148 IRS were detectable previously and suggest that at least 20%-30% of very low luminosity protostars are radio variables. Considering the variability timescale and flux level of IRAM 04191+1522 IRS and the previous detection of the circular polarization of L1014 IRS, the radio outbursts of these protostars are probably caused by magnetic flares. However, IRAM 04191+1522 IRS is too young and small to develop an internal convective dynamo. If the detected radio emission is indeed coming from magnetic flares, the discovery implies that the flares may be caused by the fossil magnetic fields of interstellar origin.

We imaged two starless molecular cloud cores, TUKH083 and TUKH122, in the Orion A giant molecular cloud in the CCS and ammonia (NH$_3$) emission with the Very Large Array. TUKH122 contains one NH$_3$ core "TUKH122-n," which is elongated and has a smooth oval boundary. Where observed, the CCS emission surrounds the NH$_3$ core. This configuration resembles that of the N$_2$H$^+$ and CCS distribution in the Taurus starless core L1544, a well-studied example of a dense prestellar core exhibiting infall motions. The linewidth of TUKH122-n is narrow (0.20 km s$^{-1}$) in the NH$_3$ emission line and therefore dominated by thermal motions. The smooth oval shape of the core boundary and narrow linewidth in NH$_3$ seem to imply that TUKH122-n is dynamically relaxed and quiescent. TUKH122-n is similar to L1544 in the kinetic temperature (10 K), linear size (0.03 pc), and virial mass ($\sim$ 2 $M_{\odot}$). Our results strongly suggest that TUKH122-n is on the verge of star formation. TUKH122-n is embedded in the 0.2 pc massive (virial mass $\sim$ 30 $M_{\odot}$) turbulent parent core, while the L1544 NH$_3$ core is embedded in the 0.2 pc less-massive (virial mass $\sim$ 10 $M_{\odot}$) thermal parent core. TUKH083 shows complicated distribution in NH$_3$, but was not detected in CCS. The CCS emission toward TUKH083 appears to be extended, and is resolved out in our interferometric observations.

The results of a maser survey toward ninety-nine protostars in the Orion molecular cloud complex are presented. The target sources are low-mass protostars identified from infrared observations. Single-dish observations were carried out in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. Most of the detected sources were mapped to determine the source positions. Five water maser sources were detected, and they are excited by HH 1-2 VLA 3, HH 1-2 VLA 1, L1641N MM1/3, NGC 2071 IRS 1/3, and an object in the OMC 3 region. The water masers showed significant variability in intensity and velocity with time scales of a month or shorter. Four methanol emission sources were detected, and those in the OMC 2 FIR 3/4 and L1641N MM1/3 regions are probably masers. The methanol emission from the other two sources in the NGC 2071 IRS 1-3 and V380 Ori NE regions are probably thermal. For the water masers, the number of detections per protostar in the survey region is about 2%, which suggests that the water masers of low-mass protostars are rarely detectable. The methanol class I maser of low-mass protostars is an even rarer phenomenon, with a detection rate much smaller than 1%.

The NGC 2024 FIR 6 region was observed in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. The water maser spectra displayed several velocity components and month-scale time variabilities. Most of the velocity components may be associated with FIR 6n, while one component was associated with FIR 4. A typical life time of the water-maser velocity-components is about 8 months. The components showed velocity fluctuations with a typical drift rate of about 0.01 km/s/day. The methanol class I masers were detected toward FIR 6. The methanol emission is confined within a narrow range around the systemic velocity of the FIR 6 cloud core. The methanol masers suggest the existence of shocks driven by either the expanding H II region of FIR 6c or the outflow of FIR 6n.

There may be non-thermal cosmic rays during big-bang nucleosynthesis (BBN) epoch (dubbed as BBNCRs). This paper investigated whether such BBNCRs can be the origin of Lithium problem or not. It can be expected that BBNCRs flux will be small in order to keep the success of standard BBN (SBBN). With favorable assumptions on the BBNCR spectrum between 0.09 -- 4 MeV, our numerical calculation showed that extra contributions from BBNCRs can account for the $^7$Li abundance successfully. However $^6$Li abundance is only lifted an order of magnitude, which is still much lower than the observed value. As the deuteron abundance is very sensitive to the spectrum choice of BBNCRs, the allowed parameter space for the spectrum is strictly constrained. We should emphasize that the acceleration mechanism for BBNCRs in the early universe is still an open question. For example, strong turbulent magnetic field is probably the solution to the problem. Whether such a mechanism can provide the required spectrum deserves further studies.

The NGC 2024 FIR 5/6 region was observed in the 6.9 mm continuum with an angular resolution of about 1.5 arcsec. The 6.9 mm continuum map shows four compact sources, FIR 5w, 5e, 6c, and 6n, as well as an extended structure of the ionization front associated with the optical nebulosity. FIR 6c has a source size of about 0.4 arcsec or 150 AU. The spectral energy distribution (SED) of FIR 6c is peculiar: rising steeply around 6.9 mm and flat around 1 mm. The possibility of a hypercompact H II region is explored. If the millimeter flux of FIR 6c comes from hot ionized gas heated by a single object at the center, the central object may be a B1 star of about 5800 solar luminosities and about 13 solar masses. The 6.9 mm continuum of FIR 6n may be a mixture of free-free emission and dust continuum emission. Archival data show that both FIR 6n and 6c exhibit water maser activity, suggesting the existence of shocked gas around them. The 6.9 mm continuum emission from FIR 5w has a size of about 1.8 arcsec or 760 AU. The SEDs suggest that the 6.9 mm emission of FIR 5w and 5e comes from dust, and the masses of the dense molecular gas are about 0.6 and 0.5 solar masses, respectively.

We obtained 350 and 850 micron continuum maps of the infrared dark cloud G049.40-00.01. Twenty-one dense clumps were identified within G049.40-00.01 based on the 350 micron continuum map with an angular resolution of about 9.6". We present submillimeter continuum maps and report physical properties of the clumps. The masses of clumps range from 50 to 600 M_sun. About 70% of the clumps are associated with bright 24 micron emission sources, and they may contain protostars. The most massive two clumps show extended, enhanced 4.5 micron emission indicating vigorous star-forming activity. The clump size-mass distribution suggests that many of them are forming high mass stars. G049.40-00.01 contains numerous objects in various evolutionary stages of star formation, from pre-protostellar clumps to HII regions.

The NGC 1333 IRAS 4A protobinary was observed in the 1.3 cm and 6.9 mm continuum and the ammonia and SiO lines, with an angular resolution of about 0.4 arcseconds. The continuum maps show the circumstellar structures of the two protostars, A1 and A2. The A1 system is brighter and more massive than the A2 system. The ratio of mass, including dense gas and protostar, is about 6. The properties of the circumstellar disks and outflows suggest that A1 may be younger than A2. The deflected part of the northeastern jet of A2 is bright in the SiO line, and the distance between the brightest peak and deflection point suggests that the enhancement of SiO takes about 100 yr after the collision with a dense core. The ammonia maps show a small structure that seems to be a part of the obstructing core. The outflow properties were studied by comparing interferometric maps of SiO, ammonia, formaldehyde, and HCN lines. Their overall structures agree well, suggesting that these species are excited by the same mechanism. However, the intensity distributions show that SiO is chemically unique. SiO may be directly linked to the primary jet while the other species may be tracing the entrained ambient gas.

The bipolar jet of the NGC 1333 IRAS 4A2 protostar shows a velocity gradient in the direction perpendicular to the jet axis. This lateral velocity gradient can be seen throughout the jet imaged in a silicon monoxide line, 2500-8700 AU from the driving source, and is consistent with the rotation of the accretion disk. If this gradient is caused by the rotation of the jet around its axis, the average specific angular momentum is about 1.5 x 10^21 cm^2 s^-1. Comparison of the kinematics between the jet and the disk suggests that the jet-launching region on the disk has a radius of about 2 AU, which supports the disk-wind models. The angular momentum transported away by the jet seems to be large enough for the protostar to accrete matter from the disk, confirming the crucial role of jets in the early phase of star formation process.

The NGC 1333 IRAS 4A protobinary was observed in the ammonia (2, 2) and (3, 3) lines with an angular resolution of 0.3 arcsec. The ammonia emission source of IRAS 4A2 is elongated in the direction perpendicular to the bipolar jet and has a size of 0.55 arcsec or 130 AU. This emission structure was interpreted as a circumstellar disk around the IRAS 4A2 protostar, and the rotation kinematics of the disk was investigated by making a position-velocity diagram along the major axis. Assuming a Keplerian rotation, the disk has a rotation velocity of 1.8 km/s at a radius of 20 AU, which implies a central object of about 0.08 solar masses. The collapse age of the protostar is about 50,000 yr. The mass, accretion rate, and age are consistent with what are expected from the standard theory of low-mass star formation. If IRAS 4A2 grows at this rate, it may become a star similar to the Sun.

We present the results of an extensive observational study of the active star-forming complex W51 that was observed in the J=2-1 transition of the 12CO and 13CO molecules over a 1.25 deg x 1.00 deg region with the University of Arizona Heinrich Hertz Submillimeter Telescope. We use a statistical equilibrium code to estimate physical properties of the molecular gas. We compare the molecular cloud morphology with the distribution of infrared (IR) and radio continuum sources, and find associations between molecular clouds and young stellar objects (YSOs) listed in Spitzer IR catalogs. The ratios of CO lines associated with HII regions are different from the ratios outside the active star-forming regions. We present evidence of star formation triggered by the expansion of the HII regions and by cloud-cloud collisions. We estimate that about 1% of the cloud mass is currently in YSOs.

It is well known that neutron stars spin down due to magnetic dipole radiation. The deconfinement phase transition of hadron matter to quark matter is expected to occur in the dense cores of the stars during spin-down. The phase transition continuously takes place inducing not only structural changes but also energy release in case of a first-order phase transition. The generation of energy increases the internal energy of the stars which is called deconfinement heating. The temperature of the stars arise when the deconfinement heating appears in the cores. We explore the deconfinement signature by studying the changes of surface temperature in the thermal evolution process of neutron stars.

We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). We use the shadow features to derive optical depth maps. We derive molecular masses for the cores and the surrounding environment; we find that the 24 micron shadow masses are always greater than or equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze--out onto dust grains. We combine this sample with two additional cores that we studied previously to bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~ 2/3 of the cores selected to have prominent 24 micron shadows are collapsing or near collapse, a result that is supported by millimeter line observations. Of this subset at least half have indications of 70 micron shadows. All cores observed to produce absorption features at 70 micron are close to collapse. We conclude that 24 micron shadows, and even more so the 70 micron ones, are useful markers of cloud cores that are approaching collapse.

We present 737 candidate Young Stellar Objects (YSOs) near the W51 Giant Molecular Cloud (GMC) over an area of 1.25 deg x 1.00 deg selected from Spitzer Space Telescope data. We use spectral energy distribution (SED) fits to identify YSOs and distinguish them from main-sequence or red giant stars, asymptotic giant branch stars, and background galaxies. Based on extinction of each YSO, we separate a total of 437 YSOs associated with the W51 region from the possible foreground sources. We identify 69 highly embedded Stage 0/I candidate YSOs in our field with masses > 5 Msun (corresponding to mid-to early-B main-sequence spectral types) 46 of which are located in the central active star forming regions of W51A and W51B. From the YSOs associated with W51, we find evidence for mass segregation showing that the most massive YSOs are concentrated on the W51 HII region complex. We find a variation in the spatial distribution of the mass function (MF) of YSOs in the mass range between 5 Msun and 18 Msun. The derived slopes of the MF are -1.26 and -2.36 in the active star-forming region and the outer region, respectively. The variation of the MF for YSOs embedded in the molecular cloud implies that the distribution of stellar masses in clusters depends on the local conditions in the parent molecular cloud.

We present Heinrich Hertz Telescope CO observations of the shell structure near the active star-forming complex W51A to investigate the process of star formation triggered by the expansion of an HII region. The CO observations confirm that dense molecular material has been collected along the shell detected in Spitzer IRAC images. The CO distribution shows that the shell is blown out toward a lower density region to the northwest. Total hydrogen column density around the shell is high enough to form new stars. We find two CO condensations with the same central velocity of 59 km/s to the east and north along the edge of the IRAC shell. We identify two YSOs in early evolutionary stages (Stage 0/I) within the densest molecular condensation. From the CO kinematics, we find that the HII region is currently expanding with a velocity of 3.4 km/s, implying that the shell's expansion age is ~1 Myr. This timescale is in good agreement with numerical simulations of the expansion of the HII region (Hosokawa et al. 2006). We conclude that the star formation on the border of the shell is triggered by the expansion of the HII region.

M17 is one of the youngest and most massive nearby star-formation regions in the Galaxy. It features a bright H II region erupting as a blister from the side of a giant molecular cloud (GMC). Combining photometry from the Spitzer GLIMPSE survey with complementary infrared (IR) surveys, we identify candidate young stellar objects (YSOs) throughout a 1.5 deg x 1 deg field that includes the M17 complex. The long sightline through the Galaxy behind M17 creates significant contamination in our YSO sample from unassociated sources with similar IR colors. Removing contaminants, we produce a highly-reliable catalog of 96 candidate YSOs with a high probability of association with the M17 complex. We fit model spectral energy distributions to these sources and constrain their physical properties. Extrapolating the mass function of 62 intermediate-mass YSOs (M >3 Msun), we estimate that >1000 stars are in the process of forming in the extended outer regions of M17. From IR survey images from IRAS and GLIMPSE, we find that M17 lies on the rim of a large shell structure ~0.5 deg in diameter (~20 pc at 2.1 kpc). We present new maps of CO and 13CO (J=2-1) emission, which show that the shell is a coherent, kinematic structure associated with M17 at v = 19 km/s. The shell is an extended bubble outlining the photodissociation region of a faint, diffuse H II region several Myr old. We provide evidence that massive star formation has been triggered by the expansion of the bubble. The formation of the massive cluster ionizing the M17 H II region itself may have been similarly triggered. We conclude that the star formation history in the extended environment of M17 has been punctuated by successive waves of massive star formation propagating through a GMC complex.

We present infrared and millimeter observations of Barnard 335, the prototypical isolated Bok globule with an embedded protostar. Using Spitzer data we measure the source luminosity accurately; we also constrain the density profile of the innermost globule material near the protostar using the observation of an 8.0 um shadow. HHT observations of 12CO 2 --> 1 confirm the detection of a flattened molecular core with diameter ~10000 AU and the same orientation as the circumstellar disk (~100 to 200 AU in diameter). This structure is probably the same as that generating the 8.0 um shadow and is expected from theoretical simulations of collapsing embedded protostars. We estimate the mass of the protostar to be only ~5% of the mass of the parent globule.

The HH 24 MMS protostellar system was observed in the 6.9 mm continuum with a high angular resolution (0.5"). HH 24 MMS was resolved into two sources. The separation between sources 1 and 2 is ~0.9" or 360 AU. The spectral energy distribution suggests that the 6.9 mm flux is almost entirely from dust. The 6.9 mm image and the spectrum suggest that HH 24 MMS may be a protostellar binary system. Total mass including the accretion disks and the inner protostellar envelope is ~1.4 Msun. Disk masses of sources 1 and 2 are 0.6 and 0.3 Msun, respectively. Both sources are highly elongated. The difference in the position angle of the two disks is ~45 degree, which means that HH 24 MMS is a highly misaligned protobinary system. The misalignment suggests that turbulent fragmentation may be the formation mechanism relevant to the binary systems with a separation of a few hundreds of AU, such as the HH 24 MMS system.

The interior of neutron stars contains nuclear matter at very high density for numerous subatomic particles compete with each other. Therefore, confirming the components and properties there is our significant task. Here we summarize the possible methods especial the way of r-mode instability to probe into the neutron star and show our some results. The KHz pulsar in XTE J1739-285 may give a significant implication

We derive the bulk viscous damping timescale of hybrid stars, neutron stars with quark matter core. The r-mode instability windows of the stars show that the theoretical results are consistent with the rapid rotation pulsar data, which may give an indication for the existence of quark matter in the interior of neutron stars. Hybrid stars instead of neutron or strange stars may lead to submillisecond pulsars.

The deconfinement phase transition will lead to the release of latent heat during spins down of neutron stars if the transition is the first-order one.We have investigated the thermal evolution of neutron stars undergoing such deconfinement phase transition. The results show that neutron stars may be heated to higher temperature.This feature could be particularly interesting for high temperature of low-magnetic field millisecond pulsar at late stage.

The NGC 1333 IRAS 4A protobinary was observed in the ammonia (2, 2) and (3, 3) lines and in the 1.3 cm continuum with a high resolution (about 1.0 arcsec). The ammonia maps show two compact sources, one for each protostar, and they are probably protostellar accretion disks. The disk associated with IRAS 4A2 is seen nearly edge-on and shows an indication of rotation. The A2 disk is brighter in the ammonia lines but dimmer in the dust continuum than its sibling disk, with the ammonia-to-dust flux ratios different by about an order of magnitude. This difference suggests that the twin disks have surprisingly dissimilar characters, one gas-rich and the other dusty. The A2 disk may be unusually active or hot, as indicated by its association with water vapor masers. The existence of two very dissimilar disks in a binary system suggests that the formation process of multiple systems has a controlling agent lacking in the isolated star formation process and that stars belonging to a multiple system do not necessarily evolve in phase with each other.