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

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 periodic activity of repeating fast radio burst (FRB) 180916.J0158+65 was recently reported by the CHIME/FRB Collaboration team. 28 bursts from this source not only show a $\sim16$-day period with an active phase of $\sim 4.0$ days but also exhibit a broken power law in differential energy distribution. In this paper, we suggest that FRB 180916.J0158+65-like periodic FRBs would provide a unique probe of extragalactic asteroid belts (EABs), based on our previously-proposed pulsar-EAB impact model, in which repeating FRBs arise from an old-aged, slowly-spinning, moderately-magnetized pulsar traveling through an EAB around another stellar-mass object. These two objects form a binary and thus the observed period is in fact the orbital period. We show that this model can be used to well interpret all of the observed data of FRB 180916.J0158+65. Furthermore, we constrain the EAB's physical properties and find that (1) the outer radius of the EAB is at least an order of magnitude smaller than that of its analogue in the solar system, (2) the differential size distribution of the EAB's asteroids at small diameters (large diameters) is shallower (steeper) than that of solar-system small objects, and (3) the two belts have a comparable mass.

It is widely believed that magnetars could be born in core-collapse supernovae (SNe), binary neutron star (BNS) or binary white dwarf (BWD) mergers, or accretion-induced collapse (AIC) of white dwarfs. In this paper, we investigate whether magnetars could also be produced from neutron star--white dwarf (NSWD) mergers, motivated by FRB 180924-like fast radio bursts (FRBs) possibly from magnetars born in BNS/BWD/AIC channels suggested by \cite{mar19}. By a preliminary calculation, we find that NSWD mergers with unstable mass transfer could result in the NS acquiring an ultra-strong magnetic field via the dynamo mechanism due to differential rotation and convection or possibly via the magnetic flux conservation scenario of a fossil field. If NSWD mergers can indeed create magnetars, then such objects could produce at least a subset of FRB 180924-like FRBs within the framework of flaring magnetars, since the ejecta, local environments, and host galaxies of the final remnants from NSWD mergers resemble those of BNS/BWD/AIC channels. This NSWD channel is also able to well explain both the observational properties of FRB 180924-like and FRB 180916.J0158+65-like FRBs within a large range in local environments and host galaxies.

We study the dynamics of plasma along the legs of an arch filament system (AFS) from the chromosphere to the photosphere, observed with high-cadence spectroscopic data from two ground-based solar telescopes: the GREGOR telescope (Tenerife) using the GREGOR Infrarred Spectrograph (GRIS) in the He I 10830 \r{A} range and the Swedish Solar Telescope (La Palma) using the CRisp Imaging Spectro-Polarimeter to observe the Ca II 8542 \r{A} and Fe I 6173 \r{A} spectral lines. The temporal evolution of the draining of the plasma was followed along the legs of a single arch filament from the chromosphere to the photosphere. The average Doppler velocities inferred at the upper chromosphere from the He I 10830 \r{A} triplet reach velocities up to 20-24~km~s$^{-1}$, in the lower chromosphere and upper photosphere the Doppler velocities reach up to 11~km~s$^{-1}$ and 1.5~km~s$^{-1}$ in the case of the Ca II 8542 \r{A} and Si I 10827 \r{A} spectral lines, respectively. The evolution of the Doppler velocities at different layers of the solar atmosphere (chromosphere and upper photosphere) shows that they follow the same LOS velocity pattern, which confirm the observational evidence that the plasma drains towards the photosphere as proposed in models of AFSs. The Doppler velocity maps inferred from the lower photospheric Ca I 10839 \r{A} or Fe I 6173 \r{A} spectral lines do not show the same LOS velocity pattern. Thus, there is no evidence that the plasma reaches the lower photosphere. The observations and the nonlinear force-free field extrapolations demonstrate that the magnetic field loops of the AFS rise with time. We found flow asymmetries at different footpoints of the AFS. The NLFFF values of the magnetic field strength give us a clue to explain these flow asymmetries.

