The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

We present a Chandra archival study of optically selected AGN pairs at a median redshift $\bar{z} \sim$ 0.1. Out of 1286 AGN pairs (with projected separations $r_p < 100$ kpc and velocity offsets $\Delta v < 600$ km s$^{-1}$) optically identified from the Sloan Digital Sky Survey Seventh Data Release, we find 67 systems with archival Chandra observations, which represent the largest sample of optically selected AGN pairs studied in the X-ray. Among the 67 AGN pairs, 21 systems have both nuclei detected in the X-rays, 36 have one nucleus detected in the X-rays, whereas 10 have no X-ray detection. The X-ray detection rate, 78/134 = 58% ($\pm 7$% 1$\sigma$ Poisson errors), is significantly higher than that (23/134 = $17\% \pm 4\%$) of a comparison sample of star-forming galaxy pairs, lending support to the optical AGN classification. In the conservative case where X-ray contamination from star formation is removed, the X-ray detection rate becomes $27\% \pm 4\%$, consistent with predictions from latest galaxy merger simulations. The 2-10 keV X-ray luminosity $L_{2-10~{\rm keV}}$ increases with decreasing projected separation in AGN pairs for $r_p \gtrsim 15$ kpc, suggesting an enhancement of black hole accretion even in early-stage mergers. On the other hand, $L_{2-10~{\rm keV}}$ appears to decrease with decreasing projected separation at $r_p \lesssim 15$ kpc, which is contradictory to predictions from merger simulations. The apparent decrease in $L_{2-10~{\rm keV}}$ of AGN pairs at $r_p \lesssim 15$ kpc may be caused by: (i) enhanced absorbing columns from merger-induced gas inflows, (ii) feedback effects from early-stage mergers, and/or (iii) small number statistics. Future X-ray studies with larger samples are needed to put our results on firmer statistical ground.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We present here the analysis of giant micropulses from the Vela pulsar. A total of 4187 giant micropulses with peak flux density $>$2.5 Jy were detected during almost 4 hours of observations carried out with the Yunnan 40-m radio telescope at 6800 MHz. Nine of the giant micropulses arrived approximately 3 to 4 ms earlier than the peak of average pulse profile, longer than that at lower frequencies. The remaining giant micropulses were clustered into three distributions which correspond to three main emission regions, including four occurring on the trailing edge of averaged profile.

Convolutional Neutral Networks have been successfully applied in searching for strong lensing systems, leading to discoveries of new candidates from large surveys. On the other hand, systematic investigations about their robustness are still lacking. In this paper, we first construct a neutral network, and apply it to $r$-band images of Luminous Red Galaxies (LRGs) of the Kilo Degree Survey (KiDS) Data Release 3 to search for strong lensing systems. We build two sets of training samples, one fully from simulations, and the other one using the LRG stamps from KiDS observations as the foreground lens images. With the former training sample, we find 48 high probability candidates after human-inspection, and among them, 27 are newly identified. Using the latter training set, about 67\% of the aforementioned 48 candidates are also found, and there are 11 more new strong lensing candidates identified. We then carry out tests on the robustness of the network performance with respect to the variation of PSF. With the testing samples constructed using PSF in the range of 0.4 to 2 times of the median PSF of the training sample, we find that our network performs rather stable, and the degradation is small. We also investigate how the volume of the training set can affect our network performance by varying it from 0.1 millions to 0.8 millions. The output results are rather stable showing that within the considered range, our network performance is not very sensitive to the volume size.

Recent observations have shown a remarkable diversity of the observational behaviors and explosion mechanisms in thermonuclear supernovae (SNe). An emerging class of peculiar type Ia SNe, called type Iax, show distinct photometric and spectroscopic behaviors. Their origin remains highly controversial, but pure turbulent deflagration of white dwarfs (WDs) has been regarded as the leading formation theory. The large population of type Iax indicates the existence of unidentified Galactic type Iax supernova remnants (SNRs). We report evidence that SNR Sgr A East in the Galactic center resulted from a pure turbulent deflagration of a Chandrasekhar-mass carbon-oxygen WD, an explosion mechanism used for type Iax SNe. Our X-ray spectroscopic study of Sgr A East using 3 Ms Chandra data shows a low ratio of intermediate-mass elements to Fe, and large Mn/Fe and Ni/Fe ratios. This abundance pattern does not accord with the core-collapse or normal type Ia models. Sgr A East is thus the first Galactic SNR for which a likely type Iax origin has been proposed, and the nearest target to study this peculiar class. We compared Sgr A East with the Fe-rich SNRs 3C 397 and W49B, which also have high Mn and Cr abundances and were claimed to result from deflagration-to-detonation explosions of Chandrasekhar-mass WDs (although with disputes). Our study shows that they have distinct abundance patterns. The X-ray spectroscopic studies of thermonuclear SNRs provide observational evidence for the theories that there are diverse explosion channels and various metal outputs for Chandrasekhar-mass WDs.

Astronomical observation data require long-term preservation, and the rapid accumulation of observation data makes it necessary to consider the cost of long-term archive storage. In addition to low-speed disk-based online storage, optical disk or tape-based offline storage can be used to save costs. However, for astronomical research that requires historical data (particularly time-domain astronomy), the performance and energy consumption of data-accessing techniques cause problems because the requested data (which are organized according to observation time) may be located across multiple storage devices. In this study, we design and develop a tool referred to as AstroLayout to redistribute the observation data using spatial aggregation. The core algorithm uses graph partitioning to generate an optimized data placement according to the original observation data statistics and the target storage system. For the given observation data, AstroLayout can copy the long-term archive in the target storage system in accordance with this placement. An efficiency evaluation shows that AstroLayout can reduce the number of devices activated when responding to data-access requests in time-domain astronomy research. In addition to improving the performance of data-accessing techniques, AstroLayout can also reduce the storage systems power consumption. For enhanced adaptability, it supports storage systems of any media, including optical disks, tapes, and hard disks.

We propose a light-weight deep convolutional neural network (CNN) to estimate the cosmological parameters from simulated 3-dimensional dark matter distributions with high accuracy. The training set is based on 465 realizations of a cubic box with a side length of $256\ h^{-1}\ \rm Mpc$, sampled with $128^3$ particles interpolated over a cubic grid of $128^3$ voxels. These volumes have cosmological parameters varying within the flat $\Lambda$CDM parameter space of $0.16 \leq \Omega_m \leq 0.46$ and $2.0 \leq 10^9 A_s \leq 2.3$. The neural network takes as an input cubes with $32^3$ voxels and has three convolution layers, three dense layers, together with some batch normalization and pooling layers. In the final predictions from the network we find a $2.5\%$ bias on the primordial amplitude $\sigma_8$ that can not easily be resolved by continued training. We correct this bias to obtain unprecedented accuracy in the cosmological parameter estimation with statistical uncertainties of $\delta \Omega_m$=0.0015 and $\delta \sigma_8$=0.0029, which are several times better than the results of previous CNN works. Compared with a 2-point analysis method using clustering region of 0-130 and 10-130 $h^{-1}$ Mpc, the CNN constraints are several times and an order of magnitude more precise, respectively. Finally, we conduct preliminary checks of the error-tolerance abilities of the neural network, and find that it exhibits robustness against smoothing, masking, random noise, global variation, rotation, reflection, and simulation resolution. Those effects are well understood in typical clustering analysis, but had not been tested before for the CNN approach. Our work shows that CNN can be more promising than people expected in deriving tight cosmological constraints from the cosmic large scale structure.

The tomographic AP method is so far the best method in separating the Alcock-Paczynski (AP) signal from the redshift space distortion (RSD) effects and deriving powerful constraints on cosmological parameters using the $\lesssim40h^{-1}\ \rm Mpc$ clustering region. To guarantee that the method can be easily applied to the future large scale structure (LSS) surveys, we study the possibility of estimating the systematics of the method using fast simulation method. The major contribution of the systematics comes from the non-zero redshift evolution of the RSD effects, which is quantified by $\hat\xi_{\Delta s}(\mu,z)$ in our analysis, and estimated using the BigMultidark exact N-body simulation and approximate COLA simulation samples. We find about 5\%/10\% evolution when comparing the $\hat\xi_{\Delta s}(\mu,z)$ measured as $z=0.5$/$z=1$ to the measurements at $z=0$. We checked the inaccuracy in the 2pCFs computed using COLA, and find it 5-10 times smaller than the intrinsic systematics of the tomographic AP method, indicating that using COLA to estimate the systematics is good enough. Finally, we test the effect of halo bias, and find $\lesssim$1.5\% change in $\hat\xi_{\Delta s}$ when varying the halo mass within the range of $2\times 10^{12}$ to $10^{14}$ $M_{\odot}$. We will perform more studies to achieve an accurate and efficient estimation of the systematics in redshift range of $z=0-1.5$.

We report new HI observations of four z$\sim$0.05 star-forming galaxies undertaken during the commissioning phase of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). FAST is the largest single-dish telescope with a 500 meter aperture and a 19-Beam receiver. Exploiting the unprecedented sensitivity provided by FAST, we aim to study the atomic gas, via the HI 21cm emission line, in low-$z$ star-forming galaxies taken from the Valpara\'iso ALMA/APEX Line Emission Survey (VALES) project. Together with previous ALMA CO($J=1-0$) observations, the HI data provides crucial information to measure the gas mass and dynamics. As a pilot HI survey, we targeted four local star-forming galaxies at $z\sim0.05$. In particular, one of them has already been detected in HI by the Arecibo Legacy Fast ALFA survey (ALFALFA), allowing a careful comparison. We use an ON-OFF observing approach that allowed us to reach an rms of 0.7mJy/beam at a 1.7km/s velocity resolution within only 20 minutes ON-target integration time. We demonstrate the great capabilities of the FAST 19-beam receiver for pushing the detectability of the HI emission line of extra-galactic sources. The HI emission line detected by FAST shows good consistency with the previous ALFALFA results. Our observations are put in context with previous multi-wavelength data to reveal the physical properties of these low-$z$ galaxies. We find that the CO($J=1-0$) and HI emission line profiles are similar. The dynamical mass estimated from the HI data is an order of magnitude higher than the baryon mass and the dynamical mass derived from the CO observations, implying that the mass probed by dynamics of HI is dominated by the dark matter halo. In one case, a target shows an excess of CO($J=1-0$) in the line centre, which can be explained by an enhanced CO($J=1-0$) emission induced by a nuclear starburst showing high velocity dispersion.

We present 1.3 mm ALMA dust polarization observations at a resolution of $\sim$0.02 pc of three massive molecular clumps, MM1, MM4, and MM9, in the infrared dark cloud G28.34+0.06. With the sensitive and high-resolution continuum data, MM1 is resolved into a cluster of condensations. The magnetic field structure in each clump is revealed by the polarized emission. We found a trend of decreasing polarized emission fraction with increasing Stokes $I$ intensities in MM1 and MM4. Using the angular dispersion function method (a modified Davis-Chandrasekhar-Fermi method), the plane-of-sky magnetic field strength in two massive dense cores, MM1-Core1 and MM4-Core4, are estimated to be $\sim$1.6 mG and $\sim$0.32 mG, respectively. \textbf{The ordered magnetic energy is found to be smaller than the turbulent energy in the two cores, while the total magnetic energy is found to be comparable to the turbulent energy.} The total virial parameters in MM1-Core1 and MM4-Core4 are calculated to be $\sim$0.76 and $\sim$0.37, respectively, suggesting that massive star formation does not start in equilibrium. Using the polarization-intensity gradient-local gravity method, we found that the local gravity is closely aligned with intensity gradient in the three clumps, and the magnetic field tends to be aligned with the local gravity in MM1 and MM4 except for regions near the emission peak, which suggests that the gravity plays a dominant role in regulating the gas collapse. Half of the outflows in MM4 and MM9 are found to be aligned within 10$^{\circ}$ of the condensation-scale ($<$0.05 pc) magnetic field, indicating that the magnetic field could play an important role from condensation to disk scale in the early stage of massive star formation. We also found that the fragmentation in MM1-Core1 cannot be solely explained by thermal Jeans fragmentation or turbulent Jeans fragmentation.

Fast radio bursts (FRBs) are a mysterious astrophysical phenomenon of bright pulses emitted at radio frequencies, and it is expected for them to be frequently detected in the future. The dispersion measures of FRBs are related to cosmological parameters, and thus FRBs have a potential to be developed into a new cosmological probe if their data can be largely accumulated in the future. In this work, we study the capability of future FRB data for improving cosmological parameter estimation in two dynamical dark energy models. We find that the simulated FRB data can break the parameter degeneracies inherent in the current cosmic microwave background (CMB) data. Therefore, the combination of the CMB and FRB data can significantly improve the constraints on the Hubble constant and dark energy parameters, compared to those using CMB or FRB alone. If ten thousand events of FRBs with known redshifts are detected in the future, they would behave better than the baryon acoustic oscillation (BAO) data in breaking the parameter degeneracies inherent in the CMB data. We also find that the combination of FRB and gravitational wave (GW) standard siren data provides an independent low-redshift probe to verify the results from CMB and BAO data. For the data combination of CMB, GW, and FRB, it is found that the main contribution to the constraints comes from the CMB and GW data, but the inclusion of FRB data still can evidently improve the constraint on the baryon density.

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.

