In 2015 October, the Be/X-ray binary 4U 0115+63 underwent a type II outburst, reaching an X-ray luminosity of $\sim$ {\DP{38}\es}. During the outburst, Nuclear Spectroscopic Telescope Array ({\it NuSTAR}) performed two Target of Opportunity observations. Using the broad-band spectra from {\it NuSTAR} (3$-$79 keV), we have detected multiple cyclotron lines of the source, i.e., $\sim$ 12, 16, 22 and 33/35 keV. Obviously, the 16 keV line is not a harmonic component of the 12 keV line. As described by the phase-dependent equivalent widths of these cyclotron lines, the 16 keV and 12 keV lines are two different fundamental lines. In our work, we apply the two-poles cyclotron line model to the observation, i.e., the two line sets are formed at the same altitude ($\sim$ 0.2 km over the NS surface) of different magnetic poles, with $\sim$ {1.1\TDP{12}} and {1.4\TDP{12}} G in two poles, respectively.

In 2015 October, the Be/X-ray binary 4U 0115+63 underwent a type II outburst, reaching an X-ray luminosity of $\sim$ {\DP{38}\es}. During the outburst, Nuclear Spectroscopic Telescope Array ({\it NuSTAR}) performed two Target of Opportunity observations. Using the broad-band spectra from {\it NuSTAR} (3$-$79 keV), we have detected multiple cyclotron lines of the source, i.e., $\sim$ 12, 16, 22 and 33/35 keV. Obviously, the 16 keV line is not a harmonic component of the 12 keV line. As described by the phase-dependent equivalent widths of these cyclotron lines, the 16 keV and 12 keV lines are two different fundamental lines. In our work, we apply the two-poles cyclotron line model to the observation, i.e., the two line sets are formed at the same altitude ($\sim$ 0.2 km over the NS surface) of different magnetic poles, with $\sim$ {1.1\TDP{12}} and {1.4\TDP{12}} G in two poles, respectively.

The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, the search for similar emission at optical wavelengths has been hampered by the high inferred extinction on the line of sight. Here, we present results from the first search for second-timescale emission from the source at near-infrared wavelengths using the Palomar Gattini-IR observing system in J-band, made possible by a recently implemented detector read-out mode that allowed for short exposure times of 0.84 s with 99.9% observing efficiency. With a total observing time of 12 hours (47728 images) on source, we place median $3\,\sigma$ limits on the second-timescale emission of $< 20$ mJy (13.1 AB mag). We present non-detection limits from epochs of four simultaneous X-ray bursts detected by the Insight-{\it HXMT} and {\it NuSTAR} telescopes during our observing campaign. The limits translate to an extinction corrected fluence limit of $< 125$ Jy ms for an estimated extinction of $A_J = 2.0$ mag. These limits provide the most stringent constraints to date on the fluence of flares at frequencies of $\sim 10^{14}$ Hz, and constrain the ratio of the near-infrared (NIR) fluence to that of coincident X-ray bursts to $R_{\rm NIR} < 2.5 \times 10^{-2}$. Our observations were sensitive enough to easily detect a near-infrared counterpart of FRB 200428 if the NIR emission falls on the same power law as that observed across its radio to X-ray spectrum. The non-detection of NIR emission around the coincident X-ray bursts constrains the fluence index of the brightest burst to be steeper than $0.35$.

Insight-HXMT is the first Chinese X-ray astronomical mission, launched successfully on June 15, 2017, from China's Jiuquan Satellite Launch Center. Insight-HXMT was designed to have a broad energy coverage in X-rays, from 1-250 keV, with excellent timing and adequate energy resolution at soft X-rays, and the largest effective area at hard X-rays. Here, we present a collection of papers of the Journal of High Energy Astrophysics on the Early Results of China's 1st X-ray Astronomy Satellite Insight-HXMT. These papers cover the in-orbit performance, the background model, and all calibration results, together with several first results on observations of binaries and details on the Galactic plane survey.

In this brief report, we investigate the existence of 4-dimensional static spherically symmetric black holes (BHs) in the Einstein-complex-scalar-Gauss-Bonnet (EcsGB) gravity with an arbitrary potential $V(\phi)$ and a coupling $f(\phi)$ between the scalar field $\phi$ and the Gauss-Bonnet (GB) term. We find that static regular BH solutions with complex scalar hairs do not exist. This conclusion does not depend on the coupling between the GB term and the scalar field, nor on the scalar potential $V(\phi)$ and the presence of a cosmological constant $\Lambda$ (which can be either positive or negative), as longer as the scalar field remains complex and is regular across the horizon.

M31N 2017-09a is a classical nova and was observed for some 160 days following its initial eruption, during which time it underwent a number of bright secondary outbursts. The light-curve is characterized by continual variation with excursions of at least 0.5 magnitudes on a daily time-scale. The lower envelope of the eruption suggests that a single power-law can describe the decline rate. The eruption is relatively long with $t_2 = 111$, and $t_3 = 153$ days.

Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

We report on detailed spectral analysis of emission from X-ray pulsar 4U 1901+03 using the data obtained by the Insight-HXMT and NuSTAR observatories during the 2019 outburst of the source. Thanks to the extensive coverage of the outburst by Insight-HXMT, we were able to investigate spectral evolution of the source as a function of flux, and compare these results to the previous reports with focus on the properties of a putative absorption feature at around 10 keV. In particular, we demonstrate that the broadband X-ray continuum of 4U 1901+03 can be well described with a two-component continuum model without absorption line at 10 keV, which casts doubt on its interpretation as a cyclotron line. The high quality of the data has also allowed us to perform both phase-averaged and phase-resolved spectral analysis as function of luminosity. Finally, we performed a detailed investigation of another absorption feature in the spectrum of the source around 30 keV recently reported in the NuSTAR data. In line with the previous findings, we show that this feature appears to be significantly detected, albeit only in the pulse-phase resolved spectra at relatively high luminosities. The non-detection of the line in the averaged spectrum is possibly associated with the observed strong variations of its parameters with the pulse phase.

We report on detailed spectral analysis of emission from X-ray pulsar 4U 1901+03 using the data obtained by the Insight-HXMT and NuSTAR observatories during the 2019 outburst of the source. Thanks to the extensive coverage of the outburst by Insight-HXMT, we were able to investigate spectral evolution of the source as a function of flux, and compare these results to the previous reports with focus on the properties of a putative absorption feature at around 10 keV. In particular, we demonstrate that the broadband X-ray continuum of 4U 1901+03 can be well described with a two-component continuum model without absorption line at 10 keV, which casts doubt on its interpretation as a cyclotron line. The high quality of the data has also allowed us to perform both phase-averaged and phase-resolved spectral analysis as function of luminosity. Finally, we performed a detailed investigation of another absorption feature in the spectrum of the source around 30 keV recently reported in the NuSTAR data. In line with the previous findings, we show that this feature appears to be significantly detected, albeit only in the pulse-phase resolved spectra at relatively high luminosities. The non-detection of the line in the averaged spectrum is possibly associated with the observed strong variations of its parameters with the pulse phase.

We discuss the potential of the eXTP X-ray telescope, in particular its Spectroscopic Focusing Array (SFA), Large Area Detector (LAD) and Wide Field Monitor (WFM) for the detection of a signal from keV-scale decaying dark matter. We show that the sensitivity of the eXTP is sufficient to improve existing constraints on the mixing angle of the neutrino Minimal extension of the Standard Model ($\nu$MSM) by a factor of 5-10 within the dark matter mass range 2-50 keV, assuming 1% level of systematic uncertainty. We assert that the eXTP will be able to probe previously inaccessible range of $\nu$MSM parameters and serve as a precursor for the Athena mission in decaying dark matter searches.

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.

