In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

The recent CHIME/FRB and STARE2 detection of an extremely bright (fluence ~MJy ms) radio burst from the Galactic magnetar SGR 1935+2154 supports the hypothesis that (some) fast radio bursts (FRBs) are emitted by magnetars at cosmological distances. With the aim of detecting more bursts from SGR 1935+2154, we performed monitoring observations of the source using (up to) four small radio telescopes observing simultaneously at complementary radio frequencies, spanning the range from 327 MHz to 8.6 GHz. In close to four weeks of almost daily monitoring (265 hrs on source) we detected two bright radio bursts that are separated in time by only 1.4 seconds, i.e. separated by 43% of the 3.25-s rotational period of this neutron star. The bursts have fluences of 112+/-22 Jy ms and 24+/-5 Jy ms, respectively, and we measure a scattering time scale of 313+/-31 us at 1.32 GHz. Given the proximity of the bursts in time, we estimate the shape parameter k and rate parameter r of an assumed Weibull distribution and find k significantly less than one. This suggests a non-Poissonian, clustered emission process and is consistent with what has been measured for FRB 121102. Together with the CHIME/FRB and STARE2 burst, as well as a faint (fluence 60 mJy ms) burst reported by FAST, these observations demonstrate that SGR 1935+2154 can produce bursts with apparent energies spanning roughly 7 orders of magnitude. This raises the question whether they all arise from similar physical processes and whether the FRB population distribution extends to very low energies (~10^30 erg, isotropic equivalent).

Are there hurricanes on exoplanets? Tidally locked terrestrial planets around M dwarfs are the main targets of space missions for finding habitable exoplanets. Whether hurricanes can form on this kind of planet is important for their climate and habitability. Using a high-resolution global atmospheric circulation model, we investigate whether there are hurricanes on tidally locked terrestrial planets under fixed sea surface temperatures (SST). The effects of planetary rotation rate, surface temperature, and bulk atmospheric compositions are examined. We find that hurricanes can form on the planets but not on all of them. For planets near the inner edge of the habitable zone of late M dwarfs, there are more and stronger hurricanes on both day and night sides. For planets in the middle and outer ranges of the habitable zone, the possibility of hurricane formation is low or even close to zero, as suggested in the study of Bin et al.(2017). Hurricane theories on Earth are applicable to tidally locked planets when atmospheric compositions are similar to Earth. However, if background atmosphere is lighter than H2O, hurricanes can hardly be produced because convection is always inhibited due to the effect of mean molecular weight, similar to that on Saturn. These results have broad implications on the precipitation, ocean mixing, climate, and atmospheric characterization of tidally locked planets. Finally, A test with a coupled slab ocean and an earth-like atmosphere in a tide-locked orbit of 10 earth days shows that there are also hurricanes in the experiment.

Active galactic nuclei (AGNs) have been attracting research attention due to their special observable properties. Specifically, a majority of AGNs are detected by Fermi-LAT missions, but not by Fermi-LAT, which raises the question of whether any differences exist between the two. To answer this issue, we compile a sample of 291 superluminal AGNs (189 FDSs and 102 non-FDSs) from available multi-wavelength radio, optical, and X-ray (or even $\gamma$-ray) data and Doppler factors and proper motion ($\mu$) (or apparent velocity ($\beta_{\rm{app}}$)); calculated the apparent velocity from their proper motion, Lorentz factor ($\Gamma$), viewing angle ($\phi$) and co-moving viewing angle ($\phi_{co}$) for the sources with available Doppler factor ($\delta$); and performed some statistical analyses for both types. Our study indicated that1. In terms of average values, FDSs have higher proper motions ($\mu$), apparent velocities ($\beta_{\rm app}$), Doppler factor ($\delta$), Lorentz factor ($\Gamma$), and smaller viewing angle ($\phi$). Nevertheless, there is no clear difference in co-moving viewing angles ($\phi_{\rm co}$).

Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$\sigma$. The ten years of data post-trial significance is 4.1$\sigma$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

Very recently, the Xenon1T collaboration has reported an intriguing electron recoil excess, which may imply for light dark matter. In order to interpret this anomaly, we proposed the the atmospheric dark matter (AMD) from the inelastic collision of cosmic rays with the atmosphere. Thanks to the boost effect of high energy CRs, we show that a ${\cal O}(10)$ keV ADM can be fast-moving and successfully fit the observed electron recoil spectra observed through the ADM-electron scattering process. Meanwhile, our ADM predict the scattering cross sections $\sigma_e \sim {\cal O}(10^{-37}- 10^{-38}$) cm$^{2}$, and thus can evade other direct detection constraints. The search for light meson rare decays, such as $\eta \to \pi + \slashed E_T$, would provide a complementary probe of our ADM in the future.

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

We present an analysis of local molecular clouds (-6 <VLSR< 30 km/s, i.e., <1.5 kpc) in the first Galactic quadrant (25.8{\deg} <l<49.7{\deg} and |b|<5{\deg}), a pilot region of the Milky Way Imaging Scroll Painting (MWISP) CO survey. Using the SCIMES algorithm to divide large molecular clouds into moderate-sized ones, we determined distances to 28 molecular clouds with the background-eliminated extinction-parallax (BEEP) method using the Gaia DR2 parallax measurements aided by AG and AV, and the distance ranges from 250 pc to about 1.5 kpc. These incomplete distance samples indicate a linear relationship between the distance and the radial velocity (VLSR) with a scatter of 0.16 kpc, and kinematic distances may be systematically larger for local molecular clouds. In order to investigate fundamental properties of molecular clouds, such as the total sample number, the linewidth, the brightness temperature, the physical area, and the mass, we decompose the spectral cube using the DBSCAN algorithm. Post selection criteria are imposed on DBSCAN clusters to remove the noise contamination, and we found that the separation of molecular cloud individuals is reliable based on a definition of independent consecutive structures in l-b-V space. The completeness of the local molecular cloud flux collected by the MWISP CO survey is about 80%. The physical area, A, shows a power-law distribution, dN/dA \propto A^{-2.20+/-0.18}, while the molecular cloud mass also follows a power-law distribution but slightly flatter, dN/dM \propto M^{-1.96+/-0.11}.

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

From near-infrared spectroscopic measurements of the MgII emission line doublet, we estimate the black hole (BH) mass of the quasar, SMSS J215728.21-360215.1, as being (3.4 +/- 0.6) x 10^10 M_sun and refine the redshift of the quasar to be z=4.692. SMSS J2157 is the most luminous known quasar, with a 3000A luminosity of (4.7 +/- 0.5) x 10^47 erg/s and an estimated bolometric luminosity of 1.6 x 10^48 erg/s, yet its Eddington ratio is only ~0.4. Thus, the high luminosity of this quasar is a consequence of its extremely large BH -- one of the most massive BHs at z > 4.

The recent measurements of the cosmological parameter $H_0$ from the direct local observations and the inferred value from the Cosmic Microwave Background show $\sim 4 \sigma$ discrepancy. We demonstrate that a keV gravitino dark matter, which has a small fraction of non-thermal component (e.g. from the late decay of NLSP bino) may reconcile this tension. Furthermore, we discuss the implied collider signatures and point out that this scenario can be tested by searching for the dilepton plus missing energy events at the LHC in near future.

Tidally locked terrestrial planets around low-mass stars are the prime targets of finding potentially habitable exoplanets. Several atmospheric general circulation models have been employed to simulate their possible climates, however, model intercomparisons showed that there are large differences in the results of the models even when they are forced with the same boundary conditions. In this paper, we examine whether model resolution contributes to the differences. Using the atmospheric general circulation model ExoCAM coupled to a 50-m slab ocean, we examine three different horizontal resolutions (440 km * 550 km, 210 km * 280 km, and 50 km * 70 km in latitude and longitude) and three different vertical resolutions (26, 51, and 74 levels) under the same dynamical core and the same schemes of radiation, convection and clouds. Among the experiments, the differences are within 5 K in global-mean surface temperature and within 0.007 in planetary albedo. These differences are from cloud feedback, water vapor feedback, and the decreasing trend of relative humidity with increasing resolution. Relatively small-scale downdrafts between upwelling columns over the substellar region are better resolved and the mixing between dry and wet air parcels and between anvil clouds and their environment are enhanced as the resolution is increased. These reduce atmospheric relative humidity and high-level cloud fraction, causing a lower clear-sky greenhouse effect, a weaker cloud longwave radiation effect, and subsequently a cooler climate with increasing model resolution. Overall, the sensitivity of the simulated climate of tidally locked aquaplanets to model resolution is small.

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

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

We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function.

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

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

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$\alpha$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

Data of $^{12}$CO/$^{13}$CO/C$^{18}$O $J=1\to0$ emission toward the Galactic plane region of $l=35^\circ$ to $45^\circ$ and $b= -5^\circ$ to $+5^\circ$ are available with the Milky Way Imaging Scroll Painting (MWISP) project. Using the data, we found a giant molecular filament (GMF) around $l\approx38\sim42^\circ$, $b\approx-3.5\sim0^\circ$, $V_{LSR} \approx 27 \sim 40$ km~s$^{-1}$, named the GMF MWISP G041-01. At a distance of 1.7 kpc, the GMF is about 160 pc long. With a median excitation temperature about 7.5 K and a median column density about $10^{21}$ cm$^{-2}$, this GMF is very cold and very diffuse compared to known GMFs. Using the morphology in the data cube, the GMF is divided into four components among which three show filamentary structure. Masses of the components are $ 10^3 \sim 10^4 M_\odot$, with a total mass for the whole filament being about $7\times10^4 M_\odot$ from the LTE method. $^{13}$CO cores inside each component are searched. Virial parameters are about 2.5 for these cores and have a power-law index of -0.34 against the mass. The mass fraction of dense cores traced by $^{13}$CO to the diffuse clouds traced by $^{12}$CO are about 7% for all components of the GMF. We found signatures of possible large scale filament-filament collision in the GMF.

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

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.

