We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $\Sigma$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$\mu$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A.

We identify 45 protostellar outflows in CO maps of the Orion A giant molecular cloud from the CARMA-NRO Orion survey. Our sample includes 11 newly detected outflows. We measure the mass and energetics of the outflows, including material at low-velocities by correcting for cloud contributions. The total momentum and kinetic energy injection rates of outflows is comparable to the turbulent dissipation rate of the cloud. We also compare the outflow position angles to the orientation of C$^{18}$O filaments. We find that the full sample of outflows is consistent with being randomly oriented with respect to the filaments. A subsample of the most reliable measurements shows a moderately perpendicular outflow-filament alignment which may reflect accretion of mass across filaments and onto the protostellar cores.

Since the start of ALMA observatory operation, new and important chemistry of infrared cold core was revealed. Molecular transitions at millimeter range are being used to identify and to characterize these sources. We have investigated the 231 GHz ALMA archive observations of the infrared dark cloud region C9, focusing on the brighter source that we called as IRDC-C9 Main. We report the existence of two sub-structures on the continuum map of this source: a compact bright spot with high chemistry diversity that we labelled as core, and a weaker and extended one, that we labelled as tail. In the core, we have identified lines of the molecules OCS(19-18), $^{13}$CS(5-4) and CH$_{3}$CH$_{2}$CN, several lines of CH$_{3}$CHO and the k-ladder emission of $^{13}$CH$_{3}$CN.We report two different temperature regions: while the rotation diagram of CH$_{3}$CHO indicates a temperature of 25 K, the rotation diagram of $^{13}$CH$_{3}$CN indicates a warmer phase at temperature of $\sim450$K. In the tail, only the OCS(19-18) and $^{13}$CS(5-4) lines were detected. We used the $Nautilus$ and the \textsc{Radex} codes to estimate the column densities and the abundances. The existence of hot gas in the core of IRDC-C9 Main suggests the presence of a protostar, which is not present in the tail.

We conducted an exploration of 12CO molecular outflows in the Orion A giant molecular cloud to investigate outflow feedback using 12CO (J = 1-0) and 13CO (J = 1-0) data obtained by the Nobeyama 45-m telescope. In the region excluding the center of OMC 1, we identified 44 12CO (including 17 newly detected) outflows based on the unbiased and systematic procedure of automatically determining the velocity range of the outflows and separating the cloud and outflow components. The optical depth of the 12CO emission in the detected outflows is estimated to be approximately 5. The total momentum and energy of the outflows, corrected for optical depth, are estimated to be 1.6 x 10 2 M km s-1 and 1.5 x 10 46 erg, respectively. The momentum and energy ejection rate of the outflows are estimated to be 36% and 235% of the momentum and energy dissipation rates of the cloud turbulence, respectively. Furthermore, the ejection rates of the outflows are comparable to those of the expanding molecular shells estimated by Feddersen et al. (2018, ApJ, 862, 121). Cloud turbulence cannot be sustained by the outflows and shells unless the energy conversion efficiency is as high as 20%.

We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degree, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in 3D. Such alignment may result from filament fragmentation or continuous mass transportation from filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields.

We investigate the relationship between turbulence and feedback in the Orion A molecular cloud using maps of $^{12}$CO(1-0), $^{13}$CO(1-0) and C$^{18}$O(1-0) from the CARMA-NRO Orion survey. We compare gas statistics with the impact of feedback in different parts of the cloud to test whether feedback changes the structure and kinematics of molecular gas. We use principal component analysis, the spectral correlation function, and the spatial power spectrum to characterize the cloud. We quantify the impact of feedback with momentum injection rates of protostellar outflows and wind-blown shells as well as the surface density of young stars. We find no correlation between shells or outflows and any of the gas statistics. However, we find a significant anti-correlation between young star surface density and the slope of the $^{12}$CO spectral correlation function, suggesting that feedback may influence this statistic. While calculating the principal components, we find peaks in the covariance matrix of our molecular line maps offset by 1-3 km s$^{-1}$ toward several regions of the cloud which may be produced by feedback. We compare these results to predictions from molecular cloud simulations.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

