GWAC will have been built an integrated FOV of 5,000 $degree^2$ and have already built 1,800 square $degree^2$. The limit magnitude of a 10-second exposure image in the moonless night is 16R. In each observation night, GWAC produces about 0.7TB of raw data, and the data processing pipeline generates millions of single frame alerts. We describe the GWAC Data Processing and Management System (GPMS), including hardware architecture, database, detection-filtering-validation of transient candidates, data archiving, and user interfaces for the check of transient and the monitor of the system. GPMS combines general technology and software in astronomy and computer field, and use some advanced technologies such as deep learning. Practical results show that GPMS can fully meet the scientific data processing requirement of GWAC. It can online accomplish the detection, filtering and validation of millions of transient candidates, and feedback the final results to the astronomer in real-time. During the observation from October of 2018 to December of 2019, we have already found 102 transients.

The outburst and spin-down glitch of magnetars are modeled from the magnetospheric point of view. We try to answer the following four questions: (1) Which pulsar on the period and peirod-derivative diagram are more likely to show magnetar outburst? (2) Which outburst will make the glitch that triggered the outburst to become a spin-down glitch? (3) Can we model the outburst and spin-down glitch in PSR J1119$-$6127 simultaneously? (4) Why the torque variation is delayed compared with the peak of the X-ray luminosity in 1E 1048.1$-$5937 and PSR J1119$-$6127? It is found that both the global and local twisted magnetic field will affect the radiation and timing behaviors of magnetars. Especially, the delay of torque variations may due to the combined effect of increasing twist in the j-bundle and untwisting of the global magnetosphere. A toy model is build for magnetar outburst and torque variations. It can catch the general trend of magnetar outburst: decaying flux, shrinking hot spot, and torque variations.

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.

In this paper, we investigate the magnetohydrodynamical structure of a jet powered by a spinning black hole, where electromagnetic fields and fluid motion are governed by the Grad-Shafranov equation and the Bernoulli equation, respectively. Assuming steady and axisymmetric jet structure, the global solution is uniquely determined with prescribed plasma loading into the jet and the poloidal shape of the outmost magnetic field line. We apply this model to the jet in the center of nearby radio galaxy M87, and we find it can naturally explain the slow flow acceleration and the flow velocity stratification within $10^5$ gravitational radii from the central black hole. In particular, we find the extremal black hole spin is disfavored by the flow velocity measurements, if the plasma loading to the jet is dominated by the electron/positron pair production at the jet base.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

It may be determined by non-parametric method if the dark energy evolves with time. We propose a method of combining PCA and biased estimation on the basis of ridge regression analysis to reconstruct parameters, meanwhile we present an interesting principal component selection criterion to avoid the arbitrariness of principal component selections, and use numerical integral by Lagrange interpolation to linearize the luminosity distance integral formula in nearly flat space to avoid instability of derivative for functional data. We get the preliminary test results that shows if $\Delta \overline w (z) = \overline {\left| {1 + w(z)} \right|} <=0.05$ included $w(z) = -1$, the probability of making a type I error for ${w_{recon}} \ne -1$ is almost zero ($1\% $) in the test; otherwise, if $\Delta \overline w (z) = \overline {\left| {1 + w(z)} \right|} > 0.05$, the probability of making a type I error for ${w_{recon}} = -1$ is not more than $10\% $. Finally, we use JLA sample to reconstruct $w(z)$, and the results reject ${\rm{w(z) }} \ne {\rm{ - 1}}$, which is agreement with $\Lambda CDM$ model.

The sources behind numerous key diffuse interstellar bands (DIBs) remain elusive, and thus evidence is presented here that a family of seven anionic hydrogen atomic clusters (H(-)2n+1, n = 1 to 7) are pertinent contributors. Configuration interaction calculations show hat these charge-induced dipole clusters are stable at temperatures characteristic of the interstellar medium, and emerge from the most abundant element in the Universe. The clusters' spectra yield 25 absorption optical lines that align with observed DIBs to within the computational uncertainties. The absorption bands are due to excitations from the ground states of the clusters to metastable states.

