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

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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

Methyl carbamate CH$_3$OC(O)NH$_2$ is an isomer of glycine. Quantum chemical analyses show that methyl carbamate is more stable isomer than glycine. Because of this, there could be a higher chance for methyl carbamte to exist in the interstellar medium as compared to glycine. Despite immense searches, till now glycine has not been detected in the ISM, therefore it is worthwhile to search its isomer methyl carbamate. In this paper, we present the constraints of methyl carbamate formation under the interstellar conditions. Large complex organic molecules are favorably produced in hot-corino environments of low mass protostars. We for the first time carried out astrochemical modeling focusing on the formation of methyl carbamate in physical conditions similar to hot-corino objects. Consequently, we examined ALMA archival data for existing spectral line observations toward hot corinos NGC1333 IRAS 4A2 and IRAS 16293B. Within the common spectral range towards these sources, we found three features are possibly related to the spectral transitions of methyl carbamate and consequently estimate the upper limit of column densities. Results of chemical modeling are consistent with the observational upper limit of estimated column density/abundance toward the sources. This may hint the validation of the proposed formation mechanism. Future observations using telescope like ngVLA may confirm the presence of MC toward the hot corinos.

Recent observations of VER J2227+608 reveal a broken power $\gamma$-ray spectrum with the spectral index increasing from $\sim 1.8$ in the GeV energy range to $\sim 2.3$ in the TeV range. Such a spectral break can be attributed to radiative energy loss of energetic electrons in the leptonic scenario for the $\gamma$-ray emission, which, in combination with characteristic age of the nearby pulsar, can be used to constrain magnetic field in the emission region. We show that the radio and X-ray observations can also be explained in such a scenario. In the hadronic scenario, the spectral break can be attributed to diffusion of energetic ions in a turbulent medium and detailed spectral measurement can be used to constrain the diffusion coefficient. These two models, however, predict drastically different spectra above 100 TeV, which will be uncovered with future high-resolution observations, such as LHAASO and/or CTA.

We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

Filaments are very common physical phenomena on the Sun and are often taken as important proxies of solar magnetic activities. The study of filaments has become a hot topic in the space weather research. For a more comprehensive understanding of filaments, especially for an understanding of solar activities of multiple solar cycles, it is necessary to perform a combined multifeature analysis by constructing a data set of multiple solar cycle data. To achieve this goal, we constructed a centennial data set that covers the H$\alpha$ data from five observatories around the world. During the data set construction, we encountered varieties of problems, such as data fusion, accurate determination of the solar edge, classifying data by quality, dynamic threshold, and so on, which arose mainly due to multiple sources and a large time span of data. But fortunately, these problems were well solved. The data set includes seven types of data products and eight types of feature parameters with which we can implement the functions of data searching and statistical analyses. It has the characteristics of better continuity and highly complementary to space observation data, especially in the wavelengths not covered by space observations, and covers many solar cycles (including more than 60 yr of high-cadence data). We expect that this new comprehensive data set as well as the tools will help researchers to significantly speed up their search for features or events of interest, for either statistical or case study purposes, and possibly help them get a better and more comprehensive understanding of solar filament mechanisms.

We study the prospect of using TianQin to detect stellar-mass binary black holes (SBBHs). We estimate the expected detection number as well as the precision of parameter estimation on SBBH inspirals, using five different population models. We note TianQin can possibly detect a few SBBH inspirals with signal to noise ratios greater than 12; lowering the threshold and combining multiple detectors can both boost the detection number. The source parameters can be recovered with good precision for most events above the detection threshold. For example, the precision of the merger time most likely occurs near 1s, making it possible to guide the detection of the ground-based detectors, the precision of the eccentricity $e_0$ most likely occurs near $10^{-4}$, making it possible to distinguish the formation channels, and the precision of the mass parameter is better than $10^{-6}$ in general and most likely occurs near $10^{-7}$. We note, in particular, that for a typical merger event, the error volume is likely to be small enough to contain only the host galaxy, which could greatly help in the study of gravitational wave cosmology and relevant studies through the multimessenger observation.

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 Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU~Orionis-type stars. Dynamic protostellar disks featuring ejections of gaseous clumps and quiescent protoplanetary disks experiencing a close encounter with an intruder star were modeled using the numerical hydrodynamics code FEOSAD. Both the gas and dust dynamics were taken into account, including dust growth and mutual friction between the gas and dust components. Only plane-of-the-disk encounters were considered. Ejected clumps produce a unique type of tail that is characterized by a bow-shock shape. Such tails originate from the supersonic motion of ejected clumps through the dense envelope that often surrounds young gravitationally unstable protostellar disks. The ejected clumps either sit at the head of the tail-like structure or disperse if their mass is insufficient to withstand the head wind of the envelope. On the other hand, close encounters with quiescent protoplanetary disks produce three types of the tail-like structure; we define these as pre-collisional, post-collisional, and spiral tails. These tails can in principle be distinguished from one another by particular features of the gas and dust flow in and around them. We find that the brown-dwarf-mass intruders do not capture circumintruder disks during the encounter, while the subsolar-mass intruders can acquire appreciable circumintruder disks with elevated dust-to-gas ratios, which can ease their observational detection. However, this is true only for prograde collisions; the retrograde intruders fail to collect appreciable amounts of gas or dust from the disk of the target. The predicted mass of dust in the model tail-like structures is higher than what was inferred for similar structures in SU~Aur, FU~Ori, and Z~CMa, making their observational detection feasible. Abridged.

Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We propose a self-similar (SS) solution to hydrodynamic non-relativistic flow behind a spherical strong blast wave (BW) passing through a homogeneous plasma with efficient relativistic particle acceleration at the shock front. The flow is described by an ideal two-fluid model with a relativistic component so that the post-shock gas has an effective SS adiabatic index $ \gamma $ varying from $ 5/3 $ to $ 4/3 $. This solution is calculated numerically and compared with the standard Sedov solution. We find that the BW center in our solution is dominated by the relativistic component with $ \gamma =4/3 $ for the divergence of expansion there, and the relativistic component dominates the interior for a moderate acceleration efficiency at the shock front. The overall efficiency of relativistic particle acceleration can be enhanced by a factor of $ 2 $ due to the slower adiabatic energy loss rate of the relativistic component during expansion. Tendency of the dominance by the relativistic component may be common in expanding astrophysical two-fluid systems such as supernova remnants, lobes of radio galaxies.

We analyzed properties of waves excited by mildly relativistic electron beams propagating along magnetic field with a ring-shape perpendicular momentum distribution in neutral and current-free solar coronal plasmas. These plasmas are subject to both the beam and the electron cyclotron maser (ECM) instabilities driven by the positive momentum gradient of the ring-beam electron distribution in the directions parallel and perpendicular to the ambient magnetic field, respectively. To explore the related kinetic processes self-consistently, 2.5-dimensional fully kinetic particle-in-cell (PIC) simulations were carried out. To quantify excited wave properties in different coronal conditions, we investigated the dependence of their energy and polarization on the ring-beam electron density and magnetic field. In general, electrostatic waves dominate the energetics of waves and nonlinear waves are ubiquitous. In weakly magnetized plasmas, where the electron cyclotron frequency $\omega_{ce}$ is lower than the electron plasma frequency $\omega_{pe}$, it is difficult to produce escaping electromagnetic waves with frequency $\omega > \omega_{pe}$ and small refractive index $|c k / \omega| < 1$ ($k$ and $c$ are the wavenumber and the light speed, respectively). Highly polarized and anisotropic escaping electromagnetic waves can, however, be effectively excited in strongly magnetized plasmas with $\omega_{ce}/\omega_{pe} \geq 1$. The anisotropy of the energy, circular polarization degree (CPD), and spectrogram of these escaping electromagnetic waves strongly depend on the number density ratio of the ring-beam electrons to the background electrons. In particular, their CPDs can vary from left-handed to right-handed with the decrease of the ring-beam density, which may explain some observed properties of solar radio bursts (e.g., radio spikes) from the solar corona.

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $\mu m$ continuum emissions. The angles between outflow axes and corresponding filaments ($\gamma_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

To deepen our understanding of the chemical properties of the Planck Galactic Cold Clump (PGCC) G168.72-15.48, we performed observations of nine molecular species, namely, \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{HC7N}, \ce{SO}, \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH}, toward two dense cores in PGCC G168.72-15.48 using the Tianma Radio Telescope and Purple Mountain Observatory Telescope. We detected \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{N2H+}, \ce{CCH}, and \ce{CH3OH} in both G168-H1 and G168-H2 cores, whereas \ce{HC7N} and \ce{CH3CCH} were detected only in G168-H1 and SO was detected only in G168-H2. Mapping observations reveal that the \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH} emissions are well coupled with the dust emission in G168-H1. Additionally, \ce{N2H+} exhibits an exceptionally weak emission in the denser and more evolved G168-H2 core, which may be attributed to the \ce{N2H+} depletion. We suggest that the \ce{N2H+} depletion in G168-H2 is dominated by \ce{N2} depletion, rather than the destruction by CO. The local thermodynamic equilibrium calculations indicate that the carbon-chain molecules of \ce{CCH}, \ce{HC5N}, \ce{HC7N}, and \ce{CH3CCH} are more abundant in the younger G168-H1 core. We found that starless core G168-H1 may have the properties of cold dark clouds based on its abundances of carbon-chain molecules. While, the prestellar core G168-H2 exhibits lower carbon-chain molecular abundances than the general cold dark clouds. With our gas-grain astrochemical model calculations, we attribute the observed chemical differences between G168-H1 and G168-H2 to their different gas densities and different evolutionary stages.

The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structure of the Galaxy. In order to study this potential link, a statistically significant sample of filaments spread throughout the Galaxy is required. In this work we present the first catalogue of $32,059$ candidate filaments automatically identified in the Hi-GAL survey of the entire Galactic Plane. For these objects we determined morphological (length, $l^{a}$, and geometrical shape) and physical (average column density, $N_{\rm H_{2}}$, and average temperature, $T$) properties. We identified filaments with a wide range of properties: 2$'$\,$\leq l^{a}\leq$\, 100$'$, $10^{20} \leq N_{\rm H_{2}} \leq 10^{23}$\,cm$^{-2}$ and $10 \leq T\leq$ 35\,K. We discuss their association with the Hi-GAL compact sources, finding that the most tenuous (and stable) structures do not host any major condensation and we also assign a distance to $\sim 18,400$ filaments for which we determine mass, physical size, stability conditions and Galactic distribution. When compared to the spiral arms structure, we find no significant difference between the physical properties of on-arm and inter-arm filaments. We compared our sample with previous studies, finding that our Hi-GAL filament catalogue represents a significant extension in terms of Galactic coverage and sensitivity. This catalogue represents an unique and important tool for future studies devoted to understanding the filament life-cycle.

We present a pilot HI survey of 17 Planck Galactic Cold Clumps (PGCCs) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). HI Narrow Self-Absorption (HINSA) is an effective method to detect cold HI being mixed with molecular hydrogen H$_2$ and improves our understanding of the atomic to molecular transition in the interstellar medium. HINSA was found in 58\% PGCCs that we observed. The column density of HINSA was found to have an intermediate correlation with that of $^{13}$CO, following $\rm log( N(HINSA)) = (0.52\pm 0.26) log(N_{^{13}CO}) + (10 \pm 4.1) $. HI abundance relative to total hydrogen [HI]/[H] has an average value of $4.4\times 10^{-3}$, which is about 2.8 times of the average value of previous HINSA surveys toward molecular clouds. For clouds with total column density N$\rm_H >5 \times 10^{20}$ cm$^{-2}$, an inverse correlation between HINSA abundance and total hydrogen column density is found, confirming the depletion of cold HI gas during molecular gas formation in more massive clouds. Nonthermal line width of $^{13}$CO is about 0-0.5 km s$^{-1}$ larger than that of HINSA. One possible explanation of narrower nonthermal width of HINSA is that HINSA region is smaller than that of $^{13}$CO. Based on an analytic model of H$_2$ formation and H$_2$ dissociation by cosmic ray, we found the cloud ages to be within 10$^{6.7}$-10$^{7.0}$ yr for five sources.

We investigate at a high angular resolution the spatial and kinematic structure of the S255IR high mass star-forming region, which demonstrated recently the first disk-mediated accretion burst in the massive young stellar object. The observations were performed with ALMA in Band 7 at an angular resolution $ \sim 0.1^{\prime\prime}$, which corresponds to $ \sim 180 $ AU. The 0.9 mm continuum, C$^{34}$S(7-6) and CCH $N=4-3$ data show a presence of very narrow ($ \sim 1000 $ AU), very dense ($n\sim 10^7$ cm$^{-3}$) and warm filamentary structures in this area. At least some of them represent apparently dense walls around the high velocity molecular outflow with a wide opening angle from the S255IR-SMA1 core, which is associated with the NIRS3 YSO. This wide-angle outflow surrounds a narrow jet. At the ends of the molecular outflow there are shocks, traced in the SiO(8-7) emission. The SiO abundance there is enhanced by at least 3 orders of magnitude. The CO(3-2) and SiO(8-7) data show a collimated and extended high velocity outflow from another dense core in this area, SMA2. The outflow is bent and consists of a chain of knots, which may indicate periodic ejections possibly arising from a binary system consisting of low or intermediate mass protostars. The C$^{34}$S emission shows evidence of rotation of the parent core. Finally, we detected two new low mass compact cores in this area (designated as SMM1 and SMM2), which may represent prestellar objects.

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios.

The spectra fits to a sample of 34 supernova remnants (Zeng et al., 2019) are updated. $\gamma$-ray spectra of 20 supernova remnants (SNRs) with a soft TeV spectrum are further analyzed. We found that 17 of them can be fitted in the hadronic scenario with a single power-law ion distribution with an index of $\sim$ 2.6, which is significantly softer than the ion distribution inferred from $\gamma$-ray observations of star-burst galaxies. If Galactic cosmic rays are mostly produced by SNRs, this result suggests that SNRs in star-burst galaxies may never reach the phase with a soft $\gamma$-ray spectra or escape of high-energy particles from SNRs before they reach this phase with a soft $\gamma$-ray spectrum dominates the contribution of SNRs to Galactic cosmic rays.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Medium Energy X-ray telescope (ME) is one of the three main telescopes on board the Insight Hard X-ray Modulation Telescope (Insight-HXMT) astronomy satellite. ME contains 1728 pixels of Si-PIN detectors sensitive in 5-30 keV with a total geometrical area of 952 cm2. Application Specific Integrated Circuit (ASIC) chips, VA32TA6, is used to achieve low power consumption and low readout noise. The collimators define three kinds of field of views (FOVs) for the telescope, 1{\deg}{\times}4{\deg}, 4{\deg}{\times}4{\deg}, and blocked ones. Combination of such FOVs can be used to estimate the in-orbit X-ray and particle background components. The energy resolution of ME is ~3 keV at 17.8 keV (FWHM) and the time resolution is 255 {\mu}s. In this paper, we introduce the design and performance of ME.

We report the detection of GeV $\gamma$-ray emission from the very-high-energy (VHE) $\gamma$-ray source VER J2227+608 associated with the "tail" region of SNR G106.3+2.7. The GeV $\gamma$-ray emission is extended and spatially coincident with molecular clouds traced by CO emission. The broadband GeV to TeV emission of VER J2227+608 can be well fitted by a single power-law function with an index of 1.90$\pm$0.04, without obvious indication of spectral cutoff toward high energies. The pure leptonic model for the $\gamma$-ray emission can be marginally ruled out by the X-ray and TeV data. In the hadronic model, the low energy content of CRs and the hard $\gamma$-ray spectrum, in combination with the center-bright source structure, suggest that VER J2227+608 may be powered by the PWN instead of shocks of the SNR. And the cutoff energy of the proton distribution needs to be higher than $\sim$ 400 TeV, which makes it an attractive PeVatron candidate. Future observations by the upcoming Large High Altitude Air Shower Observatory (LHAASO) and the Cherenkov Telescope Array in the north (CTA-North) could distinguish these models and constrain the maximum energy of cosmic rays in supernova remnants.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