We present high resolution (0.2", 1000 AU) 1.3 mm ALMA observations of massive infrared dark cloud clump, G028.37+00.07-C1, thought to harbor the early stages of massive star formation. Using $\rm N_2D^+$(3-2) we resolve the previously identified C1-S core, separating the bulk of its emission from two nearby protostellar sources. C1-S is thus identified as a massive ($\sim50\:M_\odot$), compact ($\sim0.1\:$pc diameter) starless core, e.g., with no signs of outflow activity. Being highly deuterated, this is a promising candidate for a pre-stellar core on the verge of collapse. An analysis of its dynamical state indicates a sub-virial velocity dispersion compared to a trans-Alfv\'enic turbulent core model. However, virial equilibrium could be achieved with sub-Alfv\'enic conditions involving $\sim2\:$mG magnetic field strengths.

We study the core mass function (CMF) within 32 dense clumps in seven infrared dark clouds (IRDCs) with the Atacama Large Millimeter/submillimeter Array (ALMA) via 1.3~mm continuum emission at a resolution of $\sim$1". We have identified 107 cores with the dendrogram algorithm, with a median radius of about 0.02 pc. Their masses range from 0.261 to 178 $M_{\odot}$. After applying completeness corrections, we fit the combined IRDC CMF with a power law of the form $d N / d\:{\rm log} M \propto M^{-\alpha}$ and derive an index of $\alpha\simeq0.86\pm0.11$ for $M \geq 0.79\:M_\odot$ and $\alpha\simeq0.70\pm0.13$ for $M\geq 1.26\:M_\odot$, which is a significantly more top-heavy distribution than the Salpeter stellar initial mass function (IMF) index of 1.35. We also make a direct comparison of these IRDC clump CMF results to those measured in the more evolved protocluster G286 derived with similar methods, which have $\alpha\simeq1.29\pm0.19$ and $1.08\pm0.27$ in these mass ranges, respectively. These results provide a hint that, especially for the $M\geq 1.26\:M_\odot$ range where completeness corrections are modest, the CMF in high pressure, early-stage environments of IRDC clumps may be top-heavy compared to that in the more evolved, global environment of the G286 protoclusters. However, larger samples of cores probing these different environments are needed to better establish the robustness of this potential CMF variation.

POLAR is a compact wide-field space-borne detector for precise measurements of the linear polarisation of hard X-rays emitted by transient sources in the energy range from 50 keV to 500 keV. It consists of a 40$\times$40 array of plastic scintillator bars used as a detection material. The bars are grouped in 25 detector modules. The energy range sensitivity of POLAR is optimized to match with the prompt emission photons from the gamma-ray bursts (GRBs). Polarization measurements of the prompt emission would probe source geometries, emission mechanisms and magnetic structures in GRB jets. The instrument can also detect hard X-rays from solar flares and be used for precise measurement of their polarisation. POLAR was launched into a low Earth orbit on-board the Chinese space-lab TG-2 on September 15th, 2016. To achieve high accuracies in polarisation measurements it is essential to assure both before and during the flight a precise energy calibration. Such calibrations are performed with four low activity $^{22}$Na radioactive sources placed inside the instrument. Energy conversion factors are related to Compton edge positions from the collinear annihilation photons from the sources. This paper presents main principles of the in-flight calibration, describes studies of the method based on Monte Carlo simulations and its laboratory verification and finally provides some observation results based on the in-flight data analysis.

We present the discovery of expanding spherical shells around low to intermediate-mass young stars in the Orion A giant molecular cloud using observations of $^{12}$CO (1-0) and $^{13}$CO (1-0) from the Nobeyama Radio Observatory 45-meter telescope. The shells have radii from 0.05 to 0.85 pc and expand outward at 0.8 to 5 km/s. The total energy in the expanding shells is comparable to protostellar outflows in the region. Together, shells and outflows inject enough energy and momentum to maintain the cloud turbulence. The mass-loss rates required to power the observed shells are two to three orders of magnitude higher than predicted for line-driven stellar winds from intermediate-mass stars. This discrepancy may be resolved by invoking accretion-driven wind variability. We describe in detail several shells in this paper and present the full sample in the online journal.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA to carry out a survey of dense gas and cores in the central region of the massive ($\sim10^5\:M_\odot$) Infrared Dark Cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12m-array and produces an unprecedented view of the densest structures of this IRDC. In this first paper about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the "dense gas" detection probability function (DPF), i.e., as a function of the local mass surface density, $\Sigma$, for various choices of thresholds of mm continuum emission to define "dense gas". We then estimate the dense gas mass fraction, $f_{\rm dg}$, in the central region of the IRDC and, via extrapolation with the DPF and the known $\Sigma$ probability distribution function, to the larger-scale surrounding regions, finding values of about 5\% to 15\% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of $\epsilon_{\rm ff}\sim$10\%, with approximately a factor of two systematic uncertainties.

We study the core mass function (CMF) of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array via 1.3~mm continuum emission at a resolution of 1.0\arcsec\ (2500~au). We have mapped a field of 5.3\arcmin$\times$5.3\arcmin\ centered on the protocluster clump. We measure the CMF in the central region, exploring various core detection algorithms, which give source numbers ranging from 60 to 125, depending on parameter selection. We estimate completeness corrections due to imperfect flux recovery and core identification via artificial core insertion experiments. For masses $M\gtrsim1\:M_\odot$, the fiducial dendrogram-identified CMF can be fit with a power law of the form ${\rm{d}}N/{\rm{d}}{\rm{log}}M\propto{M}^{-\alpha}$ with $\alpha \simeq1.24\pm0.17$, slightly shallower than, but still consistent with, the index of the Salpeter stellar initial mass function of 1.35. Clumpfind-identified CMFs are significantly shallower with $\alpha\simeq0.64\pm0.13$. While raw CMFs show a peak near $1\:M_\odot$, completeness-corrected CMFs are consistent with a single power law extending down to $\sim 0.5\:M_\odot$, with only a tentative indication of a shallowing of the slope around $\sim1\:M_\odot$. We discuss the implications of these results for star and star cluster formation theories.

POLAR is a new satellite-born detector aiming to measure the polarization of an unprecedented number of Gamma-Ray Bursts in the 50-500 keV energy range. The instrument, launched on-board the Tiangong-2 Chinese Space lab on the 15th of September 2016, is designed to measure the polarization of the hard X-ray flux by measuring the distribution of the azimuthal scattering angles of the incoming photons. A detailed understanding of the polarimeter and specifically of the systematic effects induced by the instrument's non-uniformity are required for this purpose. In order to study the instrument's response to polarization, POLAR underwent a beam test at the European Synchrotron Radiation Facility in France. In this paper both the beam test and the instrument performance will be described. This is followed by an overview of the Monte Carlo simulation tools developed for the instrument. Finally a comparison of the measured and simulated instrument performance will be provided and the instrument response to polarization will be presented.

We have developed ultra-low-background NaI(Tl) crystals to reproduce the DAMA results with the ultimate goal of achieving purity levels that are comparable to or better than those of the DAMA/LIBRA crystals. Even though the achieved background level does not approach that of DAMA/LIBRA, it is crucial to have a quantitative understanding of the backgrounds. We have studied background simulations toward a deeper understanding of the backgrounds and developed background models for a 9.16-kg NaI(Tl) crystal used in the test arrangement. In this paper we describe the contributions of background sources quantitatively by performing Geant4 Monte Carlo simulations that are fitted to the measured data to quantify the unknown fractions of the background compositions. In the fitted results, the overall simulated background spectrum well describes the measured data with a 9.16-kg NaI(Tl) crystal and shows that the background sources are dominated by surface $^{210}$Pb and internal $^{40}$K in the 2 to 6-keV energy interval, which produce 2.4 counts/day/keV/kg (dru) and 0.5 dru, respectively.

We present a multiwavelength investigation of a region of a nearby giant molecular cloud that is distinguished by a minimal level of star formation activity. With our new 12CO(J=2-1) and 13CO(J=2-1) observations of a remote region within the middle of the California molecular cloud, we aim to investigate the relationship between filaments, cores, and a molecular outflow in a relatively pristine environment. An extinction map of the region from Herschel Space Observatory observations reveals the presence of two 2-pc-long filaments radiating from a high-extinction clump. Using the 13CO observations, we show that the filaments have coherent velocity gradients and that their mass-per-unit-lengths may exceed the critical value above which filaments are gravitationally unstable. The region exhibits structure with eight cores, at least one of which is a starless, prestellar core. We identify a low-velocity, low-mass molecular outflow that may be driven by a flat spectrum protostar. The outflow does not appear to be responsible for driving the turbulence in the core with which it is associated, nor does it provide significant support against gravitational collapse.

Gamma-ray polarimetry is a new powerful tool to study the processes responsible for the emission from astrophysical sources and the environments in which this emission takes place. Few successful polarimetric measurements have however been performed thus far in the gamma-ray energy band due to the difficulties involved. POLAR is a dedicated polarimeter designed to perform high precision measurements of the polarization of the emission from gamma-ray burst in the 50-500 keV energy range. This new polarimeter is expected to detect approximately 50 gamma-ray bursts per year while performing high precision polarization measurements on approximately 10 bursts per year. The instrument was launched into lower earth orbit as part of the second Chinese space lab, the Tiangong-2, on September 15th 2016 and has been taking data successfully since being switched on one week after. The instrument uses a segmented scintillator array consisting of 1600 plastic scintillator bars, read out by 25 flat-panel multi-anode photomultipliers, to measure the Compton scattering angles of incoming photons. The small segmentation and relatively large uniform effective area allow the instrument to measure the polarization of a large number of transient events, such as gamma-ray bursts, with an unprecedented precision during its two year life-time. The final flight model underwent detailed calibration prior to launch as well as intensive space qualification tests, a summary of which will be presented in this paper. The instrument design will be discussed first followed by an overview of the on-ground tests, finally the in-orbit behavior as measured during the first weeks of the mission will be presented.

We carry out an ALMA $\rm N_2D^+$(3-2) and 1.3~mm continuum survey towards 32 high mass surface density regions in seven Infrared Dark Clouds with the aim of finding massive starless cores, which may be the initial conditions for the formation of massive stars. Cores showing strong $\rm N_2D^+$(3-2) emission are expected to be highly deuterated and indicative of early, potentially pre-stellar stages of star formation. We also present maps of these regions in ancillary line tracers, including C$^{18}$O(2-1), DCN(3-2) and DCO$^+$(3-2). Over 100 $\rm N_2D^+$ cores are identified with our newly developed core-finding algorithm based on connected structures in position-velocity space. The most massive core has $\sim70\:M_\odot$ (potentially $\sim170\:M_\odot$) and so may be representative of the initial conditions for massive star formation. The existence and dynamical properties of such cores constrain massive star formation theories. We measure the line widths and thus velocity dispersion of six of the cores with strongest $\rm N_2D^+$(3-2) line emission, finding results that are generally consistent with virial equilibrium of pressure confined cores.

High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitation conditions and deuterium fraction, $D_{\rm frac}^{\rm N_2H^+} \equiv [\rm N_2D^+]/[N_2H^+]$, with values of $D_{\rm frac}^{\rm N_2H^+} \simeq 0.2$--$0.7$, several orders of magnitude above the cosmic [D]/[H] ratio. Additional observations of o-H$_2$D$^+$ are also presented that help constrain the ortho-to-para ratio of $\rm H_2$, which is a key quantity affecting the degree of deuteration. We then present chemodynamical modeling of the two cores, exploring especially the implications for the collapse rate relative to free-fall, $\alpha_{\rm ff}$. In order to reach the high level of observed deuteration of $\rm N_2H^+$, we find that the most likely evolutionary history of the cores involves collapse at a relatively slow rate, $\lesssim1/10$th of free-fall.

We present ALMA follow-up observations of two massive, early-stage core candidates, C1-N & C1-S, in Infrared Dark Cloud (IRDC) G028.37+00.07, which were previously identified by their N2D+(3-2) emission and show high levels of deuteration of this species. The cores are also dark at far infrared wavelengths up to ~100 microns. We detect 12CO(2-1) from a narrow, highly-collimated bipolar outflow that is being launched from near the center of the C1-S core, which is also the location of the peak 1.3mm dust continuum emission. This protostar, C1-Sa, has associated dense gas traced by C18O(2-1) and DCN(3-2), from which we estimate it has a radial velocity that is near the center of the range exhibited by the C1-S massive core. A second outflow-driving source is also detected within the projected boundary of C1-S, but appears to be at a different radial velocity. After considering properties of the outflows, we conclude C1-Sa is a promising candidate for an early-stage massive protostar and as such it shows that these early phases of massive star formation can involve highly ordered outflow, and thus accretion, processes, similar to models developed to explain low-mass protostars.

Infrared Dark Clouds (IRDCs) are cold, dense regions that are usually found within Giant Molecular Clouds (GMCs). Ongoing star formation within IRDCs is typically still deeply embedded within the surrounding molecular gas. Characterising the properties of relatively quiescent IRDCs may therefore help us to understand the earliest phases of the star formation process. Studies of local molecular clouds have revealed that deuterated species are enhanced in the earliest phases of star formation. In this paper we test this towards IRDC G035.39-00.33. We present an 80 arcsec by 140 arcsec map of the J=2-1 transition of N2D+, obtained with the IRAM-30m telescope. We find that N2D+ is widespread throughout G035.39-00.33. Complementary observations of N2H+(1-0) are used to estimate the deuterium fraction, N(N2D+)/N(N2H+). We report a mean deuterium fraction of 0.04+-0.01, with a maximum of 0.09+-0.02. The mean deuterium fraction is ~3 orders of magnitude greater than the interstellar [D]/[H] ratio. High angular resolution observations are required to exclude beam dilution effects of compact deuterated cores. Using chemical modelling, we find that the average observed values of the deuterium fraction are in agreement with an equilibrium deuterium fraction, given the general properties of the cloud. This implies that the IRDC is at least ~3Myr old, which is ~8 times longer than the mean free-fall time of the observed deuterated region.

Gamma Ray Bursts (GRBs) are the strongest explosions in the universe which might be associated with creation of black holes. Magnetic field structure and burst dynamics may influence polarization of the emitted gamma-rays. Precise polarization detection can be an ultimate tool to unveil the true GRB mechanism. POLAR is a space-borne Compton scattering detector for precise measurements of the GRB polarization. It consists of a 40$\times$40 array of plastic scintillator bars read out by 25 multi-anode PMTs (MaPMTs). It is scheduled to be launched into space in 2016 onboard of the Chinese space laboratory TG2. We present a dedicated methodology for POLAR calibration and some calibration results based on the combined use of the laboratory radioactive sources and polarized X-ray beams from the European Synchrotron Radiation Facility. They include calibration of the energy response, computation of the energy conversion factor vs. high voltage as well as determination of the threshold values, crosstalk contributions and polarization modulation factors.

To understand massive star formation requires study of its initial conditions. Two massive starless core candidates, C1-N & C1-S, have been detected in IRDC G028.37+00.07 in $\rm N_2D^+$(3-2) with $ALMA$. From their line widths, either the cores are subvirial and are thus young structures on the verge of near free-fall collapse, or they are threaded by $\sim1$ mG $B$-fields that help support them in near virial equilibrium and potentially have older ages. We modeled the deuteration rate of $\rm N_2H^+$ to constrain collapse rates of the cores. First, to measure their current deuterium fraction, $D_{\rm frac}^{\rm N_2H^+}$ $\equiv [\rm N_2D^+]/[N_2H^+]$, we observed multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ with $CARMA$, $SMA$, $JCMT$, $NRO~45m$ and $IRAM~30m$, to complement the $ALMA$ data. For both cores we derived $D_{\rm frac}^{\rm N_2H^+}\sim0.3$, several orders of magnitude above the cosmic [D]/[H] ratio. We then carried out chemodynamical modeling, exploring how collapse rate relative to free-fall, $\alpha_{\rm ff}$, affects the level of $D_{\rm frac}^{\rm N_2H^+}$ that is achieved from a given initial condition. To reach the observed $D_{\rm frac}^{\rm N_2H^+}$, most models require slow collapse with $\alpha_{\rm ff}\sim0.1$, i.e., $\sim1/10$th of free-fall. This makes it more likely that the cores have been able to reach a near virial equilibrium state and we predict that strong $B$-fields will eventually be detected. The methods developed here will be useful for measurement of the pre-stellar core mass function.

Considering the contribution of the emission from the host galaxies of gamma-ray bursts (GRBs) to the radio afterglows, we investigate the effect of host galaxies on observations statistically. For the three types of events, e.g. low-luminosity, standard and high-luminosity GRBs, it is found that a tight correlation exists between the ratio of the radio flux (RRF) of host galaxy to the total radio peak emission and the observational frequency. Especially, toward lower frequencies, the contribution from the host increases significantly. The correlation can be used to get a useful estimate for the radio brightness of those host galaxies which only have very limited radio afterglow data. Using this prediction, we re-considered the theoretical radio afterglow light curves for four kinds of events, i.e. high-luminosity, low-luminosity, standard and failed GRBs, taking into account the contribution from the host galaxies and aiming at exploring the detectability of these events by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Lying at a typical redshift of $z=1$, most of the events can be detected by FAST easily. For the less fierce low-luminosity GRBs, their radio afterglows are not strong enough to exceed the sensitivity limit of FAST at such distances. However, since a large number of low luminosity bursts actually happen very near to us, it is expected that FAST will still be able to detect many of them.

A deep, wide-field, near-infrared imaging survey was used to construct an extinction map of the southeastern part of the California Molecular Cloud (CMC) with $\sim$ 0.5 arc min resolution. The same region was also surveyed in the $^{12}$CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1) emission lines at the same angular resolution. Strong spatial variations in the abundances of $^{13}$CO and C$^{18}$O were found to be correlated with variations in gas temperature, consistent with temperature dependent CO depletion/desorption on dust grains. The $^{13}$CO to C$^{18}$O abundance ratio was found to increase with decreasing extinction, suggesting selective photodissociation of C$^{18}$O by the ambient UV radiation field. The cloud averaged X-factor is found to be $<$X$_{\rm CO}$$>$ $=$ 2.53 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$, somewhat higher than the Milky Way average. On sub-parsec scales we find no single empirical value of the X-factor that can characterize the molecular gas in cold (T$_{\rm k}$ $\lesssim$ 15 K) regions, with X$_{\rm CO}$ $\propto$ A$_{\rm V}$$^{0.74}$ for A$_{\rm V}$ $\gtrsim$ 3 magnitudes. However in regions containing relatively hot (T$_{\rm ex}$ $\gtrsim$ 25 K) gas we find a clear correlation between W($^{12}$CO) and A$_{\rm V}$ over a large (3 $\lesssim$ A$_{\rm V}$ $\lesssim$ 25 mag) extinction range. This suggests a constant X$_{\rm CO}$ $=$ 1.5 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$ for the hot gas, a lower value than either the average for the CMC or Milky Way. We find a correlation between X$_{\rm CO}$ and T$_{\rm ex}$ with X$_{\rm CO}$ $\propto$ T$_{\rm ex}$$^{-0.7}$ suggesting that the global X-factor of a cloud may depend on the relative amounts of hot gas within it.

