Dense molecular filaments are central to the star formation process, but the detailed manner in which they fragment into prestellar cores is not yet well understood. Here, we investigate the fragmentation properties and dynamical state of several star-forming filaments in the X-shape Nebula region of the California MC, in an effort to shed some light on this issue. We used multi-wavelength far-infrared images from {\it Herschel} and the \textsl{getsources} and \textsl{getfilaments} extraction methods to identify dense cores and filaments and derive their basic properties. We also used a map of $\rm ^{13}CO (2-1)$ emission from SMT 10m submillimeter telescope to constrain the dynamical state of the filaments. We identified 10 filaments, as well as 57 dense cores. Two star-forming filaments (\# 8 and \# 10) stand out in that they harbor quasi-periodic chains of dense cores with a typical projected core spacing of $\sim$0.15 pc. These two filaments have thermally supercritical line masses and are not static. Filament~8 exhibits a prominent transverse velocity gradient, suggesting that it is accreting gas from the parent cloud gas reservoir. In both cases, the observed (projected) core spacing is similar to the filament width and significantly shorter than the canonical separation of $\sim \,$4 times the filament width predicted by classical cylinder fragmentation theory. We suggest that continuous accretion of gas onto the two star-forming filaments, as well as geometrical bending of the filaments, may account for the observed core spacing. Our findings suggest that the characteristic fragmentation lengthscale of molecular filaments is quite sensitive to external perturbations from the parent cloud, such as gravitational accretion of ambient material.

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 used two XMM-Newton and six Neutron Star Interior Composition Explorer (NICER) observations to investigate the fractional rms amplitude of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary 4U 1636-53. We studied, for the first time, the fractional rms amplitude of the mHz QPOs vs. energy in 4U 1636-53 down to 0.2 keV. We find that, as the energy increases from 0.2 keV to 3 keV, the rms amplitude of the mHz QPOs increases, different from the decreasing trend that has been previously observed above 3 keV. This finding has not yet been predicted by any current theoretical model, however, it provides an important observational feature to speculate whether a newly discovered mHz oscillation originates from the marginally stable nuclear burning process on the neutron star surface.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

We performed a multi-wavelength observation toward LkHa 101 embedded cluster and its adjacent 85arcmin*60arcmin region. The LkHa 101 embedded cluster is the first and only one significant cluster in California molecular cloud (CMC). These observations have revealed that the LkHa 101 embedded cluster is just located at the projected intersectional region of two filaments. One filament is the highest-density section of the CMC, the other is a new identified filament with a low-density gas emission. Toward the projected intersection, we find the bridging features connecting the two filaments in velocity, and identify a V-shape gas structure. These agree with the scenario that the two filaments are colliding with each other. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we measured that the RRL velocity of the LkHa 101 H II region is 0.5 km/s, which is related to the velocity component of the CMC filament. Moreover, there are some YSOs distributed outside the intersectional region. We suggest that the cloud-cloud collision together with the fragmentation of the main filament may play an important role in the YSOs formation of the cluster.

Fast radio bursts (FRBs) are extragalactic, bright pulses of emission at radio frequency with milliseconds duration. Observationally, FRBs can be divided into two classes, repeating FRBs and non-repeating FRBs. At present, twenty repeating FRBs have been discovered with unknown physical origins. Localization of the first repeating FRB 121102 and discovery of an associated persistent radio source support that FRBs are powered by young millisecond magnetars, which could be formed by core-collapses of massive stars or binary neutron stars mergers. These two formation channels can be distinguished by gravitational waves generated by binary neutron stars mergers. We first calculate the lower limit of the local formation rate of repeating FRBs observed by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) . Then we show that only a small fraction ($6\%$) of repeating FRBs is produced by young magnetars from binary neutron star mergers, basing on the gravitational wave detections by the third observing run (O3) of Advanced LIGO/Virgo gravitational-wave detectors. Therefore, we believe that repeating FRBs are more likely produced by the magnetars newborn from the core-collapses of massive stars rather than the magnetars from the binary neutron stars mergers.

Understanding large-angular-scale galactic foregrounds is crucial for future CMB experiments aiming to detect $B$-mode polarization from primordial gravitational waves. Traditionally, the dust component has been separated using its different frequency dependence. However, using non-CMB observations has potential to increase fidelity and decrease the reconstruction noise. In this exploratory paper we investigate the capability of galactic 21 cm observations to predict the dust foreground in intensity. We train a neural network to predict the dust foreground as measured by the Planck Satellite from the full velocity data-cube of galactic 21 cm emission as measured by the HI4PI survey. We demonstrate that information in the velocity structure clearly improves the predictive power over both a simple integrated emission model and a simple linear model. The improvement is significant at arc-minute scales but more modest at degree scales. This proof of principle on temperature data indicates that it might also be possible to improve foreground polarization templates from the same input data.

Four significant events of rapid $^{14}$C increase have taken place within the past several thousand years. The physical origin of these rapid increases is still a mystery but must be associated with extremely energetic cosmic processes. Pulsars are highly magnetized neutron stars that emit a beam of electromagnetic radiations. Any sudden release of the energy stored in the magnetic multipole field will trigger outbursts similar to the giant flares of magnetars. Here we show that the relativistic outflow from the outbursts of a nearby pulsar interacting with the interstellar medium generates a shock, which accelerates electrons to trillions of electron volts. The high-energy photons from synchrotron emission of the shock interact with Earth's atmosphere, producing the cosmogenic nuclide $^{14}$C, which can cause the rapid $^{14}$C increases discovered in tree rings. These same relativistic electrons can account for a significant fraction of the cosmic electron spectrum in the trillion electron volts energy range, as observed by space-borne satellites. Since these outburst events can significantly affect our environment, monitoring nearby pulsars for such outbursts may be important in the future.

In this paper, we present statistics of soft gamma repeater (SGR) bursts from SGR J1550-5418, SGR 1806-20 and SGR 1900+14 by adding new bursts from K{\i}rm{\i}z{\i}bayrak et al. (2017) detected with the Rossi X-ray Timing Explorer (RXTE). We find that the fluence distributions of magnetar bursts are well described by power-law functions with indices 1.84, 1.68, and 1.65 for SGR J1550-5418, SGR 1806-20 and SGR 1900+14, respectively. The duration distributions of magnetar bursts also show power-law forms. Meanwhile, the waiting time distribution can be described by a non-stationary Poisson process with an exponentially growing occurrence rate. These distributive features indicate that magnetar bursts can be regarded as a self-organizing critical process. We also compare these distributions with the repeating FRB 121102. The statistical properties of repeating FRB 121102 are similar with magentar bursts, combing with the large required magnetic filed ($B\geq 10^{14}$G) of neutron star for FRB 121102, which indicates that the central engine of FRB 121102 may be a magnetar.