Fast Radio Bursts (FRBs) are bright radio transients with millisecond duration at cosmological distances. Since compact dark matter/objects (COs) could act as lenses and cause split of this kind of very short duration signals, Mu$\rm{\tilde{n}}$oz et al. (2016) has proposed a novel method to probe COs with lensing of FRBs. In this Letter, we for the first time apply this method to real data and give constraints of the nature of COs with currently available FRB observations. We emphasize the information from dynamic spectra of FRBs is quite necessary for identifying any lensed signals and find no echoes in the existing data. The null search gives a constraint comparable to that from galactic wide binaries, though the methods of redshift inference from dispersion measure would impact a little. Furthermore, we make an improved forecast basing on the distributions of real data for the ongoing and upcoming telescopes. Finally, we discuss the situation where one or more lensed signals will be detected. In such a case, the parameter space of COs can be pinned down very well since the lens mass can be directly determined through the observed flux ratio and time delay between split images.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We present long-term observations of two intermittent pulsars, PSRs~J1832+0029 and J1841$-$0500 using the Parkes 64\,m radio telescope. The radio emission for these pulsars switches "off" for year-long durations. Our new observations have enabled us to improve the determination of the on-off timescales and the spin down rates during those emission states. In general our results agree with previous studies of these pulsars, but we now have significantly longer data spans. We have identified two unexpected signatures in the data. Weak emission was detected in a single observation of PSR~J1832$+$0029 during an "off" emission state. For PSR~J1841$-$0500, we identified a quasi-periodic fluctuation in the intensities of the detectable single pulses, with a modulation period between 21 and 36 pulse periods.

Bramante, Delgado, and Martin [Phys. Rev. D96, 063002(2017)., hereafter BDM17] extended the analytical formalism of dark matter (DM) capture in a very important way, which allows, in principle, the use of compact astrophysical objects, such as neutron stars (NS), as dark matter detectors. In this comment, we point out the existence of a region in the dark matter neutron scattering cross section $(\sigma_{nX})$ vs. dark matter mass $(m_{X})$ where the constraining power of this method is lost. This corresponds to a maximal temperature ($T_{crit}$) the NS has to have, in order to serve as a dark matter detector. In addition, we point out several typos and errors in BDM17 that do not affect drastically their conclusions. Moreover, we provide semi-analytical approximations for the total capture rates of dark matter particle of arbitrary mass for various limiting regimes. Those analytical approximations are used to validate our numerical results.

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

In big data era, the special data with rare characteristics may be of great significations. However, it is very difficult to automatically search these samples from the massive and high-dimensional datasets and systematically evaluate them. The DoPS, our previous work [2], provided a search method of rare spectra with double-peaked profiles from massive and high-dimensional data of LAMOST survey. The identification of the results is mainly depended on visually inspection by astronomers. In this paper, as a follow-up study, a new lattice structure named SVM-Lattice is designed based on SVM(Support Vector Machine) and FCL(Formal Concept Lattice) and particularly applied in the recognition and evaluation of rare spectra with double-peaked profiles. First, each node in the SVM-Lattice structure contains two components: the intents are defined by the support vectors trained by the spectral samples with the specific characteristics, and the relevant extents are all the positive samples classified by the support vectors. The hyperplanes can be extracted from every lattice node and used as classifiers to search targets by categories. A generalization and specialization relationship is expressed between the layers, and higher layers indicate higher confidence of targets. Then, including a SVM-Lattice building algorithm, a pruning algorithm based on association rules, and an evaluation algorithm, the supporting algorithms are provided and analysed. Finally, for the recognition and evaluation of spectra with double-peaked profiles, several data sets from LAMOST survey are used as experimental dataset. The results exhibit good consistency with traditional methods, more detailed and accurate evaluations of classification results, and higher searching efficiency than other similar methods.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

In Dec. 2011 PSR B0540-69 experienced a spin-down rate transition (SRT), after which the spin-down power of the pulsar increased by ~36%. About 1000 days after the SRT, the X-ray luminosity of the associated pulsar wind nebula (PWN) was found to brighten by 32+/-8%. After the SRT, the braking index n of PSR B0540-69 changes from n=2.12 to n=0.03 and then keeps this value for about five years before rising to n=0.9 in the following years. We find that most of the current models have difficulties in explaining the measured braking index. One exceptive model of the braking index evolution is the increasing dipole magnetic field of PSR B0540-69. We suggest that the field increase may result from some instabilities within the pulsar core that enhance the poloidal component at the price of toroidal component of the magnetic field. The increasing dipole magnetic field will result in the X-ray brightening of the PWN. We fit the PWN X-ray light curve by two models: one assumes a constant magnetic field within the PWN during the brightening and the other assumes an enhanced magnetic field proportional to the energy density of the PWN. It appears that the two models fit the data well, though the later model seems to fit the data a bit better. This provides marginal observational evidence that magnetic field in the PWN is generated by the termination shock. Future high-quality and high-cadence data are required to draw a solid conclusion.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We present the results of a ~250 square degrees CO mapping (+26d<l<+50d and -5d<b<+5d) toward the Aquila Rift region at a spatial resolution of ~50" and a grid spacing of 30". The high dynamic range CO maps with a spectral resolution of ~0.2km/s display highly structured molecular cloud (MC) morphologies with valuable velocity information, revealing complex spatial and dynamical features of the local molecular gas. In combination with the MWISP CO data and the Gaia DR2, distances of the main MC structures in the local ISM are well determined toward the Aquila Rift. We find that the total MC mass within 1 kpc is about >4.1x10^5 Msun in the whole region. In fact, the mass of the molecular gas is dominated by the W40 giant molecular cloud (GMC) at ~474 pc (~1.4x10^5 Msun) and the GMC complex G036.0+01.0 at ~560-670 pc (~2.0x10^5 Msun), while the MCs at ~220-260pc have gas masses of 10^2-10^3 Msun. Interestingly, an ~80pc long filamentary MC G044.0-02.5 at a distance of ~404 pc shows a systematic velocity gradient along and perpendicular to the major axis of the filament. The HI gas with the enhanced emission has the similar spatial morphologies and velocity features compared to the corresponding CO structure, indicating that the large-scale converging HI flows are probably responsible for the formation of the MC. Meanwhile, the long filamentary MC consists of many sub-filaments with the lengths ranging from ~0.5 pc to several pc, as well as prevalent networks of filaments in other large-scale local MCs.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

In this paper, we systematically study spherically symmetric static spacetimes in the framework of Einstein-aether theory, and pay particular attention to the existence of black holes (BHs). In the theory, one scalar and one vector two gravitational modes additionally appear. To avoid the vacuum gravi-\v{C}erenkov radiation, they must all propagate with speeds greater than or at least equal to the speed of light. In the spherical case, only the scalar mode is relevant, so BH horizons are defined by this mode. In the present studies we first clarify several subtle issues. In particular, we find that, out of the five non-trivial field equations, only three are independent, so the problem is well-posed, as now there are only three unknown functions, two, $F(r)$ and $B(r)$, describe the spacetime and one, $A(r)$, describes the aether field. In addition, the two second-order differential equations for $A$ and $F$ are independent of $B$, and once they are found, $B$ is given simply by an algebraic expression of $F,\; A$ and their derivatives. To simplify the problem further, we explore the symmetry of field redefinitions, and work first with the redefined metric and aether field, and then obtain the physical ones by the inverse transformations. These clarifications significantly simplify the computational labor, which is important, as the problem is highly involved mathematically. In fact, it is exactly because of these, we find various numerical BH solutions with an accuracy that is at least two orders higher than previous ones. More important, these BH solutions are the only ones that satisfy the self-consistent conditions and meantime are consistent with all the observational constraints obtained so far. The locations of universal horizons are also identified, together with several other observationally interesting quantities.

Charge exchange (CX) emission reveals the significant interaction between neutral and ionized interstellar medium (ISM) components of the dense, multiphase, circumnuclear region of a galaxy. We use a model including a thermal and a CX components to describe the high-resolution XMM-Newton/RGS spectrum of the diffuse emission in the central region of M51. Representative signatures of CX emission -- especially the prominent OVII forbidden line and the excess emission in the OVIII Ly$\gamma$ lines -- can be well explained by the model. Combined with the Chandra images in the OVIII and the OVII bands, we find the soft X-ray emission is dominated by the jet-driven outflow and its interaction with the ambient neutral material. The jet-driven outflow itself is likely a thermal plasma of $\sim 0.59$ keV, with mostly sub-solar abundances. It runs into the ambient neutral gas, and produces significant CX emission that accounts for one-fifth of the diffuse X-ray emission in the 7--28 {\AA} band. The effective interface area in the CX process is one order of magnitude greater than the geometrical surface area of the jet-driven outflow. The tenuous outflow driven by the nuclear star formation may also contribute a small portion to both the diffuse thermal and CX emission. The photoionization by the active galactic nuclei (AGNs) and the resonance scattering by the hot gas itself are disfavored, though the effects from past AGN events may not be ruled out.