The deuterium fraction [N$_2$D$^+$]/[N$_2$H$^+$], may provide information about the ages of dense, cold gas structures, important to compare with dynamical models of cloud core formation and evolution. Here we introduce a complete chemical network with species containing up to three atoms, with the exception of the Oxygen chemistry, where reactions involving H$_3$O$^+$ and its deuterated forms have been added, significantly improving the consistency with comprehensive chemical networks. Deuterium chemistry and spin states of H$_2$ and H$_3^+$ isotopologues are included in this primarily gas-phase chemical model. We investigate dependence of deuterium chemistry on model parameters: density ($n_{\rm H}$), temperature, cosmic ray ionization rate, and gas-phase depletion factor of heavy elements ($f_{\rm D}$). We also explore the effects of time-dependent freeze-out of gas-phase species and dynamical evolution of density at various rates relative to free-fall collapse. For a broad range of model parameters, the timescales to reach large values of $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$, observed in some low- and high-mass starless cores, are relatively long compared to the local free-fall timescale. These conclusions are unaffected by introducing time-dependent freeze-out and considering models with evolving density, unless the initial $f_{\rm D} \gtrsim$ 10. For fiducial model parameters, achieving $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$ requires collapse to be proceeding at rates at least several times slower than that of free-fall collapse, perhaps indicating a dynamically important role for magnetic fields in the support of starless cores and thus the regulation of star formation.