We report ALMA observations of NGC 1333 IRAS 4A, a young low-mass protostellar binary, referred as 4A1 and 4A2. With multiple H$_2$CO transitions and HNC (4$-$3) observed at a resolution of 0.25" ($\sim$70 au), we investigate the gas kinematics of 4A1 and 4A2. Our results show that, on the large angular scale ($\sim$10"), 4A1 and 4A2 each display a well-collimated outflow along the N-S direction, and an S-shaped morphology is discerned in the outflow powered by 4A2. On the small scale ($\sim$0.3"), 4A1 and 4A2 exhibit distinct spectral features toward the continuum centroid, with 4A1 showing simple symmetric profiles predominantly in absorption and 4A2 demonstrating rather complicated profiles in emission as well as in absorption. Based on radiative transfer modeling exercises, we find that the physical parameters inferred from earlier low-resolution observations cannot be directly extrapolated down to the 4A1 inner region. Possible reasons for the discrepancies between the observed and modelled profiles are discussed. We constrain the mass infall rate in 4A1 to be at most around 3$\times$10$^{-5}$ M$_{\odot}$ year$^{-1}$ at the layer of 75 au. For the kinematics of the 4A2 inner envelope, the absorbing dips in the H$_2$CO spectra are skewed toward the redshifted side and likely signatures of inward motion. These absorbing dips are relatively narrow. This is, like the case for 4A1, significantly slower than the anticipated inflow speed. We estimate a mass infall rate of 3.1$-$6.2 $\times$ 10$^{-5}$ M$_\odot$ year$^{-1}$ at the layer of 100 au in 4A2.

The paper presents a reliable method using deep learning to recognize solar filaments in H-alpha full-disk solar images automatically. This method cannot only identify filaments accurately but also minimize the effects of noise points of the solar images. Firstly, a raw filament dataset is set up, consisting of tens of thousands of images required for deep learning. Secondly, an automated method for solar filament identification is developed using the U-Net deep convolutional network. To test the performance of the method, a dataset with 60 pairs of manually corrected H-alpha images is employed. These images are obtained from the Big Bear Solar Observatory/Full-Disk H-alpha Patrol Telescope (BBSO/FDHA) in 2013. Cross-validation indicates that the method can efficiently identify filaments in full-disk H-alpha images.

Solar Coronal Mass Ejections (CMEs) are sometimes deflected during their propagation. This deflection may be the consequence of interaction between a CME and a coronal hole or the solar wind. We analyze 44 halo-CMEs whose deflection angle exceeds 90 degrees. The coronal magnetic field configuration is computed from daily synoptic maps of magnetic field from SOHO/MDI and SDO/HMI using a Potential Field Source Surface (PFSS) model. By comparing the ambient magnetic field configuration and the measured position angles (MPA) of the CMEs, we conclude that the deflection of 80% of the CMEs (35 of 44) are consistent with the ambient magnetic field configuration, agreeing with previous studies. Of these 35, 71% are deflected toward the heliospheric current sheet (HCS), and 29 degrees toward a pseudo-streamer (PS), the boundary between the same-polarity magnetic field regions. This implies that the ambient coronal magnetic field configuration plays an important and major role in the deflection of CMEs, and that the HCS configuration is more important than PS. If we exclude 13 CMEs having much higher uncertainty from the sample, the agreement between the deflection of CMEs and the ambient field configuration increases substantially, reaching 94% in the new sample of 31 CMEs.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

We present molecular line and dust continuum observations of a Planck-detected cold cloud, G074.11+00.11. The cloud consists of a system of curved filaments and a central star-forming clump. The clump is associated with several infrared sources and H2O maser emission. We aim to determine the mass distribution and gas dynamics within the clump, to investigate if the filamentary structure seen around the clump repeats itself on a smaller scale, and to estimate the fractions of mass contained in dense cores and filaments. The velocity distribution of pristine dense gas can be used to investigate the global dynamical state of the clump, the role of filamentary inflows, filament fragmentation and core accretion. We use molecular line and continuum observations from single dish observatories and interferometric facilities to study the kinematics of the region. The molecular line observations show that the central clump may have formed as a result of a large-scale filament collision. The central clump contains three compact cores. Assuming a distance of 2.3 kpc, based on Gaia observations and a three-dimensional extinction method of background stars, the mass of the central clump exceeds 700 solar masses, which is roughly 25% of the total mass of the cloud. Our virial analysis suggests that the central clump and all identified substructures are collapsing. We find no evidence for small-scale filaments associated with the cores. Our observations indicate that the clump is fragmented into three cores with masses in the range [10,50] solar masses and that all three are collapsing. The presence of an H2O maser emission suggests active star formation. However, the CO lines show only weak signs of outflows. We suggest that the region is young and any processes leading to star formation have just recently begun.

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

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

High-energy particle transport in pulsar wind nebulae (PWNe) plays an essential role in explaining the characteristics revealed in multiwavelength observations. In this paper, the TeV-gamma-ray-emitting electrons in the Vela X PWN are approximated to be injected impulsively when the cocoon is formed due to the interaction between the SNR reverse shock and the PWN. By solving the diffusion-loss equation analytically, we reproduce the broadband spectral energy distribution and surface brightness profile simultaneously. The diffusion coefficient of TeV electrons and positrons, which is well constrained by the spectral and spatial properties of the TeV nebula, is thus determined to be $1 \times 10^{26}$\,cm$^{2}$\,s$^{-1}$ for 10\,TeV electrons and positrons. This coefficient is more than three orders of magnitude lower than that in the interstellar medium, in agreement with a constraint recently obtained from HAWC observations of a TeV nebula associated with the Geminga pulsar. These results suggest that slow diffusion of high-energy particles might be common in PWNe.

The morphological properties of the outflowing circumstellar envelopes surrounding mass-losing stars in eccentric binary systems are presented from a set of three-dimensional hydrodynamical model simulations. Based on four template models of the envelope viewed for a range of inclination angles of the systems, we implement interpretative tools for observations at high spectral/angular resolutions (as illustrated via velocity channel maps as well as position-velocity, radius-velocity, and angle-radius diagrams). Within this framework, the image and kinematical structures can be used to place constraints on the orbital parameters of the system. Specifically, three unique characteristic patterns in the envelopes are found that distinguish these systems from those in binary systems in circular orbits. Bifurcation of the spiral pattern, asymmetry in the interarm density depression, and a concurrent spiral/ring appearance all point to a binary system with an eccentric orbit. The methodology presented in this paper is illustrated in an analysis of recent radio observations of several asymptotic giant branch stars.

We determined the chemical and kinematic properties of the Galactic thin and thick disk using a sample of 307,246 A/F/G/K-type giant stars from the LAMOST spectroscopic survey and Gaia DR2 survey. Our study found that the thick disk globally exhibits no metallicity radial gradient, but the inner disk ($R \le 8$ kpc) and the outer disk ($R>8$ kpc) have different gradients when they are studied separately. The thin disk also shows two different metallicity radial gradients for the inner disk and the outer disk, and has steep metallicity vertical gradient of d[Fe/H]/d$|z|$ $=-0.12\pm0.0007$ dex kpc$^{-1}$, but it becomes flat when it is measured at increasing radial distance, while the metallicity radial gradient becomes weaker with increasing vertical distance. Adopting a galaxy potential model, we derived the orbital eccentricity of sample stars and found a downtrend of average eccentricity with increasing metallicity for the thick disk. The variation of the rotation velocity with the metallicity shows a positive gradient for the thick disk stars and a negative one for the thin disk stars. Comparisons of our observed results with models of disk formation suggest that radial migration could have influenced the chemical evolution of the thin disk. The formation of the thick disk could be affected by more than one processes: the accretion model could play an indispensable role, while other formation mechanisms, such as the radial migration or heating scenario model could also have a contribution.

Non-standard fields are assumed to be responsible for phenomena attributed to dark energy and dark matter. Being coupled to ordinary matter, these fields modify the masses and/or charges of the elementary particles, thereby violating the Weak Equivalence Principle. Thus, values of fundamental constants such as the proton-to-electron mass ratio, mu, and/or the fine structure constant, alpha, measured in different environment conditions can be used as probes for this coupling. Here we perform differential measurements of F = mu*alpha^2 to test a non-standard coupling in the Magellanic Clouds - dwarf galaxies where the overall mass budget is dominated by dark matter. The analysis is based on [CI] and CO lines observed with the Herschel Space Observatory. Since these lines have different sensitivities to changes in mu and alpha, the combined alpha and mu variations can be evaluated through the radial velocity offsets, Delta V, between the CO and [CI] lines. Averaging over nine positions in the Magellanic Clouds, we obtain <Delta V> = -0.02+/-0.07 km/s, leading to |Delta F/F| < 2*10^-7 (1sigma), where Delta F/F = (F_obs-F_lab)/F_lab}. However, for one position observed with five times higher spectral resolution we find Delta V = -0.05+/-0.02 km/s, resulting in Delta F/F = (-1.7+/-0.7)*10^-7. Whether this offset is due to changes in the fundamental constants, due to chemical segregation in the emitting gas or merely due to Doppler noise requires further investigations.

The dependence of polarization fraction $p$ on total intensity $I$ in polarized submillimeter emission measurements is typically parameterized as $p\propto I^{-\alpha}$ $(\alpha \leq 1)$, and used to infer dust grain alignment efficiency in star-forming regions, with an index $\alpha=1$ indicating near-total lack of alignment of grains with the magnetic field. In this work we demonstrate that the non-Gaussian noise characteristics of polarization fraction may produce apparent measurements of $\alpha \sim 1$ even in data with significant signal-to-noise in Stokes $Q$, $U$ and $I$ emission, and so with robust measurements of polarization angle. We present a simple model demonstrating this behavior, and propose a criterion by which well-characterized measurements of polarization fraction may be identified. We demonstrate that where our model is applicable, $\alpha$ can be recovered by fitting the $p-I$ relationship with the mean of the Rice distribution, without statistical debiasing of polarization fraction. We apply our model to JCMT BISTRO Survey POL-2 850$\mu$m observations of three clumps in the Ophiuchus Molecular Cloud, finding that in the externally-illuminated Oph A region, $\alpha\approx 0.34$, while in the more isolated Oph B and C, despite their differing star formation histories, $\alpha \sim 0.6-0.7$. Our results thus suggest that dust grain alignment in dense gas is more strongly influenced by incident interstellar radiation field than by star formation history. We further find that grains may remain aligned with the magnetic field at significantly higher gas densities than has previously been believed, thus allowing investigation of magnetic field properties within star-forming clumps and cores.

How Supermassive Blackholes (SMBHs) are spun-up is a key issue of modern astrophysics. As an extension of the study in Wang et al. (2016), we here address the issue by comparing the host galaxy properties of nearby ($z<0.05$) radio-selected Seyfert 2 galaxies. With the two-dimensional bulge+disk decompositions for the SDSS $r$-band images, we identify a dichotomy on various host galaxy properties for the radio-powerful SMBHs. By assuming the radio emission from the jet reflects a high SMBH spin, which stems from the well-known BZ mechanism of jet production, high-mass SMBHs (i.e., $M_{\mathrm{BH}}>10^{7.9}M_\odot$) have a preference for being spun-up in classical bulges, and low-mass SMBHs (i.e., $M_{\mathrm{BH}}=10^{6-7}M_\odot$) in pseudo-bulges. This dichotomy suggests and confirms that high-mass and low-mass SMBHs are spun-up in different ways, i.e., a major "dry" merger and a secular evolution.

Protoplanets can produce structures in protoplanetary disks via gravitational disk-planet interactions. Once detected, such structures serve as signposts of planet formation. Here we investigate the kinematic signatures in disks produced by multi-Jupiter mass ($M_{\rm J}$) planets using 3D hydrodynamics and radiative transfer simulations. Such a planet opens a deep gap, and drives transonic vertical motions inside. Such motions include both a bulk motion of the entire half-disk column, and turbulence on scales comparable to and smaller than the scale height. They significantly broaden molecular lines from the gap, producing double-peaked line profiles at certain locations, and a kinematic velocity dispersion comparable to thermal after azimuthal averaging. The same planet does not drive fast vertical motions outside the gap, except at the inner spiral arms and the disk surface. Searching for line broadening induced by multi-$M_{\rm J}$ planets inside gaps requires an angular resolution comparable to the gap width, an assessment of the gap gas temperature to within a factor of 2, and a high sensitivity needed to detect line emission from the gap.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Over the past decade, $\gamma$-ray observations of supernova remnants (SNRs) and accurate cosmic-ray (CR) spectral measurements have significantly advanced our understanding of particle acceleration in SNRs. In combination with multiwavelength observations of a large sample of SNRs, it has been proposed that the highest energy particles are mostly accelerated in young remnants, and the maximum energy that middle-age and old SNRs can accelerate particles to decreases rapidly with the decrease in shock speed. If SNRs dominate the CR flux observed at Earth, a large number of particles need to be accelerated in old SNRs for the soft CR spectrum even though they cannot produce very high-energy CRs. With radio, X-ray, and $\gamma$-ray observations of seven middle-age shell-type SNRs, we derive the distribution of high-energy electrons trapped in these remnants via a simple one-zone leptonic emission model and find that their spectral evolution is consistent with such a scenario. In particular, we find that particle acceleration by shocks in middle-age SNRs with an age of $t$ can be described by a unified model with the maximum energy decreasing as $t^{-3.1}$ and the number of GeV electrons increasing as $t^{2.5}$ in the absence of escape from SNRs.

The magnetic field plays a key role in producing solar flares, so that the investigation on the relationship between the magnetic field properties and flares is significant. In this paper, based on the magnetic field extrapolated from the photospheric vector magnetograms of the active region NOAA 9077 obtained at Huairou Solar Observing Station, the magnetic field parameters including the height of field lines, force-free factor, free magnetic energy and inclination angle were studied with respect to an X-class flare in this region. We found that the magnetic field lines became lower and the ratio of number of closed field lines to those of open field lines increased after the flare. The force-free factor ($\alpha$) attained a large value before the flare and then decreased after the flare for the closed field lines, while the open field lines showed the opposite tendency. Free energy reach to maximum before flare, then decrease after flare. The magnetic inclination angles showed opposite change trends after the flare for closed and open field lines. Therefore, we may conclude that non-potential energy released by flare mostly contained in the closed magnetic field lines.

Light bridge (LB) is bright structure crossing the umbra of sunspots and associated to the breakup or assembly of sunspots. In this paper, a LB is presented and studied using the observatory data obtained by {\it Hinode} satellites. Force-free factor ($\alpha$) and the z-component of current ($J_{z}$) and tension force ($T_{z}$) are calculated basing on the vector magnetograms observed by Spectro-Polarimeter (SP) of the Solar Optical Telescope (SOT) on board {\it Hinode}. It is found that the amplitudes of $\alpha$ and $J_{z}$ of LB are generally larger than those of umbra. It is found that there are two signs of $J_{z}$ along LB, which are divided at near the middle position of LB. It is found that the amplitudes of $T_{z}$ of LB are smaller than those of umbra and there are changes of sign of $T_{z}$ between the boundary of LB and umbra. Through comparisons and investigations, it suggest that LB and umbra maybe two different magnetic systems, which is a necessary condition for interaction magnetic reconnection.

Based on the Chinese historical sunspots drawings, a data set consisting of the scanned images and all their digitized parameters from 1925 to 2015 have been constructed. In this paper, we briefly describe the developmental history of sunspots drawings in China. This paper describes the preliminary processing processes that strat from the initial data (inputing to the scanning equipment) to the parameters extraction, and finally summarizes the general features of this dataset. It is the first systematic project in Chinese solar-physics community that the historical observation of sunspots drawings were digitized. Our data set fills in an almost ninety years historical gap, which span 60 degrees from east to west and 50 degrees from north to south and have no continuous and detailed digital sunspot observation information. As a complementary to other sunspots observation in the world, our dataset provided abundant information to the long term solar cycles solar activity research.

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.

HH 212 is one of the well-studied protostellar systems, showing the first vertically resolved disk with a warm atmosphere around the central protostar. Here we report a detection of 9 organic molecules (including newly detected ketene, formic acid, deuterated acetonitrile, methyl formate, and ethanol) in the disk atmosphere, confirming that the disk atmosphere is, for HH 212, the chemically rich component, identified before at a lower resolution as a "hot-corino". More importantly, we report the first systematic survey and abundance measurement of organic molecules in the disk atmosphere within $\sim$ 40 au of the central protostar. The relative abundances of these molecules are similar to those in the hot corinos around other protostars and in Comet Lovejoy. These molecules can be either (i) originally formed on icy grains and then desorbed into gas phase or (ii) quickly formed in the gas phase using simpler species ejected from the dust mantles. The abundances and spatial distributions of the molecules provide strong constraints on models of their formation and transport in star formation. These molecules are expected to form even more complex organic molecules needed for life and deeper observations are needed to find them.

A grid of 20 millions 3-1100$\mu$m SED models is presented for synthetic young clusters embedded in dense clumps. The models depend on four primary parameters: the clump mass M$_{clump}$ and dust temperature T$_{dust}$, the fraction of mass f$_{core}$ locked in dense cores, and the age of the clump t$_{SF}$. We populate the YSO clusters using the IMF from Kroupa(2001) and the YSOs SED models grid of Robitaille et al. (2006). We conduct extensive testing of SED fitting using a simulated dataset and we find that M$_{clump}$ essentially depends on the submillimeter portion of the SED, while T$_{dust}$ is mostly determined from the shape of the SED in the 70-350$\mu$m range. Thanks to the large number of models computed we verify that the combined analysis of L/M, [8-24] and [24-70] colours removes much of the SEDs f$_{core}$-t$_{SF}$ degeneracy. The L/M values are particularly useful to diagnose f$_{core}$. L/M$\leq$1 identifies protoclusters with f$_{core}\leq$0.1 and t$_{SF} \leq 10^5$ years, while L/M$\geq$10 excludes f$_{core}\leq$0.1. We characterize lower limits of L/M where ZAMS stars are not found in models, and we also find models with L/M $\geq$10 and no ZAMS stars, in which [8-24]$\geq0.8\pm 0.1$ independently from M$_{clump}$, temperature and luminosity. This is the first set of synthesis SED models suited to model for embedded and unresolved clusters of YSOs. A set of new evolutionary tracks in the L/M diagram is also presented.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is by far the largest telescope of any kind ever built. FAST produced its first light in September 2016 and it is now under commissioning, with normal operation to commence in late 2019. During testing and early science operation, FAST has started making astronomical discoveries, particularly pulsars of various kinds, including millisecond pulsars, binaries, gamma-ray pulsars, etc. The papers in this mini-volume propose ambitious observational projects to advance our knowledge of astronomy, astrophysics and fundamental physics in many ways. Although it may take FAST many years to achieve all the goals explained in these papers, taken together they define a powerful strategic vision for the next decade.

