This White Paper summarizes the discussions from a five-day workshop, involving 50 people from 18 countries, held in Cape Town, South Africa in February 2020. Convened by the International Astronomical Union's Executive Committee Working Group on Global Coordination of Ground and Space Astrophysics and sponsored by the Kavli Foundation, we discussed existing and potential bottlenecks for transient and multi-messenger astronomy, identifying eight broad areas of concern. Some of these are very similar to the challenges faced by many astronomers engaging in international collaboration, for example, data access policies, funding, theoretical and computational resources and workforce equity. Others, including, alerts, telescope coordination and target-of-opportunity implementation, are strongly linked to the time domain and are particularly challenging as we respond to transients. To address these bottlenecks we offer thirty-five specific recommendations, some of which are simply starting points and require development. These recommendations are not only aimed at collaborative groups and individuals, but also at the various organizations who are essential to making transient collaborations efficient and effective: including the International Astronomical Union, observatories, projects, scientific journals and funding agencies. We hope those involved in transient research will find them constructive and use them to develop collaborations with greater impact and more inclusive teams.

Observationally, fast radio bursts (FRBs) can be divided into repeating and apparently non-repeating (one-off) ones. It is unclear whether all FRBs repeat and whether there are genuine non-repeating FRBs. We attempt to address these questions using Monte Carlo simulations. We define a parameter $T_c$ at which the accumulated number of non-repeating sources becomes comparable to the total number of the repeating sources, which is a good proxy to denote the intrinsic repeater fraction among FRBs. Assuming that both types of sources exist, we investigate how the {\em observed} repeater fraction evolves with time for different parameters, including $T_c$, the repeating rate distribution power-law index $q$, the minimum repetition rate $r_{\rm 0,min}$, as well as the Weibull distribution index $k$ for repeating bursts. We find that unless $T_c \rightarrow \infty$ (i.e. there is no genuine non-repeating FRB source), the observed repeater fraction should increase with time first, reaching a peak, and then declines. The peak time $T_p$ and the peak fraction $F_{\rm r,obs,p}$ depend on $T_c$ and other repeating rate parameters. With the current data, one may pose a lower limit $T_c > (0.5-25)$ d for reasonable parameter values. Future continuous monitoring of FRBs using wide-field radio telescopes such as CHIME would measure or set a more stringent lower limit on $T_c$. The detection of a peak in the observed repeater fraction would disfavor the ansatz that ``all FRB sources repeat''.

The turbulence in the diffuse intergalactic medium (IGM) plays an important role in various astrophysical processes across cosmic time, but it is very challenging to constrain its statistical properties both observationally and numerically. Via the statistical analysis of turbulence along different sightlines toward a population of fast radio bursts (FRBs), we demonstrate that FRBs provide a unique tool to probe the intergalactic turbulence. We measure the structure function (SF) of dispersion measures (DMs) of FRBs to study the multi-scale electron density fluctuations induced by the intergalactic turbulence. The SF has a large amplitude and a Kolmogorov power-law scaling with angular separations, showing large and correlated DM fluctuations over a range of length scales. Given that the DMs of FRBs are IGM dominated, our result suggests that the intergalactic turbulence has a Kolmogorov power spectrum and an outer scale on the order of $100$ Mpc.

We use K-giant stars selected from the LAMOST DR5 to study the variation of the rotational velocity of the galactic halo at different space positions. Modelling the rotational velocity distribution with both the halo and disk components, we find that the rotational velocity of the halo population decreases almost linearly with increasing vertical distance to the galactic disk plane, $Z$, at fixed galactocentric radius, $R$. The samples are separated into two parts with $6<R<12$ kpc and $12<R<20$ kpc. We derive that the decreasing rates along $Z$ for the two subsamples are $-3.07\pm0.63$ and $-1.89\pm0.37$ km s$^{-1}$ kpc$^{-1}$, respectively. Compared with the TNG simulations, we suggest that this trend is probably caused by the interaction between the disk and halo. The results from the simulations show that only the oblate halo can provide a decreasing rotational velocity with an increasing $Z$. This indicates that the Galactic halo is oblate with galactocentric radius $R<20$ kpc. On the other hand, the flaring of the disk component (mainly the thick disk) is clearly traced by this study, with $R$ between 12 and 20 kpc, the disk can vertically extend to $6\sim10$ kpc above the disk plane. What is more interesting is that, we find the Gaia-Enceladus-Sausage (GES) component has a significant contribution only in the halo with $R<12$ kpc, i.e. a fraction of 23$-$47\%. While in the outer subsample, the contribution is too low to be well constrained.

GRB 200219A is a short gamma-ray burst (GRB) with an extended emission (EE) lasting $\sim 90$s. By analyzing data observed with the {\em Swift}/BAT and {\em Fermi}/GBM, we find that a cutoff power law (CPL) model can adequately well fit the spectra of the initial short pulse with $\rm E_{p}=1387^{+232}_{-134}$ keV. More interestingly, together with the EE component and early X-ray data, it exhibits a plateau emission smoothly connected with a $\sim t^{-1}$ segment, and followed by an extremely steep decay. The short GRB composed of those three segments is unique in the {\em Swift} era, and is very difficult to explain with the standard internal/external shock model of a black hole central engine, but could be consistent with the prediction of a magnetar central engine from the merger of an NS binary. We suggest that the plateau emission followed by a $\sim t^{-1}$ decay phase is powered by the spin-down of millisecond magnetar, which loses its rotation energy via GW quadrupole radiation. Then, the abrupt drop decay is caused by magnetar collapsing into black hole before switching to the EM-dominated emission. This is the first short GRB for which the X-ray emission has such an intriguing feature powered by magnetar via GW-dominated radiation. If this is the case, one can estimate the physical parameters of a magnetar, the GW signal powered by magnetar and the merger-nova emission are also discussed.

TianQin is a proposed space-based gravitational wave observatory. It is designed to detect the gravitational wave signals in the frequency range of 0.1 mHz -- 1 Hz. At a geocentric distance of $10^5$ km, the plasma in the earth magnetosphere will contribute as the main source of environmental noises. Here, we analyze the acceleration noises that are caused by the magnetic field of space plasma for the test mass of TianQin. The real solar wind data observed by the Advanced Composition Explorer are taken as the input of the magnetohydrodynamic simulation. The Space Weather Modeling Framework is used to simulate the global magnetosphere of the earth, from which we obtain the plasma and magnetic field parameters on the detector's orbits. We calculate the time series of the residual accelerations and the corresponding amplitude spectral densities on these orbit configurations. We find that the residual acceleration produced by the interaction between the TM's magnetic moment induced by the space magnetic field and the spacecraft magnetic field ($\bm{a}_{\rm M1}$) is the dominant term, which can approach $10^{-15}$ m/s$^2$/Hz$^{1/2}$ at $f \approx$ 0.2 mHz for the nominal values of the magnetic susceptibility ($\chi_{\rm m} = 10^{-5}$) and the magnetic shielding factor ($\xi_{\rm m} = 10$) of the test mass. The ratios between the amplitude spectral density of the acceleration noise caused by the space magnetic field and the preliminary goal of the inertial sensor are 0.38 and 0.08 at 1 mHz and 10 mHz, respectively. We discuss the further reduction of this acceleration noise by decreasing $\chi_{\rm m}$ and/or increasing $\xi_{\rm m}$ in the future instrumentation development for TianQin.

Jupiter's bright persistent polar aurora and Earth's dark polar region indicate that the planets' magnetospheric topologies are very different. High-resolution global simulations show that the reconnection rate at the interface between the interplanetary and jovian magnetic fields is too slow to generate a magnetically open, Earth-like polar cap on the timescale of planetary rotation, resulting in only a small crescent-shaped region of magnetic flux interconnected with the interplanetary magnetic field. Most of the jovian polar cap is threaded by helical magnetic flux that closes within the planetary interior, extends into the outer magnetosphere and piles-up near its dawnside flank where fast differential plasma rotation pulls the field lines sunward. This unusual magnetic topology provides new insights into Jupiter's distinctive auroral morphology.

The LIGO-Virgo and Fermi collaborations recently reported a possible joint detection of a sub-threshold gravitational wave (GW) event and a sub-threshold gamma-ray burst (GRB), GBM-190816, that occurred 1.57 s after the merger. We perform an independent analysis of the publicly available data and investigate the physical implications of this potential association. By carefully studying the following properties of GBM-190816 using Fermi/GBM data, including signal-to-noise ratio, duration, f-parameter, spectral properties, energetic properties, and its compliance with some GRB statistical correlations, we confirm that this event is likely a typical short GRB. Assuming its association with the sub-threshold GW event, the inferred luminosity is $1.47_{-1.04}^{+3.40} \times 10^{49}$ erg s$^{-1}$. Based on the available information of the sub-threshold GW event, we infer the mass ratio q of the compact binary as $q=2.26_{-1.43}^{+2.75}$ according to the reported range of luminosity distance. If the heavier compact object has a mass > 3 solar masses, q can be further constrained to $q=2.26_{-0.12}^{+2.75}$. The leading physical scenario invokes an NS-BH merger system with the NS tidally disrupted. Within this scenario, we constrain the physical properties of such a system to produce a GRB. The GW data may also allow an NS-BH system with no tidal disruption of the NS or a BH-BH merger. We apply the charged compact binary coalescence (cCBC) theory (for both a constant charge and an increasing charge for the merging members) to derive the model parameters to account for GBM-190816 and found that the required parameters are extreme. Finally, we argue that the fact that the observed GW-GRB delay time scale is comparable to that of GW170817/GRB 170817A suggests that the GW-GRB time delay of these two cases is mainly defined by the time scale for the jet to propagate to the GRB emission site.

The LIGO-Virgo and Fermi collaborations recently reported a possible joint detection of a sub-threshold gravitational wave (GW) event and a sub-threshold gamma-ray burst (GRB), GBM-190816, that occurred 1.57 s after the merger. We perform an independent analysis of the publicly available data and investigate the physical implications of this potential association. By carefully studying the following properties of GBM-190816 using Fermi/GBM data, including signal-to-noise ratio, duration, f-parameter, spectral properties, energetic properties, and its compliance with some GRB statistical correlations, we confirm that this event is likely a typical short GRB. Assuming its association with the sub-threshold GW event, the inferred luminosity is $1.47_{-1.04}^{+3.40} \times 10^{49}$ erg s$^{-1}$. Based on the available information of the sub-threshold GW event, we infer the mass ratio q of the compact binary as $q=2.26_{-1.43}^{+2.75}$ according to the reported range of luminosity distance. If the heavier compact object has a mass > 3 solar masses, q can be further constrained to $q=2.26_{-0.12}^{+2.75}$. The leading physical scenario invokes an NS-BH merger system with the NS tidally disrupted. Within this scenario, we constrain the physical properties of such a system to produce a GRB. The GW data may also allow an NS-BH system with no tidal disruption of the NS or a BH-BH merger. We apply the charged compact binary coalescence (cCBC) theory (for both a constant charge and an increasing charge for the merging members) to derive the model parameters to account for GBM-190816 and found that the required parameters are extreme. Finally, we argue that the fact that the observed GW-GRB delay time scale is comparable to that of GW170817/GRB 170817A suggests that the GW-GRB time delay of these two cases is mainly defined by the time scale for the jet to propagate to the GRB emission site.

Comprehensive photometric investigation was carried out to the early K-type contact binary - IL Cnc. A few light curves from both ground-based telescopes and the Kepler space telescope were obtained (or downloaded) and then analyzed in detail. They are mostly found to be asymmetric and there are even continuously changing O'Connell effect in the light curves from Kepler K2 data, suggesting the system to be highly active. Using the Wilson-Devinney code (version 2013), photometric solutions were derived and then compared. It is found that the calculation of the mass ratio is easily affected by the spot settings. Combining the radial velocities determined from LAMOST median resolution spectral data, the mass ratio of the binary components is found to be $M_2/M_1 = 1.76\pm 0.05$. The components are in shallow contact ($f\sim 9\%$) and have a temperature difference about $T_2 - T_1 = -280\pm 20$ K. The system is demonstrated to be W-subtype, which may be a common feature of K-type contact binaries. The masses of the binary components were estimated to be $M_1\sim 0.51$ M$_{sun}$ and $M_2\sim 0.90$ M$_{sun}$. The values are in good agreement with that deduced from the parallax data of Gaia. The results suggest that the primary component lacks luminosity compared with the zero main sequence. The H$\alpha$ spectral line of the primary component is found to be peculiar. Combining newly determined minimum light times with those collected from literature, the orbital period of IL Cnc is studied. It is found that the (O$-$C)s of the primary minima show sinusoidal variation while the secondary do not. The oscillation is more likely to be caused by the starspot activities. Yet this assumption needs more data to support.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

Observations of millisecond-long radio bursts from beyond the Milky Way have revealed a repeat pattern with a roughly 16-day period -- a finding that adds to the enigma of the origin of these bursts.

If at least one of the members of a compact binary coalescence is charged, the inspiral of the two members would generate a Poynting flux with an increasing power, giving rise to a brief electromagnetic counterpart temporally associated with the chirp signal of the merger (with possibly a small temporal offset), which we term as the {\em charged Compact Binary Coalescence} (cCBC) signal. We develop a general theory of cCBC for any mass and amount of charge for each member. Neutron stars (NSs), as spinning magnets, are guaranteed to be charged, so the cCBC signal should accompany all neutron star mergers. The cCBC signal is clean in a BH-NS merger with a small mass ratio ($q \equiv m_2/m_1 < 0.2$), in which the NS plunges into the BH as a whole, and its luminosity/energy can reach that of a fast radio burst if the NS is Crab-like. The strength of the cCBC signal in Extreme Mass Ratio Inspiral Systems (EMRIs) is also estimated.

A formula is proposed to quantitatively estimate the signal emission rate of Communicating Extra-Terrestrial Intelligent civilizations (CETIs) in the Galaxy. I suggest that one possible type of CETI signal would be brief radio bursts similar to fast radio bursts (FRBs). A dedicated search for FRB-like artificial signals in the Galaxy for decades may pose a meaningful upper limit on the emission rate of these signals by CETIs. The Fermi-Hart paradox is answered in terms of not having enough observing times for this and other types of signals. Whether humans should send FRB-like signals in the far future is briefly discussed.

We investigate the relationship between the optical variability of galaxies and their distances from the centre of the Virgo Cluster using Palomar Transient Factory data. We define the ratio between the standard deviation of the galaxy brightness and the mean value of the standard deviation as a measure of a galaxy's optical variability. A sample of 814 Virgo galaxies with 230263 observations shows a monotonically decreasing trend of optical variability with increasing clustercentric distance. The variability level inside the cluster is 3.2$\sigma$ higher than the level outside. We fit the variability with a linear function and find that the data reject a distance-independent model. We examine 217 background galaxies for comparison and find no significant trend in galaxy variability. We assess the relation with Monte Carlo simulation by rebuilding the brightness of each galaxy. The simulation shows a monotonically decreasing relation for member galaxy variability and a distance-independent relation for background galaxies. Our result is consistent with the theory that the cold gas flowing inwards the cluster centre fuels AGN activity. This work is a new implementation of the method using optical variability to investigate the relation between galaxies evolution and their environment.

Binaries play key roles in determining stellar parameters and exploring stellar evolution models. We build a catalog of 88 eclipsing binaries with spectroscopic information, taking advantage of observations from both the Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) and the Palomar Transient Factory (PTF) surveys. A software pipeline is constructed to identify binary candidates by examining their light curves. The orbital periods of binaries are derived from the Lomb-Scargle method. The key distinguishing features of eclipsing binaries are recognized by a new filter \textit{Flat Test}. We classify the eclipsing binaries by applying Fourier analysis on the light curves. Among all the binary stars, 13 binaries are identified as eclipsing binaries for the first time. The catalog contains information: position, primary eclipsing magnitude and time, eclipsing depth, the number of photometry and radial velocity observations, largest radial velocity difference, binary type, the effective temperature of observable star $T_{\rm eff}$, and surface gravity of observable star log \emph{g}. The false-positive probability is calculated by using both a Monte Carlo simulation and real data from the SDSS Stripe 82 Standard Catalog. The binaries in the catalog are mostly with a period of less than one day. The period distribution shows a 0.22-day cut-off which is consistent with the low probability of an eclipsing binary rotating with such a period.

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

When a magnetar magnetosphere is trigged by crustal deformations, an electric field $E_\parallel$ parallel to the magnetic field line would appear via Alvf\'en waves in the charge starvation region. If electron-positron pair bunches pre-exist, e.g., via some possible plasma instabilities, in the magnetosphere, these pairs will undergo charge separation in the $E_\parallel$ and in the meantime emit coherent curvature radiation. Following the approach of Yang & Zhang (2018), we find that the superposed curvature radiation becomes narrower due to charge separation, with the width of spectrum depending on the separation between the electron and positron clumps. This mechanism can interpret the narrow spectra of FRBs, in particular, the Galactic FRB 200428 recently detected in association with a hard X-ray burst from the Galactic magnetar SGR J1935+2154.

Recent observations with the Chandra X-ray telescope continue to detect X-ray emission from the transient GW170817. In a total exposure of 96.6 ks, performed between March 9 and March 16 2020 (935 d to 942 d after the merger), a total of 8 photons are measured at the source position, corresponding to a significance of about 5 sigma. Radio monitoring with the Australian Telescope Compact Array (ATCA) shows instead that the source has faded below our detection threshold (<33 uJy, 3 sigma). By assuming a constant spectral index beta=0.585, we derive an unabsorbed X-ray flux of approximately 1.4E-15 erg/cm^2/s, higher than earlier predictions, yet still consistent with a simple structured jet model. We discuss possible scenarios that could account for prolonged emission in X-rays. The current dataset appears consistent both with energy injection by a long-lived central engine and with the onset of a kilonova afterglow, arising from the interaction of the sub-relativistic merger ejecta with the surrounding medium. Long-term monitoring of this source will be essential to test these different models.

We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

The discovery of a fast radio burst (FRB) in our galaxy associated with a magnetar (neutron star with strong magnetic field) has provided a critical piece of information to help us finally understand these enigmatic transients. We show that the volumetric rate of Galactic-FRB like events is consistent with the faint end of the cosmological FRB rate, and hence they most likely belong to the same class of transients. The Galactic FRB had an accompanying X-ray burst but many X-ray bursts from the same object had no radio counterpart. Their relative rates suggest that for every FRB there are roughly 100 to 1000 X-ray bursts. The radio lightcurve of the galactic FRB had two spikes separated by 30 ms in the 400-800 MHz frequency band. This is an important clue and highly constraining of the class of models where the radio emission is produced outside the light-cylinder of the magnetar. We suggest that magnetic disturbances close to the magnetar surface propagate to a distance of a few tens of neutron star radii where they damp and produce radio emission. The coincident hard X-ray spikes associated with the two FRB pulses seen in this burst and the flux ratio between the two frequency bands can be understood in this scenario. This model provides a unified picture for faint bursts like the Galactic FRB as well as the bright events seen at cosmological distances.

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

Gamma Ray Burst (GRBs) are empirically classified as long-duration GRBs (LGRBs, $>$ 2s) and short-duration GRBs (SGRBs, $<$ 2s). Physically they can be grouped into two distinct progenitor categories: those originating from collapse of massive stars (also known as Type II) and those related to mergers of compact stars (also known as Type I). Even though most LGRBs are Type II and most SGRBs are Type I, the duration criterion is not always reliable to determine the physical category of a certain GRB. Based on our previous comprehensive study of the multi-wavelength properties of long and short GRBs, here we utilize the Naive Bayes method to physically classify GRBs as Type I and Type II GRBs based on multi-wavelength criteria. It results in 0.5\% training error rate and 1\% test error rate. Moreover, there is a gap [-1.2, -0.16] in the distribution of the posterior Odds, $\log O({\rm II:I})$, the Type II to Type I probability ratio. Therefore, we propose to use ${\cal O} = \log O({\rm II:I})+0.7$ as the parameter to classify GRBs into Type I ($<0$) or Type II ($>0$). The only confirmed Type I GRB, GRB 170817A, has log $O({\rm II:I})=-10$. According to this criterion, the supernova-less long GRBs 060614 and 060505 belong to Type I, and two controversial short GRBs 090426 and 060121 belong to Type II.

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

Tidal disruption events (TDE) have been considered as cosmic-ray and neutrino sources for a decade. We suggest two classes of new scenarios for high-energy multi-messenger emission from TDEs that do not have to harbor powerful jets. First, we investigate high-energy neutrino and gamma-ray production in the core region of a supermassive black hole. In particular, we show that about 1-100 TeV neutrinos and MeV gamma-rays can efficiently be produced in hot coronae around an accretion disk. We also study the consequences of particle acceleration in radiatively inefficient accretion flows (RIAFs). Second, we consider possible cosmic-ray acceleration by disk-driven winds or interactions between tidal streams, and show that subsequent hadronuclear and photohadronic interactions inside the TDE debris lead to GeV-PeV neutrinos and sub-GeV cascade gamma-rays. We demonstrate that these models should be accompanied by soft gamma-rays or hard X-rays, which can be used for future observational tests. Although this work aims to present models of non-jetted high-energy emission, we discuss the implications of the TDE AT2019dsg that might coincide with the high-energy neutrino IceCube-191001A, by considering the corona, RIAF, hidden sub-relativistic wind, and hidden jet models. It is not yet possible to be conclusive about their physical association and the expected number of neutrinos is typically much less than unity. We find that the most optimistic cases of the corona and hidden wind models could be consistent with the observation of IceCube-191001A, whereas jet models are unlikely to explain the multi-messenger observations.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

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 the results of a detailed investigation of the prompt and afterglow emission in the HESS detected GRB 190829A. Swift and Fermi observations of the prompt phase of this GRB reveal two isolated sub-bursts or episodes, separated by a quiescent phase. The energetic and the spectral properties of the first episode are in stark contrast to the second. The first episode, which has a higher spectral peak of $\sim 120\:\text{keV}$ and a low isotropic energy $\sim 10^{50}\:\text{erg}$ is an outlier to the Amati correlation and marginally satisfies the Yonetoku correlation. However, the energetically dominant second episode has lower peak energy and is consistent with the above correlations. We compared this GRB to other low luminosity GRBs (LLGRBs). Prompt emission of LLGRBs also indicates a relativistic shock breakout origin of the radiation. For GRB 190829A, some of the properties of a shock breakout origin are satisfied. However, the absence of an accompanying thermal component and energy above the shock breakout critical limit precludes a shock breakout origin. In the afterglow, an unusual long-lasting late time flare of duration $\sim 10^4\:\text{s}$ is observed. We also analyzed the late-time \fermi-LAT emission that encapsulates the H.E.S.S. detection. Some of the LAT photons are likely to be associated with the source. All the above observational facts suggest GRB 190829A is a peculiar low luminosity GRB that is not powered by a shock breakout, and with an unusual rebrightening due to a patchy emission or a refreshed shock during the afterglow. Furthermore, our results show that TeV energy photons seem common in both high luminosity GRBs and LLGRBs.

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

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

The reported association of high-energy neutrino event IceCube-170922A and blazar TXS 0506+056 has sparkled discussion about blazars as sources of cosmic neutrinos. In this paper, we use publicly released IceCube data and blazar locations from Roma-BzCat to test spacial correlation between neutrino events and blazar locations. We also scrutinize the correlation between $\gamma$-ray flux and neutrino flux of blazars by applying a temporal filter onto the data based on Fermi monitored source list. We find no compelling evidence to prove blazars as the main source of cosmic neutrinos, as known before the detected IceCube-170922A / TXS 0506+056 association. While we do not rule out the association between IceCube-170922A and TXS 0506+056, the significance level we obtained is not high enough to claim a strong association. If such an association is real, a special physical condition is desired to allow a small fraction of blazars to become bright neutrino sources.