The inconsistent Hubble constant values derived from cosmic microwave background (CMB) observations and from local distance-ladder measurements may suggest new physics beyond the standard $\Lambda$CDM paradigm. It has been found in earlier works that, at least phenomenologically, non-standard recombination histories can reduce the $\gtrsim 4\sigma$ Hubble tension to $\sim 2\sigma$. Following this path, we vary physical and phenomenological parameters in RECFAST, the standard code to compute ionization history of the universe, to explore possible physics beyond standard recombination. We find that the CMB constraint on the Hubble constant is sensitive to the Hydrogen ionization energy and $2s \rightarrow 1s$ two-photon decay rate, both of which are atomic constants, and is insensitive to other details of recombination. Thus, the Hubble tension is very robust against perturbations of recombination history, unless exotic physics modifies the atomic constants during the recombination epoch.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

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.

In this paper, we construct a framework for investigating magnetohydrodynamical jet structure of spinning black holes (BHs), where electromagnetic fields and fluid motion are governed by the Grad-Shafranov equation and the Bernoulli equation, respectively. Assuming steady and axisymmetric jet structure, we can self-consistently obtain electromagnetic fields, fluid energy density and velocity within the jet, given proper plasma loading and boundary conditions. Specifically, we structure the two coupled governing equations as two eigenvalue problems, and develop full numerical techniques for solving them. As an example, we explicitly solve the governing equations for the split monopole magnetic field configuration and simplified plasma loading on the stagnation surface where the poloidal fluid velocity vanishes. As expected, we find the rotation of magnetic field lines is dragged down by fluid inertia, and the fluid as a whole does not contribute to energy extraction from the central BH, i.e., the magnetic Penrose process is not working. However, if we decompose the charged fluid as two oppositely charged components, we find the magnetic Penrose process does work for one of the two components when the plasma loading is low enough.

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

Motivated by the possible existence of other universes, this paper considers the evolution of massive stars with different values for the fundamental constants. We focus on variations in the triple alpha resonance energy and study its effects on the resulting abundances of $^{12}$C, $^{16}$O, and larger nuclei. In our universe, the $0^{+}$ energy level of carbon supports a resonant nuclear reaction that dominates carbon synthesis in stellar cores and accounts for the observed cosmic abundances. Here we define $\Delta{E}_R$ to be the change in this resonant energy level, and show how different values affect the cosmic abundances of the intermediate alpha elements. Using the state of the art computational package $MESA$, we carry out stellar evolution calculations for massive stars in the range $M_\ast$ = $15-40M_\odot$, and for a wide range of resonance energies. We also include both solar and low metallicity initial conditions. For negative $\Delta{E}_R$ , carbon yields are increased relative to standard stellar models, and such universes remain viable as long as the production of carbon nuclei remains energetically favorable, and stars remain stable, down to $\Delta{E}_R\approx-300$ keV. For positive $\Delta{E}_R$, carbon yields decrease, but significant abundances can be produced for resonance energy increments up to $\Delta{E}_R\approx+500$ keV. Oxygen yields tend to be anti-correlated with those of carbon, and the allowed range in $\Delta{E}_R$ is somewhat smaller. We also present yields for neon, magnesium, and silicon. With updated stellar evolution models and a more comprehensive survey of parameter space, these results indicate that the range of viable universes is larger than suggested by earlier studies.

The structure of force-free, steady and axisymmetric magnetosphere of a neutron star (NS) is governed by the Grad-Shafranov (GS) equation, which is a second-order differential equation but degrades to first-order on the light surface (LS). The key to numerically solving the GS equation is to enable magnetic field lines smoothly cross the LS, and crossing a straight LS in flat spacetime has been a well-studied problem. But the numerical algorithm implementation becomes complicate in the presence of a bent LS, e.g. in curved spacetime, since there is no suitable computation grid adapted to it. We propose to circumvent this grid-LS mismatch problem by introducing a new coordinate frame designed such that the LS in it is a straight line. As an application, we investigate the general relativistic (GR) effect in magnetosphere structure of rapidly rotating pulsars in detail, where the LS is bent towards the central NS. We split the GR effect into two parts, curvature and frame-dragging; measure each of them and examine their dependence on the NS mass and the angular velocity for pulsars embedded in aligned dipole and multipole magnetic fields. Qualitatively speaking, we find that the curvature effect compactifies the magnetic field lines near the NS, therefore reduces the open magnetic flux and the Poynting luminosity, while the frame-dragging effect contributes a minor part in shaping the magnetosphere structure but plays a role in enhancing the spacelike current generation.

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.