Using the generic hydrodynamic model of Gamma-Ray Burst (GRB) afterglows, we calculate the radio afterglow light curves of low luminosity, high luminosity, failed and standard GRBs in different observational bands of FAST's energy window. The GRBs are assumed to be located at different distances from us. Our results rank the detectability of GRBs in descending order as high luminosity, standard, failed and low luminosity GRBs. We predict that almost all types of radio afterglows except that of low luminosity GRBs could be observed by the large radio telescope as long as the domains of time and frequency are appropriate. It is important to note that FAST can detect relatively weak radio afterglows at a higher frequency of 2.5 GHz for very high redshift up to z=15 or even more. Radio afterglows of low luminosity GRBs can be detected only after the completion of the 2nd phase of FAST. FAST is expected to significantly expand the sample of GRB radio afterglows in the near future.

How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their surrounding clump medium. Their internal pressure is provided by a combination of turbulence and magnetic fields. Alternatively, the Competitive Accretion model requires strongly sub-virial initial conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate- and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we have identified four prime examples of massive (~100Msun) clumps from mid-infrared extinction mapping of infrared dark clouds (IRDCs). Fontani et al. found high deuteration fractions of N2H+ in these objects, which are consistent with them being starless. Here we present ALMA observations of these four clumps that probe the N2D+(3-2) line at 2.3" resolution. We find six N2D+ cores and determine their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the predictions of the Turbulent Core model of self-gravitating, magnetized (with Alfven Mach number m_A~1) and virialized cores that are bounded by the high pressures of their surrounding clumps. However, in the most massive cores, with masses up to ~60Msun, our results suggest that moderately enhanced magnetic fields (so that m_A~0.3) may be needed for the structures to be in virial and pressure equilibrium. Magnetically regulated core formation may thus be important in controlling the formation of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar IMF.