Recent localizations of Fast Radio Bursts and identifications of their host galaxies confirmed their extragalactic origin. While FRB 121102 resides in the bright region of a dwarf star forming galaxy, other FRBs reside in more massive galaxies and are related to older stellar populations. We compare the host galaxy properties of the five FRBs with those of several types of stellar transients: from young to old population, long duration gamma ray bursts (LGRBs), superluminous supernovae (SLSNe), Type Ib/Ic supernovae (SN Ibc), Type II supernovae (SN II), type Ia supernovae (SN Ia), and short duration gamma ray bursts (SGRBs). We find that as a whole sample, the stellar mass and star formation rate of the five FRB host galaxies prefer a medium to old population, and are against a young population similar to LGRBs and SLSNe by a null probability 0.05. Individually, the FRB 121102 host is consistent with that of young population objects, especially SLSNe; the FRB 180924 environment is similar to that of SGRBs; and the FRB 190523 environment is similar to those of SN Ia or SGRBs. These results are consistent with the magnetar engine model for FRBs, if both magnetars produced from extreme explosions (GRBs/SLSNe) and from regular channels (e.g. those producing Galactic magnetars) can produce FRBs.

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

Low-surface-brightness galaxies(LSBGs) are defined as galaxies that are fainter than dark night sky and are important for studying our universe. Particularly, edge-on galaxies are useful for the study of rotational velocity and dynamical properties of galaxies. Hence here we focus on searching for edge-on LSBGs. In order to find these edge-on dim galaxies, a series of effects caused by inclination, including the surface brightness profile, internal extinction, and scale length, have been corrected. In this work, we present a catalog of 281 edge-on LSBG candidates, which are selected from the cross-match between SDSS DR7 and the 40% ALFALFA catalog. We also present the properties of these edge-on LSBG candidates including absolute magnitude, central surface brightness, B-V color, scale length, and relative thickness. Our result suggests that the correction of inclination effects is very important for obtaining a complete sample of LSBGs.

The physical origin of fast radio bursts (FRBs) is unknown. Young magnetars born from gamma-ray bursts (GRBs) have been suggested to be a possible central engine of FRBs. We test such a hypothesis by systematically searching for GRB-FRB spatial associations from 110 FRBs and 1440 GRBs. We find that one FRB event, FRB 171209, discovered by the Parkes telescope is spatially coincident with a historical long-duration GRB 110715A at $z=0.82$. The afterglow of GRB 110715A is consistent with being powered by a millisecond magnetar. The extragalactic dispersion measure of FRB 171209 is in excess of that contributed by the intergalactic medium, which can be interpreted as being contributed by a young supernova remnant associated with the GRB. Overall, the significance of the association is $(2.28 - 2.55) \sigma$. If the association is indeed physical, our result suggests that the magnetars associated with long GRBs can be the progenitors of at least some FRBs.

The luminosity function of Fast Radio Bursts (FRBs), defined as the event rate per unit cosmic co-moving volume per unit luminosity, may help to reveal the possible origins of FRBs and design the optimal searching strategy. With the Bayesian modelling, we measure the FRB luminosity function using 46 known FRBs. Our Bayesian framework self-consistently models the selection effects, including the survey sensitivity, the telescope beam response, and the electron distributions from Milky Way / the host galaxy / local environment of FRBs. Different from the previous companion paper, we pay attention to the FRB event rate density and model the event counts of FRB surveys based on the Poisson statistics. Assuming a Schechter luminosity function form, we infer (at the 95% confidence level) that the characteristic FRB event rate density at the upper cut-off luminosity $L^*=2.9_{-1.7}^{+11.9}\times10^{44}\,\rm erg\, s^{-1}$ is $\phi^*=339_{-313}^{+1074}\,\rm Gpc^{-3}\, yr^{-1}$, the power-law index is $\alpha=-1.79_{-0.35}^{+0.31}$, and the lower cut-off luminosity is $L_0\le9.1\times10^{41}\,\rm erg\, s^{-1}$. The event rate density of FRBs is found to be $3.5_{-2.4}^{+5.7}\times10^4\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{42}\,\rm erg\, s^{-1}$, $5.0_{-2.3}^{+3.2}\times10^3\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{43}\,\rm erg\, s^{-1}$, and $3.7_{-2.0}^{+3.5}\times10^2\,\rm Gpc^{-3}\, yr^{-1}$ above $10^{44}\,\rm erg\, s^{-1}$. As a result, we find that, for searches conducted at 1.4 GHz, the optimal diameter of single-dish radio telescopes to detect FRBs is 30-40 m. The possible astrophysical implications of the measured event rate density are also discussed in the current paper.

We present results from VLBI observations of continuum and OH line emission in IRAS 02524+2046 and also arcsecond-scale radio properties of this galaxy using VLA archive data. We found that there is no significant detection of radio continuum emission from VLBI observations. The arcsecond-scale radio images of this source show no clear extended emission, the total radio flux density at L and C band are around 2.9 mJy and 1.0 mJy respectively, which indicate a steep radio spectral index between the two band. Steep spectral index, low brightness temperature and high $q$-ratio (the FIR to the radio flux density), which are three critical indicators in classification of radio activity in the nuclei of galaxies, are all consistent with the classification of this source as a starburst galaxy from its optical spectrum. The high-resolution line profile show that both of \textbf{the 1665 and 1667 MHz OH maser} line have been detected which show three and two clear components respectively. The channel maps show that the maser emission are distributed in a region $\sim$ 210 pc $\times$ 90 pc, the detected maser components at different region show similar double spectral feature, which might be an evidence that this galaxy is at a stage of major merger as seen from the optical morphology.

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

The physical origin of fast radio bursts (FRBs) is still unknown. Multiwavelength and polarization observations of an FRB source would be helpful to diagnose its progenitor and environment. So far only the first repeating source FRB 121102 appears to be spatially coincident with a persistent radio emission. Its bursts also have very large values of the Faraday rotation measure (RM) i.e., $|{\rm RM}|\sim10^5~{\rm rad~m^{-2}}$. We show that theoretically there should be a simple relation between RM and the luminosity of the persistent source of an FRB source if the observed RM mostly arises from the persistent emission region. FRB 121102 follows this relation given that the magnetic field in the persistent emission region is highly ordered and that the number of relativistic electrons powering the persistent emission is comparable to that of non-relativistic electrons that contribute to RM. The non-detections of persistent emission sources from all other localized FRB sources are consistent with their relatively small RMs ($\left|{\rm RM}\right|\lesssim{\rm a~few}\times100~{\rm rad~m^{-2}}$) according to this relation. Based on this picture, the majority of FRBs without a large RM are not supposed to be associated with bright persistent sources.

We show that the periodic FRB 180916.J0158+65 can be interpreted by invoking an interacting neutron star binary system with an orbital period of $\sim 16$ days. The FRBs are produced by a highly magnetized pulsar, whose magnetic field is ``combed'' by the strong wind from a companion star, either a massive star or a millisecond pulsar. The FRB pulsar wind retains a clear funnel in the companion's wind that is otherwise opaque to induced Compton or Raman scatterings for repeating FRB emission. The 4 day active window corresponds to the time when the funnel points toward Earth. The interaction also perturbs the magnetosphere of the FRB pulsar and may trigger emission of FRBs. We derive the physical constraints on the comb and the FRB pulsar from the observations and estimate the event rate of FRBs. In this scenario, a lower limit on the period of observable FRBs is predicted. We speculate that both the intrinsic factors (strong magnetic field and young age) and the extrinsic factor (interaction) may be needed to generate FRBs in neutron star binary systems.

The post-merger product of the first binary neutron star merger event detected in gravitational waves, GW170817, depends on neutron star equation of state (EoS) and is not well determined. We generally discuss the constraints one may pose on the maximum mass of a non-spinning neutron star, $M_{\rm TOV}$, based on the observations and some EoS-independent universal relations of rapidly-spinning neutron stars. If the merger product is a black hole after a brief hypermassive neutron star (HMNS) phase, we derive $M_{\rm TOV} < 2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04}) M_{\odot}$ at the 2$\sigma$ (1$\sigma$) level. The cases for a massive neutron star (MNS), either a supra-massive neutron star (SMNS) or even a stable neutron star (SNS), are also allowed by the data. We derive $2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04} M_{\odot}) \leq M_{\rm TOV}< 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04}) M_{\odot}$ for the SMNS case and $M_{\rm TOV} \geq 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04})M_{\odot}$ for the SNS case, at the $2\sigma$ ($1\sigma$) confidence level. In the MNS cases, we also discuss the constraints on the neutron star parameters (the dipolar magnetic field strength at the surface $B_p$ and the ellipticity $\epsilon$) that affect the spindown history, by considering different MNS survival times, e.g. 300 s, 1 d, and 155 d after the merger, as suggested by various observational arguments. We find that once an SMNS is formed, without violating the EM observational constraints, there always exist a set of ($B_p, \epsilon$) parameters that allow the SMNS to survive for 300s, 1 d, 155 d, or even longer.

In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

Fast radio bursts (FRBs) are mysterious radio transients whose physical origin is still unknown. Within a few astronomical units near an FRB source, the electric field of the electromagnetic wave is so large that the electron oscillation velocity becomes relativistic, which makes the classical Thomson scattering theory and the linear plasma theory invalid. We discuss FRBs as strong waves interacting with the ambient medium, in terms of both electron motion properties and plasma properties. Several novel features are identified. 1. The cross section of Thomson scattering is significantly enhanced for the scattering photons. 2. On the other hand, because of the nonlinear plasma properties in strong waves, the near-source plasma is more transparent and has a smaller effective dispersion measure (DM) contribution to the observed value. For a repeating FRB source, the brighter bursts would have somewhat smaller DMs contributed by the near-source plasma. 3. The radiation beam undergoes relativistic self-focusing in a dense plasma, the degree of self-focusing (or squeezing) depends on the plasma density. Such a squeezing effect would affect the collimation angle and the true event rate of FRBs. 4. When an FRB propagates in a nearby ambient plasma, a wakefield wave in the plasma will be generated by the ponderomotive force of the FRB and accelerates electrons in the ambient medium. However, such an effect is too weak to be observationally interesting.

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

Recent observations of repeating fast radio bursts (FRBs) suggest that some FRBs reside in an environment consistent with that of binary neutron star (BNS) mergers. The bursting rate for repeaters could be very high and the emission site is likely from a magnetosphere. We discuss a hypothesis of producing abundant repeating FRBs in BNS systems. Decades to centuries before a BNS system coalesces, the magnetospheres of the two neutron stars start to interact relentlessly. Abrupt magnetic reconnection accelerates particles, which emit coherent radio waves in bunches via curvature radiation. FRBs are detected as these bright radiation beams point towards Earth. This model predicts quasi-periodicity of the bursts at the rotation periods of the two merging neutron stars (tens of milliseconds and seconds, respectively) as well as the period of orbital motion (of the order of 100 s). The bursting activities are expected to elevate with time as the two neutron stars get closer. The repeating FRB sources should be gravitational wave (GW) sources for space-borne detectors such as LISA if the sources are close enough, and will be eventually detected by ground-based detectors when the two neutron stars coalesce.

Growing evidence suggests that synchrotron radiation plays a significant role in shaping the spectra of most $\gamma$-ray bursts. The relativistic jets producing them likely carry a significant fraction of energy in the form of a Poynting flux.

In order to measure the Universal Time (UT1) in real time, National Time Service Center (NTSC) has built a VGOS-like (VLBI Global Observing System) broadband VLBI network, which includes three 13-m radio telescopes located in Jilin, Sanya and Kashi, and a data analysis center in Xi'an. Each station is equipped with a highly stable hydrogen atomic clock and a self-developed VLBI backend, and is co-located with two GPS receivers. This VGOS-like VLBI network may play an important role in improving the Chinese broadband VLBI technology and making valuable contributions to domestic VLBI measurements of UT1. In this paper, we introduce the specifications of this VLBI network, and present the UT1 measurements at C-band conducted in 2018 using the Jilin-Kashi baseline of this network. The comparisons between our UT1 estimates and those provided by IERS suggest that the NTSC VLBI network is capable to determine UT1 accurate at the level of 58.8 microseconds.

Radio frequency interference (RFI) detection and excision are key steps in the data-processing pipeline of the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Because of its high sensitivity and large data rate, FAST requires more accurate and efficient RFI flagging methods than its counterparts. In the last decades, approaches based upon artificial intelligence (AI), such as codes using convolutional neural networks (CNNs), have been proposed to identify RFI more reliably and efficiently. However, RFI flagging of FAST data with such methods has often proved to be erroneous, with further manual inspections required. In addition, network construction as well as preparation of training data sets for effective RFI flagging has imposed significant additional workloads. Therefore, rapid deployment and adjustment of AI approaches for different observations is impractical to implement with existing algorithms. To overcome such problems, we propose a model called RFI-Net. With the input of raw data without any processing, RFI-Net can detect RFI automatically, producing corresponding masks without any alteration of the original data. Experiments with RFI-Net using simulated astronomical data show that our model has outperformed existing methods in terms of both precision and recall. Besides, compared with other models, our method can obtain the same relative accuracy with fewer training data, thus reducing the effort and time required to prepare the training data set. Further, the training process of RFI-Net can be accelerated, with overfittings being minimized, compared with other CNN codes. The performance of RFI-Net has also been evaluated with observing data obtained by FAST and the Bleien Observatory. Our results demonstrate the ability of RFI-Net to accurately identify RFI with fine-grained, high-precision masks that required no further modification.

The joint detection of GW 170817 and GRB 170817A indicated that at least a fraction of short gamma ray bursts (SGRBs) originate from binary neutron star (BNS) mergers. One possible remnant of a BNS merger is a rapidly rotating, strongly magnetized neutron star, which has been discussed as one possible central engine for GRBs. For a rapidly rotating magnetar central engine, the deposition of the rotation energy into the ejecta launched from the merger could lead to bright radio emission. The brightness of radio emission years after an SGRB would provide an estimation of the kinetic energy of ejecta and hence, a possible constraint on the BNS merger product. We perform a more detailed calculation on the brightness of radio emission from the interaction between the merger ejecta and circumburst medium in the magnetar scenario, invoking several important physical processes such as generic hydrodynamics, relativistic effects, and the deep Newtonian phase. We use the model to constrain the allowed parameter space for 15 SGRBs that have late radio observations. Our results show that an injected energy $E_{\rm inj} \sim 10^{52}$ erg is allowed for all the cases, which suggests that the possibility of a supra-massive or hyper-massive neutron star remnant is not disfavored by the available radio data.

The IceCube collaboration reported a $\sim 3.5\sigma$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 0506+56 during a $\sim$6 month period in 2014-2015, as well as the ($\sim3\sigma$) detection of a high-energy muon neutrino during an electromagnetic flare in 2017. We explore the possibility that the 2014-2015 neutrino excess and the 2017 multi-messenger flare are both explained in a common physical framework that relies on the emergence of a relativistic neutral beam in the blazar jet due to interactions of accelerated cosmic rays (CRs) with photons. We demonstrate that the neutral beam model provides an explanation for the 2014-2015 neutrino excess without violating X-ray and $\gamma$-ray constraints, and also yields results consistent with the detection of one high-energy neutrino during the 2017 flare. If both neutrino associations with TXS 05065+056 are real, our model requires that (i) the composition of accelerated CRs is light, with a ratio of helium nuclei to protons $\gtrsim5$, (ii) a luminous external photon field ($\sim 10^{46}$ erg s$^{-1}$) variable (on year-long timescales) is present, and (iii) the CR injection luminosity as well as the properties of the dissipation region (i.e., Lorentz factor, magnetic field, and size) vary on year-long timescales.

Two bright X-ray transients were reported from the Chandra Deep Field South archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the "free zone" where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$, $P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$, $\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent" ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers.

Two bright X-ray transients were reported from the Chandra Deep Field South archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the "free zone" where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$, $P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$, $\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent" ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers.

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

Cosmic explosions called $\gamma$-ray bursts are the most energetic bursting events in the universe. Observations of extremely high-energy emission from two $\gamma$-ray bursts provide a new way to study these gigantic explosions.

We present trigonometric parallax and proper motion measurements toward 22 GHz water and 6.7 GHz methanol masers in 16 high-mass star-forming regions. These sources are all located in the Scutum spiral arm of the Milky Way. The observations were conducted as part of the Bar and Spiral Structure Legacy (BeSSeL) survey. A combination of 14 sources from a forthcoming study and 14 sources from the literature, we now have a sample of 44 sources in the Scutum spiral arm, covering a Galactic longitude range from 0$^\circ$ to 33$^\circ$. A group of 16 sources shows large peculiar motions of which 13 are oriented toward the inner Galaxy. A likely explanation for these high peculiar motions is the combined gravitational potential of the spiral arm and the Galactic bar.

GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

The characteristics of the spectral evolution of the prompt emission of gamma-ray bursts (GRBs), which are closely related to the radiation mechanism (synchrotron or photosphere), are still an unsolved subject. Here, by performing the detailed time-resolved spectral fitting of GRB 131231A, which has a very bright and well-defined single pulse, some interesting spectral evolution features have been found. (i) Both the low-energy spectral index $\alpha$ and the peak energy $E_{\rm p}$ exhibit the "flux-tracking" pattern ("double-tracking" characteristics). (ii) The parameter relations, i.e., $F$ (the energy flux)-$\alpha$, $F$-$E_{\rm p}$, and $E_{\rm p}$-$\alpha$, along with the analogous Yonetoku $E_{\rm p}$-$L_{\gamma,\rm iso}$ relation for the different time-resolved spectra, show strong monotonous (positive) correlations, both in the rising and the decaying phases. (iii) The values of $\alpha$ do not exceed the synchrotron limit ($\alpha$= -2/3) in all slices across the pulse, favoring the synchrotron origin. We argue that the one-zone synchrotron emission model with the emitter streaming away at a large distance from the central engine can explain all of these special spectral evolution characteristics.

Low-resolution spectra are proved competitive to high-resolution spectra in determining many stellar labels at comparable precision. It is useful to consider the spectral information content when assessing the capability of a stellar spectrum in deriving precise stellar labels. In this work, we quantify the information content brought by the LAMOST-II medium-resolution spectroscopic survey (MRS) using the gradient spectra and the coefficients-of-dependence (CODs). In general, the wavelength coverage of the MRS well constrains the stellar labels but the sensitivities of different stellar labels vary with spectral types and metallicity of the stars of interest and, therefore, affect the performance of the stellar label determination from the MRS spectra. Applying the SLAM to the synthetic spectra which mimic the MRS data, we find the precision of the fundamental stellar parameters Teff, logg and [M/H] are better when combining both the blue and red bands of the MRS. This is especially important for warm stars since the H$\alpha$ line located in the red part plays a more important role in determining the effective temperature for warm stars. With blue and red parts together, we are able to reach similar performance to the low-resolution spectra except for warm stars. However, at [M/H]$\sim-2.0$ dex, the uncertainties of fundamental stellar labels estimated from MRS are substantially larger than those from low-resolution spectra. We also tested the uncertainties of Teff, logg and [M/H] of from MRS data induced from the radial velocity mismatch and find that a mismatch of about 1 km s$^{-1}$, which is typical for LAMOST MRS data, would not significantly affect the stellar label estimates. At last, reference precision limits are calculated using synthetic gradient spectra, according to which we expect abundances of at least 17 elements to be measured precisely from MRS spectra.

We investigate the long GRB140629A through multiwavelength observations, which cover optical, infrared and X-rays between 40s and 3yr after the burst, to derive the properties of the dominant jet and its host galaxy. Polarisation observations by the MASTER telescope indicate that this burst is weakly polarised. The optical spectrum contains absorption features, from which we confirm the redshift of the GRB as originating at z=2.276. We performed spectral fitting of the X-rays to optical afterglow data and find there is no strong spectral evolution. We determine the hydrogen column density to be 7.2x10^21cm^-2 along the line of sight. The afterglow in this burst can be explained by a blast wave jet with a long-lasting central engine expanding into a uniform medium in the slow cooling regime. At the end of energy injection, a normal decay phase is observed in both the optical and X-ray bands. An achromatic jet break is also found in the afterglow light curves 0.4d after trigger. We fit the multiwavelength data simultaneously with a model based on a numerical simulation and find that the observations can be explained by a narrow uniform jet in a dense environment with an opening angle of 6.7deg viewed 3.8deg off-axis, which released a total energy of 1.4x10^54erg. Using the redshift and opening angle, we find GRB 140629A follows both the Ghirlanda and Amati relations. From the peak time of the light curve, identified as the onset of the forward shock (181s after trigger), the initial Lorentz factor is constrained in the range 82-118. Fitting the host galaxy photometry, we find the host to be a low mass, star-forming galaxy with a star formation rate of logSFR=1.1^+0.9_-0.4Myr^-1. We obtain a value of the neutral hydrogen density by fitting the optical spectrum, logN(HI)=21.0+-0.3, classifying this host as a damped Lyman-alpha. High ionisation lines are also detected in the spectrum.

We compile and analyze ~200 trigonometric parallaxes and proper motions of molecular masers associated with very young high-mass stars. These measurements strongly suggest that the Milky Way is a four-arm spiral. Fitting log-periodic spirals to the locations of the masers, allows us to significantly expand our view of the structure of the Milky Way. We present an updated model for its spiral structure and incorporate it into our previously published parallax-based distance-estimation program for sources associated with spiral arms. Modeling the three-dimensional space motions yields estimates of the distance to the Galactic center, Ro = 8.15 +/- 0.15 kpc, the circular rotation speed at the Sun's position, To = 236 +/- 7 km/s, and the nature of the rotation curve. Our data strongly constrain the full circular velocity of the Sun, To + Vsun = 247 +/- 4 km/s, and its angular velocity, (To + Vsun)/Ro = 30.32 +/- 0.27 km/s/kpc. Transforming the measured space motions to a Galactocentric frame which rotates with the Galaxy, we find non-circular velocity components typically about 10 km/s. However, near the Galactic bar and in a portion of the Perseus arm, we find significantly larger non-circular motions. Young high-mass stars within 7 kpc of the Galactic center have a scale height of only 19 pc and, thus, are well suited to define the Galactic plane. We find that the orientation of the plane is consistent with the IAU-defined plane to within +/-0.1 deg., and that the Sun is offset toward the north Galactic pole by Zsun = 5.5 +/- 5.8 pc. Accounting for this offset places the central supermassive black hole, Sgr A*, in the midplane of the Galaxy. Using our improved Galactic parameters, we predict the Hulse-Taylor binary pulsar to be at a distance of 6.54 +/- 0.24 kpc, assuming its orbital decay from gravitational radiation follows general relativity.

The LAMOST survey has provided 9 million spectra in its Data Release 5 (DR5) at R$\sim$1800. Extracting precise stellar labels is crucial for such a large sample. In this paper, we report the implementation of the Stellar LAbel Machine (SLAM), which is a data-driven method based on Support Vector Regression (SVR), a robust non-linear regression technique. Thanks to the capability to model highly non-linear problems with SVR, SLAM generally can derive stellar labels over a wide range of spectral types. This gives it a unique capability compared to other popular data-driven methods. To illustrate this capability, we test the performance of SLAM on stars ranging from Teff$\sim$4000 to $\sim$8000 K trained on LAMOST spectra and stellar labels. At g-band signal-to-noise ratio (SNRg) higher than 100, the random uncertainties of Teff, logg and [Fe/H] are 50 K, 0.09 dex, and 0.07 dex, respectively. We then set up another SLAM model trained by APOGEE and LAMOST common stars to demonstrate its capability of dealing with high dimensional problems. The spectra are from LAMOST DR5 and the stellar labels of the training set are from APOGEE DR15, including Teff, logg, [M/H],[$\alpha$/M], [C/M], and [N/M]. The cross-validated scatters at SNRg$\sim$100 are 49 K, 0.10 dex, 0.037 dex,0.026 dex, 0.058 dex, and 0.106 dex for these parameters, respectively. This performance is at the same level as other up-to-date data-driven models. As a byproduct, we also provide the latest catalog of $\sim$1 million LAMOST DR5 K giant stars with SLAM-predicted stellar labels in this work.

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $\alpha$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$\alpha$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $\alpha$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream.

Recently two fast X-ray transients (XT1 and XT2) have been reported from the search in the Chandra Deep Field (CDF) data. Each transient shows an initial plateau lasting around hundreds to thousands seconds followed by a rapid decay in the light curve. In particular, CDF-S XT2 is found to be associated with a galaxy at redshift z= 0.738 and was explained as a counterpart of a binary neutron-star merger event. In this paper, motivated by the short duration and decay slopes of the two transients, we consider an alternative interpretation in which both events are accretion-driven flares from tidal disruption of white dwarfs by intermediate-mass black holes. We derive a theoretical model of the accretion rate history, and find that it fits the observed X-ray light curves well. The extremely super-Eddington peak luminosity of XT2 can be explained by the beaming effect of the system, likely in the form of a jet.