Massive stars and their associated ionized (HII) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H regions are still poorly known. To investigate the impact of HII regions on the formation of filaments, we imaged the Galactic HII region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. We used the ArT\'eMiS camera on the APEX telescope and combined the ArT\'eMiS data at 350 and 450 microns with Herschel-SPIRE/HOBYS. We studied the dense gas distribution around RCW 120 with a resolution of 8 arcsec (0.05 pc at a distance of 1.34 kpc). Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the HII region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer like profile with a deconvolved Gaussian FWHM of 0.09 pc, as observed for filaments in low-mass star-forming regions. This study suggests that compression exerted by HII regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArT\'eMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure

We present a survey of molecular outflows across the dark cloud complex in the Cygnus region, based on 46.75 deg^2 field of CO isotopologues data from Milky Way Imaging Scroll Painting (MWISP) survey. A supervised machine learning algorithm, Support Vector Machine (SVM), is introduced to accelerate our visual assessment of outflow features in the data cube of 12CO and 13CO J = 1-0 emission. A total of 130 outflow candidates are identified, of which 77 show bipolar structures and 118 are new detections. Spatially, these outflows are located inside dense molecular clouds and some of them are found in clusters or in elongated linear structures tracing the underlying gas filament morphology. Along the line of sight, 97, 31, and 2 candidates reside in the Local, Perseus, and Outer arm, respectively. Young stellar objects as outflow drivers are found near most outflows, while 36 candidates show no associated source. The clusters of outflows that we detect are inhomogeneous in their properties; nevertheless, we show that the outflows cannot inject turbulent energy on cloud scales. Instead, at best, they are restricted to affecting the so called "clump" and "core" scales, and this only on short (~0.3 Myr) estimated timescales. Combined with outflow samples in the literature, our work shows a tight outflow mass-size correlation.

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

The role of ionization feedback on high-mass (>8 Msun) star formation (HMSF) is still highly debated. Questions remain concerning the presence of nearby HII regions changes the properties of early HMSF and whether HII regions promote or inhibit the formation of high-mass stars. To characterize the role of HII regions on the HMSF, we study the properties of a sample of candidates high-mass starless clumps (HMSCs), of which about 90% have masses larger than 100 Msun. These high-mass objects probably represent the earliest stages of HMSF; we search if (and how) their properties are modified by the presence of an HII region. We took advantage of the recently published catalog of HMSC candidates. By cross matching the HMSCs and HII regions, we classified HMSCs into three categories: 1) The HMSCs associated with HII regions both in the position in the projected plane of the sky and in velocity; 2) HMSCs associated in the plane of the sky, but not in velocity; and 3) HMSCs far away from any HII regions in the projected sky plane. We carried out comparisons between associated and nonassociated HMSCs based on statistical analyses of multiwavelength data from infrared to radio. Statistical analyses suggest that HMSCs associated with HII regions are warmer, more luminous, more centrally-peaked and turbulent. We also clearly show, for the first time, that the ratio of bolometric luminosity to envelope mass of HMSCs (L/M) could not be a reliable evolutionary probe for early HMSF due to the external heating effects of the HII regions. More centrally peaked and turbulent properties of HMSCs associated with HII regions may promote the formation of high-mass stars by limiting fragmentation. High resolution interferometric surveys toward HMSCs are crucial to reveal how HII regions impact the star formation process inside HMSCs.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 ${\mu}$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files.

The Galactic plane scanning survey is one of the main scientific objectives of the Hard X-ray Modulation Telescope (known as Insight-HXMT). During the two-year operation of Insight-HXMT, more than 1000 scanning observations have been performed and the whole Galactic plane ($\rm 0^{\circ}<l<360^{\circ}$, $\rm -10^{\circ}<b<10^{\circ}$) has been covered completely. We summarize the Galactic plane scanning survey of Insight-HXMT for two years, including the characteristics of the scanning data, the data analysis process and the preliminary results of the Low-Energy telescope, the Medium-Energy telescope and the High-Energy telescope. With the light curve PSF fitting method, the fluxes of the known sources in the scanned area as well as the flux errors are obtained for each scanning observation. From the relationships of SNRs and fluxes, the $5\sigma$ sensitivities of three telescopes of Insight-HXMT are estimated as $\rm \sim7.6\times10^{-11}~erg cm^{-2}~s^{-1}$ ($\rm 3 mCrab,~1-6 keV$), $\rm \sim4.0\times10^{-10}~erg~cm^{-2}~s^{-1}$ ($\rm 20~mCrab,~7-40 keV$) and $\rm \sim2.6\times10^{-10}~erg cm^{-2}~s^{-1}$ ($\rm 18 mCrab,~25-100 keV$) for an individual scanning observation of $2-3$ hours, respectively. Up to September 2019, more than 800 X-ray sources with various types are monitored by the three telescopes and their long-term light curves with three energy bands are obtained to make further scientific analyses.

The Galactic plane scanning survey is one of the main scientific objectives of the Hard X-ray Modulation Telescope (known as Insight-HXMT). During the two-year operation of Insight-HXMT, more than 1000 scanning observations have been performed and the whole Galactic plane ($\rm 0^{\circ}<l<360^{\circ}$, $\rm -10^{\circ}<b<10^{\circ}$) has been covered completely. We summarize the Galactic plane scanning survey of Insight-HXMT for two years, including the characteristics of the scanning data, the data analysis process and the preliminary results of the Low-Energy telescope, the Medium-Energy telescope and the High-Energy telescope. With the light curve PSF fitting method, the fluxes of the known sources in the scanned area as well as the flux errors are obtained for each scanning observation. From the relationships of SNRs and fluxes, the $5\sigma$ sensitivities of three telescopes of Insight-HXMT are estimated as $\rm \sim7.6\times10^{-11}~erg cm^{-2}~s^{-1}$ ($\rm 3 mCrab,~1-6 keV$), $\rm \sim4.0\times10^{-10}~erg~cm^{-2}~s^{-1}$ ($\rm 20~mCrab,~7-40 keV$) and $\rm \sim2.6\times10^{-10}~erg cm^{-2}~s^{-1}$ ($\rm 18 mCrab,~25-100 keV$) for an individual scanning observation of $2-3$ hours, respectively. Up to September 2019, more than 800 X-ray sources with various types are monitored by the three telescopes and their long-term light curves with three energy bands are obtained to make further scientific analyses.

Gravitational accretion accumulates the original mass, and this process is crucial for us to understand the initial phases of star formation. Using the specific infall profiles in optically thick and thin lines, we searched the clumps with infall motion from the Milky Way Imaging Scroll Painting (MWISP) CO data in previous work. In this study, we selected 133 sources of them as a sub-sample for further research and identification. The excitation temperatures of these sources are between 7.0 and 38.5 K, while the H_2 column densities are between 10^21 and 10^23 cm^-2. We have observed optically thick lines HCO+ (1-0) and HCN (1-0) using the DLH 13.7-m telescope, and found 56 sources of them with blue profile and no red profile in these two lines, which are likely to have infall motions, with the detection rate of 42\%. It suggests that using CO data to restrict sample can effectively improve the infall detection rate. Among these confirmed infall sources, there are 43 associated with Class 0/I young stellar objects (YSOs), and 13 are not. These 13 sources are probably associated with the sources in earlier evolutionary stage. By comparison, the confirmed sources which are associated with Class 0/I YSOs have higher excitation temperatures and column densities, while the other sources are colder and have lower column densities. Most infall velocities of the sources we confirmed are between 10^-1 to 10^0 km s^-1, which is consistent with previous studies.

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

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.

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.