PSR B1259-63/LS 2883 is a binary system in which a 48-ms pulsar orbits around a Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. It is special for having asymmetric two-peak profiles in both the X-ray and the TeV light curves. Recently, an unexpected GeV flare was detected by $Fermi$ gamma-ray observatory several weeks after the last periastron passage. In this paper, we show that this observed GeV flare could be produced by the Doppler-boosted synchrotron emission in the bow shock tail. An anisotropic pulsar wind model, which mainly affects the energy flux injection to the termination shock in different orbital phase, is also used in this paper, and we find that the anisotropy in the pulsar wind can play a significant role in producing the asymmetric two-peak profiles in both X-ray and TeV light curves. The X-ray and TeV photons before periastron are mainly produced by the shocked electrons around the shock apex and the light curves after periastron are contributed by the emission from the shock apex and the shock tail together, which result in the asymmetric two-peak light curves.

The model of non-ballistic jet motion proposed in 2008 provides a simple explanation to the inward jet motion and bent jet. Recently, evidences of such a non-radial motion increase rapidly, and more complicated morphologies appear. On the other hand, the ballistic plus precession model likely holds in majority samples of jet motion. This paper discusses the relationship between the ballistic and non-ballistic model of jet motion, which suggests that the interaction of ejectors with ambient matter can produce knots at different stages of evolution and hence different separations to the core. And as a jet precesses, knots produced between the core and the deceleration radius result in spiral pattern expected by the model of ballistic plus precession; and knots generated at the deceleration radius display non-radial motion such as bent jet or oscillation of ridge-line. This paper develops the first non-ballistic model in four aspects. Firstly, it provides a numerical simulation to the production of multi-knot for a precessing jet. Secondly, it fits the precession behavior of multi-knot and interprets the oscillation of ridge lines like S5 1803+784. Thirdly, it gives an unified interpretation to the bent jet applicable to both multi-knot and single knot. And fourthly, the problem of very large numbers of observed outward motions as opposed to the inward ones is addressed in a new scope.