Transrelativistic supernovae (SNe), which are likely driven by central engines via jets or winds, have been among candidate sources of ultrahigh-energy cosmic rays (UHECRs). We investigate acceleration and survival of UHECR nuclei in the external reverse shock scenario. With composition models used in Zhang et al. (2018), we calculate spectra of escaping cosmic rays and secondary neutrinos. If their local rate is $\sim1$% of the core-collapse supernova rate, the observed UHECR spectrum and composition can be explained with the total cosmic-ray energy ${\mathcal E}_{\rm cr}\sim10^{51}$ erg. The maximum energy of UHECR nuclei can reach $\sim 10^{20}-{10}^{21}\rm~eV$. The diffuse flux of source neutrinos is predicted to be $\sim 10^{-11} - 10^{-10}~{\rm GeV}~{\rm cm}^{-2}~{\rm s}^{-1}~{\rm sr}^{-1}$ in the 0.1-1 EeV range, satisfying nucleus-survival bounds. The associated cosmogenic neutrino flux is calculated, and shown to be comparable or even higher than the source neutrino flux. These ultrahigh-energy neutrinos can be detected by ultimate detectors such as the Giant Radio Askaryan Neutrino Detector and Probe Of Extreme Multi-Messenger Astrophysics.

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

Early and late multiwavelength observations play an important role in determining the nature of the progenitor, circumburst medium, physical processes and emitting regions associated to the spectral and temporal features of bursts. GRB 180720B is a long and powerful burst detected by a large number of observatories in multiwavelenths that range from radio bands to sub-TeV gamma-rays. The simultaneous multiwavelength observations were presented over multiple periods of time beginning just after the trigger time and extending for more than 30 days. The temporal and spectral analysis of Fermi LAT observations suggests that it presents similar characteristics to other bursts detected by this instrument. Coupled with X-ray and optical observations, the standard external-shock model in a homogeneous medium is favored by this analysis. The X-ray flare is consistent with the synchrotron self-Compton (SSC) model from the reverse-shock region evolving in a thin shell and long-lived LAT, X-ray and optical data with the standard synchrotron forward-shock model. The best-fit parameters derived with the Markov chain Monte Carlo simulations indicate that the outflow is endowed with magnetic fields and that the radio observations are in the self-absorption regime. The SSC forward-shock model with our parameters can explain the LAT photons beyond the synchrotron limit as well as the emission recently reported by the HESS Collaboration.

The quasi-thermal components found in many Fermi gamma-ray bursts (GRBs) imply that the photosphere emission indeed contributes to the prompt emission of many GRBs. But whether the observed spectra empirically fitted by the Band function or cutoff power law, especially the spectral and peak energy ($E_{p}$) evolutions can be explained by the photosphere emission model alone needs further discussion. In this work, we investigate in detail the time-resolved spectra and $E_{p}$ evolutions of photospheric emission from a structured jet, with an inner-constant and outer-decreasing angular Lorentz factor profile. Also, a continuous wind with a time-dependent wind luminosity has been considered. We show that the photosphere spectrum near the peak luminosity is similar to the cutoff power-law spectrum. The spectrum can have the observed average low-energy spectral index $\alpha $ $ \sim -1$, and the distribution of the low-energy spectral index in our photosphere model is similar to that observed ($-2\lesssim $ $\alpha \lesssim 0$). Furthermore, the two kinds of spectral evolutions during the decay phase, separated by the width of the core ($\theta _{c}$), are consistent with the time-resolved spectral analysis results of several Fermi multi-pulse GRBs and single-pulse GRBs, respectively. Also, for this photosphere model we can reproduce the two kinds of observed $E_{p}$ evolution patterns rather well. Thus, by considering the photospheric emission from a structured jet, we reproduce the observations well for the GRBs best fitted by the cutoff power-law model for the peak-flux spectrum or the time-integrated spectrum.

The analogy of the host galaxy of the repeating fast radio burst (FRB) source FRB 121102 and those of long gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe) has led to the suggestion that young magnetars born in GRBs and SLSNe could be the central engine of repeating FRBs. We test such a hypothesis by performing dedicated observations of the remnants of six GRBs with evidence of having a magnetar central engine using the Arecibo telescope and the Robert C. Byrd Green Bank Telescope (GBT). A total of $\sim 20$ hrs of observations of these sources did not detect any FRB from these remnants. Under the assumptions that all these GRBs left behind a long-lived magnetar and that the bursting rate of FRB 121102 is typical for a magnetar FRB engine, we estimate a non-detection probability of $8.9\times10^{-6}$. Even though these non-detections cannot exclude the young magnetar model of FRBs, we place constraints on the burst rate and luminosity function of FRBs from these GRB targets.

GRB 160821B is a short duration gamma-ray burst (GRB) detected and localized by the Neil Gehrels Swift Observatory in the outskirts of a spiral galaxy at z=0.1613, at a projected physical offset of 16 kpc from the galaxy's center. We present X-ray, optical/nIR and radio observations of its counterpart and model them with two distinct components of emission: a standard afterglow, arising from the interaction of the relativistic jet with the surrounding medium, and a kilonova, powered by the radioactive decay of the sub-relativistic ejecta. Broadband modeling of the afterglow data reveals a weak reverse shock propagating backward into the jet, and a likely jet-break at 3.5 d. This is consistent with a structured jet seen slightly off-axis while expanding into a low-density medium. Analysis of the kilonova properties suggests a rapid evolution toward red colors, similar to AT2017gfo, and a low nIR luminosity, possibly due to the presence of a long-lived neutron star. The global properties of the environment, the inferred low mass (M_ej < 0.006 Msun) and velocities (v > 0.05 c) of lanthanide-rich ejecta are consistent with a binary neutron star merger progenitor.

We search for host galaxy candidates of nearby fast radio bursts (FRBs), FRB 180729.J1316+55, FRB 171020, FRB 171213, FRB 180810.J1159+83, and FRB 180814.J0422+73 (the second repeating FRB). We compare the absolute magnitudes and the expected host dispersion measure $\rm DM_{host}$ of these candidates with that of the first repeating FRB, FRB 121102, as well as those of long gamma ray bursts (LGRBs) and superluminous supernovae (SLSNe), the proposed progenitor systems of FRB 121102. We find that while the FRB 121102 host is consistent with those of LGRBs and SLSNe, the nearby FRB host candidates, at least for FRB 180729.J1316+55, FRB 171020, and FRB180814.J0422+73, either have a smaller $\rm DM_{host}$ or are fainter than FRB121102 host, as well as the hosts of LGRBs and SLSNe. In order to avoid the uncertainty in estimating $\rm DM_{host}$ due to the line-of-sight effect, we propose a galaxy-group-based method to estimate the electron density in the inter-galactic regions, and hence, $\rm DM_{IGM}$. The result strengthens our conclusion. We conclude that the host galaxy of FRB 121102 is atypical, and LGRBs and SLSNe are likely not the progenitor systems of at least most nearby FRB sources. {The recently reported two FRB hosts differ from the host of FRB 121102 and also the host candidates suggested in this paper. This is consistent with the conclusion of our paper and suggests that the FRB hosts are very diverse. }

The LIGO and Virgo scientific collaboration (LVC) alerted two binary neutron star (BNS) merger candidates, S190425z and S190426c. Fermi-GBM observed 55.6\% (for S190425z) and 100\% (for S190426c) of the probability regions of both events at the respective merger times, but no gamma-ray burst (GRB) was detected in either case. The derived luminosity upper limits suggest that a short GRB similar to GRB 170817A would not be detectable for both cases due to their larger distances than GW170817. Assuming that the jet profile obtained from the GW170817/GRB 170817A is quasi-universal for all BNS-GRB associations, we derive that the viewing angles of S190425z and S190426c should be $> (0.11-0.41)$ and $> (0.09-0.39)$, respectively. Through Monte Carlo simulations, we show that with the GRB 170817A-like jet structure, all sky gamma-ray detectors, such as GBM and GECAM, are expected to detect $\sim 4.6\%$, $3.9\%$, $1.7\%$ and $6.6\%$, $5.7\%$, $2.8\%$ BNS mergers triggered by aLIGO, aLIGO A+ and ET, respectively. The joint detection fraction would be largely reduced for Swift-BAT, SVOM-ECLAIRS and Einstein Probe, whose sensitivities are better but fields of view are smaller.

The "Breakthrough Starshot" program is planning to send transrelativistic probes to travel to nearby stellar systems within decades. Since the probe velocity is designed to be a good fraction of the light speed, \citet{zha18} recently proposed that these transrelativistic probes can be used to study astronomical objects and to test special relativity. In this work, we further propose some methods to test special relativity and constrain photon mass using the Doppler effect with the images and spectral features of astronomical objects as observed in the transrelativistic probes. We introduce more general theories to set up the framework of testing special relativity, including a parametric general Doppler effect and Doppler effect with massive photon. We find that by comparing the spectra of a certain astronomical object, one can test Lorentz invariance and constrain photon mass. Besides, using imaging and spectrograph capabilities of transrelativistic probes, one can test time dilation and constrain photon mass. For a transrelativistic probe with velocity $v\sim0.2c$, aperture $D\sim3.5~{\rm cm}$ and spectral resolution $R\sim100$ (or $1000$), we find that the probe velocity uncertainty can be constrained to $\sigma_v\sim0.01c$ (or $0.001c$), and the time dilation factor uncertainty can be constrained to $\Delta\gamma=|\hat\gamma-\gamma|\lesssim0.01$ (or $0.001$), where $\hat\gamma$ is the time dilation factor and $\gamma$ is the Lorentz factor. Meanwhile, the photon mass limit is set to $m_\gamma\lesssim10^{-33}~{\rm g}$, which is slightly lower than the energy of the optical photon.

The IceCube neutrino event detected 360 s before the trigger of LIGO/Virgo S190728q is unlikely to be physically associated with the GW event because of an energy budget reason.

Based on the early-year observations from Neil Gehrels \emph{Swift} Observatory, \cite{liang07} performed a systematic analysis for the shallow decay component of gamma-ray bursts (GRBs) X-ray afterglow, in order to explore its physical origin. Here we revisit the analysis with an updated sample (with Swift/XRT GRBs between February 2004 and July 2017). We find that with a larger sample, 1) the distributions of the characteristic properties of the shallow decay phase (e.g. $t_{b}$ , $S_{X}$, $\Gamma_{X,1}$, and $\alpha_{X,1}$) still accords with normal or lognormal distribution; 2) $\Gamma_{X,1}$ and $\Gamma_{\gamma}$ still show no correlation, but the tentative correlations of durations, energy fluences, and isotropic energies between the gamma-ray and X-ray phases still exist; 3) for most GRBs, there is no significant spectral evolution between the shallow decay segment and its follow-up segment, and the latter is usually consistent with the external-shock models; 4) assuming that the central engine has a power-law luminosity release history as $L(t)=L_{0}(\frac{t}{t_{0}})^{-q}$, we find that the value $q$ is mainly distributed between -0.5 and 0.5, with an average value of 0.16$\pm$ 0.12; 5) the tentative correlation between $E_{\rm{iso},X}$ and $t'_{b}$ disappears, so that the global 3-parameter correlation ($E_{\rm{iso},X}-E'_{p}-t'_{b}$) becomes less significant; 6) the anti-correlation between $L_{X}$ and $t'_{b}$ and the three-parameter correlation ($E_{\rm{iso},\gamma}-L_{X}-t_{b}$) indeed exist with a high confidence level. Overall, our results are generally consistent with \cite{liang07}, confirming their suggestion that the shallow decay segment in most bursts is consistent with an external forward shock origin, probably due to a continuous energy injection from a long-lived central engine.

Based on the early-year observations from Neil Gehrels \emph{Swift} Observatory, \cite{liang07} performed a systematic analysis for the shallow decay component of gamma-ray bursts (GRBs) X-ray afterglow, in order to explore its physical origin. Here we revisit the analysis with an updated sample (with Swift/XRT GRBs between February 2004 and July 2017). We find that with a larger sample, 1) the distributions of the characteristic properties of the shallow decay phase (e.g. $t_{b}$ , $S_{X}$, $\Gamma_{X,1}$, and $\alpha_{X,1}$) still accords with normal or lognormal distribution; 2) $\Gamma_{X,1}$ and $\Gamma_{\gamma}$ still show no correlation, but the tentative correlations of durations, energy fluences, and isotropic energies between the gamma-ray and X-ray phases still exist; 3) for most GRBs, there is no significant spectral evolution between the shallow decay segment and its follow-up segment, and the latter is usually consistent with the external-shock models; 4) assuming that the central engine has a power-law luminosity release history as $L(t)=L_{0}(\frac{t}{t_{0}})^{-q}$, we find that the value $q$ is mainly distributed between -0.5 and 0.5, with an average value of 0.16$\pm$ 0.12; 5) the tentative correlation between $E_{\rm{iso},X}$ and $t'_{b}$ disappears, so that the global 3-parameter correlation ($E_{\rm{iso},X}-E'_{p}-t'_{b}$) becomes less significant; 6) the anti-correlation between $L_{X}$ and $t'_{b}$ and the three-parameter correlation ($E_{\rm{iso},\gamma}-L_{X}-t_{b}$) indeed exist with a high confidence level. Overall, our results are generally consistent with \cite{liang07}, confirming their suggestion that the shallow decay segment in most bursts is consistent with an external forward shock origin, probably due to a continuous energy injection from a long-lived central engine.

We present a narrow $\rm H\alpha$-band imaging survey of 357 low surface brightness galaxies (LSBGs) that are selected from the spring sky region of the 40% Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) HI Survey. All the $\rm H\alpha$ images are obtained from the 2.16 m telescope, operated by Xinglong Observatory of the National Astronomical Observatories, Chinese Academy of Sciences. We provide the $\rm H\alpha$ fluxes and derive the global star formation rates (SFRs) of LSBGs after the Galactic extinction, internal extinction, and [NII] contamination correction. Comparing to normal star-forming galaxies, LSBGs have a similar distribution in the HI surface density ($\rm \Sigma_{HI}$), but their SFRs and star formation surface density ($\rm \Sigma_{SFR}$) are much lower. Our results show that the gas-rich LSBGs selected from the ALFALFA survey obviously deviate from the Kennicutt-Schmidt law, in the relation between the star formation surface density ($\rm \Sigma_{SFR}$) and the gas surface density ($\rm \Sigma_{gas}$). However, they follow the extended Schmidt law well when taking the stellar mass of the galaxy into consideration.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

The first gravitational wave (GW) - gamma-ray burst (GRB) association, GW170817/GRB 170817A, had an offset in time, with the GRB trigger time delayed by $\sim$1.7 s with respect to the merger time of the GW signal. We generally discuss the astrophysical origin of the delay time, $\Delta t$, of GW-GRB associations within the context of compact binary coalescence (CBC) -- short GRB (sGRB) associations and GW burst -- long GRB (lGRB) associations. In general, the delay time should include three terms, the time to launch a clean (relativistic) jet, $\Delta t_{\rm jet}$; the time for the jet to break out from the surrounding medium, $\Delta t_{\rm bo}$; and the time for the jet to reach the energy dissipation and GRB emission site, $\Delta t_{\rm GRB}$. For CBC-sGRB associations, $\Delta t_{\rm jet}$ and $\Delta t_{\rm bo}$ are correlated, and the final delay can be from 10 ms to a few seconds. For GWB-lGRB associations, $\Delta t_{\rm jet}$ and $\Delta t_{\rm bo}$ are independent. The latter is at least $\sim$10 s, so that $\Delta t$ of these associations is at least this long. For certain jet launching mechanisms of lGRBs, $\Delta t$ can be minutes or even hours long due to the extended engine waiting time $\Delta t_{\rm wait}$ to launch a jet. We discuss the cases of GW170817/GRB 170817A and GW150914/GW150914-GBM within this theoretical framework and suggest that the delay times of future GW/GRB associations will shed light into the jet launching mechanisms of GRBs.

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

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

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

For non-spinning, charged (Reissner-Nordstr\"om) black holes, the particles with an opposite sign of charge with respect to that of the black hole will be pulled into the black hole by the extra electromagnetic force. Such a hole will be quickly neutralized so that there should not exist significantly charged, non-spinning black holes in the universe. The case of spinning, charged (Kerr-Newmann) black holes is more complicated. For a given initial position and initial velocity of the particle, an oppositely charged particle does not always more easily fall into the black hole than a neutral particle. The possible existence of a magnetosphere further complicate the picture. One therefore cannot straightforwardly conclude that a charged spinning black hole will be neutralized. In this paper, we make the first step to investigate the neutralization of Kerr-Newmann black holes without introducing a magnetosphere. We track the particle trajectories under the influence of the curved spacetime and the electromagnetic field carried by the spinning, charged black hole. A statistical method is used to investigate the neutralization problem. We find a universal dependence of the falling probability into the black hole on the charge of the test particle, with the oppositely charged particles having a higher probability of falling. We therefore conclude that charged, spinning black holes without a magnetosphere should be quickly neutralized, consistent with people's intuition. The neutralization problem of Kerr-Newmann black holes with a co-rotating force-free magnetosphere is subject to further studies.

H2O (22 GHz) and SiO masers (43, 86, 129 GHz) in the bipolar proto-planetary nebula OH231.8+4.2 were simultaneously monitored using the 21-m antennas of the Korean VLBI Network in 2009-2015. Both species exhibit periodic flux variations that correlate with the central star's optical light curve, with a phase delay of up to 0.15 for the maser flux variations with respect to the optical light curve. The flux densities of SiO v = 2, J = 1-0 and H2O masers decrease with time, implying that they may disappear in 10-20 years. However, there seems to have been a transient episode of intense H2O maser emission around 2010. We also found a systematic behaviour in the velocity profiles of these masers. The velocities of the H2O maser components appear to be remarkably constant, suggesting ballistic motion for the bipolar outflow in this nebula. On the other hand, those of the SiO maser clumps show a systematic radial acceleration of the individual clumps, converging to the outflow velocity of the H2O maser clumps. Measuring the full widths at zero power of the detected lines, we estimated the expansion velocities of the compact bipolar outflow traced by H2O maser and SiO thermal line, and discussed the possibility of the expanding SiO maser region in the equatorial direction. All of our analyses support that the central host star of OH231.8 is close to the tip of the AGB phase, and that the mass-loss rate recently started to decrease because of incipient post-AGB evolution.

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

Very-high-energy (VHE; $\geq 10$ GeV) photons are expected from the nearest and brightest Gamma-ray bursts (GRBs). VHE photons, at energies higher than 300 GeV, were recently reported by the MAGIC collaboration for this burst. Immediately, GRB 190114C was followed up by a massive observational campaign covering a large fraction of the electromagnetic spectrum. In this paper, we obtain the LAT light curve of GRB 190114C and show that it exhibits similar features to other bright LAT-detected bursts; the first high-energy photon ($\geq$ 100 MeV) is delayed with the onset of the prompt phase and the flux light curve exhibits a long-lived emission (lasting much longer than the prompt phase) and a short-lasting bright peak (located at the beginning of long-lived emission). Analyzing the multi-wavelength observations, we show that the short-lasting LAT and GBM bright peaks are consistent with the synchrotron self-Compton reverse-shock model and the long-lived observations with the standard synchrotron forward-shock model that evolves from a stratified stellar-wind like medium to a uniform ISM-like medium. Given the best-fit values, a bright optical flash produced by synchrotron reverse-shock emission is expected. From our analysis we infer that the high-energy photons are produced in the deceleration phase of the outflow and some additional processes to synchrotron in the forward shocks should be considered to properly describe the LAT photons with energies beyond the synchrotron limit. Moreover, we claim that an outflow endowed with magnetic fields could describe the polarization and properties exhibited in the light curve of GRB 190114C.

Very recently \citet{XueYQ2019} reported an important detection of the X-ray transient, CDF-S XT2, whose light curve is analogous to X-ray plateau features of gamma-ray burst afterglows. They suggested that this transient is powered by a remnant stable magnetar from a binary neutron star merger since several pieces of evidence (host galaxy, location, and event rate) all point toward such an assumption. In this paper, we revisit this scenario and confirm that this X-ray emission can be well explained by the internal gradual magnetic dissipation process in an ultra-relativistic wind of the newborn magnetar. We show that both the light curve and spectral evolution of CDF-S XT2 can be well fitted by such a model. Furthermore, we can probe some key properties of the central magnetar, such as its initial spin period, surface magnetic field strength and wind saturation Lorentz factor.

We present the analysis underpinning the measurement of cosmological parameters from 207 spectroscopically classified type Ia supernovae (SNe Ia) from the first three years of the Dark Energy Survey Supernova Program (DES-SN), spanning a redshift range of 0.017<$z$<0.849. We combine the DES-SN sample with an external sample of 122 low-redshift ($z$<0.1) SNe Ia, resulting in a "DES-SN3YR" sample of 329 SNe Ia. Our cosmological analyses are blinded: after combining our DES-SN3YR distances with constraints from the Cosmic Microwave Background (CMB; Planck Collaboration 2016), our uncertainties in the measurement of the dark energy equation-of-state parameter, $w$, are .042 (stat) and .059 (stat+syst) at 68% confidence. We provide a detailed systematic uncertainty budget, which has nearly equal contributions from photometric calibration, astrophysical bias corrections, and instrumental bias corrections. We also include several new sources of systematic uncertainty. While our sample is <1/3 the size of the Pantheon sample, our constraints on $w$ are only larger by 1.4$\times$, showing the impact of the DES SN Ia light curve quality. We find that the traditional stretch and color standardization parameters of the DES SNe Ia are in agreement with earlier SN Ia samples such as Pan-STARRS1 and the Supernova Legacy Survey. However, we find smaller intrinsic scatter about the Hubble diagram (0.077 mag). Interestingly, we find no evidence for a Hubble residual step ( 0.007 $\pm$ 0.018 mag) as a function of host galaxy mass for the DES subset, in 2.4$\sigma$ tension with previous measurements. We also present novel validation methods of our sample using simulated SNe Ia inserted in DECam images and using large catalog-level simulations to test for biases in our analysis pipelines.

We present an improved measurement of the Hubble constant (H_0) using the 'inverse distance ladder' method, which adds the information from 207 Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) at redshift 0.018 < z < 0.85 to existing distance measurements of 122 low redshift (z < 0.07) SNe Ia (Low-z) and measurements of Baryon Acoustic Oscillations (BAOs). Whereas traditional measurements of H_0 with SNe Ia use a distance ladder of parallax and Cepheid variable stars, the inverse distance ladder relies on absolute distance measurements from the BAOs to calibrate the intrinsic magnitude of the SNe Ia. We find H_0 = 67.8 +/- 1.3 km s-1 Mpc-1 (statistical and systematic uncertainties, 68% confidence). Our measurement makes minimal assumptions about the underlying cosmological model, and our analysis was blinded to reduce confirmation bias. We examine possible systematic uncertainties and all are below the statistical uncertainties. Our H_0 value is consistent with estimates derived from the Cosmic Microwave Background assuming a LCDM universe (Planck Collaboration et al. 2018).

We perform a series of relativistic magnetohydrodynamics simulations to investigate how a hot magnetic jet propagates within the dynamical ejecta of a binary neutron star merger, with the focus on how the jet structure depends on the delay time of jet launching with respect to the merger time, $\Delta t_{\rm jet}$. We find that regardless of the jet launching delay time, a structured jet with an angle-dependent luminosity and Lorentz factor is always formed after the jet breaks out the ejecta. On the other hand, the jet launching delay time has an impact on the jet structure. If the jet launching delay time is relatively long, e.g., $\ge 0.5$ s, the line-of-sight material has a dominant contribution from the cocoon. On the other hand, for a relatively short jet launching delay time, the jet penetrates through the ejecta early on and develops an angular structure afterward. The line-of-sight ejecta is dominated by the structured jet itself. We discuss the case of GW170817/GRB 170817A within the framework of both long and short $\Delta t_{\rm jet}$. Future more observations of GW/GRB associations can help to differentiate between these two scenarios.

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.

With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