Site-testing is crucial for achieving the goal of scientific research and analysis of meteorological and optical observing conditions is one of the basic tasks of it. As one of three potential sites to host 12-meter Large Optical/infrared Telescope (LOT), Muztagh-ata site which is located on the Pamirs Plateau in west China's Xinjiang began its site-testing task in the spring of 2017. In this paper, we firstly start with an introduction to the site and then present a statistical analysis of the ground-level meteorological properties such as air temperature, barometric pressure, relative humidity, wind speed and direction, recorded by automatic weather station with standard meteorological sensors for two-year long. We also show the monitoring results of sky brightness during this period.

In this article, we present a detailed analysis of the statistical properties of seeing for the Muztagh-ata site which is the candidate site for hosting future Chinese Large Optical/infrared Telescope (LOT) project. The measurement was obtained with Differential Image Motion Monitor (DIMM) from April 2017 to November 2018 at different heights during different periods. The median seeing at 11 meters and 6 meters are very close but different significantly from that on the ground. We mainly analyzed the seeing at 11 meters monthly and hourly, having found that the best season for observing was from late autumn to early winter and seeing tended to improve during the night only in autumn. The analysis of the dependence on temperature inversion, wind speed, direction also was made and the best meteorological conditions for seeing is given.

The Kepler mission has discovered thousands of exoplanets around various stars with different spectral types (M, K, G, and F) and thus different masses and effective temperatures. Previous studies have shown that the planet occurrence rate, in terms of average number of planets per star, drops with increasing stellar effective temperature (Teff). In this paper, with the final Kepler Data Release (DR25) catalog, we revisit the relation between stellar effective temperature (as well as mass) and planet occurrence, but in terms of the fraction of stars with planets and the number of planets per planetary system (i.e. planet multiplicity). We find that both the fraction of stars with planets and planet multiplicity decrease with increasing stellar temperature and mass. Specifically, about 75% late-type stars (Teff<5000 K) have Kepler-like planets with an average planet multiplicity of ~2.8; while for early-type stars (Teff>6500 K), this fraction and the average multiplicity fall down to ~35% and ~1.8, respectively. The decreasing trend in the fraction of stars with planets is very significant with ${\Delta}$AIC> 30, though the trend in planet multiplicity is somewhat tentative with ${\Delta}$AIC~5. Our results also allow us to derive the dispersion of planetary orbital inclinations in relationship with stellar effective temperature. Interestingly, it is found to be similar to the well-known trend between obliquity and stellar temperature, indicating that the two trends might have a common origin.

Highly accreting quasars are quite luminous in the X-ray and optical regimes. While, they tend to become radio quiet and have optically thin radio spectra. Among the known quasars, IRAS F11119+3257 is a supercritical accretion source because it has a bolometric luminosity above the Eddington limit and extremely powerful X-ray outflows. To probe its radio structure, we investigated its radio spectrum between 0.15 and 96.15 GHz and performed very-long-baseline interferometric (VLBI) observations with the European VLBI Network (EVN) at 1.66 and 4.93 GHz. The deep EVN image at 1.66 GHz shows a two-sided jet with a projected separation about two hundred parsec and a very high flux density ratio of about 290. Together with the best-fit value of the integrated spectral index of -1.31+/-0.02 in the optically thin part, we infer that the approaching jet has an intrinsic speed at least 0.57 times of the light speed. This is a new record among the known all kinds of super-Eddington accreting sources and unlikely accelerated by the radiation pressure. We propose a scenario in which IRAS F11119+3257 is an unusual compact symmetric object with a small jet viewing angle and a radio spectrum peaking at 0.53+/-0.06 GHz mainly due to the synchrotron self-absorption.

The population of intermediate-mass black holes (IMBHs) in nearby dwarf galaxies plays an important "ground truth" role in exploring black hole formation and growth in the early Universe. In the dwarf elliptical galaxy SDSS J090613.77+561015.2 (z=0.0465), an accreting IMBH has been revealed by optical and X-ray observations. Aiming to search for possible radio core and jet associated with the IMBH, we carried out very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.66 GHz. Our imaging results show that there are two 1-mJy components with a separation of about 52 mas (projected distance 47 pc) and the more compact component is located within the 1-sigma error circle of the optical centroid from available Gaia astrometry. Based on their positions, elongated structures and relatively high brightness temperatures, as well as the absence of star-forming activity in the host galaxy, we argue that the radio morphology originates from the jet activity powered by the central IMBH. The existence of the large-scale jet implies that violent jet activity might occur in the early epochs of black hole growth and thus help to regulate the co-evolution of black holes and galaxies.

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

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

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

We present an observational study of 22 broad absorption line quasars (BAL QSOs) at $3\lesssim z \lesssim5$ based on optical/near-IR spectroscopy, aiming to investigate quasar winds and their effects. The near-IR spectroscopy covers the \hb\ and/or \mgii\ broad emission lines (BELs) for these quasars, allowing us to estimate their central black hole (BH) masses in a robust way. We found that our BAL QSOs on average do not have a higher Eddington ratio than that from non-BAL QSOs matched in redshift and/or luminosity. In a subset consisting of seven strong BAL QSOs possessing sub-relativistic BAL outflows, we see the prevalence of large \civ-BEL blueshift ($\sim$3100 km s$^{-1}$) and weak \oiii\ emission (particularly the narrow \oiii$\lambda$5007 component), indicative of nuclear outflows affecting the narrow emission-line (NEL) regions. In another subset consisting of thirteen BAL QSOs having simultaneous observations of \mgii\ and \hb, we found a strong correlation between 3000~\AA\ and 5000~\AA\ monochromatic luminosity, consistent with that from non-BAL QSOs matched in redshift and luminosity; however, there is no correlation between \mgii\ and \hb\ in FWHM, likely due to nuclear outflows influencing the BEL regions. Our spectroscopic investigations offer strong evidence that the presence of nuclear outflows plays an important role in shaping the BEL/NEL regions of these quasars and possibly, regulating the growth of central supermassive black holes (SMBHs). We propose that BEL blueshift and BAL could be different manifestations of the same outflow system viewed at different sightlines and/or phases.

The potential association between the blazar TXS 0506+056 and the neutrino event IceCube-170922A provides a unique opportunity to study the possible physical connection between the high-energy photons and neutrinos. We explore the correlated electromagnetic and neutrino emissions of blazar TXS 0506+056 by a self-consistent leptonic-hadronic model, taking into account particle stochastic acceleration and all relevant radiative processes self-consistently. The electromagnetic and neutrino spectra of blazar TXS 0506+056 are reproduced by the proton synchrotron and hybrid leptonic-hadronic models based on the proton-photon interactions. It is found that the hybrid leptonic-hadronic model can be used to better explain the observed X-ray and $\gamma$-ray spectra of blazar TXS 0506+056 than the proton synchrotron model. Moreover, the predicted neutrino spectrum of the hybrid leptonic-hadronic model is closer to the observed one compared to the proton synchrotron model. We suggest that the hybrid leptonic-hadronic model is more favored if the neutrino event IceCube-170922A is associated with the blazar TXS 0506+056.

We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well.

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.

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

The Large Optical/infrared Telescope (LOT) is a ground-based 12m diameter optical/infrared telescope which is proposed to be built in the western part of China in the next decade. Based on satellite remote sensing data, along with geographical, logistical and political considerations, three candidate sites were chosen for ground-based astronomical performance monitoring. These sites include: Ali in Tibet, Daocheng in Sichuan, and Muztagh Ata in Xinjiang. Up until now, all three sites have continuously collected data for two years. In this paper, we will introduce this site testing campaign, and present its monitoring results obtained during the period between March 2017 and March 2019.

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

The fast-rising blue optical transient AT2018cow indicated unusual early phase characteristics unlike relatively better studied explosive transients. Its afterglow may be produced by either a relativistically beamed (jetted) or intrinsically luminous (non-jetted) ejecta and carries observational signatures of the progenitor and environment. High resolution monitoring can distinguish between these scenarios and clarify the progenitor nature. We present very long baseline interferometry (VLBI) observations of AT2018cow at 5 GHz involving 21 radio telescopes from the European VLBI Network with five sessions spanning ~ 1 year. With an astrometric precision up to 25 micro-arcseconds per epoch, the rapidly fading compact milli arcsecond scale source is found to be non-jetted with a proper motion of <= 0.15 mas/yr (0.14 c). This and a dense (number density ~ 10^4 - 10^5/cm^3) magnetized environment (magnetic field strength >= 0.84 G) are characteristic of a newly formed magnetar driven central engine, originating in the successful explosion of a low-mass star.

In many particle physics models, domain wall can form during the phase transition process after discrete symmetry breaking. We study the scenario within a complex singlet extended Standard Model framework, where a strongly first order phase transition can occur depending on the hidden scalar mass and the mixing between the extra heavy Higgs and the SM Higgs mass. The gravitational wave spectrum is of a typical two-peak shape, the amplitude and the peak from the strongly first order phase transition is able to be probed by the future space-based interferometers, and the one locates around the peak from the domain wall decay is far beyond the capability of the current PTA, and future SKA.

Blazars are a sub-class of quasars with Doppler boosted jets oriented close to the line of sight, and thus efficient probes of supermassive black hole growth and their environment, especially at high redshifts. Here we report on Very Long Baseline Interferometry observations of a blazar J0906+6930 at z = 5.47, which enabled the detection of polarised emission and measurement of jet proper motion at parsec scales. The observations suggest a less powerful jet compared with the general blazar population, including lower proper motion and bulk Lorentz factor. This coupled with a previously inferred high accretion rate indicate a transition from an accretion radiative power to a jet mechanical power based transfer of energy and momentum to the surrounding gas.While alternative scenarios could not be fully ruled out, our results indicate a possibly nascent jet embedded in and interacting with a dense medium resulting in a jet bending.