The gamma-ray burst GRB 170817A associated with GW170817 is subluminous and subenergetic compared with other typical short GRBs. It may be due to a relativistic jet viewed off-axis, or a structured jet, or cocoon emission. Giant flares from magnetars may possibly be ruled out. However, the luminosity and energetics of GRB 170817A is coincident with that of magnetar giant flares. After the coalescence of the binary neutron star, a hypermassive neutron star may be formed. The hypermassive neutron star may have magnetar-strength magnetic field. During the collapse of the hypermassive neutron star, the magnetic field energy will also be released. This giant-flare-like event may explain the the luminosity and energetics of GRB 170817A. Bursts with similar luminosity and energetics are expected in future neutron star-neutron star or neutron star-black hole mergers.

We have analyzed the multipolar magnetic field structure variation at neutron star surface by means of the catastrophic eruption model, and find that the variation of the geometry of multipolar fields on the magnetar surface could result in the catastrophic rearrangement of the magnetosphere, which provides certain physical mechanism for the outburst of giant flares. The magnetospheric model we adopted consists of two assumptions: a helically twisted flux rope is suspended in an ideal force-free magnetosphere around the magnetar, and a current sheet emerges during the flux rope evolution. Magnetic energy accumulates during the flux rope's gradual evolution along with the variation of magnetar surface magnetic structure before the eruption. The two typical behaviors, either state transition or catastrophic escape, would take place once the flux rope loses equilibrium, thus tremendous accumulated energy is radiated. We have investigated the equilibrium state of the flux rope and the energy release affected by different multipolar structures, and find structures that could trigger violent eruption and provide the radiation approximately 0.5$\%$ of the total magnetic energy during the giant flare outburst. Our results provide certain multipolar structures of the neutron star's magnetic field with an energy release percentage 0.42$\%$ in the state transition and 0.51$\%$ in the catastrophic escape case, which are sufficient for the previously reported energy release from SGR 1806-20 giant flares.

The Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) began millimeter-precision ranging to the Moon in 2006. Until now, a comprehensive validation of APOLLO system range accuracy has not been possible because of centimeter-scale deficiencies in computational models of the Earth-Moon range, and because APOLLO lacked an internal timing calibration system. Here, we report on the development of a system that enables in-situ calibration of the timing response of the APOLLO apparatus, simultaneous with lunar range measurements. The system was installed in August 2016. Preliminary results show that the APOLLO system can provide lunar range measurements with millimeter accuracy.

Lunar laser ranging provides a number of leading experimental tests of gravitation -- important in our quest to unify General Relativity and the Standard Model of physics. The Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) has for years achieved median range precision at the ~2 mm level. Yet residuals in model-measurement comparisons are an order-of-magnitude larger, raising the question of whether the ranging data are not nearly as accurate as they are precise, or if the models are incomplete or ill-conditioned. This paper describes a new absolute calibration system (ACS) intended both as a tool for exposing and eliminating sources of systematic error, and also as a means to directly calibrate ranging data in-situ. The system consists of a high-repetition-rate (80 MHz) laser emitting short (< 10 ps) pulses that are locked to a cesium clock. In essence, the ACS delivers photons to the APOLLO detector at exquisitely well-defined time intervals as a "truth" input against which APOLLO's timing performance may be judged and corrected. Preliminary analysis indicates no inaccuracies in APOLLO data beyond the ~3 mm level, suggesting that historical APOLLO data are of high quality and motivating continued work on model capabilities. The ACS provides the means to deliver APOLLO data both accurate and precise below the 2 mm level.