PSR B1259-63/SS 2883 is a binary system in which a 48-ms pulsar orbits around a Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. Extensive broadband observational data are available for this system from radio band to very high energy (VHE) range. The multi-frequency emission is unpulsed and nonthermal, and is generally thought to be related to the relativistic electrons accelerated from the interaction between the pulsar wind and the stellar wind, where X-ray emission is from the synchrotron process and the VHE emission is from the inverse Compton (IC) scattering process. Here a shocked wind model with variation of the magnetic parameter $\sigma$ is developed for explaining the observations. By choosing proper param- eters, our model could reproduce two-peak profile in X-ray and TeV light curves. The effect of the disk exhibits an emission and an absorption components in the X-ray and TeV bands respectively. We suggest that some GeV flares will be produced by Doppler boosting the synchrotron spectrum. This model can possibly be used and be checked in other similar systems such as LS I+61o303 and LS 5039.

A mapping study of IRAS 05553+1631 was performed with 12CO J=3-2 and 13CO J=2-1 lines observed by the KOSMA 3 m-telescope. A core with a size of 0.65 pc and with a LTE mass of 120 M\odot was defined by the mapping with 13CO J=2-1 line. We have identified a bipolar outflow with 12CO J=3-2. For accuracy in the calculation of outflow parameters, overcoming the projection effect is important. We propose a new method to directly calculate the inclination-angle {\theta}. We establish two basic equations with the help of outflow contour diagram and finally obtain the "angle function" and the "angle equation" to derive {\theta}. We apply our method to the outflow of IRAS 05553+1631, finding that {\theta}blue is 73\circ and {\theta}red is 78\circ. Compared to the parameters initially estimated under an assumption of 45\circ inclination-angle, the newly derived parameters are changed with different factors. For instance, the timescales for the blue and the red lobes are reduced by 0.31 and 0.21, respectively. Larger influences apply to mechanical luminosity, driving force, and mass-loss rate. The comparisons between parameters before and after the correction show that the effect of the inclination-angle cannot be neglected.