Tidally locked terrestrial planets around low-mass stars are the prime targets for future atmospheric characterizations of potentially habitable systems, especially the three nearby ones--Proxima b, TRAPPIST-1e, and LHS 1140b. Previous studies suggest that if these planets have surface ocean they would be in an eyeball-like climate state: ice-free in the vicinity of the substellar point and ice-covered in the rest regions. However, an important component of the climate system--sea ice dynamics has not been well studied in previous studies. A fundamental question is: would the open ocean be stable against a globally ice-covered snowball state? Here we show that sea-ice drift cools the ocean's surface when the ice flows to the warmer substellar region and melts through absorbing heat from the ocean and the overlying air. As a result, the open ocean shrinks and can even disappear when atmospheric greenhouse gases are not much more abundant than on Earth, turning the planet into a snowball state. This occurs for both synchronous rotation and spin-orbit resonances (such as 3 to 2). These results suggest that sea-ice drift strongly reduces the open ocean area and can significantly impact the habitability of tidally locked planets.

Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an inter-comparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (20-30 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance.

How does subatomic matter organize itself? Neutron stars are cosmic laboratories uniquely poised to answer this fundamental question that lies at the heart of nuclear science. Newly commissioned rare isotope facilities, telescopes operating across the entire electromagnetic spectrum, and ever more sensitive gravitational wave detectors will probe the properties of neutron-rich matter with unprecedented precision over an enormous range of densities. Yet, a coordinated effort between observation, experiment, and theoretical research is of paramount importance for realizing the full potential of these investments. Theoretical nuclear physics provides valuable insights into the properties of neutron-rich matter in regimes that are not presently accessible to experiment or observation. In particular, nuclear density functional theory is likely the only tractable framework that can bridge the entire nuclear landscape by connecting finite nuclei to neutron stars. This compelling connection is the main scope of the present review.

Detection of non-gravitational interactions of massive dark matter with visible sector so far have given null results. The DM may communicate with the ordinary matter only through gravitational interactions. Besides, the majority of traditional direct detections have poor sensitivity of searching for light DM particle because of the small recoil energy. Thanks to the high energy cosmic rays (CRs), the light DM can be boosted by scattering with CRs and thus may be detected in those existing experiments. In this work, we first derive the exclusion limits of the cosmic ray boosted sub-GeV Dark Matter with gravitational mediator from the Xenon1T data.

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.

Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

We have conducted a large-field simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ emission toward the Orion A giant molecular cloud (GMC) with a sky coverage of $\sim$ 4.4 deg$^2$ using the PMO-13.7 m millimeter-wavelength telescope. We use the probability distribution function of the column density (N-PDF) to investigate the distribution of molecular hydrogen in the Orion A GMC. The H$_2$ column density, derived from the $^{13}$CO emission, of the GMC is dominated by log-normal distribution in the range from $\sim$4$\times10^{21}$ to $\sim$1.5$\times10^{23}$ cm$^{-2}$ with excesses both at the low-density and high-density ends. The excess of the low-density end is possibly caused by an extended and low-temperature ($\sim$10 K) component with velocities in the range of 5$-$8 km s$^{-1}$. Compared with the northern sub-regions, the southern sub-regions of the Orion A GMC contain less gas with column density in $N_{H_2} > 1.25\times 10^{22}\ \rm{cm}^{-2}$. The dispersions of the N-PDFs of the sub-regions are found to correlate with the evolutionary stages of the clouds across the Orion A GMC. The structure hierarchy of Orion A GMC is explored with the DENDROGRAM algorithm, and it is found that the GMC is composed of two branches. All structures except one in the tree have virial parameters less than 2, indicating self-gravity is important on the spatial scales from $\sim$0.3 to $\sim$4 pc. Although power-laws and departures from log-normal distributions are found at the high-density end of N-PDFs of active star-forming regions, the N-PDFs of structures in the Orion A GMC are predominantly log-normal on scales from R$\sim$0.4 to 4 pc.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

Axions are a promising cold dark matter candidate. Haloscopes, which use the conversion of axions to photons in the presence of a magnetic field to detect axions, are the basis of microwave cavity searches such as the Axion Dark Matter eXperiment (ADMX). To search for lighter, low frequency axions in the sub $2\times10^{-7}$ eV (50 MHz) range, a tunable lumped-element LC circuit has been proposed. For the first time, through ADMX SLIC (Superconducting Lc circuit Investigating Cold axions), a resonant LC circuit was used to probe this region of axion mass-coupling space. The detector used a superconducting LC circuit with piezoelectric driven capacitive tuning. The axion mass and corresponding frequency range $1.7498 -1.7519 \times10^{-7}$ eV (42.31 -- 42.36 MHz), $1.7734 - 1.7738 \times10^{-7}$ eV (42.88 -- 42.89 MHz), and $1.8007 - 1.8015 \times10^{-7}$ eV (43.54 -- 43.56 MHz) was covered at magnetic fields of 4.5 T, 5.0 T, and 7.0 T respectively. Exclusion results from the search data, for coupling below $10^{-12} \text{GeV}^{-1}$ are presented.

Using high spatial and temporal data from the New Vacuum Solar Telescope (NVST) and the Solar Dynamics Observatory (SDO), we present unambiguous observations of recurrent two-sided loop jets caused by magnetic reconnection between erupting minifilaments and nearby large filament. The observations demonstrate that three two-sided loop jets, which ejected along the large filament in opposite directions, had similar appearance and originated from the same region. We find that a minifilament erupted and drove the first jet. It reformed at the same neutral line later, and then underwent partial and total eruptions, drove the second and third jets, respectively. In the course of the jets, cool plasma was injected into the large filament. Furthermore, persistent magnetic flux cancelation occurred at the neutral line under the minifilament before its eruption and continued until the end of the observation. We infer that magnetic flux cancellation may account for building and then triggering the minifilament to erupt to produce the two-sided loop jets. This observation not only indicates that two-sided loop jets can be driven by minifilament eruptions, but also sheds new light on our understanding of the recurrent mechanism of two-sided loop jets.

We study the early evolution of a hot-channel-like magnetic flux rope (MFR) toward eruption. Combining with imaging observation and magnetic field extrapolation, we find that the hot channel possibly originated from a pre-existing seed MFR with a hyperbolic flux tube (HFT). In the precursor phase, three-dimensional tether-cutting reconnection at the HFT is most likely resulting in the heating and buildup of the hot channel. In this process, the forming hot channel was rapidly enlarged at its spatial size and slipped its feet to two remote positions. Afterward, it instantly erupted outwards with an exponential acceleration, leaving two core dimmings near its feet. We suggest that pre-flare reconnection at the HFT played a crucial role in enlarging the seed MFR and facilitating the onset of its final solar eruption. Moreover, a recently predicted drifting of MFR's footpoints was detected at both core dimmings. In particular, we find that MFR's west footpoint drift was induced by a new reconnection geometry among the erupting MFR's leg and thereby inclined arcades. As MFR's west footpoints gradually drifted to a new position, a set of newborn atypical flare loops connected into the west core dimming, causing a rapid decrease of dimmed area inside this core dimming and also generating a secondary flare ribbon at their remote feet. This reveals that core dimmings may suffer a pronounced diminishment due to the eruptive MFR's footpoint drift, implying that mapping the real footpoints of the erupting MFR down to the Sun's surface is more difficult than previously thought.

This paper reports on a cavity haloscope search for dark matter axions in the galactic halo in the mass range $2.81$-$3.31$ ${\mu}eV$. This search excludes the full range of axion-photon coupling values predicted in benchmark models of the invisible axion that solve the strong CP problem of quantum chromodynamics, and marks the first time a haloscope search has been able to search for axions at mode crossings using an alternate cavity configuration. Unprecedented sensitivity in this higher mass range is achieved by deploying an ultra low-noise Josephson parametric amplifier as the first stage signal amplifier.

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

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

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

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

About 4000 exoplanets have been confirmed since the year of 1992, and for most of the planets, the main parameters that can be measured are planetary radius and mass. Based on these two parameters, surface gravity can be estimated. In this work, we investigate the effects of varying radius and gravity on the climate of rapidly rotating terrestrial planets with assuming they have an ocean and Earth-like atmospheres (N$_2$, CO$_2$, and H$_2$O). Using a three-dimensional (3D) atmospheric general circulation model (GCM), we find that varying radius mainly influences the equator-to-pole surface temperature difference while varying gravity mainly influences the mean surface temperature. For planets of larger radii, the meridional atmospheric energy transport is weaker, which warms the tropics but cools high latitudes. For planets of larger gravities, the surface is globally cooler due to the fact that saturated vapor pressure depends on air temperature only and column water vapor mass is approximately equal to saturated vapor pressure over gravity multiplied by relative humidity. The relative humidity does not change much in our experiments. Ice albedo and water vapor feedbacks act to further amplify the effects of varying radius and gravity. These results suggest that radius and gravity are important factors for planetary climate, although they are not as critical as stellar flux, atmospheric composition and mass, and rotation rate. Future simulations are required to know how large the radius and gravity could influence the width of the habitable zone using 3D GCMs.

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

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

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

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

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

We present a new method of calculating distances of molecular clouds in the Galactic plane, using CO observations and the Gaia DR2 parallax and G-band extinction ($A_G$) measurements. Due to the complexity of dust environments in the Galactic plane, $A_G$ contains irregular variations, which is difficult to model. To overcome this difficulty, we propose that the $A_G$ of off-cloud stars (Gaia stars around molecular clouds) can be used as a baseline to calibrate the $A_G$ of on-cloud stars (Gaia stars toward molecular clouds), which removes the $A_G$ components that are unrelated to molecular clouds. The distance is subsequently inferred from the jump point in on-cloud $A_G$ with Bayesian analysis and Markov chain Monte Carlo (MCMC) sampling. We applied this baseline subtraction method to a 100 square degree region (209.75{\deg} $\leq$ l $\leq$ 219.75{\deg} and |b| $\leq$ 5{\deg}) in the third Galactic quadrant, which was mapped as part of the Milky Way Imaging Scroll Painting (MWISP) project, covering three CO isotopologue lines, and derived distances and masses for 11 molecular clouds, including the Maddalena molecular cloud and Sh 2-287. The results indicate that the distance of the Perseus Arm in this region is about 2.4 kpc and molecular clouds are present in the interarm regions.

Active galactic nuclei (AGNs) can be divided into two major classes, namely radio-loud and radio-quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff-Riley type I and type II (FRI&II) radio galaxies. Following our previous work, we present a latest sample of 966 sources with measured radio flux densities of the core and extended components. The sample includes 83 BL Lacs, 473 FSRQs, 101 Seyferts, 245 galaxies, 52 FRIs&IIs and 12 unidentified sources. We then calculate the radio core-dominance parameters and spectral indices and study their relationship. Our analysis shows that the core-dominance parameters and spectral indices are quite different for different types of sources. We also corroborate that the correlation between core-dominance parameter and radio spectral index extends over all the sources in a large sample presented.

We present observations of the CO isotopologues ($^{12}$CO, $^{13}$CO, and C$^{18}$O) toward the Galactic region with $169\arcdeg.75 \leqslant l \leqslant 174\arcdeg.75$ and $-0\arcdeg.75 \leqslant b \leqslant 0\arcdeg.5$, using the Purple Mountain Observatory 13.7~m millimeter-wavelength telescope. Based on the $^{13}$CO~($J = 1-0$) data, we find five molecular clouds within the velocity range between $-$25 and 8~km~s$^{-1}$ that are all characterized by conspicuous filamentary structures. We have identified eight filaments with a length of 6.38--28.45~pc, a mean H$_2$ column density of 0.70$\times$10$^{21}$--6.53$\times$10$^{21}$~cm$^{-2}$, and a line mass of 20.24--161.91~$M_\sun$ pc$^{-1}$, assuming a distance of $\sim$1.7~kpc. Gaussian fittings to the inner parts of the radial density profiles lead to a mean FWHM width of 1.13$\pm$0.01~pc. The velocity structures of most filaments present continuous distributions with slight velocity gradients. We find that turbulence is the dominant internal pressure to support the fragmentation of filaments instead of thermal pressure. Most filaments have virial parameters smaller than 2; thus, they are gravitationally bound. Four filaments have an LTE line mass close to the virial line mass. We further extract dense clumps using the $^{13}$CO data and find that 64$\%$ of the clumps are associated with the filaments. According to the complementary IR data, most filaments have associated Class~II young stellar objects. Class~I objects are mainly found to be located in the filaments with a virial parameter close to 1. Within two virialized filaments, $^{12}$CO outflows have been detected, indicating ongoing star-forming activity therein.

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

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

Flat-spectrum radio quasar PKS~1229$-$02 with a knotty and asymmetric radio morphology was identified as the optical and radio counterpart of a $\gamma$-ray source. In this paper, we study the properties, e.g. morphology, opacity, polarization and kinematics of the jet in PKS~1229$-$02 using radio interferometry. With our results, we find that the knotty and asymmetric morphology of this source may probably shaped by the interaction between its anterograde jet and the nonuniform dense ambient medium. By reproducing a Spectral Energy Distribution of PKS~1229$-$02 with the obtained kinematic parameters, we find that the relativistic beaming effect in PKS~1229$-$02 is not strong enough to produce the reported $\gamma$-ray emission, i.e. PKS~1229$-$02 may not be a $\gamma$-ray AGN. The misidentification may probably due to the poor spatial resolution of the $\gamma$-ray detector of the previous generation.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17.

Among more than 200 quasars known at $z\gtrsim6$, only one object, J0100+2802 (z=6.327), was found hosting a $>10^{10}M_\odot$ super-massive black hole (SMBH). In order to investigate the host galaxy properties of J0100+2802, we performed multi-band ALMA observations, aiming at mapping the dust continuum, [CII] and CO(6-5) emission lines with sub-kiloparsec scale resolution, as well as detecting high-J CO lines in CO(11-10), CO(10-9), and CO(7-6). The galaxy size is measured to be $R_{\rm major}=3.6\pm0.2$ kpc from the high resolution continuum observations. No ordered motion on kilo-parsec scales was found in both [CII] and CO(6-5) emissions. The velocity dispersion is measured to be 161$\pm$7 km s$^{-1}$, $\sim$3 times smaller than that estimated from the local M-$\sigma$ relation. In addition, we found that the CO emission is more concentrate (a factor of 1.8$\pm$0.4) than the [CII] emission. Together with CO(2-1) detected by VLA, we measured the CO Spectral Line Energy Distribution (SLED), which is best fit by a two-components model, including a cool component at $\sim24$ K with a density of $n_{\rm (H_2)}$=10$^{4.5}$ cm$^{-3}$, and a warm component at $\sim224$ K with a density of $n_{\rm (H_2)}$=10$^{3.6}$ cm$^{-3}$. We also fitted the dust continuum with a graybody model, which indicates that it has either a high dust emissivity $\beta\gtrsim2$ or a hot dust temperature $T_{\rm dust}\gtrsim60$ K, or a combination of both factors. The highly excited CO emission and hot dust temperature suggest that the powerful AGN in J0100+2802 could contribute to the gas and dust heating although future observations are needed to confirm this.

We perform a semi-automated search for strong gravitational lensing systems in the 9,000 deg$^2$ Dark Energy Camera Legacy Survey (DECaLS), part of the DESI Legacy Imaging Surveys (Dey et al.). The combination of the depth and breadth of these surveys are unparalleled at this time, making them particularly suitable for discovering new strong gravitational lensing systems. We adopt the deep residual neural network architecture (He et al.) developed by Lanusse et al. for the purpose of finding strong lenses in photometric surveys. We compile a training set that consists of known lensing systems in the Legacy Surveys and DES as well as non-lenses in the footprint of DECaLS. In this paper we show the results of applying our trained neural network to the cutout images centered on galaxies typed as ellipticals (Lang et al.) in DECaLS. The images that receive the highest scores (probabilities) are visually inspected and ranked. Here we present 335 candidate strong lensing systems, identified for the first time.

We have conducted a large-field simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ emission toward the Cassiopeia A (Cas A) supernova remnant (SNR), which covers a sky area of $3.5^{\circ}\times3.1^{\circ}$. The Cas giant molecular cloud (GMC) mainly consists of three individual clouds with masses on the order of $10^4-10^5\ M_{\odot}$. The total mass derived from the $\rm{^{13}CO}$ emission of the GMC is 2.1$\times10^{5}\ M_{\odot}$ and is 9.5$\times10^5\ M_{\odot}$ from the $\rm{^{12}CO}$ emission. Two regions with broadened (6$-$7 km s$^{-1}$) or asymmetric $^{12}$CO line profiles are found in the vicinity (within a 10$'\times10'$ region) of the Cas A SNR, indicating possible interactions between the SNR and the GMC. Using the GAUSSCLUMPS algorithm, 547 $^{13}$CO clumps are identified in the GMC, 54$\%$ of which are supercritical (i.e. $\alpha_{\rm{vir}}<2$). The mass spectrum of the molecular clumps follows a power-law distribution with an exponent of $-2.20$. The pixel-by-pixel column density of the GMC can be fitted with a log-normal probability distribution function (N-PDF). The median column density of molecular hydrogen in the GMC is $1.6\times10^{21}$ cm$^{-2}$ and half the mass of the GMC is contained in regions with H$_2$ column density lower than $3\times10^{21}$ cm$^{-2}$, which is well below the threshold of star formation. The distribution of the YSO candidates in the region shows no agglomeration.

We present a catalog of 290 `winged' or X-shaped radio galaxies (XRGs) extracted from the latest (17 Dec. 2014) data release of the `Very Large Array Faint Images of the Radio Sky at Twenty centimeter' (VLA FIRST survey). We have combined these radio images with their counterparts in the TIFR GMRT sky survey (TGSS) at 150 MHz (TGSS alternative data release, TGSS_ADR1), in an attempt to identify any low surface-brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 `strong' XRG candidates and 184 `probable' XRG candidates whose XRG designation must be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs presently known. Twenty-five of our 290 XRG candidates (9 `strong' and 16 `probable') are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (2 `strong' and 1 `probable'). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index $\alpha$ flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 GHz and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

Aims. We present a sample of 4388 AGNs with available radio core-dominance parameters defined as the ratio of the core flux densities to extended ones, namely, R = Score/Sext., which includes 630 Fermi-detected AGNs respect to the catalog of 4FGL, the fourth Fermi Large Area Telescope source catalog, and the rest of them are non-Fermi-detected AGNs. In our sample, 584 blazars are Fermi-detected and 1310 are not, and also consists of other subclasses such as Seyfert, Fanaroff-Riley I/II and normal galaxies. We investigate various different properties between Fermi-detected AGNs and non-Fermi-detected ones by using the core-dominance parameters as the previous study has shown that R is a good indication of beaming effect. Methods. We then calculate the radio spectral indices for whole sample and adopt {\gamma}-ray photon indices for Fermi AGNs from 4FGL catalog for discussing their different performances on different subclasses, and obtain the relation between core-dominance parameters and radio spectral indices for both Fermi and non-Fermi sources according to the two components model on radio band, which are consistent with our previous study. Results. We found that the core-dominance parameters and radio spectral indices are quite different for different subclasses of AGNs, not only for Fermi sources but also non-Fermi sources, particularly, R for the former ones is averagely higher than later ones. We also adopt the same relation on core-dominance parameters and {\gamma}-ray photon indices for Fermi sources by taking the same assumption with two components model on {\gamma}-ray band, and obtain the fitting results indicating that the {\gamma}-ray emissions of Fermi blazars are mainly from the core component, which is perhaps associated with the beaming effect. Therefore, Fermi blazars are beamed.

We present the results of the Extremely Luminous Quasar Survey in the $3\pi$ survey of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; PS1). This effort applies the successful quasar selection strategy of the Extremely Luminous Survey in the Sloan Digital Sky Survey footprint ($\sim12,000\,\rm{deg}^2$) to a much larger area ($\sim\rm{21486}\,\rm{deg}^2$). This spectroscopic survey targets the most luminous quasars ($M_{1450}\le-26.5$; $m_{i}\le18.5$) at intermediate redshifts ($z\ge2.8$). Candidates are selected based on a near-infrared JKW2 color cut using WISE AllWISE and 2MASS photometry to mainly reject stellar contaminants. Photometric redshifts ($z_{\rm{reg}}$) and star-quasar classifications for each candidate are calculated from near-infrared and optical photometry using the supervised machine learning technique random forests. We select 806 quasar candidates at $z_{\rm{reg}}\ge2.8$ from a parent sample of 74318 sources. After exclusion of known sources and rejection of candidates with unreliable photometry, we have taken optical identification spectra for 290 of our 334 good PS-ELQS candidates. We report the discovery of 190 new $z\ge2.8$ quasars and an additional 28 quasars at lower redshifts. A total of 44 good PS-ELQS candidates remain unobserved. Including all known quasars at $z\ge2.8$, our quasar selection method has a selection efficiency of at least $77\%$. At lower declinations $-30\le\rm{Decl.}\le0$ we approximately triple the known population of extremely luminous quasars. We provide the PS-ELQS quasar catalog with a total of 592 luminous quasars ($m_{i}\le18.5$, $z\ge2.8$). This unique sample will not only be able to provide constraints on the volume density and quasar clustering of extremely luminous quasars, but also offers valuable targets for studies of the intergalactic medium.