The structure of steady axisymmetric force-free magnetosphere of a Kerr black hole (BH) is governed by a second-order partial differential equation of $A_\phi$ depending on two "free" functions $\Omega(A_\phi)$ and $I(A_\phi)$, where $A_\phi$ is the $\phi$ component of the vector potential of the electromagnetic field, $\Omega$ is the angular velocity of the magnetic field lines and $I$ is the poloidal electric current. In this paper, we investigate the solution uniqueness. Taking asymptotically uniform field as an example, analytic studies imply that there are infinitely many solutions approaching uniform field at infinity, while only a unique one is found in general relativistic magnetohydrodynamic simulations. To settle down the disagreement, we reinvestigate the structure of the governing equation and numerically solve it with given constraint condition and boundary condition. We find that the constraint condition (field lines smoothly crossing the light surface (LS)) and boundary conditions at horizon and at infinity are connected via radiation conditions at horizon and at infinity, rather than being independent. With appropriate constraint condition and boundary condition, we numerically solve the governing equation and find a unique solution. Contrary to naive expectation, our numerical solution yields a discontinuity in the angular velocity of the field lines and a current sheet along the last field line crossing the event horizon. We also briefly discuss the applicability of the perturbation approach to solving the governing equation.

Magnetic B-type stars exhibit photometric variability due to diverse causes, and consequently on a variety of timescales. In this paper we describe interpretation of BRITE photometry and related ground-based observations of 4 magnetic B-type systems: $\epsilon$ Lupi, $\tau$ Sco, a Cen and $\epsilon$ CMa.

A sample of the LAMOST spectra of the early type M0-M3 dwarfs is compared with the Kepler observations. It is found that M dwarfs with strong chromospheric emission in $H_{\alpha}$ have large flare activity in general. The rotational periods derived from the Kepler measurements have close correlations with the sizes of the flares, the power-law distribution index and the equivalent widths of the $H_{\alpha}$ emission. A clear trend exists for higher magnetic activities being detected in faster rotating M dwarfs (rotation periods $<$ 20 day).