The origin of fast radio bursts (FRBs) is still unknown. Multi-wavelength observations during or shortly after the FRB phase would be essential to identify the counterpart of an FRB and to constrain its progenitor and environment. In this work, we investigate the brightness of the "fast optical bursts" (FOBs) associated with FRBs and the prospects of detecting them. We investigate several inverse Compton (IC) scattering processes that might produce an FOB, including both the one-zone and two-zone models. We also investigate the extension of the same mechanism of FRB emission to the optical band. We find that a detectable FOB with the current and forthcoming telescopes is possible under the IC scenarios with very special conditions. In particular, the FRB environment would need to invoke a neutron star with an extremely strong magnetic field and an extremely fast spin, or an extremely young supernova remnant (SNR) surrounding the FRB source. Furthermore, most electrons in the source are also required to have a fine-tuned energy distribution such that most of the IC energy is channeled in the optical band. We conclude that the prospect of detecting FOBs associated with FRBs is low. On the other hand, if FOBs are detected from a small fraction of FRBs, these FOBs would reveal extreme physical conditions in the FRB environments.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

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

The newly discovered second repeating fast radio burst (FRB) source, FRB 180814.J0422+73, was reported to exhibit a time-frequency downward drifting pattern, which is also seen in the first repeater FRB 121102. We propose a generic geometrical model to account for the observed downward drifting of sub-pulse frequency, within the framework of coherent curvature radiation by bunches of electrons or positrons in the magnetosphere of a neutron star. A sudden trigger event excites these coherent bunches of charged particles, which stream outwards along open field lines. As the field lines sweep across the line of sight, the bunches seen later have traveled farther into the less curved part of the magnetic field lines, thus emitting at lower frequencies. We use this model to explain the time-frequency downward drifting in two FRB generation scenarios, the transient pulsar-like sparking from the inner gap region of a slowly rotating neutron star, and the externally-triggered magnetosphere reconfiguration known as the "cosmic comb".

The merging process of binary neutron stars is the source of the quadrupole gravitational-wave (GW) radiation. Its remnant can be either a supra-massive or a stable NS. Such evidence of magnetar signature has been supported indirectly by observed X-ray plateau of some gamma-ray bursts (GRBs) afterglow. Recently, Xue et al. (2019) discovered an X-ray transient CDF-S XT2 that is claimed to be powered by a stable magnetar from the merger of double NS. In this paper, we revisit the X-ray emission of CDF-S XT2 and find that it is more consistent with a supra-massive magnetar central engine, surviving thousands of seconds to collapse black hole. We present the comparisons of the X-ray plateau luminosity, break time, and the parameters of magnetar between CDF-S XT2 and other short GRBs with internal plateau samples. By adapting the collapse time to constrain the equation of state (EOS), we find that three EOSs (GM1, DD2, and DDME2) are consistent with the observational data. On the other hand, if the most released rotation energy of magnetar is dominated by GW radiation, we also constrain the upper limit of ellipticity of NS for given EOS, and it is in a range in $[0.89-1.8]\times 10^{-3}$. Its GW signal cannot be detected by aLIGO or even for more sensitive Einstein Telescope in the future.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

The Breakthrough Starshot project aims to send centimeter-sized, gram-scale "StarChip" probes to Alpha Centauri at a speed of $\sim0.2c$. On the other hand, Zhang \& Li recently proposed that trans-relativistic cameras may be sent to any direction to study astronomical objects and test special relativity. To conduct such "relativistic astronomy", one needs to solve the motion of the probe in flight. We solve the motion of the probe (including the moving direction and the velocity) by comparing the positions of three or more point sources observed in the Earth rest frame and in the probe's comoving frame. When the positions of enough point sources are taken into account, the motion of the probe can be solved with an error which is even smaller than the diffraction limit of the transrelativistic camera. After solving the motion, when the measurement of the position of an additional point source is introduced, one can use the data to test the light aberration effect in special relativity. The upper limit of the photon mass can be placed from the deviation of aberration to slightly lower than the energy of the photon, e.g. $\sim 1 \rm {eV}$.

Enabled by the Fermi Large Area Telescope, we now know young and recycled pulsars fill the gamma-ray sky, and we are beginning to understand their emission mechanism and their distribution throughout the Galaxy. However, key questions remain: Is there a large population of pulsars near the Galactic center? Why do the most energetic pulsars shine so brightly in MeV gamma rays but not always at GeV energies? What is the source and nature of the pair plasma in pulsar magnetospheres, and what role does the polar cap accelerator play? Addressing these questions calls for a sensitive, wide-field MeV telescope, which can detect the population of MeV-peaked pulsars hinted at by Fermi and hard X-ray telescopes and characterize their spectral shape and polarization.

The conventional view of stellar-mass binary black hole (sBBH) mergers is that there should not be enough matter present to produce a detectable electromagnetic transient. However, there ARE a number of mechanisms for producing such a counterpart through accretion of matter from various reserves, charged black holes, or interactions with magnetic or exotic fields. After reviewing these mechanisms, we describe what we can learn from multi-messenger observations of sBBH mergers in the areas of stellar evolution and compact binary formation, cosmological measurements, searches for charged black holes, and tests of general relativity and fundamental physics. We end with some remarks about the need to support both observing capabilities and modeling in the next decade and beyond.

We report trigonometric parallax and proper motion measurements of 6.7-GHz CH3OH and 22-GHz H2O masers in eight high-mass star-forming regions (HMSFRs) based on VLBA observations as part of the BeSSeL Survey. The distances of these HMSFRs combined with their Galactic coordinates, radial velocities, and proper motions, allow us to assign them to a segment of the Perseus arm with ~< 70 deg. These HMSFRs are clustered in Galactic longitude from ~30 deg to ~50, neighboring a dirth of such sources between longitudes ~50 deg to ~90 deg.

The second repeating fast radio burst source, FRB 180814.J0422+73, was detected recently by the CHIME collaboration. We use the ten-year Fermi Large Area Telescope (LAT) archival data to place a flux upper limit in the energy range of 100 MeV$-$10 GeV at the position of the source, which is $\sim 10^{-11}$ erg cm$^{-2}$ s$^{-1}$ for a six-month time bin on average, and $2.35\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ for the entire ten-year time span. For the maximum redshift of $z=0.11$, the ten-year upper limit of luminosity is $7.32\times 10^{43}$ erg s$^{-1}$. We utilize these upper limits to constrain the FRB progenitor and central engine. For the rotation-powered young magnetar model, the upper limits can pose constraints on the allowed parameter space for the initial rotational period and surface magnetic field of the magnetar. We also place significant constraints on the kinetic energy of a relativistic external shock wave, ruling out the possibility that there existed a gamma-ray burst (GRB) beaming towards earth during the past ten years as the progenitor of the repeater. The case of an off-beam GRB is also constrained if the viewing angle is not much greater than the jet opening angle. All these constraints are more stringent if FRB 180814.J0422+73 is at a closer distance.

The second repeating fast radio burst source, FRB 180814.J0422+73, was detected recently by the CHIME collaboration. We use the ten-year Fermi Large Area Telescope (LAT) archival data to place a flux upper limit in the energy range of 100 MeV$-$10 GeV at the position of the source, which is $\sim 10^{-11}$ erg cm$^{-2}$ s$^{-1}$ for a six-month time bin on average, and $2.35\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ for the entire ten-year time span. For the maximum redshift of $z=0.11$, the ten-year upper limit of luminosity is $7.32\times 10^{43}$ erg s$^{-1}$. We utilize these upper limits to constrain the FRB progenitor and central engine. For the rotation-powered young magnetar model, the upper limits can pose constraints on the allowed parameter space for the initial rotational period and surface magnetic field of the magnetar. We also place significant constraints on the kinetic energy of a relativistic external shock wave, ruling out the possibility that there existed a gamma-ray burst (GRB) beaming towards earth during the past ten years as the progenitor of the repeater. The case of an off-beam GRB is also constrained if the viewing angle is not much greater than the jet opening angle. All these constraints are more stringent if FRB 180814.J0422+73 is at a closer distance.

We compare the parallaxes of stars from VLBI astrometry in the literature to those in the Gaia DR2 catalog. Our full sample contains young stellar objects, evolved AGB stars, pulsars and other radio stars. Excluding AGB stars, which show significant discrepancies between Gaia and VLBI parallaxes, and stars in binary systems, we obtain an average, systematic, parallax offset of $-75 \pm 29~\mu$as for Gaia DR2, consistent with their estimate of a parallax zero-point between $-100$ and 0 $\mu$as.

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

The first multimessenger observation of a neutron star merger was independently detected in gamma-rays by Fermi-GBM and INTEGRAL SPI-ACS and gravitational waves by Advanced LIGO and Advanced Virgo. Gravitational waves are emitted from systems with accelerating quadrupole moments, and detectable sources are expected to be compact objects. Nearly all distant astrophysical gamma-ray sources are compact objects. Therefore, serendipitous observations of these two messengers will continue to uncover the sources of gravitational waves and gamma-rays, and enable multimessenger science across the Astro2020 thematic areas. This requires upgrades to the ground-based gravitational wave network and ~keV-MeV gamma-ray coverage for observations of neutron star mergers, and broadband coverage in both gravitational waves and gamma-rays to monitor other expected joint sources.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Neutron star mergers, referring to both binary neutron star and neutron star black hole mergers, are the canonical multimessenger events. They have been detected across the electromagnetic spectrum, have recently been detected in gravitational waves, and are likely to produce neutrinos over several decades in energy. The non-thermal prompt and afterglow emission of short gamma-ray bursts and the quasi-thermal emission from the radioactively powered kilonovae provide distinct insights into the physics of neutron star mergers. When combined with direct information on coalescence from gravitational waves and neutrinos these sources may become the best understood astrophysical transients. Multimessenger observations of these cataclysmic events will determine sources of gravitational waves and astrophysical neutrinos, enable precision cosmology, and unique tests of fundamental physics, the origin of heavy elements, the behavior of relativistic jets, and the equation of state of supranuclear matter. In this white paper we present a summary of the science discoveries possible with multimessenger observations of neutron star mergers and provide recommendations to enable them in the new era of time-domain, multimessenger astronomy.