One of the most unique phenomena in the Galactic center region is the existence of numerous long and narrow filamentary structures within a few hundred parsecs of Sgr A$^{\star}$. While more than one than one hundred radio filaments have been revealed by MeerKAT, about two dozens X-ray filaments have been discovered so far. In this article, we report our analysis on the deep Chandra and NuSTAR observations of a non-thermal X-ray filament, G0.13-0.11, which is located adjacent to the Radio arc. Chandra revealed a unique morphology of G0.13-0.11, which is an elongated (0.1 pc in width and 3.2 pc in length) structure slightly bended towards the Radio arc. A pulsar candidate ($\Gamma \sim 1.4$) is detected in the middle of the filament, with a tail of diffuse non-thermal X-ray emission on one side of the filament. The filament is detected by NuSTAR up to 79 keV, with the hard X-ray centroid consistent with the pulsar candidate. We found that the X-ray intensity decays along the filament farther away from the pulsar candidate, dropping to half of its peak value at 2.2 pc away. This system is mostly likely a Pulsar Wind Nebula interacting with ambient interstellar magnetic field, where the filaments are kinetic jets from PWN as recently proposed. The nature of this filament adds to complex origin of the X-ray filaments, which serve as powerful tools to probe local and global powerful particle accelerators in the Galactic center.

Testing gravity by binary pulsars nowadays becomes a key issue. Screened modified gravity is a kind of scalar-tensor theory with screening mechanism in order to satisfy the tight Solar System tests. In this paper, we investigate how the screening mechanism affects the orbital dynamics of binary pulsars, and calculate in detail the five post-Keplerian (PK) parameters in this theory. These parameters differ from those of general relativity (GR), and the differences are quantified by the scalar charges, which lead to the dipole radiation in this theory. We combine the observables of PK parameters for the ten binary pulsars, respectively, to place the constraints on the scalar charges and possible deviations from GR. The dipole radiation in the neutron star (NS) - white dwarf (WD) binaries leads to more stringent constraints on deviations from GR. The most constraining systems for the scalar charges of NS and WD are PSR~B1913$+$16 and PSR~J1738$+$0333, respectively. The results of all tests exclude significant strong-field deviations and show good agreement with GR.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

We have carried out two-dimensional hydrodynamical simulations to study the effects of disk self-gravity and radiative cooling on the formation of gaps and spirals. (1) With disk self-gravity included, we find stronger, more tightly-wound spirals and deeper gaps in more massive disks. The deeper gaps are due to the larger Angular Momentum Flux (AMF) of the waves excited in more massive disks, as expected from the linear theory. The position of the secondary gap does not change, provided that the disk is not extremely massive ($Q \gtrsim 2$). (2) With radiative cooling included, the excited spirals become monotonically more open (less tightly-wound) as the disk's cooling timescale increases. On the other hand, the amplitude and strength of the spirals decrease when the cooling time increases from a small value to $\sim 1/\Omega$, but then the amplitude starts to increase again when the cooling time continues to increase. This indicates that radiative dissipation becomes important for waves with $T_{cool}\sim$ 1. Consequently, the induced primary gap is narrower and the secondary gap becomes significantly shallower when the cooling time becomes $\sim 1/\Omega$. When the secondary gap is present, the position of it moves to the inner disk from the fast cooling cases to the slow cooling cases. The dependence of gap properties on the cooling timescale (e.g. in AS 209) provides a new way to constrain the disk optical depth and thus disk surface density.

As smaller analogs of Active Galactic Nuclei (AGN), X-ray Binaries (XRBs) are also capable of launching jets that accelerate particles to high energies. In this work, we reexamine XRB jets as potential sources of high-energy cosmic rays (CRs) and explore whether they could provide a significant second Galactic component to the CR spectrum. In the most intriguing scenario, XRB-CRs could dominate the observed spectrum above the so-called "knee" feature at $\sim 3 \times 10^{15}$ eV, offering an explanation for several key issues in this transition zone from Galactic to extragalactic CRs. We discuss how such a scenario could be probed in the near future via multi-messenger observations of XRB jets, as well as diffuse Galactic neutrino flux measurements.

We make the in-orbit calibration to the point-spread functions (PSFs) of the collimators of the Hard X-ray Modulation Telescope with the scanning observation of the Crab. We construct the empirical adjustments to the theoretically calculated geometrical PSFs. The adjustments contain two parts: a rotating matrix to adjust the directional deviation of the collimators and a paraboloidal function to correct the inhomogeneity of the real PSFs. The parameters of the adjusting matrices and paraboloidal functions are determined by fitting the scanning data with lower scanning speed and smaller intervals during the calibration observations. After the PSF calibration, the systematic errors in source localization in the Galactic plane scanning survey are 0.010 deg, 0.015 deg, 0.113 deg for the Low-Energy Telescope (LE), the Medium-Energy telescope (ME) and the High-Energy telescope (HE), respectively; meanwhile, the systematic errors in source flux estimation are 1.8%, 1.6%, 2.7% for LE, ME and HE, respectively.

We make the in-orbit calibration to the point-spread functions (PSFs) of the collimators of the Hard X-ray Modulation Telescope with the scanning observation of the Crab. We construct the empirical adjustments to the theoretically calculated geometrical PSFs. The adjustments contain two parts: a rotating matrix to adjust the directional deviation of the collimators and a paraboloidal function to correct the inhomogeneity of the real PSFs. The parameters of the adjusting matrices and paraboloidal functions are determined by fitting the scanning data with lower scanning speed and smaller intervals during the calibration observations. After the PSF calibration, the systematic errors in source localization in the Galactic plane scanning survey are 0.010 deg, 0.015 deg, 0.113 deg for the Low-Energy Telescope (LE), the Medium-Energy telescope (ME) and the High-Energy telescope (HE), respectively; meanwhile, the systematic errors in source flux estimation are 1.8%, 1.6%, 2.7% for LE, ME and HE, respectively.

Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von K\'arm\'an crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientific goal is to obtain first active dosimetric measurements on the surface of the Moon. LND also provides observations of fast neutrons which are a result of the interaction of high-energy particle radiation with the lunar regolith and of their thermalized counterpart, thermal neutrons, which are a sensitive indicator of subsurface water content.

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $\alpha$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

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.

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.