The gases of the interstellar medium (ISM) possess orders of magnitude more mass than those of all the stars combined and are thus the prime component of the baryonic universe. With L-band surface sensitivity even better than the planned phase one Square-Kilometer-Array (SKA1), the Five-hundred-meter Aperture Spherical radio Telescope (FAST) promises unprecedented insights into two of the primary components of ISM, namely, the atomic hydrogen (HI) and the hydroxyl molecule (OH). We discuss here the evolving landscape of our understanding of ISM, particularly, its complex phases, the magnetic fields within, the so-called dark molecular gas (DMG), high velocity clouds, and the connection between local and distant ISM. We lay out, in broad strokes, several expected FAST projects, including an all northern-sky high-resolution HI survey (22,000 deg$^2$, 3\arcmin\ FWHM beam, 0.2 km/s), targeted OH mapping, searching for absorption or masing signals, and etc. Currently under commissioning, the commensal observing mode of FAST will be capable of simultaneously obtaining HI and pulsar data streams, making large-scale surveys in both science areas more efficient.

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 aim to reveal the physical properties and chemical composition of the cores in the California molecular cloud (CMC), so as to better understand the initial conditions of star formation. We made a high-resolution column density map (18.2") with Herschel data, and extracted a complete sample of the cores in the CMC with the \textsl{fellwalker} algorithm. We performed new single-pointing observations of molecular lines near 90 GHz with the IRAM 30m telescope along the main filament of the CMC. In addition, we also performed a numerical modeling of chemical evolution for the cores under the physical conditions. We extracted 300 cores, of which 33 are protostellar and 267 are starless cores. About 51\% (137 of 267) of the starless cores are prestellar cores. Three cores have the potential to evolve into high-mass stars. The prestellar core mass function (CMF) can be well fit by a log-normal form. The high-mass end of the prestellar CMF shows a power-law form with an index $\alpha=-0.9\pm 0.1$ that is shallower than that of the Galactic field stellar mass function. Combining the mass transformation efficiency ($\varepsilon$) from the prestellar core to the star of $15\pm 1\%$ and the core formation efficiency (CFE) of 5.5\%, we suggest an overall star formation efficiency of about 1\% in the CMC. In the single-pointing observations with the IRAM 30m telescope, we find that 6 cores show blue-skewed profile, while 4 cores show red-skewed profile. [$\rm {HCO}^{+}$]/[HNC] and [$\rm {HCO}^{+}$]/$\rm [N_{2}H^{+}]$ in protostellar cores are higher than those in prestellar cores; this can be used as chemical clocks. The best-fit chemical age of the cores with line observations is $\sim 5\times 10^4$~years.

Cassiopeia A, a well-observed young core-collapse supernova remnant (SNR), is considered as one of the best candidates for studying very high-energy particle acceleration up to PeV via the diffusive shock mechanism. Recently, MAGIC observations revealed a $\gamma$-ray spectral cutoff at $\sim3.5$ TeV, suggesting that if the TeV $\gamma$-rays have a hadronic origin, SNR Cas A can only accelerate particles to tens of TeV. Here, we propose a two-zone emission model for regions associated with the forward (zone 1) and inward/reverse shocks (zone 2). Given the low density in zone 1, it dominates the high-frequency radio emission, soft X-ray rim via the synchrotron process and TeV $\gamma$-ray via the inverse Comptonization. With a relatively softer particle distribution and a higher cut-off energy for electrons, emissions from zone 2 dominate the low-frequency radio, hard X-ray via the synchrotron process and GeV $\gamma$-ray via hadronic processes. There is no evidence for high-energy cutoffs in the proton distributions implying that Cas A can still be a PeVatron. Hadronic processes from zone 1 dominate very high-energy gamma-ray emission. Future observations in hundreds of TeV range can test this model.

We have investigated the linear X-ray filament, G359.944-0.052, previously identified as a likely X-ray counterpart of a parsec-scale jet from the Galactic Center super-massive black hole (SMBH), Sagittarius A* (Sgr A*), using a total of ~5.6 Ms ultra-deep Chandra observations taken from September 1999 to July 2017. This unprecedented dataset enables us to examine flux and spectral variations that might be related to intrinsic properties of the weakly accreting SMBH. We find no flux or spectral variation in G359.944-0.052 after the G2 periapsis passage around early 2014, however, a moderate flux increase of ~2 sigma significance might be associated with the periapsis passage of G1 in early 2001. The filament exhibits an unusually hard spectrum (photon-index < 1) in its portion closest to Sgr A* (i.e., near-side) and a significant spectral softening in the more distant portion, which can be interpreted as synchrotron cooling of the relativistic electrons moving along the jet path. In particular, the hard spectrum of the near-side suggests a piling up of quasi-monoenergetic electrons caused by rapid radiative cooling. The spectral and temporal properties of G359.944-0.052 strengthen the case of it being the X-ray counterpart of a jet launched by Sgr A*.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO ($J$=2--1) line emission from the protoplanetary disk around T-Tauri star SU Aurigae (hereafter SU Aur). Previous observations in optical and near infrared wavelengths find a unique structure in SU Aur. One of the highlights of the observational results is that an extended tail-like structure is associated with the disk, indicating mass transfer from or into the disk. Here we report the discovery of the counterpart of the tail-like structure in CO gas extending more than 1000 au long. Based on geometric and kinematic perspectives, both of the disk and the tail-like structure components physically connect to each other. Several theoretical studies predicted the observed tail-like structure via the following possible scenarios, 1) a gaseous stream from the molecular cloud remnant, 2) collision with a (sub)stellar intruder or a gaseous blob from the ambient cloud, and 3) ejection of a planetary or brown dwarf mass object due to gravitational instability via multi-body gravitational interaction. Since the tail-like structures associated with the SU Aur disk is a new example following RW Aurigae, some disks may experience the internal or external interaction and drastically lose mass during disk evolution.

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.

Based on the second Gaia data release (Gaia DR2) and spectroscopy from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data, we identified 20,089 halo stars kinematically and chemically. The halo streams in the solar neighborhood could be detected in the space of energy and angular momentum. We reshuffle the velocities of these stars to determine the significance of substructure. Finally, we identify 4 statistically significant substructures that are labeled GL-1 through 4. Among these substructures, GL-1 is previously known stream ("N2" stream) and the rest 3 substructures are new. These substructures may be the debris of dwarf galaxies accretion event, their dynamical and chemical information can help to understand the history of the Milky Way.

In this paper, the statistical studies on the force-freeness of photosphere magnetic fields are given. The studies are based on the vector magnetic fields observed by Solar Optical Telescope/Spectro-Polarimeter (SOT/SP) on board Hinode. Three parameters ($F_{x}/ F_{p}$, $F_{y}/ F_{p}$, and $F_{z}/ F_{p}$) are introduced to investigate the force-freeness of active regions photosphere magnetic field. Various thresholds and reductions of original resolutions are selected to calculated parameters. As for the resolutions, the reductions of original resolution data by a factor of 2, 4 and 8 are applied. While for thresholds, they are calculated from individual active region with the average of 75/77/57 G for $B_{x}$/$B_{y}$/$B_{z}$ for all original data. When the resolution are reduced by 2, 4 and 8, the corresponding averages are 72/75/55, 66/68/49 and 57/59/41 G for $B_{x}$/$B_{y}$/$B_{z}$, respectively. It is found that the forces indicated by $F_{x}/ F_{p}$ and $F_{y}/ F_{p}$ increase as thresholds/resolutions increase/decrease. While for $F_{z}/ F_{p}$ the trends become more complex. For low threshold, when the resolution decrease the trends of $F_{z}/ F_{p}$ are similar as those of $F_{x}/ F_{p}$ and $F_{y}/ F_{p}$, while for high threshold the trends of $F_{z}/ F_{p}$ are different from those of $F_{x}/ F_{p}$ and $F_{y}/ F_{p}$. For original resolution, the forces indicated by $F_{y}/ F_{p}$ increase as thresholds increase, while for others resolution they decrease as the thresholds increase contrarily.

The regular solar observations are operated at Huairou Solar Observing Station (HSOS) since 1987, which make the construction of long-term magnetic field datasets available to understand solar magnetic field and cycles. There exist some inconveniences for solar physicist to use these data, because the data storage medium and format at HSOS experienced some changes. Additionally, the processes of magnetic field calibration are not easy to deal with for who are not familiar with these data. Here shows that the magnetic field of HSOS are further processed toward international standards, in order to explore HSOS observations data for scientific research.

The energy storage and trigger mechanisms of solar flares are important for understanding of solar activity. We analyzed multi-wavelength observations of a M1.5 flare on 1 August 2014, in active region NOAA 12127 (SOL2014-08-01T18:13). There are evident large scale sunspot rotations in positive magnetic field of the main energy release region before the eruption; the rotations contain both clockwise and counter-clockwise directions. Injection of magnetic helicity from the photosphere prior to the flare. The sign of the helicity injection is reversed after the flare. It is found that both persistent larger scale ($\approx$ one day) and impulse smaller scale ($\approx$ one to two hours) magnetic-flux emergences are associated with the flare. We conclude that larger-scale flux emergence, helicity injection and sunspot rotation contribute to the energy build up, while the small-scale magnetic-flux emergence plays crucial role in triggering the flare.

The productions of data with internationally agreed standards are the common pursuits for all astronomical observation instruments both ground and spacebased undoubtedly, Since standard data to facilitate the use of scientific research. The routine solar observations are available at Huairou Solar Observing Station (HSOS) since 1987, which should be regarded as one main solar observations in the world. The data storage medium and format at HSOS experienced several changes, so there exist some inconveniences for solar physicist. This paper shows that the observations data of HSOS are further processed both for storage medium and format toward international standards, in order to explore HSOS observations data for scientific research.

Probing magnetic fields in self-gravitating molecular clouds are generally difficult even with the use of the polarimetry. Based on the properties of magneto-hydrodynamic (MHD) turbulence and turbulent reconnection, Velocity Gradient Technique (VGT) provides a new way in tracing magnetic field orientation and strength based on the spectroscopic data. Our study tests the applicability of VGT in various molecular tracers, e.g. 12CO, 13CO, and C18O. By inspecting synthetic molecular line maps of CO isotopologue generated through radiative transfer calculations, we show that the VGT method can be successfully applied in probing the magnetic field direction in the diffuse interstellar medium as well as in self-gravitating molecular clouds.

We report high angular resolution observations of NGC1333 IRAS4A, a protostellar binary including A1 and A2, at 0.84 mm with the Atacama Large Millimeter/submillimeter Array. From the continuum observations, we suggest that the dust emission from the A1 core is optically thick, and A2 is predominantly optically thin. The A2 core, exhibiting a forest of spectral lines including complex molecules, is a well known hot corino as suggested by previous works. More importantly, we report, for the first time, the solid evidence of complex organic molecules (COMs), including CH3OH, ^13CH3OH, CH2DOH, CH3CHO associated with the A1 core seen in absorption. The absorption features mostly arise from a compact region around the continuum peak position of the A1 core. Rather than originating from a larger common envelope surrounding the protobinary, the COM features are associated with individual cores A1 and A2. Considering the signatures observed in both continuum and spectral lines, we propose two possible scenarios for IRAS 4A1 - the COM absorption lines may arise from a hot-corino-like atmosphere at the surface of an optically-thick circumstellar disk around A1, or the absorption may arise from different layers of a temperature-stratified dense envelope.

HL Tau is a Class I-II protostar embedded in an infalling and rotating envelope and possibly associated with a planet forming disk, and it is co-located in a 0.1 pc molecular cloud with two nearby young stellar objects. Our ALMA observations revealed two arc-like structures on a 1000 au scale connected to the disk, and their kinematics could not be explained with any conventional model of infalling and rotational motions. In this work, we investigate the nature of these arc-like structures connected to the HL Tau disk. We conducted new observations in the 13CO and C18O (3-2; 2-1) lines with JCMT and IRAM 30m, and obtained the ACA data with the 7-m array. With the single-dish, ACA, and ALMA data, we analyzed the gas motions on both 0.1 pc and 1000 au scales in the HL Tau region. We constructed new kinematical models of an infalling and rotating envelope with the consideration of relative motion between HL Tau and the envelope. By including the relative motion between HL Tau and its protostellar envelope, our kinematical model can explain the observed velocity features in the arc-like structures. The morphologies of the arc-like structures can also be explained with an asymmetric initial density distribution in our model envelope. In addition, our single-dish results support that HL Tau is located at the edge of a large-scale (0.1 pc) expanding shell driven by the wind or outflow from XZ Tau, as suggested in the literature. The estimated expanding velocity of the shell is comparable to the relative velocity between HL Tau and its envelope in our kinematical model. These results hints that the large-scale expanding motion likely impacts the protostellar envelope around HL Tau and affects its gas kinematics. We found that the mass infalling rate from the envelope onto the HL Tau disk can be decreased by a factor of two due to this impact by the large-scale expanding shell.

The DArk Matter Particle Explorer (DAMPE) developed in China was designed to search for evidence of dark matter particles by observing primary cosmic rays and gamma rays in the energy range from 5 GeV to 10 TeV. Since its launch in December 2015, a large quantity of data has been recorded. With the data set acquired during more than a year of operation in space, a precise time-dependent calibration for the energy measured by the BGO ECAL has been developed. In this report, the instrumentation and development of the BGO Electromagnetic Calorimeter (BGO ECAL) are briefly described. The calibration on orbit, including that of the pedestal, attenuation length, minimum ionizing particle peak, and dynode ratio, is discussed, and additional details about the calibration methods and performance in space are presented.

We present our analysis of the magnetic field structures from 6000 au to 100 au scales in the Class 0 protostar B335 inferred from our JCMT POL-2 observations and the ALMA archival polarimetric data. To interpret the observational results, we perform a series of (non-)ideal MHD simulations of the collapse of a rotating non-turbulent dense core, whose initial conditions are adopted to be the same as observed in B335, and generate synthetic polarization maps. The comparison of our JCMT and simulation results suggests that the magnetic field on a 6000 au scale in B335 is pinched and well aligned with the bipolar outflow along the east-west direction. Among all our simulations, the ALMA polarimetric results are best explained with weak magnetic field models having an initial mass-to-flux ratio of 9.6. However, we find that with the weak magnetic field, the rotational velocity on a 100 au scale and the disk size in our simulations are larger than the observational estimates by a factor of several. An independent comparison of our simulations and the gas kinematics in B335 observed with the SMA and ALMA favors strong magnetic field models with an initial mass-to-flux ratio smaller than 4.8. We discuss two possibilities resulting in the different magnetic field strengths inferred from the polarimetric and molecular-line observations, (1) overestimated rotational-to-gravitational energy in B335 and (2) additional contributions in the polarized intensity due to scattering on a 100 au scale.

Distinguishing diagnostic observational signatures produced by MHD models is essential in understanding the physics for the formation of protostellar disks in the ALMA era. Developing suitable tools along with time evolution will facilitate better identification of diagnostic features. With a ray-tracing based radiative transfer code Perspective, we explore time evolution of MHD models carried out in Li, Krasnopolsky & Shang (2013) - most of which have $90^\circ$ misalignment between the rotational axis and the magnetic field. Four visible object types can be characterized, origins of which are dependent on the initial conditions. Our results show complex spiraling density, velocity and polarization structures. The systems are under constant change, but many of those distinctive features are present already early on, and they grow more visible in time, but most could not be identified from the data without examining their change in time. The results suggest that spiraling pseudodisk structures could function as an effective observation signature of the formation process, and we witness accretion in the disk with eccentric orbits which appear as spiral-like perturbation from simple circular Keplerian orbits. Magnetically aligned polarization appears purely azimuthal in the disk and magnetic field can lead to precession of the disk.

Estimating the cosmic-ray acceleration efficiency $ \epsilon $ in supernova remnants (SNRs) through observations is a challenging task in general. Based on the Rankine-Hugoniot shock conditions, we find an anticorrelation between $ \epsilon $ and the power-law spectral index $ \alpha $ of relativistic particle distribution produced via diffusive particle acceleration by nonrelativistic shocks, implying more efficient acceleration in older SNRs with harder radio spectra. Then $ \epsilon $ may be estimated from some hard radio spectral index measurements. Assuming the particle distribution in downstream of strong shocks to be a nonrelativistic Maxwellian plus a relativistic power law with a high-energy cutoff, we also find that the injection rate for relativistic particles $ \eta $ needs to $ \gtrsim 10^{-6} $ for a prominent decrease of the adiabatic index in SNRs, which implies higher compression ratio and lower values of $ \alpha $. This threshold of $ \eta $ increases with the shock speed $ u_1 $, which may explain the relatively harder radio spectra of older SNRs with lower $ u_1 $. We show that $ \eta $ and/or the relativistic cutoff momentum $ p_\text{m} $ need to be low for old SNRs, and expect a gradual increase of $ \epsilon $ as SNR evolves with gradually decreasing $ \eta $ and $ p_\text{m} $.

Supernova remnants (SNRs) have been considered as the dominant contributors to Galactic cosmic rays. However, the relation between high-energy particles trapped in SNRs and cosmic rays observed at the Earth remains obscure. In this paper, we fit the spectral energy distributions of 35 SNRs with a simple one-zone emission model and analyze correlations of model parameters to uncover the evolution of high-energy particle distribution in SNRs. We find that (1) the particle distribution in general can be described by a broken power-law function with a high-energy cutoff for all SNRs; (2) the low-energy spectrum becomes harder and the break energy decreases with aging of SNRs, (3) for most middle-age SNRs, the energy loss timescale of electrons at the high-energy cutoff is approximately equal to the age of the corresponding remnant implying quenching of very high-energy electron acceleration; for young SNRs, this energy loss timescale is shorter than the age of SNRs implying continuous electon acceleration at the cutoff energy; and for a few old age SNRs, the energy loss timescale is longer than the corresponding age which may suggest escaping of higher energy particles from SNRs. Finally, we comment on the implications of these results on the SNR origin of Galactic cosmic rays.

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