We propose that the magnetosphere reconfiguration induced by magnetic twists in the closed field line region can account for the mode-switching of intermittent pulsars. We carefully investigate the properties of axisymmetric force-free pulsar magnetospheres with magnetic twists in closed field line region around the polar caps. The magnetosphere with twisted closed lines leads to enhanced spin-down rates. The enhancement in spin-down rate depends on the size of region with twisted closed lines. Typically, it is increased by a factor of $\sim2$, which is consistent with the intermittent pulsars' spin down behavior during the `off' and `on' states. We find there is a threshold of maximal twist angle $\Delta\phi_{\rm thres}\sim1$. The magnetosphere is stable only if the closed line twist angle is less than $\Delta\phi_{\rm thres}$. Beyond this value, the magnetosphere becomes unstable and gets untwisted. The spin-down rate would reduce to its off-state value. The quasi-periodicity in spin-down rate change can be explained by long-term activities in star's crust and the untwisting induced by MHD instability. The estimated duration time of on-state is about one week, consistent with observations. Due to the MHD instability, there exists an upper limit for the spin down ratio ($f\sim3$) between the on-state and the off-state, if the Y-point remains at the light cylinder.

We report the discovery of 2 new Be stars, and re-identify one known Be star in the open cluster NGC 6830. Eleven H-alpha emitters were discovered using the H-alpha imaging photometry of the Palomar Transient Factory Survey. Stellar membership of the candidates was verified with photometric and kinematic information using 2MASS data and proper motions. The spectroscopic confirmation was carried out by using the Shane 3-m telescope at Lick observatory. Based on their spectral types, three H-alpha emitters were confirmed as Be stars with H-alpha equivalent widths > -10 Angstrom. Two objects were also observed by the new spectrograph SED-Machine on the Palomar 60 inch Telescope. The SED-Machine results show strong H-alpha emission lines, which are consistent with the results of the Lick observations. The high efficiency of the SED-Machine can provide rapid observations for Be stars in a comprehensive survey in the future.

PSR B1821$-$24 is a solitary millisecond pulsar (MSP) which radiates multi-wavelength pulsed photons. It has complex radio, X-ray and $\gamma$-ray pulse profiles with distinct peak phase-separations that challenge the traditional caustic emission models. Using the single-pole annular gap model with suitable magnetic inclination angle ($\alpha=40^\circ$) and viewing angle ($\zeta=75^\circ$), we managed to reproduce its pulse profiles of three wavebands. It is found that the middle radio peak is originated from the core gap region at high altitudes, and the other two radio peaks are originated from the annular gap region at relatively low altitudes. Two peaks of both X-ray and $\gamma$-ray wavebands are fundamentally originated from annular gap region, while the $\gamma$-ray emission generated from the core gap region contributes somewhat to the first $\gamma$-ray peak. Precisely reproducing the multi-wavelength pulse profiles of PSR B1821$-$24 enables us to understand emission regions of distinct wavebands and justify pulsar emission models.

We propose a physical mechanism to explain giant flares and radio afterglows in terms of a magnetospheric model containing both a helically twisted flux rope and a current sheet (CS). With the appearance of CS, we solve a mixed boundary value problem to get the magnetospheric field based on a domain decomposition method. We investigate properties of the equilibrium curve of the flux rope when the CS is present in background multipolar fields. In response to the variations at the magnetar surface, it quasi-statically evolves in stable equilibrium states. The loss of equilibrium occurs at a critical point and, beyond that point, it erupts catastrophically. New features show up when the CS is considered. Especially, we find two kinds of physical behaviors, i.e., catastrophic state transition and catastrophic escape. Magnetic energy would be released during state transitions. The released magnetic energy is sufficient to drive giant flares. The flux rope would go away from the magnetar quasi-statically, which is inconsistent with the radio afterglow. Fortunately, in the latter case, i.e., the catastrophic escape, the flux rope could escape the magnetar and go to infinity in a dynamical way. This is more consistent with radio afterglow observations of giant flares. We find that the minor radius of flux rope has important implications for its eruption. Flux ropes with larger minor radius are more prone to erupt. We stress that the CS provides an ideal place for magnetic reconnection, which would further enhance the energy release during eruptions.

We present an analysis of the metallicity and star formation activities of HII regions in the interacting system Arp 86, based on the first scientific observation of the multi-object spectroscopy on the 2.16m Telescope at Xinglong Observatory. We find that the oxygen abundance gradient in Arp 86 is flatter than that in normal disk galaxies, which confirms that gas inflows caused by tidal forces during encounters can flatten the metallicity distributions in galaxies. The companion galaxy NGC 7752 is currently experiencing a galaxy-wide starburst with higher surface density of star formation rate than the main galaxy NGC 7753, which can be explained that the companion galaxy is more susceptible to the effects of interaction than the primary. We also find that the galaxy 2MASX J23470758+2926531 has similar abundance and star formation properties to NGC 7753, and may be a part of the Arp 86 system.

We propose a catastrophic eruption model for magnetar's enormous energy release during giant flares, in which a toroidal and helically twisted flux rope is embedded within a force-free magnetosphere. The flux rope stays in stable equilibrium states initially and evolves quasi-statically. Upon the loss of equilibrium point is reached, the flux rope cannot sustain the stable equilibrium states and erupts catastrophically. During the process, the magnetic energy stored in the magnetosphere is rapidly released as the result of destabilization of global magnetic topology. The magnetospheric energy that could be accumulated is of vital importance for the outbursts of magnetars. We carefully establish the fully open fields and partially open fields for various boundary conditions at the magnetar surface and study the relevant energy thresholds. By investigating the magnetic energy accumulated at the critical catastrophic point, we find that it is possible to drive fully open eruptions for dipole dominated background fields. Nevertheless, it is hard to generate fully open magnetic eruptions for multipolar background fields. Given the observational importance of the multipolar magnetic fields in the vicinity of the magnetar surface, it would be worthwhile to explore the possibility of the alternative eruption approach in multipolar background fields. Fortunately, we find that flux ropes may give rise to partially open eruptions in the multipolar fields, which involve only partial opening up of background fields. The energy release fractions are greater for cases with central-arcaded multipoles than those with central-caved multipoles emerged in background fields. Eruptions would fail only when the centrally-caved multipoles become extremely strong.

We propose a catastrophic magnetospheric model for magnetar precursors and their successive giant flares. Axisymmetric models of the magnetosphere, which contain both a helically twisted flux rope and a current sheet, are established based on force-free field configurations. In this model, the helically twisted flux rope would lose its equilibrium and erupt abruptly in response to the slow and quasi-static variations at the ultra-strongly magnetized neutron star's surface. In a previous model without current sheets, only one critical point exists in the flux rope equilibrium curve. New features show up in the equilibrium curves for the flux rope when current sheets appear in the magnetosphere. The causal connection between the precursor and the giant flare, as well as the temporary re-entry of the quiescent state between the precursor and the giant flare, can be naturally explained. Magnetic energy would be released during the catastrophic state transitions. The detailed energetics of the model are also discussed. The current sheet created by the catastrophic loss of equilibrium of the flux rope provides an ideal place for magnetic reconnection. We point out the importance of magnetic reconnection for further enhancement of the energy release during eruptions.

We present the gamma-ray, X-ray, optical and radio data for GRB100814A. At the end of the slow decline phase of the X-ray and optical afterglow, a sudden and prominent rebrightening in the optical band occurs followed by a fast decay in both bands. This optical rebrightening is accompanied by possible chromatic variations. We discuss possible interpretations, such as double component scenarios and internal dissipation mechanism, with their virtues and drawbacks. We also compare GRB100814A with other Swift bursts that show optical rebrightenings with similar properties.

We use a model of an accretion flow coupled with an emergent flare to interpret the latest 1.3mm VLBI measurements for Sagittarius A*. The visibility data constrained the distances from the flare center to the black hole center as $d_{\rm EW}\lesssim20{\rm R_g}$ and $d_{\rm NS}\lesssim80{\rm R_g}$ in the East-West and North-South directions, respectively. If interpreted by the hot-spot model, the flare was preferred to pass in front of the black hole at a radius much larger than $d_{\rm EW}$. If interpreted by the episodic jet launched from a nearly edge-on hot accretion flow, the flare was preferred to be ejected with $\theta_{\rm j}\gtrsim40^\circ$ off the black hole rotating axis. This method can be generalized to help us understand future sub-millimeter VLBI observations, and study the millimeter/sub-millimeter variabilities in the vicinity of the Galactic Center supermassive black hole.

At high redshift, the universe was so young that core-collapse supernovae (SNe) are suspected to be the dominant source of dust production. However, some observations indicate that the dust production by SNe is an inefficient process, casting doubts on the existence of abundant SNe-dust in the early universe. Recently, Perley et al. (2010) reported that the afterglow of GRB 071025 - an unusually red GRB at z ~ 5 - shows evidence for the SNe-produced dust. Since this is perhaps the only high redshift GRB exhibiting compelling evidence for SNe-dust but the result could easily be affected by small systematics in photometry, we re-examined the extinction properties of GRB 071025 using our own optical/near-infrared data at a different epoch. In addition, we tested SNe-dust models with different progenitor masses and dust destruction efficiencies to constrain the dust formation mechanisms. By searching for the best-fit model of the afterglow spectral energy distribution, we confirm the previous claim that the dust in GRB 071025 is most likely to originate from SNe. We also find that the SNe-dust model of 13 or 25M_{\odot} without dust destruction fits the extinction property of GRB 071025 best, while pair-instability SNe (PISNe) models with a 170M_{\odot} progenitor poorly fit the data. Our results indicate that, at least in some systems at high redshift, SNe with intermediate masses within 10 - 30M_{\odot} were the main contributors for the dust enrichment, and the dust destruction effect due to reverse shock was negligible.

We provide precise fitting formulae for Faraday conversion and rotation coefficients in uniformly magnetized relativistic plasma. The formulae are immediately applicable to Rotation Measure and Circular Polarization (CP) production in jets and hot accretion flows. We show the recipe and results for arbitrary isotropic particle distributions, in particular thermal and power-law. The exact Faraday conversion coefficient is found to approach zero with the increasing particle energy. The non-linear corrections of Faraday conversion and rotation coefficients are found essential for reliable CP inter- pretation of Sgr A*.

The rising amount of polarized observations of relativistic sources requires the correct theory for proper model fitting. The equations for general relativistic (GR) polarized radiative transfer are derived starting from the Boltzmann equation and basic ideas of general relativity. The derivation is aimed at providing a practical guide to reproducing the synchrotron part of radio & sub-mm emission from low luminosity active galactic nuclei (LLAGNs), in particular Sgr A*, and jets. The recipe for fast exact calculation of cyclo-synchrotron emissivities, absorptivities, Faraday rotation and conversion coefficients is given for isotropic particle distributions. The multitude of physical effects influencing simulated spectrum is discussed. The application of the prescribed technique is necessary to determine the black hole (BH) spin in LLAGNs, constraining it with all observations of total flux, linear and circular polarization fractions, and electric vector position angle as functions of the observed frequency.

The vicinity of the supermassive black hole associated with the compact radio source Sagittarius (Sgr) A* is believed to dominate the observed emission at wavelengths near and shorter than $\sim$ 1 millimeter. We show that a general relativistic accretion flow, heated via the plasma wave heating mechanism, is consistent with the polarization and recent mm-VLBI observations of Sgr A* for an inclination angle of $\sim 45^\circ$, position angle of $\sim 140^\circ$, and spin $\lesssim 0.9$. Structure in visibilities produced by the black hole shadow can potentially be observed by 1.3 mm-VLBI on the existing Hawaii-CARMA and Hawaii-SMT baselines. We also consider eight additional potential mm-VLBI stations, including sites in Chile and New Zealand, finding that with these the basic geometry of the emission region can be reliably estimated.

We explore the parameter space of the two temperature pseudo-Newtonian Keplerian accretion flow model for the millimeter and shorter wavelength emission from Sagittarius A*. A general relativistic ray-tracing code is used to treat the radiative transfer of polarized synchrotron emission from the flow. The synchrotron self-Comptonization and bremsstrahlung emission components are also included. It is shown that the model can readily account for the millimeter to sub-millimeter emission characteristics with an accretion rate of ~6x10^17g.s^-1 and an inclination angle of ~40 deg. However, the corresponding model predicted near-infrared and X-ray fluxes are more than one order of magnitude lower than the observed 'quiescent' state values. While the extended quiescent-state X-ray emission has been attributed to thermal emission from the large-scale accretion flow, the NIR emission and flares are likely dominated by emission regions either within the last stable orbit of a Schwarzschild black hole or associated with outflows. With the viscous parameter derived from numerical simulations, there is still a degeneracy between the electron heating rate and the magnetic parameter. A fully general relativistic treatment with the black hole spin incorporated will resolve these issues.

We have performed monitoring observations of the 3-mm flux density toward the Galactic Center compact radio source Sgr A* with the Australia Telescope Compact Array since 2005 October. Careful calibrations of both elevation-dependent and time-dependent gains have enabled us to establish the variability behavior of Sgr A*. Sgr A* appeared to undergo a high and stable state in 2006 June session, and a low and variable state in 2006 August session. We report the results, with emphasis on two detected intra-day variation events during its low states. One is on 2006 August 12 when Sgr A* exhibited a 33% fractional variation in about 2.5 hr. The other is on 2006 August 13 when two peaks separated by about 4 hr, with a maximum variation of 21% within 2 hr, were seen. The observed short timescale variations are discussed in light of two possible scenarios, i.e., the expanding plasmon model and the sub-Keplerian orbiting hot spot model. The fitting results indicate that for the adiabatically expanding plasmon model, the synchrotron cooling can not be ignored, and a minimum mass-loss rate of 9.7*10^{-10}M_sun /yr is obtained based on parameters derived for this modified expanding plasmon model. Simultaneous multi-wavelength observation is crucial to our understanding the physical origin of rapid radio variability in Sgr A*.

In fully general relativity, we calculate the images of the radiatively inefficient accretion flow (RIAF) surrounding a Kerr black hole with arbitrary spins, inclination angles, and observational wavelengths. For the same initial conditions, such as the fixed accretion rate, it is found that the intrinsic size and radiation intensity of the images become larger, but the images become more compact in the inner region, while the size of the black hole shadow decreases with the increase of the black hole spin. With the increase of the inclination angles, the shapes of the black hole shadows change and become smaller, even disappear at all due to the obscuration by the thick disks. For median inclination angles, the radial velocity observed at infinity is larger because of both the rotation and radial motion of the fluid in the disk, which results in the luminous part of the images is much brighter. For larger inclination angles, such as the disk is edge on, the emission becomes dimmer at longer observational wavelengths (such as at 7.0mm and 3.5mm wavelengths), or brighter at shorter observational wavelengths (such as at 1.3 mm wavelength) than that of the face on case, except for the high spin and high inclination images. These complex behaviors are due to the combination of the Lorentz boosting effect and the radiative absorption in the disk. We hope our results are helpful to determine the spin parameter of the black hole in low luminosity sources, such as the Galactic center. A primary analysis by comparison with the observed sizes of Sgr A* at millimeters strongly suggests that the disk around the central black hole at Sgr A* is highly inclined or the central black hole is rotating fast.

We perform general relativistic ray-tracing calculations of the transfer of polarized synchrotron radiation through the relativistic accretion flow in Sagittarius (Sgr) A*. Based on a two-temperature magneto-rotational-instability (MRI) induced accretion mode, the birefringence effects are treated self-consistently. By fitting the spectrum and polarization of Sgr A* from millimeter to near-infrared bands, we are able to not only constrain the basic parameters related to the MRI and the electron heating rate, but also limit the orientation of the accretion torus. These constraints lead to unique polarimetric images, which may be compared with future millimeter and sub-millimeter VLBI observations. In combination with general relativistic MHD simulations, the model has the potential to test the MRI with observations of Sgr A*.

The compact dark objects with very large masses residing at the centres of galaxies are believed to be black holes. Due to the gravitational lensing effect, they would cast a shadow larger than their horizon size over the background, whose shape and size can be calculated. For the supermassive black hole candidate Sgr A*, this shadow spans an angular size of about 50 micro arc second, which is under the resolution attainable with the current astronomical instruments. Such a shadow image of Sgr A* will be observable at about 1 mm wavelength, considering the scatter broadening by the interstellar medium. By simulating the black hole shadow image of Sgr A* with the radiatively inefficient accretion flow model, we demonstrate that analyzing the properties of the visibility function can help us determine some parameters of the black hole configuration, which is instructive to the sub-millimeter VLBI observations of Sgr A* in the near future.

We report on VLBA observations of a gamma-ray bright blazar NRAO 530 at multiple frequencies (5, 8, 15, 22, 39, 43, and 45 GHz) in 1997 and 1999. These multi-epoch multi-frequency high-resolution VLBI images exhibit a consistent core-dominated morphology with a bending jet to the north of the core. The quasi-simultaneous data observed at five frequencies (5, 8, 15, 22 and 43 GHz) in February 1997 enable us to estimate the spectra of compact VLBI components in this highly variable source. Flat spectra are seen in central two components (A and B), and the most compact component A with the flattest spectral index at the south end is identified as the core.

Recent radio observations by the VLBA at 7 and 3.5 mm produced the high-resolution images of the compact radio source located at the center of our Galaxy--Sgr A*, and detected its wavelength-dependent intrinsic sizes at the two wavelengths. This provides us with a good chance of testing previously-proposed theoretical models for Sgr A*. In this {\em Letter}, we calculate the size based on the radiatively inefficient accretion flow (RIAF) model proposed by Yuan, Quataert & Narayan (2003). We find that the predicted sizes after taking into account the scattering of the interstellar electrons are consistent with the observations. We further predict an image of Sgr A* at 1.3 mm which can be tested by future observations.

We present a correlation analysis between the flux ratio of high-energy gamma-ray emission to synchrotron emission and the Doppler factor for a sample of EGRET AGNs. The results favor a model that attributes the EGRET emission of AGNs to inverse-Compton scattering on photons external to the jet.

After comparing the flux densities of a sample of active galactic nuclei detected by energetic gamma-ray experiment telescope at 90 and 230 GHz with the $\gamma$-ray emissions detected by Compton Gamma Ray Observatory and x-ray emission, a strong correlation between the emissions at the millimeter wavelength and the $\gamma$-ray emission is found. The average flux density of x-ray is almost proportional to the average flux density at the millimeter wavelength for quasars detected by energetic gamma-ray experiment telescope, which strongly supports the previous idea that the x-ray emissions of this kind sources are mainly produced by Synchrotron Self-Compton process.