The fundamental stellar atmospheric parameters T_eff and log g and 13 chemical abundances are derived for medium-resolution spectroscopy from LAMOST Medium-Resolution Survey (MRS) data sets with a deep-learning method. The neural networks we designed, named as SPCANet, precisely map LAMOST MRS spectra to stellar parameters and chemical abundances. The stellar labels derived by SPCANet are with precisions of 119 K for T_eff and 0.17 dex for log g. The abundance precision of 11 elements including [C/H], [N/H], [O/H], [Mg/H], [Al/H], [Si/H], [S/H], [Ca/H], [Ti/H], [Cr/H], [Fe/H], and [Ni/H] are 0.06~0.12 dex, while of [Cu/H] is 0.19 dex. These precisions can be reached even for spectra with signal-to-noise as low as 10. The results of SPCANet are consistent with those from other surveys such as APOGEE, GALAH and RAVE, and are also validated with the previous literature values including clusters and field stars. The catalog of the estimated parameters is available at \url{http://paperdata.china-vo.org/LAMOST/MRS_parameters_elements.csv}.

LS 5039 is a well-known $\gamma$-ray binary system which consists of an unknown compact object and a massive companion O star. It shows rather stable emissions at high energies over years and hence serves as an ideal laboratory to investigate the emission mechanism for such peculiar systems which emit prominent $\gamma$-rays. To this end, we take the orbital phase resolved energy spectrum as observed by \fermi over 10 years. We divide the orbit into four orbital phases, each with an orbital phase range of 0.25, centered at 0.00, 0.25, 0.50 and 0.75 respectively, where the phase 0.0 is the periastron and phase 0.5 is the apastron. The phases around 0.25 and 0.75 are symmetric and hence are supposed to have identical local acceleration environment. The spectral analysis shows that, the \fermi spectra are largely different from these two symmetric orbital phases: the emission from orbital phase 0.25 turns out to be significantly stronger than that from 0.75. This result does not fit a scenario that $\gamma$-rays are Doppler boosted emission from bow shock tails if LS 5039 has a shock configuration similar to PSR B1259-63, and indicates that the inverse Compton scatterings between the shock accelerated plasma and the stellar particle environment is the underline procedure. We also find that the previous report for a disappearance of the orbital modulation at 3--20 GeV is due to the similar spectral turn-over energies of the different orbital phases. The spectral properties of periastron and apastron regions are addressed in the context of the measurements in phase regions around 0.25 and 0.75.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state.

Large-scale outflows from starburst galaxies are multi-phase, multi-component fluids. Charge-exchange lines which originate from the interfacing surface between the neutral and ionised components are a useful diagnostic of the cold dense structures in the galactic outflow. From the charge-exchange lines observed in the nearby starburst galaxy M82, we conduct surface-to-volume analyses and deduce that the cold dense clumps in its galactic outflow have flattened shapes, resembling a hamburger or a pancake morphology rather than elongated shapes. The observed filamentary H$\alpha$ features are therefore not prime charge-exchange line emitters. They are stripped material torn from the slow moving dense clumps by the faster moving ionised fluid which are subsequently warmed and stretched into elongated shapes. Our findings are consistent with numerical simulations which have shown that cold dense clumps in galactic outflows can be compressed by ram pressure, and also progressively ablated and stripped before complete disintegration. We have shown that some clumps could survive their passage along a galactic outflow. These are advected into the circumgalactic environment, where their remnants would seed condensation of the circumgalactic medium to form new clumps. The infall of these new clumps back into the galaxy and their subsequent re-entrainment into the galactic outflow form a loop process of galactic material recycling.

In order to search for further observational evidence of cloud-cloud collisions in one of the promising candidates, L1188, we carried out observations of multiple molecular lines toward the intersection region of the two nearly orthogonal filamentary molecular clouds in L1188. Based on these observations, we find two parallel filamentary structures, both of which have at least two velocity components being connected with broad bridging features. We also found a spatially complementary distribution between the two molecular clouds, as well as enhanced $^{13}$CO emission and $^{12}$CO self-absorption toward their abutting regions. At the most blueshifted velocities, we unveil a 1~pc-long arc ubiquitously showing $^{12}$CO line wings. We discover two 22 GHz water masers, which are the first maser detections in L1188. An analysis of line ratios at a linear resolution of 0.2 pc suggests that L1188 is characterised by kinetic temperatures of 13--23~K and H$_{2}$ number densities of 10$^{3}$--10$^{3.6}$ cm$^{-3}$. On the basis of previous theoretical predictions and simulations, we suggest that these observational features can be naturally explained by the scenario of a cloud-cloud collision in L1188, although an additional contribution of stellar feedback from low-mass young stellar objects cannot be ruled out.

We propose a revision of the system developed by L'epine et al. (2007) for spectroscopic M subdwarf classification. Based on an analysis of subdwarf spectra and templates from Savcheva et al. (2014), we show thatthe CaH1 feature originally proposed by Gizis (1997) is important in selecting reliable cool subdwarf spectra. This index should be used in combination with the [TiO5, CaH2+CaH3] relation provided by L\'epine et al. (2007) to avoid misclassification results. In the new system, the dwarf-subdwarf separators are first derived from a sample of more than 80,000 M dwarfs and a "labeled" subdwarf subsample, these objects being all visually identified from their optical spectra. Based on these two samples, we re-fit the initial [TiO5, CaH1] relation, and propose a new [CaOH, CaH1] relation supplementing the [TiO5, CaH1] relation to reduce the impact of uncertainty in flux calibration on classification accuracy. In addition, we recalibrate the $\zeta_{TiO/CaH}$ parameter defined in L'epine et al. (2007) to enable its successful application to LAMOST spectra. Using this new system, we select candidates from LAMOST Data Release 4 and finally identify a set of 2791 new M subdwarf stars, covering the spectral sequence from type M0 to M7. This sample contains a large number of objects located at low Galactic latitudes, especially in the Galactic anti-center direction, expanding beyond previously published halo- and thick disk-dominated samples. Besides, we detect magnetic activity in 141 objects. We present a catalog for this M subdwarf sample, including radial velocities, spectral indices and errors, activity flags, with a compilation of external data (photometric and GAIA DR2 astrometric parameters). The catalog is provided on-line, and the spectra can be retrieved from the LAMOST Data Release web portal.

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

Existing analysis based on XMM-Newton/RGS spectra already shows that the G-ratio of the OVII He$\alpha$ triplet in the inner bulge of M31 is too high to be consistent with a pure optically thin thermal plasma in collisional ionization equilibrium (CIE). Different processes that may affect properties of diffuse hot plasma were proposed, such as resonance scattering (RS) and charge exchange (CX) with cold gas. To determine which physical process(es) may be responsible for this inconsistency, we present a systematic spectroscopic analysis based on 0.8 Ms XMM-Newton/RGS data, together with complementary Chandra/ACIS-S images. The combination of these data enables us to reveal multiple non-CIE spectroscopic diagnostics, including but not limited to the large G-ratios of He$\alpha$ triplets (OVII, NVI, and NeIX) and the high Lyman series line ratios (OVIII Ly$\beta$/Ly$\alpha$ and Ly$\gamma$/Ly$\alpha$, and NVII Ly$\beta$/Ly$\alpha$), which are not expected for a CIE plasma, and the high iron line ratios (FeXVIII 14.2 \AA/FeXVII~17 \AA and FeXVII~15 \AA/17 \AA), which suggest much higher temperatures than other line ratios, as well as their spatial variations. Neither CX nor RS explains all these spectroscopic diagnostics satisfactorily. Alternatively, we find that an active galactic nucleus (AGN) relic scenario provides a plausible explanation for virtually all the signatures. We estimate that an AGN was present at the center of M31 about half a million years ago and that the initial ionization parameter $\xi$ of the relic plasma is in the range of 3-4.

If captured by the gravitational field of stars or other compact objects, dark matter can self-annihilate and produce a potentially detectable particle flux. In the case of superheavy dark matter ($ m_{X} \gtrsim 10^{8} GeV $), a large number of scattering events with nuclei inside stars are necessary to slow down the dark matter particles below the escape velocity of the stars, at which point the Dark Matter (DM) particle becomes trapped, or captured. Using the recently developed analytical formalism for multiscatter capture, combined with the latest results on the constraints of dark-matter-baryon scattering cross-section, we calculate upper bounds on the capture rates for superheavy dark matter particles by the first (Pop. III) stars. Assuming that a non-zero fraction of the products of captured superheavy dark matter (SHDM) annihilations can be trapped and thermalized inside the star we find that this additional heat source could influence the evolutionary phase of Pop. III stars. Moreover, requiring that Pop. III stars shine with sub-Eddington luminosity, we find upper bounds on the masses of the Pop. III stars. This implies a DM dependent cutoff on the initial mass function (IMF) of Pop. III stars, thus opening up the intriguing possibility of constraining DM properties using the IMF of extremely metal-poor stars.

The hard X-ray transient source Swift J1845.7-0037 was discovered in 2012 by Swift/BAT. However, at that time no dedicated observations of the source were performed. On Oct 2019 the source became active again, and X-ray pulsations with a period of ~199s were detected with Swift/XRT. This triggered follow-up observations with NuSTAR. Here we report on the timing and spectral analysis of the source properties using NuSTAR and Swift/XRT. The main goal was to confirm pulsations and search for possible cyclotron lines in the broadband spectrum of the source to probe its magnetic field. Despite highly significant pulsations with period of 207.379(2) were detected, no evidence for a cyclotron line was found in the spectrum of the source. We therefore discuss the strength of the magnetic field based on the source flux and the detection of the transition to the "cold-disc" accretion regime during the 2012 outburst. Our conclusion is that, most likely, the source is a highly magnetized neutron star with B 1e13G at a large distance of d~10 kpc. The latter one consistent with the non-detection of a cyclotron line in the NuSTAR energy band.

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.

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

We present three simultaneous/quasi-simultaneous NuSTAR and Swift datasets of the black hole GRS 1716-249 in its hard intermediate state. The accretion disk in this state may have reached the innermost stable circular orbit, and the NuSTAR spectra show a broad relativistic iron line and a strong Compton hump. To measure the black hole spin, we construct a joint model consisting of a relativistic disk model kerrbb and a reflection model relxill, to fit the continuum and the reflection components, respectively. By applying this model to each dataset independently, a consistent result is obtained on the black hole spin and the disk inclination. The black hole spin is a* >~ 0.92, and the inclination angle (i) is around 40-50 degree, based on the measurements of all datasets. In the third dataset, a high black hole mass (M_BH) is strongly disfavored by the spectral fits. By unfreezing the black hole mass, we find a*>0.92, i=49.9^{+1.0}_{-1.3} degree and M_BH<8.0 Msun, at a 90% confidence level. Considering the lower limit derived from a previous optical constraint, M_BH is in a range of 4.9-8.0 Msun.

We present three simultaneous/quasi-simultaneous NuSTAR and Swift datasets of the black hole GRS 1716-249 in its hard intermediate state. The accretion disk in this state may have reached the innermost stable circular orbit, and the NuSTAR spectra show a broad relativistic iron line and a strong Compton hump. To measure the black hole spin, we construct a joint model consisting of a relativistic disk model kerrbb and a reflection model relxill, to fit the continuum and the reflection components, respectively. By applying this model to each dataset independently, a consistent result is obtained on the black hole spin and the disk inclination. The black hole spin is a* >~ 0.92, and the inclination angle (i) is around 40-50 degree, based on the measurements of all datasets. In the third dataset, a high black hole mass (M_BH) is strongly disfavored by the spectral fits. By unfreezing the black hole mass, we find a*>0.92, i=49.9^{+1.0}_{-1.3} degree and M_BH<8.0 Msun, at a 90% confidence level. Considering the lower limit derived from a previous optical constraint, M_BH is in a range of 4.9-8.0 Msun.

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.

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.

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.

The low energy (LE) X-ray telescope is one of the three main instruments of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). It is equipped with Swept Charge Device (SCD) sensor arrays with a total geometrical area of 384 cm2 and an energy band from 0.7 keV to 13 keV. In order to evaluate the particle induced X-ray background and the cosmic X-ray background simultaneously, LE adopts collimators to define four types of Field Of Views (FOVs). LE is constituted of three detector boxes (LEDs) and an electric control box (LEB) and achieves a good energy resolution of 140 eV at 5.9 keV, an excellent time resolution of 0.98 ms, as well as an extremely low pileup (<1% at 18000 cts/s). Detailed performance tests and calibration on the ground have been performed, including energy-channel relation, energy response, detection efficiency and time response.