Magnetars are a class of highly magnetized, slowly rotating neutron stars, only a small fraction of which exhibit radio emission. We propose that the coherent radio curvature emission is generated by net charge fluctuations from a twist-current-carrying bundle (the j-bundle) in the scenario of magnetar-quake. Two-photon pair production is triggered, which requires a threshold voltage not too much higher than $10^9$ V in the current-carrying bundle, and which can be regarded as the `open field lines' of a magnetar. Continued untwisting of the magnetosphere maintains change fluctuations, and hence coherent radio emission, in the progressively shrinking j-bundle, which lasts for years until the radio beam is too small to be detected. The modeled peak flux of radio emission and the flat spectrum are generally consistent with the observations. We show that this time-dependent, conal-beam, radiative model can interpret the variable radio pulsation behaviors and the evolution of the X-ray hot spot of the radio transient magnetar XTE J1810$-$197 and the high-$B$ pulsar/AXP PSR J1622$-$4950. Radio emission with luminosity of $\lesssim10^{31}\,{\rm erg\,s^{-1}}$ and high-frequency oscillations are expected to be detected for a magnetar after an X-ray outburst. Differences of radio emission between magnetars and ordinary pulsars are discussed.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

We present 22\,901 OB spectra of 16\,032 stars identified from LAMOST DR5 dataset. A larger sample of OB candidates are firstly selected from the distributions in the spectral line indices space. Then all 22\,901 OB spectra are identified by manual inspection. Based on a sub-sample validation, we find that the completeness of the OB spectra reaches about $89\pm22$\% for the stars with spectral type earlier than B7, while around $57\pm16$\% B8--B9 stars are identified. The smaller completeness for late B stars is lead to the difficulty to discriminate them from A0--A1 type stars. The sub-classes of the OB samples are determined using the software package MKCLASS. With a careful validation using 646 sub-samples, we find that MKCLASS can give fairly reliable sub-types and luminosity class for most of the OB stars. The uncertainty of the spectral sub-type is around 1 sub-type and the uncertainty of the luminosity class is around 1 level. However, about 40\% of the OB stars are failed to be assigned to any class by MKCLASS and a few spectra are significantly misclassified by MKCLASS. This is likely because that the template spectra of MKCLASS are selected from nearby stars in the solar neighborhood, while the OB stars in this work are mostly located in the outer disk and may have lower metallicity. The rotation of the OB stars may also be responsible for the mis-classifications. Moreover, we find that the spectral and luminosity classes of the OB stars located in the Galactic latitude larger than 20$^\circ$ are substantially different with those located in latitude smaller than 20$^\circ$, which may either due to the observational selection effect or hint a different origin of the high Galactic latitude OB stars.

The origin of fast radio bursts (FRBs) is still a mystery. One model proposed to interpret the only known repeating object, FRB 121102, is that the radio emission is generated from asteroids colliding with a highly magnetized neutron star (NS). With N-body simulations, we model a debris disc around a central star with an eccentric orbit intruding NS. As the NS approaches the first periastron passage, most of the comets are scattered away rather than being accreted by the NS. To match the observed FRB rate, the debris belt would have to be at least three orders of magnitude more dense than the Kuiper belt. We also consider the rate of collisions on to the central object but find that the density of the debris belt must be at least four orders of magnitude more dense than the Kuiper belt. These discrepancies in the density arise even if (1) one introduces a Kuiper-belt like comet belt rather than an asteroid belt and assume that comet impacts can also make FRBs; (2) the NS moves 2 orders of magnitude slower than their normal proper-motion velocity due to supernova kicks; and (3) the NS orbit is coplanar to the debris belt, which provides the highest rate of collisions.

We report the distances of molecular clouds at high Galactic latitudes (|b|>10$^\circ$) derived from parallax and G band extinction (A$_{\rm G}$) measurements in the second Gaia data release, Gaia DR2. Aided by Bayesian analyses, we determined distances by identifying the breakpoint in the extinction A$_{\rm G}$ towards molecular clouds and using the extinction A$_{\rm G}$ of Gaia stars around molecular clouds to confirm the breakpoint. We use nearby star-forming regions, such as Orion, Taurus, Cepheus, and Perseus, whose distances are well-known to examine the reliability of our method. By comparing with previous results, we found that the molecular cloud distances derived from this method are reliable. The systematic error in the distances is approximately 5%. In total, 52 molecular clouds have their distances well determined, most of which are at high Galactic latitudes, and we provide reliable distances for 13 molecular clouds for the first time.

The detection of the high-energy ($\sim290$ TeV) neutrino coincident with the flaring blazar TXS 0506+056, the first and only $3\sigma$ neutrino-source association to date, provides new, multimessenger tests of the weak equivalence principle (WEP) and Lorentz invariance. Assuming that the flight time difference between the TeV neutrino and gamma-ray photons from the blazar flare is mainly caused by the gravitational potential of the Laniakea supercluster of galaxies, we show that the deviation from the WEP for neutrinos and photons is conservatively constrained to have an accuracy of $10^{-6}-10^{-7}$, which is 3--4 orders of magnitude better than previous results placed by MeV neutrinos from supernova 1987A. In addition, we demonstrate that the association of the TeV neutrino with the blazar flare sets limits on the energy scales of quantum gravity for both linear and quadratic violations of Lorentz invariance (LIV) to $E_{\rm QG, 1}>3.2\times10^{15}-3.7\times10^{16}$ GeV and $E_{\rm QG, 2}>4.0\times10^{10}-1.4\times10^{11}$ GeV. These improve previous limits on both linear and quadratic LIV energy scales in neutrino propagation by 5--7 orders of magnitude.

The detection of the high-energy ($\sim290$ TeV) neutrino coincident with the flaring blazar TXS 0506+056, the first and only $3\sigma$ neutrino-source association to date, provides new, multimessenger tests of the weak equivalence principle (WEP) and Lorentz invariance. Assuming that the flight time difference between the TeV neutrino and gamma-ray photons from the blazar flare is mainly caused by the gravitational potential of the Laniakea supercluster of galaxies, we show that the deviation from the WEP for neutrinos and photons is conservatively constrained to have an accuracy of $10^{-6}-10^{-7}$, which is 3--4 orders of magnitude better than previous results placed by MeV neutrinos from supernova 1987A. In addition, we demonstrate that the association of the TeV neutrino with the blazar flare sets limits on the energy scales of quantum gravity for both linear and quadratic violations of Lorentz invariance (LIV) to $E_{\rm QG, 1}>3.2\times10^{15}-3.7\times10^{16}$ GeV and $E_{\rm QG, 2}>4.0\times10^{10}-1.4\times10^{11}$ GeV. These improve previous limits on both linear and quadratic LIV energy scales in neutrino propagation by 5--7 orders of magnitude.

As part of the BeSSeL Survey, we report trigonometric parallaxes and proper motions of molecular maser sources associated with 13 distant high mass star forming regions in the Sagittarius spiral arm of the Milky Way. In particular, we obtain improved parallax distance estimates for three well studied regions: 1.9 +0.1/-0.1 kpc for M17, 5.3 +1.3/-0.9 kpc for W51, and 7.9 +0.9/-0.7 kpc for GAL 045.5+00.0. Peculiar motions for all but one source are less than 20 km/s. We fit a log-periodic spiral to the locations and estimate an average pitch angle of 7.2+-1.9 deg. We find that the section of the arm beyond the tangent point in the first quadrant of the Milky Way appears 15 pc below the IAU-defined Galactic plane.

As shown by broad-band observations, pulsar wind nebulae (PWNe) are characterized by a broken power-law spectrum of synchrotron emission. Based on the modern magnetohydrodynamic (MHD) turbulence theories, we investigate the re-acceleration of electrons in the PWN through the adiabatic stochastic acceleration (ASA), which arises from fundamental dynamics of MHD turbulence. The ASA acts to flatten the injected energy spectrum of electrons at low energies, while synchrotron cooling results in a steep spectrum of electrons at high energies. Their dominance in different energy ranges leads to a flat radio spectrum ($F_\nu$) and a steep X-ray spectrum. Our analytical spectral shapes generally agree well with the observed synchrotron spectra of radio- and X-ray-bright PWNe. The spectral break corresponding to the balance between the ASA and synchrotron losses provides a constraint on the acceleration timescale of the ASA and the magnetic field strength in the PWN.

The spectrum of gamma-ray burst (GRB) afterglows can be studied with color indices. Here we present a large comprehensive catalogue of 70 GRBs with multi-wavelength optical transient data on which we perform a systematic study to find the temporal evolution of color indices. We categorize them into two samples based on how well the color indices are evaluated. The Golden sample includes 25 bursts mostly observed by GROND, and the Silver sample includes 45 bursts observed by other telescopes. For the Golden sample, we find that 95\% of the color indices do not vary over time. However, the color indices do vary during short periods in most bursts. The observed variations are consistent with effects of (i) the cooling frequency crossing the studied energy bands in a wind medium (43\%) and in a constant density medium (30\%), (ii) early dust extinction (12\%), (iii) transition from reverse shock to forward shock emission (5\%), or (iv) an emergent supernova emission (10\%). We also study the evolutionary properties of the mean color indices for different emission episodes. We find that 86\% of the color indices in the 70 bursts show constancy between consecutive ones. The color index variations occur mainly during the late GRB-SN bump, the flare and early reversed-shock emission components. We further perform a statistical analysis of various observational properties and model parameters (spectral index $\beta_{o}^{\rm CI}$, electron spectral indices $p^{\rm CI}$, etc.) using color indices. Overall, we conclude that $\sim$ 90\% of colors are constant in time and can be accounted for by the simplest external forward shock model, while the varying color indices call for more detailed modeling.

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a).

Multi-messenger observations of GW170817 have not conclusively established whether the merger remnant is a black hole (BH) or a neutron star (NS). We show that a long-lived magnetized NS with a poloidal field $B\approx 10^{12}$G is fully consistent with the electromagnetic dataset, when spin down losses are dominated by gravitational wave (GW) emission. The required ellipticity $\epsilon\gtrsim 10^{-5}$ can result from a toroidal magnetic field component much stronger than the poloidal component, a configuration expected from a NS newly formed from a merger. Abrupt magnetic dissipation of the toroidal component can lead to the appearance of X-ray flares, analogous to the one observed in gamma-ray burst (GRB) afterglows. In the X-ray afterglow of GW170817 we identify a low-significance ($\gtrsim 3\sigma$) temporal feature at 155 d, consistent with a sudden reactivation of the central NS. Energy injection from the NS spin down into the relativistic shock is negligible, and the underlying continuum is fully accounted for by a structured jet seen off-axis. Whereas radio and optical observations probe the interaction of this jet with the surrounding medium, observations at X-ray wavelengths, performed with adequate sampling, open a privileged window on to the merger remnant.

Recent wide-field surveys discovered new types of peculiar optical transients that showed diverse behaviors of the evolution of photospheric properties. We develop a general theory of homologous explosions with constant opacity, paying special attention on the evolution of the photospheric radius $R_{\rm ph}$. We find that regardless of the density distribution profile, $R_{\rm ph}$ always increases early on and decreases at late times. This result does not depend on the radiation and cooling processes inside the ejecta.The general rising/falling behavior of $R_{\rm ph}$ can be used to quickly diagnose whether the source originates from a supernova-like explosion. The shape of the $R_{\rm ph}$ evolution curve depends on the density profile, so the observations may be used to directly diagnose the density profile as well as the temperature profile of the ejecta. All the well-monitored supernovae show such a $R_{\rm ph}$ rising/falling behavior, which is consistent with our theory. The recently discovered peculiar transient AT2018cow showed a continuous decay of $R_{\rm ph}$, which is disfavored to be of a supernova-like explosion origin. Our result therefore supports the interpretation of this transient as a tidal disruption event.

The radio luminosity functions (RLFs) of active galactic nuclei (AGNs) are traditionally measured based on total emission, which doesn't reflect the current activity of the central black hole. The increasing interest in compact radio cores of AGNs requires determination of the RLF based on core emission (i.e., core RLF). In this work we have established a large sample (totaling 1207) of radio-loud AGNs, mainly consisting of radio galaxies (RGs) and steep-spectrum radio quasars (SSRQs). Based on the sample, we explore the relationship between core luminosity ($L_c$) and total luminosity ($L_t$) via a powerful statistical tool called "Copula". The conditional probability distribution $p(\log L_{c} \mid \log L_{t})$ is obtained. We derive the core RLF as a convolution of $p(\log L_{c} \mid \log L_{t})$ with the total RLF which was determined by previous work. We relate the separate RG and SSRQ core RLFs via a relativistic beaming model and find that SSRQs have an average Lorentz factor of $\gamma=9.84_{-2.50}^{+3.61}$, and that most are seen within $8^{\circ} \lesssim \theta \lesssim 45^{\circ}$ of the jet axis. Compared with the total RLF which is mainly contributed by extended emission, the core RLF shows a very weak luminosity-dependent evolution, with the number density peaking around $z\thicksim 0.8$ for all luminosities. Differences between core and total RLFs can be explained in a framework involving a combination of density and luminosity evolutions where the cores have significantly weaker luminosity evolution than the extended emission.

The prompt emission phase of gamma-ray bursts (GRBs) exhibits two distinct patterns of the peak-energy ($E_p$) evolution, i.e., time-resolved spectral analyses of $\nu F_{\nu}$ spectra of broad pulses reveal (1) "hard-to-soft" and (2) "flux-tracking" patterns of $E_p$ evolution in time, the physical origin of which still remains not well understood. We show here that these two patterns can be successfully reproduced within a simple physical model invoking synchrotron radiation in a bulk-accelerating emission region. We show further that the evolution patterns of the peak energy have, in fact, direct connections to the existence of two different (positive or negative) types of spectral lags, seen in the broad pulses. In particular, we predict that (1) only the positive type of spectral lags is possible for the hard-to-soft evolution of the peak energy, (2) both the positive and negative type of spectral lags can occur in the case of flux-tracking pattern of the peak energy, (3) for the flux-tracking pattern, the peak location of the flux light curve slightly lags behind the peak of the $E_p$ evolution with time if the spectral lags are positive, and (4) in the case of flux-tracking pattern, double-peaked broad pulses can appear in the light curves, the shape of which is energy-dependent.

We estimate the upper limit redshifts of known FRBs using the dispersion measure (DM) - redshift ($z$) relation and derive the upper limit peak luminosity $L_p$ and energy $E$ of FRBs within the observational band. The average $z$ upper limits range from 0.17 to 3.10, the average $L_p$ upper limits range from $1.24 \times 10^{42} \rm erg \ s^{-1}$ to $7.80 \times 10^{44} \rm erg \ s^{-1}$, and the average $E$ upper limits range from $6.91 \times 10^{39}$ erg to $1.94 \times 10^{42}$ erg. FRB 160102 with DM $=2596.1 \pm 0.3 \ {\rm pc \ cm^{-3}}$ likely has a redshift greater than 3. Assuming that its intrinsic DM contribution from the host and FRB source is ${\rm DM_{host}+DM_{scr}}\sim 100 \ {\rm pc \ cm^{-3}}$, such an FRB can be detected up to $z \sim 3.6$ by Parkes and by FAST under ideal conditions up to $z \sim 10.4$. Assuming that there exist FRBs detectable at $z\sim 15$ by sensitive telescopes such as FAST, the upper limit DM for FRB searches may be set to $\sim 9000 \ {\rm pc \ cm^{-3}}$. For single-dish telescopes, those with a larger aperture tend to detect more FRBs than those with a smaller aperture if the FRB luminosity function index $\alpha_{\rm L}$ is steeper than 2, and vice versa. In any case, large-aperture telescopes such as FAST are more capable of detecting high-$z$ FRBs, even though most of FRBs detected by them are still from relatively low redshifts.

Fast Radio bursts (FRBs), bright transients with millisecond durations at $\sim$ GHz and typical redshifts probably $>0.8$, are likely to be gravitationally lensed by intervening galaxies. Since the time delay between images of strongly lensed FRB can be measured to extremely high precision because of the large ratio $\sim10^9$ between the typical galaxy-lensing delay time $\sim\mathcal{O}$(10 days) and the width of bursts $\sim\mathcal{O}$(ms), we propose strongly lensed FRBs as precision probes of the universe. We show that, within the flat $\Lambda$CDM model, the Hubble constant $H_0$ can be constrained with a $\sim0.91\%$ uncertainty from 10 such systems probably observed with the Square Kilometer Array (SKA) in $<$ 30 years. More importantly, the cosmic curvature can be model-independently constrained to a precision of $\sim0.076$. This constraint can directly test the validity of the cosmological principle and break the intractable degeneracy between the cosmic curvature and dark energy.

The extremely high brightness temperatures of pulsars and fast radio bursts (FRBs) require their radiation mechanisms to be coherent. Coherent curvature radiation by bunches has been long discussed as the mechanism for radio pulsars and recently for FRBs. Assuming that bunches are already generated in pulsar magnetospheres, we calculate the spectrum of coherent curvature radiation under a three-dimensional magnetic field geometry. Different from previous works assuming parallel trajectories and a mono-energetic energy distribution of electrons, we consider a bunch characterized by its length, curvature radius of the trajectory family, bunch opening angle, and electron energy distribution. We find that the curvature radiation spectra of the bunches are characterized by a multi-segment broken power law, with the break frequencies depending on bunch properties and trajectory configuration. We also emphasize that in a pulsar magnetosphere only the fluctuation of net charges with respect to the background (Goldreich-Julian) outflow can make a contribution to coherent radiation. We apply this model to constrain the observed spectra of pulsars and FRBs. For a typical pulsar ($B_p=10^{12}~{\rm G}$, $P=0.1~{\rm s}$), a small fluctuation of the net charge $\delta n_{\rm GJ}\sim 0.1 n_{\rm GJ}$ can provide the observable flux. For FRBs, the fluctuating net charge may be larger due to its abrupt nature. For $\delta n_{\rm GJ}\sim n_{\rm GJ}$, a neutron star with a strong magnetic field and fast rotation is required to power an FRB in the spindown-powered model. The requirement is less stringent in the cosmic comb model thanks to the larger cross section and compressed charge density of the bunch made by the external astrophysical stream that combs the magnetosphere.

In this paper, we revisit the scenario that an internal gradual magnetic dissipation takes place within the wind from a newborn millisecond magnetar can be responsible for gamma-ray burst production. We show that a combination of two emission components in this model, i.e., the photospheric emission from the wind and the synchrotron radiation within the magnetic reconnection region, can give a reasonable fit to the observed spectrum of the prompt emission phase of GRB 160804A. We obtain the physical parameters through a Monte Carlo procedure and deduce the initial spin period and magnetic field of the central magnetar. Furthermore, the independent afterglow fitting analysis gives a consistent result, adding great credibility to this scenario. In addition, we predict a subclass of GRBs called bursts from such a Magnetar wind Internal Gradual MAgnetic Dissipation (abbreviated as "MIGMAD bursts") that have several distinctive properties.

Massive neutron star (NS) with lifetimes of at least several seconds are expected to be the result of a sizable fraction of NS mergers. We study their formation using a large set of numerical relativity simulations. We show that they are initially endowed with angular momentum that significantly exceeds the mass-shedding limit for rigidly-rotating equilibria. We find that gravitational-wave (GW) emission is not able to remove this excess angular momentum within the time over which solid body rotation should be achieved. Instead, we argue that the excess angular momentum could be carried away by massive winds. Long-lived merger remnants are also formed with larger gravitational masses than those of rigidly-rotating NSs having the same number of baryons. The excess mass is likely radiated in the form of neutrinos. The evolution of long-lived remnants on the viscous timescale is thus determined by the interplay of finite-temperature effects, mass ejection, and neutrinos with potentially dramatic consequences for the remnants' properties and stability. We also provide an empirical fit for the spin of the remnant at the end of its viscous evolution as a function of its final mass, and we discuss the implications for the magnetar model of short gamma-ray bursts (SGRBs). Finally, we investigate the possible electromagnetic signatures associated with the viscous ejecta. Massive outflows possibly resulting from the formation of long-lived remnants would power unusually bright, blue kilonova counterparts to GW events and SGRBs whose detection would provide smoking gun evidence for the formation of long-lived remnants.

We report the discovery of a new eclipsing EL CVn-type binary, consisting of a $\delta$ Sct-type pulsator and a thermally bloated low-mass pre-He white dwarf (WD). Spectroscopy and time-series $BV$ photometry of V1224 Cas were carried out. The spectroscopy reveals a spectral type of A3 for the star. Light-curve modelling indicate that V1224 Cas is a short-period detached system containing a possible low-mass WD with an effective temperature of about 9516 K. Based on the effective temperature and the surface gravity of the A-type primary star from the spectroscopic results, the absolute parameters of the components were estimated as: $M_{P}$=2.16$\pm$0.22$M_{\sun}$, $R_{P}$=3.54$\pm$0.12$R_{\sun}$, $L_{P}$=55.9$\pm$6.9$L_{\sun}$, $M_{S}$=0.19$\pm$0.02$M_{\sun}$, $R_{S}$=0.97$\pm$0.04$R_{\sun}$, and $L_{S}$=6.9$\pm$0.6$L_{\sun}$. We therefore introduce V1224 Cas as a new EL CVn-type binary candidate. The light curves in both filters all show multi-periodic pulsations, superimposed on binary effects. We performed a preliminary frequency analysis of the light residuals after removing the synthetic eclipsing curve from the original observational data. The results suggest that the rapid light variations among the light curves could be attributed to the $\delta$ Sct-type primary component. We therefore conclude that V1224 Cas is very likely a WD+$\delta$ Sct binary.

Thirty-three fast radio bursts (FRBs) had been detected by March 2018. Although the sample size is still limited, meaningful statistical studies can already be carried out. The normalised luminosity function places important constraints on the intrinsic power output, sheds light on the origin(s) of FRBs, and can guide future observations. In this paper, we measure the normalised luminosity function of FRBs. Using Bayesian statistics, we can naturally account for a variety of factors such as receiver noise temperature, bandwidth, and source selection criteria. We can also include astronomical systematics, such as host galaxy dispersion measure, FRB local dispersion measure, galaxy evolution, geometric projection effects, and Galactic halo contribution. Assuming a Schechter luminosity function, we show that the isotropic luminosities of FRBs have a power-law distribution that covers approximately three orders of magnitude, with a power-law index ranging from $-1.8$ to $-1.2$ and a cut off $\sim 2\times 10^{44}\,\rm erg\,s^{-1}$. By using different galaxy models and well-established Bayesian marginalisation techniques, we show that our conclusions are robust against unknowns, such as the electron densities in the Milky Way halo and the FRB environment, host galaxy morphology, and telescope beam response.

We measure the clustering of DES Year 1 galaxies that are intended to be combined with weak lensing samples in order to produce precise cosmological constraints from the joint analysis of large-scale structure and lensing correlations. Two-point correlation functions are measured for a sample of $6.6 \times 10^{5}$ luminous red galaxies selected using the \textsc{redMaGiC} algorithm over an area of $1321$ square degrees, in the redshift range $0.15 < z < 0.9$, split into five tomographic redshift bins. The sample has a mean redshift uncertainty of $\sigma_{z}/(1+z) = 0.017$. We quantify and correct spurious correlations induced by spatially variable survey properties, testing their impact on the clustering measurements and covariance. We demonstrate the sample's robustness by testing for stellar contamination, for potential biases that could arise from the systematic correction, and for the consistency between the two-point auto- and cross-correlation functions. We show that the corrections we apply have a significant impact on the resultant measurement of cosmological parameters, but that the results are robust against arbitrary choices in the correction method. We find the linear galaxy bias in each redshift bin in a fiducial cosmology to be $b(z$=$0.24)=1.40 \pm 0.08$, $b(z$=$0.38)=1.61 \pm 0.05$, $b(z$=$0.53)=1.60 \pm 0.04$ for galaxies with luminosities $L/L_*>$$0.5$, $b(z$=$0.68)=1.93 \pm 0.05$ for $L/L_*>$$1$ and $b(z$=$0.83)=1.99 \pm 0.07$ for $L/L_*$$>1.5$, broadly consistent with expectations for the redshift and luminosity dependence of the bias of red galaxies. We show these measurements to be consistent with the linear bias obtained from tangential shear measurements.

We present the results of a search for rapidly evolving transients in the Dark Energy Survey Supernova Programme. These events are characterized by fast light curve evolution (rise to peak in $\lesssim 10$ d and exponential decline in $\lesssim30$ d after peak). We discovered 72 events, including 37 transients with a spectroscopic redshift from host galaxy spectral features. The 37 events increase the total number of rapid optical transients by more than factor of two. They are found at a wide range of redshifts ($0.05<z<1.56$) and peak brightnesses ($-15.75>M_\mathrm{g}>-22.25$). The multiband photometry is well fit by a blackbody up to few weeks after peak. The events appear to be hot ($T\approx10000-30000$ K) and large ($R\approx 10^{14}-2\cdot10^{15}$ cm) at peak, and generally expand and cool in time, though some events show evidence for a receding photosphere with roughly constant temperature. Spectra taken around peak are dominated by a blue featureless continuum consistent with hot, optically thick ejecta. We compare our events with a previously suggested physical scenario involving shock breakout in an optically thick wind surrounding a core-collapse supernova (CCSNe), we conclude that current models for such a scenario might need an additional power source to describe the exponential decline. We find these transients tend to favor star-forming host galaxies, which could be consistent with a core-collapse origin. However, more detailed modeling of the light curves is necessary to determine their physical origin.

We present accretion disk size measurements for 15 luminous quasars at $0.7 \leq z \leq 1.9$ derived from $griz$ light curves from the Dark Energy Survey. We measure the disk sizes with continuum reverberation mapping using two methods, both of which are derived from the expectation that accretion disks have a radial temperature gradient and the continuum emission at a given radius is well-described by a single blackbody. In the first method we measure the relative lags between the multiband light curves, which provides the relative time lag between shorter and longer wavelength variations. From this, we are only able to constrain upper limits on disk sizes, as many are consistent with no lag the 2$\sigma$ level. The second method fits the model parameters for the canonical thin disk directly rather than solving for the individual time lags between the light curves. Our measurements demonstrate good agreement with the sizes predicted by this model for accretion rates between 0.3-1 times the Eddington rate. Given our large uncertainties, our measurements are also consistent with disk size measurements from gravitational microlensing studies of strongly lensed quasars, as well as other photometric reverberation mapping results, that find disk sizes that are a factor of a few ($\sim$3) larger than predictions.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

In most astrophysical processes involving synchrotron radiation, the pitch-angle distribution of the electrons is assumed to be isotropic. However, if electrons are accelerated anisotropically, e.g, in a relativistic shock wave with an ordered magnetic field or in magnetic reconnection regions, the electron pitch angles might be anisotropic. In this work, we study synchrotron radiation from electrons with a pitch-angle distribution with respect to a large-scale uniform magnetic field. Assuming that the pitch-angle distribution is normal with a scatter of $\sigma_p$ and that the viewing direction is where the pitch-angle direction peaks, we find that for electrons with a Lorentz factor $\gamma$, the observed flux satisfies $F_\nu\propto\nu^{2/3}$ for $\nu\ll\nu_{\rm cr}$ ($\nu_{\rm cr}$ is the critical frequency of synchrotron), if $\sigma_p\lesssim1/\gamma$ is satisfied. On the other hand, if $\sigma_p\gg1/\gamma$, the spectrum below $\nu_{\rm cr}$ is a broken power law with a break frequency $\nu_{\rm br}\sim2\nu_{\rm cr}/\sigma_p^3\gamma^3$, e.g., $F_\nu\propto\nu^{2/3}$ for $\nu\ll\nu_{\rm br}$ and $F_\nu\propto\nu^{1/3}$ for $\nu_{\rm br}\ll\nu\ll\nu_{\rm cr}$. Thus the ultimate synchrotron line of death is $F_\nu\propto\nu^{2/3}$. We discuss the application of this theory to blazars and gamma-ray bursts (GRBs).

The spectral evolution and spectral lag behavior of 92 bright pulses from 84 gamma-ray bursts (GRBs) observed by the Fermi GBM telescope are studied. These pulses can be classified into hard-to-soft pulses (H2S, 64/92), H2S-dominated-tracking pulses (21/92), and other tracking pulses (7/92). We focus on the relationship between spectral evolution and spectral lags of H2S and H2S-dominated-tracking pulses. %in hard-to-soft pulses (H2S, 64/92) and H2S-dominating-tracking (21/92) pulses. The main trend of spectral evolution (lag behavior) is estimated with $\log E_p\propto k_E\log(t+t_0)$ (${\hat{\tau}} \propto k_{\hat{\tau}}\log E$), where $E_p$ is the peak photon energy in the radiation spectrum, $t+t_0$ is the observer time relative to the beginning of pulse $-t_0$, and ${\hat{\tau}}$ is the spectral lag of photons with energy $E$ with respect to the energy band $8$-$25$ keV. For H2S and H2S-dominated-tracking pulses, a weak correlation between $k_{{\hat{\tau}}}/W$ and $k_E$ is found, where $W$ is the pulse width. We also study the spectral lag behavior with peak time $t_{\rm p_E}$ of pulses for 30 well-shaped pulses and estimate the main trend of the spectral lag behavior with $\log t_{\rm p_E}\propto k_{t_p}\log E$. It is found that $k_{t_p}$ is correlated with $k_E$. We perform simulations under a phenomenological model of spectral evolution, and find that these correlations are reproduced. We then conclude that spectral lags are closely related to spectral evolution within the pulse. The most natural explanation of these observations is that the emission is from the electrons in the same fluid unit at an emission site moving away from the central engine, as expected in the models invoking magnetic dissipation in a moderately-high-$\sigma$ outflow.

Electron anti-neutrinos at the Glashow resonance (GR, at $E_{\bar \nu_e} \sim 6.3$ PeV) have an enhanced probability to be detected. With three neutrinos detected by IceCube in the (1-2) PeV energy range at present, one would expect that about 1 to 4 GR $\bar\nu_e$ should have been detected. The high-energy $\sim 8.7$ PeV muon neutrino detected by IceCube may not be a GR event. If so, we expect to detect 50 to 70 GR $\bar\nu_e$, then one would have a "missing Glashow-resonance problem". This would suggest (1) that $p\gamma$ interaction rather than $pp$ interaction is the dominant channel to produce the observed IceCube high-energy neutrinos, (2) that multi-pion $p\gamma$ interactions are suppressed, and (3) that the magnetic field and photon energy density in the $p\gamma$ emission region is such that significant $\mu^+$ cooling occurs before decaying, yet $\pi^+$'s essentially do not cool before decaying.

The groundbreaking discovery of the optical transient AT2017gfo associated with GW170817 opens a unique opportunity to study the physics of double neutron star (NS) mergers. We argue that the standard interpretation of AT2017gfo as being powered by radioactive decays of r-process elements faces the challenge of simultaneously accounting for the peak luminosity and peak time of the event, as it is not easy to achieve the required high mass, and especially the low opacity of the ejecta required to fit the data. A plausible solution would be to invoke an additional energy source, which is probably provided by the merger product. We consider energy injection from two types of the merger products: (1) a post-merger black hole powered by fallback accretion; and (2) a long-lived NS remnant. The former case can only account for the early emission of AT2017gfo, with the late emission still powered by radioactive decay. In the latter case, both early- and late-emission components can be well interpreted as due to energy injection from a spinning-down NS, with the required mass and opacity of the ejecta components well consistent with known numerical simulation results. We suggest that there is a strong indication that the merger product of GW170817 is a long-lived (supramassive or even permanently stable), low magnetic field NS. The result provides a stringent constraint on the equations of state of NSs.

The favored progenitor model for Gamma-ray Bursts (GRBs) with Supernova (SN) association is the core collapse of massive stars. One possible outcome of such a collapse is a rapidly spinning, strongly magnetized neutron star ("magnetar"). We systematically analyze the multi-wavelength data of GRB/SN associations detected by several instruments before 2017 June. Twenty GRB/SN systems have been confirmed via direct spectroscopic evidence or a clear light curve bump, as well as some spectroscopic evidence resembling a GRB-SN. We derive/collect the basic physical parameters of the GRBs and the SNe, and look for correlations among these parameters. We find that the peak brightness, $\rm ^{56}Ni$ mass, and explosion energy of SNe associated with GRBs are statistically higher than other Type Ib/c SNe. A statistically significant relation between the peak energy of GRBs and the peak brightness of their associated SNe is confirmed. No significant correlations are found between the GRB energies (either isotropic or beaming-corrected) and the supernova energy. We investigate the energy partition within these systems and find that the beaming-corrected GRB energy of most systems is smaller than the SN energy, with less than 30\% of the total energy distributed in the relativistic jet. The total energy of the systems is typically smaller than the maximum available energy of a millisecond magnetar ($2\times 10^{52}$ erg), especially if aspherical SN explosions are considered. The data are consistent with-though not a proof of-the hypothesis that most, but not all, GRB/SN systems are powered by millisecond magnetars.

Context. The Gaia mission has released the second data set (Gaia DR2), which contains parallaxes and proper motions for a large number of massive, young stars. Aims. We investigate the spiral structure in the solar neighborhood revealed by Gaia DR2 and compare it with that depicted by VLBI maser parallaxes. Methods. We examined three samples with different constraints on parallax uncertainty and distance errors and stellar spectral types: (1) all OB stars with parallax errors of less than 10%; (2) only O-type stars with 0.1 mas errors imposed and with parallax distance errors of less than 0.2 kpc; and (3) only O-type stars with 0.05 mas errors imposed and with parallax distance errors of less than 0.3 kpc. Results. In spite of the significant distance uncertainties for stars in DR2 beyond 1.4 kpc, the spiral structure in the solar neighborhood demonstrated by Gaia agrees well with that illustrated by VLBI maser results. The O-type stars available from DR2 extend the spiral arm models determined from VLBI maser parallaxes into the fourth Galactic quadrant, and suggest the existence of a new spur between the Local and Sagittarius arms.

The first gravitational-wave event from the merger of a binary neutron star system (GW170817) was detected recently. The associated short gamma-ray burst (GRB 170817A) has a low isotropic luminosity ($\sim 10^{47}$ erg s$ ^{-1}$) and a peak energy $E_{p}\sim 145$ keV during the initial main emission between -0.3 and 0.4 s. The origin of this short GRB is still under debate, but a plausible interpretation is that it is due to the off-axis emission from a structured jet. We consider two possibilities. First, since the best-fit spectral model for the main pulse of GRB 170817A is a cutoff power law with a hard low-energy photon index ($\alpha =-0.62_{-0.54}^{+0.49} $), we consider an off-axis photosphere model. We develop a theory of photosphere emission in a structured jet and find that such a model can reproduce a low-energy photon index that is softer than a blackbody through enhancing high-latitude emission. The model can naturally account for the observed spectrum. The best-fit Lorentz factor along the line of sight is $\sim 20$, which demands that there is a significant delay between the merger and jet launching. Alternatively, we consider that the emission is produced via synchrotron radiation in an optically thin region in an expanding jet with decreasing magnetic fields. This model does not require a delay of jet launching but demands a larger bulk Lorentz factor along the line of sight. We perform Markov Chain Monte Carlo fitting to the data within the framework of both models and obtain good fitting results in both cases.

The first gravitational-wave event from the merger of a binary neutron star system (GW170817) was detected recently. The associated short gamma-ray burst (GRB 170817A) has a low isotropic luminosity ($\sim 10^{47}$ erg s$ ^{-1}$) and a peak energy $E_{p}\sim 145$ keV during the initial main emission between -0.3 and 0.4 s. The origin of this short GRB is still under debate, but a plausible interpretation is that it is due to the off-axis emission from a structured jet. We consider two possibilities. First, since the best-fit spectral model for the main pulse of GRB 170817A is a cutoff power law with a hard low-energy photon index ($\alpha =-0.62_{-0.54}^{+0.49} $), we consider an off-axis photosphere model. We develop a theory of photosphere emission in a structured jet and find that such a model can reproduce a low-energy photon index that is softer than a blackbody through enhancing high-latitude emission. The model can naturally account for the observed spectrum. The best-fit Lorentz factor along the line of sight is $\sim 20$, which demands that there is a significant delay between the merger and jet launching. Alternatively, we consider that the emission is produced via synchrotron radiation in an optically thin region in an expanding jet with decreasing magnetic fields. This model does not require a delay of jet launching but demands a larger bulk Lorentz factor along the line of sight. We perform Markov Chain Monte Carlo fitting to the data within the framework of both models and obtain good fitting results in both cases.

New photometric observations of NSVS 01286630 were performed and two sets of four-color (B, V , Rc, Ic) light curves (LCs) were obtained. Using the 2013 version of the Wilson-Devinney (W-D) code, we analyzed these data. The photometric solutions reveal that NSVS 01286630 is an active detached eclipsing binary (EB) with a high orbit inclination. Remarkably, the temperature of the primary component (the hotter star) is higher than the secondary one, but, the value of mass ratio q(M2/M1) for NSVS 0128663 is more than 1, which can be explained that the surface of the secondary component of NSVS 01286630 is covered with big cool-star-spots. Based on our new CCD mid-eclipse times and the data published until now, the variations of the mid-eclipse times were reanalyzed in detail using a weighted least-squares method. It is discovered that the (O-C) diagram of the system shows a cyclic oscillation with a period of 3.61 yr and an amplitude of 0.001 days. The cyclic variation may be caused by the light travel time effect (LTTE) due to the presence of a third companion, whose mass we calculated is about 0.11 solar mass. The third body may affect the orbital evolution of the central binary system by transferring the angular momentum.

We describe the derivation and validation of redshift distribution estimates and their uncertainties for the galaxies used as weak lensing sources in the Dark Energy Survey (DES) Year 1 cosmological analyses. The Bayesian Photometric Redshift (BPZ) code is used to assign galaxies to four redshift bins between z=0.2 and 1.3, and to produce initial estimates of the lensing-weighted redshift distributions $n^i_{PZ}(z)$ for bin i. Accurate determination of cosmological parameters depends critically on knowledge of $n^i$ but is insensitive to bin assignments or redshift errors for individual galaxies. The cosmological analyses allow for shifts $n^i(z)=n^i_{PZ}(z-\Delta z^i)$ to correct the mean redshift of $n^i(z)$ for biases in $n^i_{\rm PZ}$. The $\Delta z^i$ are constrained by comparison of independently estimated 30-band photometric redshifts of galaxies in the COSMOS field to BPZ estimates made from the DES griz fluxes, for a sample matched in fluxes, pre-seeing size, and lensing weight to the DES weak-lensing sources. In companion papers, the $\Delta z^i$ are further constrained by the angular clustering of the source galaxies around red galaxies with secure photometric redshifts at 0.15<z<0.9. This paper details the BPZ and COSMOS procedures, and demonstrates that the cosmological inference is insensitive to details of the $n^i(z)$ beyond the choice of $\Delta z^i$. The clustering and COSMOS validation methods produce consistent estimates of $\Delta z^i$, with combined uncertainties of $\sigma_{\Delta z^i}=$0.015, 0.013, 0.011, and 0.022 in the four bins. We marginalize over these in all analyses to follow, which does not diminish the constraining power significantly. Repeating the photo-z procedure using the Directional Neighborhood Fitting (DNF) algorithm instead of BPZ, or using the $n^i(z)$ directly estimated from COSMOS, yields no discernible difference in cosmological inferences.

Recent work has found that the metal-poor star HD 94028 shows interesting abundance features. The s-process material (e.g., Ba, La, Ce, and Pb) and r-process material (e.g., Eu, Os, Pt) are mildly overabundant while the element C is not enhanced. However, the observed supersolar ratio [As/Ge](= $0.99 \pm 0.23$) and subsolar ratio [Se/As](= $-0.16 \pm 0.22$) could not be fitted by the combination of s- and r-process abundances. In this work, adopting the abundance decomposition method, we investigate the astrophysical origins of the elements in this star. We find that the primary process of massive stars dominates the production of the light elements and iron-group elements. The lighter neutron-capture elements from Ge to Te mainly come from the weak r-process. The elements from Ba to Pb are synthesized dominantly by the main s- and main r-process. The observed abundance features of Ge, As, and Se in HD 94028 are mainly produced by the weak r-process, which possesses the features of the supersolar ratio [As/Ge] and subsolar ratio [Se/As]. Because Ge is not synthesized solely by the neutron-capture process, As should be the lightest neutron-capture element. Although the i-process has been suggested in several astrophysical environments, it should be superfluous to explaining the abundances of the lighter neutron-capture elements in HD 94028.

The abundance scatter of heavier r-process elements (Z > 56) relative to Fe ([r/Fe]) in metal-poor stars preserves excellent information of the star formation history and provides important insights into the various situations of the Galactic chemical enrichment. In this respect, the upper and lower boundaries of [r/Fe] could present useful clues for investigating the extreme situations of the star formation history and the early Galactic chemical evolution. In this paper, we investigate the formation of the upper and lower boundaries of [r/Fe] for the gas clouds. We find that, for a cloud from which metal-poor stars formed, the formation of the upper limits of [r/Fe] is mainly due to the pollution from a single main r-process event. For a cloud from which metal-poor stars formed, the formation of the lower limits of [r/Fe] is mainly due to the pollution from a single SN II event that ejects primary Fe.

Afterglows of Gamma-Ray Bursts (GRBs) are simple in the most basic model, but can show many complex features. The ultra-long duration GRB 111209A, one of the longest GRBs ever detected, also has the best-monitored afterglow in this rare class of GRBs. We want to address the question whether GRB 111209A was a special event beyond its extreme duration alone, and whether it is a classical GRB or another kind of high-energy transient. The afterglow may yield significant clues. We present afterglow photometry obtained in seven bands with the GROND imager as well as in further seven bands with the UVOT telescope on-board the Neil Gehrels Swift Observatory. The light curve is analysed by multi-band modelling and joint fitting with power-laws and broken power-laws, and we use the contemporaneous GROND data to study the evolution of the spectral energy distribution. We compare the optical afterglow to a large ensemble we have analysed in earlier works, and especially to that of another ultra-long event, GRB 130925A. We furthermore undertake a photometric study of the host galaxy. We find a strong, chromatic rebrightening event at approx 0.8 days after the GRB, during which the spectral slope becomes redder. After this, the light curve decays achromatically, with evidence for a break at about 9 days after the trigger. The afterglow luminosity is found to not be exceptional. We find that a double-jet model is able to explain the chromatic rebrightening. The afterglow features have been detected in other events and are not unique. The duration aside, the GRB prompt emission and afterglow parameters of GRB 111209A are in agreement with the known distributions for these parameters. While the central engine of this event may differ from that of classical GRBs, there are multiple lines of evidence pointing to GRB 111209A resulting from the core-collapse of a massive star with a stripped envelope.

Recent results from the Pierre Auger Collaboration have shown that the composition of ultrahigh-energy cosmic rays (UHECRs) becomes gradually heavier with increasing energy. Although gamma-ray bursts (GRBs) have been promising sources of UHECRs, it is still unclear whether they can account for the Auger results because of their unknown nuclear composition of ejected UHECRs. In this work, we revisit the possibility that low-luminosity GRBs (LL GRBs) act as the sources of UHECR nuclei, and give new predictions based on the intrajet nuclear composition models considering progenitor dependencies. We find that the nuclear component in the jet can be divided into two groups according to the mass fraction of silicon nuclei, Si-free and Si-rich. Motivated by the connection between LL GRBs and transrelativistic supernovae, we also consider the hypernova ejecta composition. Then, we discuss the survivability of UHECR nuclei in the jet base and internal shocks of the jets, and show that it is easier for nuclei to survive for typical LL GRBs. Finally, we numerically propagate UHECR nuclei ejected from LL GRBs with different composition models and compare the resulting spectra and composition to Auger data. Our results show that both the Si-rich progenitor and hypernova ejecta models match the Auger data well, while the Si-free progenitor models have more difficulty in fitting the spectrum. We argue that our model is consistent with the newly reported cross correlation between the UHECRs and starburst galaxies, since both LL GRBs and hypernovae are expected to be tracers of the star-formation activity. LL GRBs have also been suggested as the dominant origin of IceCube neutrinos in the PeV range, and the LL GRB origin of UHECRs can be critically tested by near-future multimessenger observations.

Gamma-ray burst (GRB) 120729A was detected by Swift/BAT and Fermi/GBM, and then rapidly observed by Swift/XRT, Swift/UVOT, and ground-based telescopes. It had a single long and smooth \gamma-ray emission pulse, which extends continuously to the X-rays. We report Lick/KAIT observations of the source, and make temporal and spectral joint fits of the multiwavelength light curves of GRB 120729A. It exhibits achromatic light-curve behavior, consistent with the predictions of the external shock model. The light curves are decomposed into four typical phases: onset bump (Phase I), normal decay (Phase II), shallow decay (Phase III), and post-jet break (Phase IV). The spectral energy distribution (SED) evolves from prompt \gamma-ray emission to the afterglow with photon index from $\Gamma_{\rm \gamma}=1.36$ to $\Gamma \approx 1.75$. There is no obvious evolution of the SED during the afterglow. The multiwavelength light curves from \gamma-ray to optical can be well modeled with an external shock by considering energy injection, and a time-dependent microphysics model with $\epsilon_B\propto t^{\alpha_B}$ for the emission at early times, $T < T_{\rm 0} + 157$~s. Therefore, we conclude that both the prompt \gamma-ray emission and afterglow of GRB 120729A have the same external shock physical origin. Our model indicates that the $\epsilon_B$ evolution can be described as a broken power-law function with $\alpha_{\rm B,1} = 0.18 \pm 0.04$ and $\alpha_{\rm B,2} = 0.84 \pm 0.04$. We also systematically investigate single-pulse GRBs in the Swift era, finding that only a small fraction of GRBs (GRBs 120729A, 051111, and 070318) are likely to originate from an external shock for both the prompt \gamma-ray emission and afterglow.

Gamma-ray Burst (GRB) collimation has been inferred with the observations of achromatic steepening in GRB light curves, known as jet breaks. Identifying a jet break from a GRB afterglow lightcurve allows a measurement of the jet opening angle and true energetics of GRBs. In this paper, we reinvestigate this problem using a large sample of GRBs that have an optical jet break which is consistent with being achromatic in the X-ray band. Our sample includes 99 GRBs from February 1997 to March 2015 that have optical and, for Swift GRBs, X-ray lightcurves that are consistent with the jet break interpretation. Out of 99 GRBs we have studied, 55 GRBs are found to have temporal and spectral behaviors both before and after the break consistent with the theoretical predictions of the jet break models, respectively. These include 53 long/soft (Type II) and 2 short/hard (Type I) GRBs. Only 1 GRB is classified as the candidate of a jet break with energy injection. Another 41 and 3 GRBs are classified as the candidates with the lower and upper limits of the jet break time, respectively. The typical beaming correction factor $f_b^{-1} \sim 1000$ for Type II GRBs, suggesting an even higher total GRB event rate density in the universe. Both isotropic and jet-corrected energies have a wide span in their distributions. We also investigate several empirical correlations (Amati, Frail, Ghirlanda and Liang-Zhang) previously discussed in the literature. We find that in general most of these relations are less tight than before. The existence of early jet breaks and hence small opening angle jets, which were detected in the {\em Swfit era}, is most likely the source of scatter. If one limits the sample to jet breaks later than $10^4$ s, the Liang-Zhang relation remains tight and the Ghirlanda relation still exists. These relations are derived from Type II GRBs, and Type I GRBs usually deviate from them.

We report astrometric results of VLBI phase-referencing observations of 22 GHz \hho\ masers emission toward the red hypergiant \vxsgr, one of most massive and luminous red hypergiant stars in our Galaxy, using the Very Long Baseline Array. A background source, \Jtwoze, projected 4\d4 from the target \vxsgr, was used as the phase reference. For the low declinations of these sources, such a large separation normally would seriously degrade the relative astrometry. We use a two-step method of tropospheric delay calibration, which combines the VLBI geodetic-block (or GPS) calibration with an image-optimization calibration, to obtain a trigonometric parallax of $0.64\pm0.04$ mas, corresponding to a distance of 1.56$^{+0.11}_{-0.10}$ kpc. The measured proper motion of \vxsgr\ is $0.36\pm0.76$ and $-2.92\pm0.78$ \masy\ in the eastward and northward directions. The parallax and proper motion confirms that \vxsgr\ belong to the Sgr OB1 association. Rescaling bolometric luminosities in the literature to our parallax distance, we find the luminosity of \vxsgr~is $(1.95 \pm 0.62) \times 10^5$ \Lsun, where the uncertainty is dominated by differing photometry measurements.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

Rapid response and short time latency are very important for Time Domain Astronomy, such as the observations of Gamma-ray Bursts (GRBs) and electromagnetic (EM) counterparts of gravitational waves (GWs). Based on the near real-time Fermi/GBM data, we developed a low-latency pipeline to automatically calculate the temporal and spectral properties of GRBs. With this pipeline, some important parameters can be obtained, such as T90 and fluence, within ~20 minutes after the GRB trigger. For ~90% GRBs, T90 and fluence are consistent with the GBM catalog results within 2 sigma errors. This pipeline has been used by the Gamma-ray Bursts Polarimeter (POLAR) and the Insight Hard X-ray Modulation Telescope (Insight-HXMT) to follow up the bursts of interest. For GRB 170817A, the first EM counterpart of GW events detected by Fermi/GBM and INTEGRAL/SPI-ACS, the pipeline gave T90 and spectral information in 21 minutes after the GBM trigger, providing important information for POLAR and Insight-HXMT observations.

We present new criteria for selecting HII regions from the Infrared Astronomical Satellite (IRAS) Point Source catalogue (PSC), based on an HII region catalogue derived manually from the all-sky Wide-field Infrared Survey Explorer (WISE). The criteria are used to augment the number of HII region candidates in the Milky Way. The criteria are defined by the linear decision boundary of two samples: IRAS point sources associated with known HII regions, which serve as the HII region sample, and IRAS point sources at high Galactic latitudes, which serve as the non-HII region sample. A machine learning classifier, specifically a support vector machine (SVM), is used to determine the decision boundary. We investigate all combinations of four IRAS bands and suggest that the optimal criterion is log(F$_{\rm 60}$/F$_{\rm 12}$)$\ge$(-0.19$\times$log(F$_{\rm 100}$/F$_{\rm 25}$)+ 1.52), with detections at 60 and 100 micron. This selects 3041 HII region candidates from the IRAS PSC. We find that IRAS HII region candidates show evidence of evolution on the two-colour diagram. Merging the WISE HII catalogue with IRAS HII region candidates, we estimate a lower limit of approximately 10200 for the number of HII regions in the Milky Way.

The origin of fast radio bursts (FRBs) remains mysterious. Recently, the only repeating FRB source, FRB 121102, was reported to possess an extremely large and variable rotation measure (RM). The inferred magnetic field strength in the burst environment is comparable to that in the vicinity of the super-massive black hole Sagittarius A* of our Galaxy. Here we show that all the observational properties of FRB 121102 (including the high RM and its evolution, the high linear polarization degree, an invariant polarization angle across each burst and other properties previously known) can be interpreted within the cosmic comb model invoking a regular neutron star whose magnetosphere is sporadically and marginally combed by a variable outflow from a nearby low-luminosity accreting super-massive black hole in the host galaxy. We propose three falsifiable predictions of the model and discuss other FRBs within the context of the cosmic comb model as well as the challenges encountered by other repeating FRB models in light of the new observations.

The repeating FRB 121102 (the "repeater") shows {repetitive bursting activities} and was localized in a host galaxy at $z=0.193$. On the other hand, despite dozens of hours of telescope time spent on follow-up observations, no other FRBs have been observed to repeat. Yet, it has been speculated that the repeater is the prototype of FRBs, and that other FRBs should show similar repeating patterns. Using the published data, we compare the repeater with other FRBs in the observed time interval ($\Delta t$) - flux ratio ($S_i / S_{i+1}$) plane. We find that whereas other FRBs occupy the upper (large $S_i / S_{i+1}$) and right (large $\Delta t$) regions of the plane due to the non-detections of other bursts, some of the repeater bursts fall into the lower-left region of the plot (short interval and small flux ratio) excluded by the non-detection data of other FRBs. The trend also exists even if one only selects those bursts detectable by the Parkes radio telescope. If other FRBs were similar to the repeater, our simulations suggest that the probability that none of them have been detected to repeat with the current searches would be $\sim (10^{-4}-10^{-3})$. We suggest that the repeater is not representative of the entire FRB population and that there is strong evidence of more than one population of FRBs.

Limited by the sensitivities of the current gravitational wave (GW) detectors, the central remnant of the binary neutron star (NS) merger associated with GW170817 remains an open question. Considering the relatively large total mass, it is generally proposed that the merger of GW170817 would lead to a shortly lived hypermassive NS or directly produce a black hole (BH). There is no clear evidence to support or rule out a long-lived NS as the merger remnant. Here we utilize the GW and electromagnetic (EM) signals to comprehensively investigate the parameter space that allows a long-lived NS to survive as the merger remnant of GW170817. We find that for some stiff equations of state, the merger of GW170817 could, in principle, lead to a massive NS, which has a millisecond spin period. The post-merger GW signal could hardly constrain the ellipticity of the NS. If the ellipticity reaches $10^{-3}$, in order to be compatible with the multi-band EM observations, the dipole magnetic field of the NS ($B_p$) is constrained to the magnetar level of $\sim10^{14}$ G. If the ellipticity is smaller than $10^{-4}$, $B_p$ is constrained to the level of $\sim10^{10}-10^{12}\,$G. These conclusions weakly depend on the adoption of equations of state.

Recently, the LIGO-Virgo collaboration reported their first detection of gravitational wave (GW) signals from a low mass compact binary merger GW170817, which is most likely due to a double neutron star (NS) merger. With the GW signals only, the chirp mass of the binary is precisely constrained to $1.188^{+0.004}_{-0.002}~\rm{M_{\odot}}$, but the mass ratio is loosely constrained in the range $0.4-1$, so that a very rough estimation of the individual NS masses ($0.86~{\rm M_{\odot}}<M_1<1.36~\rm{M_{\odot}}$ and $1.36~{\rm M_{\odot}}<M_2<2.26~\rm{M_{\odot}}$) was obtained. Here we propose that if one can constrain the dynamical ejecta mass through performing kilonova modeling of the optical/IR data, by utilizing an empirical relation between the dynamical ejecta mass and the mass ratio of NS binaries, one may place a more stringent constraint on the mass ratio of the system. For instance, considering that the red "kilonova" component is powered by the dynamical ejecta, we reach a tight constraint on the mass ratio in the range of $0.46-0.59$. Alternatively, if the blue "kilonova" component is powered by the dynamical ejecta, the mass ratio would be constrained in the range of $0.53-0.67$. Overall, such a multi-messenger approach could narrow down the mass ratio of GW170817 system to the range of $0.46-0.67$, which gives a more precise estimation of the individual NS mass than pure GW signal analysis, i.e. $0.90~{\rm M_{\odot}}<M_1<1.16~{\rm M_{\odot}}$ and $1.61~{\rm M_{\odot}}<M_2<2.11~{\rm M_{\odot}}$.

The remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions.

We determined the \fodd\ values, $0.46\pm0.08$, $0.51\pm0.09$, $0.50\pm0.13$, $0.48\pm0.12$, which correspond to the r-contribution 100\% for four r-II stars, \cs, \hen, \hes\, and \het, respectively. Our results suggest that almost all of the heavy elements (in the range from Ba to Pb) in r-II stars have a common origin, that is, from a single r-process (the main r-process). We found that the \fodd\ has a intrinsic nature, and should keep constant value of about 0.46 in the main r-process yields, which is responsible for the heavy element enhancement of r-II stars and of our Galaxy chemical enhancement. In addition, except the abundance ratio [Ba/Eu] the \fodd\ is also an important indicator, which can be used to study the relative contributions of the r- and s-process during the chemical evolution history of the Milky Way and the enhancement mechanism in stars with peculiar abundance of heavy elements.

The "Breakthrough Starshot" aims at sending near-speed-of-light cameras to nearby stellar systems in the future. Due to the relativistic effects, a trans-relativistic camera naturally serves as a spectrograph, a lens, and a wide-field camera. We demonstrate this through a simulation of the optical-band image of the nearby galaxy M51 in the rest frame of the trans-relativistic camera. We suggest that observing celestial objects using a trans-relativistic camera may allow one to study the astronomical objects in a special way, and to perform unique tests on the principles of special relativity. We outline several examples that trans-relativistic cameras may make important contributions to astrophysics and suggest that the Breakthrough Starshot cameras may be launched in any direction to serve as a unique astronomical observatory.

We perform Monte Carlo simulations of transrelativistic shear acceleration dedicated to a jet-cocoon system of active galactic nuclei. A certain fraction of galactic cosmic rays in a halo is entrained, and sufficiently high-energy particles can be injected to the reacceleration process and further accelerated up to 100 EeV. We show that the shear reacceleration mechanism leads to a hard spectrum of escaping cosmic rays, $dL_E/dE\propto E^{-1}-E^0$, distinct from a conventional $E^{-2}$ spectrum. The supersolar abundance of ultrahigh-energy nuclei is achieved due to injections at TeV-PeV energies. As a result, we find that the highest-energy spectrum and mass composition can be reasonably explained by our model without contradictions with the anisotropy data.

We consider a simple model for gamma-ray bursts induced by magnetic reconnection in turbulent media. The magnetic field in a jet is subject to kink instabilities, which distort the regular structure of the spiral magnetic field, drive turbulence, and trigger reconnection. The resulting reconnection takes place in a high Reynolds number medium, where turbulence is further enhanced and in turn accelerates the reconnection process. This boot-strap reconnection gives rise to bursts of reconnection events, through which the free energy of magnetic field is transformed into a gamma-ray burst. The efficiency of magnetic reconnection and magnetic energy dissipation is not constrained by microphysical properties of plasmas.

We aim to ascertain the physical parameters of a propagating wave over the solar disk detected by the Interface Region Imaging Spectrograph (IRIS). Using imaging data from the IRIS and the Solar Dynamic Observatory (SDO), we tracked bright spots to determine the parameters of a propagating transverse wave in active region (AR) loops triggered by activation of a filament. Deriving the Doppler velocity of Si IV line from spectral observations of IRIS, we have determined the rotating directions of active region loops which are relevant to the wave. On 2015 December 19, a filament was located on the polarity inversion line of the NOAA AR 12470. The filament was activated and then caused a C 1.1 two-ribbon flare. Between the flare ribbons, two rotation motions of a set of bright loops were observed to appear in turn with opposite directions. Following the end of the second rotation, a propagating wave and an associated transverse oscillation were detected in these bright loops. In 1400 A channel, there was bright material flowing along the loops in a wave-like manner, with a period of ~128 s and a mean amplitude of ~880 km. For the transverse oscillation, we tracked a given loop and determine the transverse positions of the tracking loop in a limited longitudinal range. In both of 1400 A and 171 A channels, approximately four periods are distinguished during the transverse oscillation. The mean period of the oscillation is estimated as ~143 s and the displacement amplitude as between ~1370 km and ~690 km. We interpret these oscillations as a propagating kink wave and obtain its speed of ~1400 km s-1. Our observations reveal that a flare associated with filament activation could trigger a kink propagating wave in active region loops over the solar disk.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

In this paper, CCD photometric light curves for the short-period eclipsing binary 1SWASP J140533.33+114639.1 (hereafter J1405) in the $BVR$ bands are presented and analyzed using the 2013 version of the Wilson-Devinney (W-D) code. It is discovered that the J1405 is a W-subtype shallow contact binary with a contact degree of $f$ = 7.9$\pm$ 0.5\% and a mass ratio of $q$ = 1.55 $\pm$ 0.02. In order to explain the asymmetric light curves of the system, a cool star-spot on the more massive component was employed. This shallow contact eclipsing binary may be formed from a short-period detached system through the orbital shrinkage due to angular momentum loss. Based on $(O-C)$ method, the variation of the orbital period was studied using all the available times of the minimum light. The $(O-C)$ diagram reveals that the period is increasing continuously at a rate of $dP/dt=+2.09\times{10^{-7}}$days yr$^{-1}$, which can be explained by mass transfer from the less massive component to the more massive one.

Mergers of neutron star-neutron star (NS-NS) or neutron star-black hole (NS-BH) binaries are candidate sources of gravitational waves (GWs). At least a fraction of the merger remnant should be a stellar mass BH with a sub-relativistic ejecta. A collimated jet is launched via Blandford-Znajek mechanism from the central BH to trigger a short gamma-ray burst (SGRB). At the same time, a near-isotropic wind can be driven by the Blandford-Payne mechanism (BP). In previous work, additional energy injection to the ejecta from the BP mechanism was ignored, and radioactive decay has long been thought as the main source of energy for "kilonova" in the BH scenario. In this letter, we propose that the wind driven by BP mechanism from the new-born BH-disk can heat up and push the ejecta during the prompt emission phase or even at late time when there is mass fall-back. Such a BP-powered merger-nova could be bright at optical band even for a low-luminosity SGRB. The detection of a GW event with a merger product of BH, and accompanied by a bright merger-nova, would be a robust test of our model.

`Internal plateau' followed by a sharp decay is commonly seen in short gamma-ray burst (GRB) light curves. The plateau component is usually interpreted as the dipole emission from a supra-massive magnetar, and the sharp decay may imply the collapse of the magnetar to a black hole (BH). Fall-back accretion onto the new-born BH could produce long-lasting activities via the Blandford-Znajek (BZ) process. The magnetic flux accumulated near the BH would be confined by the accretion disks for a period of time. As the accretion rate decreases, the magnetic flux is strong enough to obstruct gas infall, leading to a magnetically-arrested disk (MAD). Within this scenario, we show that the BZ process could produce two types of typical X-ray light curves: type I exhibits a long-lasting plateau, followed by a power-law decay with slopes ranging from 5/3 to 40/9; type II shows roughly a single power-law decay with slope of 5/3. The former requires low magnetic filed strength, while the latter corresponds to relatively high values. We search for such signatures of the new-born BH from a sample of short GRBs with an internal plateau, and find two candidates: GRB 101219A and GRB 160821B, corresponding to type II and type I light curve, respectively. It is shown that our model can explain the data very well.