Using the Purple Mountain Observatory Delingha 13.7 m telescope, we conducted a large-scale $^{12}$CO $\left(J=1 \rightarrow 0\right)$ outflow survey (over $\sim$ 110 deg$^2$) toward the W3/4/5 complex and its surroundings. In total, 459 outflow candidates were identified. Approximately 62\% (284) were located in the Perseus arm, including W3/4/5 complex and its surroundings, while $\sim$ 35\% (162) were located in the Local arm, $\sim$ 1\% (5) in the Outer arm, and $\sim$ 2\% (8) in two interarm regions. This result indicated that star formation was concentrated in the Galactic spiral arms. The detailed spatial distribution of the outflow candidates showed that the Perseus arm presented the most active star formation among the study regions. The W3/4/5 complex is a great region to research massive star formation in a triggered environment. A key region, which has been well-studied by other researches, is in the eastern high-density W3 complex that neighbors the W4 complex. Two shell-like structures in the Local arm contain candidates that can be used to study the impact on star formation imposed by massive or intermediate-mass stars in relatively isolated systems. The majority of outflow candidates in the two interarm regions and the Outer arm are located in filamentary structures.

We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

Recent studies have revealed that solar jets involving minifilament eruptions may be initiated under the well-known magnetic-breakout mechanism. Before or just at the onset of those jets, there should be a current sheet, where breakout magnetic reconnection takes place, between open fields and the outside of the jet-base arcade carrying minifilament in its core. In this paper we present a jet produced by eruption of two minifilaments lying at the jet base. A current sheet is directly detected near the jet base before the onset of the eruption, suggesting the magnetic-breakout mechanism. However, we further find that the current sheet undergoes a transition. The current sheet first shortens to zero in length, but then lengthens towards an orthogonal direction relative to its initial orientation. The change of the current sheet gives rise to a reversal of the breakout reconnection, as the inflow and outflow regions before the transition become the outflow and inflow regions after the transition, respectively. We therefore propose that this observation provides evidence for the so-called oscillatory reconnection which is defined by a series of reconnection reversals but not yet proved to exist in real plasma environment of the solar atmosphere.

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.

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High spatial resolution measurements of the source size and displacement can discriminate between these scenarios. We present Very Long Baseline Interferometry observations, performed 207.4 days after the merger, using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milliarcseconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