Cosmic voids, the underdense regions in the universe, are particularly sensitive to diffuse density components such as cosmic neutrinos. This sensitivity is enhanced by the match between void sizes and the free-streaming scale of massive neutrinos. Using the massive neutrino simulations \texttt{MassiveNuS}, we investigate the effect of neutrino mass on dark matter halos as a function of environment. We find that the halo mass function depends strongly on neutrino mass and that this dependence is more pronounced in voids than in high-density environments. An observational program that measured the characteristic mass of the most massive halos in voids should be able to place novel constraints on the sum of the masses of neutrinos $\sum m_\nu$. The neutrino mass effect in the simulations is quite strong: In a 512$^3$ $h^{-3}$ Mpc$^3$ survey, the mean mass of the 1000 most massive halos in the void interiors is $(4.82 \pm 0.11) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.6$ eV and $(8.21 \pm 0.13) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.1$ eV. Subaru (SuMIRe), Euclid and WFIRST will have both spectroscopic and weak lensing surveys. Covering volumes at least 50 times larger than our simulations, they should be sensitive probes of neutrino mass through void substructure.

We investigate the frequency and amplitude of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary (NS LMXB) 4U 1636-53 using Rossi X-ray Timing Explorer observations. We find that no mHz QPOs appear when the source is in the hard spectral state. We also find that there is no significant correlation between the frequency and the fractional rms amplitude of the mHz QPOs. Notwithstanding, for the first time, we find that the absolute RMS amplitude of the mHz QPOs is insensitive to the parameter Sa, which measures the position of the source in the colour-colour diagram and is usually assumed to be an increasing function of mass accretion rate. This finding indicates that the transition from marginally stable burning to stable burning or unstable burning could happen very rapidly since, before the transition, the mHz QPOs do not gradually decay as the rate further changes.

We present for the neutron-star low-mass X-ray binary 4U 1636$-$53, and for the first time for any source of kilohertz quasi-periodic oscillations (kHz QPOs), the two-dimensional behaviour of the fractional rms amplitude of the kHz QPOs in the parameter space defined by QPO frequency and photon energy. We find that the rms amplitude of the lower kHz QPO increases with energy up to $\sim12$ keV and then decreases at higher energies, while the rms amplitude of the upper kHz QPO either continues increasing or levels off at high energies. The rms amplitude of the lower kHz QPO increases and then decreases with frequency, peaking at $\sim 760$ Hz, while the amplitude of the upper kHz QPO decreases with frequency, with a local maximum at around $\sim 770$ Hz, and is consistent with becoming zero at the same QPO frequency, $\sim1400$ Hz, in all energy bands, thus constraining the neutron-star mass at $M_{NS} \leq 1.6 M_{\odot}$, under the assumption that this QPO reflects the Keplerian frequency at the inner edge of the accretion disc. We show that the slope of the rms energy spectrum is connected to the changing properties of the kHz QPOs in different energy bands as its frequencies change. Finally, we discuss a possible mechanism responsible for the radiative properties of the kHz QPOs and, based on a model in which the QPO arises from oscillations in a Comptonising cloud of hot electrons, we show that the properties of the kHz QPOs can constrain the thermodynamic properties of the inner accretion flow.

Fast radio bursts (FRBs) are milliseconds radio transients with large dispersion measures (DMs). An outstanding question is the relation between repeating FRBs and those with a single burst. In this paper, we study the energy distribution of the repeating FRB 121102. For a power-law distribution of energy $dN/dE\propto E^{-\alpha_E}$, we show that the value of $\alpha_E$ is in a narrow range $1.6-1.8$ for bursts observed by different telescopes at different frequencies, which indicates a universal energy distribution for FRB 121102. Interestingly, similar power-law index of energy distribution for non-repeating FRBs observed by Parkes and ASKAP is also found. However, if low-energy bursts below completeness threshold of Arecibo are discarded for FRB 121102, the slope could be up to 2.2. Implications of such a universal energy distribution are discussed.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

In this paper, we investigate the energy function, formation rate and environment of fast radio bursts (FRBs) using Parkes sample and Australian Square Kilometer Array Pathfinder (ASKAP) sample. For the first time, the metallicity effect on the formation rate is considered. If FRBs are produced by the mergers of compact binaries, the formation rate of FRBs should have a time delay relative to cosmic star formation rate (CSFR). We get the time delay is about 3-5 Gyr and the index of differential energy function $\gamma$ ($dN/dE\propto E^{-\gamma}$) is between 1.6 and 2.0 from redshift cumulative distribution. The value of $\gamma$ is similar to that of FRB 121102, which indicates single bursts may share the same physical mechanism with the repeaters. In another case, if the formation rate of FRB is proportional to the SFR without time delay, the index $\gamma$ is about 2.3. In both cases, we find that FRBs may prefer to occur in low-metallicity environment with $ 12 +\log(\rm{O/H}) \simeq 8.40$, which is similar to those of long gamma-ray bursts (GRBs) and hydrogen-poor superluminous supernovae (SLSNe-I).

We aim to investigate the impact of the ionized radiation from the M16 HII region on the surrounding molecular cloud and on its hosted star formation. To present comprehensive multi-wavelength observations towards the M16 HII region, we used new CO data and existing infrared, optical, and submillimeter data. The 12CO J=1-0, 13CO J=1-0, and C18O J=1-0 data were obtained with the Purple Mountain Observatory (PMO) 13.7m radio telescope. To trace massive clumps and extract young stellar objects (YSOs) associated with the M16 HII region, we used the ATLASGAL and GLIMPSE I catalogs, respectively. From CO data, we discern a large-scale filament with three velocity components. Because these three components overlap with each other in both velocity and space, the filament may be made of three layers. The M16 ionized gas interacts with the large-scale filament and has reshaped its structure. In the large-scale filament, we find 51 compact cores from the ATLASGAL catalog, 20 of them being quiescent. The mean excitation temperature of these cores is 22.5 K, while this is 22.2 K for the quiescent cores. This high temperature observed for the quiescent cores suggests that the cores may be heated by M16 and do not experience internal heating from sources in the cores. Through the relationship between the mass and radius of these cores, we obtain that 45% of all the cores are massive enough to potentially form massive stars. Compared with the thermal motion, the turbulence created by the nonthermal motion is responsible for the core formation. For the pillars observed towards M16, the H II region may give rise to the strong turbulence.

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.