A study on cosmogenic activation in germanium was carried out to evaluate the cosmogenic background level of natural and $^{70}$Ge depleted germanium detectors. The production rates of long-lived radionuclides were calculated with Geant4 and CRY. Results were validated by comparing the simulated and experimental spectra of CDEX-1B detector. Based on the validated codes, the cosmogenic background level was predicted for further tonne-scale CDEX experiment. The suppression of cosmogenic background level could be achieved by underground germanium crystal growth and high-purity germanium detector fabrication to reach the sensitivity requirement for direct detection of dark matter. With the low cosmogenic background, new physics channels, such as solar neutrino research and neutrinoless double-beta decay experiments, were opened and the corresponding simulations and evaluations were carried out.

The discovery of thousands of exoplanets has shown that the birth of planets is a very efficient process in nature. Several physical mechanisms have been proposed to describe the assembly of planets in disks surrounding young stars. However, observational constraints have been sparse on account of insufficient sensitivity and resolution, especially for imaging the inner $10 - 20$ au from the star, where most of planets are observed. Thanks to its unprecedented angular resolution and sensitivity at wavelengths where the emission from the circumstellar material is optically thin, the ngVLA has the potential to transform our understanding of planet formation. In this chapter, we present state-of-the-art theoretical models of planetary systems in the act of forming that demonstrate ngVLA capabilities of imaging and follow the temporal evolution of young solar system analogues up to the distance of the Orion Nebula. These images will unveil how planets form and interact with their parental disks, and will shed light on the diverse properties of exoplanetary systems.

HESS J1640-465 is an extended TeV $\gamma$-ray source and its $\gamma$-ray emission whether from the shell of a supernova remnant (SNR) or a pulsar wind nebula (PWN) is still under debate. We reanalyze the GeV $\gamma$-ray data in the field of HESS J1640-465 using eight years of Pass 8 data recorded by the Fermi Large Area Telescope. An extended GeV $\gamma$-ray source positionally coincident with HESS J1640-465 is found. Its photon spectrum can be described by a power-law with an index of $1.42\pm0.19$ in the energy range of 10-500 GeV, and smoothly connects with the TeV spectrum of HESS J1640-465. The broadband spectrum of HESS J1640-465 can be well fit by a leptonic model with a broken power-law spectrum of electrons with an exponential cut-off at $\sim$ 300 TeV. The spectral properties of HESS J1640-465 are broadly consistent with the characteristics of other sources identified as PWNe, such as the correlations between high-energy luminosity ratios and the physical parameters of pulsar, including spin-down luminosity $\dot{E}$ and characteristic age $\tau_c$. All these pieces of evidence support that the $\gamma$-ray emission of HESS J1640-465 may originate from the PWN powered by PSR J1640-4631 rather than the shell of the SNR G338.3-0.0.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

This paper reports the first results of a direct dark matter search with the DEAP-3600 single-phase liquid argon (LAr) detector. The experiment was performed 2 km underground at SNOLAB (Sudbury, Canada) utilizing a large target mass, with the LAr target contained in a spherical acrylic vessel of 3600 kg capacity. The LAr is viewed by an array of PMTs, which would register scintillation light produced by rare nuclear recoil signals induced by dark matter particle scattering. An analysis of 4.44 live days (fiducial exposure of 9.87 tonne-days) of data taken with the nearly full detector during the initial filling phase demonstrates the detector performance and the best electronic recoil rejection using pulse-shape discrimination in argon, with leakage $<1.2\times 10^{-7}$ (90% C.L.) between 16 and 33 keV$_{ee}$. No candidate signal events are observed, which results in the leading limit on WIMP-nucleon spin-independent cross section on argon, $<1.2\times 10^{-44}$ cm$^2$ for a 100 GeV/c$^2$ WIMP mass (90% C.L.).

The search for telluric extrasolar planets with the Radial Velocity (RV) technique is intrinsically limited by the stellar jitter due to the activity of the star, because stellar surface inhomogeneities, including spots, plages and convective granules, induce perturbations hiding or even mimicking the planetary signal. This kind of noise is poorly understood in all the stars, but the Sun, due to their unresolved surfaces. For these reasons, the effects of the surface inhomogeneities on the measurement of the RV are very difficult to characterize. On the other hand, a better knowledge of these phenomena can allow us a step forward in our understanding of solar and stellar RV noise sources. This will allow to develop more tools for an optimal activity correction leading to more precise stellar RVs. Due to the high spatial resolution with which the Sun is observed, this noise is well known for it. Despite this, a link is lacking between the single observed photospheric phenomena and the behavior of the Sun observed as a star. LOCNES (Low Cost NIR Extended Solar Telescope) will allow to gather time series of RVs in order to disentangle the different contributions to the stellar (i.e., solar) RV jitter. Since July 2015, a Low Cost Solar Telescope (LCST) has been installed outside the TNG dome to feed solar light to the HARPS-N spectrograph (0.38-0.69 $\mu$m; R=115000). The refurbishment of the Near Infrared (NIR) High Resolution Spectrograph GIANO (now GIANO-B) and the new observing mode GIARPS at TNG (simultaneous observations in visible with HARPS-N and in NIR with GIANO-B) is a unique opportunity to extend the wavelength range up to 2.4 $\mu$m for measuring the RV time series of the Sun as a star. This paper outlines the LOCNES project and its scientific drivers.

Temporal photometric variations at near infrared to submillimeter wavelengths have been found in low-mass young stellar objects. These phenomena are generally interpreted as accretion events of star-disk systems with varying accretion rates. There is growing evidence suggesting that similar luminosity flaring also occurs in high-mass star/cluster-forming regions. We report in this Letter the rise and fall of the 900 ${\mu}$m continuum emission and the newly found 349.1 GHz methanol maser emission in the massive star forming region S255IR~SMA1 observed with the Submillimeter Array and the Atacama Large Millimeter/submillimeter Array. The level of flux variation at a factor of $\sim$ 2 at the submillimeter band and the relatively short 2-year duration of this burst suggest that the event is probably similar to those milder and more frequent minor bursts seen in 3D numerical simulations.

We present a visual analytics tool, based on the VisIVO suite, to exploit a combination of all new-generation surveys of the Galactic Plane to study the star formation process of the Milky Way. The tool has been developed within the VIALACTEA project, founded by the 7th Framework Programme of the European Union, that creates a common forum for the major new-generation surveys of the Milky Way Galactic Plane from the near infrared to the radio, both in thermal continuum and molecular lines. Massive volumes of data are produced by space missions and ground-based facilities and the ability to collect and store them is increasing at a higher pace than the ability to analyze them. This gap leads to new challenges in the analysis pipeline to discover information contained in the data. Visual analytics focuses on handling these massive, heterogeneous, and dynamic volumes of information accessing the data previously processed by data mining algorithms and advanced analysis techniques with highly interactive visual interfaces offering scientists the opportunity for in-depth understanding of massive, noisy, and high-dimensional data.

Based on the first Gaia data release and spectroscopy from the LAMOST Data Release 4, we study the kinematics and chemistry of the local halo stars. The halo stars are identified kinematically with a relative speed of at least 220 km s$^{-1}$ with respect to the local standard of rest. In total, 436 halo stars are identified. From this halo sample, 16 high velocity (HiVel) stars are identified. We studied the metallicity and [$\alpha$/Fe] distribution of these HiVel stars. Though most of HiVel stars are metal-poor, there are several stars that have metallicity above $-0.5$ dex. To understand the origin of high velocity stars, we evolve the trajectory of the star backwards along the orbit in our adopted Galaxy potential model to determine the orbital parameters and assess whether the star could have originated in the Galactic center. We found that some high velocity stars could have originated from the Galactic center, but other stars were probably kicked up from the Galactic disk.

We study the metallicity distribution and kinematic properties of 4,680 A/F/G/K-type giant stars with $|z|>$ 5 kpc selected from the LAMOST spectroscopic survey. The metallicity distribution of giant stars with 5 $<|z|\leqslant$ 15 kpc can be described by a three-peak Gaussian model with peaks at [Fe/H] $\sim-0.6\pm0.1$, $-1.2\pm0.3$ and $-2.0\pm0.2$, corresponding to the ratio of 19$\%$, 74$\%$ and 7$\%$, respectively. The $\rm{[\alpha/Fe]}$ is used to associate the three peaks with the thick disk, inner-halo and outer-halo components of the Galaxy. The metallicity distribution of these giant stars, which is fit with Gaussians corresponding to the three components, show a growing fraction of inner-halo component and declining fraction of the thick-disk component with increasing distance from the Galactic plane. Adopting a galaxy potential model, we also derive the orbital parameters of the sample stars, such as orbit eccentricity and rotation velocity. The peak values of derived orbital eccentricity for stars covering different metallicity regions maintain $e\sim$ 0.75, independent of height above the plane, within the range 5$<|z|<$ 15 kpc. By comparing the MDFs of stars in different rotation velocity intervals, we find that the majority of the retrograde stars are more metal-poor than the prograde stars.

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

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

We observed a C-class flare at the outer boundary of the super-penumbra of a sunspot. The flare was triggered by an emerging magnetic bipolar region that was obliquely oriented with respect to the super-penumbral fibrils. The flare started due to the low height magnetic reconnection of emerging magnetic flux with super-penumbral field resulting hot multi-temperature plasma flows in the inverse Evershed flow channel and its overlying atmosphere. The inverse Evershed flows in the chromosphere start from super penumbra towards sunspot that end at the outer boundary of the penumbra. The hot plasma flow towards the sunspot in the inverse Evershed channels show about 10 km s$^{-1}$ higher velocity in H$\alpha$ wavelengths compared to the plasma emissions at various temperatures as seen in different AIA filters. Even though these velocities are about seven times higher than the typical inverse-Evershed flow speeds, the flow is diminished at the outer boundary of the sunspot's penumbra. This suggests that the super-penumbral field lines that carry the inverse Evershed flows, are discontinued at the boundary where the penumbral field lines dive into the sun and these two sets of field lines are completely distinct. The discontinuity in the typical magnetic field and plasma properties at the adjoining of these two sets of field lines further leads the discontinuity in the characteristic magnetoacoustic and Alfv\'en speeds, therefore, stopping the plasma flows further on. The multi-temperature plasma in the inverse Evershed channels exhibits \textbf{possible} longitudinal oscillations initially during the onset of the flare, and later flows towards the sunspot. In the multi-temperature view, the different layers above the flare region have the mixture of supersonic as well as subsonic flows.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

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

The results of a detailed analysis of SMA, VLA, and IRAM observations of the region of massive star formation S255N in CO(2---1), \nh, \nhh, \co and some other lines is presented. Combining interferometer and single-dish data has enabled a more detailed investigation of the gas kinematics in the moleclar core on various spatial scales. There are no signs of rotation or isotropic compression on the scale of the region as whole. The largest fragments of gas ($\approx$0.3 pc) are located near the boundary of the regions of ionized hydrogen S255 and S257. Some smaller-scale fragments are associated with protostellar clumps. The kinetic temperatures of these fragments lie in the range 10---80 K. A circumstellar torus with inner radius R$_{in}$ $\approx$ 8000 AU and outer radius R$_{out}$ 12 000 AU has been detected around the clump SMA1. The rotation profile indicates the existence of a central object with mass $\approx$ 8.5/ sin 2 (i) M$_\odot$ . SMA1 is resolved into two clumps, SMA1---NE and SMA1---SE, whose temperatures are $\approx$150 K and $\approx$25 K, respectively. To all appearances, the torus is involved in the accretion of surrounding gas onto the two protostellar clumps.

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

The Dark matter Experiment using Argon Pulse-shape discrimination (DEAP) has been designed for a direct detection search for particle dark matter using a single-phase liquid argon target. The projected cross section sensitivity for DEAP-3600 to the spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons is $10^{-46}~\rm{cm}^{2}$ for a 100 GeV/$c^2$ WIMP mass with a fiducial exposure of 3 tonne-years. This paper describes the physical properties and construction of the DEAP-3600 detector.

Recent Atacama Large Millimeter and Submillimeter Array (ALMA) observations of the protoplanetary disk around the Herbig Ae star HD 163296 revealed three depleted dust gaps at 60, 100 and 160 au in the 1.3 mm continuum as well as CO depletion in the middle and outer dust gaps. However, no CO depletion was found in the inner dust gap. To examine the planet--disk interaction model, we present results of two-dimensional two fluid (gas + dust) hydrodynamic simulations coupled with three-dimensional radiative transfer simulations. In order to fit the high gas-to-dust ratio of the first gap, we find the Shakura--Sunyaev viscosity parameter $\alpha$ must be very small ($\lesssim 10^{-4}$) in the inner disk. On the other hand, a relatively large $\alpha$ ($\sim 7.5\times 10^{-3}$) is required to reproduce the dust surface density in the outer disk. We interpret the variation of $\alpha$ as an indicator of the transition from an inner dead zone to the outer magnetorotational instability (MRI) active zone. Within $\sim 100$ au, the HD 163296 disk's ionization level is low, and non-ideal magnetohydrodynamic (MHD) effects could suppress the MRI, so the disk can be largely laminar. The disk's ionization level gradually increases toward larger radii, and the outermost disk ($r > 300$ au) becomes turbulent due to MRI. Under this condition, we find that the observed dust continuum and CO gas line emissions can be reasonably fit by three half-Jovian-mass planets (0.46, 0.46 and 0.58 $M_\textrm{J}$) at 59, 105 and 160 au, respectively.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

The conservative Post-Newtonian (PN) Hamiltonian formulation of spinning compact binaries has six integrals of motion including the total energy, the total angular momentum and the constant unit lengths of spins. The manifold correction method can effectively eliminate the integration errors accumulation in a long time. In this paper, the accelerated manifold correction method based on graphics processing unit (GPU) is designed to simulate the dynamic evolution of spinning compact binaries. The feasibility and the efficiency of parallel computation on GPU for spinning compact binaries have been confirmed by various numerical experiments. The numerical comparisons show that the accuracy on GPU execution of manifold corrections method has a good agreement with the execution of codes on merely central processing unit (CPU-based) method. The acceleration ability when the codes are implemented on GPU can increase enormously through the use of shared memory and register optimization techniques without additional hardware costs, implying that the speedup is nearly 13 times as compared with the codes executed on CPU for phase space scan (including orbits). In addition, GPU-accelerated manifold correction method is used to numerically study how dynamics are affected by the spin-induced quadrupole-monopole interaction for black hole binary system.

We describe the general structure of the well known S255IR high mass star forming region, as revealed by our recent ALMA observations. The data indicate a physical relation of the major clumps SMA1 and SMA2. The driving source of the extended high velocity well collimated bipolar outflow is not the most pronounced disk-like SMA1 clump harboring a 20 M$_\odot$ young star (S255 NIRS3), as it was assumed earlier. Apparently it is the less evolved SMA2 clump, which drives the outflow and contains a compact rotating structure (probably a disk). At the same time the SMA1 clump drives another outflow, with a larger opening angle. The molecular line data do not show an outflow from the SMA3 clump (NIRS1), which was suggested by IR studies of this region.

HH212, a nearby (400 pc) object in Orion, is a Class 0 protostellar system with a Keplerian disk and collimated bipolar SiO jets. Deuterated water, HDO and a deuterated complex molecule, methanol (CH2DOH) have been reported in the source. Here, we report the HDCO (deuterated formaldehyde) line observation from ALMA data to probe the inner region of HH212. We compare HDCO line with other molecular lines to understand the possible chemistry and physics of the source. The distribution of HDCO emission suggests it may be associated with the base of the outflow. The emission also shows a rotation but it is not associated with the Keplerian rotation of disk or the rotating infalling envelope, rather it is associated with the outflow as previously seen in C 34 S. From the possible deuterium fractionation, we speculate that the gas phase formation of deuterated formaldehyde is active in the central hot region of the low-mass protostar system, HH212.

We investigate the potential of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Next Generation Very Large Array (ngVLA) to observe substructures in nearby young disks which are due to the gravitational interaction between disk material and planets close to the central star. We simulate the gas and dust dynamics in the disk using the LA-COMPASS hydrodynamical code. We generate synthetic images for the dust continuum emission at sub-millimeter to centimeter wavelengths and simulate ALMA and ngVLA observations. We explore the parameter space of some of the main disk and planet properties that would produce substructures that can be visible with ALMA and the ngVLA. We find that ngVLA observations with an angular resolution of 5 milliarcsec at 3 mm can reveal and characterize gaps and azimuthal asymmetries in disks hosting planets with masses down to $\approx~5 M_{\oplus}$ $\approx 1 - 5$ au from a Solar-like star in the closest star forming regions, whereas ALMA can detect gaps down to planetary masses of $\approx 20~M_{\oplus}$ at 5 au. Gaps opened by super-Earth planets with masses $\approx 5 - 10 M_{\oplus}$ are detectable by the ngVLA in the case of disks with low viscosity ($\alpha \sim 10^{-5}$) and low pressure scale height ($h \approx 0.025$ au at 5 au). The ngVLA can measure the proper motion of azimuthal asymmetric structures associated with the disk-planet interaction, as well as possible circumplanetary disks on timescales as short as one to a few weeks for planets at $1 - 5$ au from the star.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

