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

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, in synergy with the SUDARE survey, is a deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). The observations for the CDFS field comprise about 4.38 deg$^2$ down to $r\sim26$ mag. The total on-sky time spans over four years in this field, distributed over four adjacent sub-fields. In this paper, we use the multi-epoch $r$-band imaging data to measure the variability of the detected objects and search for transients. We perform careful astrometric and photometric calibrations and point spread function (PSF) modeling. A new method, referring to as differential running-average photometry, is proposed to measure the light curves of the detected objects. With the method, the difference of PSFs between different epochs can be reduced, and the background fluctuations are also suppressed. Detailed uncertainty analysis and detrending corrections on the light curves are performed. We visually inspect the light curves to select variable objects, and present some objects with interesting light curves. Further investigation of these objects in combination with multi-band data will be presented in our forthcoming paper.

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

The jet composition and energy dissipation mechanism of Gamma-ray bursts (GRBs) and Blazars are fundamental questions which remain not fully understood. One plausible model is to interpret the $\gamma$-ray emission of GRBs and optical emission of blazars as synchrotron radiation of electrons accelerated from the collision-induced magnetic dissipation regions in Poynting-flux-dominated jets. The polarization observation is an important and independent information to test this model. Based on our recent 3D relativistic MHD simulations of collision-induced magnetic dissipation of magnetically dominated blobs, here we perform calculations of the polarization properties of the emission in the dissipation region and apply the results to model the polarization observational data of GRB prompt emission and blazar optical emission. We show that the same numerical model with different input parameters can reproduce well the observational data of both GRBs and blazars, especially the $90^{\circ}$ polarization angle (PA) change in GRB 100826A and the $180^{\circ}$ PA swing in Blazar 3C279. This supports a unified model for GRB and blazar jets, suggesting that collision-induced magnetic reconnection is a common physical mechanism to power the relativistic jet emission from events with very different black hole masses.

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $\sigma_0=B^2/(4\pi n(m_i+m_e)c^2)$, the magnetization for each species is $\sigma_i \sim \sigma_0$ and $\sigma_e \sim (m_i/m_e) \sigma_0$, respectively. We have studied the magnetically dominated regime by varying $\sigma_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $\gamma_{be} \sim \sigma_e$ and $\gamma_{bi} \sim \sigma_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

The optical radiation and polarization signatures in blazars are known to be highly variable during flaring activities. It is frequently argued that shocks are the main driver of the flaring events. However, the spectral variability modelings generally lack detailed considerations of the self-consistent magnetic field evolution modeling, thus so far the associated optical polarization signatures are poorly understood. We present the first simultaneous modeling of the optical radiation and polarization signatures based on 3D magnetohydrodynamic simulations of relativistic shocks in the blazar emission environment, with the simplest physical assumptions. By comparing the results with observations, we find that shocks in a weakly magnetized environment will largely lead to significant changes in the optical polarization signatures, which are seldom seen in observations. Hence an emission region with relatively strong magnetization is preferred. In such an environment, slow shocks may produce minor flares with either erratic polarization fluctuations or considerable polarization variations, depending on the parameters; fast shocks can produce major flares with smooth PA rotations. In addition, the magnetic fields in both cases are observed to actively revert to the original topology after the shocks. All these features are consistent with observations. Future observations of the radiation and polarization signatures will further constrain the flaring mechanism and the blazar emission environment.

We perform 3D relativistic ideal MHD simulations to study the collisions between high-$\sigma$ (Poynting-flux-dominated) blobs which contain both poloidal and toroidal magnetic field components. This is meant to mimic the interactions inside a highly variable Poynting-flux-dominated jet. We discover a significant electromagnetic field (EMF) energy dissipation with an Alfv\'enic rate with the efficiency around 35\%. Detailed analyses show that this dissipation is mostly facilitated by the collision-induced magnetic reconnection. Additional resolution and parameter studies show a robust result that the relative EMF energy dissipation efficiency is nearly independent of the numerical resolution or most physical parameters in the relevant parameter range. The reconnection outflows in our simulation can potentially form the multi-orientation relativistic mini-jets as needed for several analytical models. We also find a linear relationship between the $\sigma$ values before and after the major EMF energy dissipation process. Our results give support to the proposed astrophysical models that invoke significant magnetic energy dissipation in Poynting-flux-dominated jets, such as the internal collision-induced magnetic reconnection and turbulence (ICMART) model for GRBs, and reconnection triggered mini-jets model for AGNs.