A small cluster of massive stars residing in the Galactic center, collectively known as IRS13E, is of special interest due to its close proximity to Sgr A* and the possibility that an embedded intermediate-mass black hole (IMBH) binds its member stars. It has been suggested that colliding winds from two member stars, both classified as Wolf-Rayet type, are responsible for the observed X-ray, infrared and radio emission from IRS13E. We have conducted an in-depth study of the X-ray spatial, temporal and spectral properties of IRS13E, based on 5.6 Ms of ultra-deep Chandra observations obtained over 20 years. These X-ray observations show no significant evidence for source variability. We have also explored the kinematics of the cluster members, using Keck near-infrared imaging and spectroscopic data on a 14-yr baseline that considerably improve the accuracy of stars' proper motions. The observations are interpreted using 3-dimensional hydrodynamical simulations of colliding winds tailored to match the physical conditions of IRS13E, leading us to conclude that the observed X-ray spectrum and morphology can be well explained by the colliding wind scenario, in the meantime offering no support for the presence of a putative IMBH. An IMBH more massive than a few $10^3{\rm~M_\odot}$ is also strongly disfavored by the stellar kinematics.

An increasing number of young circumstellar disks show strikingly ordered (sub)millimeter polarization orientations along the minor axis, which is strong evidence for polarization due to scattering by ~0.1 mm sized grains. To test this mechanism further, we model the dust continuum and polarization data of HD 163296, one of the best observed disks with prominent rings and gaps, using the RADMC-3D radiative transfer code. We find that scattering by grains with a maximum size of 90$\mu$m can simultaneously reproduce the polarization observed at 0.87 mm (ALMA Band 7) and the unusually low spectral index of $\alpha$ ~ 1.5 between 0.87 and 1.25 mm (ALMA Band 6) in the optically thick inner disk as a result of more efficient scattering at a shorter wavelength. The relatively low spectral index of $\alpha$ ~ 2.5 inferred for the optically thin gaps is also reproduced by the same (relatively small) grains, as a result of telescope beam averaging of the gaps (with an intrinsic $\alpha$ ~ 4) and their adjacent optically thick rings (where $\alpha$ << 2). In this case, the long-standing tension between the grain sizes inferred from polarization and spectral index disappears because the relatively low $\alpha$ values are illusory and do not require large mm-sized grains. In addition, the polarization fraction has a unique pattern of azimuthal variation: higher along the major axis than the minor axis in the gaps but higher along the minor axis than the major axis in the rings. We find a rapidly declining polarization spectrum (with the fraction $p \propto \lambda^{-3}$ approximately) in the gaps, which becomes flattened or even inverted towards short wavelengths in the optically thick rings. These contrasting behaviors in the rings and gaps provide further tests of scattering-induced polarization that can be tested via multi-wavelength observations that resolve the disk substructure.

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

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.

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.

Quenched massive spiral galaxies have attracted great attention recently, as more data is available to constrain their environment and cold gas content. However, the quenching mechanism is still uncertain, as it depends on the mass range and baryon budget of the galaxy. In this letter, we report the identification of a rare population of very massive, quenched spiral galaxies with stellar mass $\gtrsim10^{11}{\rm~M_\odot}$ and halo mass $\gtrsim10^{13}{\rm~M_\odot}$ from the Sloan Digital Sky Survey at redshift $z\sim0.1$. Our CO observations using the IRAM-30m telescope show that these galaxies contain only a small amount of molecular gas. Similar galaxies are also seen in the state-of-the-art semi-analytical models and hydro-dynamical simulations. It is found from these theoretical models that these quenched spiral galaxies harbor massive black holes, suggesting that feedback from the central black holes has quenched these spiral galaxies. This quenching mechanism seems to challenge the popular scenario of the co-evolution between massive black holes and massive bulges.

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.

Quenched massive spiral galaxies have attracted great attention recently, as more data is available to constrain their environment and cold gas content. However, the quenching mechanism is still uncertain, as it depends on the mass range and baryon budget of the galaxy. In this letter, we report the identification of a rare population of very massive, quenched spiral galaxies with stellar mass $\gtrsim10^{11}{\rm~M_\odot}$ and halo mass $\gtrsim10^{13}{\rm~M_\odot}$ from the Sloan Digital Sky Survey at redshift $z\sim0.1$. Our CO observations using the IRAM-30m telescope show that these galaxies contain only a small amount of molecular gas. Similar galaxies are also seen in the state-of-the-art semi-analytical models and hydro-dynamical simulations. It is found from these theoretical models that these quenched spiral galaxies harbor massive black holes, suggesting that feedback from the central black holes has quenched these spiral galaxies. This quenching mechanism seems to challenge the popular scenario of the co-evolution between massive black holes and massive bulges.

In this paper, we present a numerical method to study the predicted lightcurves as a function of viewing angle. We extrapolate the fitting formulae for the mass and velocity of tidal dynamical ejecta across a wide range of mass ratio validated with 66 simulations and use them in the calculations of the kilonova lightcurves. The calculated peak luminosity of a BH-NS merger kilonova is typically about a few times $10^{41}\ {\rm erg\ s^{-1}}$, which is always $\lesssim4.5\times10^{41}\ {\rm erg\ s^{-1}}$. This corresponds to the AB absolute magnitudes fainter than $\sim -15\ {\rm mag}$ in optical and $\sim -16\ {\rm mag}$ in infrared. Since the projected photosphere area of the dynamical ejecta is much larger than that of the disk wind outflows, the dynamical ejecta usually contribute to the majority of the kilonova emission from BH-NS mergers. The fitted blackbody temperature and the shape of the observed multi-band lightcurves are insensitive to the line of sight. The peak time of the observed multi-band lightcurves, affected by the light propagation effect, is related to the relative motion direction between the dynamical ejecta and the observer. The observed luminosity varies with the projected photosphere area determined by the viewing angles. However, the predicted peak luminosity only varies by a factor of $\sim (2 - 3)$ (or by $\sim1\ {\rm mag}$) for different viewing angles. When the short-duration gamma-ray burst afterglow is taken into account, for an on-axis geometry, the kilonova emission is usually outshone by the afterglow emission and can be only observed in the redder bands, especially in the $K$-band at late times. Compared with GW170817/AT2017gfo, the BH-NS merger kilonovae are optically dim but possibly infrared bright. At the same epoch after the merger, the blackbody fitting temperature of the BH-NS merger kilonovae is lower than that of GW170817/AT2017gfo.

We perform an extensive review of the numerous studies and methods used to determine the total mass of the Milky Way. We group the various methods into seven broad classes, including: i) estimating Galactic escape velocity using high velocity objects; ii) measuring the rotation curve through terminal and circular velocities; iii) modeling halo stars, globular clusters and satellite galaxies with the Spherical Jeans equation and iv) with phase-space distribution functions; v) simulating and modeling the dynamics of stellar streams and their progenitors; vi) modeling the motion of the Milky Way, M31 and other distant satellites under the framework of Local Group timing argument; and vii) measurements made by linking the brightest Galactic satellites to their counterparts in simulations. For each class of methods, we introduce their theoretical and observational background, the method itself, the sample of available tracer objects, model assumptions, uncertainties, limits and the corresponding measurements that have been achieved in the past. Both the measured total masses within the radial range probed by tracer objects and the extrapolated virial masses are discussed and quoted. We also discuss the role of modern numerical simulations in terms of helping to validate model assumptions, understanding systematic uncertainties and calibrating the measurements. While measurements in the last two decades show a factor of two scatters, recent measurements using \textit{Gaia} DR2 data are approaching a higher precision. We end with a detailed discussion of future developments, especially as the size and quality of the observational data will increase tremendously with current and future surveys. In such cases, the systematic uncertainties will be dominant and thus will necessitate a much more rigorous testing and characterization of the various mass determination methods.

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.

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

It has been recently suggested that (i) nuclear rings in barred galaxies (including our own Milky Way) form at the radius where the shear parameter of the rotation curve reaches a minimum; (ii) the acoustic instability of Montenegro et al. is responsible for driving the turbulence and angular momentum transport in the central regions of barred galaxies. Here we test these suggestions by running simple hydrodynamical simulations in a logarithmic barred potential. Since the rotation curve of this potential is scale-free, the shear minimum theory predicts that no ring should form. We find that in contrast to this prediction, a ring does form in the simulation, with morphology consistent with that of nuclear rings in real barred galaxies. This proves that the presence of a shear-minimum is not a necessary condition for the formation of a ring. We also find that perturbations that are predicted to be acoustically unstable wind up and eventually propagate off to infinity, so that the system is actually stable. We conclude that (i) the shear-minimum theory is an unlikely mechanism for the formation of nuclear rings in barred galaxies; (ii) the acoustic instability is a spurious result and may not be able to drive turbulence in the interstellar medium, at least for the case without self-gravity. The question of the role of turbulent viscosity remains open.

Two new fast single radio bursts FRB 180924 and FRB 190523 well localized to massive galaxies have opened a new window to probe and characterize how cosmic baryons are allocated between galaxies, their surroundings and intergalactic medium. We are motivated by testing Einstein's weak equivalence principle with these two cosmic transients which have accurate redshifts. Using photons with different energies emitted by FRB 180924, we obtain, so far, the most stringent bound $\Delta\gamma<2.16\times10^{-10}$ for non-repeating FRBs with accurate redshifts when only considering the gravitational potential of the Milk Way. If using the gravitational potential of the Laniakea supercluster instead of the Milk Way one, we also obtain the strictest bound $\Delta\gamma<1.06\times10^{-14}$ to date. In light of rapid progress of FRB cosmology, towards the next two decades, we give an universal limitation $\Delta\gamma<8.24\times10^{-22}$ from photons with different energies emitted by single FRBs with accurate redshifts. Moreover, we analyze detailedly the effects of various astrophysical parameters on the precision of weak equivalence principle. We also estimate the abilities of single FRBs with known redshifts to test the validity of swampland criterion, and to distinguish which value of $H_0$ is preferred.

Five fast radio bursts (FRBs), including three apparently non-repeating ones FRB 180924, FRB 181112, and FRB 190523, and two repeaters, FRB 121102 and FRB 180916.J0158+65, have already been localized so far. We apply a method developed recently by us (Li et al. 2019) to these five localized FRBs to give a cosmology-insensitive estimate of the fraction of baryon mass in the IGM, $f_{\rm IGM}$. Using the measured dispersion measure (DM) and luminosity distance $d_{\rm L}$ data (inferred from the FRB redshifts and $d_{\rm L}$ of type Ia supernovae at the same redshifts) of the five FRBs, we constrain the local $f_{\rm IGM} = 0.84^{+0.16}_{-0.22}$ with no evidence of redshift dependence. This cosmology-insensitive estimate of $f_{\rm IGM}$ from FRB observations is in excellent agreement with previous constraints using other probes. Moreover, using the three apparently non-repeating FRBs only we get a little looser but consistent result $f_{\rm IGM} = 0.74^{+0.24}_{-0.18}$. In these two cases, reasonable estimations for the host galaxy DM contribution (${\rm DM_{host}}$) can be achieved by modelling it as a function of star formation rate. The constraints on both $f_{\rm IGM}$ and ${\rm DM_{host}}$ are expected to be significantly improved with the rapid progress in localizing FRBs.

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M=1.23_{-0.18}^{+0.21}\times 10^{12} M_\odot$ and a concentration of $c=9.4_{ -2.1}^{ +2.8}$ with the prior based on the $M$-$c$ relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. By varying the disruption rate of such satellites in this simulation, we estimate a $\sim 5\%$ scatter in the inferred MW halo mass among hydrodynamics-based simulations.

Polarized continuum emission from aligned grains in disks around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively or mechanically. We show that it can also be used to probe another key disk property -- the temperature gradient -- along sight lines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth tau at the tau=1 surface. The formula is then validated using models of stellar irradiated disks with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I disks with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212 and HH 111. We find, in particular, that the sight lines passing through the near-side of a highly inclined disk trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disk may be such a case. Our technique for probing the disk temperature gradient through dust polarization can complement other methods, particularly those using molecular lines.

The co-evolution between black holes and galaxies suggests that feedback of active galactic nuclei influence host galaxies through ejecting radiative and kinetic energies to surroundings. Larger scale outflow in local universe are frequently observed by spatially resolved spectroscopy, while smaller scale outflow cannot be directly resolved by current observations. At the scale of the dusty torus, radiative and kinetic energies ejected from the central active nucleus interact with the materials. However, observations of such outflow are rarely reported due to the lack detection of unambiguously gas emission. Here we report the detection of clear and rich emission lines origin from the scale of dusty tours in an partially obscured quasar. The lines share a common intermediate width with full width at half maximum about 1900 \kmps\ and are shown in two systems: a major system is unshifted and a minor system has a blue-shifts of about 2600 \kmps. The line intensity ratios, combining photo-ionization simulations, indicates an ultradense line-emitting region with the density as high as $\sim$ $10^{13}~\rm cm^{-3}$. We interpret this as the lines being excited by a shock induced by the high-density and high-temperature gases at the scale of dusty torus, rather than photo-ionized by the central accretion disk. We speculate that the outflow, launched from the accretion disk, collides onto the inner wall of the dusty torus and shock-heat the gases to cause the major emission lines. The outflowing gases may also collide onto surrounding isolated clouds, and give rise to blue-shifted minor emission lines.