Recent studies have shown that ocean dynamics can have a significant warming effect on the permanent night sides of 1 to 1 tidally locked terrestrial exoplanets with Earth-like atmospheres and oceans in the middle of the habitable zone. However, the impact of ocean dynamics on the habitable zone's boundaries (inner edge and outer edge) is still unknown and represents a major gap in our understanding of this type of planets. Here we use a coupled atmosphere-ocean global climate model to show that planetary heat transport from the day to night side is dominated by the ocean at lower stellar fluxes and by the atmosphere near the inner edge of the habitable zone. This decrease in oceanic heat transport (OHT) at high stellar fluxes is mainly due to weakening of surface wind stress and a decrease in surface shortwave energy deposition. We further show that ocean dynamics have almost no effect on the observational thermal phase curves of planets near the inner edge of the habitable zone. For planets in the habitable zone's middle range, ocean dynamics moves the hottest spot on the surface eastward from the substellar point. These results suggest that future studies of the inner edge may devote computational resources to atmosphere-only processes such as clouds and radiation. For studies of the middle range and outer edge of the habitable zone, however, fully coupled ocean-atmosphere modeling will be necessary. Note that due to computational resource limitations, only one rotation period (60 Earth days) has been systematically examined in this study; future work varying rotation period as well as other parameters such as atmospheric mass and composition is required.

We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.

The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed.

The $\mu$eV axion is a well-motivated extension to the standard model. The Axion Dark Matter eXperiment (ADMX) collaboration seeks to discover this particle by looking for the resonant conversion of dark-matter axions to microwave photons in a strong magnetic field. In this Letter, we report results from a pathfinder experiment, the ADMX "Sidecar," which is designed to pave the way for future, higher mass, searches. This testbed experiment lives inside of and operates in tandem with the main ADMX experiment. The Sidecar experiment excludes masses in three widely spaced frequency ranges (4202-4249, 5086-5799, and 7173-7203 MHz). In addition, Sidecar demonstrates the successful use of a piezoelectric actuator for cavity tuning. Finally, this publication is the first to report data measured using both the TM$_{010}$ and TM$_{020}$ modes.

The Milky Way Imaging Scroll Painting (MWISP) project is an unbiased Galactic plane CO survey for mapping regions of l=-10deg to +250deg and |b|<5.2deg with the 13.7 m telescope of the Purple Mountain Observatory. The legacy survey aims to observe the 12CO, 13CO, and C18O (J=1-0) lines simultaneously with full sampling using the nine-beam Superconducting SpectroScopic Array Receiver (SSAR) system with an instantaneous bandwidth of 1 GHz. In this paper, the completed 250 deg^2 data from l=+25.8deg to +49.7deg are presented with a grid spacing of 30" and a typical rms noise level of ~0.5 K for 12CO at the channel width of 0.16 km/s and ~0.3 K for 13CO and C18O at 0.17 km/s. The high quality data with moderate resolution (50"), uniform sensitivity, and high spatial dynamic range, allow us to investigate the details of molecular clouds (MCs) traced by the three CO isotope lines. Three interesting examples are briefly investigated, including distant Galactic spiral arms traced by CO emission with VLSR<0 km/s, the bubble-like dense gas structure near the HII region W40, and the MCs distribution perpendicular to the Galactic plane.

We have designed the Extremely Luminous Quasar Survey (ELQS) to provide a highly complete census of unobscured UV-bright quasars during the cosmic noon, $z=2.8-5.0$. Here we report the discovery of 70 new quasars in the ELQS South Galactic Cap (ELQS-S) quasar sample, doubling the number of known extremely luminous quasars in $4,237.3\,\rm{deg}^2$ of the SDSS footprint. These observations conclude the ELQS and we present the properties of the full ELQS quasar catalog, containing 407 quasars over $11,838.5\,\rm{deg}^2$. Our novel ELQS quasar selection strategy resulted in unprecedented completeness at the bright end and allowed us to discover 109 new quasars in total. This marks an increase of $\sim36\%$ (109/298) to the known population at these redshifts and magnitudes, while we further are able to retain a selection efficiency of $\sim80\%$. On the basis of 166 quasars from the full ELQS quasar catalog, who adhere to the uniform criteria of the 2MASS point source catalog, we measure the bright-end quasar luminosity function (QLF) and extend it one magnitude brighter than previous studies. Assuming a single power law with exponential density evolution for the functional form of the QLF, we retrieve the best fit parameters from a maximum likelihood analysis. We find a steep bright-end slope of $\beta\approx-4.1$ and we can constrain the bright-end slope to $\beta\leq-3.4$ with $99\%$ confidence. The density is well modeled by the exponential redshift evolution, resulting in a moderate decrease with redshift ($\gamma\approx-0.4$).

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

This is the fourth paper in a series of publications aiming at discovering quasars at the epoch of reionization. In this paper, we expand our search for $z\sim 7$ quasars to the footprint of the Dark Energy Survey (DES) Data Release One (DR1), covering $\sim 5000$ deg$^2$ of new area. We select $z\sim 7$ quasar candidates using deep optical, near-infrared (near-IR) and mid-IR photometric data from the DES DR1, the VISTA Hemisphere Survey (VHS), the VISTA Kilo-degree Infrared Galaxy (VIKING) survey, the UKIRT InfraRed Deep Sky Surveys -- Large Area Survey (ULAS) and the unblurred coadds from the Wide-field Infrared Survey Explore ($WISE$) images (unWISE). The inclusion of DES and unWISE photometry allows the search to reach $\sim$ 1 magnitude fainter, comparing to our $z \gtrsim 6.5$ quasar survey in the northern sky (Wang et al. 2018). We report the initial discovery and spectroscopic confirmation of six new luminous quasars at $z>6.4$, including an object at $z=7.02$, the fourth quasar yet known at $z>7$, from a small fraction of candidates observed thus far. Based on the recent measurement of $z \sim 6.7 $ quasar luminosity function using the quasar sample from our survey in the northern sky, we estimate that there will be $\gtrsim$ 55 quasars at $z > 6.5$ at $M_{1450} < -24.5$ in the full DES footprint.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

Active galactic nuclei (AGNs) can be divided into two major classes, namely radio loud and radio quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff & Riley type I/II (FRI/II) radio galaxies. Following our previous work (Fan et al. 2011), we present a sample of 2400 sources with measured radio flux densities of the core and extended components. The sample contains 250 BL Lacs, 520 quasars, 175 Seyferts, 1178 galaxies, 153 Fanaroff & Riley type I and type II galaxies and 104 unidentified sources. We then calculate the radio core-dominance parameters and spectral indices and study their relationship. Our analysis shows that the core-dominance parameters and spectral indices are quite different for different types of sources. We also confirm that the correlation between core-dominance parameter and spectral index exists for a large sample presented in this work.

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.

