In an effort to probe the origin of surface brightness profile (SBP) breaks widely observed in nearby disk galaxies, we carry out a comparative study of stellar population profiles of 635 disk galaxies selected from the MaNGA spectroscopic survey. We classify our galaxies into single exponential (TI), down-bending (TII) and up-bending (TIII) SBP types, and derive their spin parameters and radial profiles of age/metallicity-sensitive spectral features. Most TII (TIII) galaxies have down-bending (up-bending) star formation rate (SFR) radial profiles, implying that abrupt radial changes of SFR intensities contribute to the formation of both TII and TIII breaks. Nevertheless, a comparison between our galaxies and simulations suggests that stellar migration plays a significant role in weakening down-bending $\Sigma_{\star}$ profile breaks. While there is a correlation between the break strengths of SBPs and age/metallicity-sensitive spectral features for TII galaxies, no such correlation is found for TIII galaxies, indicating that stellar migration may not play a major role in shaping TIII breaks, as is evidenced by a good correspondence between break strengths of $\Sigma_{\star}$ and surface brightness profiles of TIII galaxies. We do not find evidence for galaxy spin being a relevant parameter for forming different SBP types, nor do we find significant differences between the asymmetries of galaxies with different SBP types, suggesting that environmental disturbances or satellite accretion in the recent past do not significantly influence the break formation. By dividing our sample into early and late morphological types, we find that galaxies with different SBP types follow nearly the same tight stellar mass-$R_{25}$ relation, which makes the hypothesis that stellar migration alone can transform SBP types from TII to TI and then to TIII highly unlikely.

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 propose a new non-thermal mechanism of dark matter production based on vacuum misalignment. A global $X$-charge asymmetry is generated at high temperatures, under which both the will-be Higgs and the dark matter are charged. At lower energies, the vacuum changes alignment and breaks the $U(1)_X$, leading to the emergence of the Higgs and of a fraction of charge asymmetry stored in the stable dark matter relic. This mechanism can be present in a wide variety of models based on vacuum misalignment, and we demonstrate it in a composite Higgs template model, where all the necessary ingredients are naturally present. A light pseudo-scalar $\eta$ is always predicted, with interesting implications for cosmology, future supernova observations and exotic $Z \to \gamma \eta$ decays.

Compact steep-spectrum sources (CSSs) likely represent a population of young radio-loud active galactic nuclei (AGNs) and have been identified as gamma-ray emitting sources. We present a comprehensive analysis of their gamma-ray emission observed with Fermi/LAT and establish their broadband spectral energy distributions (SEDs). We derive their jet properties by the SED fits with a two-zone leptonic model for radiations from the compact core and large-scale extended region, and explore the possible signature of a unification picture of jet radiation among subclasses of AGNs. We show that the observed gamma-rays of CSSs with significant variability are contributed by the radiation of their compact cores via the inverse Compton process of the torus photons. The derived power-law distribution index of the radiating electrons is p_1~1.5-1.8, magnetic field strength is B~0.15-0.6 G, and Doppler boosting factor is \delta~2.8-8.9. Assuming that the jet is composed of electron-positron pairs, the compact cores of CSSs are magnetized and have a high radiation efficiency, similar to that of flat spectrum radio quasars. The six CSSs on average have higher Eddington ratio and black hole mass than those non-GeV-detected CSSs, and they follow the correlation between the jet power in units of Eddington luminosity (P_jet/L_Edd) and Eddington ratio (R_Edd) with other sub-classes of AGNs, P_jet/L_Edd~R_Edd^0.52, indicating that R_Edd would be a key physical driver for the unification scheme of AGN jet radiation.

Owing to the detection of an extremely bright fast radio burst (FRB) 200428 associated with a hard X-ray counterpart from the magnetar soft gamma-ray repeater (SGR) 1935+2154, the distance of SGR 1935+2154 potentially hosted in the supernova remnant (SNR) G57.2+0.8 can be revisited. Under the assumption that the SGR and the SNR are physically related, in this Letter, by investigating the dispersion measure (DM) of the FRB contributed by the foreground medium of our Galaxy and the local environments and combining with other observational constraints, we find that the distance of SGR 1935+2154 turns out to be $9.0\pm2.5\,$kpc and the SNR radius falls into $10$ to $18\,$pc since the local DM contribution is as low as $0-18\,$pc cm$^{-3}$. These results are basically consistent with the previous studies. In addition, an estimate for the Faraday rotation measure of the SGR and SNR is also carried out.

Star-forming galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into $\gamma$-rays by colliding with the interstellar medium (ISM). Several nearby star-forming galaxies have been detected in GeV-TeV $\gamma$-rays. It is also found that the $\gamma$-ray luminosities in 0.1-100 GeV correlate well with indicators of star formation rates of the galaxies, such as the total infrared (IR) luminosity. In this paper, we report a systematic search for possible $\gamma$-ray emission from galaxies in the IRAS Revised Bright Galaxies Sample, using 11.4 years of $\gamma$-ray data taken by the Fermi Large Area Telescope (LAT). Two new galaxies, M33 and Arp 299, are detected significantly. The two galaxies are consistent with the empirical correlation between the $\gamma$-ray luminosity and total infrared luminosity, suggesting that their $\gamma$-ray emissions should mainly originate from CRs interacting with ISM. Nevertheless, there is a tentative evidence that the flux of the $\gamma$-ray emission from Arp~299 is variable. If the variability is true, part of the emission from Arp 299 should originate from the obscured AGN in this interacting galaxy system. In addition, we find that the $\gamma$-ray excess from M33 is located at the northeast region of the galaxy, where a supergiant H II region, NGC604, resides. This indicates that some bright star-forming regions in spiral galaxies could play a dominant role in the galaxy in producing $\gamma$-ray emission.

AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose formation, nature and role in galaxy evolution are yet to be understood. In this work we first present a new automated method based on deep neural networks to detect clumps in galaxy images and show that it is more sensitive and faster than previous proposed methods. We then use this method to systematically detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA forward modeled zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair and direct comparisons between observations and simulations. After correcting for these effects with a simple mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. We find that between $20-40\%$ of $\sim10^{10}$ solar mass star-forming galaxies contain at least one massive ($M_c>10^{7}M_\odot$) off-centered clump but only $\sim2-5\%$ of the total galaxy stellar mass is in those clumps. We also show indications that the contribution of clumps to the stellar mass is more important in extended and low mass galaxies. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function when confronted under the same conditions although tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ through violent disk instabilities.

We report measurements of the sky-projected spin-orbit angle for AU\,Mic\,b, a Neptune-size planet orbiting a very young ($\sim20$\,Myr) nearby pre-main sequence M dwarf star which also hosts a bright, edge-on, debris disk. The planet was recently discovered from preliminary analysis of radial velocity observations and confirmed to be transiting its host star from photometric data from the NASA's \textit{TESS} mission. We obtained radial velocity measurements of the Rossiter--McLaughlin effect for AU\,Mic\,b over the course of two partially observable transits and one full transit from high-resolution spectroscopic observations made with the {\textsc{Minerva}}-Australis telescope array. We find that the planet's orbit is aligned with the spin-axis of its host star with a sky-projected spin-orbit angle of $\lambda = 5{^{+16}_{-15}}^{\circ}$. This measurement is in excellent agreement with observations carried out on other instruments around the same time. AU\,Mic is the youngest exoplanetary system for which the projected spin-orbit angle has been measured, making it a key data point in the study of the formation and migration of exoplanets -- particularly given that the system is also host to a bright debris disk. Given that we find this system (star, planet, and disk) to be in near spin-orbit alignment, it seems likely that the planet formed farther out beyond the ice-line and then migrated inwards to its current location via quiescent migration through the disk, rather than through mechanisms that would excite the inclination of the planet's orbit.