The Local Group dwarf elliptical galaxy M32 hosts one of the nearest and most under-luminous super-massive black holes (SMBHs) ever known, offering a rare opportunity to study the physics of accreting SMBHs at the most quiescent state. Recent Very Large Array (VLA) observations have detected a radio source at the nucleus of M32, which is suggested to be the radio counterpart of the SMBH. To further investigate the radio properties of this nuclear source, we have conducted follow-up, high-resolution VLA observations in four epochs between 2015--2017, each with dual frequencies. At 6 GHz, the nuclear source is resolved under an angular resolution of $\sim$0\farcs4, exhibiting a coreless, slightly lopsided morphology with a detectable extent of $\sim$2.5 \arcsec ($\sim$10 parsec). No significant variability can be found among the four epochs. At 15 GHz, no significant emission can be detected within the same region, pointing to a steep intrinsic radio spectrum (with a 3\,$\sigma$ upper limit of -1.46 for the spectral index). We discuss possible scenarios for the nature of this nuclear source and conclude that a stellar origin, in particular planetary nebulae, X-ray binaries, supernova remnants or diffuse ionized gas powered by massive stars, can be ruled out.Instead, the observed radio properties can be explained by synchrotron radiation from a hypothetical wind driven by the weakly accreting SMBH.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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.

We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of cosmic muons. It's able to image the overburden density length from the surface of overburden to the detector along the muon tracks. The open sky muon flux measurement outside the tunnel has a good agreement with the modified Gassier Formula model. The distributions of muon flux at three positions inside the tunnel are very similar to that of open sky. Assuming the average density of overburden compact sandstone is $2.65 g/cm^3$, the overburden thickness can be obtained from the density length derived from the difference of muon flux inside and outside the tunnel. Moreover, for known penetrated lengths (i.e., topography of overburden), the density anomalies of the overburden can also been obtained. This study suggests a potential application for imaging and detecting subsurface structures in civil engineering, tunnels or caverns with the cosmic ray muon telescope.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

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.

The point spread function reflects the state of an optical telescope and it is important for data post-processing methods design. For wide field small aperture telescopes, the point spread function is hard to model, because it is affected by many different effects and has strong temporal and spatial variations. In this paper, we propose to use the denoising autoencoder, a type of deep neural network, to model the point spread function of wide field small aperture telescopes. The denoising autoencoder is a pure data based point spread function modelling method, which uses calibration data from real observations or numerical simulated results as point spread function templates. According to real observation conditions, different levels of random noise or aberrations are added to point spread function templates, making them as realizations of the point spread function, i.e., simulated star images. Then we train the denoising autoencoder with realizations and templates of the point spread function. After training, the denoising autoencoder learns the manifold space of the point spread function and can map any star images obtained by wide field small aperture telescopes directly to its point spread function, which could be used to design data post-processing or optical system alignment methods.

We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows.

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.

Double Degenerate systems (DDs) are supposed to be significant gravitational wave (GW) sources for future space-based gravitational-wave detectors, e.g., Laser Interferometer Space Antenna (LISA). Recently, one type of DDs with Extremely low-mass WD (ELM WD; $\leq 0.30\; M_\odot$) companions has been largely found in the ELM Survey. They have very short orbital periods and are therefore important sources for LISA detection. Besides, due to the thick envelope of ELM WDs compared with massive WDs (e.g. CO WDs), they are much easier to be found by the combination of electromagnetic (EM) and GW observations. In this paper, we first obtain the population of ELM WDs in DDs with considering the detailed evolutionary tracks of ELM WDs, and then analyse the GW radiation of these systems. We found that about $6\times10^3$ sources could be solely detected by LISA, including $\sim2\times10^3$ chirping sources, and $\sim13$ ($\sim107$) more sources are expected to be detected by both LISA and ELM Survey (Gaia).

The symmetry energy obtained with the effective Skyrme energy density functional is related to the values of isoscalar effective mass and isovector effective mass, which is also indirectly related to the incompressibility of symmetric nuclear matter. In this work, we analyze the values of symmetry energy and its related nuclear matter parameters in five-dimensional parameter space by describing the heavy ion collision data, such as isospin diffusion data at 35 MeV/u and 50 MeV/u, neutron skin of $^{208}$Pb, and tidal deformability and maximum mass of neutron star. We obtain the parameter sets which can describe the isospin diffusion, neutron skin, tidal deformability and maximum mass of neutron star, and give the incompressibility $K_0$=250.23$\pm$20.16 MeV, symmetry energy coefficient $S_0$=31.35$\pm$2.08 MeV, the slope of symmetry energy $L$=59.57$\pm$10.06 MeV, isoscalar effective mass $m_s^*/m$=0.75$\pm$0.05 and quantity related to effective mass splitting $f_I$=0.005$\pm$0.170. At two times normal density, the symmetry energy we obtained is in 35-55 MeV. To reduce the large uncertainties of $f_I$, more critical works in heavy ion collisions at different beam energies are needed.

IGR J17591-2342 is a recently INTEGRAL discovered accreting millisecond X-ray pulsar that went into outburst around July 21, 2018. To better understand the physics acting in these systems during the outburst episode we performed detailed temporal-, timing- and spectral analyses across the 0.3-300 keV band using data from NICER, XMM-Newton, NuSTAR and INTEGRAL. The hard X-ray 20-60 keV outburst profile is composed of four flares. During the maximum of the last flare we discovered a type-I thermonuclear burst in INTEGRAL JEM-X data. We derived a distance of 7.6+/-0.7 kpc, adopting Eddington luminosity limited photospheric radius expansion burst emission and assuming anisotropic emission. In the timing analysis using all NICER 1-10 keV monitoring data we observed a rather complex behaviour starting with a spin-up period, followed by a frequency drop, a episode of constant frequency and concluding with irregular behaviour till the end of the outburst. The 1-50 keV phase distributions of the pulsed emission, detected up to $\sim$ 120 keV using INTEGRAL ISGRI data, was decomposed in three Fourier harmonics showing that the pulsed fraction of the fundamental increases from ~10% to ~17% going from ~1.5 to ~4 keV, while the harder photons arrive earlier than the soft photons for energies <10 keV. The total emission spectrum of IGR J17591-2342 across the 0.3-150 keV band could adequately be fitted in terms of an absorbed compPS model yielding as best fit parameters a column density of N_H=(2.09+/-0.05) x 10^{22} /cm2, a blackbody seed photon temperature kT_bb,seed of 0.64+/- 0.02 keV, electron temperature kT_e=38.8+/-1.2 keV and Thomson optical depth Tau_T=1.59+/-0.04. The fit normalisation results in an emission area radius of 11.3+/-0.5 km adopting a distance of 7.6 kpc. Finally, the results are discussed within the framework of accretion physics- and X-ray thermonuclear burst theory.

We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbinary disk in a highly anisotropic manner. Stream-like SO emission around the circumbinary disk reveals the complicated shocks caused by the interactions between the disk, accretion flows and outflows. A wall-like structure is discovered south of VLA1623B. The discovery of two outflow cavity walls at the same position traveling at different velocities suggests the two outflows from both VLA1623A and VLA1623B overlays on top of each other in the plane of sky. Our detailed flat and flared disk modeling shows that Cycle 2 C18O J = 2-1 data is inconsistent with the combined binary mass of 0.2 Msun as suggested by early Cycle 0 studies. The combined binary mass for VLA1623A should be modified to 0.3 ~ 0.5 Msun.

Using archival {\it Chandra} observations with a total effective exposure of 323 ks, we derive an updated catalog of point sources in the bulge globular cluster M28. The catalog contains 502 X-ray sources within an area of $\sim475\, \rm arcmin^{2}$, and more than $90\%$ of these sources are first detected in this cluster. We find significant dips in the radial distribution profiles of X-ray sources in M28, with the projected distance and width of the distribution dip for bright ($L_{X} \gtrsim 4.5\times 10^{30} {\rm\ erg\ \,s^{-1}}$) X-ray sources are larger than the faint ($L_{X} \lesssim 4.5\times 10^{30} {\rm\ erg\ \,s^{-1}}$) sources. The "generalized King model" fitting give a slightly larger average mass for the bright sources ($1.30\pm0.15\,M_{\odot}$) than the faint sources ($1.09\pm0.14\,M_{\odot}$), which support a universal mass segregation delay between heavy objects in GCs. Compared with 47 Tuc and Terzan 5, we show that the dynamical age of M28 is comparable to Terzan 5 and much smaller than 47 Tuc, but it is evolving more fast (i.e., with smaller two-body relaxation timescale) than 47 Tuc. These features may suggest an acceleration effect of cluster dynamical evolution by tidal shock in M28. Besides, we find an abnormal deficiency of X-ray sources in the central region ($R \lesssim 1.5 \rm~arcmin$) of M28 than its outskirts, which indicate that M28 may have suffered an early phase of primordial binary disruption within its central region, and mass segregation effect will erase such a phenomenon as cluster evolve to older dynamical age.

In the literature, compact binary coalescences (CBCs) have been proposed as one of the main scenarios to explain the origin of some non-repeating fast radio bursts (FRBs). The large discrepancy between the FRB and CBC event rate densities suggest their associations, if any, should only apply at most for a small fraction of FRBs. Through a Bayesian estimation method, we show how a statistical analysis of the coincident associations of FRBs with CBC gravitational wave (GW) events may test the hypothesis of these associations. We show that during the operation period of advanced LIGO, the detection of $\sim100$ ($\sim1000$) GW-less FRBs with dispersion measure (DM) values smaller than 500 ${\rm pc~cm^{-3}}$ could reach the constraint that less than $10\%$ (or $1\%$) FRBs are related to binary black hole (BBH) mergers. The same number of FRBs with DM values smaller than 100 ${\rm pc~cm^{-3}}$ is required to reach the same constraint for binary neutron star (BNS) mergers. With the upgrade of GW detectors, the same constraints for BBH and BNS mergers can be reached with less FRBs or looser requirements for the DM values. It is also possible to pose constraints on the fraction of each type of CBCs that are able to produce observable FRBs based on the event density of FRBs and CBCs. This would further constrain the dimensionless charge of black holes in binary BH systems.

We study the statistical properties of 320 bulges of disk galaxies in the Carnegie-Irvine Galaxy Survey, using robust structural parameters of galaxies derived from image fitting. We apply the Kormendy relation to classify classical and pseudo bulges and characterize bulge dichotomy with respect to bulge structural properties and physical properties of host galaxies. We confirm previous findings that pseudo bulges on average have smaller S\'{e}rsic indices, smaller bulge-to-total ratios, and fainter surface brightnesses when compared with classical bulges. Our sizable sample statistically shows that pseudo bulges are more intrinsically flattened than classical bulges. Pseudo bulges are most frequent (incidence $\gtrsim 80\%$) in late-type spirals (later than Sc). Our measurements support the picture in which pseudo bulges arose from star formation induced by inflowing gas, while classical bulges were born out of violent processes such as mergers and coalescence of clumps. We reveal differences with the literature that warrant attention: (1) the bimodal distribution of S\'{e}rsic indices presented by previous studies is not reproduced in our study; (2) classical and pseudo bulges have similar relative bulge sizes; and (3) the pseudo bulge fraction is considerably smaller in early-type disks compared with previous studies based on one-dimensional surface brightness profile fitting. We attribute the above differences to our improved image quality, more robust bulge-to-disk decomposition technique, and different classification criteria applied. Moreover, we find that barred galaxies do not host more pseudo bulges or more prominent pseudo bulges than unbarred galaxies. Various implications of these findings are discussed.

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.

One of the most important checks for systematic errors in CMB studies is the cross correlation of maps made by independent experiments. In this paper we report on the cross correlation between maps from the Atacama B-mode Search (ABS) and Atacama Cosmology Telescope (ACT) experiments in both temperature and polarization. These completely different measurements have a clear correlation with each other and with the Planck satellite in both the EE and TE spectra at $\ell<400$ over the roughly $1100$ deg$^2$ common to all three. The TB, EB, and BB cross spectra are consistent with noise. Exploiting such cross-correlations will be important for future experiments operating in Chile that aim to probe the $30<\ell<8,000$ range.

Taking a large amount of images, the Cassini Imaging Science Subsystem (ISS) has been routinely used in astrometry. In ISS images, disk-resolved objects often lead to false detection of stars that disturb the camera pointing correction. The aim of this study was to develop an automated processing method to remove the false image stars in disk-resolved objects in ISS images. The method included the following steps: extracting edges, segmenting boundary arcs, fitting circles and excluding false image stars. The proposed method was tested using 200 ISS images. Preliminary experimental results show that it can remove the false image stars in more than 95% of ISS images with disk-resolved objects in a fully automatic manner, i.e. outperforming the traditional circle detection based on Circular Hough Transform (CHT) by 17%. In addition, its speed is more than twice as fast as that of the CHT method. It is also more robust (no manual parameter tuning is needed) when compared with CHT. The proposed method was also applied to a set of ISS images of Rhea to eliminate the mismatch in pointing correction in automatic procedure. Experiment results showed that the precision of final astrometry results can be improve by roughly 2 times than that of automatic procedure without the method. It proved that the proposed method is helpful in the astrometry of ISS images in fully automatic manner.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

Massive neutrinos influence the background evolution of the Universe as well as the growth of structure. Being able to model this effect and constrain the sum of their masses is one of the key challenges in modern cosmology. Weak gravitational lensing by the large-scale structure has proven to be a powerful tool to achieve this goal, and its importance to precision cosmology is borne out in the scientific results of galaxy surveys such as the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES) and Hyper SuprimeCam (HSC). Weak-lensing cosmological constraints will also soon reach higher levels of precision with next-generation surveys like LSST, WFIRST and Euclid. In this context, we use the MassiveNus simulations to derive constraints on the sum of neutrino masses $M_{\nu}$, the present-day total matter density ${\Omega}_m$, and the primordial power spectrum normalization $A_s$ in a tomographic setting. Along with the lensing power spectrum, we use peak counts as higher-order statistics. We find that the combination of the two statistics does not add any relevant information over considering just the peaks alone. We also compare two different filtering techniques on noisy maps and show that a starlet (wavelet) filter produces tighter constraints than Gaussian filtering by 50%, 25%, and 38% on $M_{\nu}$, ${\Omega}_m$, and $A_s$, respectively. For the starlet case, we further find a minimum resolution that allows us to obtain constraints comparable to what is achieved with the full wavelet decomposition, and we show that the information contained in the coarse map cannot be neglected.

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