The low dust temperatures (<14 K) of Planck Galactic Cold Clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. TOP-SCOPE is a joint survey program targeting ~2000 PGCCs in J=1-0 transitions of CO isotopologues and ~1000 PGCCs in 850 micron continuum emisison. The objective of the TOP-SCOPE survey and the joint surveys (SMT 10-m, KVN 21-m and NRO 45-m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, the length and the mean line-mass (M/L) of the G26 filament are ~6200 Msun, ~12 pc and ~500 Msun/pc, respectively. Ten massive clumps including eight starless ones are found along the filament. The most massive Clump as a whole may be still in global collapse while its denser part seems to be undergoing expansion due to outflow feedback. The fragmentation in G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, non-magnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index ($\beta$) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed due to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.

Timing analysis provides information about the dynamics of matter accreting on to neutron stars and black holes, and hence is crucial for studying the physics of the accretion flow around these objects. It is difficult, however, to associate the different variability components with each of the spectral components of the accretion flow. We apply several new methods to two Rossi X-ray Timing Explorer observations of the black hole binary GRS 1915+105 during its heartbeat state to explore the origin of the X-ray variability and the interactions of the accretion-flow components. We offer a promising window into the disc--corona interaction through analysing the formation regions of the disc aperiodic variabilities with different time-scales via comparing the corresponding transition energies of the amplitude-ratio spectra. In a previous paper, we analysed the Fourier power density as a function of energy and frequency to study the origin of the aperiodic variability, and combined that analysis with the phase lag as a function of frequency to derive a picture of the disc--corona interaction in this source. We here, for the first time, investigate the phase lag as a function of energy and frequency, and display some interesting details of the disc--corona interaction. Besides, the results from the shape of amplitude-ratio spectrum and from several other aspects suggest that the quasi-periodic oscillation originates from the corona.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

We present preliminary results of the high resolution $ (0.10^{\prime\prime} \times 0.15^{\prime\prime}) $ observations of the high mass star forming region S255IR with ALMA in several spectral windows from ~ 335 GHz to ~ 350 GHz. The main target lines were C$^{34}$S(7-6), CH$_3$CN($19_K-18_K$), CO(3-2) and SiO(8-7), however many other lines of various molecules have been detected, too. We present sample spectra and maps, discuss briefly the source structure and kinematics. A new, never predicted methanol maser line has been discovered.

The routinely flaring events from \sgras\ trace dynamic, high-energy processes in the immediate vicinity of the supermassive black hole. We statistically study temporal and spectral properties, as well as fluence and duration distributions, of the flares detected by the \chandra\ X-ray Observatory from 1999 to 2012. The detection incompleteness and bias are carefully accounted for in determining these distributions. We find that the fluence distribution can be well characterized by a power-law with a slope of $1.73^{+0.20}_{-0.19}$, while the durations ($\tau$ in seconds) by a log-normal function with a mean $\log(\tau)=3.39^{+0.27}_{-0.24}$ and an intrinsic dispersion $\sigma=0.28^{+0.08}_{-0.06}$. No significant correlation between the fluence and duration is detected. The apparent positive correlation, as reported previously, is mainly due to the detection bias (i.e., weak flares can be detected only when their durations are short). These results indicate that the simple self-organized criticality model has difficulties in explaining these flares. We further find that bright flares usually have asymmetric lightcurves with no statistically evident difference/preference between the rising and decaying phases in terms of their spectral/timing properties. Our spectral analysis shows that although a power-law model with a photon index of $2.0\pm0.4$ gives a satisfactory fit to the joint spectra of strong and weak flares, there is weak evidence for a softer spectrum of weaker flares. This work demonstrates the potential to use statistical properties of X-ray flares to probe their trigger and emission mechanisms, as well as the radiation propagation around the black hole.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

Aims. We aimed at investigating the structure and kinematics of the gaseous disk and outflows around the massive YSO S255 NIRS3 in the S255IR-SMA1 dense clump. Methods. Observations of the S255IR region were carried out with ALMA at two epochs in the compact and extended configurations. Results. We serendipitously detected a new, never predicted, bright maser line at about 349.1 GHz, which most probably represents the CH$_3$OH $14_{1} - 14_{0}$ A$^{- +}$ transition. The emission covers most of the 6.7 GHz methanol maser emission area of almost 1$^{\prime\prime}$ in size and shows a velocity gradient in the same sense as the disk rotation. No variability was found on the time interval of several months. It is classified as Class II maser and probably originates in a ring at a distance of several hundreds AU from the central star.

Stellar feedback from high-mass stars (e.g., H{\sc ii} regions) can strongly influence the surrounding interstellar medium and regulate star formation. Our new ALMA observations reveal sequential high-mass star formation taking place within one sub-virial filamentary clump (the G9.62 clump) in the G9.62+0.19 complex. The 12 dense cores (MM 1-12) detected by ALMA are at very different evolutionary stages, from starless core phase to UC H{\sc ii} region phase. Three dense cores (MM6, MM7/G, MM8/F) are associated with outflows. The mass-velocity diagrams of outflows associated with MM7/G and MM8/F can be well fitted with broken power laws. The mass-velocity diagram of SiO outflow associated with MM8/F breaks much earlier than other outflow tracers (e.g., CO, SO, CS, HCN), suggesting that SiO traces newly shocked gas, while the other molecular lines (e.g., CO, SO, CS, HCN) mainly trace the ambient gas continuously entrained by outflow jets. Five cores (MM1, MM3, MM5, MM9, MM10) are massive starless core candidates whose masses are estimated to be larger than 25 M$_{\sun}$, assuming a dust temperature of $\leq$ 20 K. The shocks from the expanding H{\sc ii} regions ("B" \& "C") to the west may have great impact on the G9.62 clump through compressing it into a filament and inducing core collapse successively, leading to sequential star formation. Our findings suggest that stellar feedback from H{\sc ii} regions may enhance the star formation efficiency and suppress the low-mass star formation in adjacent pre-existing massive clumps.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

We have conducted an interferometric line survey in the 0.8 mm band toward the young high-mass protostar candidate NGC 2264 CMM3 with ALMA. CMM3 is resolved into the two continuum peaks, CMM3A and CMM3B, at an angular separation of 0.9". Thus, CMM3 is found to be a binary system candidate. We have detected molecular outflows associated with CMM3A and CMM3B each, indicating active star formation. In addition to the two peaks, six faint continuum peaks are detected around CMM3A and CMM3B, most of which are thought to be evolved low-mass protostars. CMM3A is found to be rich in molecular line emission including complex organic molecules such as HCOOCH3 and CH3OCH3. The emission of complex organic molecules is distributed within a compact region around the continuum peak of CMM3A. Hence, CMM3A apparently harbors a hot core. On the other hand, CMM3B is deficient in molecular line emission, although its continuum flux is almost comparable to that of CMM3A. Possible origins of the striking difference between CMM3A and CMM3B are discussed.

The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.

Recent precise measurements of cosmic ray (CR) spectra show that the energy distribution of protons is softer than those of heavier nuclei, and there are spectral hardenings for all nuclear compositions above $\sim$200 GV. Models proposed for these anomalies generally assume steady-state solutions of the particle acceleration process. We show that, if the diffusion coefficient has a weak dependence on the particle rigidity near shock fronts of supernova remnants (SNRs), time-dependent solutions of the linear diffusive shock acceleration at two stages of SNR evolution can naturally account for these anomalies. The high-energy component of CRs is dominated by acceleration in the free expansion and adiabatic phases with enriched heavy elements and a high shock speed. The low energy component may be attributed to acceleration by slow shocks propagating in dense molecular clouds with low metallicity in the radiative phase. Instead of a single power-law distribution, the spectra of time-dependent solutions soften gradually with the increase of energy, which may beresponsible for the "knee" of CRs.

We present the results from a series of ground-based radio observations toward a Planck Galactic Cold Clump (PGCC), PGCC G108.84-00.81, which is located in one curved filamentary cloud in the vicinity of an extended HII region Sh2-152 and SNR G109.1-1.0. PGCC G108.84-00.81 is mainly composed of two clumps, "G108-N" and "G108-S". In the 850 micron dust continuum emission map, G108-N is shown as one component while G108-S is fragmented into four components. There is no infrared source associated with G108-N while there are two infrared sources (IRS 1 and IRS 2) associated with G108-S. The total mass of G108-N is larger than the jeans mass, suggesting that G108-N is gravitationally unstable and a potential place for a future star formation. The clump properties of G108-N and G108-S such as the gas temperature and the column density, are not distinctly different. However, G108-S is slightly more evolved than G108-N, in the consideration of the CO depletion factor, molecular abundances, and association with infrared sources. G108-S seems to be affected by the compression from Sh2-152, while G108-N is relatively protected from the external effect

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

Puppis A is a very famous and extensively studied supernova remnant (SNR) that shows strong evidence of shock-cloud interaction. We re-analyze the GeV $\gamma$-ray emission of it using seven years Pass 8 data recorded by the Fermi Large Area Telescope (Fermi-LAT). The morphology of the $\gamma$-ray emission is more compatible with that of the thermal X-ray and IR emissions than the radio image, which suggests a possible correlation between the gamma-ray emitting region and dense clouds. The $\gamma$-ray spectrum in the energy range of 1-500 GeV shows a break at $7.92\pm1.91$ GeV with the photon indices of $1.81\pm0.08$ below the break and $2.53\pm0.12$ above the break, which can naturally explain the lack of TeV $\gamma$-ray emission from Puppis A. The multi-wavelength observations favor a hadronic origin for the $\gamma$-ray emission.

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.

Preplanetary nebulae and planetary nebulae are evolved, mass-losing stellar objects that show a wide variety of morphologies. Many of these nebulae consist of outer structures that are nearly spherical (spiral/shell/arc/halo) and inner structures that are highly asymmetric (bipolar/multipolar). The coexistence of such geometrically distinct structures is enigmatic because it hints at the simultaneous presence of both wide and close binary interactions, a phenomenon that has been attributed to stellar binary systems with eccentric orbits. Here we report new high-resolution molecular-line observations of the circumstellar spiral-shell pattern of AFGL 3068, an asymptotic giant branch star transitioning to the preplanetary nebula phase. The observations clearly reveal that the dynamics of the mass loss is influenced by the presence of an eccentric-orbit binary. This quintessential object opens a new window on the nature of deeply embedded binary stars through the circumstellar spiral-shell patterns that reside at distances of several thousand Astronomical Units from the stars.

We present a detailed study of an earth-directed coronal mass ejection (Full halo CME) event happened on 2011 February 15 making use of white light observations by three coronagraphs and radio observations by Wind/WAVES. We applied three different methods to reconstruct the propagation direction and traveling distance of the CME and its driven shock. We measured the kinematics of the CME leading edge from white light images observed by STEREO A and B, tracked the CME-driven shock using the frequency drift observed by Wind/WAVES together with an interplanetary density model, and obtained the equivalent scattering centers of the CME by Polarization Ratio(PR) method. For the first time, we applied PR method to different features distinguished from LASCO/C2 polarimetric observations and calculated their projections onto white light images observed by STEREO A and B. By combining the graduated cylindrical shell (GCS) forward modeling with the PR method, we proposed a new GCS-PR method to derive 3D parameters of a CME observed from a single perspective at Earth. Comparisons between different methods show a good degree of consistence in the derived 3D results.

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 report the detection of a GeV $\gamma$-ray source which is spatially overlapping and thus very likely associated with the unidentified very-high-energy (VHE) $\gamma$-ray source HESS J1427-608 with the Pass 8 data recorded by the Fermi Large Area Telescope. The photon spectrum of this source is best described by a power-law with an index of $1.85\pm0.17$ in the energy range of $3-500$ GeV, and the measured flux connects smoothly with that of HESS J1427-608 at a few hundred GeV. This source shows no significant extension and time variation. The broadband GeV-TeV emission over four decades of energies can be well fitted by a single power-law function with an index of 2.0, without obvious indication of spectral cutoff toward high energies. Such a result implies that HESS J1427-608 may be a PeV particle accelerator. We discuss possible nature of HESS J1427-608 according to the multi-wavelength spectral fittings. Given the relatively large errors, either a leptonic or a hadronic model can explain the multi-wavelength data from radio to VHE $\gamma$-rays. The inferred magnetic field strength is a few $\mu$G, which is smaller than typical values of supernova remnants (SNRs), and is consistent with some pulsar wind nebulae (PWNe). On the other hand, the flat $\gamma$-ray spectrum is slightly different from typical PWNe but similar to that of some known SNRs.

Multi-wavelength observations of mature supernova remnants (SNRs), especially with recent advances in gamma-ray astronomy, make it possible to constrain energy distribution of energetic particles within these remnants. In consideration of the SNR origin of Galactic cosmic rays and physics related to particle acceleration and radiative processes, we use a simple one-zone model to fit the nonthermal emission spectra of three shell-type SNRs located within 2 degrees on the sky: RX J1713.7-3946, CTB 37B, and CTB 37A. Although radio images of these three sources all show a shell (or half-shell) structure, their radio, X-ray, and gamma-ray spectra are quite different, offering an ideal case to explore evolution of energetic particle distribution in SNRs. Our spectral fitting shows that 1) the particle distribution becomes harder with aging of these SNRs, implying a continuous acceleration process, and the particle distributions of CTB 37A and CTB 37B in the GeV range are harder than the hardest distribution that can be produced at a shock via the linear diffusive shock particle acceleration process, so spatial transport may play a role; 2) the energy loss timescale of electrons at the high-energy cutoff due to synchrotron radiation appears to be always a bit (within a factor of a few) shorter than the age of the corresponding remnant, which also requires continuous particle acceleration; 3) double power-law distributions are needed to fit the spectra of CTB 37B and CTB 37A, which may be attributed to shock interaction with molecular clouds.

We present a timing analysis of two Rossi X-ray Timing Explorer observations of the microquasar GRS 1915+105 during the heartbeat state. The phase-frequency-power maps show that the intermediate-frequency aperiodic X-ray variability weakens as the source softens in the slow rise phase, and when the quasi-periodic oscillation disappears in the rise phase of the pulse of the double-peaked class its sub-harmonic is still present with a hard phase lag. In the slow rise phase, the energy-frequency-power maps show that most of the aperiodic variability is produced in the corona, and may also induce the aperiodic variability observed at low energies from an accretion disk, which is further supported by the soft phase lag especially in the intermediate-frequency range (with a time delay up to 20 ms). In the rise phase of the pulse, the low-frequency aperiodic variability is enhanced significantly and there is a prominent hard lag (with a time delay up to 50 ms), indicating that the variability is induced by extension of the disk toward small radii as implied by the increase in flux and propagates into the corona. However, during the hard pulse of the double-peaked class, the variability shows no significant lag, which may be attributed to an optically thick corona. These timing results are generally consistent with the spectral results presented by Neilsen et al. (2011, 2012) which indicated that the slow rise phase corresponds to a local Eddington limit and the rise phase of the pulse corresponds to a radiation pressure instability in the disk.

Power-law frequency distributions of the peak flux of solar flare X-ray emission have been studied extensively and attributed to a system of self-organized criticality (SOC). In this paper, we first show that, so long as the shape of the normalized light curve is not correlated with the peak flux, the flux histogram of solar flares also follows a power-law distribution with the same spectral index as the power-law frequency distribution of the peak flux, which may partially explain why power-law distributions are ubiquitous in the Universe. We then show that the spectral indexes of the histograms of soft X-ray fluxes observed by GOES satellites in two different energy channels are different: the higher energy channel has a harder distribution than the lower energy channel, which challenges the universal power-law distribution predicted by SOC models and implies a very soft distribution of thermal energy content of plasmas probed by the GOES. The temperature ($T$) distribution, on the other hand, approaches a power-law distribution with an index of 2 for high values of $T$. Application of SOC models to statistical properties of solar flares needs to be revisited.

The VIALACTEA project has a work package dedicated to Tools and Infrastructure and, inside it, a task for the Database and Virtual Observatory Infrastructure. This task aims at providing an infrastructure to store all the resources needed by the, more purposely, scientific work packages of the project itself. This infrastructure includes a combination of: storage facilities, relational databases and web services on top of them, and has taken, as a whole, the name of VIALACTEA Knowledge Base (VLKB). This contribution illustrates the current status of this VLKB. It details the set of data resources put together; describes the database that allows data discovery through VO inspired metadata maintenance; illustrates the discovery, cutout and access services built on top of the former two for the users to exploit the data content.