Recent observations by the Fermi satellite suggest that a photosphere emission component is contributing to the observed spectrum of many GRBs. One important question is whether the photosphere component can interpret the typical ``Band'' function of GRBs with a typical low energy photon spectral index $\alpha \sim -1$. We perform a detailed study of the photosphere emission spectrum by progressively introducing several physical ingredients previously not fully incorporated, including the probability distribution of the location of a dynamically evolving photosphere, superposition of emission from an equal-arrival-time ``volume'' in a continuous wind, the evolution of optical depth of a wind with finite but evolving outer boundary, as well as the effect of different top-hat wind luminosity ($L_w$) profiles. By assuming a co-moving blackbody spectrum emerging from the photosphere, we find that for an outflow with a constant or increasing $L_w$, the low-energy spectrum below the peak energy ($E_{p}$), can be modified to $F_\nu \sim \nu^{1.5}$ ($\alpha \sim +0.5$). A softer ($-1<\alpha<+0.5$) or flat ($\alpha=-1$) spectrum can be obtained during the $L_w$ decreasing phase or high-latitude-emission-dominated phase. We also study the evolution of $E_{p}$ as a function of wind and photosphere luminosity in this photosphere model. An $E_p-L$ tracking pattern can be reproduced if a certain positive dependence between the dimensionless entropy $\eta$ and $L_w$ is introduced. However, the hard-to-soft evolution pattern cannot be reproduced unless a contrived condition is invoked. In order to interpret the Band spectrum, a more complicated photosphere model or a different energy dissipation and radiation mechanism are needed.

If a small fraction of Fast Radio Bursts (FRBs) are associated with Gamma-Ray Bursts (GRBs), as recently suggested by Zhang, the combination of redshift measurements of GRBs and dispersion measure (DM) measurements of FRBs opens a new window to study cosmology. At $z<2$ where the universe is essentially fully ionized, detections of FRB/GRB pairs can give an independent measurement of the intergalactic medium portion of the baryon mass fraction, $\Omega_b f_{\rm IGM}$, of the universe. If a good sample of FRB/GRB associations are discovered at higher redshifts, the free electron column density history can be mapped, which can be used to probe the reionization history of both hydrogen and helium in the universe. We apply our formulation to GRBs 101011A and 100704A that each might have an associated FRB, and constrained $\Omega_b f_{\rm IGM}$ to be consistent with the value derived from other methods. The methodology developed here is also applicable, if the redshifts of FRBs not associated with GRBs can be measured by other means.

Multiple rebrightenings have been observed in the multiband afterglow of GRB 030329. Especially, a marked and quick rebrightening occurred at about t ~ 1.2 * 10^5 s. Energy injection from late and slow shells seems to be the best interpretation for these rebrightenings. Usually it is assumed that the energy is injected into the whole external shock. However, in the case of GRB 030329, the rebrightenings are so quick that the usual consideration fails to give a satisfactory fit to the observed light curves. Actually, since these late/slow shells coast freely in the wake of the external shock, they should be cold and may not expand laterally. The energy injection then should only occur at the central region of the external shock. Considering this effect, we numerically re-fit the quick rebrightenings observed in GRB 030329. By doing this, we were able to derive the beaming angle of the energy injection process. Our result, with a relative residual of only 5% - 10% during the major rebrightening, is better than any previous modeling. The derived energy injection angle is about 0.035. We assume that these late shells are ejected by the central engine via the same mechanism as those early shells that produce the prompt gamma-ray burst. The main difference is that their velocities are much slower, so that they catch up with the external shock very lately and manifest as the observed quick rebrightenings. If this were true, then the derived energy injection angle can give a good measure of the beaming angle of the prompt gamma-ray emission. Our study may hopefully provide a novel method to measure the beaming angle of gamma-ray bursts.

We report studies of ultra-high energy cosmic ray composition via analysis of depth of airshower maximum (Xmax), for airshower events collected by the High Resolution Fly's Eye (HiRes) observatory. The HiRes data are consistent with a constant elongation rate d<Xmax>/d(log(E)) of 47.9 +- 6.0 (stat.) +- 3.2 (syst.)g/cm^2/decade for energies between 1.6 EeV and 63 EeV, and are consistent with a predominantly protonic composition of cosmic rays when interpreted via the QGSJET01 and QGSJET-II high-energy hadronic interaction models. These measurements constrain models in which the galactic-to-extragalactic transition is the cause of the energy spectrum "ankle' at 4 EeV.

Stereo data collected by the HiRes experiment over a six year period are examined for large-scale anisotropy related to the inhomogeneous distribution of matter in the nearby Universe. We consider the generic case of small cosmic-ray deflections and a large number of sources tracing the matter distribution. In this matter tracer model the expected cosmic ray flux depends essentially on a single free parameter, the typical deflection angle theta. We find that the HiRes data with threshold energies of 40 EeV and 57 EeV are incompatible with the matter tracer model at a 95% confidence level unless theta is larger than 10 degrees and are compatible with an isotropic flux. The data set above 10 EeV is compatible with both the matter tracer model and an isotropic flux.