We present 2D-spectroscopic observations from Centro Astron\'omico Hispano Alem\'an (CAHA) 3.5 m telescope and the millimetre observation from NOrthern Extended Millimeter Array (NOEMA) of the nearby S0 galaxy PGC 26218, which shows central star-formation activity and post-starburst outside in the disk. We estimate the star formation rate (SFR = $0.28\pm0.01$ $M_{\odot} \rm yr^{-1}$) and molecular gas mass ($\rm log\ $$M_{\rm H_{2}}=7.60\pm0.15\ M_{\odot}$) of PGC 26218 based on the extinction-corrected H$\alpha$ emission line and the CO-$\rm H_{2}$ conversion factor ($\alpha_{\rm CO}$) of the Milky Way, respectively. We find that PGC 26218 follows the star forming main sequence (SFMS) and the Kennicutt-Schmidt law. Comparing the kinematics of CO($J$=1-0), stars and H$\alpha$, we find that the rotational axis of CO($J$=1-0) is 45$^{\circ}$ different from that of H$\alpha$. In addition, the profile of the CO($J$=1-0) emission line shows asymmetry and has an inflow component of $\sim$ 46 $\rm km\ s^{-1}$. With the decomposition of the optical image, we confirm that PGC 26218 shows multiple nuclear structures. The projected offset between the most luminous optical center and the center of CO($J$=1-0) is $5.2\arcsec$ ($\sim$ 0.6 kpc) and the latter overlaps with one of the optical cores. These results support that PGC 26218 may have experienced a gas-rich minor merger, extending its star formation and locating it in the SFMS.

The on-going vertical phase mixing, manifesting itself as a snail shell in the $Z-V_{Z}$ phase space, has been discovered with the Gaia DR2 data. To better understand the origin and properties of the phase mixing process, we study the vertical phase-mixing signatures in arches (including the classical ``moving groups'') of the $V_{R}-V_{\phi}$ phase space near the Solar circle. Interestingly, the phase space snail shell exists only in the arches with $|V_{\phi} - V_{\rm LSR}| \lesssim 30$ km/s, i.e., stars on dynamically ``colder'' orbits. The snail shell becomes much weaker and eventually disappears for increasingly larger radial action ($J_{R}$), quantifying the ``hotness'' of orbits. Thus one should pay closer attention to the colder orbits in future phase mixing studies. We also confirm that the Hercules stream has two branches (at fast and slow $V_{\phi}$), which may not be explained by a single mechanism, since only the fast branch shows the prominent snail shell feature. The hotter orbits may have phase-wrapped away already due to the much larger dynamical range in radial variation to facilitate faster phase mixing. To explain the lack of a well-defined snail shell in the hotter orbits, the disk should have been perturbed at least $500$ Myr ago. Our results offer more support to the recent satellite-disk encounter scenario than the internal bar buckling perturbation scenario as the origin of the phase space mixing. Origin of the more prominent snail shell in the $V_{\phi}$ color-coded phase space is also discussed.

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

The synchrotron radiation from secondary electrons and positrons (SEPs) generated by hadronic interactions in the shock of supernova remnant (SNR) could be a distinct evidence of cosmic ray (CR) production in SNR shocks. Here we provide a method where the observed gamma-ray flux from SNRs, created by pion decays, is directly used to derive the SEP distribution and hence the synchrotron spectrum. We apply the method to three gamma-ray bright SNRs. In the young SNR RX J1713.7-3946, if the observed GeV-TeV gamma-rays are of hadronic origin and the magnetic field in the SNR shock is $B\gtrsim 0.5$mG, the SEPs may produce a spectral bump at $10^{-5}-10^{-2}$eV, exceeding the predicted synchrotron component of the leptonic model, and a soft spectral tail at $\gtrsim 100$keV, distinct from the hard spectral slope in the leptonic model. In the middle-aged SNRs IC443 and W44, if the observed gamma-rays are of hadronic origin, the SEP synchrotron radiation with $B\sim 400 - 500 \mu$G can well account for the observed radio flux and spectral slopes, supporting the hadronic origin of gamma-rays. Future microwave to far-infrared and hard X-ray (>100keV) observations are encouraged to constraining the SEP radiation and the gamma-ray origin in SNRs.

The fuzzy dark matter model (FDM, also known as quantum wave dark matter model) argues that light bosons with a mass of $\sim10^{-22}{\;{\rm eV}}$ are a possible candidate for dark matter in the Universe. One of the most important predictions of FDM is the formation of a soliton core instead of a density cusp at the center of galaxies. If FDM is the correct theory of dark matter, then the predicted soliton core can help to form the Central Molecular Zone (CMZ) in the Milky Way. We present high-resolution hydrodynamical simulations of gas flow patterns to constrain the properties of the soliton core based on a realistic Milky Way potential. We find that a dense center is required to form a reasonable CMZ. The size and kinematics of the CMZ offer a relatively strong constraint on the inner enclosed mass profile of the Galaxy. If a soliton core is not considered, a compact nuclear bulge alone with a radially varying mass-to-light ratio can match the observed size and kinematics of the CMZ. A soliton core model with a mass of $\approx4.0\times10^8{\; {\rm M}_{\odot}}$ and a core radius of $\approx0.05{\;{\rm kpc}}$, together with a less massive nuclear bulge with a constant mass-to-light ratio, also agrees nicely with the current data. Such a FDM soliton core corresponds to a boson mass of $\sim2-7\times10^{-22}{\;{\rm eV}}$, which could be further constrained by the improved determination of the mass-to-light ratio in the Galactic center.

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 explore the survival of cool clouds in multi-phase circum-galactic media. We revisit the "cloud crushing problem" in a large survey of simulations including radiative cooling, self-shielding, self-gravity, magnetic fields, and anisotropic Braginskii conduction and viscosity (with saturation). We explore a wide range of parameters including cloud size, velocity, ambient temperature and density, as well as a variety of magnetic field configurations and cloud turbulence. We find that realistic magnetic fields and turbulence have weaker effects on cloud survival; the most important physics is radiative cooling and conduction. Self-gravity and self-shielding are important for clouds which are initially Jeans-unstable, but largely irrelevant otherwise. Non-self-gravitating, realistically magnetized clouds separate into four regimes: (1) At low column densities, clouds evaporate rapidly via conduction. (2) A "failed pressure confinement" regime, where the ambient hot gas cools too rapidly to provide pressure confinement for the cloud. (3) An "infinitely long-lived" regime, in which the cloud lifetime becomes longer than the cooling time of gas swept up in the leading bow shock, so the cloud begins to accrete and grow. (4) A "classical cloud destruction" regime, where clouds are eventually destroyed by instabilities. In the final regime, the cloud lifetime can exceed the naive cloud-crushing time owing to conduction-induced compression. However, small and/or slow-moving clouds can also evaporate more rapidly than the cloud-crushing time. We develop simple analytic models that explain the simulated cloud destruction times in this regime.

Context: Subdwarf B stars (sdBs) play a crucial role in stellar evolution, asteroseismology, and far-UV radiation of early-type galaxies, and have been intensively studied with observation and theory. It has theoretically been predicted that sdBs with neutron star (NS) companions exist in the Galaxy, but none have been discovered yet. This remains a puzzle in this field. In a previous study (hereafter Paper I), we have studied the formation channels of sdB+NS binaries from main-sequence (MS) stars plus NS binaries by establishing a model grid, but it is still unclear how these binaries consisting of MS stars and NS binaries came to be in the first place. Aims: We systematically study the formation of sdB+NS binaries from their original zero-age main-sequence progenitors. We bridge the gap left by our previous study in this way. We obtain the statistical population properties of sdB+NS binaries and provide some guidance for observational efforts. Methods: We first used Hurley's rapid binary evolution code BSE to evolve 10^7 primordial binaries to the point where the companions of NS+MS, NS+Hertzsprung gap (HG) star, and NS+Giant Branch (GB) star binaries have just filled their Roche lobes. Next, we injected these binaries into the model grid we developed in Paper I to obtain the properties of the sdB+NS populations. We adopted two prescriptions of NS natal kicks. Different values of common-envelope ejection efficiency were chosen to examine the effect of common-envelope evolution on the results. Conclusions: Most sdB+NS binaries are located in the Galactic disk with small RV semi-amplitudes. SdB+NS binaries with large RV semi-amplitudes are expected to be strong GWR sources, some of which could be detected by LISA in the future.

We report the discovery of SDSS J0849+1114 as the first known triple Type 2 Seyfert nucleus. It represents three active black holes that are identified from new spatially resolved optical slit spectroscopy using the Dual Imaging Spectrograph on the 3.5 m telescope at the Apache Point Observatory. We also present new complementary observations including the Hubble Space Telescope Wide Field Camera 3 U- and Y-band imaging, Chandra Advanced CCD Imaging Spectrometer S-array X-ray 0.5--8 keV imaging spectroscopy, and NSF Karl G. Jansky Very Large Array radio 9.0 GHz imaging in its most extended A configuration. These comprehensive multiwavelength observations, when combined together, strongly suggest that all three nuclei are active galactic nuclei. While they are now still at kiloparsec-scale separations, where the host-galaxy gravitational potential dominates, the black holes may evolve into a bound triple system in $\lesssim$2 Gyr. These triple merger systems may explain the overly massive stellar cores that have been observed in some elliptical galaxies such as M87, which are expected to be unique gravitational wave sources. Similar systems may be more common in the early universe, when galaxy mergers are thought to have been more frequent.

Galaxy clustering provides insightful clues to our understanding of galaxy formation and evolution, as well as the universe. The redshift assignment for the random sample is one of the key steps to measure the galaxy clustering accurately. In this paper, by virtue of the mock galaxy catalogs, we investigate the effect of two redshift assignment methods on the measurement of galaxy two-point correlation functions (hereafter 2PCFs), the Vmax method and the "shuffled" method. We found that the shuffled method significantly underestimates both of the projected 2PCFs and the two-dimensional 2PCFs in redshift space. While the Vmax method does not show any notable bias on the 2PCFs for volume-limited samples. For flux-limited samples, the bias produced by the Vmax method is less than half of the shuffled method on large scales. Therefore, we strongly recommend the Vmax method to assign redshifts to random samples in the future galaxy clustering analysis.

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

Model-independent measurements for the cosmic spatial curvature, which is related to the nature of cosmic space-time geometry, plays an important role in cosmology. On the basis of the Distance Sum Rule in the Friedmann-Lema{\^i}tre-Robertson-Walker metric, (distance ratio) measurements of strong gravitational lensing (SGL) systems together with distances from type Ia supernovae observations have been proposed to directly estimate the spatial curvature without any assumptions for the theories of gravity and contents of the universe. However, previous studies indicated that a spatially closed universe was strongly preferred. In this paper, we re-estimate the cosmic curvature with the latest SGL data which includes 163 well-measured systems. In addition, possible factors, e.g. combination of SGL data from different surveys and stellar mass of the lens galaxy, which might affect estimations for the spatial curvature, are considered in our analysis. We find that, except the case where only SGL systems from the Sloan Lens ACS Survey are considered, a spatially flat universe is consistently favored at very high confidence level by the latest observations. It is suggested that the increasing number of well-measured strong lensing events might significantly reduce the bias of estimation for the cosmic curvature.

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.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

We propose a versatile and accurate method to estimate the halo mass and concentration from the kinematics of satellite galaxies. We construct the 6D phase-space distribution function of satellites from a cosmological simulation based on the similarity of internal dynamics for different halos. Within the Bayesian statistical framework, not only can we infer the halo mass and concentration efficiently, but also treat various observational effects, including the selection function, incomplete data, and measurement errors, in a rigorous and straightforward manner. Through tests with mock samples, we show that our method is valid and accurate, and more precise than pure steady-state methods. It can constrain the halo mass to within ~ 20% using only 20 tracers and has a small intrinsic uncertainty of ~ 10%. In addition to the clear application to the Milky Way and similar galaxies, our method can be extended to galaxy groups or clusters.

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.

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 POLAR detector is a space based Gamma-Ray Burst (GRB) polarimeter sensitive in the 15-500 keV energy range. Apart from its main scientific goal as a Gamma-Ray Burst polarimeter it is also able to detect photons from pulsars in orbit. By using the six-months in-orbit observation data, significant pulsation from the PSR B0531+21 (Crab pulsar) was obtained. In this work, we present the precise timing analysis of the Crab pulsar, together with a phase-resolved spectroscopic study using a joint-fitting method adapted for wide field of view instruments like POLAR. By using single power law fitting over the pulsed phase, we obtained spectral indices ranging from 1.718 to 2.315, and confirmed the spectral evolution in a reverse S shape which is homogenous with results from other missions over broadband. We will also show, based on the POLAR in-orbit performance and Geant4 Monte-Carlo simulation, the inferred capabilities of POLAR-2, the proposed follow-up mission of POLAR on board the China Space Station (CSS), for pulsars studies.

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

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.