The low energy (LE) X-ray telescope is one of the three main instruments of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). It is equipped with Swept Charge Device (SCD) sensor arrays with a total geometrical area of 384 cm2 and an energy band from 0.7 keV to 13 keV. In order to evaluate the particle induced X-ray background and the cosmic X-ray background simultaneously, LE adopts collimators to define four types of Field Of Views (FOVs). LE is constituted of three detector boxes (LEDs) and an electric control box (LEB) and achieves a good energy resolution of 140 eV at 5.9 keV, an excellent time resolution of 0.98 ms, as well as an extremely low pileup (<1% at 18000 cts/s). Detailed performance tests and calibration on the ground have been performed, including energy-channel relation, energy response, detection efficiency and time response.

The Hard X-ray Modulation Telescope ($Insight$-HXMT) was successfully launched on June 15th, 2017. It performs broad band X-ray scan survey of the Galactic Plane to detect new black holes and other objects in active states. It also observes X-ray binaries to study their X-ray variabilities. Here we will introduce the Science Operations of $Insight$-HXMT, which is responsible for collecting and evaluating observation proposals, scheduling observations, and monitoring the working status of the payloads.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

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.

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.

IRAS F11119+3257 is a quasar-dominated Ultra-Luminous InfraRed Galaxy, with a partially obscured narrow-line seyfert 1 nucleus. In this paper, we present the NIR spectroscopy of F11119+3257, in which we find unusual Paschen emission lines, and metastable He I* $\lambda$10830 absorption associated with the previously reported atomic sodium and molecular OH mini-BAL (Broad Absorption Line) outflow. Photo-ionization diagnosis confirms previous findings that the outflows are at kilo-parsec scales. Such large-scale outflows should produce emission lines. We indeed find that high-ionization emission lines ([O III], [Ne III], and [Ne V]) are dominated by blueshifted components at similar speeds to the mini-BALs. The blueshifted components are also detected in some low-ionization emission lines, such as [O II] $\lambda$3727 and some Balmer lines (H$\alpha$, H$\beta$, and H$\gamma$), even though their cores are dominated by narrow ($FWHM_{\rm NEL} = 570\pm40$km s$^{-1}$) or broad components at the systemic redshift of $z=0.18966\pm0.00006$. The mass flow rate (230-730$~M_\odot \rm yr^{-1}$) and the kinetic luminosity ($\dot{E}_k \sim 10^{43.6-44.8} $erg s$^{-1}$) are then inferred jointly from the blueshifted emission and absorption lines. In the NIR spectrum of F11119+3257, we also find that the Paschen emission lines are unique, in which a very narrow ($FWHM=260\pm20~$km s$^{-1}$) component is shown in only Pa$\alpha$. This narrow component most probably comes from heavily obscured star formation. Based on the Pa$\alpha$ and Pa$\beta$ emissions, we obtain an extinction at the $H$ band, $A_H~>~2.1$ (or a reddenning of $E_{B-V}~>~$3.7), and a star formation rate of $SFR~>~130\rm M_\odot yr^{-1}$ that resembles the estimates inferred from the FIR emissions ($SFR_{\rm FIR} = 190\pm90$ $M_\odot$ yr$^{-1}$).

We present spectropolarimetric and spectrophotometric observations of the peculiar active galactic nucleus (AGN) SDSS J120300.19+162443.7 (hereafter J1203+1624) at $z=0.1656$. Its optical total flux spectra clearly show broad emission lines (BELs) in H$\alpha$ and H$\beta$. After removal of narrow emission lines (NELs), the full width at half maximum of the lorentzian BEL is $FWHM_{\rm BEL} \approx 1,270$ km s$^{-1}$, fulfilling the conventional definition of a narrow line Seyfert 1 (NLS1) galaxy. However, its NELs are unprecedentedly strong when compared to type 1AGNs. This, together with its large MIR excess ($g-W_4 = 13.172$ mag), implies that the observer and the narrow emission line region (NELR) might see a different ionization continuum. Our optical spectropolarimetry confirms its type 2 nature by detecting a polarized blue continuum and Balmer BELs ($FWHM_{\rm Polarized~BEL} \approx 1,183$ km s$^{-1}$), with a high polarization degree of $> 20$% in the blue wing. The mass and Eddington rate of the central black hole is estimated based on both transmitted and scattered AGN radiation, which is $M_{\bullet} < 2.9 \times 10^7 M_\odot$ and $L_{\rm bol}/L_{\rm Edd} > 1.5$. Severe extinction of the AGN emission also enables a clear view of the compact host galaxy. Discovery of J1203+1624 suggests that NLS1 follows AGN unfication scheme, studying of its analogs could blaze a new trail for exploring the connection between black hole growth and star formation in the host galaxy. Interesting features of J1203+1624 like the peculiar NELs, inflowing scattering clouds within sublimation radius are worth detailed follow-ups in the future.

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning 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 introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning 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 introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

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.

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.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) defined by the collimators, and CsI(Na) is used as an active shield detector to NaI(Tl) by pulse shape discrimination. CsI(Na) is also used as an omnidirectional gamma-ray monitor. The HE collimators have a diverse FOV: 1.1{\deg}x 5.7{\deg} (15 units), 5.7{\deg}x 5.7{\deg} (2 units) and blocked (1 unit), thus the combined FOV of HE is about 5.7{\deg}x 5.7{\deg}. Each HE detector has a diameter of 190 mm, resulting in the total geometrical area of about 5100 cm_2. The energy resolution is ~15% at 60 keV. The timing accuracy is better than 10 {\mu}s and dead-time for each detector is less than 10 {\mu}s. HE is devoted to observe the spectra and temporal variability of X-ray sources in the 20-250 keV band either by pointing observations for known sources or scanning observations to unveil new sources, and to monitor the gamma-ray sky in 0.2-3 MeV. This paper presents the design and performance of the HE instruments. Results of the on-ground calibration experiments are also reported.