We observed 146 Galactic clumps in HCN (4-3) and CS (7-6) with the Atacama Submillimeter Telescope Experiment (ASTE) 10-m telescope. A tight linear relationship between star formation rate and gas mass traced by dust continuum emission was found for both Galactic clumps and the high redshift (z>1) star forming galaxies (SFGs), indicating a constant gas depletion time of ~100 Myr for molecular gas in both Galactic clumps and high z SFGs. However, low z galaxies do not follow this relation and seem to have a longer global gas depletion time. The correlations between total infrared luminosities (L_TIR) and molecular line luminosities (L'_mol) of HCN (4-3) and CS (7-6) are tight and sublinear extending down to clumps with LTIR~10^{3} L_sun. These correlations become linear when extended to external galaxies. A bimodal behavior in the LTIR--L'mol correlations was found for clumps with different dust temperature, luminosity-to-mass ratio, and sigmaline/sigmavir. Such bimodal behavior may be due to evolutionary effects. The slopes of LTIR--L'mol correlations become more shallow as clumps evolve. We compared our results with lower J transition lines in wu et al. (2010). The correlations between clump masses and line luminosities are close to linear for low effective excitation density tracers but become sublinear for high effective excitation density tracers for clumps with LTIR larger than LTIR~10^4.5 Lsun. High effective excitation density tracers cannot linearly trace the total clump masses, leading to a sublinear correlations for both Mclump-L'mol and LTIR-L'mol relations.

The merger-driven Gamma-ray Bursts (GRBs) and their associated gravitational wave (GW) radiation, if both successfully detected, have some far-reaching implications, including for instance: (i) The statistical comparison of the physical properties of the short/long-short GRBs with and without GW detection can test the general origin model; (ii) Revealing the physical processes taking place at the central engine; (iii) Measuring the velocity of the Gravitational wave directly/accurately. In this work we discuss these implications in the case of possible association of GW150914/ GBM transient 150914. We compared GBM transient 150914 with other SGRBs and found that such an event {may be} a distinct outlier in some statistical diagrams, possibly due to its specific binary-black-hole merger origin. However, the presence of a "new" group of SGRBs with "unusual" physical parameters is also possible. If the outflow of GBM transient 150914 was launched by the accretion onto the nascent black hole, the magnetic activity rather than the neutrino process is likely responsible for the energy extraction and the accretion disk mass is estimated to be $\sim 10^{-5}~M_\odot$. The GW150914/GBM transient 150914 association, {if confirmed, would} provide the first opportunity to directly measure the GW velocity and its departure from the speed of the light {should be within} a factor of $\sim 10^{-17}$.

Aims. We study a sample of eight B[e] stars with uncertain evolutionary status to shed light on the origin of their circumstellar dust. Methods. We performed a diagnostic analysis on the spectral energy distribution beyond infrared wavelengths, and conducted a census of neighboring region of each target to ascertain its evolutionary status. Results. In comparison to pre-main sequence Herbig stars, these B[e] stars show equally substantial excess emission in the near-infrared, indicative of existence of warm dust, but much reduced excess at longer wavelengths, so the dusty envelopes should be compact in size. Isolation from star-forming regions excludes the possibility of their pre-main sequence status. Six of our targets, including HD 50138, HD 45677, CD-245721, CD-49 3441, MWC 623, and HD 85567, have been previously considered as FS CMa stars, whereas HD 181615/6 and HD 98922 are added to the sample by this work. We argue that the circumstellar grains of these isolated B[e] stars, already evolved beyond the pre-main sequence phase, should be formed in situ. This is in contrast to Herbig stars, which inherit large grains from parental molecular clouds. It has been thought that HD 98922, in particular, is a Herbig star because of its large infrared excess, but we propose it being in a more evolved stage. Because dust condenses out of stellar mass loss in an inside-out manner, the dusty envelope is spatially confined, and anisotropic mass flows, or anomalous optical properties of tiny grains, lead to the generally low line-of-sight extinction toward these stars.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

We present new spectral line observations of the CH3CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array that for the first time measure the disk density, temperature, and rotational velocity with sufficient resolution (0.37", equivalent to ~600 AU) to assess the gravitational stability of the disk through the Toomre-Q parameter. Our observations resolve the central 2000 AU region that shows steeper velocity gradients with increasing upper state energy, indicating an increase in the rotational velocity of the hotter gas nearer the star. Such spin-up motions are characteristics of an accretion flow in a rotationally supported disk. We compare the observed data with synthetic image cubes produced by three-dimensional radiative transfer models describing a thin flared disk in Keplerian motion enveloped within the centrifugal radius of an angular-momentum-conserving accretion flow. Given a luminosity of 1.3x10^4 Lsun, the optimized model gives a disk mass of 1.5 Msun and a radius of 858 AU rotating about a 12.0 Msun protostar with a disk mass accretion rate of 3.9x10^{-5} Msun/yr. Our study finds that, in contrast to some theoretical expectations, the disk is hot and stable to fragmentation with Q > 2.8 at all radii which permits a smooth accretion flow. These results put forward the first constraints on gravitational instabilities in massive protostellar disks, which are closely connected to the formation of companion stars and planetary systems by fragmentation.

With the laws of mass conservation, momentum conservation and energy conservation, incorporating the processes of neutral gas ionization and ion diffusion, we develop a self-consistent model for the bright ribbon --- the most prominent feature in Io's plasma torus. The model parameters are well constrained by earlier {\it in situ} observations with the Galileo and Voyager spacescrafts. Our model calculation indicates that the total power dissipated inside the torus is 3.6 times bigger than the total power transported to Jovian ionosphere via Birkeland current. The power dissipation inside the torus is relatively uniform. Most of the power transportation associated with the Birkeland current, however, is localized near the flux tube of Io. With a height-intergrated conductivity of 0.15$\,$mho in Jovian ionosphere, consistent with earlier aeronomy models, the model gives a reasonable fit to the recent observations of the FUV Io tail on Jupiter. Extra mass loading near Io is required in the model. This excess of mass injection into the plasma torus may explain the Io-correlated energy source revealed by earlier EUV observations with the Voyager spacecrafts. Further investigation of the model implications is warranted.

Supernova remnants (SNRs) are the most attractive candidates for the acceleration sites of Galactic cosmic rays. We report the detection of GeV $\gamma$-ray emission with the Pass 8 events recorded by Fermi Large Area Telescope (Fermi-LAT) in the vicinity of the shell type SNR CTB 37B that is likely associated with the TeV $\gamma-$ray source HESS J1713-381. The photon spectrum of CTB 37B is consistent with a power-law with an index of $1.89\pm0.08$ in the energy range of $0.5-500$ GeV, and the measured flux connects smoothly with that of HESS J1713-381 at a few hundred GeV. No significant spatial extension and time variation are detected. The multi-wavelength data can be well fitted with either a leptonic model or a hadronic one. However, parameters of both models suggest more efficient particle acceleration than typical SNRs. Meanwhile, the X-ray and $\gamma$-ray spectral properties of CTB 37B show that it is an interesting source bridging young SNRs dominated by non-thermal emission and old SNRs interacting with molecular clouds.

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.

The VIALACTEA project aims at building a predictive model of star formation in our galaxy. We present the innovative integrated framework and the main technologies and methodologies to reach this ambitious goal.

This paper describes several applications in astronomy and cosmology that are addressed using probabilistic modelling and statistical inference.

Mini-Neptunes and volatile-poor super-Earths coexist on adjacent orbits in proximity to host stars such as Kepler-36 and Kepler-11. Several post-formation processes have been proposed for explaining the origin of the compositional diversity: the mass loss via stellar XUV irradiation, degassing of accreted material, and in-situ accumulation of the disk gas. Close-in planets are also likely to experience giant impacts during the advanced stage of planet formation. This study examines the possibility of transforming volatile-rich super-Earths / mini-Neptunes into volatile-depleted super-Earths through giant impacts. We present the results of three-dimensional giant impact simulations in the accretionary and disruptive regimes. Target planets are modeled with a three-layered structure composed of an iron core, silicate mantle and hydrogen/helium envelope. In the disruptive case, the giant impact can remove most of the H/He atmosphere immediately and homogenize the refractory material in the planetary interior. In the accretionary case, the planet can retain more than half of the gaseous envelope, while a compositional gradient suppresses efficient heat transfer as its interior undergoes double-diffusive convection. After the giant impact, a hot and inflated planet cools and contracts slowly. The extended atmosphere enhances the mass loss via both a Parker wind induced by thermal pressure and hydrodynamic escape driven by the stellar XUV irradiation. As a result, the entire gaseous envelope is expected to be lost due to the combination of those processes in both cases. We propose that Kepler-36b may have been significantly devolatilized by giant impacts, while a substantial fraction of Kepler-36c's atmosphere may remain intact. Furthermore, the stochastic nature of giant impacts may account for the large dispersion in the mass--radius relationship of close-in super-Earths and mini-Neptunes.

CIT 6 is a carbon star in the transitional phase from the asymptotic giant branch (AGB) to the protoplanetary nebulae (pPN). Observational evidences of two point sources in the optical, circumstellar arc segments in an HC$_3$N line emission, and a bipolar nebula in near-infrared provide strong support for the presence of a binary companion. Hence, CIT 6 is very attractive for studying the role of companions in the AGB-pPN transition. We have carried out high resolution $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ observations of CIT 6 with the Submillimeter Array combined with the Submillimeter Telescope (single-dish) data. The $^{12}$CO channel maps reveal a spiral-shell pattern connecting the HC$_3$N segments in a continuous form, and an asymmetric outflow corresponding to the near-infrared bipolar nebula. Rotation of the $^{12}$CO channel peak position may be related to the inner spiral winding and/or the bipolar outflow. An eccentric orbit binary is suggested for the presences of an anisotropic mass loss to the west and a double spiral pattern. The lack of interarm emission to the west may indicate a feature corresponding to the periastron passage of a highly eccentric orbit of the binary. Spatially-averaged radial and spectral profiles of $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ are compared with simple spherical radiative transfer models, suggesting a change of $^{12}$CO/$^{13}$CO abundance ratio from $\sim30$ to $\sim50$ inward in the CSE of CIT 6. The millimeter continuum emission is decomposed into extended dust thermal emission (spectral index $\sim-2.4$) and compact emission from radio photosphere (spectral index $\sim-2.0$).

The Large High Altitude Air Shower Observatory (LHAASO) is designed for high energy gamma ray and cosmic ray detection. A Water Cherenkov Detector Array which is sensitive to gamma ray showers above a few hundred GeV is proposed to survey gamma ray sources. The WCDA consists of 3600 PhotoMultiplier Tubes (PMT) which collect the Cherenkov light produced by the shower particles in water. Both high precision time and charge measurement are required over a large dynamic range from 1 photo electron (P.E.) to 4000 P.E. Prototype of an analogue front-end Application Specific Integrated Circuit (ASIC) fabricated in Chartered 0.35 {\mu}m CMOS technology is designed to read out PMT signal in the WCDA. This ASIC employs leading edge discrimination and RC4 shaping structure; combined with the following Time-to-Digital Converter (TDC) and Analog-to-Digital Converter (ADC), both the arrival time and charge of the PMT signal can be measured. Initial test results indicate that time resolution is better than 350 ps and charge resolution is better than 10% at 1 P.E. and better than 1% with large input signals (300 P.E. to 4000 P.E.). Besides, this ASIC has a small channel-to-channel crosstalk and low ambient temperature dependency.

With the accurate cosmic ray (CR) electron and positron spectra (denoted as $\Phi_{\rm e^{-}}$ and $\Phi_{\rm e^{+}}$, respectively) measured by AMS-02 collaboration, the difference between the electron and positron fluxes (i.e., $\Delta \Phi=\Phi_{\rm e^{-}}-\Phi_{\rm e^{+}}$), dominated by the propagated primary electrons, can be reliably inferred. In the standard model, the spectrum of propagated primary CR electrons at energies $\geq 30$ GeV softens with the increase of energy. The absence of any evidence for such a continuous spectral softening in $\Delta \Phi$ strongly suggests a significant `excess' of primary CR electrons and at energies of $100-400$ GeV the identified excess component has a flux comparable to that of the observed positron excess. Middle-age but `nearby' supernova remnants (e.g., Monogem and Geminga) are favored sources for such an excess.

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

We report the results of our observations of the S255IR area with the SMA at 1.3 mm in the very extended configuration and at 0.8 mm in the compact configuration as well as with the IRAM-30m at 0.8 mm. The best achieved angular resolution is about 0.4 arcsec. The dust continuum emission and several tens of molecular spectral lines are observed. The majority of the lines is detected only towards the S255IR-SMA1 clump, which represents a rotating structure (probably disk) around the young massive star. The achieved angular resolution is still insufficient for conclusions about Keplerian or non-Keplerian character of the rotation. The temperature of the molecular gas reaches 130-180 K. The size of the clump is about 500 AU. The clump is strongly fragmented as follows from the low beam filling factor. The mass of the hot gas is significantly lower than the mass of the central star. A strong DCN emission near the center of the hot core most probably indicates a presence of a relatively cold ($\lesssim 80$ K) and rather massive clump there. High velocity emission is observed in the CO line as well as in lines of high density tracers HCN, HCO+, CS and other molecules. The outflow morphology obtained from combination of the SMA and IRAM-30m data is significantly different from that derived from the SMA data alone. The CO emission detected with the SMA traces only one boundary of the outflow. The outflow is most probably driven by jet bow shocks created by episodic ejections from the center. We detected a dense high velocity clump associated apparently with one of the bow shocks. The outflow strongly affects the chemical composition of the surrounding medium.

Spectral line survey observations are conducted toward the high-mass protostar candidate NGC 2264 CMM3 in the 4 mm, 3 mm, and 0.8 mm bands with the Nobeyama 45 m telescope and the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope. In total, 265 emission lines are detected in the 4 mm and 3 mm bands, and 74 emission lines in the 0.8 mm band. As a result, 36 molecular species and 30 isotopologues are identified. In addition to the fundamental molecular species, many emission lines of carbon-chain molecules such as HC5N, C4H, CCS, and C3S are detected in the 4 mm and 3 mm bands. Deuterated molecular species are also detected with relatively strong intensities. On the other hand, emission lines of complex organic molecules such as HCOOCH3, and CH3OCH3 are found to be weak. For the molecules for which multiple transitions are detected, rotation temperatures are derived to be 7-33 K except for CH3OH. Emission lines with high upper-state energies (Eu > 150 K) are detected for CH3OH, indicating existence of a hot core. In comparison with the chemical composition of the Orion KL, carbon-chain molecules and deuterated molecules are found to be abundant in NGC 2264 CMM3, while sulfur-bearing species and complex organic molecules are deficient. These characteristics indicate chemical youth of NGC 2264 CMM3 in spite of its location at the center of the cluster forming core, NGC 2264 C.

With the standard deviation for the logarithm of the re-scaled range $\langle |F(t+\tau)-F(t)|\rangle$ of simulated fractal Brownian motions $F(t)$ given in a previous paper \cite{q14}, the method of least squares is adopted to determine the slope, $S$, and intercept, $I$, of the log$(\langle |F(t+\tau)-F(t)|\rangle)$ vs $\rm{log}(\tau)$ plot to investigate the limitation of this procedure. It is found that the reduced $\chi^2$ of the fitting decreases with the increase of the Hurst index, $H$ (the expectation value of $S$), which may be attributed to the correlation among the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters $S$ and $I$ are usually smaller than their corresponding standard deviations. These results show the limitation of using the simple least square method to determine the dimension of a fractal time series. Nevertheless, they may be used to reinterpret the fitting results of the least square method to determine the dimension of fractal Brownian motions more self-consistently. The currency exchange rate between Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 1.511 is obtained for spans greater than 30 transactions.

Our analysis of archival VLBI data of PSR 0834+06 revealed that its scintillation properties can be precisely modelled using the inclined sheet model (Pen & Levin 2014), resulting in two distinct lens planes. These data strongly favour the grazing sheet model over turbulence as the primary source of pulsar scattering. This model can reproduce the parameters of the observed diffractive scintillation with an accuracy at the percent level. Comparison with new VLBI proper motion results in a direct measure of the ionized ISM screen transverse velocity. The results are consistent with ISM velocities local to the PSR 0834+06 sight-line (through the Galaxy). The simple 1D structure of the lenses opens up the possibility of using interstellar lenses as precision probes for pulsar lens mapping, precision transverse motions in the ISM, and new opportunities for removing scattering to improve pulsar timing. We describe the parameters and observables of this double screen system. While relative screen distances can in principle be accurately determined, a global conformal distance degeneracy exists that allows a rescaling of the absolute distance scale. For PSR B0834+06, we present VLBI astrometry results that provide (for the fist time) a direct measurement of the distance of the pulsar. For targets where independent distance measurements are not available, which are the cases for most of the recycled millisecond pulsars that are the targets of precision timing observations, the degeneracy presented in the lens modelling could be broken if the pulsar resides in a binary system.

We numerically study the detailed evolutionary features of the wave-like disturbance and its propagation in the eruption. This work is a follow-up to Wang et al., using significantly upgraded new simulations. We focus on the contribution of the velocity vortices and the fast shock reflection and refraction in the solar corona to the formation of the EUV waves. Following the loss of equilibrium in the coronal magnetic structure, the flux rope exhibits rapid motions and invokes the fast-mode shock forward of the rope, which then produces the type II radio burst. The expansion of the fast shock, which is associated with outward motion, takes place in various directions, and the downward expansion shows the reflection and the refraction as a result of the non-uniform background plasma. The reflected component of the fast shock propagates upward and the refracted component propagates downward. As the refracted component reaches the boundary surface, a weak echo is excited. The Moreton wave is invoked as the fast shock touches the bottom boundary, so the Moreton wave lags the type II burst. A secondary echo occurs in the area where reflection of the fast shock encounters the slow-mode shock, and the nearby magnetic field lines are further distorted because of the interaction between the secondary echo and the velocity vortices. Our results indicate that the EUV wave may arise from various processes that are revealed in the new simulations.

Aims: We investigate the electron acceleration in convective electric fields of cascading magnetic reconnection in a flaring solar corona and show the resulting hard X-ray (HXR) radiation spectra caused by Bremsstrahlung for the coronal source. Methods: We perform test particle calculation of electron motions in the framework of a guiding center approximation. The electromagnetic fields and their derivatives along electron trajectories are obtained by linearly interpolating the results of high-resolution adaptive mesh refinement (AMR) MHD simulations of cascading magnetic reconnection. Hard X-ray (HXR) spectra are calculated using an optically thin Bremsstrahlung model. Results: Magnetic gradients and curvatures in cascading reconnection current sheet accelerate electrons: trapped in magnetic islands, precipitating to the chromosphere and ejected into the interplanetary space. The final location of an electron is determined by its initial position, pitch angle and velocity. These initial conditions also influence electron acceleration efficiency. Most of electrons have enhanced perpendicular energy. Trapped electrons are considered to cause the observed bright spots along coronal mass ejection CME-trailing current sheets as well as the flare loop-top HXR emissions.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

There is a growing appreciation that the environmental conditions that we call space weather impact the technological infrastructure that powers the coupled economies around the world. With that comes the need to better shield society against space weather by improving forecasts, environmental specifications, and infrastructure design. [...] advanced understanding of space weather requires a coordinated international approach to effectively provide awareness of the processes within the Sun-Earth system through observation-driven models. This roadmap prioritizes the scientific focus areas and research infrastructure that are needed to significantly advance our understanding of space weather of all intensities and of its implications for society. Advancement of the existing system observatory through the addition of small to moderate state-of-the-art capabilities designed to fill observational gaps will enable significant advances. Such a strategy requires urgent action: key instrumentation needs to be sustained, and action needs to be taken before core capabilities are lost in the aging ensemble. We recommend advances through priority focus (1) on observation-based modeling throughout the Sun-Earth system, (2) on forecasts more than 12 hrs ahead of the magnetic structure of incoming coronal mass ejections, (3) on understanding the geospace response to variable solar-wind stresses that lead to intense geomagnetically-induced currents and ionospheric and radiation storms, and (4) on developing a comprehensive specification of space climate, including the characterization of extreme space storms to guide resilient and robust engineering of technological infrastructures. The roadmap clusters its implementation recommendations by formulating three action pathways, and outlines needed instrumentation and research programs and infrastructure for each of these. [...]

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

We propose a valid scheme to measure the Hubble parameter $H(z)$ at high redshifts by detecting the Sandage-Loeb signal (SL signal) which can be realized by the next generation extremely large telescope. It will largely extend the current observational Hubble parameter data (OHD) towards the redshift region of $z \in [2.0,5.0]$, the so-called "redshift desert", where other dark energy probes are hard to provide useful information of the cosmic expansion. Quantifying the ability of this future measurement by simulating observational data for a CODEX (COsmic Dynamics and EXo-earth experiment)-like survey and constraining various cosmological models, we find that the SL signal scheme brings the redshift upper-limit of OHD from $z_\mathrm{max}=2.3$ to $z_\mathrm{max}\simeq 5.0$, provides more accurate constraints on different dark energy models, and greatly changes the degeneracy direction of the parameters. For the $\Lambda$CDM case, the accuracy of $\Omega_m$ is improved by $58\%$ and the degeneracy between $\Omega_m$ and $\Omega_ {\Lambda}$ is rotated to the vertical direction of $\Omega_k = 0$ line strongly; for the $w$CDM case, the accuracy of $w$ is improved by $15\%$. The Fisher matrix forecast on different time-dependent $w(z)$ is also performed.

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

We consider the origin of compact, short-period, Jupiter-mass planets. We propose that their diverse structure is caused by giant impacts of embryos and super-Earths or mergers with other gas giants during the formation and evolution of these hot Jupiters. Through a series of numerical simulations, we show that typical head-on collisions generally lead to total coalescence of impinging gas giants. Although extremely energetic collisions can disintegrate the envelope of gas giants, these events seldom occur. During oblique and moderately energetic collisions, the merger products retain higher fraction of the colliders' cores than their envelopes. They can also deposit considerable amount of spin angular momentum to the gas giants and desynchronize their spins from their orbital mean motion. We find that the oblateness of gas giants can be used to infer the impact history. Subsequent dissipation of stellar tide inside the planets' envelope can lead to runaway inflation and potentially a substantial loss of gas through Roche-lobe overflow. The impact of super-Earths on parabolic orbits can also enlarge gas giant planets' envelope and elevates their tidal dissipation rate over $\sim $ 100 Myr time scale. Since giant impacts occur stochastically with a range of impactor sizes and energies, their diverse outcomes may account for the dispersion in the mass-radius relationship of hot Jupiters.

Aims. Two-dimensional MHD simulations are used to model the emission properties of TeV-bright shell-type supernova remnants (SNRs) and to explore their nature. Methods. In the leptonic scenario for the TeV emission, the $\gamma$-ray emission is produced via Inverse Compton scattering of background soft photons by high-energy electrons accelerated by the shocks of the SNRs. The TeV emissivity is proportional to the magnetic field energy density and MHD simulations can be used to model the TeV structure of such remnants directly. 2D MHD simulations for SNRs are then performed under the assumption that the ambient interstellar medium is turbulent with the magnetic field and density fluctuations following a Kolmogorov-like power-law spectrum. Results. (1) As expected, these simulations confirm early 1D and 2D modelings of these sources, namely the hydrodynamical evolution of the shock waves and amplification of magnetic field by Rayleigh-Taylor convective flows and by shocks propagating in a turbulent medium; (2) We reproduce rather complex morphological structure for $\gamma$-rays, suggesting intrinsic variations of the source morphology not related to the structure of the progenitor and environment; (3)The observed radial profile of several remnants are well reproduced with an ambient medium density of $0.1-1$ cm$^{-3}$. An even lower ambient density leads to a sharper drop of the TeV brightness with radius than what is observed near the outer edge of these remnants. Conclusions. In a turbulent background medium, we can reproduce the observed characteristics of several shell-type TeV SNRs with reasonable parameters except for a higher ambient density than that inferred from X-ray observations.

We carried out detailed analysis of X-ray and radio observations of a simple flare loop that occurred on 12th August 2002, with the impulsive hard X-ray (HXR) light curves dominated by a single pulse. The emission spectra of the early impulsive phase are consistent with an isothermal model in the coronal loop with a temperature reaching several keVs. A power-law high-energy spectral tail is evident near the HXR peak time, in accordance with the appearance of footpoints at high energies, and is well correlated with the radio emission. The energy content of the thermal component keeps increasing gradually after the disappearance of this nonthermal component. These results suggest that electron acceleration only covers a central period of a longer and more gradual energy dissipation process and that the electron transport within the loop plays a crucial role in the formation of the inferred power-law electron distribution. The spectral index of power-law photons shows a very gradual evolution indicating a quasi-steady state of the electron accelerator, which is confirmed by radio observations. These results are consistent with the theory of stochastic electron acceleration from a thermal background. Advanced modeling with coupled electron acceleration and spatial transport processes is needed to explain these observations more quantitatively, which may reveal the dependence of the electron acceleration on the spatial structure of the acceleration region.

The inorganic CsI(Tl) crystal scintillator is a candidate served as an anti-compton detector for China Dark matter EXperiment (CDEX). Studying the intrinsic radiopurity of CsI(Tl) crystal is an issue of major importance. The timing, energy and spatial correlations, as well as the capability of pulse shape discrimination provide powerful methods for the measurement of intrinsic radiopurities. The experimental design, detector performance and event-selection algorithms are described. A total of 359*3 kg-days data from three prototypes of CsI(Tl) crystals were taken at China Jinping underground laboratory where offers a good shielding environment. The contamination levels of internal isotopes from Cs-137, Th-232 and U-238 series, as well as the upper bounds of U-235 series are reported. Identification of the whole alpha8 peaks from U/Th decay chains and derivation of those corresponding quenching factors are achieved.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

We investigate equilibrium height of the flux rope, and its internal equilibrium in a realistic plasma environment by carrying out numerical simulations of the evolution of systems including a current-carrying flux rope. We find that the equilibrium height of the flux rope is approximately a power-law function of the relative strength of the background field. Our simulations indicate that the flux rope can escape more easily from a weaker background field. This further confirms the catastrophe in the magnetic configuration of interest can be triggered by decrease of strength of the background field. Our results show that it takes some time to reach internal equilibrium depending on the initial state of the flux rope. The plasma flow inside the flux rope due to the adjustment for the internal equilibrium of the flux rope remains small and does not last very long when the initial state of the flux rope commences from the stable branch of the theoretical equilibrium curve. This work also confirms the influence of the initial radius of flux rope on its evolution, the results indicate that the flux rope with larger initial radius erupts more easily. In addition, by using the realistic plasma environment and much higher resolution in our simulations, we notice some different characters compared to previous studies in Forbes (1990).

The Sky Brightness Monitor (SBM) is an important instrument to measure the brightness level for the sky condition, which is a critical parameter for judging a site for solar coronal observations. In this paper we present an automatic method for the processing of SBM data in large quantity, which can separate the regions of the Sun and the nearby sky as well as recognize the regions of the supporting arms in the field of view. These processes are implemented on the data acquired by more than one SBM instruments during our site survey project in western China. An analysis applying the result from our processes has been done for the assessment of the scattered-light levels by the instrument. Those results are considerably significant for further investigations and studies, notably to derive a series of the other important atmospheric parameters such as extinctions, aerosol content and precipitable water vapor content for candidate sites. Our processes also provide a possible way for full-disk solar telescopes to track the Sun without an extra guiding system.

Using $5.4$ year Fermi-LAT data, we report the detection of GeV $\gamma$-ray emission from the shell-type supernova remnant RCW 86 (G315.4-2.3) with a significance of $\sim5.1\sigma$. The data slightly favors an extended emission of this supernova remnant. The spectral index of RCW 86 is found to be very hard, $\Gamma\sim1.4$, in the $0.4$ to $300$ GeV range. A one zone leptonic model can well fit the multi-wavelength data from radio to very high energy $\gamma$-rays. The very hard GeV $\gamma$-ray spectrum and the inferred low gas density seem to disfavor the hadronic origin of the $\gamma$-rays. The $\gamma$-ray behavior of RCW 86 is very similar to several other TeV shell-type supernova remnants, e.g., RX J1713.7-3946, RX J0852.0-4622, SN 1006 and HESS J1731-347.

The DAMPE (DArk Matter Particle Explorer) is a scientific satellite being developed in China, aimed at cosmic ray study, gamma ray astronomy, and searching for the clue of dark matter particles, with a planned mission period of more than 3 years and an orbit altitude of about 500 km. The BGO Calorimeter, which consists of 308 BGO (Bismuth Germanate Oxid) crystal bars, 616 PMTs (photomultiplier tubes) and 1848 dynode signals, has approximately 32 radiation lengths. It is a crucial sub-detector of the DAMPE payload, with the functions of precisely measuring the energy of cosmic particles from 5 GeV to 10TeV, distinguishing positrons/electrons and gamma rays from hadron background, and providing trigger information for the whole DAMPE payload. The dynamic range for a single BGO crystal is about 2?105 and there are 1848 detector signals in total. To build such an instrument in space, the major design challenges for the readout electronics come from the large dynamic range, the high integrity inside the very compact structure, the strict power supply budget and the long term reliability to survive the hush environment during launch and in orbit. Currently the DAMPE mission is in the end of QM (Qualification Model) stage. This paper presents a detailed description of the readout electronics for the BGO calorimeter.

By adopting the differential age method, we utilize selected 17832 luminous red galaxies (LRGs) from Sloan Digital Sky Survey Data Release Seven (SDSS DR7) covering redshift $0<z<0.4$ to measure Hubble parameters. Using a full spectrum fitting package UlySS, these spectra are reduced with single stellar population (SSP) models and optimal age information of our selected sample are derived. With the decreasing age-redshift relation, four new observational $H(z)$ data (OHD) points are obtained, which are $H(z)=69.0\pm19.6$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.07$, $H(z)=68.6\pm26.2$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.12$, $H(z)$=$72.9\pm29.6$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.2$ and $H(z)$=$88.8\pm36.6$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.28$, respectively. Combined with other 21 available OHD data points, a performance of constraint on both flat and non-flat $\Lambda$CDM model is presented.

Context: HESS J1731-347 has been identified as one of the few TeV-bright shell-type supernova remnants (SNRs). These remnants are dominated by nonthermal emission, and the nature of TeV emission has been continuously debated for nearly a decade. Aims: We carry out the detailed modeling of the radio to gamma-ray spectrum of HESS J1731-347 to constrain the magnetic field and energetic particles sources, which we compare with those of the other TeV-bright shell-type SNRs explored before. Methods: Four years of data from Fermi Large Area Telescope (LAT) observations for regions around this remnant are analyzed, leading to no detection correlated with the source discovered in the TeV band. The Markov Chain Monte Carlo method is used to constrain parameters of one-zone models for the overall emission spectrum. Results: Based on the 99.9% upper limits of fluxes in the GeV range, one-zone hadronic models with an energetic proton spectral slope greater than 1.8 can be ruled out, which favors a leptonic origin for the gamma-ray emission, making this remnant a sibling of the brightest TeV SNR RX J1713.7-3946, the Vela Junior SNR RX J0852.0-4622, and RCW 86. The best-fit leptonic model has an electron spectral slope of 1.8 and a magnetic field of about 30 muG, which is at least a factor of 2 higher than those of RX J1713.7-3946 and RX J0852.0-4622, posing a challenge to the distance estimate and/or the energy equipartition between energetic electrons and the magnetic field of this source. A measurement of the shock speed will address this challenge and has implications on the magnetic field evolution and electron acceleration driven by shocks of SNRs.

With the SDO/AIA instrument, continuous and intermittent plasma outflows are observed on the boundaries of an active region along two distinct open coronal loops. %with the speed of 30$\sim$200 km s$^{-1}$. By investigating the temporal sequence magnetograms obtained from HMI/SDO, it is found that small-scale magnetic reconnection probably plays an important role in the generation of the plasma outflows in the coronal loops. It is found that the origin of the plasma outflows coincides with the locations of the small-scale magnetic fields with mixed polarities, which suggests that the plasma outflows along coronal loops probably results from the magnetic reconnection between the small-scale close emerging loops and the large-scale open active region coronal loops.

Magnetic fields dominate most solar activities, there exist direct relations between solar flare and the distributions of magnetic field, and also its corresponding magnetic energy. In this paper, the statistical results about the relationships between the spatial magnetic field and solar flare are given basing on vector magnetic field observed by the Solar Magnetic Field Telescope (SMFT) at Huairou Solar Observing Station (HSOS). The spatial magnetic fields are obtained by extrapolated photosphere vector magnetic field observed by SMFT. There are 23 active regions with flare eruption are chosen as data samples, which were observed from 1997 to 2007. The results are as follows: 1. Magnetic field lines become lower after flare for 16(69%) active regions; 2. The free energy are decreased after flare for 17 (74%) active regions. It can conclude that for most active regions the changes of magnetic field after solar flare are coincident with the previous observations and studies.

Cosmological observables could be used to construct cosmological models, however, a fixed number of observables limited on the light cone is not enough to uniquely determine a certain model. A reconstructed spherically symmetric, inhomogeneous model that share the same angular-diameter-distance-redshift relationship $d_A(z)$ and Hubble parameter $H(z)$ besides $\Lambda\textrm{CDM}$ model (which we call LTB-$\Lambda\textrm{CDM}$ model in this paper), may provide another solution. Cosmic age, which is off the light cone, could be employed to distinguish these two models. We derive the formulae for age calculation with origin conditions. From the data given by 9-year WMAP measurement, we compute the likelihood of the parameters in these two models respectively by using the Distance Prior method and do likelihood analysis by generating Monte Carlo Markov Chain for the purpose of breaking the degeneracy of $\Omega_m$ and $H_0$ (the parameters that we use for calculation). The results yield to be: $t_{\Lambda\textrm{CDM}} =13.76 \pm 0.09 ~\rm Gyr$, $t_{\rm {LTB}-\Lambda\textrm{CDM}} =11.38 \pm 0.15 ~\rm Gyr$, both in $1\sigma$ agreement with the constraint of cosmic age given by metal-deficient stars. The cosmic age method that is set in this paper enables us to distinguish between the inhomogeneous model and $\Lambda\textrm{CDM}$ model.

We have built an electron-phonon coupling model to describe the behavior of the NbxSi1-x transition edge sensor (TES) bolometers, fabricated by electron-beam coevaporation and photolithography techniques on a 2-inch silicon wafer. The resistance versus temperature curves of several sensors with different thickness are measured with different bias currents, ranging from 200 nA to 10 micro A, and the electron-phonon coupling coefficient and the electron-phonon thermal conductance are calculated herein. Our values are quite comparable with those in metallic TES samples of other groups using different measurement methods, while we are using the transition region of our TES sample to calculate the electron-phonon coupling interaction.

We present the various science cases for building Band 1 receivers as part of ALMA's ongoing Development Program. We describe the new frequency range for Band 1 of 35-52 GHz, a range chosen to maximize the receiver suite's scientific impact. We first describe two key science drivers: 1) the evolution of grains in protoplanetary disks and debris disks, and 2) molecular gas in galaxies during the era of re-ionization. Studies of these topics with Band 1 receivers will significantly expand ALMA's Level 1 Science Goals. In addition, we describe a host of other exciting continuum and line science cases that require ALMA's high sensitivity and angular resolution. For example, ALMA Band 1 continuum data will probe the Sunyaev-Zel'dovich Effect in galaxy clusters, Very Small Grains and spinning dust, ionized jets from young stars, spatial and flaring studies of Sgr A*, the acceleration sites of solar flares, pulsar wind nebulae, radio supernovae, and X-ray binaries. Furthermore, ALMA Band 1 line data will probe chemical differentiation in cloud cores, complex carbon chain molecules, extragalactic radio recombination lines, masers, magnetic fields through Zeeman effect measurements, molecular outflows from young stars, the co-evolution of star formation and active galactic nuclei, and the molecular content of galaxies at z ~ 3. ALMA provides similar to better sensitivities than the JVLA over 35-50 GHz, with differences increasing with frequency. ALMA's smaller antennas and shorter baselines, greater number of baselines, and single-dish capabilities, however, give it a significant edge for observing extended emission, making wide-field maps (mosaics), or attaining high image fidelity, as required by the described science cases.