The prompt emission spectrum of Gamma-ray bursts (GRB) is characterized by a smoothly joint broken power-law spectrum known as the Band function. The typical low-energy photon index is $\sim -1$, which poses challenge to the standard synchrotron radiation models. We investigate the electron energy spectrum as a result of the interplay among adiabatic stochastic acceleration (ASA), particle injection, and synchrotron cooling. In the ASA-dominated low-energy range, ASA enables an efficient hardening of the injected energy spectrum to approach a spectral index $-1$. In the synchrotron cooling-dominated high-energy range, the injected high-energy electrons undergo fast synchrotron cooling and have a softer photon spectrum. With the energy range of the injected electrons broadly covering both the ASA- and synchrotron cooling-dominated ranges, the resulting photon number spectrum has the low- and high-energy indices as $\alpha_s \sim -1$ and $\beta_s \sim -p/2-1$, respectively. The break energy is of the order of $\sim 100 $ keV, depending on the turbulence properties.

The low-energy spectra of gamma-ray bursts' (GRBs) prompt emission are closely related to the energy distribution of electrons, which is further regulated by their cooling processes. We develop a numerical code to calculate the evolution of the electron distribution with given initial parameters, in which three cooling processes (i.e., adiabatic, synchrotron and inverse Compton cooling) and the effect of decaying magnetic field are coherently considered. A sequence of results are presented by exploring the plausible parameter space for both the fireball and the Poynting-flux-dominated regime. Different cooling patterns for the electrons can be identified and they are featured by a specific dominant cooling mechanism. Our results show that the hardening of the low-energy spectra can be attributed to the dominance of synchrotron self-Compton cooling within the internal shock model, or to decaying synchrotron cooling within the Poynting-flux-dominated jet scenario. These two mechanisms can be distinguished by observing the hard low-energy spectra of isolated short pulses in some GRBs. The dominance of adiabatic cooling can also lead to hard low-energy spectra when the ejecta is moving at an extreme relativistic speed. The information from the time-resolved low-energy spectra can help to probe the physical characteristics of the GRB ejecta via our numerical results.