Two accretion columns have been argued to form over the surface of a newborn millisecond magnetar for an extremely high accretion rate $\gtrsim1.8\times10^{-2}M_\odot\ {\rm s^{-1}}$ that may occur in the core-collapse of a massive star. In this paper, we investigate the characteristics of these accretion columns and their gravitational wave (GW) radiation. For a typical millisecond magnetar (surface magnetic field strength $B\sim10^{15}$ G and initial spin period $P\sim1$ ms), we find (1) its accretion columns are cooled via neutrinos and can reach a height $\sim1$ km over the stellar surface; (2) its column-induced characteristic GW strain is comparable to the sensitivities of the next generation ground-based GW detectors within a horizon $\sim1$ Mpc; (3) the magnetar can survive only a few tens of seconds; (4) during the survival timescale, the height of the accretion columns increases rapidly to the peak and subsequently decreases slowly; (5) the column mass, characteristic GW strain, and maximum GW luminosity have simultaneous peaks in a similar rise-fall evolution. In addition, we find that the magnetar's spin evolution is dominated by the column accretion torque. A possible association with failed supernova is also discussed.

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

For a binary composed of a spinning black hole (BH) (with mass $\gtrsim 7M_\odot$) and a strongly magnetized neutron star (NS) (with surface magnetic field strength $B_{\rm S,NS}\gtrsim10^{12}$\,G and mass $\sim 1.4M_\odot$), the NS as a whole will possibly eventually plunge into the BH. During the inspiral phase, the spinning BH could be charged to the Wald charge quantity $Q_{\rm W}$ until merger in an electro-vacuum approximation. During the merger, if the spinning charged BH creates its own magnetosphere due to an electric field strong enough for pair cascades to spark, the charged BH would transit from electro-vacuum to force-free cases and could discharge in a time $\gtrsim1~{\rm ms}$. As the force-free magnetosphere is full of a highly conducting plasma, the magnetic flux over the NS's caps would be retained outside the BH's event horizon under the frozen-in condition. Based on this scenario, we here investigate three possible energy-dissipation mechanisms that could produce electromagnetic (EM) counterparts in a time interval of the BH's discharge post a BH-NS merger-induced gravitational wave event: (1) magnetic reconnection at the BH's poles would occur, leading to a millisecond bright EM signal, (2) a magnetic shock in the zone of closed magnetic field lines due to the detachment and reconnection of the entire BH magnetic field would probably produce a bright radio emission, e.g., a fast radio burst, and (3) the Blandford-Znajek mechanism would extract the BH's rotational energy, giving rise to a millisecond-duration luminous high-energy burst. We also calculate the luminosities due to these mechanisms as a function of BH's spin for different values of $B_{\rm S,NS}$.

The emergence of active regions (ARs) closely relates to the solar dynamo and the dynamical atmospheric phenomena.With high-resolution and long-lasting observations from the New Vacuum Solar Telescope, we report a new dynamic activity phenomenon named "fan-shaped activity (FSA)" in the emerging phase of NOAA AR 12700. The FSAs are clearly observed at Ha wavelength and are closely related to the dynamics of the adjacent arch filament system (AFS), including threads deformation and materials downward motions. On 2018 February 26, the two most representative FSAs appeared around 05:21 UT and 06:03 UT, respectively, and they firstly ascended and then decayed in around 10 minutes. At the ascending phase, accompanied by the uplifting of an adjacent AFS, each FSA launches up at one end of the AFS and extends for up to 11 Mm. At the decaying phase, the FSA gradually vanishes, and materials downflows towards the other end of the AFS are detected. After checking the evolution of the magnetic fields of AR 12700, we find that each FSA is located between the end of an AFS and an adjacent magnetic patch with the same polarity and launches at the onset of the collision and compression between these two magnetic patches. We propose that the collision lifts up the AFS, and then the initially compact AFS laterally expands, resulting in the formation of FSA. A cartoon model is proposed to depict the activities.

Gamma-ray burst (GRB) central engines and jet production mechanisms are still open questions. Assuming that the shallow decay segments of canonical X-ray afterglow lightcurves of {\em Swift} GRBs are attributed to the magnetic dipole (MD) radiations of newly-born magnetars, we derive the parameters of the magnetars and explore their possible relations to jet and MD wind emission. We show that the magnetar initial spin period ($P_0$) are tightly correlated with the jet energy ($E_{\rm jet}$), which is almost proportional to the wind energy ($E_{\rm wind}$). Our least square fits yield $P_0\propto E^{-0.36\pm 0.03}_{\rm jet}$ and $E_{\rm wind}\propto E^{0.91\pm 0.07}_{\rm jet}$. These relations may imply that a magnetar with faster rotating speed can power a more energetic GRB, and energy partition between the jet and wind may be quasi-universal. Although the $P_0-E_{\rm jet}$ relation is driven by a few sub-energetic GRBs in our sample, our Monte Carlo simulation analysis shows that sample selection biases from instrumental flux limits and contaminations of the bright jet afterglows cannot make this correlation. Within this jet-wind paradigm, we propose that GRB 101225A-like ultra-long GRBs, whose prompt gamma-ray/X-ray lightcurves are featured as a long-lasting plateau with a sharp drop, may be the orphan MD wind emission being due to misalignment of their jet axis to the light of sight. Brief discussion on the orphan MD wind emission and its association with the gravitational wave radiation of newly-born magnetars is presented.