To model a given time series $F(t)$ with fractal Brownian motions (fBms), it is necessary to have appropriate error assessment for related quantities. Usually the fractal dimension $D$ is derived from the Hurst exponent $H$ via the relation $D=2-H$, and the Hurst exponent can be evaluated by analyzing the dependence of the rescaled range $\langle|F(t+\tau)-F(t)|\rangle$ on the time span $\tau$. For fBms, the error of the rescaled range not only depends on data sampling but also varies with $H$ due to the presence of long term memory. This error for a given time series then can not be assessed without knowing the fractal dimension. We carry out extensive numerical simulations to explore the error of rescaled range of fBms and find that for $0<H<0.5$, $|F(t+\tau)-F(t)|$ can be treated as independent for time spans without overlap; for $0.5<H<1$, the long term memory makes $|F(t+\tau)-F(t)|$ correlated and an approximate method is given to evaluate the error of $\langle|F(t+\tau)-F(t)|\rangle$. The error and fractal dimension can then be determined self-consistently in the modeling of a time series with fBms.

With the advent of high-resolution high-sensitivity observations, spiral patterns have been revealed around several asymptotic giant branch (AGB) stars. Such patterns can provide possible evidence for the existence of central binary stars embedded in outflowing circumstellar envelopes. Here, we suggest the viability of explaining the previously observed incomplete ring-like patterns with the spiral-shell structure due to the motion of (unknown) binary components viewed at an inclination with respect to the orbital plane. We describe a method of extracting such spiral-shells from an incomplete ring-like pattern to place constraints on the characteristics of the central binary stars. The use of gas kinematics is essential in facilitating a detailed modeling for the three-dimensional structure of the circumstellar pattern. We show that a hydrodynamic radiative transfer model can reproduce the structure of the HC3N molecular line emission of the extreme carbon star, CIT 6. This method can be applied to other sources in the AGB phase and to the outer ring-like patterns of pre-planetary nebulae for probing the existence of embedded binary stars, which are highly anticipated with future observations using the Atacama Large Millimeter/submillimeter Array.