Fast radio bursts (FRBs) are bright milliseconds radio transients with large dispersion measures (DMs), and can be used as potential astrophysical and cosmological tools. Many models for FRB progenitors have been proposed. However, none of them can explain all the observational properties, meaning that FRBs remain one of the most intriguing mysteries in astronomy. The statistical properties of radio bursts can unveil the underlying physics. Here, we report statistical results of the repeating FRB 121102, and show that FRB 121102 and solar type III radio bursts share four statistical properties: power-law frequency distributions for energies, fluxes, durations and waiting times. All of the distributions can be explained by avalanche models of self-organized criticality (SOC) systems. It is well known that solar type III radio burst arises from the nonlinear conversion of Langmuir waves generated by two-stream instability of electron beams, which are accelerated by bursty magnetic reconnections. The similarities support that repeating FRBs are coherent emissions from magnetic-reconnection-driven beams in the magnetospheres of magnetars.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present optical, UV and X-ray monitoring of the short orbital period black hole X-ray binary candidate Swift J1753.5-0127, focusing on the final stages of its 12$-$year long outburst that started in 2005. From September 2016 onward, the source started to fade and within three months, the optical flux almost reached the quiescent level. Soon after that, using a new proposed rebrightening classification method we recorded a mini-outburst and a reflare in the optical light curves, peaking in February (V$\rm\sim$17.0) and May (V$\rm\sim$17.9) 2017, respectively. Remarkably, the mini-outburst has a peak flux consistent with the extrapolation of the slow decay before the fading phase preceding it. The following reflare was fainter and shorter. We found from optical colors that the temperature of the outer disk was $\sim 11$,000 K when the source started to fade rapidly. According to the disk instability model, this is close to the critical temperature when a cooling wave is expected to form in the disk, shutting down the outburst. The optical color could be a useful tool to predict decay rates in some X-ray transients. We notice that all X-ray binaries that show mini-outbursts following a main outburst are short orbital period systems ($<$ 7 h). In analogy with another class of short period binaries showing similar mini-outbursts, the cataclysmic variables of the RZ LMi type, we suggest mini-outbursts could occur if there is a hot inner disk at the end of the outburst decay.

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.

Disc instability models predict that for X-ray binaries in quiescence, there should be a brightening of the optical flux prior to an X-ray outburst. Tracking the X-ray variations of X-ray binaries in quiescence is generally not possible, so optical monitoring provides the best means to measure the mass accretion rate variability between outbursts, and to identify the beginning stages of new outbursts. With our regular Faulkes Telescope/Las Cumbres Observatory (LCO) monitoring we are routinely detecting the optical rise of new X-ray binary outbursts before they are detected by X-ray all-sky monitors. We present examples of detections of an optical rise in X-ray binaries prior to X-ray detection. We also present initial optical monitoring of the new black hole transient MAXI J1820+070 (ASASSN-18ey) with the Faulkes, LCO telescopes and Al Sadeem Observatory in Abu Dhabi, UAE. Finally, we introduce our new real-time data analysis pipeline, the 'X-ray Binary New Early Warning System (XB-NEWS)' which aims to detect and announce new X-ray binary outbursts within a day of first optical detection. This will allow us to trigger X-ray and multi-wavelength campaigns during the very early stages of outbursts, to constrain the outburst triggering mechanism.

The Near-infrared Sky Brightness Monitor (NIRBM) aims to measure the middle infrared sky background in Antarctica. The NIRBM mainly consists of an InGaAs detector, a chopper, a reflector, a cooler and a black body. The reflector can rotate to scan the sky with a field of view ranging from 0{\deg} to 180{\deg}. Electromechanical control and weak signal readout functions are accomplished by the same circuit, whose core chip is a STM32F407VG microcontroller. Considering the environment is harsh for humans in Antarctica, a multi-level remote control software system is designed and implemented. A set of EPICS IOCs are developed to control each hardware module independently via serial port communication with the STM32 microcontroller. The tornado web framework and PyEpics are introduced as a combination where PyEpics is used to monitor or change the EPICS Process Variables, functioning as a client for the EPICS framework. Tornado is responsible for the specific operation process of inter-device collaboration, and expose a set of interfaces to users to make calls. Considering the high delay and low bandwidth of the network environment, the tornado back-end is designed as a master-and-agent architecture to improve domestic user experience. The master node is deployed in Antarctic while multiple agent nodes can be deployed domestic. The master and agent nodes communicate with each other through the WebSocket protocol to exchange latest information so that bandwidth is saved. The GUI is implemented in the form of single-page application based on the Vue framework which communicates with tornado through WebSocket and AJAX requests. The web page integrates device control, data curve drawing, alarm display, auto observation and other functions together.

We describe an improved nonequilibrium ionization (NEI) method that we have developed as an optional module for the FLASH magnetohydrodynamic simulation code. The method employs an eigenvalue approach rather than the earlier iterative ordinary differential equation approach to solve the stiff differential equations involved in NEI calculations. The new code also allows the atomic data to be easily updated from the AtomDB database. We compare both the updated atomic data and the methods separately. The new atomic data are shown to make a significant difference in some circumstances, although the general trends remain the same. Additionally, the new method also allows simultaneous calculation of the nonequilibrium radiative cooling, which is not included in the original method. The eigenvalue method improves the calculation efficiency overall with no loss of accuracy. We explore some common ways to present the NEI state with a sample simulation and find that using the average ionic charge difference from the equilibrium tends to be the clearest method.

The mixed morphology class of supernova remnants (MMSNRs) comprises a substantial fraction of observed remnants and yet there is as yet no consensus on their origin. A clue to their nature is the presence of regions that show X-ray evidence of recombining plasmas. Recent calculations of remnant evolution in a cloudy interstellar medium (ISM) that included thermal conduction but not non-equilibrium ionization (NEI) showed promise in explaining observed surface brightness distributions but could not determine if recombining plasmas were present. In this paper, we present numerical hydrodynamical models of MMSNRs in 2D and 3D including explicit calculation of NEI effects. Both the spatial ionization distribution and temperature-density diagrams show that recombination occurs inside the simulated MMSNR, and that both adiabatic expansion and thermal conduction cause recombination, albeit in different regions. Features created by the adiabatic expansion stand out in the spatial and temperature-density diagrams, but thermal conduction also plays a role. Thus thermal conduction and adiabatic expansion both contribute significantly to the cooling of high-temperature gas. Realistic observational data are simulated with both spatial and spectral input from various regions. We also discuss the possibility of analyzing the sources of recombination and dominant hydrodynamical processes in observations using temperature-density diagrams and spatial maps.

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

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We studied five XMM-Newton observations of the neutron-star binary 4U 1728$-$34 covering the hard, intermediate and soft spectral states. By jointly fitting the spectra with several reflection models, we obtained an inclination angle of 25$-$53$\deg$ and an iron abundance up to 10 times the solar. From the fits with reflection models, we found that the fluxes of the reflection and the Comptonised components vary inconsistently; since the latter is assumed to be the illuminating source, this result possibly indicates the contribution of the neutron star surface/boundary layer to the disc reflection. As the source evolved from the relatively soft to the intermediate state, the disc inner radius decreased, opposite to the prediction of the standard accretion disc model. We also explore the possible reasons why the supersolar iron abundance is required by the data and found that this high value is probably caused by the absence of the hard photons in the XMM-Newton data.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