The emergence of active regions (ARs) leads to various dynamic activities. Using high-resolution and long-lasting H$\alpha$ observations from the New Vacuum Solar Telescope, we report the dynamics of NOAA AR 12700 in its emerging phase on 26 February 2018 in detail. In this AR, constant interchange reconnections between emerging fibrils and preexisting ones were detected. Driven by the flux emergence, small-scale fibrils observed in H$\alpha$ wavelength continuously emerged at the center of the AR and reconnected with the ambient preexisting fibrils, forming new longer fibrils. We investigate three scenarios of such interchange reconnection in two hours. Specially, the third scenario of reconnection resulted in the formation of longer fibrils that show pronounced rotation motion. To derive the evolution of the magnetic structure during the reconnections, we perform nonlinear force-free field extrapolations. The extrapolated three-dimensional magnetic fields clearly depict a set of almost potential emerging loops, two preexisting flux ropes at 03:00 UT before the second reconnection scenario, and a set of newly formed loops with less twist at 03:48 UT after the third reconnection scenario. All of these extrapolated structures are consistent with the fibrils detected in H$\alpha$ wavelength. The aforementioned observations and extrapolation results suggest that the constant interchange reconnections resulted in that the magnetic twist was redistributed from preexisting flux ropes towards the newly-formed system with longer magnetic structure and weaker twist.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