The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) defined by the collimators, and CsI(Na) is used as an active shield detector to NaI(Tl) by pulse shape discrimination. CsI(Na) is also used as an omnidirectional gamma-ray monitor. The HE collimators have a diverse FOV: 1.1{\deg}x 5.7{\deg} (15 units), 5.7{\deg}x 5.7{\deg} (2 units) and blocked (1 unit), thus the combined FOV of HE is about 5.7{\deg}x 5.7{\deg}. Each HE detector has a diameter of 190 mm, resulting in the total geometrical area of about 5100 cm_2. The energy resolution is ~15% at 60 keV. The timing accuracy is better than 10 {\mu}s and dead-time for each detector is less than 10 {\mu}s. HE is devoted to observe the spectra and temporal variability of X-ray sources in the 20-250 keV band either by pointing observations for known sources or scanning observations to unveil new sources, and to monitor the gamma-ray sky in 0.2-3 MeV. This paper presents the design and performance of the HE instruments. Results of the on-ground calibration experiments are also reported.

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.

We present the observational results from a detailed timing analysis of the black hole candidate Swift J1658.2-4242 during its 2018 outburst with the observations of Hard X-ray Modulation Telescope (Insight-HXMT), Neutron Star Interior Composition Explorer (NICER) and AstroSat in 0.1-250keV. The evolution of intensity, hardness and integrated fractional root mean square (rms) observed by Insight-HXMT and NICER are presented in this paper. Type-C quasi-periodic oscillations (QPOs) observed by NICER (0.8-3.5Hz) and Insight-HXMT (1-1.6Hz) are also reported in this work. The features of the QPOs are analysed with an energy range of 0.5-50keV. The relations between QPO frequency and other characteristics such as intensity, hardness and QPO rms are carefully studied. The timing and spectral properties indicate that Swift J1658.2-4242 is a black hole binary system. Besides, the rms spectra of the source calculated from the simultaneous observation of Insight-HXMT, NICER and AstroSat support the Lense-Thirring origin of the QPOs. The relation between QPO phase lag and the centroid frequency of Swift J1658.2-4242 reveals a near zero constant when < 4Hz and a soft phase lag at 6.68 Hz. This independence follows the same trend as the high inclination galactic black hole binaries such as MAXI J1659-152.

We present the observational results from a detailed timing analysis of the black hole candidate Swift J1658.2-4242 during its 2018 outburst with the observations of Hard X-ray Modulation Telescope (Insight-HXMT), Neutron Star Interior Composition Explorer (NICER) and AstroSat in 0.1-250keV. The evolution of intensity, hardness and integrated fractional root mean square (rms) observed by Insight-HXMT and NICER are presented in this paper. Type-C quasi-periodic oscillations (QPOs) observed by NICER (0.8-3.5Hz) and Insight-HXMT (1-1.6Hz) are also reported in this work. The features of the QPOs are analysed with an energy range of 0.5-50keV. The relations between QPO frequency and other characteristics such as intensity, hardness and QPO rms are carefully studied. The timing and spectral properties indicate that Swift J1658.2-4242 is a black hole binary system. Besides, the rms spectra of the source calculated from the simultaneous observation of Insight-HXMT, NICER and AstroSat support the Lense-Thirring origin of the QPOs. The relation between QPO phase lag and the centroid frequency of Swift J1658.2-4242 reveals a near zero constant when < 4Hz and a soft phase lag at 6.68 Hz. This independence follows the same trend as the high inclination galactic black hole binaries such as MAXI J1659-152.

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 long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

In this paper, we study the gravitational waveforms, polarizations and radiation powers of the first relativistic triple systems PSR J0337 + 1715, observed in 2014, by using the post-Newtonian approximations to their lowest order. Although they cannot be observed either by current or next generation of the detectors, they do provide useful information to test different theories of gravity. In particular, we carry out the studies in three different theories, general relativity (GR), Brans-Dicke (BD) gravity, and Einstein-aether ($\ae$) theory. The tensor modes $h_{+}$ and $h_{\times}$ exist in all three theories and have almost equal amplitudes. Their frequencies are all peaked at two locations, $ f^{(+, \times)}_1 \simeq 0.0686656 \mu$Hz and $ f^{(+, \times)}_2 \simeq 14.2138 \mu$Hz, which are about twice of the outer and inner orbital frequencies of the triple system, as predicted in GR. In $\ae$-theory, all the six polarization modes are different from zero, but the breathing $h_b$ and longitudinal $h_L$ modes are not independent and also peaked at two frequencies, but at the frequencies, $ f^{(b, L, \ae)}_{1} \simeq 0.0457771 \mu$Hz and $f^{(b, L, \ae)}_{ 2} \simeq 7.09545 \mu$Hz. A similar phenomenon is also observed in BD gravity, in which only the three modes $h_{+},\; h_{\times}$ and $h_{b}$ exit. We also study the radiation powers, and find that the quadrupole emission in each of the three theories has almost the same amplitude, but the dipole emission can be as big as the quadrupole emission in $\ae$-theory. This can provide a very promising window to obtain severe constraints on $\ae$-theory by the multi-band gravitational wave astronomy.

We extract the multiple-epoch Balmer-line profiles of the heavily obscured quasar 3C 445 from the spectral curves in the literature, and analyze the emission-line profiles of the H$\alpha$ and H$\beta$ lines and the profile variability in the H$\alpha$ line in the large time interval of more than three decades. The profile comparison between the H$\alpha$ and H$\beta$ lines shows that both Balmer lines share the profile with the same form, while the blue system of the H$\beta$ line is seriously weaker than that of the H$\alpha$ line. Moreover, the blue system of the H$\alpha$ line suddenly disappeared completely and then did not appear again, however the other two components did not exhibit significant variation in the velocity or the amplitude. These findings suggest that the blue system of 3C 445, as with SDSS J153636.22+044127.0 and its analogs, is probable the result of the shock-heated outflowing gases. The observation angle of almost edge-on which the previous studies suggested can easily produce the high-speed and high-temperature shock in the collision between the massive outflow and the inner surface of the dusty torus.

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 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.

The interstellar medium (ISM) is a complex non-linear system governed by gravity and magneto-hydrodynamics, as well as radiative, thermodynamical, and chemical processes. Our understanding of it mostly progresses through observations and numerical simulations, and a quantitative comparison between these two approaches requires a generic and comprehensive statistical description. The goal of this paper is to build such a description, with the purpose to permit an efficient comparison independent of any specific prior or model. We start from the Wavelet Scattering Transform (WST), a low-variance statistical description of non-Gaussian processes, developed in data science, that encodes long-range interactions through a hierarchical multiscale approach based on the Wavelet transform. We perform a reduction of the WST through a fit of its angular dependencies, allowing to gather most of the information it contains into a few components whose physical meanings are identified, and that describe, e.g., isotropic and anisotropic behaviours. The result of this paper is the Reduced Wavelet Scattering Transform (RWST), a statistical description with a small number of coefficients that characterizes complex structures arising from non-linear phenomena, free from any specific prior. The RWST coefficients encode moments of order up to four, have reduced variances, and quantify the couplings between scales. To show the efficiency and generality of this description, we apply it successfully to three kinds of processes: fractional Brownian motions, MHD simulations, and Herschel observations in a molecular cloud. With fewer than 100 coefficients when probing 6 scales and 8 angles on 256*256 maps, we were able with the RWST to perform quantitative comparisons, to infer relevant physical properties, and to produce realistic synthetic fields.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