More than one hundred galaxy-scale strong gravitational lens systems have been found by searching for the emission lines coming from galaxies with redshifts higher than the lens galaxies. Based on this spectroscopic-selection method, we introduce the deep Residual Networks (ResNet, a kind of deep Convolutional Neural Networks) to search for the galaxy-Ly$\alpha$ emitter (LAE) lens candidates by recognizing the Ly$\alpha$ emission lines coming from high redshift galaxies ($2 < z < 3$) in the spectra of early-type galaxies (ETGs) at middle redshift ($z\sim 0.5$). The spectra of the ETGs come from the Data Release 12 (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey \uppercase\expandafter{\romannumeral3} (SDSS-\uppercase\expandafter{\romannumeral3}). In this paper, we first build a 28 layers ResNet model, and then artificially synthesize 150,000 training spectra, including 140,000 spectra without Ly$\alpha$ lines and 10,000 ones with Ly$\alpha$ lines, to train the networks. After 20 training epochs, we obtain a near-perfect test accuracy at 0.9954. The corresponding loss is 0.0028 and the completeness is 93.6\%. We finally apply our ResNet model to our predictive data with 174 known lens candidates. We obtain 1232 hits including 161 of the 174 known candidates (92.5\% discovery rate). Apart from the hits found in other works, our ResNet model also find 536 new hits. We then perform several subsequent selections on these 536 hits and present 5 most believable lens candidates.

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.

MARST (Multi-channel Antarctic Solar Telescope) aims to observe the Sun in multiple wavebands in Antarctica and will be China's first solar telescope in Antarctica. The telescope has two tubes, corresponding to Photosphere observation which uses 11 filters and Chromosphere observation in H{\alpha} waveband. The two tubes need to be coordinated to observe at the same time. The telescope will observe the Sun for a long time, so a self-guiding module is needed to improve sun tracking. Besides, performing solar specific flat-field exposure is necessary for analyzing. EPICS is introduced to control each hardware and an autonomous observation system based on RTS2 is designed under such demands. EPICS application modules are implemented for each device: telescope mount & focuser, filter wheel, H{\alpha} filter, dome with webcams, Andor CCD and PI CCD. We also integrate EPICS modules into RTS2 framework with an XML format configuration. To control these applications autonomously we have developed a RTS2 executor module where two plan classes are instantiated to control two sets of filters and CCDs, and to ensure only one could control the mount at the same time. Different types of observation plans are designed to describe different series of processes with different priorities. To improve sun tracking, we calculate the centroid of each image to get the offset, then apply the correction to the telescope during observation process. For frontend users, a GUI based on PyQt5 and QML is implemented and connects to rts2-httpd and rts2-proxy modules so that users can control devices, check images and get logs.

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.

Diffuse soft X-ray line emission is commonly used to trace the thermal and chemical properties of the hot interstellar medium, as well as its content, in nearby galaxies. Although resonant line scattering complicates the interpretation of the emission, it also offers an opportunity to measure the kinematics of the medium. We have implemented a direct Monte Carlo simulation scheme that enables us to account for resonant scattering effect in the medium, in principle, with arbitrary spatial, thermal, chemical, and kinematic distributions. Here we apply this scheme via dimensionless calculation to an isothermal, chemically uniform, and spherically symmetric medium with a radial density distribution characterized by a $\beta$-model. This application simultaneously account for both optical depth-dependent spatial distortion and intensity change of the resonant line emission due to the scattering, consistent with previous calculations. We further apply the modeling scheme to the OVII and OVIII emission line complex observed in the XMM-Newton RGS spectrum of the M31 bulge. This modeling, though with various limitations due to its simplicity, shows that the resonant scattering could indeed account for much of the spatial distortion of the emission, as well as the relative strengths of the lines, especially the large forbidden to resonant line ratio of the OVII He$\alpha$ triplet. We estimate the isotropic turbulence Mach number of the medium in M31 as $\sim0.17$ for the first time and the line-emitting gas temperature as $\sim2.3\times10^6$ K. We conclude that the resonant scattering may in general play an important role in shaping the soft X-ray spectra of diffuse hot gas in normal galaxies.

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.

Automated spectral classification is an active research area in astronomy at the age of data explosion. While new generation of sky survey telescopes (e.g. LAMOST and SDSS) produce huge amount of spectra, automated spectral classification is highly required to replace the current model fitting approach with human intervention. Galaxies, and especially active galactic nucleus (AGNs), are important targets of sky survey programs. Efficient and automated methods for galaxy spectra classification is the basis of systematic study on physical properties and evolution of galaxies. To address the problem, in this paper we carry out an experiment on Alibaba Cloud AI plaform (PAI) to explore automated galaxy spectral classification using machine learning approach. Supervised machine learning algorithms (Logistic Regression, Random Forest and Linear SVM) were performed on a dataset consist of ~ 10000 galaxy spectra of SDSS DR14, and the classification results of which are compared and discussed. These galaxy spectra each has a subclass tag (i.e. AGNs, Starburst, Starforming, and etc.) that we use as training labels.

In this paper, we study the luminosity function and formation rate of short gamma-ray bursts (sGRBs). Firstly, we derive the $E_p-L_p$ correlation using 16 sGRBs with redshift measurements and determine the pseudo redshifts of 284 Fermi sGRBs. Then, we use the Lynden-Bell c$^-$ method to study the luminosity function and formation rate of sGRBs without any assumptions. A strong evolution of luminosity $L(z)\propto (1+z)^{4.47}$ is found. After removing this evolution, the luminosity function is $ \Psi (L) \propto L_0 ^ {- 0.29 \pm 0.01} $ for dim sGRBs and $ \psi (L) \propto L_0 ^ {- 1.07 \pm 0.01} $ for bright sGRBs, with the break point $8.26 \times 10^{50} $ erg s$^{-1}$. We also find that the formation rate decreases rapidly at $z<1.0$, which is different with previous works. The local formation rate of sGRBs is 7.53 events Gpc$^{-3}$ yr$^{-1}$. Considering the beaming effect, the local formation rate of sGRBs including off-axis sGRBs is $ 203.31^{+1152.09}_{-135.54} $ events Gpc$^{-3}$ yr$^{-1}$. We also estimate that the event rate of sGRBs detected by the advanced LIGO and Virgo is $0.85^{+4.82}_{-0.56} $ events yr$^{-1}$ for NS-NS binary.

The origin of the high-frequency radio emission detected from several magnetars is poorly understood. In this paper, we report the ~40 GHz properties of SGR J1745-29 as measured using Jansky Very Large Array (JVLA) and Robert C. Byrd Green Bank Telescope (GBT) observations between 2013 October 26 and 2014 May 31. Our analysis of a Q-band (45 GHz) GBT observation on 2014 April 10 resulted in the earliest detection of pulsed radio emission at high frequencies (>20 GHz); we found that the average pulse has a singly peaked profile with width ~75 ms (~2% of the 3.764 s pulse period) and an average pulsed flux density of ~100 mJy. We also detected very bright, short (<10 ms) single pulses during ~70% of this neutron star's rotations, and the peak flux densities of these bright pulses follow the same log-normal distribution as measured at 8.5 GHz. Additionally, our analysis of contemporaneous JVLA observations suggest that its 41/44 GHz flux density varied between ~1-4 mJy during this period, with a ~2x change observed on ~20 minute timescales during a JVLA observation on 2014 May 10. Such a drastic change over short time-scales is inconsistent with the radio emission resulting from a shock powered by the magnetar's supersonic motion through the surrounding medium, and instead is dominated by pulsed emission generated in its magnetosphere.

We performed a multiwavelength study towards HII region Sh2-104. New maps of 12CO J=1-0 and 13CO J=1-0 were obtained from the Purple Mountain Observatory (PMO) 13.7 m radio telescope. Sh2-104 displays a double-ring structure. The outer ring with a radius of 4.4 pc is dominated by 12 um, 500 um, 12CO J=1-0, and 13CO J=1-0 emission, while the inner ring with a radius of 2.9 pc is dominated by 22 um and 21 cm emission. We did not detect CO emission inside the outer ring. The north-east portion of the outer ring is blueshifted, while the south-west portion is redshifted. The present observations have provided evidence that the collected outer ring around Sh2-104 is a two-dimensional structure. From the column density map constructed by the Hi-GAL survey data, we extract 21 clumps. About 90\% of all the clumps will form low-mass stars. A power-law fit to the clumps yields M=281Msun(r/pc)^1.31. The selected YSOs are associated with the collected material on the edge of Sh2-104. The derived dynamical age of Sh2-104 is 1.6*10^6 yr. Compared the Sh2-104 dynamical age with the YSOs timescale and the fragmentation time of the molecular ring, we further confirm that collect-and-collapse process operates in this region, indicating a positive feedback from a massive star for surrounding gas.

The candidate black hole X-ray binary Swift J1753.5-0127 faded to quiescence in 2016 November, after a prolonged outburst that was discovered in 2005. Nearly three months later the system displayed renewed activity that lasted through 2017 July. Here, we present radio and X-ray monitoring over ~3 months of the renewed activity to study the coupling between the jet and the inner regions of the disk/jet system. Our observations cover low X-ray luminosities that have not historically been well-sampled (Lx~2e33 - 1e36 erg/s; 1-10 keV), including time periods when the system was both brightening and fading. At these low luminosities Swift J1753.5-0127 occupies a parameter space in the radio/X-ray luminosity plane that is comparable to "canonical" systems (e.g., GX 339-4), regardless of whether the system was brightening or fading, even though during its >11-year outburst Swift J1753.5-0127 emitted less radio emission from its jet than expected. We discuss implications for the existence of a single radio/X-ray luminosity correlation for black hole X-ray binaries at the lowest luminosities (Lx < 1e35 erg/s), and we compare to supermassive black holes. Our campaign includes the lowest luminosity quasi-simultaneous radio/X-ray detection to date for a black hole X-ray binary during its rise out of quiescence, thanks to early notification from optical monitoring combined with fast responses from sensitive multiwavelength facilities.

We investigate the relation between the parameters of the energy spectrum and the frequency and amplitude of the kilohertz quasi-periodic oscillations (kHz QPOs) in the low-mass X-ray binary 4U 1636$-$53. We fit the $3-180$-keV spectrum of this source with a model that includes a thermal Comptonisation component. We show that the frequencies of both kHz QPOs follow the same relation as a function of the parameters of this spectral component, except for a systematic frequency shift, whereas the rms fractional amplitude of each QPO follows a different relation with respect to those same parameters. This implies that, while the dynamical mechanism that sets the frequencies of the QPO can be the same for both kHz QPOs, the radiative mechanisms that set the amplitudes of the lower and the upper kHz QPO are likely different. We discuss the implications of these results to the modelling of the kHz QPOs and the possibility that the lower kHz QPO reflects a resonance between the Comptonising medium and the photons from the accretion disc and/or the neutron star surface.

We present the results of our analysis of new optical, ultraviolet (UV) and X-ray observations of a highly variable source $-$ AX J1549.8$-$5416. Both the detection of several fast rise, exponential decay outbursts in the optical light curve and the lack of He II emission lines in the optical spectra suggest AX J1549.8$-$5416 is a cataclysmic variable of the dwarf nova (DN) type. The multiwavelength analysis of three mini-outbursts and one normal outburst represent one of the most complete multiwavelength studies of a DN and help to refine the relationship between the X-ray, UV and optical emission in this system. We find that the UV emission is delayed with respect to the optical by $1.0-5.4$ days during the rising phase of the outburst. The X-ray emission is suppressed during the peak of the optical outburst and recovers during the end of the outburst. From our analysis of archival Swift, Chandra and XMM-Newton observations of AX J1549.8$-$5416, we estimate this DN has a high duty cycle ($\sim50\%$), suggesting a quiescent X-ray luminosity larger than $10^{32}$ erg/s. We also find the X-ray and UV flux are roughly anti-correlated. Furthermore, we find that, at low X-ray fluxes, the X-ray spectrum is well described by a single temperature thermal plasma model, while at high X-ray fluxes, an isobaric cooling flow model also works. We find that the maximum temperature of the plasma in quiescence is significantly higher than that in outburst.

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.

Sometimes the early star formation can be found in cold and dense molecular clouds, such as infrared dark cloud (IRDC). Considering star formation often occurs in clustered condition, HII regions may be triggering a new generation of star formation, so we can search for initial stage of massive star formation around HII regions. Based on that above, this work is to introduce one method of how to search for initial stage of massive star formation around HII regions. Towards one sample of the HII region G18.2-0.3, multiwavelength observations are carried out to investigate its physical condition. In contrast and analysis, we find three potential initial stages of massive star formation, suggesting that it is feasible to search for initial stage of massive star formation around HII regions.

Kes 79 (G33.6+0.1) is an aspherical thermal composite supernova remnant (SNR) observed across the electromagnetic spectrum and showing an unusual highly-structured morphology, in addition to harboring a central compact object (CCO). Using the CO J=1-0, J=2-1, and J=3-2 data, we provide the first direct evidence and new morphological evidence to support the physical interaction between the SNR and the molecular cloud at $V_LSR\sim 105$ km s$^{-1}$. We revisit the 380 ks XMM-Newton observations and perform a dedicated spatially resolved X-ray spectroscopic study with careful background subtraction. The overall X-ray-emitting gas is characterized by an under-ionized ($\tau_c \sim 6\times 10^{11}$ cm^${-3}$) cool ($kT_c \approx 0.20$ keV) plasma with solar abundances, plus an under-ionized ($\tau_h\sim 8\times 10^{10}$ cm$^{-3}$) hot ($kT_h\approx 0.80$ keV) plasma with elevated Ne, Mg, Si, S and Ar abundances. Kes 79 appears to have a double-hemisphere morphology viewed along the symmetric axis. Projection effect can explain the multiple shell structures and the thermal composite morphology. The X-ray filaments, spatially correlated with the 24 um IR filaments, are suggested to be due to the SNR shock interaction with dense gas, while the halo forms from SNR breaking out into a tenuous medium. The high-velocity, hot ($kT_h\sim 1.4$--1.6 keV) ejecta patch with high metal abundances, together with the non-uniform metal distribution across the SNR, indicating an asymmetric SN explosion of Kes 79. We refine the Sedov age to 4.4--6.7 kyr and the mean shock velocity to 730 km s$^{-1}$. Our multi-wavelength study suggests a progenitor mass of $\sim 15$--20 solar masses for the core-collapse explosion that formed Kes 79 and its CCO, PSR J1852+0040.

We fitted the $3-180$-keV spectrum of all the observations of the neutron-star low-mass X-ray binary 4U 1636$-$53 taken with the {\it Rossi X-ray Timing Explorer} using a model that includes a thermal Comptonisation component. We found that in the low-hard state the power-law index of this component, $\Gamma$, gradually increases as the source moves in the colour-colour diagram. When the source undergoes a transition from the hard to the soft state $\Gamma$ drops abruptly; once the source is in the soft state $\Gamma$ increases again and then decreases gradually as the source spectrum softens further. The changes in $\Gamma$, together with changes of the electron temperature, reflect changes of the optical depth in the corona. The lower kilohertz quasi-periodic oscillation (kHz QPO) in this source appears only in observations during the transition from the hard to the soft state, when the optical depth of the corona is high and changes depends strongly upon the position of the source in the colour-colour diagram. Our results are consistent with a scenario in which the lower kHz QPO reflects a global mode in the system that results from the resonance between, the disc and/or the neutron-star surface, and the Comptonising corona.

We investigated the convexity of all type I X-ray bursts with millihertz quasi-periodic oscillations (mHz QPOs) in 4U 1636-53 using archival observations with the Rossi X-ray Timing Explorer. We found that, at a 3.5 $\sigma$ confidence level, in all 39 cases in which the mHz QPOs disappeared at the time of an X-ray burst, the convexity of the burst is positive. The convexity measures the shape of the rising part of the burst light curve and, according to recent models, it is related to the ignition site of bursts on the neutron star surface. This finding suggests that in 4U 1636-53 these 39 bursts and the marginally-stable nuclear burning process responsible for the mHz QPOs take place at the neutron-star equator. This scenario could explain the inconsistency between the high accretion rate required for triggering mHz QPOs in theoretical models and the relatively low accretion rate derived from observations.

While there are many dynamical mechanisms and models that try to explain the origin and phenomenology of the quasi-periodic oscillations (QPOs) seen in the X-ray light curves of low-mass X-ray binaries, few of them address how the radiative processes occurring in these extreme environments give rise to the rich set of variability features actually observed in these light curves. A step towards this end comes from the study of the energy and frequency dependence of the phase lags of these QPOs. Here we used a methodology that allowed us to study, for the first time, the dependence of the phase lags of all QPOs in the range of 1 Hz to 1300 Hz detected in the low-mass X-ray binary 4U 1636-53 upon energy and frequency as the source changes its states as it moves through the colour-colour diagram. Our results suggest that within the context of models of up-scattering Comptonization, the phase lags dependencies upon frequency and energy can be used to extract size scales and physical conditions of the medium that produces the lags.

For an unattended telescopes in Antarctic, the remote operation, autonomous observation and control are essential. An EPICS (Experimental Physics and Industrial Control System) and RTS2(Remote Telescope System, 2nd Version) based autonomous observation and control system with remoted operation is introduced in this paper. EPICS is a set of Open Source software tools, libraries and applications developed collaboratively and used worldwide to create distributed soft real-time control systems for scientific instruments while RTS2 is an open source environment for control of a fully autonomous observatory. Using the advantage of EPICS and RTS2 respectively, a combined integrated software framework for autonomous observation and control is established that use RTS2 to fulfill the function of astronomical observation and use EPICS to fulfill the device control of telescope. A command and status interface for EPICS and RTS2 is designed to make the EPICS IOC (Input/Output Controller) components integrate to RTS2 directly. For the specification and requirement of control system of telescope in Antarctic, core components named Executor and Auto-focus for autonomous observation is designed and implemented with remote operation user interface based on Browser-Server mode. The whole system including the telescope is tested in Lijiang Observatory in Yunnan Province for practical observation to complete the autonomous observation and control, including telescope control, camera control, dome control, weather information acquisition with the local and remote operation.

Coherent oscillations and the evolution of the X-ray spectrum during thermonuclear X-ray bursts in accreting neutron-star X-ray binaries have been studied intensively but separately. We analysed all the X-ray bursts of the source 4U 1728-34 with the Rossi X-ray Timing Explorer. We found that the presence of burst oscillations can be used to predict the behaviour of the blackbody radius during the cooling phase of the bursts. If a burst shows oscillations, during the cooling phase the blackbody radius remains more or less constant for ~2 - ~8s, whereas in bursts that do not show oscillations the blackbody radius either remains constant for more than ~2 - ~8s or it shows a rapid (faster than ~2s) decrease and increase. Both the presence of burst oscillations and the time-dependent spectral behaviour of the bursts are affected by accretion rate. We also found that the rise time and convexity of the bursts' light curve are different in bursts with and without oscillations in 4U 1728--34. Bursts with oscillations have a short rise time (~0.5s) and show both positive and negative convexity, whereas bursts without oscillations have a long rise time (~1s) and mostly positive convexity. This is consistent with the idea that burst oscillations are associated with off-equator ignition.

We investigate the spectral and timing properties of the millihertz quasi-periodic oscillations (mHz QPOs) in neutron-star low-mass X-ray binary 4U 1636-53 using XMM-Newton and Rossi X-ray Timing Explorer (RXTE) observations. The mHz QPOs in the XMM-Newton/RXTE observations show significant frequency variation and disappear right before type I X-ray bursts. We find no significant correlation between the mHz QPO frequency and the temperature of the neutron-star surface, which is different from theoretical predictions. For the first time we observed the full lifetime of a mHz QPO lasting 19 ks. Besides, we also measure a frequency drift timescale around 15 ks, we speculate that this is the cooling timescale of a layer deeper than the burning depth, possibly heated by the previous burst. Moreover, the analysis of all X-ray bursts in this source shows that all type I X-ray bursts associated with the mHz QPOs are short, bright and energetic, suggesting a potential connection between mHz QPOs and He-rich X-ray bursts.

Hadron emission from supernova remnant (SNR)-molecular cloud (MC) association systems has been widely regarded as a probe of the shock-accelerated cosmic-ray protons. Here, we report on the detection of a {\gamma}-ray emission source, with a significance of 24{\sigma} in 0.2-300 GeV, projected to lie to the northwest of the thermal composite SNR Kesteven 41, using 5.6 years of Fermi-Large Area Telescope (LAT) observation data. No significant long-term variability in the energy range 0.2--300 GeV is detected around this source. The 3{\sigma} error circle, 0.09 degree; in radius, covers the 1720MHz OH maser and is essentially consistent with the location of the V_{LSR} ~-50 km/s MC with which the SNR interacts. The source emission has an exponential cutoff power-law spectrum with a photon index of 1.9+/-0.1 and a cutoff energy of 4.0+/-0.9 GeV, and the corresponding 0.2-300 GeV luminosity is ~1.3*10^36 erg/s at a distance of 12 kpc. There is no radio pulsar in the 3{\sigma} circle responsible for the high {\gamma}-ray luminosity. While the inverse Compton scattering scenario would lead to a difficulty in the electron energy budget, the source emission can naturally be explained by the hadronic interaction between the relativistic protons accelerated by the shock of SNR Kesteven 41 and the adjacent northwestern MC. In this paper, we present a list of Galactic thermal composite SNRs detected at GeV {\gamma}-ray energies by Fermi-LAT is presented.

We present an analysis of the energy and frequency dependence of the Fourier time lags of the hectoHertz quasi-periodic oscillations (QPOs) and of the QPOs at the frequency at which the power density spectrum shows a break in the neutron-star low-mass X-ray binary 4U 1636-53, using a large data set obtained with the Rossi X-ray Timing Explorer. We found that: (i) For the break frequency QPO: for low frequencies, in general the time lag is positive, but it is decreasing with increasing frequency, reaching zero lag at 20 Hz. Between 20 and 35 Hz there is a small fluctuation around zero, from where the time lags become positive again and increase slightly above zero up to 65 Hz. (ii) For the hHz QPO: we see that when the frequency is 100 Hz the time lag is negative, but it increases to zero already at 110 Hz, being consistent with this value up to 130 Hz from where it increases to 0.5 msec at around 140 Hz. From 140 Hz the time lag decreases sharply, being strongly negative for hHz greater than 220 Hz. We compare these lags with our previous results for the lags of the kiloHertz QPOs in this same source and discuss possible scenarios for producing the lags in this system in the context of reflection off the accretion disc or (up-)down-scattering in a hot medium close to the neutron star.

The Active Particle-induced X-ray Spectrometer (APXS) is one of the payloads on board the Yutu rover of Chang'E-3 mission. In order to assess the instrumental performance of APXS, a ground verification test was done for two unknown samples (basaltic rock, mixed powder sample). In this paper, the details of the experiment configurations and data analysis method are presented. The results show that the elemental abundance of major elements can be well determined by the APXS with relative deviations < 15 wt. % (detection distance = 30 mm, acquisition time = 30 min). The derived detection limit of each major element is inversely proportional to acquisition time and directly proportional to detection distance, suggesting that the appropriate distance should be < 50mm.

The physical nature of thermal composite supernova remnants (SNRs) remains controversial. We have revisited the archival XMM-Newton and Chandra data of the thermal composite SNR Kesteven 41 (Kes 41 or G337.8-0.1) and performed a millimeter observation toward this source in the $^{12}$CO, $^{13}$CO, and C$^{18}$O lines. The X-ray emission, mainly concentrated toward the southwestern part of the SNR, is characterized by distinct S and Ar He-like lines in the spectra. The X-ray spectra can be fitted with an absorbed nonequilibrium ionization collisional plasma model at a temperature of 1.3-2.6 keV and an ionization timescale of 0.1-1.2$\times$10$^{12}$ cm$^{-3}$ s. The metal species S and Ar are overabundant, with 1.2-2.7 and 1.3-3.8 solar abundances, respectively, which strongly indicate the presence of a substantial ejecta component in the X-ray-emitting plasma of this SNR. Kes 41 is found to be associated with a giant molecular cloud (MC) at a systemic local standard of rest velocity of -50 km s$^{-1}$ and confined in a cavity delineated by a northern molecular shell, a western concave MC that features a discernible shell, and an HI cloud seen toward the southeast of the SNR. The birth of the SNR in a preexisting molecular cavity implies a mass of $\gtrsim$18 M$_{\odot}$ for the progenitor if it was not in a binary system. Thermal conduction and cloudlet evaporation seem to be feasible mechanisms to interpret the X-ray thermal composite morphology, and the scenario of gas reheating by the shock reflected from the cavity wall is quantitatively consistent with the observations. An updated list of thermal composite SNRs is also presented in this paper.

Stars are born in dense cores of molecular clouds. The core mass function (CMF), which is the mass distribution of dense cores, is important for understanding the stellar initial mass function (IMF). We obtained 350 $\mu$m dust continuum data using the SHARC-II camera at the Caltech Submillimeter Observatory (CSO) telescope. A 350 $\mu$m map covering 0.25 ${deg}^{2}$ of the Ophiuchus molecular cloud was created by mosaicing 56 separate scans. The CSO telescope had an angular resolution of 9 $^{\prime\prime}$, corresponding to $1.2\times {10}^{3}\ $AU at the distance of the Ophiuchus molecular cloud (131 pc). The data was reduced using the Comprehensive Reduction Utility for SHARC-II (CRUSH). The flux density map was analyzed using the GaussClumps algorithm, within which 75 cores has been identified. We used the Spitzer c2d catalogs to separate the cores into 63 starless cores and 12 protostellar cores. By locating Jeans instabilities, 55 prestellar cores (a subcategory of starless cores) were also identified. The excitation temperatures, which were derived from FCRAO ${}^{12}$CO data, help to improve the accuracy of the masses of the cores. We adopted a Monte Carlo approach to analyze the CMF with two types of functional forms; power law and log-normal. The whole and prestellar CMF are both well fitted by a log-normal distribution, with $\mu =-1.18\pm0.10,\ \sigma =0.58\pm0.05$ and $\mu =1.40\pm0.10,\ \sigma =0.50\pm0.05$ respectively. This finding suggests that turbulence influences the evolution of the Ophiuchus molecular cloud.

We analysed all archival data of the low-mass X-ray binary 4U 1636--53 with the Rossi X-ray Timing Explorer (1490 observations). We found a total of 336 type-I X-ray bursts from this source. In the time-resolved spectra of 69 of these bursts, close to the peak of the burst, the best-fitting blackbody radius shows the sharp increase and decrease that is typical of photospheric radius-expansion (PRE) bursts. We found that in 17 of these 69 PRE bursts, after the touchdown point, the blackbody radius increases again quickly after about 1 second, and from then on the radius decreases slightly or it remains more or less constant. In the other 52 PRE bursts, after touchdown, the radius of the blackbody stays more or less constant for $\sim 2 - 8$ seconds, and after that it increases slowly. Interestingly, those PRE bursts in which the blackbody radius remains more or less constant for $\simmore 2$ seconds show coherent oscillations in the tail of the burst, whereas those PRE bursts in which the blackbody radius changes rapidly after touchdown show no coherent oscillations in the tail of the burst. We found that the distribution of durations of the post touchdown phase between these two groups of PRE bursts is significantly different; the Kolmogorov-Smirnov probability that the two groups of PRE bursts come from the same parent populations is only $3.5 \times 10^{-7}$. This is the first time that the presence of burst oscillations in the tail of X-ray bursts is associated with a systematic behaviour of the spectral parameters in that phase of the bursts. This result is consistent with predictions of models that associate the oscillations in the tail of X-ray bursts with the propagation of a cooling wake in the material on the neutron-star surface during the decay of the bursts.

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.

We analysed the complete set of bursts from the neutron-star low-mass X-ray binary 4U 1820-30 detected with the Rossi X-ray Timing Explorer (RXTE). We found that all are photospheric radius expansion bursts, and have similar duration, peak flux and fluence. From the analysis of time-resolved spectra during the cooling phase of the bursts, we found that the relation between the bolometric flux and the temperature is very different from the canonical F - T^4 relation that is expected if the apparent emitting area on the surface of the neutron star remains constant. The flux-temperature relation can be fitted using a broken power law, with indices 2.0$\pm$0.3 and 5.72$\pm$0.06. The departure from the F - T^4 relation during the cooling phase of the X-ray bursts in 4U 1820-30 could be due to changes in the emitting area of the neutron star while the atmosphere cools-down, variations in the colour-correction factor due to chemical evolution, or the presence of a source of heat, e.g. residual hydrogen nuclear burning, playing an important role when the burst emission ceases.

Time-resolved spectra during the cooling phase of thermonuclear X-ray bursts in low-mass X-ray binaries (LMXBs) can be used to measure the radii and masses of neutron stars. We analyzed ~ 300 bursts of the LMXB 4U 1636-53 using data from the Rossi X-ray Timing Explorer. We divided the bursts in three groups, photospheric radius expansion (PRE), hard non-PRE and soft non-PRE bursts, based on the properties of the bursts and the state of the source at the time of the burst. For the three types of bursts, we found that the average relation between the bolometric flux and the temperature during the cooling phase of the bursts is significantly different from the canonical $F \propto T^4$ relation that is expected if the apparent emitting area on the surface of the neutron star remains constant as the flux decreases during the decay of the bursts. We also found that a single power law cannot fit the average flux-temperature relation for any of the three types of bursts, and that the flux-temperature relation for the three types of bursts is significantly different. Finally, for the three types of bursts, the temperature distribution at different flux levels during the decay of the bursts is significantly different. From the above we conclude that hard non-PRE bursts ignite in a hydrogen-rich atmosphere, whereas for soft non-PRE and PRE bursts the fuel is helium-rich. We further conclude that the metal abundance in the neutron star atmosphere decreases as the bursts decay, probably because the heavy elements sink faster in the atmosphere than H and He.

Recently, the neutron star X-ray binary EXO 0748-676 underwent a transition to quiescence. We analyzed an XMM-Newton observation of this source in quiescence, where we fitted the spectrum with two different neutron-star atmosphere models. From the fits we constrained the allowed parameter space in the mass-radius diagram for this source for an assumed range of distances to the system. Comparing the results with different neutron-star equations of state, we constrained the distance to EXO 0748-676. We found that the EOS model 'SQM1' is rejected by the atmosphere model fits for the known distance, and the 'AP3' and 'MS1' is fully consistent with the known distance.

We have discovered a triple-peaked X-ray burst from the low-mass X-ray binary (LMXB) 4U 1636-53 with the Rossi X-ray Timing Explorer (RXTE). This is the first triple-peaked burst reported from any LMXB using RXTE, and it is only the second burst of this kind observed from any source. (The previous one was also from 4U 1636-53, and was observed with EXOSAT.) From fits to time-resolved spectra, we find that this is not a radius-expansion burst, and that the same triple-peaked pattern seen in the X-ray light curve is also present in the bolometric light curve of the burst. Similar to what was previously observed in double-peaked bursts from this source, the radius of the emitting area increases steadily during the burst, with short periods in between during which the radius remains more or less constant. The temperature first increases steeply, and then decreases across the burst also showing three peaks. The first and last peak in the temperature profile occur, respectively, significantly before and significantly after the first and last peaks in the X-ray and bolometric light curves. We found no significant oscillations during this burst. This triple-peaked burst, as well as the one observed with EXOSAT and the double-peak bursts in this source, all took place when 4U 1636-53 occupied a relatively narrow region in the colour-colour diagram, corresponding to a relatively high (inferred) mass-accretion rate. No model presently available is able to explain the multiple-peaked bursts.

We report the detections of the anti-correlated soft and hard X-rays, and the time lags of $\sim$ hecto-second from the neutron star low-mass X-ray binary Cyg X-2, a well-known Z-type luminous source. Both the anti-correlation and the positive correlation were detected during the low-intensity states, while only the latter showed up during high-intensity states. Comparing with the lower part of normal branch and flaring branch, more observations located on the horizontal and the upper normal branches are accompanied with the anti-correlation, implying the occurrence of the anti-correlation under circumstance of a low mass accretion rate. So far the anti-correlated hard lag of thousand-second timescale are only reported from the Galactic black hole candidates in their hard states. Here we provide the first evidence that a similar feature can also establish in a neutron-star system like Cyg X-2. Finally, the possible origins of the observed time lags are discussed under the current LMXB models.

We have carried out detailed analysis on the black hole candidate (BHC) X- ray transient Swift J1753.5-0127 observed by the Rossi X-Ray Timing Explorer (RXTE) during its outburst in 2005 {2006. The spectral analysis shows that the emissions are dominated by the hard X-rays, thus revealing the low/hard state of the source during the outburst. The peak luminosity is found lower than the typical value of balancing the mass flow and evaporation of the inner edge of disk (Meyer-Hofmeister 2004). As a result, the disk is prevented from extending inward to produce strong soft X-rays, corresponding to the so-called high/soft state. These are the typical characteristics for a small subset of BHCs, i.e. those soft X-ray transients stay at the low/hard state during the outburst. In most observational time, the QPO frequencies are found to vary roughly linearly with the fluxes and the spectral indices, while the deviation from this relationship at the peak luminosity might provide the first observational evidence of a partially evaporated inner edge of the accretion disk. The anti-correlation between the QPO frequency and spectral color suggests that the global disk oscillation model proposed by Titarchuk & Osherovich (2000) is not likely at work.

Using the SU($N$) representation of the group theory, we derive the general form of the spin swapping operator for the quantum Heisenberg spin-$s$ systems. We further prove that such a spin swapping operator is equal to the spin singlet pairing operator under the partial transposition. For SU(2) invariant states, it is shown that the expectation value of the spin swapping operator and its generalizations, the permutations, can be used as an entanglement witness, especially, for the formulation of observable conditions of entanglement.