The most luminous quasars at high redshift harbour the fastest-growing and most massive black holes in the early Universe. They are exceedingly rare and hard to find. Here, we present our search for the most luminous quasars in the redshift range from $z=4.5$ to $5$ using data from SkyMapper, Gaia and WISE. We use colours to select likely high-redshift quasars and reduce the stellar contamination of the candidate set with parallax and proper motion data. In $\sim$12,500~deg$^2$ of Southern sky, we find 92 candidates brighter than $R_p=18.2$. Spectroscopic follow-up has revealed 21 quasars at $z\ge 4$ (16 of which are within $z=[4.5,5]$), as well as several red quasars, BAL quasars and objects with unusual spectra, which we tentatively label OFeLoBALQSOs at redshifts of $z\approx 1$ to $2$. This work lifts the number of known bright $z\ge 4.5$ quasars in the Southern hemisphere from 10 to 26 and brings the total number of quasars known at $R_p<18.2$ and $z\ge 4.5$ to 42.

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.

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

We present multi-component decomposition of high-quality $R$-band images of 312 disk galaxies from the Carnegie-Irvine Galaxy Survey. In addition to bulges and disks, we successfully model nuclei, bars, disk breaks, nuclear/inner lenses, and inner rings. Our modeling strategy treats nuclear rings and nuclear bars as part of the bulge component, while other features such as spiral arms, outer lenses, and outer rings are omitted from the fits because they are not crucial for accurate bulge measurements. The error budget of bulge parameters includes the uncertainties from sky level measurements and model assumptions. Comparison with multi-component decomposition from the \textit{Spitzer} Survey of Stellar Structure in Galaxies reveals broad agreement for the majority of the overlapping galaxies, but for a considerable fraction of galaxies there are significant differences in bulge parameters caused by different strategies in model construction. We confirm that on average bulge prominence decreases from early to late-type disk galaxies, although the large scatter of bulge-to-total ratios in each morphological bin limits the application of Hubble type as an accurate predictor of bulge-to-total ratio. In contrast with previous studies claiming that barred galaxies host weaker bulges, we find that barred and unbarred spiral galaxies have similar bulge prominence.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

During non-flaring times, the radio flux of the Sun at the wavelength of a few centimeters to several tens of centimeters mostly originates from the thermal bremsstrahlung emission, very similar to the EUV radiation. Owing to such a proximity, it is feasible to investigate the relationship between the EUV emission and radio emission in a quantitative way. In this paper, we reconstruct the radio images of the Sun through the differential emission measure obtained from the multi-wavelength EUV images of the Atmospheric Imaging Assembly on board Solar Dynamic Observatory. Through comparing the synthetic radio images at 6 GHz with those observed by Siberian Radioheliograph, we find that the predicted radio flux is qualitatively consistent with the observed value, confirming thermal origin of the coronal radio emission during non-flaring times. The results further show that the predicted radio flux is closer to the observations in the case of including the contribution of the plasma with temperatures above 3 MK than in the case of only involving the low temperature plasma as was usually done in the era of pre-SDO. We also discuss the applications of the method and uncertainties of the results.

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.

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

Tremaine and Weinberg (TW) proposed a conceptually simple procedure relying on long-slit spectroscopy to measure the pattern speeds of bars ($\Omega_{\rm p}$) in disk galaxies. Using a simulated galaxy, we investigate the potential biases and uncertainties of TW measurements using increasingly popular integral-field spectrographs (IFSs), for which multiple pseudo-slits (and thus independent measurements) can be constructed with a single observation. Most important, to establish the spatial coverage required and ensure the validity of the measurements, the inferred $\Omega_{\rm p}$ must asymptotically converge as the (half-)length of each pseudo-slit used is increased. The requirement for our simulation is to reach $\approx1.3$ times the half-light radius, but this may vary from galaxy to galaxy. Only those slits located within the bar region yield accurate measurements. We confirm that the position angle of the disk is the dominant source of systematic error in TW $\Omega_{\rm p}$ measurements, leading to under/overestimates of tens of percents for inaccuracies of even a few degrees. Recasting the data so that the data grid aligns with the disk major axis leads to slightly reduced uncertainties. Accurate measurements are obtained only for well-defined ranges of the bar angle (relative to the galaxy major axis) $\phi_{\rm bar}$ and the inclination angle $i$, here $10\lesssim\phi_{\rm bar}\lesssim75^{\circ}$ and $105\lesssim\phi_{\rm bar}\lesssim170^{\circ}$ and $15\lesssim i\lesssim70^{\circ}$. The adopted (pseudo-)slit widths, spatial resolution, and (unless extremely aggressive) spatial binning of IFS data have no significant impact on the measurements. Our results thus provide useful guidelines for reliable and accurate direct $\Omega_{\rm p}$ measurements with IFS observations.

We present Chandra ACIS-S X-ray imaging spectroscopy for five dual active galactic nucleus (AGN) candidates. Our targets were drawn from a sample of 1286 [O III]-selected AGN pairs systematically selected from the Sloan Digital Sky Survey Seventh Data Release. Each of the targets contains two nuclei separated by $\sim$3-9 kpc in projection, both of which are optically classified as Type 2 (obscured) AGNs based on diagnostic ratios of the narrow emission lines. Combined with independent, empirical star formation rate estimates based on the host-galaxy stellar continua, the new Chandra X-ray observations allow us to evaluate the dual-AGN hypothesis for each merging system. We confirm two (SDSS J0907+5203 and SDSS J1544+0446) of the five targets as bona-fide dual AGNs. For the other three targets, the existing data are consistent with the dual-AGN scenario, but we cannot rule out the possibility of stellar/shock heating and/or one AGN ionizing both gaseous components in the merger. The average X-ray-to-[O III] luminosity ratio in our targets seems to be systematically smaller than that observed in single AGNs but is higher than that seen in dual AGNs selected from AGNs with double-peaked narrow emission lines. We suggest that the systematically smaller X-ray-to-[O III] luminosity ratio observed in dual AGNs than in single AGNs is due to a high nuclear gas column likely from strong merger-induced inflows. Unlike double-peaked-[O III]-selected dual AGNs, the new sample selected from resolved galaxy pairs are not subject to the orientation bias caused by the double-peak line-of-sight velocity splitting selection, which also contributes to lowering the X-ray-to-[O III] luminosity ratio.

Spatially resolved (sub)millimeter polarization has been detected by ALMA in an increasing number of disks around young stellar objects. The majority of the observations show polarization patterns that are consistent with self-scattering, especially at ALMA Band 7. The inferred sizes of the grains are typically of order 100 microns, which is very different from the millimeter size commonly inferred from the dust opacity index beta. In an effort to resolve this discrepancy, we first introduced the so-called "Coplanar Isotropic Radiation Field" approximation, which enables the computation of the (signed) polarization fraction analytically. With an oft-adopted dust composition (used, e.g., by DSHARP), we find that models with big dust grains produce very small polarization with reversed orientation, which hasn't been observed. The semi-analytic results are validated through Monte Carlo radiative transfer simulations. In these models, the "correct" polarization orientation and the small beta index are mutually exclusive. To resolve this tension, we explore a wide range of dust models, parameterized by their complex refractive index. We find that both the fraction of the polarization and whether it is reversed or not depend on the refractive index in a complex way, and this dependence is mapped out on an n-k plane for a representative dust size distribution (MRN with amax=3mm) and wavelength of 870 microns. In particular, refractory organics dust grains produce reversed polarization, whereas grains of the same size but made of absorptive carbonaceous materials can produce a percent-level polarization that is not reversed; the latter may alleviate the tension between the grain sizes inferred from the scattering-induced polarization and the opacity index beta. We conclude that scattering-induced polarization has the potential to probe not only the sizes of the grains but also their composition.