We use the Submillimeter Array to observe, at 1.4 mm, the blue-lobe of the L1157 outflow at high spatial resolution (~ 3"). We detected SiO, H_2CO, and CH_3OH lines from several molecular clumps that constitute the outflow. All three molecules were detected along the wall of the inner cavity that is supposedly related with the later ejection event. On the other hand, no emission was detected towards positions related to an old ejection episode, likely due to space filtering from the interferometer. The H_2CO and CH_3OH emission is detected only at velocities close to the systemic velocity. The spatial distributions of the H_2CO and CH_3OH are similar. These emission lines trace the U-shaped structure seen in the mid-infrared image. In contrast, the SiO emission is detected in wider velocity range with a peak at ~14 km s/s blue-shifted from the systemic velocity. The SiO emission is brightest at the B1 position, which corresponds to the apex of the U-shaped structure. There are two compact SiO clumps along the faint arc-like feature to the east of the U-shaped structure. At the B1 position, there are two velocity components; one is a compact clump with a size of ~1500 AU seen in the high-velocity and the other is an extended component with lower velocities. The kinematic structure at the B1 position is different from that expected in a single bow shock. It is likely that the high-velocity SiO clump at the B1 position is kinetically independent from the low-velocity gas. The line ratio between SiO (5--4) and SiO (2--1) suggests that the high velocity SiO clumps consist of high density gas of n ~ 10^5 - 10^6 cm^-3, which is comparable to the density of the bullets in the extremely high velocity (EHV) jets. It is likely that the high-velocity SiO clumps in L1157 have the same origin as the EHV bullets.

We discuss 94 Extended Green Objects (EGOs) in the northern Galactic plane cataloged by Cyganowski et al, based on near-infrared narrowband H2 (2.122 {\mu}m and continuum observations from the United Kingdom Infrared Telescope. This data set is three times larger than our previous study, and is unbiased by preselection. As discussed in the previous paper, the morphologies of the 4.5 {\mu}m emission generally resemble those of the near-infrared continuum, but are different from those of the H2 emission. Of our sample, only 28% of EGOs with H2 emission show similar morphologies between 4.5 {\mu}m and H2 emission. These results suggest that the 4.5 {\mu}m emission mainly comes from scattered continuum from the embedded young stellar object (YSO), and partially from H2 emission. About half of EGOs are associated with H2 outflows, if the H2 outflow incompleteness is considered. The H2 outflow detection rate for EGOs with K-band detections (61%) is significantly higher than for those without K-band detections (36%). This difference may be due to the fact that both H2 and K-band emissions are associated with outflows, i.e., H2 emission and K-band continuum are associated with shocks and outflow cavities, respectively. We also compared the correlation between the H2 outflows and Class I 44 GHz methanol masers from literature. The methanol masers can be located upstream or downstream of the H2 outflows and some bright H2 spots or outflows are not associated with methanol masers, suggesting that methanol masers and H2 emission trace different excitation conditions.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The LHAASO (Large High Altitude Air Shower Observatory) experiment is proposed for very high energy gamma ray source survey, in which the WCDA (Water Cherenkov Detector Array) is the one of the major components. In the WCDA, a total of 3600 PMTs are placed under water in four ponds, each with a size of 150 m x 150 m. Precise time and charge measurement is required for the PMT signals, over a large signal amplitude range from single P.E. (photo electron) to 4000 P.E. To fulfill the high requirement of signal measurement in so many front end nodes scattered in a large area, special techniques are developed, such as multiple gain readout, hybrid transmission of clocks, commands, and data, precise clock phase alignment, and new trigger electronics. We present the readout electronics architecture for the WCDA and several prototype modules, which are now under test in the laboratory.

Current observations of supernova remnant (SNR) RX J1713.7-3946 favor the leptonic scenario for the TeV emission, where the radio to X-ray emission is produced via the synchrotron process and the $\gamma$-ray emission is produced via the inverse Comptonization of soft background photons, and the electron distribution can be inferred from the observed $\gamma$-ray spectrum with a spectral inversion method. It is shown that the observed correlation between the X-ray and $\gamma$-ray brightness of SNR RX J1713.7-3946 can be readily explained with the assumption that the energy density of energetic electrons is proportional to that of the magnetic field in such a scenario. A 2D magnetohydrodynamic simulation is then carried out to model the overall emission spectrum. It is found that the total energy of electrons above $\sim 1$ GeV is equal to that of the magnetic field. This is the first piece of observational evidence for energy equipartition between energetic electrons and magnetic field in the downstream of strong collision-less astrophysical shocks of SNRs.

China JinPing underground Laboratory (CJPL) is the deepest underground laboratory presently running in the world. In such a deep underground laboratory, the cosmic ray flux is a very important and necessary parameter for rare event experiments. A plastic scintillator telescope system has been set up to measure the cosmic ray flux. The performance of the telescope system has been studied using the cosmic ray on the ground laboratory near CJPL. Based on the underground experimental data taken from November 2010 to December 2011 in CJPL, which has effective live time of 171 days, the cosmic ray muon flux in CJPL is measured to be (2.0+-0.4)*10^(-10)/(cm^2)/(s). The ultra-low cosmic ray background guarantees CJPL's ideal environment for dark matter experiment.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

Understanding the relationship among different emission components plays an essential role in the study of particle acceleration and energy conversion in solar flares. In flares where gradual and impulsive emission components can be readily identified the impulsive emission has been attributed to non-thermal particles. We carry out detailed analysis of H\alpha\ and X-ray observations of a GOES class B microflare loop on the solar disk. The impulsive hard X-ray emission, however, is found to be consistent with a hot, quasi-thermal origin, and there is little evidence of emission from chromospheric footpoints, which challenges conventional models of flares and reveals a class of microflares associated with dense loops. H\alpha\ observations indicate that the loop lies very low in the solar corona or even in the chromosphere and both emission and absorption materials evolve during the flare. The enhanced H\alpha\ emission may very well originate from the photosphere when the low-lying flare loop heats up the underlying chromosphere and reduces the corresponding H\alpha\ opacity. These observations may be compared with detailed modeling of flare loops with the internal kink instability, where the mode remains confined in space without apparent change in the global field shape, to uncover the underlying physical processes and to probe the structure of solar atmosphere.

Ultra-luminous infrared galaxies (ULIRGs) are the most luminous and intense starburst galaxies in the Universe. Both their star-formation rate (SFR) and gas surface mass density are very high, implying a high supernovae rate and an efficient energy conversion of energetic protons. A small fraction of these supernovae is the so-called hypernovae with a typical kinetic energy ~1e52 erg and a shock velocity >=1e9 cm/s. The strong shocks driven by hypernovae are able to accelerate cosmic ray protons up to 1e17 eV. These energetic protons lose a good fraction of their energy through proton-proton collision when ejected into very dense interstellar medium, and as a result, produce high energy neutrinos (<=5 PeV). Recent deep infrared surveys provide solid constraints on the number density of ULIRGs across a wide redshift range 0<z<2.3, allowing us to derive the flux of diffuse neutrinos from hypernovae. We find that at PeV energies, the diffuse neutrinos contributed by ULIRGs are comparable with the atmosphere neutrinos with the flux of 2e-9GeV cm^-2/s/sr, by assuming the injected cosmic ray power law spectrum with an index of -2.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

A large population of planetary candidates in short-period orbits have been found through transit searches. Radial velocity surveys have also revealed several Jupiter-mass planets with highly eccentric orbits. Measurements of the Rossiter-McLaughlin effect indicate some misaligned planetary systems. This diversity could be induced by post-formation dynamical processes such as planet-planet scattering, the Kozai effect, or secular chaos which brings planets to the vicinity of their host stars. In this work, we propose a novel mechanism to form close-in super-Earths and Neptune-like planets through the tidal disruption of giant planets as a consequence of these dynamical processes. We model the core-envelope structure of giant planets with composite polytropes. Using three-dimensional hydrodynamical simulations of close encounters between planets and their host stars, we find that the presence of a core with a mass more than ten Earth masses can significantly increase the fraction of envelope which remains bound to it. After the encounter, planets with cores are more likely to be retained by their host stars in contrast with previous studies which suggested that coreless planets are often ejected. As a substantial fraction of their gaseous envelopes is preferentially lost while the dense incompressible cores retain most of their original mass, the resulting metallicity of the surviving planets is increased. Our results suggest that some gas giant planets can be effectively transformed into either super-Earths or Neptune-like planets after multiple close stellar passages. Finally, we analyze the orbits and structure of known planets and Kepler candidates and find that our model is capable producing some of the shortest-period objects.

Shocks of supernova remnants (SNRs) are important (and perhaps the dominant) agents for production of the Galactic cosmic rays. Recent $\gamma$-ray observations of several SNRs have made this case more compelling. However, these broadband high-energy measurements also reveal a variety of spectral shape demanding more comprehensive modeling of emissions from SNRs. According to the locally observed fluxes of cosmic ray protons and electrons, the electron-to-proton number ratio is known to be about 1%. Assuming such a ratio is universal for all SNRs and identical spectral shape for all kinds of accelerated particles, we propose a unified model that ascribes the distinct $\gamma$-ray spectra of different SNRs to variations of the medium density and the spectral difference between cosmic ray electrons and protons observed at Earth to transport effects. For low density environments, the $\gamma$-ray emission is inverse-Compton dominated. For high density environments like systems of high-energy particles interacting with molecular clouds, the $\gamma$-ray emission is $\pi^0$-decay dominated. The model predicts a hadronic origin of $\gamma$-ray emission from very old remnants interacting mostly with molecular clouds and a leptonic origin for intermediate age remnants whose shocks propagate in a low density environment created by their progenitors via e.g., strong stellar winds. These results can be regarded as evidence in support of the SNR-origin of the Galactic cosmic rays.

We present a set of compelling science cases for the ALMA Band 1 receiver suite. For these cases, we assume in tandem the updated nominal Band 1 frequency range of 35-50 GHz with a likely extension up to 52 GHz; together these frequencies optimize the Band 1 science return. The scope of the science cases ranges from nearby stars to the re-ionization edge of the Universe. Two cases provide additional leverage on the present ALMA Level One Science Goals and are seen as particularly powerful motivations for building the Band 1 Receiver suite: (1) detailing the evolution of grains in protoplanetary disks, as a complement to the gas kinematics, requires continuum observations out to ~35 GHz (~9mm); and (2) detecting CO 3-2 line emission from galaxies like the Milky Way during the epoch of re-ionization, i.e., 6 < z < 10, also requires Band 1 receiver coverage. The range of Band 1 science is wide, however, and includes studies of very small dust grains in the ISM, pulsar wind nebulae, radio supernovae, X-ray binaries, the Galactic Center (i.e., Sgr A*), dense cloud cores, complex carbon-chain molecules, masers, magnetic fields in the dense ISM, jets and outflows from young stars, distant galaxies, and galaxy clusters (i.e., the Sunyaev-Zel'dovich Effect). A comparison of ALMA and the Jansky VLA (JVLA) at the same frequencies of Band 1 finds similar sensitivity performance at 40-50 GHz, with a slight edge for ALMA at higher frequencies (e.g., within a factor of 2 for continuum observations). With its larger number of instantaneous baselines, however, ALMA Band 1data will have greater fidelity than those from the JVLA at similar frequencies.

We report limits on annual modulation of the low-energy event rate from the Cryogenic Dark Matter Search (CDMS II) experiment at the Soudan Underground Laboratory. Such a modulation could be produced by interactions from Weakly Interacting Massive Particles (WIMPs) with masses ~10 GeV/c^2. We find no evidence for annual modulation in the event rate of veto-anticoincident single-detector interactions consistent with nuclear recoils, and constrain the magnitude of any modulation to <0.06 event [keVnr kg day]^-1 in the 5-11.9 keVnr energy range at the 99% confidence level. These results disfavor an explanation for the reported modulation in the 1.2-3.2 keVee energy range in CoGeNT in terms of nuclear recoils resulting from elastic scattering of WIMPs at >98% confidence. For events consistent with electron recoils, no significant modulation is observed for either single- or multiple-detector interactions in the 3.0-7.4 keVee range.

It is an indisputable fact that solar magnetic fields are force-free in the corona, where force free fields means that current and magnetic fields are parallel and there is no Lorentz force in the fields. While the force-free extent of photospheric magnetic fields remains open. In this paper, the statistical results about it is given. The vector magnetograms (namely, $B_{x}$, $B_{y}$ and $B_{z}$ in heliocentric coordinates) are employed, which are deduced and calibrated from Stokes spectra, observed by Solar Magnetic Field Telescope (SMFT) at Huairou Solar Observing Station (HSOS) are used. We study and calibrated 925 magnetograms calibrated by two sets of calibration coefficients, that indicate the relations between magnetic fields and the strength of Stokes spectrum and can be calculated either theoretically or empirically. The statistical results show that the majority of active region magnetic fields are not consistent with the force-free model.

We present observational data for two main components (S255IR and S255N) of the S255 high mass star forming region in continuum and molecular lines obtained at 1.3 mm and 1.1 mm with the SMA, at 1.3 cm with the VLA and at 23 and 50 cm with the GMRT. The angular resolution was from ~ 2" to ~ 5" for all instruments. With the SMA we detected a total of about 50 spectral lines of 20 different molecules (including isotopologues). About half of the lines and half of the species (in particular N2H+, SiO, C34S, DCN, DNC, DCO+, HC3N, H2CO, H2CS, SO2) have not been previously reported in S255IR and partly in S255N at high angular resolution. Our data reveal several new clumps in the S255IR and S255N areas by their millimeter wave continuum emission. Masses of these clumps are estimated at a few solar masses. The line widths greatly exceed expected thermal widths. These clumps have practically no association with NIR or radio continuum sources, implying a very early stage of evolution. At the same time, our SiO data indicate the presence of high-velocity outflows related to some of these clumps. In some cases, strong molecular emission at velocities of the quiescent gas has no detectable counterpart in the continuum. We discuss the main features of the distribution of NH3, N2H+, and deuterated molecules. We estimate properties of decimeter wave radio continuum sources and their relationship with the molecular material.

We present the first interferometric polarization map of the W3(OH) massive star-forming region observed with the Submillimeter Array (SMA) at 878 mum with an angular resolution of 1.5 (about 3 \times 10 AU). Polarization is detected in the W3(H2O) hot core, an extended emission structure in the north-west of W3(H2O), and part of the W3(OH) ultracompact HII region. The W3(H2O) hot core is known to be associated with a synchrotron jet along the east-west direction. In this core, the inferred magnetic field orientation is well aligned with the synchrotron jet and close to the plane of sky. Using the Chandrasekhar-Fermi method with the observed dispersion in polarization angle, we estimate a plane-of-sky magnetic field strength of 17.0 mG. Combined with water maser Zeeman measurements, the total magnetic field strength is estimated to be 17.1 mG, comparable to the field strength estimated from the synchrotron model. The magnetic field energy dominates over turbulence in this core. In addition, the depolarization effect is discerned in both SMA and JCMT measurements. Despite the great difference in angular resolutions and map extents, the polarization percentage shows a similar power-law dependence with the beam averaged column density. We suggest that the column density may be an important factor to consider when interpreting the depolarization effect.

Stochastic acceleration of charged particles due to their interactions with plasma waves may be responsible for producing superthermal particles in a variety of astrophysical systems. This process can be described as a diffusion process in the energy space with the Fokker-Planck equation. In this paper, a time-dependent numerical code is used to solve the reduced Fokker-Planck equation involving only time and energy variables with general forms of the diffusion coefficients. We also propose a self-similar model for particle acceleration in Sedov explosions and use the TeV SNR RX J1713.7-3946 as an example to demonstrate the model characteristics. Markov Chain Monte Carlo method is utilized to constrain model parameters with observations.