Multiple rebrightenings have been observed in the multiband afterglow of GRB 030329. Especially, a marked and quick rebrightening occurred at about t ~ 1.2 * 10^5 s. Energy injection from late and slow shells seems to be the best interpretation for these rebrightenings. Usually it is assumed that the energy is injected into the whole external shock. However, in the case of GRB 030329, the rebrightenings are so quick that the usual consideration fails to give a satisfactory fit to the observed light curves. Actually, since these late/slow shells coast freely in the wake of the external shock, they should be cold and may not expand laterally. The energy injection then should only occur at the central region of the external shock. Considering this effect, we numerically re-fit the quick rebrightenings observed in GRB 030329. By doing this, we were able to derive the beaming angle of the energy injection process. Our result, with a relative residual of only 5% - 10% during the major rebrightening, is better than any previous modeling. The derived energy injection angle is about 0.035. We assume that these late shells are ejected by the central engine via the same mechanism as those early shells that produce the prompt gamma-ray burst. The main difference is that their velocities are much slower, so that they catch up with the external shock very lately and manifest as the observed quick rebrightenings. If this were true, then the derived energy injection angle can give a good measure of the beaming angle of the prompt gamma-ray emission. Our study may hopefully provide a novel method to measure the beaming angle of gamma-ray bursts.