Employing Solar Dynamic Observatory/Atmosphertic Imaging Assembly (AIA) multi-wavelength images, we study an eruption of two crossing filaments, and firstly report the current sheets (CSs) connecting the lower flare ribbons and the upper erupting filaments. On July 8, 2014, two crossing filaments are observed in the NOAA active region (AR) 12113. The lower-lying filament rises first, and then meets the higher-lying one. Thereafter, both of them erupt together. The filament eruption draws the overlying magnetic field lines upward, leading to the approach of two legs, with opposite magnetic polarities, of the overlying field lines. Two sets of bright CSs form at the interface of these two legs, and magnetic reconnection takes place in the CSs producing the underneath flare ribbons and post-flare loops. Several bright plasmoids appear in the CSs, and propagate along the CSs bi-directionally. The CSs and plasmoids are observed in AIA multi-wavelength channels, indicating that both of them have been heated during the reconnection process, with hot and warm components. Employing the differential emission measure (EM) analysis, we find that both the temperature and EM of the CSs decrease from the flare arcades outward to the erupting filaments, and those of the plasmoids are significantly larger than the regions where no plasmoid is detected.

In this paper, we study the onset process of a solar eruption on 21 February 2015, focusing on its unambiguous precursor phase. With multi-wavelength imaging observations from the Atmospheric Imaging Assembly (AIA), definitive tether-cutting (TC) reconnection signatures, i.e., flux convergence and cancellation, bidirectional jets, as well as topology change of hot loops, were clearly observed below the pre-eruption filament. As TC reconnection progressed between the sheared arcades that enveloped the filament, a channel-like magnetic flux rope (MFR) arose in multi-wavelength AIA passbands wrapping around the main axis of the filament. With the subsequent ascent of the newborn MFR, the filament surprisingly split into three branches. After a 7-hour slow rise phase, the high-lying branch containing by the MFR abruptly accelerated causing a two-ribbon flare; while the two low-lying branches remained stable forming a partial eruption. Complemented by kinematic analysis and decay index calculation, we conclude that TC reconnection played a key role in building up the eruptive MFR and triggering its slow rise. The onset of the torus instability may have led the high-lying branch into the standard eruption scenario in the fashion of a catastrophe.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

We present the final results from our survey of luminous $z \sim $ 5.5 quasars. This is the first systematic quasar survey focusing on quasars at $z \sim$ 5.5, during the post-reionization epoch. It has been challenging to select quasars at $5.3 < z < 5.7$ using conventional color selections, due to their similar optical colors to those of late-type stars, especially M dwarfs. We developed a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR, and mid-IR photometry, using data from the Sloan Digital Sky Survey (SDSS), PanSTARR1 (PS1), the UKIRT Infrared Deep Sky Surveys - Large Area Survey, the UKIRT Hemisphere Survey, the VISTA Hemisphere Survey, and the Wide Field Infrared Survey Explorer ($WISE$), covering $\sim$ 11000 deg$^2$ of high galactic latitude sky. In this paper, we present the discovery of 15 new quasars at $z\sim 5.5$. Together with results from Yang et al. (2017), our survey provides a complete, flux-limited sample of 31 quasars at $5.3 \le z \le 5.7$. We measure the quasar spatial density at $z \sim 5.5$ and $M_{1450} < -26.2$. Our result is consistent with the rapid decline of quasar spatial density from $z = 5$ to 6, with $k=-0.66 \pm 0.05$ ($\rho(z) \propto 10^{kz}$). In addition, we present a new survey using optical colors only from the full PS1 area for luminous quasars at $z = 5.0 - 5.5$, beyond the SDSS footprint, and report the preliminary results from this survey, including 51 new quasars discovered at $4.61\le z \le5.71$.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

We present nearly simultaneous NuSTAR and XMM-Newton observations of the nearby (832 kpc) ultraluminous X-ray source (ULX) M33 X-8. M33 X-8 has a 0.3-10 keV luminosity of LX ~ 1.4 x 10^39 erg/s, near the boundary of the "ultraluminous" classification, making it an important source for understanding the link between typical Galactic X-ray binaries and ULXs. Past studies have shown that the 0.3-10 keV spectrum of X-8 can be characterized using an advection-dominated accretion disk model. We find that when fitting to our NuSTAR and XMM-Newton observations, an additional high-energy (>10 keV) Comptonization component is required, which allows us to rule out single advection-dominated disk and classical sub-Eddington models. With our new constraints, we analyze XMM-Newton data taken over the last 17 years to show that small (~30%) variations in the 0.3-10 keV flux of M33 X-8 result in spectral changes similar to those observed for other ULXs. The two most likely phenomenological scenarios suggested by the data are degenerate in terms of constraining the nature of the accreting compact object (i.e., black hole versus neutron star). We further present a search for pulsations using our suite of data; however, no clear pulsations are detected. Future observations designed to observe M33 X-8 at different flux levels across the full 0.3-30 keV range would significantly improve our constraints on the nature of this important source.

The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9{\deg} for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1{\sigma} and with the Auger result to within 1.4{\sigma}, meaning that it is not capable to discriminate between these two hypotheses.

Active galactic nuclei (AGN) accreting at rates close to the Eddington limit can host radiatively driven mildly relativistic outflows. Some of these X-ray absorbing but powerful outflows may produce strong shocks resulting in a significant non-thermal emission. This outflow-driven radio emission may be detectable in the radio-quiet quasar PDS 456 since it has a bolometric luminosity reaching the Eddington limit and a relativistic wide-aperture X-ray outflow with a kinetic power high enough to quench the star formation in its host galaxy. To investigate this possibility, we performed very-long-baseline interferometric (VLBI) observations of the quasar with the European VLBI Network (EVN) at 5 GHz. The EVN image with the full resolution reveals two faint and diffuse radio components with a projected separation of about 20 pc and an average brightness temperature of around two million Kelvin. In relation to the optical sub-mas-accuracy position measured by the Gaia mission, the two components are very likely on opposite sides of an undetected radio core. The VLBI structure at the deca-pc scale can thus be either a young jet or a bidirectional radio-emitting outflow, launched in the vicinity of a strongly accreting central engine. Two diffuse components at the hecto-pc scale, likely the relic radio emission from the past AGN activity, are tentatively detected on each side in the low-resolution EVN image.

BL Lac objects are known to have compact jets inclined to our line of sight at a small angle, showing prominent radio emission. Two radio-weak BL Lac candidates with no counterparts in current radio surveys were recently reported by Massaro et al. (2017). Both sources were selected as candidate low-energy counterparts of unassociated Fermi {\gamma}-ray sources. We carried out very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 5 GHz to explore their radio properties at milli-arcsecond (mas) scale. One target, J1410+7405, is clearly detected with the EVN. Its measured 5-GHz flux density, 2.4 mJy, is consistent with recent interferometric measurements with the Karl G. Jansky Very Large Array, suggesting that the radio emission is confined to the inner < ~10 mas region. J1410+7405 is therefore identified as a radio-loud jetted active galactic nucleus, and its brightness temperature exceeds ~10^9 K. Its properties are similar to those of other {\gamma}-ray detected BL Lac objects. On the other hand, the second target, J0644+6031, remains undetected with the EVN with a 6-{\sigma} brightness upper limit of 0.12 mJy/beam. This source is thus radio-quiet, confirming its peculiarity, or possibly questioning its BL Lac nature.

We present the results of CS J=2-1 mapping observations towards 39 massive star-forming regions selected from the previous CO line survey of cold IRAS sources with high-velocity CO flows along the Galactic plane (Yang et al. 2002). All sources are detected in CS J=2-1 showing the existence of CS clumps around the IRAS sources. However, one-third of the sources are not deeply embedded in the dense clumps by comparison of the central powering IRAS sources and the morphologies of CS clumps. Physical parameters of the dense molecular clumps are presented. We have identified 12 high-mass protostellar object (HMPO) candidates by checking the association between the dense cores and the IRAS sources, the detection of water maser, and the radio properties towards the IRAS sources. We find that the HMPO sources are characterized by low FIR luminosity to virial mass ratios since they are in very early evolutionary stages when the massive protostars have not reached their full luminosities, which are typical for zero-age main sequence stars. Large turbulent motion in the HMPO sources may be largely due to the large kinetic energy ejected by the central protostars formed in the dense clumps. However, alternative means or undetected outflows may also be responsible for the turbulence in the clumps.

An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional (1D) radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G- and M-stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 um) and in the region between 0.2 and 1.5 um. Differences in outgoing longwave radiation increase with surface temperature and reach 10-20 Wm^-2; differences in shortwave reach up to 60 Wm^-2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ~10% of modern Earth's solar constant (i.e., ~34 Wm^-2 in global mean) among band models and ~3% between the two line-by-line models. These comparisons show that future work is needed focusing on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

Motivated by the important role of the ocean in the Earth climate system, here we investigate possible scenarios of ocean circulations on exoplanets using a one-layer shallow water ocean model. Specifically, we investigate how planetary rotation rate, wind stress, fluid eddy viscosity and land structure (a closed basin vs. a reentrant channel) influence the pattern and strength of wind-driven ocean circulations. The meridional variation of the Coriolis force, arising from planetary rotation and the spheric shape of the planets, induces the western intensification of ocean circulations. Our simulations confirm that in a closed basin, changes of other factors contribute to only enhancing or weakening the ocean circulations (e.g., as wind stress decreases or fluid eddy viscosity increases, the ocean circulations weaken, and vice versa). In a reentrant channel, just as the Southern Ocean region on the Earth, the ocean pattern is characterized by zonal flows. In the quasi-linear case, the sensitivity of ocean circulations characteristics to these parameters is also interpreted using simple analytical models. This study is the preliminary step for exploring the possible ocean circulations on exoplanets, future work with multi-layer ocean models and fully coupled ocean-atmosphere models are required for studying exoplanetary climates.

Ongoing and future space missions aim to identify potentially habitable planets in our Solar System and beyond. Planetary habitability is determined not only by a planet's current stellar insolation and atmospheric properties, but also by the evolutionary history of its climate. It has been suggested that icy planets and moons become habitable after their initial ice shield melts as their host stars brighten. Here we show from global climate model simulations that a habitable state is not achieved in the climatic evolution of those icy planets and moons that possess an inactive carbonate-silicate cycle and low concentrations of greenhouse gases. Examples for such planetary bodies are the icy moons Europa and Enceladus, and certain icy exoplanets orbiting G and F stars. We find that the stellar fluxes that are required to overcome a planet's initial snowball state are so large that they lead to significant water loss and preclude a habitable planet. Specifically, they exceed the moist greenhouse limit, at which water vapour accumulates at high altitudes where it can readily escape, or the runaway greenhouse limit, at which the strength of the greenhouse increases until the oceans boil away. We suggest that some icy planetary bodies may transition directly to a moist or runaway greenhouse without passing through a habitable Earth-like state.