With the published data of apparent axis ratios for 1109 ultra-diffuse galaxies (UDGs) located in 17 low-redshift (z~ 0.020 - 0.063) galaxy clusters and 84 UDGs in 2 intermediate-redshift (z~ 0.308 - 0.348) clusters, we take advantage of a Markov Chain Monte Carlo approach and assume a ubiquitous triaxial model to investigate the intrinsic morphologies of UDGs. In contrast to the conclusion of Burkert (2017), i.e., the underlying shapes of UDGs are purely prolate ($C=B<A$), we find that the data favor the oblate-triaxial models ($C<B\lesssim A$) over the nearly prolate ones. We also find that the intrinsic morphologies of UDGs are relevant to their stellar masses/luminosities, environments, and redshifts. First, for the low-redshift UDGs in the same environment, the more-luminous ones are always thicker than the less-luminous counterparts, possibly due to the more voilent internal supernovae feedback or external tidal interactions for the progenitors of the more-luminous UDGs. The UDG thickness dependence on luminosity is distinct from that of the typical quiescent dwarf ellipticals (dEs) and dwarf spheroidals (dSphs) in the local clusters and groups, but resembles that of massive galaxies; in this sense, UDGs may not be simply treated as an extension of the dE/dSph class with similar evolutionary histories. Second, for the low-redshift UDGs within the same luminosity range, the ones with smaller cluster-centric distances are more puffed-up, probably attributed to tidal interactions. Finally, the intermediate-redshift cluster UDGs are more flattened, which plausibly suggests a `disky' origin for high-redshift, initial UDGs.

The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.

The tilt angle, current helicity and twist of solar magnetic fields can be observed in solar active regions. We carried out estimates of these parameters by two ways. Firstly, we consider the model of turbulent convective cells (super-granules) which have a loop floating structure towards the surface of the Sun. Their helical properties are attained during the rising process in the rotating stratified convective zone. The other estimate is obtained from a simple mean-field dynamo model that accounts magnetic helicity conservation. The both values are shown to be capable to give important contributions to the observable tilt, helicity and twist.

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the \startrack\ binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for \gaia\ and \lamost\ surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$--$340$ nBHBs to be observed by \gaia, although the numbers may drop to $\sim10$--$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For \lamost\ we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy.

We present $\gtrsim 15,000$ metal-rich (${\rm [Fe/H]}>-0.2$dex) A and F stars whose surface abundances deviate strongly from Solar abundance ratios and cannot plausibly reflect their birth material composition. These stars are identified by their high [Ba/Fe] abundance ratios (${\rm [Ba/Fe]}>1.0$dex) in the LAMOST DR5 spectra analyzed by Xiang et al. (2019). They are almost exclusively main sequence and subgiant stars with $T_{\rm eff}\gtrsim6300$K. Their distribution in the Kiel diagram ($T_{\rm eff}$--$\log g$) traces a sharp border at low temperatures along a roughly fixed-mass trajectory (around $1.4M_\odot)$ that corresponds to an upper limit in convective envelope mass fraction of around $10^{-4}$. Most of these stars exhibit distinctly enhanced abundances of iron-peak elements (Cr, Mn, Fe, Ni) but depleted abundances of Mg and Ca. Rotational velocity measurements from GALAH DR2 show that the majority of these stars rotate slower than typical stars in an equivalent temperature range. These characteristics suggest that they are related to the so-called Am/Fm stars. Their abundance patterns are qualitatively consistent with the predictions of stellar evolution models that incorporate radiative acceleration, suggesting they are a consequence of stellar internal evolution particularly involving the competition between gravitational settling and radiative acceleration. These peculiar stars constitute 40% of the whole population of stars with mass above 1.5$M_\odot$, affirming that "peculiar" photospheric abundances due to stellar evolution effects are a ubiquitous phenomenon for these intermediate-mass stars. This large sample of Ba-enhanced chemically peculiar A/F stars with individual element abundances provides the statistics to test more stringently the mechanisms that alter the surface abundances in stars with radiative envelopes.

We present a sample of $\sim$ 140,000 primary red clump (RC) stars of spectral signal-to-noise ratios higher than 20 from the LAMOST Galactic spectroscopic surveys, selected based on their positions in the metallicity-dependent effective temperature--surface gravity and color--metallicity diagrams, supervised by high-quality $Kepler$ asteroseismology data. The stellar masses and ages of those stars are further determined from the LAMOST spectra, using the Kernel Principal Component Analysis method, trained with thousands of RCs in the LAMOST-$Kepler$ fields with accurate asteroseismic mass measurements. The purity and completeness of our primary RC sample are generally higher than 80 per cent. For the mass and age, a variety of tests show typical uncertainties of 15 and 30 per cent, respectively. Using over ten thousand primary RCs with accurate distance measurements from the parallaxes of Gaia DR2, we re-calibrate the $K_{\rm s}$ absolute magnitudes of primary RCs by, for the first time, considering both the metallicity and age dependencies. With the the new calibration, distances are derived for all the primary RCs, with a typical uncertainty of 5--10 per cent, even better than the values yielded by the Gaia parallax measurements for stars beyond 3--4 kpc. The sample covers a significant volume of the Galactic disk of $4 \leq R \leq 16$ kpc, $|Z| \leq 5$ kpc, and $-20 \leq \phi \leq 50^{\circ}$. Stellar atmospheric parameters, line-of-sight velocities and elemental abundances derived from the LAMOST spectra and proper motions of Gaia DR2 are also provided for the sample stars. Finally, the selection function of the sample is carefully evaluated in the color-magnitude plane for different sky areas. The sample is publicly available.

RV variable stars are important in astrophysics. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic survey has provided ~ 6.5 million stellar spectra in its Data Release 4 (DR4). During the survey, ~ 4.7 million unique sources were targeted and ~ 1 million stars observed repeatedly. The probabilities of stars being RV variables are estimated by comparing the observed radial velocity variations with the simulated ones. We build a catalog of 80,702 RV variable candidates with probability greater than 0.60 by analyzing the duplicate-observed multi-epoch sources covered by the LAMOST DR4. Simulations and cross-identifications show that the purity of the catalog is higher than 80%. The catalog consists of 77% binary systems and 7% pulsating stars as well as 16% pollution by single stars. 3,138 RV variables are classified through cross-identifications with published results in literatures. By using the 3,138 sources common to both LAMOST and a collection of published RV variable catalogs we are able to analyze LAMOST's RV variable detection rate. The efficiency of the method adopted in this work relies not only on the sampling frequency of observations but also periods and amplitudes of RV variables. With the progress of LAMOST, Gaia and other surveys, more and more RV variables would will be confirmed and classified. This catalog is valuable for other large-scale surveys, especially for RV variable searches. The catalog will be released according to the LAMOST Data Policy via http://dr4.lamost.org.

We report the detection of high-energy gamma-ray signal towards the young star-forming region, W40. Using 10-year Pass 8 data from the Fermi Large Area Telescope (Fermi-LAT), we extracted an extended gamma-ray excess region with a significance of about 18sigma. The radiation has a spectrum with a photon index of 2.49 +/- 0.01. The spatial correlation with the ionized gas content favors the hadronic origin of the gamma-ray emission. The total cosmic-ray (CR) proton energy in the gamma-ray production region is estimated to be the order of 10^47 erg. However, this could be a small fraction of the total energy released in cosmic rays (CRs) by local accelerators, presumably by massive stars, over the lifetime of the system. If so, W40, together with earlier detections of gamma-rays from Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and 30 Dor C, supports the hypothesis that young star clusters are effective CR factories. The unique aspect of this result is that the gamma-ray emission is detected, for the first time, from a stellar cluster itself, rather than from the surrounding "cocoons".

This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a single electron-hole pair. The energy spectrum is reported from a blind analysis with 1.2 gram-days of exposure acquired in an above-ground laboratory. With charge carrier trapping and impact ionization effects incorporated into the dark matter signal models, the dark matter-electron cross section $\bar{\sigma}_{e}$ is constrained for dark matter masses from 0.5--$10^{4} $MeV$/c^{2}$; in the mass range from 1.2--50 eV$/c^{2}$ the dark photon kinetic mixing parameter $\varepsilon$ and the axioelectric coupling constant $g_{ae}$ are constrained. The minimum 90% confidence-level upper limits within the above mentioned mass ranges are $\bar{\sigma}_{e}\,=\,8.7\times10^{-34}$ cm$^{2}$, $\varepsilon\,=\,3.3\times10^{-14}$, and $g_{ae}\,=\,1.0\times10^{-9}$.

We revisit the diagonal ridge feature (diagonal distributions in the $R, v_{\phi}$ plane) found in $Gaia$ and present timing analysis for it between Galactocentric distances of $R=7.5$ and 12 \,kpc, using Main-Sequence-Turn-Off and OB stars selected from the LAMOST Galactic spectroscopic surveys. We recover the ridge pattern in the $R$-$v_{\phi}$ plane color coded by stellar number density and mean radial velocity and find this feature is presented from very young (OB stars, few hundred \,Myr) to very old populations ($\tau > $ 9 \,Gyr), accompanying with the north-south asymmetrical ridges. One of the ridge patterns, with constant angular momentum per unit mass, shows clear variations with different age populations compared. However, the remaining two ones are very stable, clearly implying there might have two kinds of ridge patterns with different dynamical origins and evolution. For the vertical velocity and height, there are no ridge features corresponding to the in-plane Galactoseismogy features, but the metallicity and [$\alpha/Fe$] have some weak signals. This first temporal chemo-kinematical analysis for ridge strongly supports that the ridge pattern might be originated from the phase mixing of disrupting spiral arms implying different ridges have different physical processes and no coupling phenomenon, although other possible scenarios such as Sgr perturbations could not be ruled out.

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

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.

In the era of big data astronomy, next generation telescopes and large sky surveys produce data sets at the TB or even PB level. Due to their large data volumes, these astronomical data sets are extremely difficult to transfer and analyze using personal computers or small clusters. In order to offer better access to data, data centers now generally provide online science platforms that enable analysis close to the data. The Chinese Virtual Observatory (China-VO) is one of the member projects in the International Virtual Observatory Alliance and it is dedicated to providing a research and education environment where globally distributed astronomy archives are simple to find, access, and interoperate. In this study, we summarize highlights of the work conducted at the China-VO, as well the experiences and lessons learned during the full life-cycle management of astronomical data. Finally, We discuss the challenges and future trends for astronomical science platforms.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

We present a study on the radial profiles of D4000,luminosity-weighted stellar ages $\tau_L$,and luminosity-weighted stellar metallicities $[Z/H]_L$ of 3654 nearby galaxies($0.01<z<0.15$)using the IFU spectroscopic data from the MaNGA survey available in the SDSS DR15,in an effort to explore the connection between median stellar population radial gradients($\nabla$D4000,$\nabla\tau_L,\nabla[Z/H]_L$)out to~$1.5R_e$ and various galaxy properties,including stellar mass($M_\star$),specific star formation rate(sSFR),morphologies,and local environment. We find that $M_\star$ is the single most predictive physical property for$\nabla$D4000 and$\nabla[Z/H]_L$. The most predictive properties for $\nabla\tau_L$ are sSFR,and to a lesser degree,$M_\star$. The environmental parameters,including local galaxy overdensities and central-satellite division,have virtually no correlation with stellar population radial profiles for the whole sample,but the $\nabla$D4000 of star-forming satellite galaxies with$M_\star\lesssim 10^{10}M_\odot$exhibit a significant positive correlation with galaxy overdensities. Galaxies with lower sSFR have on average steeper negative stellar population gradients,and this sSFR dependence is stronger for more massive star-forming galaxies. The negative correlation between the median stellar population gradients and$M_\star$ are best described largely as segmented relationships, whereby median gradients of galaxies with$\log M_\star\lesssim 10$(with the exact value depending on sSFR)have much weaker mass dependence than galaxies with higher$M_\star$. While the dependence of the radial gradients of ages and metallicities on T-Types and central stellar mass surface densities are generally not significant,galaxies with later T-Types or lower central mass densities tend to have significantly lower D4000,younger$\tau_L$ and lower$[Z/H]_L$ across the radial ranges probed in this study.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $\gamma ({\rm age})$ decreasing with age and its first derivative $\dot{\gamma} ({\rm age})$ is different from zero. The azimuth of line of nodes $\phi _w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.

Clues to the formation and evolution of Nuclear Star Clusters (NSCs) lie in their stellar populations. However, these structures are often very faint compared to their host galaxy, and spectroscopic analysis of NSCs is hampered by contamination of light from the rest of the system. With the introduction of wide-field IFU spectrographs, new techniques have been developed to model the light from different components within galaxies, making it possible to cleanly extract the spectra of the NSCs and study their properties with minimal contamination from the light of the rest of the galaxy. This work presents the analysis of the NSCs in a sample of 12 dwarf galaxies in the Fornax Cluster observed with MUSE. Analysis of the stellar populations and star-formation histories reveal that all the NSCs show evidence of multiple episodes of star formation, indicating that they have built up their mass further since their initial formation. The NSCs were found to have systematically lower metallicities than their host galaxies, which is consistent with a scenario for mass-assembly through mergers with infalling globular clusters, while the presence of younger stellar populations and gas emission in the core of two galaxies is indicative of in-situ star formation. We conclude that the NSCs in these dwarf galaxies likely originated as globular clusters that migrated to the core of the galaxy which have built up their mass mainly through mergers with other infalling clusters, with gas-inflow leading to in-situ star formation playing a secondary role.

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

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$\alpha$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$\alpha$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

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.

We present the morphology and stellar population of 27 extremely metal poor galaxies (EMPGs) at $z\sim0$ with metallicities of $0.01-0.1$ Z$_{\odot}$. We conduct multi-component surface brightness (SB) profile fitting for the deep Subaru/HSC $i$-band images of the EMPGs with the Galfit software, carefully removing the SB contributions of potential associated galaxies (PAGs). We find that the EMPGs with a median stellar mass of $\log(M_{*}/{\rm M}_{\odot})=5.77$ have a median S\'ersic index of $n=1.08$ and a median effective radius of $r_{\rm e}=176$ pc, suggesting that typical EMPGs are very compact disk galaxies. We compare the EMPGs with $z\sim 6$ galaxies and local galaxies on the size-mass ($r_{\rm e}$--$M_*$) diagram, and identify that the majority of the EMPGs have a $r_{\rm e}$--$M_*$ relation similar to $z\sim0$ star-forming galaxies rather than $z\sim6$ galaxies. Not every EMPG is a local analog of high-$z$ young galaxies in the $r_{\rm e}$--$M_*$ relation. We also study the PAGs of our galaxies, and find that 23 out of the 27 EMPGs show detectable PAGs within the projected distance of 10 kpc. The PAGs have median values of $n=0.93$, $r_{\rm e}=1.41$ kpc, and $\log(M_{*}/{\rm M}_{\odot})=7.47$ that are similar to those of local dwarf irregulars and ultra-diffuse galaxies. A spectrum of one EMPG-PAG system, so far available, indicates that the PAG is dynamically related to the EMPG with a median velocity difference of $\Delta V=96.1$ km s$^{-1}$. This moderately-large $\Delta V$ cannot be explained by the dynamics of the PAG, but likely by the infall on the PAG. EMPGs may form in infalling gas, and become part of PAG disks undergoing further metal enrichment.

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 have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates from the helical motion jet in a supermassive binary black hole (SMBBH) system undergoing major merger, we estimate the primary black hole mass M~7.7*10^{9} M_sun. To explore the origin of the gamma-ray, we investigated the optical-gamma-ray correlations using discrete correlation function (DCF) method, and found that the correlation between the two bands is very significant. This strong correlation tends to imply lepton self-synchro-Compton (LSSC) model to produce the gamma-ray.

The Search for Extraterrestrial Intelligence (SETI) attempts to address the possibility of the presence of technological civilizations beyond the Earth. Benefiting from high sensitivity, large sky coverage, an innovative feed cabin for China's Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed the SETI first observations with FAST's newly commisioned 19-beam receiver; we report preliminary results in this paper. Using the data stream produced by the SERENDIP VI realtime multibeam SETI spectrometer installed at FAST, as well as its off-line data processing pipelines, we identify and remove four kinds of radio frequency interference(RFI): zone, broadband, multi-beam, and drifting, utilizing the Nebula SETI software pipeline combined with machine learning algorithms. After RFI mitigation, the Nebula pipeline identifies and ranks interesting narrow band candidate ET signals, scoring candidates by the number of times candidate signals have been seen at roughly the same sky position and same frequency, signal strength, proximity to a nearby star or object of interest, along with several other scoring criteria. We show four example candidates groups that demonstrate these RFI mitigation and candidate selection. This preliminary testing on FAST data helps to validate our SETI instrumentation techniques as well as our data processing pipeline.

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.

Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at 43 GHz) fluxes in these sources. The gamma-ray and radio fluxes of 3C 454.3 and PKS 1633+382 are correlated without any time-lag, suggesting that they are from the same radiation region. Similar correlation is also observed in 3C 273 and PKS 1222+216, but the radio flux is lag behind the gamma-ray flux approximately ~160 days and ~290 days, respectively. This likely suggests that their gamma-ray emission regions are located at the upstream of their radio cores at 43 GHz. The gamma-ray and radio fluxes of the other four blazars are not correlated, implying that the gamma-ray and radio emission may be from different regions in their jets. The gamma-ray light curves of the eight blazars can be decomposed into some long timescale variability components and fast spike flares. We propose that they may be attributed to the central engine activity and the magnetic reconnection process or turbulence in the local emission region, respectively.

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored.

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

We continue our empirical study of the emission line flux originating in the cool ($T\sim10^4$ K) gas that populates the halos of galaxies and their environments. Specifically, we present results obtained for a sample of galaxy pairs with a range of projected separations, {\bf $10 < {S_p/\rm kpc} < 200$}, and mass ratios $<$ 1:5, intersected by 5,443 SDSS lines of sight at projected radii of 10 to 50 kpc from either or both of the two galaxies. We find significant enhancement in H$\alpha$ emission and a moderate enhancement in [N {\small II}]6583 emission for low mass pairs (mean stellar mass per galaxy, $\overline{\rm M}_*, <10^{10.4} {\rm M}_\odot$) relative to the results from a control sample. This enhanced H$\alpha$ emission comes almost entirely from sight lines located between the galaxies, consistent with a short-term, interaction-driven origin for the enhancement. We find no enhancement in H$\alpha$ emission, but significant enhancement in [N {\small II}]6583 emission for high mass ($\overline{\rm M}_* >10^{10.4}{\rm M}_\odot$) pairs. Furthermore, we find a dependence of the emission line properties on the galaxy pair mass ratio such that those with a mass ratio below 1:2.5 have enhanced [N {\small II}]6583 and those with a mass ratio between 1:2.5 and 1:5 do not. In all cases, departures from the control sample are only detected for close pairs ($S_p <$ 100 kpc). Attributing an elevated [N {\small II}]6583/H$\alpha$ ratio to shocks, we infer that shocks play a role in determining the CGM properties for close pairs that are among the more massive and have mass ratios closer to 1:1.

Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate nonthermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anti-correlated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the period of the QPOs is associated to the kink growth time scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities.

The absorption features of olivine in visible and near-infrared (VNIR) reflectance spectra are the key spectral parameters in its mineralogical studies. Generally, these spectral parameters can be obtained by exploiting the Modified Gaussian Model (MGM) with a proper continuum removal. However, different continua may change the deconvolution results of these parameters. This paper investigates the diagnostic spectral features of olivine with diverse chemical compositions. Four different continuum removal methods with MGM for getting the deconvolution results are presented and the regression equations for predicting the Mg-number (Fo#) are introduced. The results show that different continua superimposed on the mineral absorption features will make the absorption center shift, as well as the obvious alterations in shape, width, and strength of the absorption band. Additionally, it is also found that the logarithm of a second-order polynomial continuum can match the overall shape of the spectrum in logarithmic space, and the improved regression equations applied to estimate the chemical composition of olivine-dominated spectra also have a better performance. As an application example, the improved approach is applied to pulse laser irradiated olivine grains to simulate and study the space weathering effects on olivine diagnostic spectral features. The experiments confirm that space weathering can make the absorption band center shift toward longer wavelength. Therefore, the Fo# estimated from remote sensing spectra may be less than its actual chemical composition. These results may provide valuable information for revealing the difference between the spectra of olivine grains and olivine-dominated asteroids.

It has long been speculated that many starburst or compact dwarf galaxies are resulted from dwarf-dwarf galaxy merging, but unequivocal evidence for this possibility has rarely been reported in the literature. We present the first study of deep optical broadband images of a gas-dominated blue compact dwarf galaxy (BCD) VCC848 (Mstar=2e8Msun) which hosts extended stellar shells and thus is confirmed to be a dwarf-dwarf merger. VCC848 is located in the outskirts of the Virgo Cluster. By analyzing the stellar light distribution, we found that VCC848 is the result of a merging between two dwarf galaxies with a primary-to-secondary mass ratio < ~ 5 for the stellar components and < ~ 2 for the presumed dark matter halos. The secondary progenitor galaxy has been almost entirely disrupted. The age-mass distribution of photometrically selected star cluster candidates in VCC848 implies that the cluster formation rate (CFR, proportional to star formation rate) was enhanced by a factor of ~ 7 - 10 during the past 1 Gyr. The merging-induced enhancement of CFR peaked near the galactic center a few hundred Myr ago and has started declining in the last few tens of Myr. The current star formation activities, as traced by the youngest clusters, mainly occur at large galactocentric distances (> ~ 1 kpc). The fact that VCC848 is still (atomic) gas-dominated after the period of most violent collision suggests that gas-rich dwarf galaxy merging can result in BCD-like remnants with extended atomic gas distribution surrounding a blue compact center, in general agreement with previous numerical simulations.

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

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

We present the first systematic study of the stellar populations of ultra-diffuse galaxies (UDGs) in the field, integrating the large area search and characterization of UDGs by the SMUDGes survey with the twelve-band optical photometry of the S-PLUS survey. Based on Bayesian modeling of the optical colors of UDGs, we determine the ages, metallicities and stellar masses of 100 UDGs distributed in an area of $\sim 330$ deg$^2$ in the Stripe 82 region. We find that the stellar masses and metallicities of field UDGs are similar to those observed in clusters and follow the trends previously defined in studies of dwarf and giant galaxies. However, field UDGs have younger luminosity-weighted ages than do UDGs in clusters. We interpret this result to mean that field UDGs have more extended star formation histories, including some that continue to form stars at low levels to the present time. Finally, we examine stellar population scaling relations that show that UDGs are, as a population, similar to other low-surface brightness galaxies.

We present a catalog of 5,290 RR Lyrae stars (RRLs) with metallicities estimated from spectra of the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) surveys. Nearly 70 per cent of them (3,642 objects) also have systemic radial velocities measured. Given the pulsating nature of RRLs, metallicity estimates are based on spectra of individual exposures, by matching them with the synthetic templates. The systemic radial velocities are measured by fitting the observed velocity as a function of phase assuming an empirical pulsating velocity template curve. Various tests show that our analyses yield metallicities with a typical precision of 0.20\,dex and systemic radial velocities with uncertainties ranging from 5 to 21\,km\,s$^{-1}$ (depending on the number of radial velocity measurements available for a given star). Based on the well calibrated near-infrared $PM_{W1}Z$ or $PM_{K_{\rm s}}Z$, and $M_{V}$-[Fe/H] relations, precise distances are derived for these RRLs. Finally, we include Gaia DR2 proper motions in our catalog. The catalog should be very useful for various Galactic studies, especially of the Galactic halo.

In this paper we study the effects of hemispheric imbalance of magnetic helicity density on breaking the equatorial reflection symmetry of the dynamo generated large-scale magnetic field. Our study employs the axisymmetric dynamo model which takes into account the nonlinear effect of magnetic helicity conservation. We find that the evolution of the net magnetic helicity density, in other words, the magnetic helicity imbalance, on the surface follows the evolution of the parity of the large-scale magnetic field. Random fluctuations of the $\alpha$-effect and the helicity fluxes can inverse the causal relationship, i.e., the magnetic helicity imbalance or the imbalance of magnetic helicity fluxes can drive the magnetic parity breaking. We also found that evolution of the net magnetic helicity of the small-scale fields follows the evolution of the net magnetic helicity of the large-scale fields with some time lag. We interpret this as an effect of the difference of the magnetic helicity fluxes out of the Sun from the large and small scales.

We present a catalogue of 3,305 H$\alpha$ emission-line point sources observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33 during September 2011 and January 2016. The catalogue contains 1,487 emission-line stars, 532 emission-line nebulae including 377 likely planetary nebulae (PNe), 83 H~{\textsc{ii}} regions candidates and 20 possible supernovae remnants (SNRs) and 1,286 unknown objects. Among them, 24 PN candidates, 19 H~{\sc ii} region candidates, 10 SNR candidates and 1 symbiotic star candidate are new discoveries. Radial velocities and fluxes estimated from the H$\alpha$ line and those quantities of seven other major emission lines including H$\beta$, [O~{\textsc{iii}}]~$\lambda$4959, [O~{\textsc{iii}}]~$\lambda$5007, [N~{\textsc{ii}}]~$\lambda$6548, [N~{\textsc{ii}}]~$\lambda$6583, [S~{\textsc{ii}}]~$\lambda$6717 and [S~{\textsc{ii}}]~$\lambda$6731 lines of all the catalogued sources yielded from the LAMOST spectra are also presented in our catalogue. Our catalogue is an ideal starting point to study the chemistry properties and kinematics of M31 and M33.

We present a large, homogeneous catalogue of molecular clouds within 4 kpc from the Sun at low Galactic latitudes ($|b|$ $<$ 10\degr) with unprecedented accurate distance determinations. Based on the three-dimensional dust reddening map and estimates of colour excesses and distances of over 32 million stars presented in Chen et al, we have identified 567 dust/molecular clouds with a hierarchical structure identification method and obtained their distance estimates by a dust model fitting algorithm. The typical distance uncertainty is less than 5 per cent. As far as we know, this is the first large catalogue of molecular clouds in the Galactic plane with distances derived in a direct manner. The clouds are seen to lie along the Sagittarius, Local and Perseus Arms. In addition to the known structures, we propose the existence of a possible {\it spur}, with a pitch angle of about 34\degr, connecting the Local and the Sagittarius Arms in the fourth quadrant. We have also derived the physical properties of those molecular clouds. The distribution of cloud properties in different parameter spaces agrees grossly with the previous results. Our cloud sample is an ideal starting point to study the concentration of dust and gas in the solar vicinity and their star formation activities.

We report the discovery of a warm sub-Saturn, TOI-257b (HD 19916b), based on data from NASA's Transiting Exoplanet Survey Satellite (TESS). The transit signal was detected by TESS and confirmed to be of planetary origin based on radial-velocity observations with the Minerva-Australis telescope array. An analysis of the TESS photometry, the Minerva-Australis, FEROS, and HARPS radial velocities, and the asteroseismic data of the stellar oscillations reveals that TOI-257b has a mass of $M_P=0.134^{+0.023}_{-0.022}$$\rm{M_J}$ ($42.6^{+7.3}_{-7.0}$$\rm{M_{\oplus}}$), a radius of $R_P=0.626^{+0.013}_{-0.012}$$\rm{R_J}$ ($7.02^{+0.15}_{-0.13}$$\rm{R_{\oplus}}$), and an orbit with eccentricity $0.242^{+0.040}_{-0.065}$ and period $18.38827\pm0.00072$$\rm{days}$. TOI-257b orbits a bright ($\mathrm{V}=7.570$mag) somewhat evolved late F-type star with $M_*=1.390\pm0.046$$\rm{M_{\odot}}$, $R_*=1.888\pm0.033$$\rm{R_{\odot}}$, $T_{\rm eff}=6075\pm90$$\rm{K}$, and $v\sin{i}=11.3\pm0.5$km/s. Additionally, we statistically validate a second non-transiting sub-Saturn mass planet on a $\sim71$day orbit using the radial velocity data. This system joins the ranks of a small number of exoplanet host stars that have been characterized with asteroseismology. Warm sub-Saturns are rare in the known sample of exoplanets, and thus the discovery of TOI-257b is important in the context of future work studying the formation and migration history of similar planetary systems.

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 the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons.

The weak-field approximation implying linear relationship between Stokes $V/I$ and longitudinal magnetic field, $B_{\Vert}$, often suffers from the saturation observed in strong magnetic field regions such as sunspot umbra. In this work we intend to improve the magnetic field observations carried out by the \textit{Solar Magnetic Field Telescope} (SMFT) at the Huairou Solar Observing Station, China. We propose to use non-linear relationship between Stokes $V/I$ and $B_{\Vert}$ to derive magnetic field. To determine the form of the relationship we perform cross-calibration of SMFT data and magnetograms provided by the \textit{Helioseismic and Magnetic Imager} on board the \textit{Solar Dynamics Observatory}. The algorithm of the magnetic field deriving is described in details. We show that using non-linear relationship between Stokes $V/I$ and $B_{\Vert}$ allows eliminating the magnetic field saturation inside sunspot umbra.

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.

In this paper, we briefly introduce the basic questions in the measurements of solar magnetic fields and the possible error sources due to the approximation of the theory of radiation transfer of spectral lines in the solar atmosphere. We introduce some basic research progress in magnetic field measurement at Huairou Solar Observing Station of National Astronomical Observatories of the Chinese Academy of Sciences, especially concerning the non-potentiality in solar active regions, such as the magnetic shear, current and helicity. We also discuss some basic questions for the measurements of the magnetic fields and corresponding challenges for the future studies.

We present the discoveries of KELT-25b (TIC 65412605, TOI-626.01) and KELT-26b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A-stars. The transit signals were initially detected by the KELT survey, and subsequently confirmed by \textit{TESS} photometry. KELT-25b is on a 4.40-day orbit around the V = 9.66 star CD-24 5016 ($T_{\rm eff} = 8280^{+440}_{-180}$ K, $M_{\star}$ = $2.18^{+0.12}_{-0.11}$ $M_{\odot}$), while KELT-26b is on a 3.34-day orbit around the V = 9.95 star HD 134004 ($T_{\rm eff}$ =$8640^{+500}_{-240}$ K, $M_{\star}$ = $1.93^{+0.14}_{-0.16}$ $M_{\odot}$), which is likely an Am star. We have confirmed the sub-stellar nature of both companions through detailed characterization of each system using ground-based and \textit{TESS} photometry, radial velocity measurements, Doppler Tomography, and high-resolution imaging. For KELT-25, we determine a companion radius of $R_{\rm P}$ = $1.64^{+0.039}_{-0.043}$ $R_{\rm J}$, and a 3-sigma upper limit on the companion's mass of $\sim64~M_{\rm J}$. For KELT-26b, we infer a planetary mass and radius of $M_{\rm P}$ = $1.41^{+0.43}_{-0.51}$ $M_{\rm J}$ and $R_{\rm P}$ = $1.940^{+0.060}_{-0.058}$ $R_{\rm J}$. From Doppler Tomographic observations, we find KELT-26b to reside in a highly misaligned orbit. This conclusion is weakly corroborated by a subtle asymmetry in the transit light curve from the \textit{TESS} data. KELT-25b appears to be in a well-aligned, prograde orbit, and the system is likely a member of a cluster or moving group.

Using the $f\sigma_8(z)$ redshift space distortion (RSD) data, the $\sigma_8^0-\Omega_m^0$ tension is studied utilizing a parameterization of growth rate $f(z) = \Omega_m(z)^\gamma$. Here, $f(z)$ is derived from the expansion history $H(z)$ which is reconstructed from the observational Hubble data applying the Gaussian Process method. It is found that different priors of $H_0$ have a great influence on the evolution curve of $H(z)$ and the constraint of $\sigma_8^0-\Omega_m^0$. When using a larger $H_0$ prior, the low redshifts $H(z)$ deviate significantly from that of the $\Lambda$CDM model, which indicates that a dark energy model different from the cosmological constant can help to relax the $H_0$ tension problem. The tension between our best-fit values of $\sigma_8^0-\Omega_m^0$ and that of the \textit{Planck} 2018 $\Lambda$CDM (PLA) will disappear (less than $1\sigma$) when taking a prior for $H_0$ obtained from PLA. Moreover, the tension exceeds $2\sigma$ level when applying the prior $H_0 = 73.52 \pm 1.62$ resulted from the Hubble Space Telescope photometry.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

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

In this work, we examine the possibility of using the diffuse supernova neutrino background (DSNB) to test the Chaplygin gas (CG) models of the Universe. With a typical supernova rate $R_{\mathrm{SN}}(z)$ and supernova neutrino spectrum $\mathrm{d} N(E_\nu)/\mathrm{d} E_\nu$, the DSNB flux spectrum $n(E_\nu)$ in three categories of CG models, the generalized CG (GCG), modified CG (MCG) and extended CG (ECG) models, are studied. It is found that generally the flux spectra take a form similar to a Fermi-Dirac distribution with a peak centered around 3.80-3.97 MeV. The spectrum shape and peak positions are primarily determined by $R_{\mathrm{SN}}(z)$ and $\mathrm{d} N(E_\nu)/\mathrm{d} E_\nu$ and only slightly affected by the CG models. However, the height of the spectra in each category of the CG models can vary dramatically for different models, with variances of 13.2%, 23.6% and 14.9% for GCG, MCG and ECG categories respectively. The averaged total flux in each category are also different, with the ECG model average 10.0% and 12.7% higher than that of the GCG and MCG models. These suggest that the DSNB flux spectrum height and total flux can be used to constrain the CG model parameters, and if the measured to a sub-10% accuracy, might be used to rule out some models.

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

In this letter we propose that the X-ray and the TeV observations in the vicinity of Geminga can be understood in the framework of anisotropic diffusion of injected electrons/positrons. This interpretation only requires the turbulence in the vicinity of Geminga to be sub-Alfv\'enic with the local mean magnetic field direction approximately aligned with our line of sight towards Geminga, without invoking extreme conditions for the environment, such as an extremely small diffusion coefficient and a weak magnetic field of strength $<1\mu$G as suggested in previous literature.

We continue our empirical study of the emission line flux originating in the cool ($T\sim10^4$ K) gas that populates the halos of galaxies and their environments. Specifically, we present results obtained for a sample of nearly half a million individual galaxies, groups, and clusters of galaxies, intersected by more than two million SDSS lines of sight at projected separations of up to a quarter of the virial radius. Adopting simple power law relationships between the circumgalactic (CGM) cool gas fraction and either the halo or stellar mass, we present expressions for the CGM cool gas fraction as a function of either halo or stellar mass, $f_{\rm cool}(M_h) = (0.25^{+0.07}_{-0.06}) \times (M_h/10^{12}M_\odot)^{(-0.39^{+0.06}_{-0.07})}$ or $f_{\rm cool}(M_{*}) = (0.28^{+0.07}_{-0.06}) \times (M_{\rm *}/10^{10.0}M_\odot)^{(-0.33\pm 0.06)}$. Where we can compare, our results are consistent with previous constraints from absorption line studies, our own previous emission line work, and simulations. The cool gas can be the dominant baryonic CGM component, comprising a fraction as high as $> 90\%$ of halo gaseous baryons, in low mass halos, $M_h\sim$ $10^{10.5} M_\odot$, and a minor fraction, $<$ 5\%, in groups and clusters, $M_h > 10^{14} M_\odot$.

We perform analysis of the three-dimensional kinematics of Milky Way disk stars in mono-age populations. We focus on stars between Galactocentric distances of $R=6$ and 14 \,kpc, selected from the combined LAMOST DR4 red clump giant stars and Gaia DR2 proper motion catalogue. We confirm the 3D asymmetrical motions of recent works and provide time tagging of the Galactic outer disk asymmetrical motions near the anticenter direction out to Galactocentric distances of 14\,kpc. Radial Galactocentric motions reach values up to 10 km s$^{-1}$, depending on the age of the population, and present a north-south asymmetry in the region corresponding to density and velocity substructures that were sensitive to the perturbations in the early 6 \,Gyr. After that time, the disk stars in this asymmetrical structure have become kinematically hotter, and are thus not sensitive to perturbations, and we find the structure is a relatively younger population. With quantitative analysis, we find stars both above and below the plane at $R\gtrsim 9$ kpc that exhibit bending mode motions of which the sensitive duration is around 8 \,Gyr. We speculate that the in-plane asymmetries might not be mainly caused by a fast rotating bar, intrinsically elliptical outer disk, secular expansion of the disk, or streams. Spiral arm dynamics, out-of-equilibrium models, minor mergers or others are important contributors. Vertical motions might be dominated by bending and breathing modes induced by complicated inner or external perturbers. It is likely that many of these mechanisms are coupled together.

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.

During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.

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.

Starburst galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into gamma-rays by colliding with the interstellar medium (ISM). The produced GeV gamma-ray emission is temporally stable and the flux is found to correlate well with indicators of star formation rates, such as the total infrared (IR) luminosity $L_{\rm IR}$ and monochromatic radio continuum luminosity at 1.4 GHz $L_{\rm 1.4\ GHz}$, i.e., following $L_{\gamma}-L_{\rm IR}$ and $L_{\gamma}-L_{\rm 1.4\ GHz}$ relations. Recently, gamma-ray excesses are reported to be spatially coincident with two starburst galaxies NGC 3424 and UGC 11041 in the fourth \textsl{Fermi} Large Area Telescope (LAT) source catalog (4FGL). Different from other starburst galaxies detected by \textsl{Fermi}-LAT, we find that the gamma-ray emission associated with NGC 3424 and UGC 11041 show significant flux variability. With relatively weak infrared and radio emission, NGC 3424 and UGC 11041 appear as outliers of the $L_{\gamma}-L_{\rm IR}$ and $L_{\gamma}-L_{\rm 1.4\ GHz}$ relations of starburst galaxies. These results suggest that NGC 3424 and UGC 11041 may harbor obscured active galactic nuclei (AGNs) and the AGN activities provide the dominant contribution to the gamma-ray emission as compared to that provided by the starburst activities.

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

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

The occurrence of planetary nebulae (PNe) in globular clusters (GCs) provides an excellent chance to study low-mass stellar evolution in a special (low-metallicity, high stellar density) environment. We report a systematic spectroscopic survey for the [O{\sc iii}] 5007 emission line of PNe in 1469 Virgo GCs and 121 Virgo ultra-compact dwarfs (UCDs), mainly hosted in the giant elliptical galaxies M87, M49, M86, and M84. We detected zero PNe in our UCD sample and discovered one PN ($M_{5007} = -4.1$ mag) associated with an M87 GC. We used the [O{\sc iii}] detection limit for each GC to estimate the luminosity-specific frequency of PNe, $\alpha$, and measured $\alpha$ in the Virgo cluster GCs to be $\alpha \sim 3.9_{-0.7}^{+5.2}\times 10^{-8}\mathrm{PN}/L_\odot$. $\alpha$ in Virgo GCs is among the lowest values reported in any environment, due in part to the large sample size, and is 5--6 times lower than that for the Galactic GCs. We suggest that $\alpha$ decreases towards brighter and more massive clusters, sharing a similar trend as the binary fraction, and the discrepancy between the Virgo and Galactic GCs can be explained by the observational bias in extragalactic surveys toward brighter GCs. This low but non-zero efficiency in forming PNe may highlight the important role played by binary interactions in forming PNe in GCs. We argue that a future survey of less massive Virgo GCs will be able to determine whether PN production in Virgo GCs is governed by internal process (mass, density, binary fraction), or is largely regulated by external environment.

The Minerva-Australis telescope array is a facility dedicated to the follow-up, confirmation, characterisation, and mass measurement of bright transiting planets discovered by the Transiting Exoplanet Survey Satellite (TESS) -- a category in which it is almost unique in the southern hemisphere. It is located at the University of Southern Queensland's Mount Kent Observatory near Toowoomba, Australia. Its flexible design enables multiple 0.7m robotic telescopes to be used both in combination, and independently, for high-resolution spectroscopy and precision photometry of TESS transit planet candidates. Minerva-Australis also enables complementary studies of exoplanet spin-orbit alignments via Doppler observations of the Rossiter-McLaughlin effect, radial velocity searches for non-transiting planets, planet searches using transit timing variations, and ephemeris refinement for TESS planets. In this first paper, we describe the design, photometric instrumentation, software, and science goals of Minerva-Australis, and note key differences from its Northern hemisphere counterpart -- the Minerva array. We use recent transit observations of four planets--WASP-2b, WASP-44b, WASP-45b, and HD 189733b to demonstrate the photometric capabilities of Minerva-Australis.

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.

Imposing an external magnetic field in short-pulse intense laser-plasma interaction is of broad scientific interest in related plasma research areas. We propose a simple method using a virtual current layer by introducing an extra current density term to simulate the external magnetic field, and demonstrate it with three-dimensional particle-in-cell simulations. The field distribution and its evolution in sub-picosecond time scale are obtained. The magnetization process takes a much longer time than that of laser-plasma interaction due to plasma diamagnetism arising from collective response. The long-time evolution of magnetic diffusion and diamagnetic current can be predicted based on a simplified analytic model in combination with simulations.

We study a simple class of dark matter models with N_f copies of electroweak fermionic multiplets, stabilized by O(N_F) global symmetry. Unlike conventional minimal dark matter which usually suffers from Landau poles, in these models the gauge coupling g_2 has a non-trivial ultraviolet fixed point, and thus is asymptotically safe as long as N_F is large enough. These fermionic n-plet models have only two free parameters: N_F and a common mass M_DM, which relate to dark matter relic abundance. We find that the mass of triplet fermionic dark matter with N_F being dozens of flavors can be several hundred GeV, which is testable on LHC. A benefit of large N_F is that DM pair annihilation rate in dwarf galaxies is effectively suppressed by 1/N_F, and thus they can evade the constraint from gamma-ray continuous spectrum observation. For the case of triplets, we find that the models in the range 3 <= N_F <= 20 are consistent with all current experiments. However, for N_F quintuplets, even with large N_F they are still disfavored by the gamma-ray continuous spectrum.

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

The hot massive luminous blue variables (LBVs) represent an important evolutionary phase of massive stars. Here, we report the discovery of a new LBV -- LAMOST J0037+4016 in the distant outskirt of the Andromeda galaxy. It is located in the south-western corner (a possible faint spiral arm) of M31 with an unexpectedly large projection distance of $\sim$ 22 kpc from the center. The optical light curve shows a 1.2 mag variation in $V$ band and its outburst and quiescence phases both last over several years. The observed spectra indicate an A-type supergiant at epoch close to the outburst phase and a hot B-type supergiant with weak [Fe II] emission lines at epoch of much dimmer brightness. The near-infrared color-color diagram further shows it follows the distribution of Galactic and M31 LBVs rather than B[e] supergiants. All the existing data strongly show that LAMOST J0037+4016 is an LBV. By spectral energy distribution fitting, we find it has a luminosity ($4.42 \pm 1.64$)$\times 10^5$ $L_{\odot}$ and an initial mass $\sim 30$ $M_{\odot}$, indicating its nature of less luminosity class of LBV.

We investigate the three-dimensional asymmetrical kinematics and present time stamps of the Milky Way disk between Galactocentric distances of $R=12$ and 15 \,kpc, using red clump stars selected from the LAMOST Galactic survey, also with proper motion measurements provided by the Gaia DR2. We discover velocity substructure above the Galactic plane corresponding to a density dip found recently ("south-middle opposite" density structure[R $\sim$ 12-15 \,kpc, Z $\sim$ 1.5 \,kpc] discovered in \citet{Wang2018b, Wang2018c}) in the radial and azimuthal velocity. For the vertical velocity, we detect clear vertical bulk motions or bending mode motions, which has no clear north-south asymmetry corresponding to the in-plane asymmetrical features. In the subsample of stars with different ages, we find that there is little temporal evolution of the in-plane asymmetry from 0$-$14 \,Gyr, which means the structure is sensitive to the perturbations in almost cosmic time possibly. We propose that the possible scenario of this asymmetric velocity structure is caused by the mechanisms generated in-plane, rather than vertical perturbations.

Recently, very high-energy photons above 100 GeV were reported to be detected from GRB~190114C and GRB~180720B at, respectively, 100-1000 s and 10 hours after the burst. We model the available broad-band data of both GRBs with the synchrotron plus synchrotron self-Compton (SSC) emission of the afterglow shocks. We find that the sub-TeV emission of GRB~180720B can be interpreted as the SSC emission from afterglow shocks expanding in a constant density circum-burst medium. The SSC emission of GRB~190114C dominates over the synchrotron component from GeV energies at $\sim100$ s, which can explain the possible hard spectrum of the GeV emission at this time. The extrapolated flux of this SSC component to sub-TeV energies can explain the high-significance detection of GRB~190114C by the MAGIC telescope. The parameter values (such as the circum-burst density and shock microphysical parameters) in the modeling are not unusual for both GRBs, implying that the detection of sub-TeV photons from these two bursts should be attributed to their large burst energies and low redshifts.

The timing argument provides a lower limit on the mass of the Milky Way. We find, using a sample of 32 stars at $R > 60$ kpc drawn from the H3 Spectroscopic Survey and mock catalogs created from published numerical simulations, that M$_{200} > 0.91\times 10^{12}$ M$_\odot$ with 90% confidence. We recommend using this limit to refine the allowed prior mass range in more complex and sophisticated statistical treatments of Milky Way dynamics. The use of such a prior would have significantly reduced many previously published uncertainty ranges. Our analysis suggests that the most likely value of M$_{200}$ is $\sim 1.4 \times 10^{12}$ M$_\odot$, but establishing this as the Milky Way mass requires a larger sample of outer halo stars and a more complete analysis of the inner halo stars in H3. The imminent growth in the sample of outer halo stars due to ongoing and planned surveys will make this possible.

We extract 18 candidate short gamma-ray bursts (SGRBs) with precursors from 660 SGRBs observed by {\em Fermi} and {\em Swift} satellites, and carry out a comprehensive analysis on their temporal and spectral features. We obtain the following results: (1) For a large fraction of candidates, the main burst durations are longer than their precursor durations, comparable to their quiescent times from the end of precursors to the beginning of their main bursts. (2) The average flux of precursors tends to increase as their main bursts brighten. (3) As seen from the distributions of hardness ratio and spectral fitting, the precursors are slightly spectrally softer with respect to the main bursts. Moreover, a significant portion of precursors and all main bursts favor a non-thermal spectrum. (4) The precursors might be a probe of the progenitor properties of SGRBs such as the magnetic field strength and the crustal equation of state if they arise from some processes before mergers of binary compact objects rather than post-merger processes.

We present the determination of stellar parameters and individual elemental abundances for 6 million stars from $\sim$8 million low-resolution ($R\sim1800$) spectra from LAMOST DR5. This is based on a modeling approach that we dub $The$ $Data$--$Driven$ $Payne$ ($DD$--$Payne$), which inherits essential ingredients from both {\it The Payne} \citep{Ting2019} and $The$ $Cannon$ \citep{Ness2015}. It is a data-driven model that incorporates constraints from theoretical spectral models to ensure the derived abundance estimates are physically sensible. Stars in LAMOST DR5 that are in common with either GALAH DR2 or APOGEE DR14 are used to train a model that delivers stellar parameters ($T_{\rm eff}$, $\log g$, $V_{\rm mic}$) and abundances for 16 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, Co, Ni, Cu, and Ba) when applied to LAMOST spectra. Cross-validation and repeat observations suggest that, for ${\rm S/N}_{\rm pix}\ge 50$, the typical internal abundance precision is 0.03--0.1\,dex for the majority of these elements, with 0.2--0.3\,dex for Cu and Ba, and the internal precision of $T_{\rm eff}$ and $\log g$ is better than 30\,K and 0.07\,dex, respectively. Abundance systematics at the $\sim$0.1\,dex level are present in these estimates, but are inherited from the high-resolution surveys' training labels. For some elements, GALAH provides more robust training labels, for others, APOGEE. We provide flags to guide the quality of the label determination and to identify binary/multiple stars in LAMOST DR5. The abundance catalogs are publicly accessible via \href{url}{http://dr5.lamost.org/doc/vac}.

HISPEC (High-resolution Infrared Spectrograph for Exoplanet Characterization) is a proposed diffraction-limited spectrograph for the W.M. Keck Observatory, and a pathfinder for the MODHIS facility project (Multi-Object Diffraction-limited High-resolution Infrared Spectrograph) on the Thirty Meter Telescope. HISPEC/MODHIS builds on diffraction-limited spectrograph designs which rely on adaptively corrected single-mode fiber feeds. Seeing-limited high-resolution spectrographs, by virtue of the conservation of beam etendue, grow in volume following a D^3 power law (D is the telescope diameter), and are subject to daunting challenges associated with their large size. Diffraction-limited spectrographs fed by single mode fibers are decoupled from the telescope input, and are orders of magnitude more compact and have intrinsically stable line spread functions. Their efficiency is directly proportional to the performance of the adaptive optics (AO) system. AO technologies have matured rapidly over the past two decades and are baselined for future extremely large telescopes. HISPEC/MODHIS will take R>100,000 spectra of a few objects in a 10" field-of-view sampled at the diffraction limit (~10-50 mas), simultaneously from 0.95 to 2.4 microns (y-K). The scientific scope ranges from exoplanet infrared precision radial velocities, spectroscopy of transiting, close-in, and directly imaged exoplanets (atmospheric composition and dynamics, RM effect, spin measurements, Doppler imaging), brown dwarf characterization, stellar physics/chemistry, proto-planetary disk kinematics/composition, Solar system, extragalactic science, and cosmology. HISPEC/MODHIS features a compact, cost-effective design optimized to fully exploit the existing Keck-AO and future TMT-NFIRAOS infrastructures and boost the scientific reach of Keck Observatory and TMT soon after first light.

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.

Using the photometric data from the Next Generation Fornax Survey, we find a significant radial alignment signal among the Fornax dwarf galaxies. For the first time, we report that the radial alignment signal of nucleated dwarfs is stronger than that of non-nucleated ones at 2.4$\sigma$ confidence level, and the dwarfs located in the outer region ($R>R_{\rm{vir}}/3$; $R_{\rm{vir}}$ is the Fornax virial radius) show slightly stronger radial alignment signal than those in the inner region ($R<R_{\rm{vir}}/3$) at $1.5\sigma$ level. We also find that the significance of radial alignment signal is independent of the luminosities or sizes of the dwarfs.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

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.

Transiting Exoplanet Survey Satellite (TESS) observations have revealed a compact multi-planet system around the sixth-magnitude star HR 858 (TIC 178155732, TOI 396), located 32 parsecs away. Three planets, each about twice the size of Earth, transit this slightly-evolved, late F-type star, which is also a member of a visual binary. Two of the planets may be in mean motion resonance. We analyze the TESS observations, using novel methods to model and remove instrumental systematic errors, and combine these data with follow-up observations taken from a suite of ground-based telescopes to characterize the planetary system. The HR 858 planets are enticing targets for precise radial velocity observations, secondary eclipse spectroscopy, and measurements of the Rossiter-McLaughlin effect.

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.

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work.

Precision measurements of charged cosmic rays have recently been carried out by space-born (e.g. AMS-02), or ground experiments (e.g. HESS). These measured data are important for the studies of astro-physical phenomena, including supernova remnants, cosmic ray propagation, solar physics and dark matter. Those scenarios usually contain a number of free parameters that need to be adjusted by observed data. Some techniques, such as Markov Chain Monte Carlo and MultiNest, are developed in order to solve the above problem. However, it is usually required a computing farm to apply those tools. In this paper, a genetic algorithm for finding the optimum parameters for cosmic ray injection and propagation is presented. We find that this algorithm gives us the same best fit results as the Markov Chain Monte Carlo but consuming less computing power by nearly 2 orders of magnitudes.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We measure differences in the emission line flux from the circumgalactic medium, CGM, of galaxies in different environments. Such differences could be a critical clue in explaining a range of galaxy properties that depend on environment. Using large samples of stacked archival spectra from the Sloan Digital Sky Survey, we find that the H$\alpha$ + [N {\small II}] emission line flux from the CGM within 50 kpc of $\sim$ L$^*$ galaxies is lower both for galaxies that lie within a projected distance of $\sim$ 500 kpc from a massive ($M_*>10^{11} M_\odot$) galaxy and for galaxies in richer/denser environments. The environmental differences are statistically significant even after we control for galaxy mass and morphology. We interpret these observations as a direct signature of environmentally-caused strangulation. We present a simple, heuristic model for the effect of a massive parent galaxy. In this model, the CGM cool gas fraction within 50 kpc is significantly decreased for galaxies that lie within 700 kpc of a massive galaxy, with about 80\% of the cool gas removed even when the galaxy is at a distance of 500 kpc from its massive parent. However, we discuss alternative physical causes for the observed behavior and discuss ways forward in addressing open questions.

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$\alpha$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$\alpha$ emission of the LABs. The Ly$\alpha$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$\alpha$ luminosity relation of the Ly$\alpha$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

Pulsar Wind Nebula (PWN) DA 495 (G65.7+1.2) was detected in TeV gamma-rays by the High Altitude Water Cherenkov Observatory (HAWC) in 2017 (2HWC J1953+294). Follow-up observations by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) confirmed the association between 2HWC J1953+294 and DA 495 and found the TeV emission to be spatially coincident with the radio emission first reported in 1968. The detection of TeV gamma-rays from DA 495, along with past X-ray detection up to 10 keV, prompted high energy X-ray observations as part of the NuSTAR Galactic Legacy Survey. We present the results of these NuSTAR observations, combined with archival Chandra and XMM-Newton observations, and confirm the previous X-ray photon index of $\Gamma_{2-20 \rm\ keV} = 2.0 \pm 0.1$. We find no spectral cutoff up to 20 keV. With the spectral information for DA 495 extended to TeV gamma-rays, we were able to perform analytical modeling to test leptonic and hadronic emission scenarios. The leptonic models can explain the broadband emission, but also imply a diffuse X-ray nebula of similar extent to the radio and TeV nebulae, which cannot be confirmed by our observations. The hadronic models can simultaneously explain the spectrum and the spatial extent in all wavelengths; however, we need a very high magnetic field strength pervading the radio and TeV nebulae and a surprisingly high particle kinetic energy. These requirements deepen the mystery of the physical nature of DA 495. Future observations in radio to infrared bands and spatially resolved $\gamma$-rays can further constrain the physical conditions and radiation mechanisms in DA 495.

Accurate determinations of atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log $g$ and metallicity [Fe/H]) and distances for large complete samples are of vital importance for various Galactic studies. We have developed a photometric method to select red giant stars and estimate their atmospheric parameters from the photometric colors provided by the SkyMapper Southern Survey (SMSS) data release (DR) 1.1, using stars in common with the LAMOST Galactic spectroscopic surveys as a training set. Distances are estimated with two different approaches: one based on the Gaia DR2 parallaxes for nearby ($d \leq 4.5$ kpc) bright stars and another based on the absolute magnitudes predicted by intrinsic color $(g-i)_0$ and photometric metallicity [Fe/H] for distant ($d > 4.5$ kpc) faint stars. Various tests show that our method is capable of delivering atmospheric parameters with a precision of $\sim$80 K for $T_{\rm eff}$, $\sim$0.18 dex for [Fe/H] and $\sim$0.35 dex for log $g$, but with a significant systematic error at log $g \sim$ 2.3. For distances delivered from $(g-i)_0$ and photometric [Fe/H], our test with the member stars of globular clusters show a median uncertainty of 16 per cent with a negligible zero-point offset. Using this method, atmospheric parameters and distances of nearly one million red giant stars are derived from SMSS DR1.1. Proper motion measurements from Gaia DR2 are available for almost all of the red giant stars, and radial velocity measurements from several large spectroscopic surveys are available for 44 per cent of these. This sample will be accessible online at https://yanghuang0.wixsite.com/yangh/research .

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 equipartition of magnetic and thermal energy in the interstellar medium (ISM) indicates the magneto-hydrodynamic nature of the interstellar turbulence, which can be decomposed into three wave modes: Alfv\'en, fast and slow magnetosonic modes\cite{Herlofson50,LG01,CL03}. Even for studies performed in the case of subsonic turbulence, the magnetosonic modes have been revealed to be more effective than the Alfv\'enic modes for processes such as cosmic ray (CR) transport and acceleration \cite{YL04,Lynn14}. The multiphase nature of ISM and diversity of driving mechanisms give rise to spatial variation of turbulence properties. Nevertheless, the employed model of magneto-hydrodynamic turbulence is often oversimplified being assumed to be only Alfv\'enic or even hydrodynamic due to a lack of observational evidence. Here we report the employment of our novel method, the signature from polarization analysis (SPA), on unveiling the plasma modes in interstellar turbulence. Its application leads to the first discovery of magnetosonic modes in the Cygnus X region. It is found that the magnetosonic modes overlap to a high degree with Fermi cocoon, consistent with theoretical expectations. Moreover, through comparison with the spectrum at other wavelengths, the plasma modes of turbulence are unveiled in active star formation zones. The SPA casts light on the plasma modes composition of the Galactic turbulence, and marks the onset of a new era in the study of interstellar turbulence and accordingly our understandings of relevant processes including cosmic ray transport and star formation.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

The circumgalactic medium (CGM), which harbors > 50% of all the baryons in a galaxy, is both the reservoir of gas for subsequent star formation and the depository of chemically processed gas, energy, and angular momentum from feedback. As such, the CGM obviously plays a critical role in galaxy evolution. We discuss the opportunity to image this component using recombination line emission, beginning with the early results coming from recent statistical detection of this emission to the final goal of realizing spectral-line images of the CGM in individual nearby galaxies. Such work will happen in the next decade and provide new insights on the galactic baryon cycle.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present a new, model-independent method to reconstruct the temporal evolution of the speed of light $c(z)$ using astronomical observations. After validating our pipeline using mock datasets, we apply our method to the latest BAO and supernovae observations, and reconstruct $c(z)$ in the redshift range of $z\in[0,1.5]$. We find no evidence of a varying speed of light, although we see some interesting features of $\Delta c(z)$, the fractional difference between $c(z)$ and $c_0$ (the speed of light in SI), \eg, $\Delta c(z)<0$ and $\Delta c(z)>0$ at $0.2\lesssim z\lesssim0.5$ and $0.8\lesssim z\lesssim1.3$, respectively, although the significance of these features is currently far below statistical importance.

With the photometric data from the SDSS survey, the spectroscopic data from the SDSS/SEGUE and the LAMOST surveys, and the astrometric data from the Gaia DR2, we have identified 71 highly-probable member stars of the GD-1 cold stellar stream spread along almost its entire length (i.e. from 126 to 203 degree in Right Ascension). With the accurate spectroscopic (i.e. metallicity and line-of-sight velocity) and astrometric (i.e. proper motions) information, the position-velocity diagrams, i.e. $\phi_{1}$-$\mu_{\alpha}$, $\phi_{1}$-$\mu_{\delta}$ and $\phi_{1}$-${v_{\rm gsr}}$, of the GD-1 stream are well mapped. The stream has an average metallicity [Fe/H] $= -1.95 \pm 0.23$. The rich information of member stars of the stream now available allow one not only to model its origin, but also to place strong constraints on the mass distribution and the gravitational potential of the Milky Way.

We have investigated the distributions of stellar azimuthal and radial velocity components $V_{\Phi}$ and $V_{R}$ in the vertical position-velocity plane $Z$-$V_{Z}$ across the Galactic disc of $6.34 \lesssim R \lesssim 12.34$\,kpc and $|\Phi| \lesssim 7.5^{\circ}$ using a Gaia and Gaia-LAMOST sample of stars. As found in previous works, the distributions exhibit significant spiral patterns. The $V_{R}$ distributions also show clear quadrupole patterns, which are the consequence of the well-known tilt of the velocity ellipsoid. The observed spiral and quadrupole patterns in the phase space plane vary strongly with radial and azimuthal positions. The phase spirals of $V_{\Phi}$ become more and more relaxed as $R$ increases. The spiral patterns of $V_{\Phi}$ and $V_{R}$ and the quadrupole patterns of $V_{R}$ are strongest at $-2^{\circ} < \Phi < 2^{\circ}$ but negligible at $4^{\circ} < \Phi < 6^{\circ}$ and $-6^{\circ} < \Phi < -4^{\circ}$. Our results suggest an external origin of the phase spirals. In this scenario, the intruder, most likely the previously well-known Sagittarius dwarf galaxy, passed through the Galactic plane in the direction towards either Galactic center or anti-center. The azimuthal variations of the phase spirals also help us constrain the passage duration of the intruder. A detailed model is required to reproduce the observed radial and azimuthal variations of the phase spirals of $V_{\Phi}$ and $V_{R}$.

Sub-Gauss magnetic fields are crucial for many physical processes from Galactic $kpc$ scale to stellar $Au$ scale, but are nontrivial to detect. All present magnetic diagnoses trace only one component of a magnetic field at best. Here we report the first observational results that unveil 3D topology as well as the strength of a sub-milliGauss magnetic field beyond our solar system. We found that two weak neutral iron absorption lines from the ground state in the atmosphere of $89~Her$ produced counterintuitive high-amplitude polarizations and consistent polarization angles, exactly in line with the theoretical prediction from ground state alignment (GSA). Our analysis reveals the first sub-AU scale magnetic field on $89~Her$, that is $1.3kpc$ away, has a 3D orientation aligned to the stellar outflow axis and a strength of $70\mu{G} \lesssim B \lesssim150\mu{G}$, thus substantially improving the accuracy by five orders of magnitude compared to the previous $10G$ upper limit set by non-detection of Zeeman effect. This long-awaited discovery is anticipated to usher in a new era of precise measurement of sub-gauss magnetic field in the Universe.

We present optical luminosity functions (LFs) of galaxies in the $^{0.1}g$, $^{0.1}r$, and $^{0.1}i$ bands, calculated using data in $\sim 40$ $deg^{2}$ sky area of LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in Southern Galactic Cap. Redshifts for galaxies brighter $r = 18.1$ were obtained mainly with LAMOST. In each band, LFs derived using both parametric and non-parametric maximum likelihood methods agree well with each other. In the $^{0.1}r$ band, our fitting parameters of the Schechter function are $\phi_{*}=(1.65\pm0.36)\times10^{-2}h^{3}Mpc^{-3}$, $M_{*}=-20.69\pm0.06$ mag, and $\alpha=-1.12\pm0.08$, in agreements with previous studies. Separate LFs are also derived for emission line galaxies and absorption line galaxies, respectively. The LFs of absorption line galaxies show a dip at $^{0.1}r \sim 18.5$ and can be well fitted by a double-Gaussian function, suggesting a bi-modality in passive galaxies.

The Commensal Radio Astronomy Five-hundred-meter Aperture Spherical radio Telescope (FAST) Survey (CRAFTS) utilizes the novel drift-scan commensal survey mode of FAST and can generate billions of pulsar candidate signals. The human experts are not likely to thoroughly examine these signals, and various machine sorting methods are used to aid the classification of the FAST candidates. In this study, we propose a new ensemble classification system for pulsar candidates. This system denotes the further development of the pulsar image-based classification system (PICS), which was used in the Arecibo Telescope pulsar survey, and has been retrained and customized for the FAST drift-scan survey. In this study, we designed a residual network model comprising 15 layers to replace the convolutional neural networks (CNNs) in PICS. The results of this study demonstrate that the new model can sort >96% of real pulsars to belong the top 1% of all candidates and classify >1.6 million candidates per day using a dual--GPU and 24--core computer. This increased speed and efficiency can help to facilitate real-time or quasi-real-time processing of the pulsar-search data stream obtained from CRAFTS. In addition, we have published the labeled FAST data used in this study online, which can aid in the development of new deep learning techniques for performing pulsar searches.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

Understanding the sources of lunar water is crucial for studying the history of lunar evolution, and also the solar wind interaction with the Moon and other airless bodies. Recent observations revealed lunar hydration is very likely a surficial dynamic process driven by solar wind. Solar wind is shielded over a period of 3-5 days as the Moon passes through the Earth's magnetosphere, during which a significant loss of hydration is expected from previous works.Here we study lunar hydration inside the magnetosphere using orbital spectral data, which unexpectedly found that the polar surficial OH/H2O abundance remains at the same level when in the solar wind and in the magnetosphere. We suggest that particles from the magnetosphere (Earth wind, naturally different from solar wind) contribute to lunar hydration. From lunar orbital plasma observations, we find the existence of optimal energy ranges, other than 1 keV as previously thought, for surface hydration formation. These optimal energy ranges deduced from space observations may provide strong implications for laboratory experiments simulating lunar hydration processes.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and the neutrino flux of hadronic models, and predict their X-ray to $\gamma$-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic (RMHD) simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multi-wavelength polarization variability of the proton synchrotron model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and $\gamma$-ray bands, and to infer the magnetic field strength in the emission region. Additionally, the temporal multi-wavelength polarization signatures in the proton synchrotron model shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multi-wavelength polarimetry facilities such as {\it IXPE} and {\it AMEGO} together with neutrino telescopes such as {\it IceCube} can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.

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.

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

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

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

It is known that de Rham-Gabadadze-Tolley (dRGT) massive gravity does not permit a homogeneous and isotropic universe with flat or spherical spatial metrics. We demonstrate that a singular reference metric solves this problem in an economic and straightforward way. In the dRGT massive gravity with a singular reference metric, there are sound homogeneous and isotropic cosmological solutions. We investigate cosmologies with the static and dynamical singular reference metrics, respectively. The term like dark energy appears naturally and the universe accelerates itself in some late time evolution. The term simulating dark matter also naturally emerges. We make a preliminary constraint on the parameters in the dRGT massive gravity in frame of the present cosmological model by using the data of supernovae, cosmic microwave back ground radiations, and baryonic acoustic oscillations.

Protoplanetary disks (PPDs) are characterized by different kinds of gas dynamics and chemistry, which are coupled via ionization, heating and cooling processes, as well as advective and turbulent transport. However, directly coupling gas dynamics with time-dependent chemistry is prohibitively computationally expensive when using comprehensive chemical reaction networks. In this paper, we evaluate the utility of a species-based network reduction method in different disk environments to produce small chemical networks that reproduce the abundances of major species found in large gas-phase chemistry networks. We find that the method works very well in disk midplane and surfaces regions, where approximately 20-30 gas phase species, connected by $\sim$50-60 gas phase reactions, are sufficient to reproduce the targeted ionization fraction and chemical abundances. Most species of the reduced networks, including major carriers of oxygen, carbon and nitrogen, also have similar abundances in the reduced and complete network models. Our results may serve as an initial effort for future hydrodynamic/magnetohydrodynamic simulations of PPDs incorporating time-dependent chemistry in appropriate regions. Accurately modeling the abundances of major species at intermediate disk heights, however, will require much more extended network incorporating gas-grain chemistry and are left for future studies.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

We present accurate fits for the remnant properties of generically precessing binary black holes, trained on large banks of numerical-relativity simulations. We use Gaussian process regression to interpolate the remnant mass, spin, and recoil velocity in the 7-dimensional parameter space of precessing black-hole binaries with mass ratios $q\leq2$, and spin magnitudes $\chi_1,\chi_2\leq0.8$. For precessing systems, our errors in estimating the remnant mass, spin magnitude, and kick magnitude are lower than those of existing fitting formulae by at least an order of magnitude (improvement is also reported in the extrapolated region at high mass ratios and spins). In addition, we also model the remnant spin and kick directions. Being trained directly on precessing simulations, our fits are free from ambiguities regarding the initial frequency at which precessing quantities are defined. We also construct a model for remnant properties of aligned-spin systems with mass ratios $q\leq8$, and spin magnitudes $\chi_1,\chi_2\leq0.8$. As a byproduct, we also provide error estimates for all fitted quantities, which can be consistently incorporated into current and future gravitational-wave parameter-estimation analyses. Our model(s) are made publicly available through a fast and easy-to-use Python module called surfinBH.

We take advantage of the red clump stars to build the relation of the optical extinction (AV) and distance in each direction of supernova remnants (SNRs) with known extinction in the fourth Galactic quadrant. The distances of 9 SNRs are well determined by this method. Their uncertainties range from 10% to 30%, which is significantly improved for 8 SNRs, G279.0+1.1, G284.3-1.8, G296.1-0.5, G299.2-2.9, G308.4-1.4, G309.2-0.6, G309.8- 2.6, G332.4-0.4. In addition, SNR G284.3-1.8 with the new distance of 5.5 kpc is not likely associated with the PSR J1016-5857 at 3 kpc.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

The absence of large-angle correlations in the map of cosmic microwave background temperature fluctuations is among the well-established anomalies identified in full-sky and cut-sky maps over the past three decades. Suppressed large-angle correlations are rare statistical flukes in standard inflationary cosmological models. One natural explanation could be that the underlying primordial density perturbations lack correlations on large distance scales. To test this idea, we replace Fourier modes by a wavelet basis with compact spatial support. While the angular correlation function of perturbations can readily be suppressed, the observed monopole and dipole-subtracted correlation function is not generally suppressed. This suggests that suppression of large-angle temperature correlations requires a mechanism that has both real-space and harmonic-space effects.

We present new three-dimensional (3D) interstellar dust reddening maps of the Galactic plane in three colours, E(G-Ks), E(Bp-Rp) and E(H-Ks). The maps have a spatial angular resolution of 6 arcmin and covers over 7000 deg$^2$ of the Galactic plane for Galactic longitude 0 deg $<$ $l$ $<$ 360 deg and latitude $|b|$ $<$ $10$ deg. The maps are constructed from robust parallax estimates from the Gaia Data Release 2 (Gaia DR2) combined with the high-quality optical photometry from the Gaia DR2 and the infrared photometry from the 2MASS and WISE surveys. We estimate the colour excesses, E(G-Ks), E(Bp-Rp) and E(H-Ks), of over 56 million stars with the machine learning algorithm Random Forest regression, using a training data set constructed from the large-scale spectroscopic surveys LAMOST, SEGUE and APOGEE. The results reveal the large-scale dust distribution in the Galactic disk, showing a number of features consistent with the earlier studies. The Galactic dust disk is clearly warped and show complex structures possibly spatially associated with the Sagittarius, Local and Perseus arms. We also provide the empirical extinction coefficients for the Gaia photometry that can be used to convert the colour excesses presented here to the line-of-sight extinction values in the Gaia photometric bands.

Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses.

We present a homogeneous catalog of 275 large (effective radius $\gtrsim$ 5.3 arcsec) ultra-diffuse galaxy (UDG) candidates lying within an $\approx$ 290 square degree region surrounding the Coma cluster. The catalog results from our automated postprocessing of data from the Legacy Surveys, a three-band imaging survey covering 14,000 square degrees of the extragalactic sky. We describe a pipeline that identifies UDGs and provides their basic parameters. The survey is as complete in these large UDGs as previously published UDG surveys of the central region of the Coma cluster. We conclude that the majority of our detections are at roughly the distance of the Coma cluster, implying effective radii $\ge 2.5$ kpc, and that our sample contains a significant number of analogs of DF 44, where the effective radius exceeds 4 kpc, both within the cluster and in the surrounding field. The $g-z$ color of our UDGs spans a large range, suggesting that even large UDGs may reflect a range of formation histories. A majority of the UDGs are consistent with being lower stellar mass analogs of red sequence galaxies, but we find both red and blue UDG candidates in the vicinity of the Coma cluster and a relative overabundance of blue UDG candidates in the lower density environments and the field. Our eventual processing of the full Legacy Surveys data will produce the largest, most homogeneous sample of large UDGs.

We carry out a project to independently measure the distances of supernova remnants (SNRs) in the first quadrant of the Galaxy. In this project, red clump (RC) stars are used as standard candles and extinction probes to build the optical extinction (A$_V$) - distance(D) relation in each direction of extinction-known SNRs. 15 SNRs' distances are well determined. Among them, the distances of G65.8-0.5, G66.0-0.0 and G67.6+0.9 are given for the first time. We also obtain 32 upper/lower limits of distances, and the distances to G5.7-0.1, G15.1-1.6, G28.8+1.5 and G78.2+2.1 are constrained. Most of the distances measured by the RC method are consistent with previous results. The RC method provides an independent access to the distances of SNRs.

We report the photometric properties of 16 dwarf galaxies, 15 of which are newly identified, in the Western halo of the nearby giant elliptical galaxy NGC5128. All candidates are found at projected distances $\sim\!100\!-\!225$kpc from their giant host, with luminosities $-10.82\!\leq\!M_V/{\rm mag}\!\leq\!-7.42$ and effective radii $4''\!\leq\!r_{\rm eff}\!\leq\!17''$ (or $75\!\leq\!r_{\rm eff}/{\rm pc}\!\leq\!300$ at the distance of NGC5128). We compare to other low-mass dwarf galaxies in the local universe and find them to populate the faint/compact extension of the size-luminosity relation that was previously not well-sampled by dwarf galaxies in the Centaurus A system, with optical colors similar to compact stellar systems like globular clusters and ultra-compact dwarf galaxies despite having much more diffuse morphologies. From optical $u'g'r'i'z'$ photometry, stellar masses are estimated to be $5.17\!\leq\!\log{\cal M}_\star/M_\odot\!\leq\!6.48$, with colors that show them to fall redward of the dwarf galaxy mass-metallicity relation. These colors suggest star formation histories that require some mechanism that would give rise to extra metal enrichment such as primordial formation within the halos of their giant galaxy hosts, non-primordial star formation from previously enriched gas, or extended periods of star formation leading to self-enrichment. We also report the existence of at least two sub-groups of dwarf candidates, each subtending $10-20'$ on the sky, corresponding to projected physical separations of $10\!-\!20$kpc. True physical associations of these groups, combined with their potentially extended star formation histories, would imply that they may represent dwarf galaxy groups in the early stage of interaction upon infall into a giant elliptical galaxy halo in the very nearby universe.

We have investigated the metallicity distributions of mono-age stellar populations across the disc of $6 \lesssim R \lesssim 12$\,kpc and $|Z| \lesssim 2$\,kpc using samples selected from the main-sequence turn-off and sub-giant (MSTO-SG) stars targeted by the LAMOST Galactic Spectroscopic surveys. Both the mean values and the profiles of the distributions exhibit significant variations with age and position. We confirm that the oldest ($>11$\,Gyr) stars have nearly flat radial [Fe/H] gradients at all heights above the disc but show negative vertical [Fe/H] gradients. For stars younger than 11\,Gyr, the radial [Fe/H] gradients flatten with $|Z|$, while the vertical [Fe/H] gradients flatten with $R$. Stars of 4--6\,Gyr exhibit steeper negative radial [Fe/H] gradients than those of either younger or older ages. Values of [$\alpha$/Fe] of mono-age stellar populations also show significant radial and vertical gradients, with patterns varying with age. The [Fe/H] distribution profiles of old ($>8$\,Gyr) stars vary little with $R$, while those of younger stars exhibit strong radial variations, probably a consequence of significant radial migration. The [$\alpha$/Fe] radial distribution profiles show opposite patterns of variations with age compared to those of [Fe/H]. We have also explored the impacts of stellar mixing by epicycle motions (blurring) on the [Fe/H] and [$\alpha$/Fe] distributions, and found that blurring mainly change the widths of the distribution profiles. Our results suggest that the disc may have experienced a complex assemblage history, in which both the "inside-out" and "upside-down" formation processes may have played an important role.

The LAMOST-Kepler (LK-) project was initiated to use the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) to make spectroscopic follow-up observations for the targets in the field of the Kepler mission. The Kepler field is divided into 14 subfields that are adapted to the LAMOST circular field with a diameter of 5 degrees. During the regular survey phase of LAMOST, the LK-project took data from 2012 June to 2017 June and covered all the 14 subfields at least twice. In particular, we describe in this paper the second Data Release of the LK-project, including all spectra acquired through 2015 May to 2017 June together with the first round observations of the LK-project from 2012 June to 2014 September. The LK-project now counts 227 870 spectra of 156 390 stars, among which we have derived atmospheric parameters (log g, T eff and [Fe/H]) and heliocentric radial velocity (RV) for 173 971 spectra of 126 172 stars. These parameters were obtained with the most recent version of the LAMOST Stellar Parameter Pipeline v 2.9.7. Nearly one half, namely 76 283 targets, are observed both by LAMOST and Kepler telescopes. These spectra, establishing a large spectroscopy library, will be useful for the entire astronomical community, particularly for planetary science and stellar variability on Kepler targets.

SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $\gamma$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $\gamma$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss.

We find significantly different diagnostic emission line ratios for the circumgalactic gas associated with galaxies of stellar masses above and below $10^{10.4}$ M$_\odot$ using SDSS spectroscopy. Specifically, in a sample of 17,393 galaxies, intersected by 18,535 lines of sight at projected radii between 10 and 50 kpc, we stack measured fluxes for nebular strong emission lines, [O {\small III}] $\lambda$5007, H$\alpha$ and [N {\small II}] $\lambda6583$, and find that the gas surrounding the lower mass galaxies exhibits similar line ratios to those of gas ionized by star formation and that surrounding the higher mass galaxies similar to those of gas ionized by AGN or shocks. This finding highlights yet another characteristic of galaxies that is distinctly different above and below this stellar mass threshold, but one that is more closely connected to the gas accretion behavior hypothesized to be responsible for this dichotomy.

The composition of the astrophysical relativistic jets remains uncertain. By kinetic particle-in-cell simulations, we show that the baryon component in the jet, or the so-called baryon loading effect (BLE), heavily affects relativistic jets transport dynamics in the interstellar medium. On the one hand, with the BLE, relativistic jets can transport in a much longer distance, because jet electrons draw a significant amount of energy from jet baryons via the Buneman-induced electrostatic waves and the Weibel-mediated collisionless shock; on the other hand, the jet electron phase space distribution may transform from a bottom-wide-single-peak structure to a center-wide-multiple-peak one by increasing the BLE, which largely influences the observed jet morphology. Implications for related astrophysical studies are also discussed.

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

Taking the neutrino oscillation data into consideration, a dimensionless parameter $\Delta = (m_3-m_1)/(m_3+m_1)$ is adopted to parameterize the three neutrino mass eigenstates and the normal (positive $\Delta$) or inverted (negative $\Delta$) mass hierarchies in three typical cosmological models. Using the currently available cosmic observational data, several Markov Chain Monte Carlo chains are obtained with uniform priors on the free parameters at first. Applying importance sampling the results are compared with three new priors, i.e., logarithmic prior on $|\Delta|$, linear and logarithmic priors on $\Sigma m_\nu$. It turns out that the three new priors increase the upper limits of neutrino mass, but do not change the tendency towards different model's preference for different hierarchies, i.e., the normal hierarchy tends to be favored by $\Lambda$CDM and $w$CDM, which, however, disappears in the $w_0 w_a$CDM model. In addition, the almost symmetrical contours in the $w-\Delta$, $w_0-\Delta$, $w_a-\Delta$ planes indicate that the normal and inverted hierarchy have strong degeneracy. Finally, we perform a Bayesian model comparison analysis, finding that flat linear prior on $\Delta$ and $w_0 w_a$CDM are the most preferred prior and model, respectively.

4C +49.22 is a gamma-ray flat spectrum radio quasar with a bright and knotty jet. We investigate the properties of the core and large-scale knots by using their spectral energy distributions (SEDs). Analyzing its Fermi/LAT data in the past 8 years, a long-term steady gamma-ray emission component is found besides bright outbursts. For the core region, the gamma-ray emission together with the simultaneous emission in the low-energy bands at different epochs is explained with the single-zone leptonic model. The derived magnetization parameters and radiation efficiencies of the radio-core jet decrease as gamma-ray flux decays, likely indicating that a large part of the magnetic energy is converted to the kinetic energy of particles in pc-scale. For the large-scale knots, their radio-optical-X-ray SEDs can be reproduced with the leptonic model by considering the inverse Compton scattering of cosmic microwave background photons. The sum of the predicted gamma-ray fluxes of these knots is comparable to that observed with LAT at 10^{24} Hz of the steady gamma-ray component, indicating that the steady gamma-ray emission may be partially contributed by these large-scale knots. This may conceal the flux variations of the low-level gamma-ray emission from the radio-core. The derived bulk Lorentz factors of the knots decrease along the distance to the core, illustrating as deceleration of jet in large-scale. The powers of the core and knots are roughly in the same order, but the jet changes from highly magnetized at the core region into particle-dominated at the large-scale knots.

A population of so-called "blueberry" starbursting dwarf galaxies with extremely blue colors, low-metallicities, and enormous ionization ratios, has recently been found by Yang et al. (2017). Yet we still do not know their detailed properties, such as morphologies, AGN occupations, massive star contents, infrared emission, dust properties, etc. As a pilot study of the blueberries, we investigate the spectroscopic and photometric properties of two blueberry candidates, RGG B and RGG 5, for which Hubble Space Telescope high-resolution images are available. We find that RGG B and RGG 5 perhaps are likely to be two merging dwarf galaxy systems. RGG B may have a close merging companion; yet the current evidence still cannot exclude the possibility that RGG B is just disturbed by in-situ star formation through, e.g., outflows, rather than undergoing a merger. RGG 5 presents stellar shells in the outskirt which can be the powerful evidence of galaxy merging. We also find that, all of the blueberries, including RGG B and RGG 5, are located close to the theoretical maximum-starburst-line in the BPT diagram, have very high ionization parameters, and relatively low hardness ionizing radiation fields, exhibit nitrogen overabundances, and show extremely red mid-IR colors, and reside in the so-called "ULIRGs/LINERs/Obscured AGN" region. The blueberry galaxies may not harbor AGN.

During the evolution of proto-planetary disc, photo-evaporations of both central and external stars play important roles. Considering the complicated radiation surroundings in the clusters, where the star formed, the proto-planetary discs survive in different lifetimes due to flyby events. In this paper, we mainly focus on the disc around a T Tauri star, which encounters with another main-sequence star with different temperatures in hyperbolic orbits with different peri-center distances, eccentricities and inclinations. We find the criterion for gap-opening due to photo-evaporation of central star after the flyby event. A gap is opened in the late stage of gas disc, and induce that the gap only influence the planet formation and migration limitedly. If the flyby orbit has a moderate value of peri-center distance, which weakly depends on the eccentricity and inclination, the external photo-evaporation lead to a maximum mass loss during the flyby event. Flyby stars in orbits with smaller eccentricities or larger inclinations induce larger mass loss. Adopting a simple multiple flyby models, we conclude that in open clusters, gas discs usually survive in typical lifetimes between 1 and 10 Myr, except there are many massive stars in dense open clusters. In globular clusters, discs disperse very quickly and hardly produce the gas giant planets. The fast-depleted discs are probably responsible for the null detection of giant planets in globular clusters.

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

We develop a formalism that treats the calculation of solar kinetic energy and cross helicity spectra in an equal manner to that of magnetic energy and helicity spectra. The magnetic helicity spectrum is shown to be equal to the vertical part of the current helicity spectrum divided by the square of the wavenumber. For the cross helicity, we apply the recently developed two-scale approach globally over an entire active region to account for the sign change between the two polarities. Using vector magnetograms and Dopplergrams of NOAA 11158 and 12266, we show that kinetic and magnetic energy spectra have similar slopes at intermediate wavenumbers, where the contribution from the granulation velocity has been removed. At wavenumbers around 0.3 Mm$^{-1}$, the magnetic helicity is found to be close to its maximal value. The cross helicity spectra are found to be within about 10% of the maximum possible value. Using the two-scale method for NOAA 12266, the global cross helicity spectrum is found to be particularly steep, similarly to what has previously been found in theoretical models of spot generation. In the quiet Sun, by comparison, the cross helicity spectrum is found to be small.

This study expands the coverage and improves the homogeneity of the distribution of MILES template stars in the parameter space, as well as extends the wavelength coverage of the template spectra to the far red beyond the Ca\,{\sc{ii}} triplet. To achieve this we have carried out a major observational campaign using the OMR long-slit spectrograph mounted on the NAOC 2.16\,m telescope and the YFOSC long-slit spectrograph mounted on the YNAO 2.4\,m telescope. The original sample is based on the MILES library, supplemented by 918 stars selected from PASTEL database. In total, 822 OMR and 1,324 YFOSC spectra, covering respectively the wavelength ranges $\lambda\lambda$3800--5180 and $\lambda\lambda$5150--9000, have been collected and reduced. The spectra have a mean resolution FWHM (full-width at half-maximum) of $\sim 3.3$\,{\AA} and are wavelength- and flux-calibrated to an accuracy of $\sim 20$\,km\,s$^{-1}$ and $\sim 5$\,per\,cent, respectively. The spectra are further corrected for systematic errors in the wavelength calibration to an accuracy of $\sim 4$\,km\,s$^{-1}$ by cross-correlating with the theoretical spectra. Almost all the spectra have an average signal to noise ratio (SNR) better than 100 per pixel. Combined with the MILES spectra, there are now 1,731, 1,542, 1,324 and 1,273 stars with spectra covering respectively $\lambda\lambda$3800--5180, $\lambda\lambda$3800--7500, $\lambda\lambda$5150--9000 and $\lambda\lambda$3800--9000. This paper describes our template star selection, the observation and data reduction, and presents the reduced spectra collected hitherto.

As a key process that refreshes the interstellar medium, the dynamics and radiative properties of the supernova remnant (SNR) expansion front not only reflect the physical environment of the old interstellar medium (ISM) surrounding the supernova, but they also provide information about the refreshed ISM. However the expansion dynamics of SNRs cannot be simply explained by the conventional law of spherical shock wave propagation; on the other hand, the high energy radiation requires an additional electron acceleration mechanism in the shock front beyond thermal collision. We consider herein the detonation wave description of the SNR expansion, in which magnetic reconnection follows the shock front and transfers the SNR magnetic field energy to both fluid thermal energy and particle kinetic energy. The structure of the magnetic reconnection detonation (MRD) is identified based on scaling analysis in this paper. By applying the MRD description of the SNR expansion shock to the example of the Crab Nebula, this paper shows that the MRD description can explain both the accelerative expansion of the nebula as well as the origin of the luminous expanding shell.

We analyse the spatial distribution of current helicity in solar active regions. A comparison of current helicity maps derived from three different instruments (Helioseismic and Magnetc Imager on board the Solar Dynamics Observatory, SDO/HMI, Spectro-Polarimeter on board the Hinode, and Solar Magnetic Field Telescope at the Huairou Solar Observing Station, China, HSOS/SMFT) is carried out. The comparison showed an excellent correlation between the maps derived from the spaceborne instruments and moderate correlation between the maps derived from SDO/HMI and HSOS/SMFT vector magnetograms. The results suggest that the obtained maps characterize real spatial distribution of current helicity over an active region. To analyse intermittency of current helicity, we traditionally use the high-order structure functions and flatness function approach. The slope of a flatness function within some range of scales - the flatness exponent - is a measure of the degree of intermittency. SDO/HMI vector magnetograms for 3 ARs (NOAA 11158, 12494, and 12673) were used to calculate the flatness exponent time variations. All three ARs exhibited emergence of a new magnetic flux during the observational interval. The flatness exponent indicated the increase of intermittency 12-20 hours before the emergence of a new flux. We suppose that this behaviour can indicate subphotospheric fragmentation or distortion of the pre-existed electric current system by emerging magnetic flux.

The association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. We perform semi-analytical spectral fitting of the multi-wavelength emission to obtain estimates of the jet physical parameters, and find that the multi-wavelength emission can be explained by either a proton synchrotron scenario or an electron inverse Compton scattering scenario. In the proton synchrotron scenario, a strong magnetic field of $10-100$~G is required, implying that the particle acceleration is likely driven by magnetic energy dissipation such as magnetic reconnection events. The inverse Compton scenario implies a magnetic field of $0.1-1$~G. Thus the particle acceleration is likely driven by the kinetic energy dissipation such as shocks. We also discuss the neutrino production in the context of single-zone and multi-zone models based on the above two scenarios. We demonstrate that the variability and optical polarization signatures can be used to distinguish the two scenarios due to their drastically different magnetic field. Specifically, the proton synchrotron scenario may show orphan fast variability in the low-energy spectral component on top of the active state, with an optical polarization degree $\lesssim 10\%$ throughout the active state. The inverse Compton scattering scenario instead predicts co-variability of the low- and high-energy components on both short and long time scales, as well as a strongly variable optical polarization degree that can reach $\gtrsim 20\%$. Our results suggest that optical polarization measurements and well-resolved multi-wavelength light curves can be used to understand the electromagnetic and high-energy neutrino emissions by TXS~0506+056 and similar events in the future.

Polar Crown Filaments (PCFs) form above the magnetic polarity inversion line, which separates the unipolar polar fields and the nearest dispersed fields from trailing part of active regions with opposite polarity. The statistical properties of PCFs are correlated with the solar cycle. Therefore, study of PCFs plays an important role in understanding the variation of solar cycle, especially the prolonged cycle 23 and the current "abnormal" solar cycle 24. In this study, we investigate PCFs using full disk H{\alpha} data from 1973 to early 2018, recorded by Kanzelhohe Solar Observatory (KSO) and Big Bear Solar Observatory (BBSO), in digital form from 1997 to 2018 and in 35 mm film (digitized) from 1973 to 1996. PCFs are identified manually because their segmented shape and close-to-limb location were not handled well by automatical detections in several previous studies. Our results show that the PCFs start to move poleward at the beginning of each solar cycle. When the PCFs approach to the maximum latitude, the polar field strength reduces to zero followed by a reversal. The migration rates are about 0.4 to 0.7 degree per Carrington rotation, with clear N-S asymmetric pattern. In cycles 21 and 23, the PCFs in the northern hemisphere migrate faster than those in the southern hemisphere. However, in the "abnormal" cycle 24, the southern PCFs migrate faster, which is consistent with other observations of magnetic fields and radio emission. In addition, there are more days in cycle 23 and 24 without PCFs than in the previous cycles.