Conventionally, long GRBs are thought to be caused by the core collapses of massive stars. During the lifetime of a massive star, a stellar wind bubble environment should be produced. Furthermore, the microphysics shock parameters may vary along with the evolution of the fireball. Here we investigate the variation of the microphysics shock parameters under the condition of wind bubble environment, and allow the microphysics shock parameters to be discontinuous at shocks in the ambient medium. It is found that our model can acceptably reproduce the rebrightenings observed in GRB afterglows, at least in some cases. The effects of various model parameters on rebrightenings are investigated. The rebrightenings observed in both the R-band and X-ray afterglow light curves of GRB 060206, GRB 070311 and GRB 071010A are reproduced in this model.

Extensive multi-band afterglow data are available for GRB 980703. Especially, its radio afterglow was very bright and was monitored until more than 1000 days after the trigger time. Additionally, there is no obvious special feature, i.e., no rebrightenings, no plateau, and no special steep decay or slow decay in the multi-band afterglow light curves. All these conditions make GRB 980703 a precious sample in gamma-ray burst research. Here we use the observational data of GRB 980703 to test the standard fireball model in depth. It is found that the model can give a satisfactory explanation to the multi-band and overall afterglow light curves. The beaming angle of GRB 980703 is derived as ~0.23 radian, and the circum-burst medium density is ~ 27 cm-3. The total isotropic equivalent kinetic energy of the ejecta is ~ 3.8E52 ergs. A rest-frame extinction of Av ~ 2.5 mag in the host galaxy is also derived.

According to popular progenitor models of gamma-ray bursts, twin jets should be launched by the central engine, with a forward jet moving toward the observer and a receding jet (or the counter jet) moving backwardly. However, in calculating the afterglows, usually only the emission from the forward jet is considered. Here we present a detailed numerical study on the afterglow from the receding jet. Our calculation is based on a generic dynamical description, and includes some delicate ingredients such as the effect of the equal arrival time surface. It is found that the emission from the receding jet is generally rather weak. In radio bands, it usually peaks at a time of $t \geq 1000$ d, with the peak flux nearly 4 orders of magnitude lower than the peak flux of the forward jet. Also, it usually manifests as a short plateau in the total afterglow light curve, but not as an obvious rebrightening as once expected. In optical bands, the contribution from the receding jet is even weaker, with the peak flux being $\sim 8$ orders of magnitude lower than the peak flux of the forward jet. We thus argue that the emission from the receding jet is very difficult to detect. However, in some special cases, i.e., when the circum-burst medium density is very high, or if the parameters of the receding jet is quite different from those of the forward jet, the emission from the receding jet can be significantly enhanced and may still emerge as a marked rebrightening. We suggest that the search for receding jet emission should mostly concentrate on nearby gamma-ray bursts, and the observation campaign should last for at least several hundred days for each event.