The High Resolution Fly's Eye experiment has measured the flux of ultrahigh energy cosmic rays using the stereoscopic air fluorescence technique. The HiRes experiment consists of two detectors that observe cosmic ray showers via the fluorescence light they emit. HiRes data can be analyzed in monocular mode, where each detector is treated separately, or in stereoscopic mode where they are considered together. Using the monocular mode the HiRes collaboration measured the cosmic ray spectrum and made the first observation of the Greisen-Zatsepin-Kuzmin cutoff. In this paper we present the cosmic ray spectrum measured by the stereoscopic technique. Good agreement is found with the monocular spectrum in all details.

We investigate the time average mean square displacement $\overline{\delta^2}(x(t))=\int_0^{t-\Delta}[x(t^\prime+\Delta)-x(t^\prime)]^2 dt^\prime/(t-\Delta)$ for fractional Brownian and Langevin motion. Unlike the previously investigated continuous time random walk model $\overline{\delta^2}$ converges to the ensemble average $<x^2 > \sim t^{2 H}$ in the long measurement time limit. The convergence to ergodic behavior is however slow, and surprisingly the Hurst exponent $H=3/4$ marks the critical point of the speed of convergence. When $H<3/4$, the ergodicity breaking parameter ${EB} = {Var} (\overline{\delta^2}) / < \overline{\delta^2} >^2\sim k(H) \cdot\Delta\cdot t^{-1}$, when $H=3/4$, ${EB} \sim (9/16)(\ln t) \cdot\Delta \cdot t^{-1}$, and when $3/4<H <1, {EB} \sim k(H)\Delta^{4-4H} t^{4H-4}$. In the ballistic limit $H \to 1$ ergodicity is broken and ${EB} \sim 2$. The critical point $H=3/4$ is marked by the divergence of the coefficient $k(H)$. Fractional Brownian motion as a model for recent experiments of sub-diffusion of mRNA in the cell is briefly discussed and comparison with the continuous time random walk model is made.

We have searched for correlations between the pointing directions of ultrahigh energy cosmic rays observed by the High Resolution Fly's Eye experiment and Active Galactic Nuclei (AGN) visible from its northern hemisphere location. No correlations, other than random correlations, have been found. We report our results using search parameters prescribed by the Pierre Auger collaboration. Using these parameters, the Auger collaboration concludes that a positive correlation exists for sources visible to their southern hemisphere location. We also describe results using two methods for determining the chance probability of correlations: one in which a hypothesis is formed from scanning one half of the data and tested on the second half, and another which involves a scan over the entire data set. The most significant correlation found occurred with a chance probability of 24%.

Air-fluorescence detectors such as the High Resolution Fly's Eye (HiRes) detector are very sensitive to upward-going, Earth-skimming ultrahigh energy electron-neutrino-induced showers. This is due to the relatively large interaction cross sections of these high-energy neutrinos and to the Landau-Pomeranchuk-Migdal (LPM) effect. The LPM effect causes a significant decrease in the cross sections for bremsstrahlung and pair production, allowing charged-current electron-neutrino-induced showers occurring deep in the Earth's crust to be detectable as they exit the Earth into the atmosphere. A search for upward-going neutrino-induced showers in the HiRes-II monocular dataset has yielded a null result. From an LPM calculation of the energy spectrum of charged particles as a function of primary energy and depth for electron-induced showers in rock, we calculate the shape of the resulting profile of these showers in air. We describe a full detector Monte Carlo simulation to determine the detector response to upward-going electron-neutrino-induced cascades and present an upper limit on the flux of electron-neutrinos.

We report results based on the monitoring of the BL Lac object Mrk 501 in the optical (B, V and R) passbands from March to May 2000. Observations spread over 12 nights were carried out using 1.2 meter Mount Abu Telescope, India and 61 cm Telescope at Sobaeksan Astronomy Observatory, South Korea. The aim is to study the intra-day variability (IDV), short term variability and color variability in the low state of the source. We have detected flux variation of 0.05 mag in the R-band in time scale of 15 min in one night. In the B and V passbands, we have less data points and it is difficult to infer any IDVs. Short term flux variations are also observed in the V and R bands during the observing run. No significant variation in color (B$-$R) has been detected but (V$-$R) shows variation during the present observing run. Assuming the shortest observed time scale of variability (15 min) to represent the disk instability or pulsation at a distance of 5 Schwarschild radii from the black hole (BH), mass of the central BH is estimated $\sim$ 1.20 $\times$ 10$^{8} M_{\odot}$.