We present the first results of high-spectral resolution (0.023 km/s) N$_2$H$^+$ observations of dense gas dynamics at core scales (~0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position-velocity (PV) space for all 5 targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming disks nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disk formation.

Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr.

The formation of Low mass X-ray binaries (LMXB) is favored within dense stellar systems such as Globular Clusters (GCs). The connection between LMXB and Globular Clusters has been extensively studied in the literature, but these studies have always been restricted to the innermost regions of galaxies. We present a study of LMXB in GCs within the central 1.5 deg^2 of the Fornax cluster with the aim of confirming the existence of a population of LMXB in intra-cluster GCs and understand if their properties are related to the host GCs, to the environment or/and to different formation channels.

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.

Testing deviations from the $\Lambda$CDM model using the Cosmic Microwave Background (CMB) power spectra requires a pristine understanding of instrumental systematics. In this work we discuss the properties of a new observable ${\cal R}^{TE}_{\ell}$, the correlation coefficient of temperature and E modes. We find that this observable is mostly unaffected by systematics introducing multiplicative biases such as errors in calibration, polarisation efficiency, beam and transfer function measurements. We discuss the dependency of this observable on the cosmological model and derive its statistical properties. We then compute the T-E correlation coefficients of Planck legacy data and compare them with expectations from the Planck best-fit $\Lambda$CDM and foreground model.

Hydrodynamical simulations show that a pair of spiral arms can form in the disk around a rapidly-growing young star and that the arms are crucial in transporting angular momentum as the disk accretes material from the surrounding envelope. Here we report the detection of a pair of symmetric spiral structures in a protostellar disk, supporting the formation of spiral arms in the disk around a forming star. The HH 111 VLA 1 source is a young Class I source embedded in a massive infalling protostellar envelope and is actively accreting, driving the prominent HH 111 jet. Previous observations showed a ring of shock emission around the disk's outer edge, indicating accretion of the envelope material onto the disk at a high rate. Now with ALMA observations of thermal emission from dust particles, we detect a pair of spiral arms extending from the inner region to the disk's outer edge, similar to that seen in many simulations. Additionally, the disk is massive, with Toomre's Q parameter near unity in the outer parts where the spiral structures are detected, supporting the notion that envelope accretion is driving the outer disk gravitationally unstable. In our observations, another source, HH 111 VLA 2, is spatially resolved for the first time, showing a disk-like structure with a diameter of ~ 26 au and an orientation nearly orthogonal to that of the HH 111 VLA 1 disk.

It is well known that the polarized continuum emission from magnetically aligned dust grains is determined to a large extent by local magnetic field structure. However, the observed significant anticorrelation between polarization fraction and column density may be strongly affected, perhaps even dominated by variations in grain alignment efficiency with local conditions, in contrast to standard assumptions of a spatially homogeneous grain alignment efficiency. Here we introduce a generic way to incorporate heterogeneous grain alignment into synthetic polarization observations of molecular clouds, through a simple model where the grain alignment efficiency depends on the local gas density as a power-law. We justify the model using results derived from radiative torque alignment theory. The effects of power-law heterogeneous alignment models on synthetic observations of simulated molecular clouds are presented. We find that the polarization fraction-column density correlation can be brought into agreement with observationally determined values through heterogeneous alignment, though there remains degeneracy with the relative strength of cloud-scale magnetized turbulence and the mean magnetic field orientation relative to the observer. We also find that the dispersion in polarization angles-polarization fraction correlation remains robustly correlated despite the simultaneous changes to both observables in the presence of heterogeneous alignment.

The inference of the Milky Way halo mass requires modelling the phase space structure of dynamical tracers, with different tracers following different models and having different levels of sensitivity to the halo mass. For steady-state models, phase correlations among tracer particles lead to an irreducible stochastic bias. This bias is small for satellite galaxies and dark matter particles, but as large as a factor of 2 for halo stars. This is consistent with the picture that satellite galaxies closely trace the underlying phase space distribution of dark matter particles, while halo stars are less phase-mixed. As a result, the use of only $\sim 100$ satellite galaxies can achieve a significantly higher accuracy than that achievable with a much larger sample of halo stars.

With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 {\mu}m, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in Serpens Emb8(N). We use 3 mm and 1 mm molecular tracers to investigate outflow and dense gas properties and their correlation with the polarization. These observations allow us to investigate the physical processes involved in the Radiative Alignment Torques (RATs) acting on dust grains along the outflow cavity walls, which experience irradiation from accretion processes and outflow shocks. The inner core of SMM1-a presents a polarization pattern with a poloidal magnetic field at the bases of the two lobes of the bipolar outflow. To the south of SMM1-a we see two polarized filaments, one of which seems to trace the redshifted outflow cavity wall. The other may be an accretion streamer of material infalling onto the central protostar. We propose that the polarized emission we see at millimeter wavelengths along the irradiated cavity walls can be reconciled with the expectations of RAT theory if the aligned grains present at '<' 500 au scales in Class 0 envelopes have grown larger than the 0.1 {\mu}m size of ISM dust grains. Our observations allow us to constrain the star-forming sources magnetic field morphologies within the central cores, along the outflow cavity walls, and in possible accretion streamers.

Disks are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disk formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically-induced flattened pseudodisk that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disk at early times with or without ambipolar diffusion, but such a disk remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong ambipolar diffusion is also present. The AD-enabled disks in laminar simulations tend to fragment gravitationally. The disk fragmentation is suppressed by initial turbulence. The ambipolar diffusion facilitates the disk formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and ambipolar diffusion complement each other in promoting disk formation. The disks formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-$\beta$ of order a few tens or smaller, which is 2-3 orders of magnitude lower than the values commonly adopted in MHD simulations of protoplanetary disks. To resolve this potential tension, longer-term simulations of disk formation and evolution with increasingly more realistic physics are needed.

Extragalactic dispersion measures (DMs) obtained from observations of fast radio bursts (FRBs) are an excellent tool for probing intergalactic medium (IGM) and for conducting cosmography. However, the DM contribution from the IGM (${\rm DM_{IGM}}$) depends on the fraction of baryon mass in the IGM, $f_{\rm IGM}$, which is not properly constrained. As $f_{\rm IGM}(z)$ is geometrically related to the Hubble parameter $H(z)$ and ${\rm DM_{IGM}}(z)$, here we propose that combining two independent measurements of FRBs and $H(z)$ in similar redshift ranges provides a novel and cosmology-free method to constrain the evolution of $f_{\rm IGM}(z)$. Under the assumption that $f_{\rm IGM}$ is evolving with redshift in a functional form, we forecast that the evolution of $f_{\rm IGM}(z)$ can be well inferred in a combined analysis of $\sim3000$ ${\rm DM_{IGM}}(z)$ derived from FRBs and $\sim50$ $H(z)$ derived from the Hubble parameter data. Though the efficiency of our method is not as good as that of the other model-independent method involving the joint measurements of DM and luminosity distance of FRBs, our method offers a new model-independent way to constrain $f_{\rm IGM}(z)$.

The tomographic Alcock-Paczynski (AP) method utilizes the redshift evolution of the AP distortion to place constraints on cosmological parameters. It has proved to be a robust method that can separate the AP signature from the redshift space distortion (RSD) effect, and deliver powerful cosmological constraints using the $\lesssim 40h^{-1}\ \rm Mpc$ clustering region. In previous works, the tomographic AP method was performed via the anisotropic 2-point correlation function statistic. In this work we consider the feasibility of conducting the analysis in the Fourier domain and examine the pros and cons of this approach. We use the integrated galaxy power spectrum (PS) as a function of direction, $\hat P_{\Delta k}(\mu)$, to quantify the magnitude of anisotropy in the large-scale structure clustering, and use its redshift variation to do the AP test. The method is tested on the large, high resolution Big-MultiDark Planck (BigMD) simulation at redshifts $z=0-1$, using the underlying true cosmology $\Omega_m=0.3071,\ w=-1$. Testing the redshift evolution of $\hat P_{\Delta k}(\mu)$ in the true cosmology and cosmologies deviating from the truth with $\delta \Omega_m=0.1,\ \delta w=0.3$, we find that the redshift evolution of the AP distortion overwhelms the effects created by the RSD by a factor of $\sim1.7-3.6$. We test the method in the range of $k\in(0.2,1.8)\ h\ \rm Mpc^{-1}$, and find that it works well throughout the entire regime. We tune the halo mass within the range $2\times 10^{13}$ to $10^{14}\ M_{\odot}$, and find that the change of halo bias results in $\lesssim 5 \%$ change in $\hat P_{\Delta k}(\mu)$, which is less significant compared with the cosmological effect. Our work shows that it is feasible to conduct the tomographic AP analysis in the Fourier space.

In Cassini ISS (Imaging Science Subsystem) images, contour detection is often performed on disk-resolved object to accurately locate their center. Thus, the contour detection is a key problem. Traditional edge detection methods, such as Canny and Roberts, often extract the contour with too much interior details and noise. Although the deep convolutional neural network has been applied successfully in many image tasks, such as classification and object detection, it needs more time and computer resources. In the paper, a contour detection algorithm based on H-ELM (Hierarchical Extreme Learning Machine) and DenseCRF (Dense Conditional Random Field) is proposed for Cassini ISS images. The experimental results show that this algorithm's performance is better than both traditional machine learning methods such as SVM, ELM and even deep convolutional neural network. And the extracted contour is closer to the actual contour. Moreover, it can be trained and tested quickly on the general configuration of PC, so can be applied to contour detection for Cassini ISS images.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

Using archival {\it Chandra} observations with a total effective exposure of 734 ks, we derive an updated catalog of point sources in the massive globular cluster Terzan 5. Our catalog covers an area of $58.1\, \rm arcmin^{2}$ ($R\leq 4.3 \, \rm arcmin$) with 489 X-ray sources, and more than $75\%$ of these sources are first detected in this cluster. We find significant dips in the radial distribution profiles of X-ray sources in Terzan 5, with the projected distance and width of the distribution dips for bright ($L_{X} \gtrsim 9.5\times 10^{30} {\rm\ erg\ \,s^{-1}}$) X-ray sources are larger than that of the faint ($L_{X} \lesssim 9.5\times 10^{30} {\rm\ erg\ \,s^{-1}}$) sources. By fitting the radial distribution of the X-ray sources with a"generalized King model", we estimated an average mass of $1.48\pm0.11\,M_{\odot}$ and $1.27\pm0.13\,M_{\odot}$ for the bright and faint X-ray sources, respectively. These results are in agreement with that observed in 47 Tuc, which may suggest a universal mass segregation effect for X-ray sources in GCs. Compared with 47 Tuc, we show that the two-body relaxation timescale of Terzan 5 is much smaller, but its dynamical age is significantly younger than 47 Tuc. These features suggest that the evolution of Terzan 5 is not purely driven by two-body relaxation, and tidal stripping effect also plays an important role in accelerating the dynamical evolution of this cluster.

Simultaneous $\gamma$-ray measurements of gamma-ray burst (GRB) spectra and polarization offer a unique way to determine the underlying emission mechanism(s) in these objects as well as probing the particle acceleration mechanism(s) that lead to the observed $\gamma$-ray emission. Herein we examine the jointly-observed data from POLAR and GBM of GRB 170114A to determine its spectral and polarization properties and seek to understand the emission processes that generate these observations. We aim to develop an extensible and statistically sound framework for these types of measurements applicable to other instruments. We leverage the existing 3ML analysis framework to develop a new analysis pipeline for simultaneously modeling the spectral and polarization data. We derive the proper Poisson likelihood for $\gamma$-ray polarization measurements in the presence of background. The developed framework is publicly available for similar measurements with other $\gamma$-ray polarimeters. The data are analyzed within a Bayesian probabilistic context and the spectral data from both instruments are simultaneously modeled with a physical, numerical synchrotron code. The spectral modeling of the data is consistent with a synchrotron photon model as has been found in a majority of similarly analyzed single-pulse GRBs. The polarization results reveal a slight trend of growing polarization in time reaching values of ~30% at the temporal peak of the emission. Additionally, it is observed that the polarization angle evolves with time throughout the emission. These results hint at a synchrotron origin of the emission but further observations of many GRBs are required to verify these evolutionary trends. Furthermore, we encourage the development of time-resolved polarization models for the prompt emission of GRBs as the current models are not predictive enough to enable a full modeling of our current data.

The Chinese Space Station Optical Survey (CSS-OS) is a planned full sky survey operated by the Chinese Space Station Telescope (CSST). It can simultaneously perform the photometric imaging and spectroscopic slitless surveys, and will probe weak and strong gravitational lensing, galaxy clustering, individual galaxies and galaxy clusters, active galactic nucleus (AGNs), and so on. It aims to explore the properties of dark matter and dark energy and other important cosmological problems. In this work, we focus on two main CSS-OS scientific goals, i.e. the weak gravitational lensing (WL) and galaxy clustering surveys. We generate the mock CSS-OS data based on the observational COSMOS and zCOSMOS catalogs. We investigate the constraints on the cosmological parameters from the CSS-OS using the Markov Chain Monte Carlo (MCMC) method. The intrinsic alignments, galaxy bias, velocity dispersion, and systematics from instrumental effects in the CSST WL and galaxy clustering surveys are also included, and their impacts on the constraint results are discussed. We find that the CSS-OS can improve the constraints on the cosmological parameters by a factor of a few (even one order of magnitude in the optimistic case), compared to the current WL and galaxy clustering surveys. The constraints can be further enhanced when performing joint analysis with the WL, galaxy clustering, and galaxy-galaxy lensing data. Therefore, the CSS-OS is expected to be a powerful survey for exploring the Universe. Since some assumptions may be still optimistic and simple, it is possible that the results from the real survey could be worse. We will study these issues in details with the help of simulations in the future.

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.

We investigate the effects of extended multi-year light curves (9-year photometry and 5-year spectroscopy) on the detection of time lags between the continuum variability and broad-line response of quasars at z>~1.5, and compare with the results using 4-year photometry+spectroscopy presented in a companion paper. We demonstrate the benefits of the extended light curves in three cases: (1) lags that are too long to be detected by the shorter-duration data but can be detected with the extended data; (2) lags that are recovered by the extended light curves but missed in the shorter-duration data due to insufficient light curve quality; and (3) lags for different broad line species in the same object. These examples demonstrate the importance of long-term monitoring for reverberation mapping to detect lags for luminous quasars at high-redshift, and the expected performance of the final dataset from the Sloan Digital Sky Survey Reverberation Mapping project that will have 11-year photometric and 7-year spectroscopic baselines.

We present our piggyback search for fast radio bursts using the Nanshan 26m Radio Telescope and the Kunming 40m Radio Telescope. The observations are performed in the L-band from 1380 MHz to 1700 MHz at Nanshan and S-band from 2170 MHz to 2310 MHz at Kunming. We built the \textsc{Roach2}-based FFT spectrometer and developed the real-time transient search software. We introduce a new radio interference mitigation technique named \emph{zero-DM matched filter} and give the formula of the signal-to-noise ratio loss in the transient search. Though we have no positive detection of bursts in about 1600 and 2400 hours data at Nanshan and Kunming respectively, an intriguing peryton was detected at Nanshan, from which hundreds of bursts were recorded. Perytons are terrestrial radio signals that mimic celestial fast radio bursts. They were first reported at Parkes and identified as microwave oven interferences later. The bursts detected at Nanshan show similar frequency swept emission and have double-peaked profiles. They appeared in different sky regions in about tens of minutes observations and the dispersion measure index is not exactly 2, which indicates the terrestrial origin. The peryton differs drastically from the known perytons detected at Parkes, because it appeared in a precise period of $p=1.71287\pm 0.00004$ s. Its origin remains unknown.

Here we present new ALMA observations of polarized dust emission from six of the most massive clumps in W43-Main. The clumps MM2, MM3, MM4, MM6, MM7, and MM8, have been resolved into two populations of fragmented filaments. From these two populations we extracted 81 cores (96 with the MM1 cores) with masses between 0.9 \Msun\ to 425 \Msun\ and a mass sensitivity of 0.08 M$_{\odot}$. The MM6, MM7, and MM8 clumps show significant fragmentation, but the polarized intensity appears to be sparse and compact. The MM2, MM3, and MM4 population shows less fragmentation, but with a single proto-stellar core dominating the emission at each clump. Also, the polarized intensity is more extended and significantly stronger in this population. From the polarized emission, we derived detailed magnetic field patterns throughout the filaments which we used to estimate field strengths for 4 out of the 6 clumps. The average field strengths estimations were found between 500 $\mu$G to 1.8 mG. Additionally, we detected and modeled infalling motions towards MM2 and MM3 from single dish HCO$^{+}(J=4 \rightarrow 3)$ and HCN$(J=4 \rightarrow 3)$ data resulting in mass infall rates of $\dot{\mathrm{M}}_{\mathrm{MM2}} = 1.2 \times 10^{-2}$ \Msun\ yr$^{-1}$ and $\dot{\mathrm{M}}_{\mathrm{MM3}} = 6.3 \times 10^{-3}$ \Msun\ yr$^{-1}$. By using our estimations, we evaluated the dynamical equilibrium of our cores by computing the total virial parameter $\alpha_{\mathrm{total}}$. For the cores with reliable field estimations, we found that 71\% of them appear to be gravitationally bound while the remaining 29\% are not. We concluded that these unbound cores, also less massive, are still accreting and have not yet reached a critical mass. This also implies different evolutionary time-scales, which essentially suggests that star-formation in high mass filaments is not uniform.

Fast radio bursts (FRBs) are highly dispersed and probably extragalactic radio flashes with millisecond-duration. Recently, the Canadian Hydrogen Intensity Mapping Experiment (using the CHIME/FRB instrument) has reported detections of 13 FRBs during a pre-commissioning phase. It is more exciting that one of the 13 FRBs is a second source of repeaters which suggests that CHIME/FRB and other wide-field sensitive radio telescopes will find a substantial population of repeating FRBs. We have proposed strongly lensed repeating FRBs as a precision cosmological probe, e.g. constraining the Hubble constant and model-independently estimating the cosmic curvature. Here, we study complementary constraints on the equation of state of dark energy from strongly lensed FRBs to currently available popular probes. It is found that, in the framework of Chevalier-Polarski-Linder parametrization, adding time delay measurement of 30 strongly lensed FRB systems to cosmic microwave background radiation and type Ia supernovae can improve the dark energy figure of merit by a factor 2. In the precision cosmology era, this improvement is of great significance for studying the nature of dark energy.

Previous X-ray observations toward the Nuclear Star Cluster (NSC) at the Galactic center have discovered thousands of point sources, most of which were believed to be cataclysmic variables (CVs), i.e., a white dwarf (WD) accreting from a low-mass companion. However, the population properties of these CVs remain unclear, which otherwise contain important information about the evolutionary history of the NSC. In this work we utilize ultradeep archival \textit{Chandra} observations to study the spectral properties of the NSC CVs, in close comparison with those in the Solar vicinity. We find that the NSC CVs have strong Fe XXV and Fe XXVI lines (both of which show equivalent widths $\sim200-300$ eV), indicating metal-rich companions. Moreover, their Fe XXVI to Fe XXV line flux ratio is used to diagnose the characteristic white dwarf mass ($M_{\rm WD}$) of NSC CVs. The results show that the CVs with $L_{\rm 2-10 keV}>6\times10^{31}$ erg s$^{-1}$ have a mean $M_{\rm WD}$ of $\sim0.6/1.0\,M_{\odot}$ if they are magnetic/non-magnetic CVs; while those with $L_{\rm 2-10 keV}$ between $1-6\times10^{31}$ erg s$^{-1}$ have a mean $M_{\rm WD}$ of $\sim0.8/1.2\,M_{\odot}$ if they are magnetic/non-magnetic CVs. All these \textit{Chandra}-detected CVs collectively contribute $\sim$30-50\% of the unresolved 20-40 keV X-ray emission from the NSC. The CV population with massive (i.e., $M_{\rm WD}\sim1.2M_{\odot}$) WDs have not been observed in the Solar vicinity or the Galactic bulge, and they might have been formed via dynamical encounters in the NSC.

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 Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

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.

We study the properties of hot beta-stable nuclear matter using equations of state derived within the Brueckner-Hartree-Fock approach at finite temperature including consistent three-body forces. Simple and accurate parametrizations of the finite-temperature equations of state are provided. The properties of hot neutron stars are then investigated within this framework, in particular the temperature dependence of the maximum mass. We find very small temperature effects and analyze the interplay of the different contributions.

Fast radio bursts (FRBs) are short duration ($\sim$millisecond) radio transients with cosmological origin. The simple sharp features of the FRB signal have been utilized to probe two fundamental laws of physics, namey, testing Einstein's weak equivalence principle and constraining the rest mass of the photon. Recently, \cite{hessels18} found that after correcting for dispersive delay, some of the bursts in FRB 121102 have complex time-frequency structures that include sub-pulses with a time-frequency downward drifting property. Using the delay time between sub-pulses in FRB 121102, here we show that the parameterized post-Newtonian parameter $\gamma$ is the same for photons with different energies to the level of $\left|\gamma_1-\gamma_2\right|<2.5\times 10^{-16}$, which is 1000 times better than previous constraints from FRBs using similar methods. We also obtain a stringent constraint on the photon mass, $m_{\gamma} < 5.1\times10^{-48}$ g, which is 10 times smaller than previous best limits on the photon mass derived through the velocity dispersion method.

We present the revised and high quality CCD time series observations in 2017 and 2018 of the W UMa-type contact binary BH~Cas. The combination of our new determinations of minimum times and those collected from the literatures reveal a long-term period increase with a rate of $\textit{dP/dt}\,=\,+3.27 \times 10^{-7} \textrm{d}\,\textrm{yr}^{-1}$, which can be accounted for by angular momentum transfer as the less-massive component loses its mass to the more massive one. One cyclical oscillation ($A'=0.00300$~days and $P'=20.09$~yrs) superimposed on the long-term period increase tendency is attributed to one additional object that orbits BH Cas. The low-resolution spectra of this source are observed at phase $\sim$0.0 and $\sim$0.5 in 2018 and 2019, and the spectral types are identified as K3V$\,\pm\,$1 and G8V$\,\pm\,$2. According to the luminosity contributions and spectra at these phases, the spectral type of the primary star is close to K3V$\,\pm\,$1. We infer a cool spot on the hotter and less-massive component to fit the asymmetry of light curves. In addition, changes in the location, temperature, and area of spots on the secondary star in different observations may indicate that the magnetic field activity of the secondary is more active than that of the primary. Such strong spot activities and optical flare are both detected at the same time from a W UMa-type star for the first time. With the \emph{XMM-Newton} observations taken in 2003 and 2012, the X-ray light curve is obtained and the X-ray spectra can be well described with a black-body or a bremsstrahlung model. The soft X-ray luminosity of BH Cas is estimated as 8.8\,--\,9.7$\,\times\,10^{30}$\,erg\,s$^{-1}$.

We derive dust masses ($M_{\rm dust}$) from the spectral energy distributions of 58 post-starburst galaxies (PSBs). There is an anticorrelation between specific dust mass ($M_{\rm dust}$/$M_{\star}$) and the time elapsed since the starburst ended, indicating that dust was either destroyed, expelled, or rendered undetectable over the $\sim$1 Gyr after the burst. The $M_{\rm dust}$/$M_{\star}$ depletion timescale, 205$^{+58}_{-37}$ Myr, is consistent with that of the CO-traced $M_{\rm H_2}/M_{\star}$, suggesting that dust and gas are altered via the same process. Extrapolating these trends leads to the $M_{\rm dust}/M_{\star}$ and $M_{\rm H_2}/M_{\star}$ values of early-type galaxies (ETGs) within 1-2 Gyr, a timescale consistent with the evolution of other PSB properties into ETGs. Comparing $M_{\rm dust}$ and $M_{\rm H_2}$ for PSBs yields a calibration, log $M_{\rm H_2}$ = 0.45 log $M_{\rm dust}$ + 6.02, that allows us to place 33 PSBs on the Kennicutt-Schmidt (KS) plane, $\Sigma \rm SFR-\Sigma M_{\rm H_2}$. Over the first $\sim$200-300 Myr, the PSBs evolve down and off of the KS relation, as their star formation rate (SFR) decreases more rapidly than $M_{\rm H_2}$. Afterwards, $M_{\rm H_2}$ continues to decline whereas the SFR levels off. These trends suggest that the star-formation efficiency bottoms out at 10$^{-11}\ \rm yr^{-1}$ and will rise to ETG levels within 0.5-1.1 Gyr afterwards. The SFR decline after the burst is likely due to the absence of gas denser than the CO-traced H$_2$. The mechanism of the $M_{\rm dust}/M_{\star}$ and$M_{\rm H_2}/M_{\star}$ decline, whose timescale suggests active galactic nucleus (AGN) or low-ionization nuclear emission-line region (LINER) feedback, may also be preventing the large CO-traced molecular gas reservoirs from collapsing and forming denser star forming clouds.

In this paper, we investigate the statistical signal-processing algorithm to measure the instant local clock jump from the timing data of multiple pulsars. Our algorithm is based on the framework of Bayesian statistics. In order to make the Bayesian algorithm applicable with limited computational resources, we dedicated our efforts to the analytic marginalization of irrelevant parameters. We found that the widely used parameter for pulsar timing systematics, the `Efac' parameter, can be analytically marginalized. This reduces the Gaussian likelihood to a function very similar to the Student's $t$-distribution. Our iterative method to solve the maximum likelihood estimator is also explained in the paper. Using pulsar timing data from the Yunnan Kunming 40m radio telescope, we demonstrate the application of the method, where 80-ns level precision for the clock jump can be achieved. Such a precision is comparable to that of current commercial time transferring service using satellites. We expect that the current method could help developing the autonomous pulsar time scale.

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

We present polarimetric observations of four Class II protoplanetary disks (DG Tau, Haro 6-13, RY Tau, and MWC 480) taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at 3 mm. The polarization morphologies observed fall into two distinct categories: azimuthal polarization (DG Tau and Haro 6-13) and polarization parallel to the disk minor axis (RY Tau and MWC 480). The mechanisms responsible for disk polarization at millimeter and submillimeter wavelengths are currently under debate. In this Letter, we investigate two mechanisms capable of producing polarized dust emission in disks: self-scattering and grain alignment to the radiation anisotropy. The polarization morphologies of DG Tau and Haro 6-13 are broadly consistent with that expected from radiation alignment (though radiative alignment still does not account for all of the features seen in these disks), while RY Tau and MWC 480 are more consistent with self-scattering. Such a variation in the polarized morphology may provide evidence of dust grain size differences between the sources.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns.

We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity and high resolution ($\sim 100$ au) observations show a clear hourglass magnetic field morphology centered on the protostellar system, but the central pattern is consistent with a toroidal field indicative of a circumstellar disk, although other interpretations are also possible, including field lines dragged by an equatorial accretion flow into a configuration parallel to the midplane. If a relatively large disk does exist, it would suggest that the magnetic braking catastrophe is averted in this system, not through a large misalignment between the magnetic and rotation axes, but rather through some other mechanisms, such as non-ideal magneto-hydrodynamic effects and/or turbulence. We have also found a relationship of decreasing polarization fractions with intensities and the various slopes of this relationship can be understood as multiple polarization mechanisms and/or depolarization from a changing field morphology. In addition, we found a prominent clumpy depolarization strip crossing the center perpendicular to the bipolar outflow. Moreover, a rough estimate of the magnetic field strength indicates that the field is strong enough to hinder formation of a rotationally supported disk, which is inconsistent with the feature of a central toroidal field.

The sensitivity of radial velocity (RV) surveys for exoplanet detection are extending to increasingly long orbital periods, where companions with periods of several years are now being regularly discovered. Companions with orbital periods that exceed the duration of the survey manifest in the data as an incomplete orbit or linear trend, a feature that can either present as the sole detectable companion to the host star, or as an additional signal overlain on the signatures of previously discovered companion(s). A diagnostic that can confirm or constrain scenarios in which the trend is caused by an unseen stellar, rather than planetary, companion is the use of high-contrast imaging observations. Here, we present RV data from the Anglo-Australian Planet Search (AAPS) for twenty stars that show evidence of orbiting companions. Of these, six companions have resolved orbits, with three that lie in the planetary regime. Two of these (HD~92987b and HD~221420b) are new discoveries. Follow-up observations using the Differential Speckle Survey Instrument (DSSI) on the Gemini South telescope revealed that five of the twenty monitored companions are likely stellar in nature. We use the sensitivity of the AAPS and DSSI data to place constraints on the mass of the companions for the remaining systems. Our analysis shows that a planetary-mass companion provides the most likely self-consistent explanation of the data for many of the remaining systems.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

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.

The inner disc of the local group galaxy M33 appears to be in settled rotational balance, and near IR images reveal a mild, large-scale, two-arm spiral pattern with no strong bar. We have constructed N-body models that match all the extensive observational data on the kinematics and surface density of stars and gas in the inner part of M33. We find that currently favoured models are unstable to the formation of a strong bar of semi-major axis 2 < a_B < 3 kpc within 1 Gyr, which changes the dynamical properties of the models to become inconsistent with the current, apparently well-settled, state. The formation of a bar is unaffected by how the gas component is modelled, by increasing the mass of the nuclear star cluster, or by making the dark matter halo counter-rotate, but it can be prevented by either reducing the mass-to-light ratio of the stars to Upsilon_V ~ 0.6 or Upsilon_K ~ 0.23 in solar units or by increasing the random motions of the stars. Also a shorter and weaker bar results when the halo is rigid and unresponsive. However, all three near-stable models support multi-arm spirals, and not the observed large-scale bi-symmetric spiral. A two-arm spiral pattern could perhaps be tidally induced, but such a model would require an implausibly low mass disc to avoid a bar and there is no visible culprit. Thus the survival of the current state of this exceptionally well-studied galaxy is not yet understood. We also suspect that many other unbarred galaxies present a similar puzzle.

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

It has been shown that time delay of strong gravitational lensing is not only an effective cosmological probe to constrain the Hubble constant, the matter density and the curvature of the universe but also useful for breaking the degeneracy of the dark energy equation of state and thus provide complementarity with other popular probes, such as type Ia supernovae, baryon acoustic oscillations, and cosmic microwave background radiation. Interestingly, compared to traditional strong lensing systems where quasars act as sources, strongly lensed gravitational waves (GWs) from compact binary coalescence and their electromagnetic (EM) counterparts systems have been recently proposed to be more powerful for studying cosmology since GWs and their EM counterparts are transients. Short durations of GWs and their EM counterparts are very advantageous to time delay measurement and lens profile modeling. Here, in the framework of Chevalier-Polarski-Linder (CPL) parametrization, we investigate improvement of constraining power on dark energy equation of state due to including time delay measurements of strong lensed GW systems. It is suggested that, on the basis of the third generation ground-based detector, e.g. Einstein Telescope, adding time delay of only 30 strong lensed GW systems to type Ia supernovae and cosmic microwave background radiation can improve the dark energy figure of merit by a factor 2. For the traditional standard siren method where the uncertainties of luminosity distances of GWs are $\sim10\%$, a few $\times~10^4$ events are expected to present similar constraints. In the precision cosmology era, this progress is of great significance for studying the nature of dark energy.

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

The excessive dispersion measure (DM) of fast radio bursts (FRBs) has been proposed to be a powerful tool to study intergalactic medium (IGM) and to perform cosmography. One issue is that the fraction of baryons in the IGM, $f_{\rm IGM}$, is not properly constrained. Here we propose a method of estimating $f_{\rm IGM}$ using a putative sample of FRBs with the measurements of both DM and luminosity distance $d_{\rm L}$. The latter can be obtained if the FRB is associated with a distance indicator (e.g. a gamma-ray burst or a gravitational wave event), or the redshift $z$ of the FRB is measured and $d_{\rm L}$ at the corresponding $z$ is available from other distance indicators (e.g. type Ia supernovae) at the same redshift. Since $d_{\rm L}/{\rm DM}$ essentially does not depend on cosmological parameters, our method can determine $f_{\rm IGM}$ independent of cosmological parameters. We parameterize $f_{\rm IGM}$ as a function of redshift and model the DM contribution from a host galaxy as a function of star formation rate. Assuming $f_{\rm IGM}$ has a mild evolution with redshift with a functional form and by means of Monte Carlo simulations, we show that an unbiased and cosmology-independent estimate of the present value of $f_{\rm IGM}$ with a $\sim 12\%$ uncertainty can be obtained with 50 joint measurements of $d_{\rm L}$ and DM. In addition, such a method can also lead to a measurement of the mean value of DM contributed from the local host galaxy.

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.

Knowledge of protostellar evolution has been revolutionized with the advent of surveys at near-infrared to submillimeter wavelengths. This has enabled the bolometric luminosities and bolometric temperatures (traditional protostellar evolution diagnostics) to be measured for large numbers of protostars. However, further progress is difficult without knowing the masses of the central protostars. Protostar masses can be most accurately determined via molecular line kinematics from millimeter interferometers (i.e., ALMA). Theoretical investigations have predicted the protostellar mass function (PMF) for various protostellar mass accretion models, and it is now imperative to observationally constrain its functional form. While ALMA has enabled protostellar mass measurements, samples approaching 100 sources are necessary to constrain the functional form of the PMF, and upgrades to ALMA and/or a new mm/cm facility will increase the feasibility of measuring such a large number of protostar masses. The masses of protostars will enable their stellar structure (radius and intrinsic luminosity), evolution, and accretion histories to be better understood. This is made more robust when effective temperatures and accretion rates can be measured via ground/space-based near to mid-infrared spectroscopy. Furthermore, access to supercomputing facilities is essential to fit the protostar masses via radiative transfer modeling and updated theoretical/numerical modeling of stellar structure may also be required.

We combine constraints from linear and nonlinear scales, for the first time, to study the interaction between dark matter and dark energy. We devise a novel N-body simulation pipeline for cosmological models beyond $\Lambda$CDM. This pipeline is fully self-consistent and opens a new window to study the nonlinear structure formation in general phenomenological interacting dark energy models. By comparing our simulation results with the SDSS galaxy-galaxy weak lensing measurements, we are able to constrain the strength of interaction between dark energy and dark matter. Compared with the previous studies using linear examinations, we point to plausible improvements on the constraints of interaction strength by using small scale information from weak lensing. This improvement is mostly due to the sensitivity of weak lensing measurements on nonlinear structure formation at low redshift. With this new pipeline, it is possible to look for smoking gun signatures of dark matter-dark energy interaction.

Gravitational wave (GW) echoes, if they exist, would be a probe to the near-horizon quantum structure of black hole (BH), which has motivated the searching for the echo signals in GW data. We point out that the echo phenomenology related with the potential structure might be not so simple as expected. In particular, if the near-horizon regime of BH is modelled as a multiple-barriers filter, the late-time GW ringdown waveform will exhibit the mixing of echoes, even the superpositions. As a result, the amplitudes of successive echoes might not drop sequentially.

Radial velocity (RV) searches for exoplanets have surveyed many of the nearest and brightest stars for long-term velocity variations indicative of a companion body. Such surveys often detect high-amplitude velocity signatures of objects that lie outside the planetary mass regime, most commonly those of a low-mass star. Such stellar companions are frequently discarded as false-alarms to the main science goals of the survey, but high-resolution imaging techniques can be employed to either directly detect or place significant constraints on the nature of the companion object. Here, we present the discovery of a compact companion to the nearby star HD~118475. Our Anglo-Australian Telescope (AAT) RV data allow the extraction of the full Keplerian orbit of the companion, found to have a minimum mass of 0.445~$M_\odot$. Follow-up speckle imaging observations at the predicted time of maximum angular separation rule out a main sequence star as the source of the RV signature at the 3.3$\sigma$ significance level, implying that the companion must be a low-luminosity compact object, most likely a white dwarf. We provide an isochrone analysis combined with our data that constrain the possible inclinations of the binary orbit. We discuss the eccentric orbit of the companion in the context of tidal circularization timescales and show that non-circular orbit was likely inherited from the progenitor. Finally, we emphasize the need for utilizing such an observation method to further understand the demographics of white dwarf companions around nearby stars.

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

Using archival {\it Chandra} observations with a total exposure of 510 ks, we present an updated catalog of point sources for Globular Cluster 47 Tucanae. Our study covers an area of $\sim 176.7$ arcmin$^{2}$ (i.e., with $R\lesssim7.5\arcmin$) with 537 X-ray sources. We show that the surface density distribution of X-ray sources in 47 Tuc is highly peaked in cluster center, rapidly decreases at intermediate radii, and finally rises again at larger radii, with two distribution dips at $R\sim 100\arcsec$ and $R\sim 170\arcsec$ for the faint ($L_{X}\lesssim 5.0\times 10^{30} {\rm\ erg\,s^{-1}}$) and bright ($L_{X}\gtrsim 5.0\times 10^{30} {\rm\ erg\,s^{-1}}$) groups of X-ray sources, separately. These distribution features are similar to those of Blue Straggler Stars (BSS), where the distribution dip is located at $R\sim 200\arcsec$ \citep{ferraro2004}. By fitting the radial distribution of each group of sources with a "generalized King model", we estimated an average mass of $1.51\pm0.17\ M_{\odot}$, $1.44\pm0.15\ M_{\odot}$ and $1.16\pm0.06\ M_{\odot}$ for the BSS, bright and faint X-ray sources, respectively. These results are consistent with the mass segregation effect of heavy objects in GCs, where more massive objects drop to the cluster center faster and their distribution dip propagates outward further. Besides, the peculiar distribution profiles of X-ray sources and BSS are also consistent with the mass segregation model of binaries in GCs, which suggests that in addition to the dynamical formation channel, primordial binaries are also a significant contributor to the X-ray source population in GCs.

Young stars accrete mass from a circumstellar disk, but at the same time disk and star eject outflows and jets. These outflows have an onion-like structure where the innermost and fastest layers are surrounded by increasingly lower velocity components. The outer layers are probably photo-evaporative and magnetocentrifugally launched disk winds, but the nature of the inner winds is still uncertain. Since the fastest components carry only a small fraction of the mass, they are best observed at high-energies (X-ray and UV) as the slower, more massive components do not reach plasma temperatures sufficient for relevant X-ray or UV emission. Outflows are the most likely way in which a star or its disk can shed angular momentum and allow accretion to proceed; thus we cannot understand the accretion and the rotation rate of young stars if we cannot solve the origin of the inner jet components. Stellar jets share characteristics with their counterparts in more massive astrophysical objects, such as stellar mass black holes and AGN, with the added benefit that young stars are found at much closer distances and thus scales not accessible in other types of objects can be resolved. To understand the origin and impact of the inner jets, sub-arcsecond imaging and spectroscopy in the UV and X-rays is required, together with theory and modelling to interpret existing and future observations.

Background: low-mass stars are the dominant product of the star formation process, and they trace star formation over the full range of environments, from isolated globules to clusters in the central molecular zone. In the past two decades, our understanding of the spatial distribution and properties of young low-mass stars and protostars has been revolutionized by sensitive space-based observations at X-ray and IR wavelengths. By surveying spatial scales from clusters to molecular clouds, these data provide robust measurements of key star formation properties. Goal: with their large numbers and their presence in diverse environments, censuses of low mass stars and protostars can be used to measure the dependence of star formation on environmental properties, such as the density and temperature of the natal gas, strengths of the magnetic and radiation fields, and the density of stars. Here we summarize how such censuses can answer three basic questions: i.) how is the star formation rate influenced by environment, ii.) does the IMF vary with environment, and iii.) how does the environment shape the formation of bound clusters? Answering these questions is an important step toward understanding star and cluster formation across the extreme range of environments found in the Universe. Requirements: sensitivity and angular resolution improvements will allow us to study the full range of environments found in the Milky Way. High spatial dynamic range (< 1arcsec to > 1degree scales) imaging with space-based telescopes at X-ray, mid-IR, and far-IR and ground-based facilities at near-IR and sub-mm wavelengths are needed to identify and characterize young stars.

The majority of the ultimate main-sequence mass of a star is assembled in the protostellar phase, where a forming star is embedded in an infalling envelope and encircled by a protoplanetary disk. Studying mass accretion in protostars is thus a key to understanding how stars gain their mass and ultimately how their disks and planets form and evolve. At this early stage, the dense envelope reprocesses most of the luminosity generated by accretion to far-infrared and submillimeter wavelengths. Time-domain photometry at these wavelengths is needed to probe the physics of accretion onto protostars, but variability studies have so far been limited, in large part because of the difficulty in accessing these wavelengths from the ground. We discuss the scientific progress that would be enabled with far-infrared and submillimeter programs to probe protostellar variability in the nearest kiloparsec.

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process.

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes.

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

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

Extremely low-mass white dwarfs (ELM WDs) are helium WDs with a mass less than $\sim$$0.3\rm\;M_\odot$. Most ELM WDs are found in double degenerates (DDs) in the ELM Survey led by Brown and Kilic. These systems are supposed to be significant gravitational-wave sources in the mHz frequency. In this paper, we firstly analyzed the observational characteristics of ELM WDs and found that there are two distinct groups in the ELM WD mass and orbital period plane, indicating two different formation scenarios of such objects, i.e. a stable Roche lobe overflow channel (RL channel) and common envelope ejection channel (CE channel). We then systematically investigated the formation of ELM WDs in DDs by a combination of detailed binary evolution calculation and binary population synthesis. Our study shows that the majority of ELM WDs with mass less than $0.22\rm\;M_\odot$ are formed from the RL channel. The most common progenitor mass in this way is in the range of $1.15-1.45\rm\;M_\odot$ and the resulting ELM WDs have a peak around $0.18\rm\;M_\odot$ when selection effects are taken into account, consistent with observations. The ELM WDs with a mass larger than $0.22\rm\;M_\odot$ are more likely to be from the CE channel and have a peak of ELM WD mass around $0.25\rm\;M_\odot$ which needs to be confirmed by future observations. By assuming a constant star formation rate of 2$\rm\;M_\odot yr^{-1}$ for a Milky Way-like galaxy, the birth rate and local density are $5\times10^{-4}\rm\;yr^{-1}$ and $1500\rm\;kpc^{-3}$, respectively, for DDs with an ELM WD mass less than $0.25\rm\;M_\odot$.

The complete three-dimensional structure of the magnetic field within molecular clouds has eluded determination despite its high value in determining controlling factors in the star formation process, as it cannot be directly probed observationally. Considering that inclination of the magnetic field relative to the plane of sky is one of the major sources of depolarization of thermal emission from dust in molecular clouds, we propose here a new method to estimate the inclination angle of the cloud-scale magnetic field based on the statistical properties of the observed polarization fraction. We test this method using a series of Monte Carlo experiments, and find that the method works well provided that deviations of magnetic field direction from the averaged values are small. When applied to synthetic observations of numerical simulations of star-forming clouds, our method gives fairly accurate measurements of the mean magnetic field inclination angle (within $10^\circ-25^\circ$), which can further be improved if we restrict our technique to regions of low dispersion in polarization angles ${\cal S}$. We tested our method on the BLASTPol polarimetric observations of the Vela C molecular cloud complex, which suggests that the magnetic field of Vela C has a high inclination angle ($\sim 60^\circ$), consistent with previous analyses.

Although many high-energy neutrinos detected by the IceCube telescope are believed to have anextraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedenteddiscovery of a high-energy muon neutrino event coincident with a multiwavelength flare from ablazar, TXS 0506+056, shed some light on the origin of the neutrinos. It is usually believed that ablazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH).Here we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclearinteractions between the accelerated cosmic-ray protons in the relativistic jet and the dense gasclouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in thejet, which could be much smaller than that required in the conventional hadronic model whichinstead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare fromthe optical to gamma-ray band can be well explained by invoking a second radiation zone in thejet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time thanthe multiwavelength flare so the neutrino event is not necessarily correlated with the flare but it is probably accompanied by a spectrum hardening above a few GeV.

Based on archival {\it Chandra} observations with a total exposure of 1.3 Ms, we study X-ray point sources in the Fornax cluster of galaxies, with the primary aim of searching for intra-cluster X-ray source populations. We detect 1177 point sources out to a projected radius of $\sim$30 arcmin ($\sim$180 kpc) from the cluster center and down to a limiting 0.5--8 keV luminosity of $\sim3\times10^{37}{\rm~erg~s^{-1}}$. We construct source surface density profile, after excluding sources associated with foreground stars, known globular clusters, ultra-compact dwarfs and galactic nuclei. From this profile we statistically identify $\sim$183 excess sources that are not associated with the bulk stellar content of the individual member galaxies of Fornax, nor with the cosmic X-ray background. Taking into account Poisson error and cosmic variance, the cumulative significance of this excess is at $\gtrsim 2\,\sigma$ level (with a maximum of 3.6\,$\sigma$) outside three effective radii of the central giant elliptical, NGC\,1399. The luminosity function of the excess sources is found to be significantly steeper than that of the GC-hosting sources (presumably low-mass X-ray binaries [LMXBs]), disfavoring the possibility that unidentified GCs are primarily responsible for the excess. We show that a large fraction of the excess can be related to the extended stellar halo of NGC\,1399 and/or the diffuse intra-cluster light, thus providing strong evidence for the presence of intra-cluster X-ray sources in Fornax, the second unambiguous case for a galaxy cluster after Virgo. Other possible origins of the excess, including supernova-kicked LMXBs and stripped nucleated dwarf galaxies are discussed.

In this paper we analyse the implications of the latest cosmological data sets to test the Etherington's distance duality relation (DDR), which connects the luminosity distance $D_L$ and angular diameter distance $D_A$ at the same redshift. For $D_L$ we consider the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), which can be considered as standard candles (or standard siren), while the angular diameter distances $D_A$ are derived from the newly-compiled sample of compact radio quasars observed by very-long-baseline interferometry (VLBI), which represents a type of new cosmological standard ruler. Alleviating the absorption and scattering effects of dust in the Universe, this will create a valuable opportunity to directly test DDR at much higher precision with the combination of gravitational wave (GW) and electromagnetic (EM) signals. Our results show that, with the combination of the current radio quasar observations, the duality-distance relation can be verified at the precision of $10^{-2}$. Moreover, the Einstein Telescope ET would produce more robust constraints on the validity of such distance duality relation (at the precision of $10^{-3}$), with a larger sample of compact milliarcsecond radio quasars detected in future VLBI surveys.