Using the PMO-13.7m millimeter telescope at Delingha in China, we have conducted a large-scale simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=1-0 emission toward the Rosette molecular cloud (RMC) region with a sky coverage of 3.5 $\times$ 2.5 deg$^2$. The majority of the emission in the region comes from the RMC complex with velocities lying in the range from -2 km s$^{-1}$ to 20.5 km s$^{-1}$. Beyond this velocity range, 73 molecular clumps are identified with kinematic distances from 2.4 kpc to 11 kpc. Based on the spatial and velocity distribution, nine individual clouds, C1-C9, have been identified for the RMC complex. It appears that the C3 cloud is different from other clouds in the RMC complex in view of its characteristic velocity, excitation temperature, and velocity dispersion. Most of the young stellar clusters in the region are located in positions of both high column density and high excitation temperature. Seven new molecular filaments are discovered in the RMC complex. Evidence for cloud-cloud collision is found in the region of young stellar clusters REFL9 and PouF, showing that these young stellar clusters probably result from a cloud-cloud collision. The abundance ratios of $^{13}$CO to C$^{18}$O in the region have a mean value of 13.7 which is 2.5 times larger than the solar system value, showing that UV photons from the nearby OB clusters have strong influence on the chemistry of clouds in the RMC complex.

We present our deep $Y$-band imaging data of a two square degree field within the F22 region of the VIMOS VLT Deep Survey. The observations were conducted using the WIRCam instrument mounted at the Canada--France--Hawaii Telescope (CFHT). The total on-sky time was 9 hours, distributed uniformly over 18 tiles. The scientific goals of the project are to select faint quasar candidates at redshift $z>2.2$, and constrain the photometric redshifts for quasars and galaxies. In this paper, we present the observation and the image reduction, as well as the photometric redshifts that we derived by combining our $Y$-band data with the CFHTLenS $u^*g'r'i'z'$ optical data and UKIDSS DXS $JHK$ near-infrared data. With $J$-band image as reference total $\sim$80,000 galaxies are detected in the final mosaic down to $Y$-band $5\sigma$ point source limiting depth of 22.86 mag. Compared with the $\sim$3500 spectroscopic redshifts, our photometric redshifts for galaxies with $z<1.5$ and $i'\lesssim24.0$ mag have a small systematic offset of $|\Delta{z}|\lesssim0.2$, 1$\sigma$ scatter $0.03<\sigma_{\Delta z} < 0.06$, and less than 4.0% of catastrophic failures. We also compare to the CFHTLenS photometric redshifts, and find that ours are more reliable at $z\gtrsim0.6$ because of the inclusion of the near-infrared bands. In particular, including the $Y$-band data can improve the accuracy at $z\sim 1.0-2.0$ because the location of the 4000\AA-break is better constrained. The $Y$-band images, the multi-band photometry catalog and the photometric redshifts are released at \url{http://astro.pku.edu.cn/astro/data/DYI.html}.

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.

Recently, a rotating black hole solution was found in Vector-Tensor Galileons modified gravity that has some significant differences from the Kerr black hole solution. We study the iron line shape that is part of the reflection spectrum of accretion disks around black holes in this new black hole solution. We simulate and compare the iron lines of this solution with those of the Kerr solution to see if the technique of iron line spectroscopy can be used as a tool to test Vector-Tensor Galileons modified gravity. Our analysis shows that current X-ray facilities can, in principle, be used to test and place constraints on Vector-Tensor Galileons modified gravity in the strong coupling and ultraspinning regimes.

Using the 135-second cadence of the photospheric vector data provided by the Helioseismic and Magnetic Imager telescope on board the Solar Dynamic Observatory, we examined the time-evolution of magnetic helicity fluxes across the photosphere during 16 flares with the energy class lower than M5.0. During the flare in 4 out of 16 events, we found impulsive changes in the helicity fluxes. This indicates that even the flare with less energy could be associated with anomalistic transportation of the magnetic helicity across the photosphere. Accompanying the impulsive helicity fluxes, the poynting fluxes across the photosphere evolved from positive to negative. As such, the transportations of magnetic energy across the photosphere were toward solar interior during these flares. In each of the 4 events, the impulsive change in the helicity flux was always mainly contributed by abrupt change in horizontal velocity field on a sunspot located near the flaring polarity inversion line. The velocity field on each sunspot shows either an obvious vortex patten or an shearing patten relative to the another magnetic polarity, which tended to relax the magnetic twist or shear in the corona. During these flares, abrupt change in the Lorentz force acting on these sunspots were found. The rotational motions and shearing motions of these sunspots always had the same directions with the resultant Lorentz forces. These results support the view that the impulsive helicity transportation during the flare could be driven by the change in the Lorentz force applied on the photosphere.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

We present the North Galactic Cap sample of the Extremely Luminous Quasar Survey (ELQS-N), which targets quasars with $M_{1450}<-27$ at $2.8 \leq z < 5$ in an area of $\sim7600\,\rm{deg}^2$ of the Sloan Digital Sky Survey (SDSS) footprint with $90\text{\textdegree}<\rm{RA}<270\text{\textdegree}$. Based on a near-infrared/infrared \textit{JKW2} color cut, the ELQS selection efficiently uses random forest methods to classify quasars and to estimate photometric redshifts; this scheme overcomes some of the difficulties of pure optical quasar selection at $z\approx3$. As a result, we retain a completeness of $>70\%$ over $z\sim3.0-5.0$ at $m_{i}\lesssim17.5$, limited toward fainter magnitudes by the depth of the Two Micron All Sky Survey (2MASS). The presented quasar catalog consists of a total of 270 objects, of which 39 are newly identified in this work with spectroscopy obtained at the Vatican Advanced Technology Telescope and the MMT $6.5\,\rm{m}$ telescope. In addition to the high completeness, which allowed us to discover new quasars in the already well-surveyed SDSS North Galactic Cap, the efficiency of our selection is relatively high at $\sim79\%$. Using 120 objects of this quasar sample we are able to extend the previously measured optical quasar luminosity function (QLF) by one magnitude toward the bright end at $2.8 \leq z \leq 4.5$. A first analysis of the QLF suggests a relatively steep bright-end slope of $\beta\approx-4$ for this sample. This result contrasts with previous results in the same redshift range, which find a much flatter slope around $\beta\sim-2.5$, but agrees with recent measurements of the bright-end slope at lower and higher redshifts. Our results constrain the bright-end slope at $z=2.8-4.5$ to $\beta<-2.94$ with a 99\% confidence.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

With new high-resolution CO and HI data, we revisited the large-scale interstellar medium (ISM) environment toward the SS 433/W50 system. We find that two interesting molecular cloud (MC) concentrations, G39.315-1.155 and G40.331-4.302, are well aligned along the precession cone of SS 433 within a smaller opening angle of ~7deg. The kinematic features of the two MCs at ~73--84 km/s, as well as those of the corresponding atomic-gas counterparts, are consistent with the kinematic characteristics of SS 433. That is, the receding gas from SS 433 jet is probably responsible for the redshifted feature of G39.315-1.155 near the Galactic plane and the approaching one may power the blueshifted gas of G40.331-4.302 toward the observer. Moreover, the HI emission at VLSR~70--90 km/s displays the morphological resemblance with the radio nebula W50. We suggest that the VLSR=77\pm5 km/s gas is physically associated with SS 433/W50, leading to a near kinematic distance of 4.9\pm0.4 kpc for the system. The observed gas features, which are located outside the current radio boundaries of W50, are probably the fossil record of jet-ISM interactions at ~10^5 years ago. The energetic jets of the unique microquasar have profound effects on its ISM environment, which may facilitate the formation of molecular gas on the timescale of <0.1 Myr for the ram pressure of ~2x10^6 K cm-3.

We present the photometric calibration of the Beijing-Arizona Sky Survey (BASS) and Mayall z-band Legacy Survey (MzLS), which are two of the three wide-field optical legacy imaging surveys to provide the baseline targeting data for the Dark Energy Spectroscopic Instrument (DESI) project. The method of our photometric calibration is subdivided into the external and internal processes. The former utilizes the point-source objects of Pan-STARRS1 survey (PS1) as the reference standards to achieve the zero points of the absolute flux for individual exposures. And then the latter revise the zero points to make them consistent across the survey based on multiple tilings and large offset overlaps. Our process achieves a homogeneous photometric calibration over most of the sky with precision better than 10 mmag for g and r bands, 15 mmag for z band. The accuracy of the calibration is better than 1% at the bright end (16-18 mag) over most of the survey area.

Using high spatial and temporal data from the \emph{Solar Dynamics Observatory} (\emph{SDO}) and the \emph{Interface Region Imaging Spectrograph} (\emph{IRIS}), several observational signatures of magnetic reconnection in the course of magnetic flux cancellation are presented, including two loop-loop interaction processes, multiple plasma blob ejections, and a sheet-like structure that appeared above the flux cancellation sites with a Y-shaped and an inverted Y-shaped ends. The \emph{IRIS} 1400 \AA\ observations show that the plasma blobs were ejected from the tip of the Y-shaped ends of the sheet-like structure. Obvious photospheric magnetic flux cancellation occurred after the first loop-loop interaction and continued until the end of the observation. Complemented by the nonlinear force-free field extrapolation, we found that two sets of magnetic field lines, which reveal an X-shaped configuration, align well with the interacted coronal loops. Moreover, a magnetic null point is found to be situated at about $0.9$ Mm height right above the flux cancellation sites and located between the two sets of magnetic field lines. These results suggest that the flux cancellation might be a result of submergence of magnetic field lines following magnetic reconnection that occurs in the lower atmosphere of the Sun, and the ejected plasma blobs should be plasmoids created in the sheet-like structure due to the tearing-mode instability. This observation reveals detailed magnetic field structure and dynamic process above the flux cancellation sites and will help us to understand magnetic reconnection in the lower atmosphere of the Sun.