We present our optical photometric and spectroscopical observations of GRB 140629A. A redshift of $z=2.275\pm0.043$ is measured through the metal absorption lines in our spectroscopic data. Using our photometric data and multiple observational data observed with other telescopes, we show that its optical light curve is well interpreted with the standard forward shock models in the thin shell case. Its optical-X-ray afterglow spectrum is jointly fitted with a single power-law function, yielding a photon index of $-1.90\pm 0.05$. The optical extinction and neutral hydrogen absorption of the GRB host galaxy are negligible. The fit to the light curve with the standard models shows that the ambient density is $60\pm 9$ cm$^{-3}$ and the GRB radiating efficiency is as low as $\sim 0.24\%$, likely indicating a baryonic-dominated ejecta of this GRB. This burst agrees well with the $L_{\rm p, iso}-E_p^{'}-\Gamma_0$ relation, but confidently violates those empirical relations involving geometric corrections (or jet break time). This gives rise to an issue of possible selection effect on these relations since the jet opening angle of this GRB is extremely narrow (0.04 rad).

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.

We study the rates of tidal disruption of stars by intermediate-mass to supermassive black holes on bound to unbound orbits by using high-accuracy direct N-body experiments. The approaching stars from the star cluster to the black hole can take three types of orbit: eccentric, parabolic, and hyperbolic orbits. Since the mass fallback rate shows a different variability depending on these orbital types, we can classify tidal disruption events (TDEs) into three main categories: eccentric, parabolic, and hyperbolic TDEs. Respective TDEs are characterized by two critical values of the orbital eccentricity: the lower critical eccentricity is the one below which the stars on eccentric orbits cause the finite, intense accretion, and the higher critical eccentricity above which the stars on hyperbolic orbits cause no accretion. Moreover, we find that the parabolic TDEs are divided into three subclasses: precisely parabolic, marginally eccentric, and marginally hyperbolic TDEs. We analytically derive that the mass fallback rate of the marginally eccentric TDEs can be flatter and slightly higher than the standard fallback rate proportional to $t^{-5/3}$, whereas it can be flatter and lower for the marginally hyperbolic TDEs. We confirm by N-body experiments that only few eccentric, precisely parabolic, and hyperbolic TDEs can occur in a spherical stellar system with a single intermediate-mass to supermassive black hole. A substantial fraction of the stars approaching to the black hole would cause the marginally eccentric or marginally hyperbolic TDEs.

The conservative Post-Newtonian (PN) Hamiltonian formulation of spinning compact binaries has six integrals of motion including the total energy, the total angular momentum and the constant unit lengths of spins. The manifold correction method can effectively eliminate the integration errors accumulation in a long time. In this paper, the accelerated manifold correction method based on graphics processing unit (GPU) is designed to simulate the dynamic evolution of spinning compact binaries. The feasibility and the efficiency of parallel computation on GPU for spinning compact binaries have been confirmed by various numerical experiments. The numerical comparisons show that the accuracy on GPU execution of manifold corrections method has a good agreement with the execution of codes on merely central processing unit (CPU-based) method. The acceleration ability when the codes are implemented on GPU can increase enormously through the use of shared memory and register optimization techniques without additional hardware costs, implying that the speedup is nearly 13 times as compared with the codes executed on CPU for phase space scan (including orbits). In addition, GPU-accelerated manifold correction method is used to numerically study how dynamics are affected by the spin-induced quadrupole-monopole interaction for black hole binary system.

GRB 160625B is an extremely bright GRB with three distinct emission episodes. By analyzing its data observed with the GBM and LAT on board the {\em Fermi} mission, we find that a multi-color black body (mBB) model can be used to fit the spectra of initial short episode (Episode I) very well within the hypothesis of photosphere emission of a fireball model. The time-resolved spectra of its main episode (Episode II), which was detected with both GBM and LAT after a long quiet stage ($\sim 180$ seconds) of the initial episode, can be fitted with a model composing of an mBB component plus a cutoff power-law (CPL) component. This GRB was detected again in the GBM and LAT bands with a long extended emission (Episode III) after a quiet period of $\sim 300$ seconds. The spectrum of Episode III is adequately fitted with a CPL plus a single power-law models, and no mBB component is required. These features may imply that the emission of three episodes are dominated by distinct physics process, i.e., Episode I is possible from cocoon emission surrounding the relativistic jet, Episode II may be from photosphere emission and internal shock of relativistic jet, and Episode III is contributed by internal and external shocks of relativistic jet. On the other hand, both X-ray and optical afterglows are consistent with standard external shocks model.

GRB 160821B is a short gamma-ray burst (GRB) at redshift $z=0.16$, with a duration less than 1 second and without detection of any "extended emission" up to more than 100 seconds in both {\em Swift}/BAT and {\em Fermi}/GBM bands. An X-ray plateau with a sharp drop 180 seconds after the BAT trigger was observed with {\em Swift}/XRT. No supernova or kilo-nova signature was detected. Assuming the central engine of this SGRB is a recently born supra-massive magnetar, we can explain the SGRB as jet radiation and its X-ray plateau as the internal energy dissipation of the pulsar wind as it spins down. We constrain its surface magnetic field as $B_{\rm p}<3.12\times 10^{16}$ G and initial spin period as $P_0< 8.5\times 10^{-3}$ seconds. Its equation of state is consistent with the GM1 model with $M_{\rm TOV} \sim 2.37 M_\odot$ and ellipticity $\epsilon<0.07$. Its gravitational wave (GW) radiation may be detectable with the future Einstein Telescope, but is much weaker than the current detectability limit of advanced-LIGO. The GW radiation of such an event would be detectable by advanced-LIGO if it occurred at a distance of 100 Mpc ($z=0.023$).

With the application of advanced astronomical technologies, equipments and methods all over the world, astronomy covers from radio, infrared, visible light, ultraviolet, X-ray and gamma ray band, and enters into the era of full wavelength astronomy. How to effectively integrate data from different ground- and space-based observation equipments, different observers, different bands, different observation time, requires the data fusion technology. In this paper we introduce the cross-match tool that is developed by the Python language and based on the PostgreSQL database and uses Q3C as the core index, facilitating the cross-match work of massive astronomical data. It provides four different cross-match functions, namely: I) cross-match of custom error range; II) cross-match of catalog error; III) cross-match based on the elliptic error range; IV) cross-match of the nearest algorithm. The cross-match result set provides good foundation for subsequent data mining and statistics based on multiwavelength data. The most advantage of this tool is a user-oriented tool applied locally by users. By means of this tool, users may easily create their own databases, manage their own data and cross-match databases according to their requirements. In addition, this tool is also able to transfer data from one database into another database. More importantly, the tool is easy to get started and used by astronomers without writing any code.