Long-duration gamma-ray bursts (LGRBs) and ultra-LGRBs (ULGRBs) originate from collapsars, in the center of which a newborn rotating stellar-mass black hole (BH) surrounded by a massive accretion disk may form. In the scenario of BH hyperaccretion inflow-outflow model and Blandford-Znajek (BZ) mechanism to trigger gamma-ray bursts (GRBs), the real accretion rate to power a BZ jet is far lower than the mass supply rate from the progenitor star. The characteristics of the progenitor stars can be constrained by GRB luminosity observations, and the results exceed usual expectations. LGRBs lasting from several seconds to tens of seconds in the rest frame may originate from solar-metallicity ($Z \sim 1~ Z_\odot$, where $Z$ and $Z_\odot$ are the metallicities of progenitor stars and the Sun), massive ($M \gtrsim 34 ~M_\odot$, where $M$ and $M_\odot$ are the masses of progenitor stars and the Sun) stars or some zero-metallicity ($Z \sim 0$) stars. A fraction of low-metallicity ($Z \lesssim 10^{-2}~Z_\odot$) stars, including Population III stars, can produce ULGRBs such as GRB 111209A. The fraction of LGRBs lasting less than tens of seconds in the rest frame is more than 40$\%$, which cannot conform to the fraction of the demanded type of progenitor star. It possibly implies that the activity timescale of central engine may be much longer than the observed timescale of prompt emission phase, as indicated by X-ray late-time activities. Alternatively, LGRBs and ULGRBs may be powered by a millisecond magnetar central engine.

We conducted multi-epoch VLBA phase reference observations of LS I +61 303 in order to study its precessing radio jet. Compared to similar observations in 2006, we find that the observed elliptical trajectory of emission at 8.4 GHz repeats after the 9-year gap. The accurate alignment of the emission patterns yields a precession period of 26.926 +- 0.005 d, which is consistent with that determined by Lomb-Scargle analysis of the radio light curve. We analytically model the projection on the sky plane of the peak position of a precessing, synchrotron-emitting jet, which traces an elliptical trajectory on the sky. Comparing the simulation with the VLBA astrometry we improve our knowledge of the geometry of the system.We measure the LS I +61 303 absolute proper motion to be -0.150 +- 0.006 mas/yr eastward and -0.264 +- 0.006 mas/yr northward. Removing Galactic rotation, this reveals a small, < 20 km/s, non-circular motion, which indicates a very low kick velocity when the black hole was formed.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

A hyperaccreting stellar-mass black hole has been proposed as the candidate central engine of gamma-ray bursts (GRBs). The rich observations of GRBs by \textit{Fermi} and \textit{Swift} make it possible to constrain the central engine model by comparing the model predications against data. This paper is dedicated to studying the temporal evolution of central engine parameters for both prompt emission and afterglow phases. We consider two jet launching mechanisms, i.e., $\nu\bar{\nu}$ annihilations and the Blandford-Znajek (BZ) processe, and obtain analytical solutions to these two models. We then investigate the black hole central engine parameters, such as the jet power, the dimensionless entropy $\eta$, and the central engine parameter $\mu_0=\eta (1+\sigma_0)$ (where $\sigma_0$ is the initial magnetization of the engine) at the base of the jet. The black hole may be spun up by accretion, or spun down by the BZ process, leaving imprints in GRB lightcurves. Usually, a BZ jet is more powerful and is likely responsible for the late time central engine activities. However, an initially non-spinning black hole central engine may first launch a thermal "fireball" via neutrino annihilations, and then launch a Poynting-flux-dominated jet via the BZ process. Multiple flares, giant bumps, and plateaus in GRB afterglows can be well produced as the result of late time accretion onto the black hole.

We conducted an SiO maser survey towards 221 O-rich AGB stars with the aim of identifying maser emission associated with the Sagittarius stellar stream. In this survey, maser emission was detected in 44 targets, of which 35 were new detections. All of these masers are within 5 kpc of the Sun. We also compiled a Galactic SiO maser catalogue including ~2300 SiO masers from the literature. The distribution of these SiO masers give a scale height of 0.40 kpc, while 42 sources deviate from the Galactic plane by more than 1.2 kpc, half of which were found in this survey. Regarding SiO masers in the disc, we found both the rotational speeds and the velocity dispersions vary with the Galactic plane distance. Assuming Galactic rotational speed $\Theta$0 = 240 km/s , we derived the velocity lags are 15 km/s and 55 km/s for disc and off-plane SiO masers respectively. Moreover, we identified three groups with significant peculiar motions (with 70% confidence). The most significant group is in the thick disc that might trace stream/peculiar motion of the Perseus arm. The other two groups are mainly made up of off-plane sources. The northern and southern off-plane sources were found to be moving at ~33 km/s and 54 km/s away from the Galactic plane, respectively. Causes of these peculiar motions are still unclear. For the two off-plane groups, we suspect they are thick disc stars whose kinematics affected by the Sgr stellar stream or very old Sgr stream debris.

Tidal disruption events (TDEs) by supermassive or intermediate mass black holes have been suggested as candidate sources of ultrahigh-energy cosmic rays (UHECRs) and high-energy neutrinos. Motivated by the recent measurements from the Pierre Auger Observatory, which indicates a metal-rich cosmic-ray composition at ultrahigh energies, we investigate the fate of UHECR nuclei loaded in TDE jets. First, we consider the production and survival of UHECR nuclei at internal shocks, external forward and reverse shocks, and nonrelativistic winds. Based on the observations of Swift J1644+57, we show that the UHECRs can survive for external reverse and forward shocks, and disk winds. On the other hand, UHECR nuclei are significantly disintegrated in internal shocks, although they could survive for low-luminosity TDE jets. Assuming that UHECR nuclei can survive, we consider implications of different composition models of TDEs. We find that the tidal disruption of main sequence stars or carbon-oxygen white dwarfs does not successfully reproduce UHECR observations, namely the observed composition or spectrum. The observed mean depth of the shower maximum and its deviation could be explained by oxygen-neon-magnesium white dwarfs, but they may be too rare to be the sources of UHECRs.

We use a sample of 1338 spectroscopically confirmed and photometrically classified Type Ia Supernovae (SNe Ia), sourced from the CSP, CfA, SDSS-II, and SNLS supernova samples, to examine the relationships between SNe Ia and the galaxies that host them. Our results provide confirmation with improved statistical significance that SNe Ia, after standardization, are on average more luminous in massive hosts (significance $\rm > 5 \sigma$), and decline more rapidly in massive hosts (significance $\rm > 9\sigma$) and in hosts with low specific star formation rates (significance $\rm > 8\sigma$). We study the variation of these relationships with redshift and detect no evolution. We split SNe Ia into pairs of subsets that are based on the properties of the hosts, and fit cosmological models to each subset. Including both systematic and statistical uncertainties, we do not find any significant shift in the best-fit cosmological parameters between the subsets. Among different SN Ia subsets, we find that SNe Ia in hosts with high specific star formation rates have the least intrinsic scatter ($\rm \sigma_{int}=0.08\pm0.01$) in luminosity after standardization.

GRB 160625B, one of the brightest bursts in recent years, was simultaneously observed by Fermi and Swift satellites, and ground-based optical telescopes in three different events separated by long periods of time. In this paper the non-thermal multiwavelength observations of GRB 160625B are described and a transition phase from wind-type-like medium to interstellar medium between the early (event II) and the late (event III) afterglow is found. The multiwavelength observations of the early afterglow are consistent with the afterglow evolution starting at $\sim$ 150 s in a stellar wind medium whereas the observations of the late afterglow are consistent with the afterglow evolution in interstellar medium (ISM). The wind-to-ISM transition is calculated to be at $\sim 8\times 10^3$ s when the jet has decelerated, at a distance of $\sim$ 1 pc from the progenitor. Using the standard external shock model, the synchrotron and synchrotron self-Compton emission from reverse shock is required to model the GeV $\gamma$-ray and optical observations in the early afterglow, and synchrotron radiation from the adiabatic forward shock to describe the X-ray and optical observations in the late afterglow. The derived values of the magnetization parameter, the slope of the fast decay of the optical flash and the inferred magnetic fields suggest that Poynting flux-dominated jet models with arbitrary magnetization could account for the spectral properties exhibited by GRB 160625B.

GRB 160625B, one of the brightest bursts in recent years, was simultaneously observed by Fermi and Swift satellites, and ground-based optical telescopes in three different events separated by long periods of time. In this paper the non-thermal multiwavelength observations of GRB 160625B are described and a transition phase from wind-type-like medium to interstellar medium between the early (event II) and the late (event III) afterglow is found. The multiwavelength observations of the early afterglow are consistent with the afterglow evolution starting at $\sim$ 150 s in a stellar wind medium whereas the observations of the late afterglow are consistent with the afterglow evolution in interstellar medium (ISM). The wind-to-ISM transition is calculated to be at $\sim 8\times 10^3$ s when the jet has decelerated, at a distance of $\sim$ 1 pc from the progenitor. Using the standard external shock model, the synchrotron and synchrotron self-Compton emission from reverse shock is required to model the GeV $\gamma$-ray and optical observations in the early afterglow, and synchrotron radiation from the adiabatic forward shock to describe the X-ray and optical observations in the late afterglow. The derived values of the magnetization parameter, the slope of the fast decay of the optical flash and the inferred magnetic fields suggest that Poynting flux-dominated jet models with arbitrary magnetization could account for the spectral properties exhibited by GRB 160625B.

1SWASP\,J162117.36$+$441254.2 was originally classified as an EW-type binary with a period of 0.20785\,days. However, it was detected to have undergone a stellar outburst on June 3, 2016. Although the system was latter classified as a cataclysmic variable (CV) and the event was attributed as a dwarf-nova outburst, the physical reason is still unknown. This binary has been monitored photometrically since April 19, 2016 and many light curves were obtained before, during and after the outburst. Those light and color curves observed before the outburst indicate that the system is a special CV. The white dwarf is not accreting material from the secondary and there are no accretion disks surrounding the white dwarf. By comparing the light curves obtained from April 19 to September 14, 2016, it was found that magnetic activity of the secondary is associated with the outburst. We show strong evidence that the $L_1$ region on the secondary was heavily spotted before and after the outburst and thus quench the mass transfer, while the outburst is produced by a sudden mass accretion of the white dwarf. These results suggest that J162117 is a good astrophysical laboratory to study stellar magnetic activity and its influences on CV mass transfer and mass accretion.

We use numerical simulations to characterize the performance of a clustering-based method to calibrate photometric redshift biases. In particular, we cross-correlate the weak lensing (WL) source galaxies from the Dark Energy Survey Year 1 (DES Y1) sample with redMaGiC galaxies (luminous red galaxies with secure photometric redshifts) to estimate the redshift distribution of the former sample. The recovered redshift distributions are used to calibrate the photometric redshift bias of standard photo-$z$ methods applied to the same source galaxy sample. We apply the method to three photo-$z$ codes run in our simulated data: Bayesian Photometric Redshift (BPZ), Directional Neighborhood Fitting (DNF), and Random Forest-based photo-$z$ (RF). We characterize the systematic uncertainties of our calibration procedure, and find that these systematic uncertainties dominate our error budget. The dominant systematics are due to our assumption of unevolving bias and clustering across each redshift bin, and to differences between the shapes of the redshift distributions derived by clustering vs photo-$z$'s. The systematic uncertainty in the mean redshift bias of the source galaxy sample is $\Delta z \lesssim 0.02$, though the precise value depends on the redshift bin under consideration. We discuss possible ways to mitigate the impact of our dominant systematics in future analyses.

We introduce a non-resonant acceleration mechanism arising from the second adiabatic invariant in magnetic turbulence and apply it to study the prompt emission spectra of gamma-ray bursts (GRBs). The mechanism contains both the first- and second-order Fermi acceleration, originating from the interacting turbulent reconnection and dynamo processes. It leads to a hard electron energy distribution up to a cutoff energy at the balance between the acceleration and synchrotron cooling. The sufficient acceleration rate ensures a rapid hardening of any initial energy distribution to a power-law distribution with the index $p \sim 1$, which naturally produces a low-energy photon index $\alpha \sim -1$ via the synchrotron radiation. For typical GRB parameters, the synchrotron emission can extend to a characteristic photon energy on the order of $\sim 100$ keV.

The repeating fast radio burst (FRB) 121102 was recently localized in a dwarf galaxy at a cosmological distance. The dispersion measure (DM) derived for each burst from FRB 121102 so far has not shown significant evolution, even though an apparent increase was recently seen with newly detected VLA bursts. It is expected that more repeating FRB sources may be detected in the future. In this work, we investigate a list of possible astrophysical processes that might cause DM variation of a particular FRB source. The processes include (1) the cosmological scale effects such as Hubble expansion and large-scale structure fluctuations; (2) the FRB local effects such as gas density fluctuation, expansion of a supernova remnant, a pulsar wind nebula, and an HII region; and (3) the propagation effect due to plasma lensing. We find that the DM variations contributed by the large-scale structure are extremely small, and any observable DM variation is likely caused by the plasma local to the FRB source. Besides mechanisms that produce decreasing DM with time, we suggest that an FRB source in an expanding supernova remnant around a nearly neutral ambient medium during the deceleration (Sedov-Taylor and snowplow) phases or in a growing HII region can introduce DM increasing. Some effects (e.g. an FRB source moving in an HII region or plasma lensing) can give either positive or negative DM variations. Future observations of DM variations of FRB 121102 and other repeating FRB sources can bring important clues for the physical origin of these sources.

We present results for the first three years of OzDES, a six-year programme to obtain redshifts for objects in the Dark Energy Survey (DES) supernova fields using the 2dF fibre positioner and AAOmega spectrograph on the Anglo-Australian Telescope. OzDES is a multi-object spectroscopic survey targeting multiple types of targets at multiple epochs over a multi-year baseline, and is one of the first multi-object spectroscopic surveys to dynamically include transients into the target list soon after their discovery. At the end of three years, OzDES has spectroscopically confirmed almost 100 supernovae, and has measured redshifts for 17,000 objects, including the redshifts of 2,566 supernova hosts. We examine how our ability to measure redshifts for targets of various types depends on signal-to-noise, magnitude, and exposure time, finding that our redshift success rate increases significantly at a signal-to-noise of 2 to 3 per 1-Angstrom bin. We also find that the change in signal-to-noise with exposure time closely matches the Poisson limit for stacked exposures as long as 10 hours. We use these results to predict the redshift yield of the full OzDES survey, as well as the potential yields of future surveys on other facilities such as the 4m Multi-Object Spectroscopic Telescope (4MOST), the Subaru Prime Focus Spectrograph (PFS), and the Maunakea Spectroscopic Explorer (MSE). This work marks the first OzDES data release, comprising 14,693 redshifts. OzDES is on target to obtain over a yield of approximately 5,700 supernova host-galaxy redshifts.

The jet compositions, central engines, and progenitors of gamma-ray bursts (GRBs) remain open questions in GRB physics. Applying broadband observations, including GRB prompt emission and afterglow properties derived from {\em Fermi} and {\em Swift} data, as well as from Keck host-galaxy observations, we address these questions for the peculiar, bright GRB 110731A. By using the pair-opacity method, we derive $\Gamma_{0}>190$ during the prompt emission phase. Alternatively, we derive $\Gamma_{0} \approx 580$ and $\Gamma_{0} \approx 154$ by invoking the early-afterglow phase within the homogeneous density and wind cases, respectively. On the other hand, nondetection of a thermal component in the spectra suggests that the prompt emission is likely powered by dissipation of a Poynting-flux-dominated jet leading to synchrotron radiation in an optically thin region. The nondetection of a jet break in the X-ray and optical bands allows us to place a lower limit on the jet opening angle $\theta_{\rm j}>5.5^{\circ}$. Within a millisecond magnetar central engine scenario, we derive the period $P_{0}$ and polar magnetic field strength $B_{\rm p}$, which have extreme (but still allowed) values. The moderately short observed duration (7.3\,s) and relatively large redshift ($z=2.83$) places the burst as a "rest-frame short" GRB, so the progenitor of the burst is subject to debate. Its relatively large $f_{{\rm eff}, z}$ parameter (ratio of the 1\,s peak flux of a pseudo-GRB and the background flux) and a large physical offset from a potential host galaxy suggest that the progenitor of GRB 110731A may be a compact-star merger.

The early optical afterglow emission of several gamma-ray bursts (GRBs) shows a high linear polarization degree (PD) of tens of percent, suggesting an ordered magnetic field in the emission region. The light curves are consistent with being of a reverse shock (RS) origin. However, the magnetization parameter, $\sigma$, of the outflow is unknown. If $\sigma$ is too small, an ordered field in the RS may be quickly randomized due to turbulence driven by various perturbations so that the PD may not be as high as observed. Here we use the "Athena++" relativistic MHD code to simulate a relativistic jet with an ordered magnetic field propagating into a clumpy ambient medium, with a focus on how density fluctuations may distort the ordered magnetic field and reduce PD in the RS emission for different $\sigma$ values. For a given density fluctuation, we discover a clear power-law relationship between the relative PD reduction and the $\sigma$ value of the outflow. Such a relation may be applied to estimate $\sigma$ of the GRB outflows using the polarization data of early afterglows.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) are plausible candidates for the central engines of gamma-ray bursts (GRBs). NDAFs are hyperaccretion disks with accretion rates in the range of around 0.001-10 $M_\odot ~\rm s^{-1}$, which have high density and temperature and therefore are extremely optically thick and geometrically slim or even thick. We review the theoretical progresses in studying the properties of NDAFs as well as their applications to the GRB phenomenology. The topics include: the steady radial and vertical structure of NDAFs and the implications for calculating neutrino luminosity and annihilation luminosity, jet power due to neutrino-antineutrino annihilation and Blandford-Znajek mechanism and their dependences on parameters such as BH mass, spin, and accretion rate, time evolution of NDAFs, effect of magnetic fields, applications of NDAF theories to the GRB phenomenology such as lightcurve variability, extended emission, X-ray flares, kilonovae, etc., as well as probing NDAFs using multi-messenger signals such as MeV neutrinos and gravitational waves.

We present a study of the three-dimensional structure of the molecular clouds in the Galactic Centre (GC) using CO emission and OH absorption lines. Two CO isotopologue lines, $^{12}$CO ($J$=1$\rightarrow$0) and $^{13}$CO ($J$=1$\rightarrow$0), and four OH ground-state transitions, surveyed by the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH), contribute to this study. We develop a novel method to calculate the OH column density, excitation temperature, and optical depth precisely using all four OH lines, and we employ it to derive a three-dimensional model for the distribution of molecular clouds in the GC for six slices in Galactic latitude. The angular resolution of the data is 15.5 arcmin, which at the distance of the GC (8.34 kpc) is equivalent to 38 pc. We find that the total mass of OH in the GC is in the range 2400-5100 Solar mass . The face-on view at a Galactic latitude of b = 0{\deg} displays a bar-like structure with an inclination angle of 67.5 $\pm$ 2.1{\deg} with respect to the line of sight. No ring-like structure in the GC is evident in our data, likely due to the low spatial resolution of the CO and OH maps.

Neutron stars may sustain a non-axisymmetric deformation due to magnetic distortion and are potential sources of continuous gravitational waves (GWs) for ground-based interferometric detectors. With decades of searches using available GW detectors, no evidence of a GW signal from any pulsar has been observed. Progressively stringent upper limits of ellipticity have been placed on Galactic pulsars. In this work, we use the ellipticity inferred from the putative millisecond magnetars in short gamma-ray bursts (SGRBs) to estimate their detectability by current and future GW detectors. For $\sim 1$ ms magnetars inferred from the SGRB data, the detection horizon is $\sim 30$ Mpc and $\sim 600$ Mpc for advanced LIGO (aLIGO) and Einstein Telescope (ET), respectively. Using the ellipticity of SGRB millisecond magnetars as calibration, we estimate the ellipticity and gravitational wave strain of Galactic pulsars and magnetars assuming that the ellipticity is magnetic-distortion-induced. We find that the results are consistent with the null detection results of Galactic pulsars and magnetars with the aLIGO O1. We further predict that the GW signals from these pulsars/magnetars may not be detectable by the currently designed aLIGO detector. The ET detector may be able to detect some relatively low frequency signals ($<50$ Hz) from some of these pulsars. Limited by its design sensitivity, the eLISA detector seems not suitable for detecting the signals from Galactic pulsars and magnetars.

[Zhang 2014] propose to redefine the true GRB central engine activity duration, $t_{\rm burst}$, by considering the contributions from the prompt $\gamma$-ray emission, X-ray flare and internal plateau features. With a comprehensive study on a large sample of Swift GRBs, it is shown that the $t_{\rm burst}$ distribution in the observer frame consists of a bimodal feature, suggesting the existence of a new population of ultra-long GRBs. In this work, we make a series of further studies on $t_{\rm burst}$: we update the Swift GRB sample up to June 2016; we investigate the properties of $t_{\rm burst}$ distribution in the rest frame; we redefine $t_{\rm burst}$ by involving external plateau contributions; we make a multiple parameter analysis to investigate whether the bursts within the ultra-long population being statistically different in sense of other features besides the duration distribution. We find that for all situations, the distribution of $t_{\rm burst}$ requires two normal distributions in logarithmic space to provide a good fit, both in observer frame and rest frame. Considering the observational gap effect would not completely erase the bimodal distribution feature. However the bursts within the ultra-long population may have no statistically different in sense of other features besides the duration term. We thus suggest that if the ultra-long population of GRBs indeed exists, their central engine mechanism and radiation mechanism should be similar to the normal population, but they have longer central engine activity timescale.

The merger of a double neutron star (NS-NS) binary may result in a rapidly rotating massive NS with an extremely strong magnetic field (i.e., a millisecond magnetar). In this case, the magnetic spin-down of the NS remnant provides an additional source of sustained energy injection, which would continuously power the merger ejecta. The thermal emission from the merger ejecta would give rise to a bright optical "magnetar-powered merger-nova". In this work, we carry out a complete search for magnetar-powered merger-nova from Swift short gamma-ray burst (SGRB) sample. We focus on short GRBs with extended emission or internal plateau, which may signify the presence of magnetars as the central engine. We eventually find three candidates of magnetar-powered merger-nova from the late observations of GRB 050724, GRB 070714B and GRB 061006. With standard parameter values, the magnetar remnant scenario could well interpret the multi-band data of all three bursts, including the extended emission and their late chromatic features in the optical and X-ray data. The peak luminosities of these merger-novae reach several times $10^{42}~{\rm erg~s^{-1}}$, more than one order of magnitude brighter than the traditional "kilo-novae" with peak luminosity of $\sim 10^{41}~{\rm erg~s^{-1}}$. Intense, multi-color late time observations of short GRBs are encouraged to identify more merger-novae in the future.