We present observations of the formation process of a small-scale filament on the quiet Sun during 5-6 February 2016 and investigate its formation cause. Initially, a small dipole emerged and its associated arch filament system was found to reconnect with overlying coronal fields accompanied by numerous EUV bright points. When bright points faded out, many elongated dark threads formed bridging the positive magnetic element of dipole and external negative network fields. Interestingly, an anti-clockwise photospheric rotational motion (PRM) set in within the positive endpoint region of newborn dark threads following the flux emergence and lasted for more than 10 hours. Under the drive of the PRM, these dispersive dark threads gradually aligned along the north-south direction and finally coalesced into an inverse S-shaped filament. Consistent with the dextral chirality of the filament, magnetic helicity calculations show that an amount of negative helicity was persistently injected from the rotational positive magnetic element and accumulated during the formation of the filament. These observations suggest that twisted emerging fields may lead to the formation of the filament via reconnection with pre-existing fields and release of its inner magnetic twist. The persistent PRM might trace a covert twist relaxation from below photosphere to the low corona.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

We report the detection of the radio afterglow of a long gamma-ray burst (GRB) 111005A at 5-345 GHz, including the very long baseline interferometry observations with the positional error of 0.2 mas. The afterglow position is coincident with the disk of a galaxy ESO 580-49 at z= 0.01326 (~1" from its center), which makes GRB 111005A the second closest GRB known to date, after GRB 980425. The radio afterglow of GRB 111005A was an order of magnitude less luminous than those of local low-luminosity GRBs, and obviously than those of cosmological GRBs. The radio flux was approximately constant and then experienced an unusually rapid decay a month after the GRB explosion. Similarly to only two other GRBs, we did not find the associated supernovae (SN), despite deep near- and mid-infrared observations 1-9 days after the GRB explosion, reaching ~20 times fainter than other SNe associated with GRBs. Moreover, we measured twice solar metallicity for the GRB location. The low gamma-ray and radio luminosities, rapid decay, lack of a SN, and super-solar metallicity suggest that GRB 111005A represents a different rare class of GRBs than typical core-collapse events. We modelled the spectral energy distribution of the GRB 111005A host finding that it is a dwarf, moderately star-forming galaxy, similar to the host of GRB 980425. The existence of two local GRBs in such galaxies is still consistent with the hypothesis that the GRB rate is proportional to the cosmic star formation rate (SFR) density, but suggests that the GRB rate is biased towards low SFRs. Using the far-infrared detection of ESO 580-49, we conclude that the hosts of both GRBs 111005A and 980425 exhibit lower dust content than what would be expected from their stellar masses and optical colours.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

We report the evidence for the anti-correlation between pulsed fraction (PF) and luminosity of the X-ray pulsar SXP 1323, found for the first time in a luminosity range $10^{35}$--$10^{37}$ erg s$^{-1}$ from observations spanning 15 years. The phenomenon of a decrease in X-ray PF when the source flux increases has been observed in our pipeline analysis of other X-ray pulsars in the Small Magellanic Cloud (SMC). It is expected that the luminosity under a certain value decreases as the PF decreases due to the propeller effect. Above the propeller region, an anti-correlation between the PF and flux might occur either as a result of an increase in the un-pulsed component of the total emission or a decrease of the pulsed component. Additional modes of accretion may also be possible, such as spherical accretion and a change in emission geometry. At higher mass accretion rates, the accretion disk could also extend closer to the neutron star (NS) surface, where a reduced inner radius leads to hotter inner disk emission. These modes of plasma accretion may affect the change in the beam configuration to fan-beam dominant emission.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

To date, PMN J2134-0419 (at a redshift z=4.33) is the second most distant quasar known with a milliarcsecond-scale morphology permitting direct estimates of the jet proper motion. Based on two-epoch observations, we constrained its radio jet proper motion using the very long baseline interferometry (VLBI) technique. The observations were conducted with the European VLBI Network (EVN) at 5 GHz on 1999 November 26 and 2015 October 6. We imaged the central 10-pc scale radio jet emission and modeled its brightness distribution. By identifying a jet component at both epochs separated by 15.86 yr, a proper motion of mu=0.035 +- 0.023 mas/yr is found. It corresponds to an apparent superluminal speed of beta_a=4.1 +- 2.7 c . Relativistic beaming at both epochs suggests that the jet viewing angle with respect to the line of sight is smaller than 20 deg, with a minimum bulk Lorentz factor Gamma=4.3. The small value of the proper motion is in good agreement with the expectations from the cosmological interpretation of the redshift and the current cosmological model. Additionally we analyzed archival Very Large Array observations of J2143-0419 and found indication of a bent jet extending to ~30 kpc.

We present a 45 ks Chandra observation of the quasar ULAS J1342+0928 at z=7.54. We detect 14.0^{+4.8}_{-3.7} counts from the quasar in the observed-frame energy range 0.5-7.0 keV (6-sigma detection), representing the most distant non-transient astronomical source identified in X-rays to date. The present data are sufficient only to infer rough constraints on the spectral parameters. We find an X-ray hardness ratio of HR = -0.51^{+0.26}_{-0.28} between the 0.5-2.0 keV and 2.0-7.0 keV ranges and derive a power-law photon index of Gamma = 1.95^{+0.55}_{-0.53}. Assuming a typical value for high-redshift quasars of Gamma = 1.9, ULAS J1342+0928 has a 2-10 keV rest-frame X-ray luminosity of L_{2-10} = 11.6^{+4.3}_{-3.5} x 10^{44} erg/s. Its X-ray-to-optical power-law slope is alpha_{OX}=-1.67^{+0.16}_{-0.10}, consistent with the general trend indicating that the X-ray emission in the most bolometrically powerful quasars is weaker relative to their optical emission.

This is the second installment for the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) Quasar Survey, which includes quasars observed from September 2013 to June 2015. There are 9024 confirmed quasars in DR2 and 10911 in DR3. After cross-match with the SDSS quasar catalogs and NED, 12126 quasars are discovered independently. Among them 2225 quasars were released by SDSS DR12 QSO catalogue in 2014 after we finalised the survey candidates. 1801 sources were identified by SDSS DR14 as QSOs. The remaining 8100 quasars are considered as newly founded, and among them 6887 quasars can be given reliable emission line measurements and the estimated black hole masses. Quasars found in LAMOST are mostly located at low-to-moderate redshifts, with a mean value of 1.5. The highest redshift observed in DR2 and DR3 is 5. We applied emission line measurements to H$\alpha$, H$\beta$, Mg{\sc ii} and C{\sc iv}. We deduced the monochromatic continuum luminosities using photometry data, and estimated the virial black hole masses for the newly discovered quasars. Results are compiled into a quasar catalog, which will be available online.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

Using a simplified framework, we attempt to explain the recent DAMPE cosmic $e^+ + e^-$ flux excess by leptophilic Dirac fermion dark matter (LDM). The scalar ($\Phi_0$) and vector ($\Phi_1$) mediator fields connecting LDM and Standard Model particles are discussed. Under constraints of DM relic density, gamma-rays, cosmic-rays and Cosmic Microwave Background (CMB), we find that the couplings $P \otimes S$, $P \otimes P$, $V \otimes A$ and $V \otimes V$ can produce the right bump in $e^+ + e^-$ flux for a DM mass around 1.5 TeV with a natural thermal annihilation cross-section $<\sigma v> \sim 3 \times 10^{-26} cm^3/s$ today. Among them, $V \otimes V$ coupling is tightly constrained by PandaX-II data (although LDM-nucleus scattering appears at one-loop level) and the surviving samples appear in the resonant region, $m_{\Phi_1} \simeq 2m_{\chi}$. We also study the related collider signatures, such as dilepton production $pp \to \Phi_1 \to \ell^+\ell^-$, and muon $g-2$ anomaly. Finally, we present a possible $U(1)_X$ realization for such leptophilic dark matter.

We present a large scale survey of CO outflows in the Gem OB1 molecular cloud complex and its surroundings using the Purple Mountain Observatory Delingha 13.7 m telescope. A total of 198 outflow candidates were identified over a large area ($\sim$ 58.5 square degrees), of which 193 are newly detected. Approximately 68% (134/198) are associated with the Gem OB1 molecular cloud complex, including clouds GGMC 1, GGMC 2, BFS 52, GGMC 3 and GGMC 4. Other regions studied are: Local Arm (Local Lynds, West Front), Swallow, Horn, and Remote cloud. Outflow candidates in GGMC 1, BFS 52, and Swallow are mainly located at ring-like or filamentary structures. To avoid excessive uncertainty in distant regions ($\gtrsim$ 3.8 kpc), we only estimated the physical parameters for clouds in the Gem OB1 molecular cloud complex and in the Local arm. In those clouds, the total kinetic energy and the energy injection rate of the identified outflow candidates are $\lesssim$ 1% and $\lesssim$ 3% of the turbulent energy and the turbulent dissipation rate of each cloud, indicating that the identified outflow candidates cannot provide enough energy to balance turbulence of their host cloud at the scale of the entire cloud (several pc to dozens of pc). The gravitational binding energy of each cloud is $\gtrsim$ 135 times the total kinetic energy of the identified outflow candidates within the corresponding cloud, indicating that the identified outflow candidates cannot cause major disruptions to the integrity of their host cloud at the scale of the entire cloud.

We propose to interpret the DAMPE electron excess at 1.5 TeV through scalar or Dirac fermion dark matter (DM) annihilation with doubly charged scalar mediators that have lepton-specific Yukawa couplings. Hierarchy of such lepton-specific Yukawa couplings is generated through the Froggatt-Nielsen mechanism, so that the dark matter annihilation products can be dominantly electrons. Stringent constraints from LEP2 on intermediate vector boson production can be evaded in our scenarios. In the case of scalar DM, we discuss one scenario with DM annihilating directly to leptons and the other scenario with DM annihilating to scalar mediators followed by their decays. We also discuss the Breit-Wigner resonant enhancement and the Sommerfeld enhancement in case that the s-wave annihilation process is small or helicity suppressed. With both types of enhancement, constraints on the parameters can be relaxed and new ways for model building will be open in explaining the DAMPE results.