We report our observations of very bright prompt optical and reverse shock (RS) optical emission of GRB 140512A and analyze its multi-wavelength data observed with the {\em Swift} and {\em Fermi} missions. It is found that the joint optical-X-ray-gamma-ray spectrum with our first optical detection (R=13.09 mag) at $T_0+136$ seconds during the second episode of the prompt gamma-rays can be fit by a single power-law with index $-1.32\pm 0.01$. Our empirical fit to the afterglow lightcurves indicates that the observed bright optical afterglow with R=13.00 mag at the peak time is consistent with predictions of the RS and forward shock (FS) emission of external shock models. Joint optical-X-ray afterglow spectrum is well fit with an absorbed single power-law, with an index evolving with time from $-1.86\pm 0.01$ at the peak time to $-1.57\pm 0.01$ at late epoch, which could be due to the evolution of the ratio of the RS to FS emission fluxes. We fit the lightcurves with standard external models, and derive the physical properties of the outflow. It is found that the ratio $R_B\equiv\epsilon_{\rm B, r}/\epsilon_{\rm B, f}$ is 8187, indicating a high magnetization degree in the RS region. Measuring the relative radiation efficiency with $R_e\equiv\epsilon_{\rm e, r}/\epsilon_{\rm e, f}$, we have $R_e= 0.02$, implying the radiation efficiency of the RS is much lower than that in FS. We also show that the $R_B$ of GRBs 990123, 090102, and 130427A are similar to that of GRB 140512A and their apparent difference may be mainly attributed to the difference of the jet kinetic energy, initial Lorentz factor, and medium density among them.

A rapidly spinning, strongly magnetized neutron star is invoked as the central engine for some Gamma-ray bursts (GRBs), especially, the $"$internal plateau$"$ feature of X-ray afterglow. However, for these $"$internal plateau$"$ GRBs, how to produce their prompt emission remains an open question. Two different physical process have been proposed in the literature, (1) a new-born neutron star is surrounded by a hyper-accreting and neutrino cooling disk, the GRB jet can be powered by neutrino annihilation aligning the spin axis; (2) a differentially rotating millisecond pulsar was formed due to different angular velocity between the interior core and outer shell parts of the neutron star, which can power an episodic GRB jet. In this paper, by analyzing the data of one peculiar GRB 070110 (with internal plateau), we try to test which model being favored. By deriving the physical parameters of magnetar with observational data, the parameter regime for initial period ($P_{0\rm }$) and surface polar cap magnetic field ($B_{\rm p}$) of the central NS are $(0.96\sim 1.2 )~\rm ms$ and $(2.4\sim 3.7)\times 10^{14}~\rm G$, respectively. The radiative efficiency of prompt emission is about $\eta_{\gamma} \sim 6\%$. However, the radiative efficiency of internal plateau ($\eta_{\rm X}$) is larger than $31\%$ assuming the $M_{\rm NS}\sim1.4 M_{\odot}$ and $P_{0\rm }\sim1.2 ~\rm ms$. The clear difference between the radiation efficiencies of prompt emission and internal plateau implies that they maybe originated from different components (e.g. prompt emission from the relativistic jet powered by neutrino annihilation, while the internal plateau from the magnetic outflow wind).

We report our very early optical observations of GRB 110530A and investigate its jet properties together with its X-ray afterglow data. A peculiar broad onset bump followed by a plateau is observed in its early R band afterglow lightcurve. The optical data in the other bands and the X-ray data are well consistent with the temporal feature of the R band lightcurve. Our joint spectral fits of the optical and X-ray data show that they are in the same regime, with a photon index of $\sim 1.70$. The optical and X-ray afterglow lightcurves are well fitted with the standard external shock model by considering a delayed energy injection component. Based on our modeling results, we find that the radiative efficiency of the GRB jet is $\sim 1\%$ and the magnetization parameter of the afterglow jet is $<0.04$ with the derived extremely low $\epsilon_B$ (the fraction of shock energy to magnetic field) of $(1.64\pm 0.25)\times 10^{-6}$. These results indicate that the jet may be matter dominated. Discussion on delayed energy injection from accretion of late fall-back material of its pre-supernova star is also presented.

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 sec to about 50 hours after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data observed with {\em Swift} and {\em Fermi} missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features as low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow lightcurves with the standard afterglow model via the Markov Chain Monte Carlo (MCMC) technique shows that $\epsilon_e=(6.9\pm 0.3)\times 10^{-2}$, $\epsilon_B=(7.73\pm 0.62)\times 10^{-6}$, $E_{\rm K}=(6.32\pm 0.86)\times 10^{53}\rm erg$, $n=0.100\pm 0.014$ cm$^{-3}$. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. Achromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The $E_p$ of its time-integrated prompt gamma-ray spectrum is $\sim 26$ KeV. Using the initial Lorentz factor ($\Gamma_0=476^{+225}_{-237}$) derived from our afterglow model fit, it is found that GRB 111228A satisfies the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