Using the archive data from the Rossi X-ray Timing Explorer ({\sl RXTE}), we have studied the evolution of the X-ray profile of the Crab pulsar in a time span of 11 years. The X-ray profiles, as characterized by a few parameters, changed slightly but significantly in these years: the separation of the two peaks increased with a rate $0.88\pm0.20\,\textordmasculine$\,per century, the flux ratio of the second pulse to the first pulse decreased with $(3.64\pm0.86)\times10^{-2}$\,per century, and the pulse widths of the two pulses descended with $1.44\pm0.15\,\textordmasculine$\, and $1.09\pm0.73\,\textordmasculine$\,per century, respectively. The evolutionary trends of the above parameters are similar to the radio results, but the values are different. We briefly discussed the constraints of these X-ray properties on the geometry of the emission region of this pulsar.

It is believed that an isolated pulsar loses its rotational energy mainly through a relativistic wind consisting of electrons, positrons and possibly Poynting flux\cite{Pacini1973,Rees1974,Kennel1984}. As it expands, this wind may eventually be terminated by a shock, where particles can be accelerated to energies of X-ray synchrotron emission, and a pulsar wind nebula (PWN) is usually detectable surrounding a young energetic pulsar\cite{Pacini1973,Rees1974,Kennel1984}. However, the nature and/or energetics of these physical processes remain very uncertain, largely because they typically cannot be studied in a time-resolved fashion. Here we show that the X-ray PWN around the young pulsar PSR B0540--69 brightens gradually up to 32$\pm8\%$ over the mean previous flux, after a sudden spin-down rate ($\dot{\nu}$) transition (SRT) by $\sim36\%$\ in December 2011, which has very different properties from a traditional pulsar glitch\cite{Marshall2015}. No evidence is seen for any change in the pulsed X-ray emission. We conclude that the SRT results from a sudden change in the pulsar magnetosphere that increases the pulsar wind power and hence the PWN X-ray emission. The X-ray light curve of the PWN suggests a mean life time of the particles of $397\pm374$\,days, corresponding to a magnetic field strength of $0.78_{-0.28}^{+4.50}$\,mG in the PWN.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

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.

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.

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

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 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.

POLAR is a compact space-borne detector initially designed to measure the polarization of hard X-rays emitted from Gamma-Ray Bursts in the energy range 50-500keV. This instrument was launched successfully onboard the Chinese space laboratory Tiangong-2 (TG-2) on 2016 September 15. After being switched on a few days later, tens of gigabytes of raw detection data were produced in-orbit by POLAR and transferred to the ground every day. Before the launch date, a full pipeline and related software were designed and developed for the purpose of quickly pre-processing all the raw data from POLAR, which include both science data and engineering data, then to generate the high level scientific data products that are suitable for later science analysis. This pipeline has been successfully applied for use by the POLAR Science Data Center in the Institute of High Energy Physics (IHEP) after POLAR was launched and switched on. A detailed introduction to the pipeline and some of the core relevant algorithms are presented in this paper.

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.

The Parkes Pulsar Timing Array (PPTA) project monitors two dozen millisecond pulsars (MSPs) in order to undertake a variety of fundamental physics experiments using the Parkes 64m radio telescope. Since June 2017 we have been undertaking commensal searches for fast radio bursts (FRBs) during the MSP observations. Here, we report the discovery of four FRBs (171209, 180309, 180311 and 180714). The detected events include an FRB with the highest signal-to-noise ratio ever detected at the Parkes observatory, which exhibits unusual spectral properties. All four FRBs are highly polarized. We discuss the future of commensal searches for FRBs at Parkes.

We present a new method to investigate the effective magnetic field decay of isolated neutron stars, from the analysis of the long-term timing data of a large sample of radio pulsars \citep{2010MNRAS.402.1027H}. There are some differences between the distributions of frequency's second derivatives of the pulsar spins with different effective field decay timescales. Kolmogorov-Smirnov tests are performed to reexamine the consistency of distributions of the simulated and reported data for a series values of decay timescales. \textbf{We show that the timescale of the effective field decay exceeds $\sim 5~\rm {Myr}$ for pulsars with spin-down age $\tau_{\rm C}< 10^{7}~{\rm yr}$ or $\sim 100~\rm {Myr}$ for pulsars with $10^{7}<\tau_{\rm C}< 10^{9}~{\rm yr}$ in the sample. The result does not depend on any specific theories of the field evolution, the inclination decay or the variation in the moment of inertia. It is also found that the extent of the closed line region of the magnetic field is close to the light cylinder $r_{\rm lc}$, i.e., the corotating radius $r_{\rm c}\approx r_{\rm lc}$ is a good approximation for the observed pulsar population.

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 studied the long-term evolution of the spectral-temporal correlated properties of the black hole candidate Swift J1753.5-0127 from the onset of its outburst until 2011 with the Rossi X-ray Timing Explorer (RXTE). The source stayed most of its lifetime during hard state, with occasionally transitioned to the hard intermediate state. Similar to typical black hole transients, Swift J1753.5-0127 traces a clear hard line in absolute rms - intensity diagram during the low hard state, with expected highest absolute rms, while shows a clear turn during the hard intermediate state, accompanied by lower absolute rms. Different from Cyg X-1, we found that frequency-dependent time lag increased significantly in the 0.02 -- 3.2 Hz band during state transition in this source. The X-ray time lags in 0.02 -- 3.2 Hz can therefore be used as indicators of state transition in this source. Type-C QPO frequency is positively related with its fractional rms and X-ray photon index, suggesting a moving inwards disk/corona scenario. We discussed the physical interpretation of our results in our paper.

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.

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.

Large High Altitude Air Shower Observatory(LHAASO) is a composite cosmic ray observatory consisting of three detector arrays: kilometer square array (KM2A) which includes the electromagnetic detector array and muon detector array, water Cherenkov detector array (WCDA) and wide field of view Cherenkov telescope array (WFCTA). One of the main scientific objectives of LHAASO is to precisely measure the cosmic rays energy spectrum of individual components from 1014 eV to 1018 eV. The hybrid observation will be employed by LHAASO experiment, in which the lateral and longitudinal distributions of the extensive air shower can be observed simultaneously. Thus many kinds of parameters can be used for primary nuclei identification. In this paper, high purity cosmic ray simulation samples of light nuclei component are obtained through Multi-Variable Analysis. The apertures of 1/4 LHAASO array for pure proton and mixed proton and helium (H&He) samples are 900 m2Sr and 1800 m2Sr respectively. A prospect of proton and H&He spectra from 100 TeV to 4 PeV is discussed.

The binary emission-line system, particularly the binary broad-line emission system, is considered the most effective indicator of the black hole binary. A plausible example of such a system, SDSS J153636.22+044127.0, was reported as the first known object with two hydrogen Balmer broad-line systems, which are interpreted to be the result of broad-line regions around a pair of black holes (Boroson \& Lauer 2009). Here, we show the follow-up optical and near-infrared spectral observations of SDSS J153636.22+044127.0 and its analogues. In these objects, the broad hydrogen Balmer and Paschen, He I and Mg II lines share the same peculiar emission-line profile (including a blue system, a red system and a double-peaked disk-line component); however, the invariance in the large time interval, the absence of the blue system in He I $\lambda$10830 profile and the abnormally strong emission of the hydrogen Pa$\beta$ blue system oppose the binary proposal. We suggest that these unique broad lines arise from the AGN emission-line region and the shock-heated outflowing gases rather than a binary system of two active black holes.

Gravitational wave emissions (GWEs) of pulsars could not only make them promising targets for continuous gravitational wave searches but also leave imprints in their timing data. We interpret the measured braking index of PSR J1640-4631 with a model involving both the GWE and dipole magnetic field decay. Combining the timing data of PSR J1640-4631 and the theory of magnetic field decay, we propose a new approach of constraining the number of precession cycles, $\xi$, which is highly uncertain currently but can be tightly related to the interior physics of a neutron star and its GWE. We suggest that future observation of the tilt angle $\chi$ of PSR J1640-4631 would not merely help to constrain $\xi$ but also possibly provide information about the internal magnetic field configuration of this pulsar. We find that $\xi$ would be larger than previous estimates unless a tiny angle $\chi\lesssim5^\circ$ is observed. Furthermore, a measured angle $\chi\gtrsim12^\circ$ would indicate $\xi\gtrsim 10^6$, which is at least ten times larger than that suggested previously.