We carry out a systematical study of the spectral lag properties of 50 single-pulsed Gamma-Ray Bursts (GRBs) detected by Fermi/GBM. By dividing the light curves into multiple consecutive energy channels we provide a new measurement of the spectral lag which is independent on energy channel selections. We perform a detailed statistical study of our new measurements. We find two similar power-law energy dependencies of both the pulse arrival time and pulse width. Our new results on the power-law indices would favor the relativistic geometric effects for the origin of spectral lag. However, a complete theoretical framework that can fully account for the diverse energy dependencies of both arrival time and pulse width revealed in this work is still missing. We also study the spectral evolution behaviors of the GRB pulses. We find that the GRB pulse with negligible spectral lag would usually have a shorter pulse duration and would appear to have a "hardness-intensity tracking" (HIT) behavior and the GRB pulse with a significant spectral lag would usually have a longer pulse duration and would appear to have a "hard-to-soft" (HTS) behavior.

Presently a ${>}3\sigma$ tension exists between values of the Hubble constant $H_0$ derived from analysis of fluctuations in the Cosmic Microwave Background by Planck, and local measurements of the expansion using calibrators of type Ia supernovae (SNe Ia). We perform a blinded reanalysis of Riess et al. 2011 to measure $H_0$ from low-redshift SNe Ia, calibrated by Cepheid variables and geometric distances including to NGC 4258. This paper is a demonstration of techniques to be applied to the Riess et at. 2016 data. Our end-to-end analysis starts from available CfA3 and LOSS photometry, providing an independent validation of Riess et al. 2011. We obscure the value of $H_0$ throughout our analysis and the first stage of the referee process, because calibration of SNe Ia requires a series of often subtle choices, and the potential for results to be affected by human bias is significant. Our analysis departs from that of Riess et al. 2011 by incorporating the covariance matrix method adopted in SNLS and JLA to quantify SN Ia systematics, and by including a simultaneous fit of all SN Ia and Cepheid data. We find $H_0 = 72.5 \pm 3.1$ (stat) $\pm 0.77$ (sys) km s$^{-1}$ Mpc$^{-1}$ with a three-galaxy (NGC 4258+LMC+MW) anchor. The relative uncertainties are 4.3% statistical, 1.1% systematic, and 4.4% total, larger than in Riess et al. 2011 (3.3% total) and the Efstathiou 2014 reanalysis (3.4% total). Our error budget for $H_0$ is dominated by statistical errors due to the small size of the supernova sample, whilst the systematic contribution is dominated by variation in the Cepheid fits, and for the SNe Ia, uncertainties in the host galaxy mass dependence and Malmquist bias.

The BeSSeL Survey is mapping the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with high-mass star formation. While parallax techniques for water masers at high frequency (22 GHz) have been well documented, recent observations of methanol masers at lower frequency (6.7 GHz) have revealed astrometric issues associated with signal propagation through the ionosphere that could significantly limit parallax accuracy. These problems displayed as a "parallax gradient" on the sky when measured against different background quasars. We present an analysis method in which we generate position data relative to an "artificial quasar" at the target maser position at each epoch. Fitting parallax to these data can significantly mitigate the problems and improve parallax accuracy.

The cosmological evolution of primordial black holes (PBHs) is considered. A comprehensive view of the accretion and evaporation histories of PBHs across the entire cosmic history is presented, with focus on the critical mass holes. The critical mass of a PBH for current era evaporation is $M_{cr}\sim 5.1\times10^{14}$ g. Across cosmic time such a black hole will not accrete radiation or matter in sufficient quantity to hasten the inevitable evaporation, if the black hole remains within an average volume of the universe. The accretion rate onto PBHs is most sensitive to the mass of the hole, the sound speed in the cosmological fluid, and the energy density of the accreted components. It is not easy for a PBH to accrete the average cosmological fluid to reach $30M_\odot$ by $z\sim0.1$, the approximate mass and redshift of the merging BHs that were the sources of the gravitational wave events GW150914 and GW151226. A PBH located in an overdense region can undergo enhanced accretion leading to the possibility of growing by many orders of magnitude across cosmic history. Thus, two merging PBHs are a plausible source for the observed gravitational wave events. However, it is difficult for isolated PBHs to grow to supermassive black holes (SMBHs) at high redshift with masses large enough to fit observational constraints.

A hypothetical photon mass, $m_\gamma$, gives an energy-dependent light speed in a Lorentz-invariant theory. Such a modification causes an additional time delay between photons of different energies when they travel through a fixed distance. Fast radio bursts (FRBs), with their short time duration and cosmological propagation distance, are excellent astrophysical objects to constrain $m_\gamma$. Here for the first time we develop a Bayesian framework to study this problem with a catalog of FRBs. Those FRBs with and without redshift measurement are both useful in this framework, and can be combined in a Bayesian way. A catalog of 21 FRBs (including 20 FRBs without redshift measurement, and one, FRB 121102, with a measured redshift $z=0.19273 \pm 0.00008$) give a combined limit $m_\gamma \leq 8.7 \times 10^{-51}\, {\rm kg}$, or equivalently $m_\gamma \leq 4.9 \times 10^{-15}\, {\rm eV}/c^2$ ($m_\gamma \leq 1.5\times10^{-50} \, {\rm kg}$, or equivalently $m_\gamma \leq 8.4 \times 10^{-15} \,{\rm eV}/c^2$) at 68% (95%) confidence level, which represents the best limit that comes purely from kinematics. The framework proposed here will be valuable when FRBs are observed daily in the future. Increment in the number of FRBs, and refinement in the knowledge about the electron distributions in the Milky Way, the host galaxies of FRBs, and the intergalactic median, will further tighten the constraint.

The repeating fast radio burst (FRB) source that produced FRB 121102 was recently localized in a star forming galaxy at $z=0.193$, which is associated with an extended radio source at the burst location. One possibility is that the repeating FRBs are produced by a new-born magnetar, which also powers the radio nebula. If so, the magnetar may produce $\gamma$-ray emission due to magnetic dipolar spin-down. The luminosity depends on the magnetar spin parameters and age. We process the eight-year Fermi LAT data at the position of FRB 121102 and place an energy flux upper limit of $\sim 10^{-11} \ {\rm erg \ cm^{-2} \ s^{-1}}$ in time bins with six-month intervals, and an accumulated energy flux upper limit of $\sim 4\times 10^{-12} \ {\rm erg \ cm^{-2} \ s^{-1}}$ over the eight year span. The corresponding $\gamma$-ray luminosity upper limits are $\sim 10^{45} \ {\rm erg \ s^{-1}}$ and $\sim 4\times 10^{44} \ {\rm erg \ s^{-1}}$ for the time-resolved and time-integrated analyses, respectively. We discuss the implications of these limits on the young magnetar model.

The repeating fast radio burst (FRB) source that produced FRB 121102 was recently localized in a star forming galaxy at $z=0.193$, which is associated with an extended radio source at the burst location. One possibility is that the repeating FRBs are produced by a new-born magnetar, which also powers the radio nebula. If so, the magnetar may produce $\gamma$-ray emission due to magnetic dipolar spin-down. The luminosity depends on the magnetar spin parameters and age. We process the eight-year Fermi LAT data at the position of FRB 121102 and place an energy flux upper limit of $\sim 10^{-11} \ {\rm erg \ cm^{-2} \ s^{-1}}$ in time bins with six-month intervals, and an accumulated energy flux upper limit of $\sim 4\times 10^{-12} \ {\rm erg \ cm^{-2} \ s^{-1}}$ over the eight year span. The corresponding $\gamma$-ray luminosity upper limits are $\sim 10^{45} \ {\rm erg \ s^{-1}}$ and $\sim 4\times 10^{44} \ {\rm erg \ s^{-1}}$ for the time-resolved and time-integrated analyses, respectively. We discuss the implications of these limits on the young magnetar model.

Violations of Lorentz invariance can lead to an energy-dependent vacuum dispersion of light, which results in arrival-time differences of photons arising with different energies from a given transient source. In this work, direction-dependent dispersion constraints are obtained on nonbirefringent Lorentz-violating effects, using the observed spectral lags of the gamma-ray burst GRB 160625B. This burst has unusually large high-energy photon statistics, so we can obtain constraints from the true spectral time lags of bunches of high-energy photons rather than from the rough time lag of a single highest-energy photon. Also, GRB 160625B is the only burst to date having a well-defined transition from positive lags to negative lags, which provides a unique opportunity to distinguish Lorentz-violating effects from any source-intrinsic time lag in the emission of photons of different energy bands. Our results place comparatively robust two-sided constraints on a variety of isotropic and anisotropic coefficients for Lorentz violation, including first bounds on Lorentz-violating effects from operators of mass dimension ten in the photon sector.

Pulsars are magnetized rotating compact objects. They spin down due to magnetic dipole radiation and wind emission. If photon has a nonzero mass, the spin down rate would be smaller than the zero mass case. We show that an upper limit of the photon mass, i.e. $m_\gamma\lesssim h/Pc^2$, may be placed if a pulsar with period $P$ is observed to spin down. Recently, a white dwarf (WD) --- M dwarf binary, AR Scorpii was discovered to emit pulsed broadband emission with pulses. The spin-down luminosity of the WD can comfortably power the non-thermal radiation from the system. Applying our results to the WD pulsar with $P=117~\rm{s}$, we obtain a stringent upper limit of the photon mass between $m_\gamma<6.3\times10^{-50}~\rm{g}$ assuming a vacuum dipole spindown, and $m_\gamma<9.6\times10^{-50}~\rm{g}$ assuming a spindown due to a fully developed pulsar wind.

Multi-color photometric light curves of NO Cam in $V$ $R_C$ and $I_C$ bands are obtained and analyzed simultaneously using the Wilson-Devinney (W-D) program. The solutions suggest that NO Cam is an A-subtype overcontact binary with its mass ratio to be $q = 0.439$ and a contact degree of $f = 55.5\,\%$. The small temperature difference ($\Delta T = 44K$) between its two components indicates that the system is under thermal contact. The high orbital inclination ($i = 84.5^{\circ}$) strengthens our confidence in the parameters determined from the light curves. All available times of minimum light are collected and period variations are analyzed for the first time. The $O - C$ curve reveals that its period is increasing continuously at a rate of $dP/dt=+1.46\times{10^{-9}}$, which can be explained by mass transfer from the less massive component to the more massive one. After the upward parabolic variation is subtracted, the residuals suggest that there may be a cyclic variation with a period of 2.23 years and an amplitude of $A_3 = 0.00153$ days, which may due to the light-travel time effect (LTTE) arising from the gravitational influence of a close-in tertiary component. The close-in companion reveals that early dynamic interaction among triple system may have played very important role in the formation of the W UMa type binaries.

Zhang (2013) proposed a type of GRB-less X-ray transient associated with double neutron star (NS-NS) mergers under the conjecture of a rapidly-spinning magnetar merger product with the line of sight off the short GRB jet. We investigate possible light curves of these transients by considering different observer's viewing angles. We perform Monte Carlo simulations to calculate the peak luminosity function (LF) and event rate density of these X-ray transients. By considering that a fraction of massive neutron stars may be supra-massive and later collapse into black holes after spinning down, we investigate how the predicted LF depends on the equation of state (EoS) of the central object and the geometry of the system. In general, the LF can be fit by two log-normal distributions peaking around $10^{46.4}$ and $10^{49.6}$ $\rm erg\,s^{-1}$, corresponding to the trapped and free zones, respectively. For the majority of the EoS models, the current non-detection is consistent with having a free zone solid angle at most a few times of the solid angle of the short GRB jet. The event rate density of these X-ray transients is around a few tens of $\rm Gpc^{-3}yr^{-1}$ for luminosity above $10^{45}$ $\rm erg\,s^{-1}$. We predict that future X-ray telescopes (such as Einstein Probe) with sensitivity $\sim 10^{-11}$ $\rm erg\,s^{-1}\,cm^{-2}$ would detect as many as several tens of such transients per year per steradian. Within 200 Mpc, the aLIGO average range for NS-NS mergers, the estimated event rate of these transients is about 1 transient per year all sky.

Fast radio bursts (FRBs) have excessive dispersion measures (DMs) and an all-sky distribution, which point toward an extragalactic or even a cosmological origin. We develop a method to extract the mean host galaxy DM ($\left\langle{\rm DM_{HG,loc}}\right\rangle$) and the characterized luminosity ($L$) of FRBs using the observed DM-Flux data, based on the assumption of a narrow luminosity distribution. Applying Bayesian inference to the data of 21 FRBs, we derive a relatively large mean host DM, i.e. $\left\langle{\rm DM_{HG,loc}}\right\rangle \sim 270~\rm{pc~cm^{-3}}$ with a large dispersion. A relatively large ${\rm DM_{HG}}$ of FRBs is also supported by the millisecond scattering times of some FRBs and the relatively small redshift $z=0.19273$ of FRB 121102 (which gives ${\rm DM_{HG,loc}} \sim 210 ~\rm{pc~cm^{-3}}$). The large host galaxy DM may be contributed by the ISM or a near-source plasma in the host galaxy. If it is contributed by the ISM, the type of the FRB host galaxies would not be Milky Way (MW)-like, consistent with the detected host of FRB 121102. We also discuss the possibility of having a near-source supernova remnant (SNR), pulsar wind nebula (PWN) or HII region that gives a significant contribution to the observed ${\rm DM_{HG}}$.

IceCube has detected a cumulative flux of PeV neutrinos, which origin is unknown. Blazars, active galactic nuclei with relativistic jets pointing to us, are long and widely expected to be one of the strong candidates of high energy neutrino sources. The neutrino production depends strongly on the cosmic ray power of blazar jets, which is largely unknown. The recent null results in stacking searches of neutrinos for several blazar samples by IceCube put upper limits on the neutrino fluxes from these blazars. Here we compute the cosmic ray power and PeV neutrino flux of Fermi-LAT blazars, and find that the upper limits for known blazar sources give stringent constraint on the cosmic ray loading factor of blazar jets (i.e., the ratio of the cosmic ray to bolometric radiation luminosity of blazar jets), $\xi_{\rm cr}\la(2-10)\zeta^{-1}$ (with $\zeta\la1$ the remained fraction of cosmic ray energy when propagate into the blazar broad line region) for flat cosmic ray spectrum, and that the cumulative PeV neutrino flux contributed by all-sky blazars is a fraction $\la(10-50)\%$ of the IceCube detected flux.

The SkyMapper 1.3 m telescope at Siding Spring Observatory has now begun regular operations. Alongside the Southern Sky Survey, a comprehensive digital survey of the entire southern sky, SkyMapper will carry out a search for supernovae and other transients. The search strategy, covering a total footprint area of ~2000 deg2 with a cadence of $\leq 5$ days, is optimised for discovery and follow-up of low-redshift type Ia supernovae to constrain cosmic expansion and peculiar velocities. We describe the search operations and infrastructure, including a parallelised software pipeline to discover variable objects in difference imaging; simulations of the performance of the survey over its lifetime; public access to discovered transients; and some first results from the Science Verification data.

Recent observations of fast radio bursts (FRBs) indicate a perplexing, inconsistent picture. We propose a unified scenario to interpret diverse FRBs observed. A regular pulsar, otherwise unnoticeable at a cosmological distance, may produce a bright FRB if its magnetosphere is suddenly "combed" by a nearby, strong plasma stream towards the anti-stream direction. If the Earth is to the night side of the stream, the combed magnetic sheath would sweep across the direction of Earth and make a detectable FRB. The stream could be an AGN flare, a GRB or supernova blastwave, a tidal disruption event, or even a stellar flare. Since it is the energy {\em flux} received by the pulsar rather than the luminosity of the stream origin that defines the properties of the FRB, this model predicts a variety of counterparts of FRBs, including a possible connection between FRB 150418 and an AGN flare, a possible connection between FRB 131104 and a weak GRB, a steady radio nebula associated with the repeating FRB 121102, and probably no bright counterparts for some FRBs.

We constrain the X-ray properties of the nearby $(360\,{\rm pc})$, old ($5\,{\rm Myr}$) pulsar B1133+16 with $\sim 100\,{\rm ks}$ effective exposure time by {\it XMM-Newton}. The observed pulsar flux in the 0.2-3 keV energy range is $\sim 10^{-14} \, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$, which results in the recording of $\sim 600$ source counts with the EPIC pn and MOS detectors. The X-ray radiation is dominated by nonthermal radiation and is well described by both a single power-law model (PL) and a sum of blackbody and power-law emission (BB+PL). The BB+PL model results in a spectral photon index $\Gamma=2.4^{+0.4}_{-0.3}$ and a nonthermal flux in the 0.2-3 keV energy range of $(7\pm 2) \times 10^{-15}\, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$. The thermal emission is consistent with the blackbody emission from a small hot spot with a radius of $R_{\rm pc} \approx 14^{+7}_{-5} \, {\rm m}$ and a temperature of $T_{\rm s} = 2.9^{+0.6}_{-0.4} \, {\rm MK}$. Assuming that the hot spot corresponds to the polar cap of the pulsar, we can use the magnetic flux conservation law to estimate the magnetic field at the surface $B_{\rm s} \approx 3.9 \times 10^{14} \, {\rm G}$. The observations are in good agreement with the predictions of the partially screened gap model, which assumes the existence of small-scale surface magnetic field structures in the polar cap region.

Possible violations of Lorentz invariance (LIV) have been investigated for a long time using the observed spectral lags of gamma-ray bursts (GRBs). However, these generally have relied on using a single photon in the highest energy range. Furthermore, the search for LIV lags has been hindered by our ignorance concerning the intrinsic time lag in different energy bands. GRB 160625B, the only burst so far with a well-defined transition from $positive$ lags to $negative$ lags provides a unique opportunity to put new constraints on LIV. Using multi-photon energy bands we consider the contributions to the observed spectral lag from both the intrinsic time lag and the lag by LIV effects, and assuming the intrinsic time lag to have a positive dependence on the photon energy, we obtain robust limits on LIV by directly fitting the spectral lag data of GRB 160625B. Here we show that these robust limits on the quantum gravity energy scales are $E_{\rm QG,1}\geq0.5\times10^{16}$ GeV for the linear, and $E_{\rm QG,2}\geq1.4\times10^{7}$ GeV for the quadratic LIV effects, respectively. In addition, we give for the first time a reasonable formulation of the intrinsic energy-dependent time lag.

A $\gamma$-ray transient, Swift J0644.5-5111, has been claimed to be associated with FRB 131104. The $\gamma$-ray energy output is estimated as $E_\gamma \approx 5\times 10^{51}$\,erg at the nominal $z\approx 0.55$ redshift implied by the dispersion measure of FRB 131104. However, a long-term radio imaging follow-up observations only place an upper limit on the radio afterglow flux of Swift J0644.5-5111. Applying the external shock model, we make a detailed constraint on the afterglow parameters for the FRB 131104/Swift J0644.5-5111 system. We find that for the commonly used microphysics shock parameters (e.g., $\epsilon_e=0.1$, $\epsilon_B=0.01$ and $p=2.3$), if the redshift value inferred from the DM value is correct to order of magnitude (i.e., $z>0.1$), the ambient medium number density should be $\leq 10^{-3}~\rm{cm^{-3}}$, which is the typical value for a compact binary merger environment but disfavors a massive star origin. Assuming a typical ISM density, one would require that the redshift of the FRB much smaller than the value inferred from DM value ($z\ll0.1$), implying a non-cosmological origin of DM. The constraints are much looser if one adopts smaller $\epsilon_B$ and $\epsilon_e$ values, as observed in some GRB afterglows. The FRB 131104/Swift J0644.5-5111 association remains plausible. We critically discuss possible progenitor models for the system.

We work on a GRB sample whose initial Lorentz factors ($\Gamma_0$) are constrained with the afterglow onset method and the jet opening angles ($\theta_{\rm j}$) are determined by the jet break time. We confirm the $\Gamma_0$ - $E_{\gamma,\rm iso}$ correlation by Liang et al. (2010), and the $\Gamma_0$ - $L_{\gamma,\rm iso}$ correlation by L{\"u} et al. (2012). Furthermore, we find correlations between $\Gamma_0$ and the beaming corrected $\gamma$-ray energy ($E_{\gamma}$) and mean $\gamma$-ray luminosity ($L_{\gamma}$). By also including the kinetic energy of the afterglow, we find rough correlations (with larger scatter) between $\Gamma_0$ and the total ($\gamma$-ray plus kinetic) energy and the total mean luminosity, both for isotropic values and beaming corrected values: these correlations allow us to test the data with GRB central engine models. Limiting our sample to the GRBs that likely have a black hole central engine, we compare the data with theoretical predictions of two types of jet launching mechanisms from BHs, i.e. the non-magnetized $\nu \bar\nu$ - annihilation mechanism, and the strongly magnetized Blandford-Znajek (BZ) mechanism. We find that the data are more consistent with the latter mechanism, and discuss the implications of our findings for GRB jet composition.

Observations reveal a uniform Kolmogorov turbulence throughout the diffuse ionized interstellar medium (ISM) and supersonic turbulence preferentially located in the Galactic plane. Correspondingly, we consider the Galactic distribution of electron density fluctuations consisting of not only a Kolmogorov density spectrum but also a short-wave-dominated density spectrum with the density structure formed at small scales due to shocks. The resulting dependence of the scatter broadening time on the dispersion measure (DM) naturally interprets the existing observational data for both low and high-DM pulsars. According to the criteria that we derive for a quantitative determination of scattering regimes over wide ranges of DMs and frequencies $\nu$, we find that the pulsars with low DMs are primarily scattered by the Kolmogorov turbulence, while those at low Galactic latitudes with high DMs undergo more enhanced scattering dominated by the supersonic turbulence, where the corresponding density spectrum has a spectral index $\approx 2.6$. Besides, by considering a volume filling factor of the density structures with the dependence on $\nu$ as $\propto \nu^{1.4}$ in the supersonic turbulence, our model can also explain the observed shallower $\nu$ scaling of the scattering time than the Kolmogorov scaling for the pulsars with relatively large DMs. The comparison between our analytical results and the scattering measurements of pulsars in turn makes a useful probe of the properties of the large-scale ISM turbulence, e.g., an injection scale of $\sim 100$ pc, and also characteristics of small-scale density structures.

Einstein's weak equivalence principle (EEP) can be tested through the arrival time delay between photons with different frequencies. Assuming that the arrival time delay is solely caused by the gravitational potential of the Milky Way, we show that a "nano-shot" giant pulse with a time delay between energies corrected for all known effects, e.g. $\Delta t<0.4~\rm{ns}$, from the Crab pulsar poses a new upper limit on the deviation from EEP, i.e. $\Delta\gamma < (0.6-1.8)\times 10^{-15}$. This result provides the hitherto most stringent constraint on the EEP, improving by at least 2 to 3 orders of magnitude from the previous results based on fast radio bursts.

We summarize basic observational properties of gamma-ray bursts (GRBs), including prompt emission properties, afterglow properties, and classification schemes. We also briefly comment on the current physical understanding of these properties.

We investigated the physical properties of molecular clouds and star formation processes around infrared bubbles which are essentially expanding HII regions. We performed observations of 13 galactic infrared bubble fields containing 18 bubbles. Five molecular lines, 12CO (J=1-0), 13CO (J=1-0), C18O(J=1-0), HCN (J=1-0), and HCO+ (J=1-0), were observed, and several publicly available surveys, GLIMPSE, MIPSGAL, ATLASGAL, BGPS, VGPS, MAGPIS, and NVSS, were used for comparison. We find that these bubbles are generally connected with molecular clouds, most of which are giant. Several bubble regions display velocity gradients and broad shifted profiles, which could be due to the expansion of bubbles. The masses of molecular clouds within bubbles range from 100 to 19,000 solar mass, and their dynamic ages are about 0.3-3.7 Myr, which takes into account the internal turbulence pressure of surrounding molecular clouds. Clumps are found in the vicinity of all 18 bubbles, and molecular clouds near four of these bubbles with larger angular sizes show shell-like morphologies, indicating that either collect-and-collapse or radiation-driven implosion processes may have occurred. Due to the contamination of adjacent molecular clouds, only six bubble regions are appropriate to search for outflows, and we find that four of them have outflow activities. Three bubbles display ultra-compact HII regions at their borders, and one of them is probably responsible for its outflow. In total, only six bubbles show star formation activities in the vicinity, and we suggest that star formation processes might have been triggered.