Tidal Disruption of stars by supermassive central black holes from dense rotating star clusters is modelled by high-accuracy direct N-body simulation. As in a previous paper on spherical star clusters we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, now according to several classes of orbits which only occur in axisymmetric systems (short axis tube and saucer). Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss-cone in phase space which is raised from the fact that total angular momentum $\bf J$ is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals interesting features. Inside the influence radius of the supermassive black hole the angular distribution of disrupted stars has a conspicuous bimodal structure with a local minimum near the equatorial plane. Outside the influence radius this dependence is weak. We show that the bimodal structure of orbital parameters can be explained by the presence of two families of regular orbits, namely short axis tube and saucer orbits. Also the consequences of our results for the loss cone in axisymmetric galactic nuclei are presented.

We present a study of the waiting time distributions (WTDs) of solar energetic particle (SEP) events observed with the spacecraft $WIND$ and $GOES$. Both the WTDs of solar electron events (SEEs) and solar proton events (SPEs) display a power-law tail $\sim \Delta t^{-\gamma}$. The SEEs display a broken power-law WTD. The power-law index is $\gamma_{1} =$ 0.99 for the short waiting times ($<$70 hours) and $\gamma_{2} =$ 1.92 for large waiting times ($>$100 hours). The break of the WTD of SEEs is probably due to the modulation of the corotating interaction regions (CIRs). The power-law index $\gamma \sim$ 1.82 is derived for the WTD of SPEs that is consistent with the WTD of type II radio bursts, indicating a close relationship between the shock wave and the production of energetic protons. The WTDs of SEP events can be modeled with a non-stationary Poisson process which was proposed to understand the waiting time statistics of solar flares (Wheatland 2000; Aschwanden $\&$ McTiernan 2010). We generalize the method and find that, if the SEP event rate $\lambda = 1/\Delta t$ varies as the time distribution of event rate $f(\lambda) = A \lambda^{-\alpha}exp(-\beta \lambda)$, the time-dependent Poisson distribution can produce a power-law tail WTD $\sim \Delta t^{\alpha - 3}$, where $0 \leq \alpha < 2$.

Tidal Disruption of stars by super massive central black holes from dense star clusters is modeled by high-accuracy direct $N$-body simulation. The time evolution of the stellar tidal disruption rate, the effect of tidal disruption on the stellar density profile and for the first time the detailed origin of tidally disrupted stars are carefully examined and compared with classic papers in the field. Up to 128k particles are used in simulation to model the star cluster around the super massive black hole, we use the particle number and the tidal radius of black hole as free parameters for a scaling analysis. The transition from full to empty loss-cone is analyzed in our data, the tidal disruption rate scales with the particle number $N$ in the expected way for both cases. For the first time in numerical simulations (under certain conditions) we can support the concept of a critical radius of Frank & Rees (1976), which claims that most stars are tidally accreted on highly eccentric orbits originating from regions far outside the tidal radius. Due to the consumption of stars moving on radial orbits, a velocity anisotropy is founded inside the cluster. Finally we make an estimation for the real galactic center based on our simulation results and the scaling analysis.

We present direct astrophysical N-body simulations with up to a few million bodies using our parallel MPI/CUDA code on large GPU clusters in China, Ukraine and Germany, with different kinds of GPU hardware. These clusters are directly linked under the Chinese Academy of Sciences special GPU cluster program in the cooperation of ICCS (International Center for Computational Science). We reach about the half the peak Kepler K20 GPU performance for our phi-GPU code [2], in a real application scenario with individual hierarchically block time-steps with the high (4th, 6th and 8th) order Hermite integration schemes and a real core-halo density structure of the modeled stellar systems. The code and hardware are mainly used to simulate star clusters [23, 24] and galactic nuclei with supermassive black holes [20], in which correlations between distant particles cannot be neglected.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.