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

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

Using numerical simulations we show that low-amplitude Alfv\'en waves from a magnetar quake propagate to the outer magnetosphere and convert to "plasmoids" (closed magnetic loops) which accelerate from the star, driving blast waves into the magnetar wind. Quickly after its formation, the plasmoid becomes a thin relativistic pancake. It pushes out the magnetospheric field lines, and they gradually reconnect behind the pancake, generating a variable wind far stronger than the normal spindown wind of the magnetar. Repeating ejections drive blast waves in the amplified wind. We suggest that these ejections generate the simultaneous X-ray and radio bursts detected from SGR 1935+2154. A modest energy budget of the magnetospheric perturbation $\sim 10^{40}$ erg is sufficient to produce the observed bursts with the energy ratio $\mathcal{E}_X/\mathcal{E}_{\rm radio}\sim 10^5$. Our simulation predicts a narrow (a few ms) X-ray spike from the magnetosphere, arriving simultaneously with the radio burst emitted far outside the magnetosphere. This timing is caused by the extreme relativistic motion of the ejecta.

We study the equivalence principle, regarded as the cornerstone of general relativity, by analyzing the deformation observable of black hole shadows. Such deformation can arise from new physics and may be expressed as a phenomenological violation of the equivalence principle. Specifically, we assume that there is an additional background vector field that couples to the photons around the supermassive black hole. This type of coupling yields impact on the way the system depends on initial conditions, and affects the black hole shadow at different wavelengths by a different amount, and therefore observations of the shadow in different wavelengths could constrain such couplings. This can be tested by future multi-band observations. Adopting a specific form of the vector field, we obtain constraints on model parameters from Event Horizon Telescope observations and measurements of gas/stellar orbits.

The exact relationship between the long gamma-ray burst (LGRB) rate and the cosmic star formation rate (CSFR) is essential for using LGRBs as cosmological probes. In this work, we collect a large sample composed of 371 Swift LGRBs with known redshifts and prompt emission properties. We first compare the rest-frame prompt properties of these bursts in different redshift bins, finding negligible redshift evolution of the luminosity of LGRBs with $L_{\mathrm{iso}}\gtrsim10^{51}\,\mathrm{erg\, s^{-1}}$ between $z\sim1$ and $z\sim4$. Then, by utilizing the CSFR obtained from the large-scale cosmological hydrodynamical simulation, the Illustris simulation, we calculate the cumulative redshift distribution of LGRBs under different metallicity thresholds. After comparing with our sample, we find that the predictions with a moderate threshold between $0.3\,Z_{\odot}\leqslant Z_{\mathrm{th}}\leqslant1.0\,Z_{\odot}$ are consistent with the sample between redshift $0<z<3$, while at higher redshifts between $3<z<5$, all metallicity thresholds fit the data well. When changing to an empirical model based on observations, the predictions show similar results as well. After comparing with the metallicity distribution of the observed LGRB host galaxies between $0<z<1$, we confirm that the production of LGRBs in galaxies with super-solar metallicity is suppressed. Nevertheless, considering that a significant fraction of stars are born in sub-solar metallicity environments at $z\gtrsim3$, we suggest that, as a first approximation, LGRBs can be used as direct tracers of the CSFR in this redshift range.

The exact relationship between the long gamma-ray burst (LGRB) rate and the cosmic star formation rate (CSFR) is essential for using LGRBs as cosmological probes. In this work, we collect a large sample composed of 371 Swift LGRBs with known redshifts and prompt emission properties. We first compare the rest-frame prompt properties of these bursts in different redshift bins, finding negligible redshift evolution of the luminosity of LGRBs with $L_{\mathrm{iso}}\gtrsim10^{51}\,\mathrm{erg\, s^{-1}}$ between $z\sim1$ and $z\sim4$. Then, by utilizing the CSFR obtained from the large-scale cosmological hydrodynamical simulation, the Illustris simulation, we calculate the cumulative redshift distribution of LGRBs under different metallicity thresholds. After comparing with our sample, we find that the predictions with a moderate threshold between $0.3\,Z_{\odot}\leqslant Z_{\mathrm{th}}\leqslant1.0\,Z_{\odot}$ are consistent with the sample between redshift $0<z<3$, while at higher redshifts between $3<z<5$, all metallicity thresholds fit the data well. When changing to an empirical model based on observations, the predictions show similar results as well. After comparing with the metallicity distribution of the observed LGRB host galaxies between $0<z<1$, we confirm that the production of LGRBs in galaxies with super-solar metallicity is suppressed. Nevertheless, considering that a significant fraction of stars are born in sub-solar metallicity environments at $z\gtrsim3$, we suggest that, as a first approximation, LGRBs can be used as direct tracers of the CSFR in this redshift range.

We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $\tau(\lambda)$, generally follow the relation $\tau(\lambda) \propto \lambda^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $\tau(\lambda) \propto \lambda^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_\nu\propto\nu^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate.

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.

Theoretical and observational arguments suggest that there is a large amount of hot ($\sim 10^6$ K), diffuse gas residing in the Milky Way's halo, while its total mass and spatial distribution are still unclear. In this work, we present a general model for the gas density distribution in the Galactic halo, and investigate the gas evolution under radiative cooling with a series of 2D hydrodynamic simulations. We find that the mass inflow rate in the developed cooling flow increases with gas metallicity and the total gas mass in the halo. For a fixed halo gas mass, the spatial gas distribution affects the onset time of the cooling catastrophe, which starts earlier when the gas distribution is more centrally-peaked, but does not substantially affect the final mass inflow rate. The gravity from the Galactic bulge and disk affects gas properties in inner regions, but has little effect on the final inflow rate either. We confirm our results by investigating cooling flows in several density models adopted from the literature, including the Navarro-Frenk-White (NFW) model, the cored-NFW model, the Maller & Bullock model, and the $\beta$ model. Typical mass inflow rates in our simulations range from $\sim 5 M_{\odot}$ yr$^{-1}$ to $\sim 60 M_{\odot}$ yr$^{-1}$, and are much higher than the observed star formation rate in our Galaxy, suggesting that stellar and active galactic nucleus feedback processes may play important roles in the evolution of the Milky Way (MW) and MW-type galaxies.

Recent modeling of Neutron Star Interior Composition Explorer(NICER) observations of the millisecond pulsar PSR J0030+0451 suggests that the magnetic field of the pulsar is non-dipolar. We construct a magnetic field configuration where foot points of the open field lines closely resemble the hotspot configuration from NICER observations. Using this magnetic field as input, we perform force-free simulations of the magnetosphere of PSR J0030+0451, showing the three-dimensional structure of its plasma-filled magnetosphere. Making simple and physically motivated assumptions about the emitting regions, we are able to construct the multi-wavelength lightcurves that qualitatively agree with the corresponding observations. The agreement suggests that multipole magnetic structures are the key to modeling this type of pulsars, and can be used to constrain the magnetic inclination angle and the location of radio emission.

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$\beta$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$\beta$, He II, He I, and Fe II lines with respect to the varying 5100\AA continuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

We present the two-dimensional inflow-outflow solutions of radiation hydrodynamic (RHD) equations of supercritical accretion flows. Compared with prior studies, we include all components of the viscous stress tensor. We assume steady state flow and use self-similar solutions in the radial direction to solve the equations in $ r-\theta $ domain of the spherical coordinates. The set of differential equations have been integrated from the rotation axis to the equatorial plane. We find that the self-similarity assumption requires that the radial profile of density is described by $ \rho(r) \propto r^{-0.5} $. Correspondingly, the radial profile of the mass inflow rate decreases with decreasing radii as $ \dot{M}_\mathrm{in} \propto r $. Inflow-outflow structure has been found in our solution. In the region $ \theta > 65^{\circ} $ there exist inflow while above that flow moves outward and outflow could launch. The driving forces of the outflow are analyzed and found that the radiation force is dominant and push the gas particles outwards with poloidal velocity $ \sim 0.25 c $. The properties of outflow are also studied. The results show that the mass flux weighted angular momentum of the inflow is lower than that of outflow, thus the angular momentum of the flow can be transported by the outflow. We also analyze the convective stability of the supercritical disk and find that in the absence of the magnetic field, the flow is convectively unstable. Our analytical results are fully consistent with the previous numerical simulations of the supercritical accretion flow.

We study the long-term variability in the optical monitoring database of Ark~120, a nearby radio-quiet active galactic nucleus (AGN) at a distance of 143 Mpc (z=0.03271). We compiled the historical archival photometric and spectroscopic data since 1974 and conducted a new two-year monitoring campaign in 2015-2017, resulting in a total temporal baseline over four decades. The long-term variations in the optical continuum exhibit a wave-like pattern and the Hbeta integrated flux series varies with a similar behavior. The broad Hbeta profiles have asymmetric double peaks, which change strongly with time and tend to merge into a single peak during some epochs. The period in the optical continuum determined from various period-search methods is about 20 yr and the estimated false alarm probability with null hypothesis simulations is about 1*10^-3. The overall variations of the broad Hbeta profiles also follow the same period. However, the present database only covers two cycles of the suggested period, which strongly encourages continued monitoring to track more cycles and confirm the periodicity. Nevertheless, in light of the possible periodicity and the complicated Hbeta profile, Ark~120 is one candidate of the nearest radio-quiet AGNs with possible periodic variability, and it is thereby a potential candidate host for a sub-parsec supermassive black hole binary.

This is the second paper in a series where we examine the physics of pair producing gaps in low-luminosity accreting supermassive black hole systems. In this paper, we carry out time-dependent self-consistent fully general relativistic 1D PIC simulations of the gap, including full inverse Compton scattering and photon tracking. Similar to the previous paper, we find a highly time-dependent solution where a macroscopic vacuum gap can open quasi-periodically, producing bursts of $e^\pm$ pairs and high energy radiation. We present the light curve, particle and photon spectra from this process. Using an empirical scaling relation, we rescale the parameters to the inferred values at the base of the jet in M87, and find that the observed TeV flares could potentially be explained by this model under certain parameter assumptions.

Inverse Compton-pair cascades are initiated when gamma-rays are absorbed on an ambient soft photon field to produce relativistic pairs, which in turn up-scatter the same soft photons to produce more gamma-rays. If the Compton scatterings take place in the deep Klein-Nishina regime, then triplet pair production ($e\gamma_b \rightarrow ee^{+}e^{-}$) becomes relevant and may even regulate the development of the cascade. We investigate the properties of pair-Compton cascades with triplet pair production in accelerating gaps, i.e., regions with an unscreened electric field. Using the method of transport equations for the particle evolution, we compute the growth rate of the pair cascade as a function of the accelerating electric field in the presence of black-body and power-law ambient photon fields. Informed by the numerical results, we derive simple analytical expressions for the peak growth rate and the corresponding electric field. We show that for certain parameters, which can be realized in the vicinity of accreting supermassive black holes at the centers of active galactic nuclei, the pair cascade may well be regulated by inverse Compton scattering in the deep Klein-Nishina regime and triplet pair production. We present indicative examples of the escaping gamma-ray radiation from the gap, and discuss our results in application to the TeV observations of radio galaxy M87.

Proton-proton interaction and photo-hadronic interaction in cosmic accelerators are the two main channels for the production of cosmic ultra-high energy photons and neutrinos (TeV-PeV). In this Letter, we use FWW approach to obtain the production of cosmic ultra-high energy photons and neutrinos from the collision between UHE proton with magnetic field which could be considered as the virtualphoton in the rest frame of UHE proton. We name this as $pB$ process. The threshold for the occurrence of the $pB$ process is that the combination of the Lorentz factor of proton and the strength of the magnetic field is about $\gamma_p B \simeq 5\times 10^{18}$Gauss. Beyond this threshold, the rate of energy loss of proton due to the $pB$ process is about three orders higher than that due to the synchrotron radiation of proton in the same magnetic field. The $pB$ process might potentially happen in the atmosphere of white dwarfs, neutron stars or even that of stellar massive black holes.

In some Seyfert Galaxies, the hard X-rays that produce fluorescent emission lines are thought to be generated in a hot corona that is compact and located at only a few gravitational radii above the supermassive black hole. We consider the possibility that this X-ray source may be powered by small scale magnetic flux tubes attached to the accretion disk near the black hole. We use three dimensional, time dependent force-free simulations in a simplified setting to study the dynamics of such flux tubes as they get continuously twisted by the central compact star/black hole. We find that, the dynamical evolution of the flux tubes connecting the central compact object and the accretion disk is strongly influenced by the confinement of the surrounding field. Although differential rotation between the central object and the disk tends to inflate the flux tubes, strong confinement from surrounding field quenches the formation of a jet-like outflow, as the inflated flux tube becomes kink unstable and dissipates most of the extracted rotational energy relatively close to the central object. Such a process may be able to heat up the plasma and produce strong X-ray emission. We estimate the energy dissipation rate and discuss its astrophysical implications.

We compiled a sample of 73 active galactic nuclei (AGNs) with reverberation mapping (RM) observations from RM campaigns including our ongoing campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs). This sample covers a large range of black hole (BH) mass $(M_{\bullet}=10^{6-9}~M_\odot)$, dimensionless accretion rates $(\dot{\mathscr{M}}=10^{-2.7}-10^{2.7})$ and 5100~\AA~luminosity $(L_{5100}=10^{42-46}~\rm erg~s^{-1})$, allowing us to systematically study the AGN variability and their relations with BH mass, accretion rates, and optical luminosity. We employed the damped random walk (DRW) model to delineate the optical variability of continuum at 5100~\AA~and obtained damped variability timescale ($\tau_{\rm d}$) and amplitude ($\sigma_{\rm d}$) using a Markov Chain Monte Carlo (MCMC) method. We also estimated the traditional variability amplitudes ($F_{\rm var}$), which provide a model-independent measure and therefore are used to test the DRW results. We found that AGN variability characteristics are generally correlated with $(M_{\bullet},\dot{\mathscr{M}},L_{5100})$. These correlations are smooth from sub-Eddington to super-Eddington accretion AGNs, probably implying that the AGN variability may be caused by the same physical mechanism.

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

We report results of the first reverberation mapping campaign of I Zwicky 1 during $2014$-$2016$, which showed unambiguous reverberations of the broad H$\beta$ line emission to the varying optical continuum. From analysis using several methods, we obtain a reverberation lag of $\tau_{\rm H\beta}=37.2^{+4.5}_{-4.9}\,$ days. Taking a virial factor of $f_{_{\rm BLR}}=1$, we find a black hole mass of $M_{\bullet}=9.30_{-1.38}^{+1.26}\times 10^6 M_{\odot}$ from the mean spectra. The accretion rate is estimated to be $203.9_{-65.8}^{+61.0}\,L_{\rm Edd}c^{-2}$, suggesting a super-Eddington accretor, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. By decomposing {\it Hubble Space Telescope} images, we find that the stellar mass of the bulge of its host galaxy is $\log (M_{\rm bulge}/M_{\odot}) = \rm 10.92\pm 0.07$. This leads to a black hole to bulge mass ratio of $\sim 10^{-4}$, which is significantly smaller than that of classical bulges and elliptical galaxies. After subtracting the host contamination from the observed luminosity, we find that I Zw 1 follows the empirical $R_{\rm BLR}\propto L_{5100}^{1/2}$ relation.

In this work, we compiled a sample of 202 Fermi/LAT blazars with available broad line emissions. Out of the 202 sources, 66 have known Doppler factors. The correlation between ($\gamma$-ray) and broad-line emission, and that between radio and broad-line emission are investigated by removing the effects of redshift and beaming boosting for the whole sample and the subclasses, flat spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs) respectively. We obtained a strong positive correlation between ($\gamma$-ray) and broad-line emission and between radio and broad-line emission for the 202 blazars; It is worth noting that the correlation still exists after removing redshift effect. For the 66 sources with Doppler factors, there is also a positive correlation between ($\gamma$-ray) and broad-line emission after removing the Doppler factors, as well as that between radio and broad-line emission. Our analysis suggest that 1. There are strong correlations between the ($\gamma$-ray) and the broad line emission for the whole blazar sample and their subclasses. The correlations exist when the redshift effect is removed for the whole sample and their subclasses, confirming the results by Ghisellini et al. (2014) and Chen (2018). 2. For the 66 blazars with available Doppler factors, a strong correlation between the broad line emission and the Doppler factor is found. The correlation between the ($\gamma$-ray) and the broad line emission exists after the Doppler factor effect is removed. Similar results for 1 and 2 also obtained between radio and broad-line emission. 3. Our analysis suggests a robust connection between the accretion process and the jet.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

There is observational evidence that the X-ray continuum source that creates the broad fluorescent emission lines in some Seyfert Galaxies may be compact and located at a few gravitational radii above the black hole. We consider the possibility that this compact source may be powered by small scale flux tubes near the black hole that are attached to the orbiting accretion disk. As a first step, this paper investigates the salient features of black hole magnetospheres that contain small scale, disk-hole linking "closed" flux tubes, using the force-free approximation in an axisymmetric setting. We find that the extent of the closed zone is a result of the balance between the black hole spin induced twist in the closed zone and the confinement pressure of the external (open) field of the disk. The maximal extent of the closed zone, for a typical external confinement, is usually a few gravitational radii. The pressure competition between the closed zone and the external confinement could in principle lead to interesting dynamics and dissipation relevant for the compact X-ray corona.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

Most results of the reverberation monitoring of active galaxies showed a universal scaling of the time delay of the Hbeta emission region with the monochromatic flux at 5100 A, with very small dipersion. Such a scaling favored the dust-based formation mechanism of the Broad Line Region (BLR). Recent reverberation measurements showed that actually a significant fraction of objects exhibits horter lags than the previously found scaling. Here we demonstrate that these shorter lags can be explained by the old concept of scaling of the BLR size with the ionization parameter. Assuming a universal value of this parameter and universal value of the cloud density reproduces the distribution of observational points in the time delay vs. monochromatic flux plane, provided that a range of black hole spins is allowed. However, a confirmation of the new measurements for low/moderate Eddington ratio sources is strongly needed before the dust-based origin of the BLR can be excluded.

This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring season. In this approach, continuum variations are reconstructed using a damped random walk process, and BLR structure is delineated using a flexible disk-like geometry, in which BLR clouds move around the central BH with Keplerian orbits or inflow/outflow motion. The approach also includes the possibilities of anisotropic emission of BLR clouds, non-linear response of the line emission to the continuum, and different long-term trends in the continuum and emission-line variations. We implement the approach in a Bayesian framework that is apt for parallel computation and use a Markov Chain Monte Carlo technique to recover the parameters and uncertainties for the modeling, including mass of the central BH. We apply three BLR models with different prescriptions of BLR clouds distributions and find that the best model for fitting the data of Mrk 142 is a two-zone BLR model, consistent with the theoretical BLR model surrounding slim accretion disks. The best model yields a BH mass of $\log (M_\bullet/M_\odot)=6.23_{-0.45}^{+0.26}$, resulting in a virial factor of $\log f=-0.36_{-0.54}^{+0.33}$ for the full width at half maximum of the H$\beta$ line measured from the mean spectrum. The virial factors for the other measures of the H$\beta$ line width are also presented.

NGC 5548 has been well spectroscopically monitored for reverberation mapping of the central kinematics by 19 campaigns. Using the maximum entropy method in this Letter, we build up a high-quality velocity-delay map of the H$\beta$ emission line in the light curves of the continuum and the line variations observed between 2015-2016. The map shows the response strength and lags of the velocity fields of the H$\beta$ emitting regions. The velocity-delay structure of the map is generally symmetric, with strong red and blue wings at time lag $\tau \leq 15$ days, a narrower velocity distribution at $\tau \geq 15$ days, and a deficit of response in the core. This is suggestive of a disk geometry of the broad-line region (BLR). The relatively weaker H$\beta$ response at the longer lags in the red side indicates anisotropic emission from the outer part of the BLR. We also recover the velocity-delay maps of NGC 5548 from the historical data of 13 years to investigate the long-term variability of its BLR. In general, the BLR of NGC 5548 was switching between the inflow and virialized phases in the past years. The resultant maps of seven years reveal inflow signatures and show decreasing lags, indicating that the changes in the BLR size are related to the infalling BLR gas. The other four maps show potential disk signatures which are similar to our map.

As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their H\beta\ and H\gamma\ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of H\beta\ and H\gamma\ emission lines of Mrk 335 indicate that the gas of its broad-line region (BLR) is infalling. For Mrk 142, its H\beta\ and H\gamma\ lines show signatures of outflow. The H\beta\ and H\gamma\ maps of Mrk 1044 demonstrate complex kinematics -- a virialized motion accompanied by an outflow signature, and the H\beta\ map of IRAS F12397+3333 is consistent with a disk or a spherical shell. The H\beta\ maps of Mrk 486 and MCG +06-26-012 suggest the presence of an inflow and outflow, respectively. These super-Eddington accretors show diverse geometry and kinematics. Brief discussions of their BLRs are provided for each individual object.

The spectral flux density of stars can indicate their atmospheric physical properties. A detector can obtain any band flux density at the design stage. However, the band flux density is confirmed and fixed in the process of operation because of the restriction of filters. Other band flux densities cannot be obtained through the same detector. In this study, a computational model of stellar spectral flux density is established based on basic physical parameters which are effective temperature and angular parameter. The stochastic particle swarm optimization algorithm is adopted to address this issue with appropriately chosen values of the algorithm parameters. Four star catalogues are studied and consist of the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST), Wide-field Infrared Survey Explorer (WISE), Midcourse Space Experiment (MSX), and Two Micron All Sky Survey (2MASS). The given flux densities from catalogues are input parameters. Stellar effective temperatures and angular parameters are inversed using the given flux densities according to SPSO algorithm. Then the flux density is calculated according to Planck's law on the basis of stellar effective temperatures and angular parameters. The calculated flux density is compared with the given value from catalogues. It is found that the inversion results are in good agreement for all bands of the MSX and 2MASS catalogues, whereas they do not agree well in some bands of the LAMOST and WISE catalogues. Based on the results, data from the MSX and 2MASS catalogues can be used to calculate the spectral flux density at different wavelengths of given wavelength ranges. The stellar flux density is obtained and can provide data support and an effective reference for detection and recognition of stars.

Galaxy clusters host a large reservoir of diffuse plasma with radially-varying temperature profiles. The efficiency of thermal conduction in the intracluster medium (ICM) is complicated by the existence of turbulence and magnetic fields, and has received a lot of attention in the literature. Previous studies suggest that the magnetothermal instability developed in outer regions of galaxy clusters would drive magnetic field lines preferentially radial, resulting in efficient conduction along the radial direction. Using a series of spherically-symmetric simulations, here we investigate the impact of thermal conduction on the observed temperature distributions in outer regions of three massive clusters, and find that thermal conduction substantially modifies the ICM temperature profile. Within 3 Gyr, the gas temperature at a representative radius of $0.3r_{500}$ typically decreases by ~10 - 20% and the average temperature slope between $0.3r_{500}$ and $r_{500}$ drops by ~ 30 - 40%, indicating that the observed ICM would not stay in a long-term equilibrium state in the presence of thermal conduction. However, X-ray observations show that the outer regions of massive clusters have remarkably similar radially-declining temperature profiles, suggesting that they should be quite stable. Our study thus suggests that the effective conductivity along the radial direction must be suppressed below the Spitzer value by a factor of 10 or more, unless additional heating sources offset conductive cooling and maintain the observed temperature distributions. Our study provides a smoking-gun evidence for the suppression of parallel conduction along magnetic field lines in low-collisionality plasmas by kinetic mirror or whistler instabilities.

First results are presented from kinetic numerical simulations of relativistic collisionless magnetic reconnection in pair plasma that include radiation reaction from both synchrotron and inverse Compton (IC) processes, motivated by non-thermal high-energy astrophysical sources, including in particular blazars. These simulations are initiated from a configuration known as 'ABC fields' that evolves due to coalescence instability and generates thin current layers in its linear phase. Global radiative efficiencies, instability growth rates, time-dependent radiation spectra, lightcurves, variability statistics and the structure of current layers are investigated for a broad range of initial parameters. We find that the IC radiative signatures are generally similar to the synchrotron signatures. The luminosity ratio of IC to synchrotron spectral components, the Compton dominance, can be modified by more than one order of magnitude with respect to its nominal value. For very short cooling lengths, we find evidence for modification of the temperature profile across the current layers, no systematic compression of plasma density, and very consistent profiles of E.B. We decompose the profiles of E.B with the use of the Vlasov momentum equation, demonstrating a contribution from radiation reaction at the thickness scale consistent with the temperature profile.

In some low-luminosity accreting supermassive black hole systems, the supply of plasma in the funnel region can be a problem. It is believed that a local region with unscreened electric field can exist in the black hole magnetosphere, accelerating particles and producing high energy gamma-rays that can create $e^{\pm}$ pairs. We carry out time-dependent self-consistent 1D PIC simulations of this process, including inverse Compton scattering and photon tracking. We find a highly time-dependent solution where a macroscopic gap opens quasi-periodically to create $e^{\pm}$ pairs and high energy radiation. If this gap is operating at the base of the jet in M87, we expect an intermittency on the order of a few $r_g/c$, which coincides with the time scale of the observed TeV flares from the same object. For Sagittarius A* the gap electric field can potentially grow to change the global magnetospheric structure, which may explain the lack of a radio jet at the center of our galaxy.

We investigated the relationship between the accretion process and jet properties by ultilizing the VLBA and mid-infrared (MIR) data for a sample of 45 3CRR radio galaxies selected with a flux density at 178 MHz $>16.4$ Jy, 5 GHz VLA core flux density $\geq$ 7 mJy, and MIR observations. The pc-scale radio structure at 5 GHz are presented by using our VLBA observations for 21 sources in February, 2016, the analysis on the archival data for 16 objects, and directly taking the measurements for 8 radio galaxies available in literatures. The accretion mode is constrained from the Eddington ratio with a dividing value of 0.01, which is estimated from the MIR-based bolometric luminosity and the black hole masses. While most FRII radio galaxies have higher Eddington ratio than FRIs, we found that there is indeed no single correspondence between the FR morphology and accretion mode with eight FRIIs at low accretion and two FRIs at high accretion rate. There is a significant correlation between the VLBA core luminosity at 5 GHz and the Eddington ratio. Various morphologies are found in our sample, including core only, single-sided core-jet, and two-sided core-jet structures. We found that the higher accretion rate may be more likely related with the core-jet structure, thus more extended jet. These results imply that the higher accretion rates are likely able to produce more powerful jets. There is a strong correlation between the MIR luminosity at 15 $\mu$m and VLBA 5 GHz core luminosity, in favour of the tight relation between the accretion disk and jets. In our sample, the core brightness temperature ranges from $10^{9}$ to $10^{13.38}$ K with a median value of $10^{11.09}$ K indicating that systematically the beaming effect may not be significant....

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

As one of the series of papers reporting on a large reverberation mapping campaign of super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs), we present the results of 10 SEAMBHs monitored spectroscopically during 2015-2017. Six of them are observed for the first time, and have generally higher 5100 \AA\ luminosities than the SEAMBHs monitored in our campaign from 2012 to 2015; the remaining four are repeat observations to check if their previous lags change. Similar to the previous SEAMBHs, the H$\beta$ time lags of the newly observed objects are shorter than the values predicted by the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation of sub-Eddington AGNs, by factors of $\sim2-6$, depending on the accretion rate. The four previously observed objects have lags consistent with previous measurements. We provide linear regressions for the $R_{\mathrm{H\beta}}$-$L_{5100}$ relation, solely for the SEAMBH sample and for low-accretion AGNs. We find that the relative strength of Fe II and the profile of the H$\beta$ emission line can be used as proxies of accretion rate, showing that the shortening of H$\beta$ lags depends on accretion rates. The recent SDSS-RM discovery of shortened H$\beta$ lags in AGNs with low accretion rates provides compelling evidence for retrograde accretion onto the black hole. These evidences show that the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation holds only in AGNs with moderate accretion rates. At low accretion rates, it should be revised to include the effects of black hole spin, whereas the accretion rate itself becomes a key factor in the regime of high accretion rates.

The evolution of main sequence binaries that reside in the Galactic Centre can be heavily influenced by the central supermassive black hole (SMBH). Due to these perturbative effects, the stellar binaries in dense environments are likely to experience mergers, collisions or ejections through secular and/or non-secular interactions. More direct interactions with the central SMBH are thought to produce hypervelocity stars (HVSs) and tidal disruption events (TDEs). In this paper, we use N-body simulations to study the dynamics of stellar binaries orbiting a central SMBH primary with an outer SMBH secondary orbiting this inner triple. The effects of the secondary SMBH on the event rates of HVSs, TDEs and stellar mergers are investigated, as a function of the SMBH-SMBH binary mass ratio. Our numerical experiments reveal that, relative to the isolated SMBH case, the TDE and HVS rates are enhanced for, respectively, the smallest and largest mass ratio SMBH-SMBH binaries. This suggests that the observed event rates of TDEs and HVSs have the potential to serve as a diagnostic of the mass ratio of a central SMBH-SMBH binary. The presence of a secondary SMBH also allows for the creation of hypervelocity binaries. Observations of these systems could thus constrain the presence of a secondary SMBH in the Galactic Centre.

We solved the set of two-dimensional magnetohydrodynamic (MHD) equations for optically thin black hole accretion flows incorporating toroidal component of magnetic field. Following global and local MHD simulations of black hole accretion disks, the magnetic field inside the disk is decomposed into a large scale field and a fluctuating field. The effects of the fluctuating magnetic field in transferring the angular momentum and dissipating the energy are described through the usual $ \alpha $ description. We solved the MHD equations by assuming steady state and radially self-similar approximation in $ r-\theta $ plane of spherical coordinate system. We found that as the amount of magnetic field at the equatorial plane increases, the heating by the viscosity decreases. In addition, the maximum amount of the heating by the viscous dissipation is produced at the mid-plane of the disk, while that of the heating by the magnetic field dissipation is produced at the surface of the disk. Our main conclusion is that in terms of the no-outflow solution, thermal equilibrium still exists for the strong magnetic filed at the equatorial plane of the disk.

The cooling flow problem is one of the central problems in galaxy clusters, and active galactic nucleus (AGN) feedback is considered to play a key role in offsetting cooling. However, how AGN jets heat and suppress cooling flows remains highly debated. Using an idealized simulation of a cool-core cluster, we study the development of central cooling catastrophe and how a subsequent powerful AGN jet event averts cooling flows, with a focus on complex gasdynamical processes involved. We find that the jet drives a bow shock, which reverses cooling inflows and overheats inner cool core regions. The shocked gas moves outward in a rarefaction wave, which rarefies the dense core and adiabatically transports a significant fraction of heated energy to outer regions. As the rarefaction wave propagates away, inflows resume in the cluster core, but a trailing outflow is uplifted by the AGN bubble, preventing gas accumulation and catastrophic cooling in central regions. Inflows and trailing outflows constitute meridional circulations in the cluster core. At later times, trailing outflows fall back to the cluster centre, triggering central cooling catastrophe and potentially a new generation of AGN feedback. We thus envisage a picture of cool cluster cores going through cycles of cooling-induced contraction and AGN-induced expansion. This picture naturally predicts an anti-correlation between the gas fraction (or X-ray luminosity) of cool cores and the central gas entropy, which may be tested by X-ray observations.

We mapped the kinetic temperature structure of the Orion molecular cloud 1 with para-H2CO(303-202, 322-221, and 321-220) using the APEX 12m telescope. This is compared with the temperatures derived from the ratio of the NH3(2,2)/(1,1) inversion lines and the dust emission. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured averaged line ratios of para-H2CO 322-221/303-202 and 321-220/303-202. The gas kinetic temperatures derived from the para-H2CO line ratios are warm, ranging from 30 to >200 K with an average of 62 K at a spatial density of 10$^5$ cm$^{-3}$. These temperatures are higher than those obtained from NH3(2,2)/(1,1) and CH3CCH(6-5) in the OMC-1 region. The gas kinetic temperatures derived from para-H2CO agree with those obtained from warm dust components measured in the mid infrared (MIR), which indicates that the para-H2CO(3-2) ratios trace dense and warm gas. The cold dust components measured in the far infrared (FIR) are consistent with those measured with NH3(2,2)/(1,1) and the CH3CCH(6-5) line series. With dust at MIR wavelengths and para-H2CO(3-2) on one side and dust at FIR wavelengths, NH3(2,2)/(1,1), and CH3CCH(6-5) on the other, dust and gas temperatures appear to be equivalent in the dense gas of the OMC-1 region, but provide a bimodal distribution, one more directly related to star formation than the other. The non-thermal velocity dispersions of para-H2CO are positively correlated with the gas kinetic temperatures in regions of strong non-thermal motion (Mach number >2.5) of the OMC-1, implying that the higher temperature traced by para-H2CO is related to turbulence on a 0.06 pc scale. Combining the temperature measurements with para-H2CO and NH3(2,2)/(1,1) line ratios, we find direct evidence for the dense gas along the northern part of the OMC-1 10 km s$^{-1}$ filament heated by radiation from the central Orion nebula.

We study the secular dynamics of compact binaries (consisting of white dwarfs, neutron stars or black holes) with tertiary companions in hierarchical triple systems. As the inner binary (with initially negligible eccentricity) undergoes orbital decay due to gravitational radiation, its eccentricity can be excited by gravitational forcing from the tertiary. This excitation occurs when the triple system passes through an "apsidal precession resonance", when the precession rate of the inner binary, driven by the gravitational perturbation of the external companion and general relativity, matches the precession rate of the outer binary. The eccentricity excitation requires the outer companion to be on an eccentric orbit, with the mutual inclination between the inner and outer orbits less than $\sim 40^\circ$. Gravitational wave (GW) signals from the inner binary can be significantly modified as the system evolves through the apsidal precession resonance. For some system parameters (e.g., a white dwarf binary with a brown dwarf tertiary), the resonance can happen when the binary emits GWs in the $10^{-4}-10^{-1}$~Hz range (the sensitivity band of LISA).

Using the 1.56m telescope at the Shanghai Observatory (ShAO), China, we monitored two sources, BL Lac object S5 0716+714 and Flat Spectrum Radio Quasar (FSRQ) 3C 273. For S5 0716+714, we report 4969 sets of CCD (Charge-coupled Device) photometrical optical observations (1369 for V band, 1861 for R band and 1739 for I band) in the monitoring time from Dec.4, 2000 to Apr.5, 2014. For 3C 273, we report 460 observations (138 for V band, 146 for R band and 176 for I band) in the monitoring time from Mar. 28, 2006 to Apr. 9, 2014. The observations provide us with a large amount of data to analyze the short-term and long-term optical variabilities. Based on the variable timescales, we can estimate the central black hole mass and the Doppler factor. An abundance of multi-band observations can help us to analyze the relations between the brightness and spectrum. We use Gaussian fitting to analyze the intra-day light curves and obtain the intra-day variability (IDV) timescales. We use the discrete correlation function (DCF) method and Jurkevich method to analyze the quasi-periodic variability. Based on the VRI observations, we use the linear fitting to analyze the relations between brightness and spectrum. The two sources both show IDV properties for S5 0716+714. The timescales are in the range from 17.3 minutes to 4.82 hours; for 3C273, the timescale is 35.6 minutes. Based on the periodic analysis methods, we find the periods P(V) = 24.24 days, P(R)=24.12 days, P(I)=24.82 days for S5 0716+714, and P = 12.99, 21.76 yr for 3C273. The two sources displayed the "bluer-when-brighter" spectral evolution properties. S5 0716+714 and 3C 273 are frequently studied objects. The violent optical variability and IDV may come from the jet. Gaussian fitting can be used to analyze IDVs. The relations between brightness (flux density) and spectrum are strongly influenced by the frequency.

Pulsar Wind Nebulae, Blazars, Gamma Ray Bursts and Magnetars all contain regions where the electromagnetic energy density greatly exceeds the plasma energy density. These sources exhibit dramatic flaring activity where the electromagnetic energy distributed over large volumes, appears to be converted efficiently into high energy particles and gamma-rays. We call this general process magnetoluminescence. Global requirements on the underlying, extreme particle acceleration processes are described and the likely importance of relativistic beaming in enhancing the observed radiation from a flare is emphasized. Recent research on fluid descriptions of unstable electromagnetic configurations are summarized and progress on the associated kinetic simulations that are needed to account for the acceleration and radiation is discussed. Future observational, simulation and experimental opportunities are briefly summarized.

In this contribution we review the recent progress in the modeling of Pulsar Wind Nebulae (PWN). We start with a brief overview of the relevant physical processes in the magnetosphere, the wind-zone and the inflated nebula bubble. Radiative signatures and particle transport processes obtained from 3D simulations of PWN are discussed in the context of optical and X-ray observations. We then proceed to consider particle acceleration in PWN and elaborate on what can be learned about the particle acceleration from the dynamical structures called "wisps" observed in the Crab nebula. We also discuss recent observational and theoretical results of gamma-ray flares and the inner knot of the Crab nebula, which had been proposed as the emission site of the flares. We extend the discussion to GeV flares from binary systems in which the pulsar wind interacts with the stellar wind from a companion star. The chapter concludes with a discussion of solved and unsolved problems posed by PWN.

Variability is one of the extreme observational properties of BL Lacertae objects. AO 0235+164 is a well studied BL Lac through the whole electro-magnetic wavebands. In the present work, we show its optical R band photometric observations carried out during the period of Nov, 2006 to Dec. 2012 using the Ap6E CCD camera attached to the primary focus of the $\rm 70-cm$ meniscus telescope at Abastumani Observatory, Georgia. It shows a large variation of $\Delta R$ = 4.88 mag (14.19 - 19.07 mag) and a short time scale of $\Delta T_v$ = 73.5 min during our monitoring period. During the period of Dec. 2006 to Nov. 2009, we made radio observations of the source using the 25-m radio telescope at Xinjiang Astronomical Observatory. When a discrete correlation function (DCF) is adopted to the optical and radio observations, we found that the optical variation leads the radio variation by 23.2$\pm$12.9 days.

In this work, which is a continuation of Castell\'o-Mor et al. (2016), we present new X-ray and infrared (IR) data for a sample of active galactic nuclei (AGN) covering a wide range in Eddington ratio over a small luminosity range. In particular, we rigorously explore the dependence of the optical-to-X-ray spectral index $\alpha_{OX}$ and the IR-to-optical spectral index on the dimensionless accretion rate, $\dot{\mathcal{M}}=\dot{m}/\eta$ where $\dot{m}=L_{AGN}/L_{Edd}$ and $\eta$ is the mass-to-radiation conversion efficiency, in low and high accretion rate sources. We find that the SED of the faster accreting sources are surprisingly similar to those from the comparison sample of sources with lower accretion rate. In particular: i) the optical-to-UV AGN SED of slow and fast accreting AGN can be fitted with thin AD models. ii) The value of $\alpha_{OX}$ is very similar in slow and fast accreting systems up to a dimensionless accretion rate $\dot{\mathcal{M}}_{c}\sim$10. We only find a correlation between $\alpha_{OX}$ and $\dot{\mathcal{M}}$ for sources with $\dot{\mathcal{M}} > \dot{\mathcal{M}}_{c}$. In such cases, the faster accreting sources appear to have systematically larger $\alpha_{OX}$ values. iii) We also find that the torus in the faster accreting systems seems to be less efficient in reprocessing the primary AGN radiation having lower IR-to-optical spectral slopes. These findings, failing to recover the predicted differences between the SEDs of slim and thin ADs within the observed spectral window, suggest that additional physical processes or very special geometry act to reduce the extreme UV radiation in fast accreting AGN. This may be related to photon trapping, strong winds, and perhaps other yet unknown physical processes.

The evolution of main sequence binaries resided in the galactic centre is influenced a lot by the central super massive black hole (SMBH). Due to this perturbation, the stars in a dense environment are likely to experience mergers or collisions through secular or non-secular interactions. In this work, we study the dynamics of the stellar binaries at galactic center, perturbed by another distant SMBH. Geometrically, such a four-body system is supposed to be decomposed into the inner triple (SMBH-star-star) and the outer triple (SMBH-stellar binary-SMBH). We survey the parameter space and determine the criteria analytically for the stellar mergers and the tidal disruption events (TDEs). For a relative distant and equal masses SMBH binary, the stars have more opportunities to merge as a result from the Lidov-Kozai(LK) oscillations in the inner triple. With a sample of tight stellar binaries, our numerical experiments reveal that a significant fraction of the binaries, ~70 per cent, experience merger eventually. Whereas the majority of the stellar TDEs are likely to occur at a close periapses to the SMBH, induced by the outer Kozai effect. The tidal disruptions are found numerically as many as ~10 per cent for a close SMBH binary that is enhanced significantly than the one without the external SMBH. These effects require the outer perturber to have an inclined orbit (>=40 degree) relatively to the inner orbital plane and may lead to a burst of the extremely astronomical events associated with the detection of the SMBH binary.

The Large Magellanic Cloud (LMC), the closest star forming galaxy with low metallicity, provides an ideal laboratory to study star formation in such an environment. The classical dense molecular gas thermometer NH3 is rarely available in a low metallicity environment because of photoionization and a lack of nitrogen atoms. Our goal is to directly measure the gas kinetic temperature with formaldehyde toward six star-forming regions in the LMC. Three rotational transitions of para-H2CO near 218 GHz were observed with the APEX 12m telescope toward six star forming regions in the LMC. Those data are complemented by C18O 2-1 spectra. Using non-LTE modeling with RADEX, we derive the gas kinetic temperature and spatial density, using as constraints the measured para-H2CO 321-220/303-202 and para-H2CO 303-202/C18O 2-1 ratios. Excluding the quiescent cloud N159S, where only one para-H2CO line could be detected, the gas kinetic temperatures derived from the preferred para-H2CO 321-220/303-202 line ratios range from 35 to 63 K with an average of 47 K. Spatial densities of the gas derived from the paraH2CO 303-202/C18O 2-1 line ratios yield 0.4-2.9x10^5 cm^-3 with an average of 1.5x10^5 cm^-3. Temperatures derived from the para-H2CO line ratio are similar to those obtained with the same method from Galactic star forming regions and agree with results derived from CO in the dense regions of the LMC. A comparison of kinetic temperatures derived from para-H2CO with those from the dust also shows good agreement. This suggests that the dust and para-H2CO are well mixed in the studied star forming regions. A correlation between the gas kinetic temperatures derived from para-H2CO and infrared luminosity, represented by the 250um flux, suggests that the kinetic temperatures traced by para-H2CO are correlated with the ongoing massive star formation in the LMC.

For a general understanding of the physics involved in the star formation process, measurements of physical parameters such as temperature and density are indispensable. The chemical and physical properties of dense clumps of molecular clouds are strongly affected by the kinetic temperature. Therefore, this parameter is essential for a better understanding of the interstellar medium. Formaldehyde, a molecule which traces the entire dense molecular gas, appears to be the most reliable tracer to directly measure the gas kinetic temperature.We aim to determine the kinetic temperature with spectral lines from formaldehyde and to compare the results with those obtained from ammonia lines for a large number of massive clumps.Three 218 GHz transitions (JKAKC=303-202, 322-221, and 321-220) of para-H2CO were observed with the 15m James Clerk Maxwell Telescope (JCMT) toward 30 massive clumps of the Galactic disk at various stages of high-mass star formation. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured para-H2CO 322-221/303-202and 321-220/303-202 ratios. The gas kinetic temperatures derived from the para-H2CO (321-220/303-202) line ratios range from 30 to 61 K with an average of 46 K. A comparison of kinetic temperature derived from para-H2CO, NH3, and the dust emission indicates that in many cases para-H2CO traces a similar kinetic temperature to the NH3 (2,2)/(1,1) transitions and the dust associated with the HII regions. Distinctly higher temperatures are probed by para-H2CO in the clumps associated with outflows/shocks. Kinetic temperatures obtained from para-H2CO trace turbulence to a higher degree than NH3 (2,2)/(1,1) in the massive clumps. The non-thermal velocity dispersions of para-H2CO lines are positively correlated with the gas kinetic temperature. The massive clumps are significantly influenced by supersonic non-thermal motions.

Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques and data sets used, it is essentially impossible to compare the results from the literature. This problem is exacerbated by the low event rates of large solar flares. The challenges of forecasting rare events have long been recognized in the meteorology community, but have yet to be fully acknowledged by the space weather community. During the interagency workshop on "all clear" forecasts held in Boulder, CO in 2009, the performance of a number of existing algorithms was compared on common data sets, specifically line-of-sight magnetic field and continuum intensity images from MDI, with consistent definitions of what constitutes an event. We demonstrate the importance of making such systematic comparisons, and of using standard verification statistics to determine what constitutes a good prediction scheme. When a comparison was made in this fashion, no one method clearly outperformed all others, which may in part be due to the strong correlations among the parameters used by different methods to characterize an active region. For M-class flares and above, the set of methods tends towards a weakly positive skill score (as measured with several distinct metrics), with no participating method proving substantially better than climatological forecasts.

(Abridged) In this paper, multi-wavelength data are compiled for a sample of 1425 Fermi blazars to calculate their spectral energy distributions (SEDs). A parabolic function, $\log(\nu F_{\nu}) = P_1(\log\nu - P_2)^2 + P_3,$ is used for SED fitting. Synchrotron peak frequency ($\log \nu_p$), spectral curvature ($P_1$), peak flux ($\nu_{\rm p}F_{\rm \nu_p}$), and integrated flux ($\nu F_{\nu}$) are successfully obtained for 1392 blazars (461 flat spectrum radio quasars-FSRQs, 620 BL Lacs-BLs and 311 blazars of uncertain type-BCUs, 999 sources have known redshifts). Monochromatic luminosity at radio 1.4 GHz, optical R band, X-ray at 1 keV and $\gamma$-ray at 1 GeV, peak luminosity, integrated luminosity and effective spectral indexes of radio to optical ($\alpha_{\rm RO}$), and optical to X-ray ($\alpha_{\rm OX}$) are calculated. The "Bayesian classification" is employed to log$\nu_{\rm p}$ in the rest frame for 999 blazars with available redshift and the results show that 3 components are enough to fit the $\log\nu_{\rm p}$ distribution, there is no ultra high peaked subclass. Based on the 3 components, the subclasses of blazars using the acronyms of Abdo et al. (2010a) are classified, and some mutual correlations are also studied. Conclusions are finally drawn as follows: (1) SEDs are successfully obtained for 1392 blazars. The fitted peak frequencies are compared with common sources from samples available (Sambruna et al. 1996, Nieppola et al. 2006, 2008, Abdo et al. 2010a). (2) Blazars are classified as low synchrotron peak sources (LSPs) if $\log\nu_{\rm p}$(Hz) $\leq 14.0$, intermediate synchrotron peak sources (ISPs) if $14.0 < \log\nu_{\rm p}$(Hz) $\leq 15.3$, and high synchrotron peak sources (HSPs) if $\log\nu_{\rm p}$(Hz) $> 15.3$. (3) $\gamma$-ray emissions are strongly correlated with radio emissions. (...)

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short time scales. These are likely due to rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reaction. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The "flares" are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

In this work, we aim to characterise high-mass clumps with infall motions. We selected 327 clumps from the Millimetre Astronomy Legacy Team 90-GHz (MALT90) survey, and identified 100 infall candidates. Combined with the results of He et al. (2015), we obtained a sample of 732 high-mass clumps, including 231 massive infall candidates and 501 clumps where infall is not detected. Objects in our sample were classified as pre-stellar, proto-stellar, HII or photo-dissociation region (PDR). The detection rates of the infall candidates in the pre-stellar, proto-stellar, HII and PDR stages are 41.2%, 36.6%, 30.6% and 12.7%, respectively. The infall candidates have a higher H$_{2}$ column density and volume density compared with the clumps where infall is not detected at every stage. For the infall candidates, the median values of the infall rates at the pre-stellar, proto-stellar, HII and PDR stages are 2.6$\times$10$^{-3}$, 7.0$\times$10$^{-3}$, 6.5$\times$10$^{-3}$ and 5.5$\times$10$^{-3}$ M$_\odot$ yr$^{-1}$, respectively. These values indicate that infall candidates at later evolutionary stages are still accumulating material efficiently. It is interesting to find that both infall candidates and clumps where infall is not detected show a clear trend of increasing mass from the pre-stellar to proto-stellar, and to the HII stages. The power indices of the clump mass function (ClMF) are 2.04$\pm$0.16 and 2.17$\pm$0.31 for the infall candidates and clumps where infall is not detected, respectively, which agree well with the power index of the stellar initial mass function (2.35) and the cold Planck cores (2.0).

Considering features of stellar spectral radiation and survey explorers, we established a computational model for stellar effective temperatures, detected angular parameters, and gray rates. Using known stellar flux data in some band, we estimated stellar effective temperatures and detected angular parameters using stochastic particle swarm optimization (SPSO). We first verified the reliability of SPSO, and then determined reasonable parameters that produced highly accurate estimates under certain gray deviation levels. Finally, we calculated 177,860 stellar effective temperatures and detected angular parameters using the Midcourse Space Experiment (MSX) catalog data. These derived stellar effective temperatures were accurate when we compared them to known values from literatures. This research made full use of catalog data and presented an original technique for studying stellar characteristics. It proposed a novel method for calculating stellar effective temperatures and detected angular parameters, and provided theoretical and practical data support for finding radiation flow information for any band.

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$\beta$ lags ($\tau_{_{\rm H\beta}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm H\beta}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$\beta$ lags tightly depends on accretion rates, $\tau_{_{\rm H\beta}}/\tau_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $\tau_{_{R-L}}$ is the H$\beta$ lag from the normal $R_{_{\rm H\beta}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm H\beta}}=\alpha_1\ell_{44}^{\beta_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-\gamma_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 \AA\ continuum luminosity, and coefficients of $\alpha_1=(29.6_{-2.8}^{+2.7})$ lt-d, $\beta_1=0.56_{-0.03}^{+0.03}$, $\gamma_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

We investigate state transitions in black hole X-ray binaries through different parameters by using two-dimensional axisymmetric hydrodynamical simulation method. For radiative cooling in hot accretion flow, we take into account the bremsstrahlung, synchrotron and synchrotron-self Comptonization self-consistently in the dynamics. Our main result is that the state transitions occur when the accretion rate reaches a critical value $\dot M \sim 3\alpha\ \dot M_{\rm Edd}$, above which cold and dense clumpy/filamentary structures are formed, embedded within the hot gas. We argued this mode likely corresponds to the proposed two-phase accretion model, which may be responsible for the intermediate state of black hole X-ray binaries. When the accretion rate becomes sufficiently high, the clumpy/filamentary structures gradually merge and settle down onto the mid-plane. Eventually the accretion geometry transforms to a disc-corona configuration. In summary our results are consistent with the truncated accretion scenario for the state transition.

Previous works have clearly shown the existence of winds from black hole hot accretion flow and investigated their detailed properties. In extremely low accretion rate systems, the collisional mean-free path of electrons is large compared with the length-scale of the system, thus thermal conduction is dynamically important. When the magnetic field is present, the thermal conduction is anisotropic and energy transport is along magnetic field lines. In this paper, we study the effects of anisotropic thermal conduction on the wind production in hot accretion flows by performing two-dimensional magnetohydrodynamic simulations. We find that thermal conduction has only moderate effects on the mass flux of wind. But the energy flux of wind can be increased by a factor of $\sim 10$ due to the increase of wind velocity when thermal conduction is included. The increase of wind velocity is because of the increase of driving forces (e.g. gas pressure gradient force and centrifugal force) when thermal conduction is included. This result demonstrates that thermal conduction plays an important role in determining the properties of wind.

We present the results of particle-in-cell numerical pair plasma simulations of relativistic 2D magnetostatic equilibria known as the 'ABC' fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of 'X-point collapse' first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration time scale. We describe in detail the evolving structure of the dynamical current layers, and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$\beta$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$\beta$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

Jets associated with Active Galactic Nuclei (AGN) have been observed for almost a century, initially at optical and radio wavelengths. They are now widely accepted as "exhausts" produced electromagnetically by the central, spinning, massive black hole and its orbiting, accreting gas. Observations at X-ray and, especially, gamma-ray energies have transformed our understanding of how these jets evolve dynamically, accelerate electrons (and positrons) and radiate throughout the entire electromagnetic spectrum. Some new approaches to modeling the powerful and rapidly variable TeV emission observed from many blazars are sketched. Observations at the highest TeV energies, to which the High Altitude Water Cherenkov Gamma-Ray Observatory (HAWC) will contribute, promise crucial discrimination between rival models of AGN jets.

The peaks of the spectra of the accretion disks surrounding massive black holes in quasars are in the far-UV or soft X-ray band, which are usually not observed. However, in the disk corona model, the soft photons from the disk are Comptonized to high energy in the hot corona, and the hard X-ray spectra (luminosity and spectral shape) contain the information of the incident spectra from the disk. The values of black hole spin parameter $a_{\ast}$ are inferred from the spectral fitting, which spread over a large range, $\sim -0.94$ to $0.998$. We find that the inclination angles and mass accretion rates are well determined by the spectral fitting, while the results are sensitive to the accuracy of black hole mass estimates. No tight constraints on the black hole spins are achieved, if the uncertainties of black hole mass measurements are a factor of four, which are typical for the single-epoch reverberation mapping method. Recently, the accuracy of black hole mass measurement has been significantly improved to $0.2-0.4$~dex with velocity resolved reverberation mapping method (Pancoast et al., 2014). The black hole spin can be well constrained if the mass measurement accuracy is $\la 50$\%. In the accretion disk corona scenario, a fraction of power dissipated in the disk is transported into the corona, and therefore the accretion disk is thinner than a bare disk for the same mass accretion rate, because the radiation pressure in the disk is reduced. We find that the thin disk approximation, $H/R\la 0.1$, is still valid if $0.3<\dot{m}<0.5$, provided a half of the dissipated power is radiated in the corona above the disk.

First stars can only form in structures that are suitably dense, which can be parametrized by the minimum dark matter halo mass $M_{\rm min}$. $M_{\rm min}$ must plays an important role in star formation. The connection of long gamma-ray bursts (LGRBs) with the collapse of massive stars has provided a good opportunity for probing star formation in dark matter halos. We place some constraints on $M_{\rm min}$ using the latest $Swift$ LGRB data. We conservatively consider that LGRB rate is proportional to the cosmic star formation rate (CSFR) and an additional evolution parametrized as $(1+z)^{\alpha}$, where the CSFR model as a function of $M_{\rm min}$. Using the $\chi^{2}$ statistic, the contour constraints on the $M_{\rm min}$--$\alpha$ plane show that at the $1\sigma$ confidence level, we have $M_{\rm min}<10^{10.5}$ $\rm M_{\odot}$ from 118 LGRBs with redshift $z<4$ and luminosity $L_{\rm iso}>1.8\times10^{51}$ erg $\rm s^{-1}$. We also find that adding 12 high-\emph{z} $(4<z<5)$ LGRBs (consisting of 104 LGRBs with $z<5$ and $L_{\rm iso}>3.1\times10^{51}$ erg $\rm s^{-1}$) could result in much tighter constraints on $M_{\rm min}$, for which, $10^{7.7}\rm M_{\odot}<M_{\rm min}<10^{11.6}\rm M_{\odot}$ ($1\sigma$). Through Monte Carlo simulations, we estimate that future five years of Sino-French spacebased multiband astronomical variable objects monitor (\emph{SVOM}) observations would tighten these constraints to $10^{9.7}\rm M_{\odot}<M_{\rm min}<10^{11.3}\rm M_{\odot}$. The strong constraints on $M_{\rm min}$ indicate that LGRBs are a new promising tool for investigating star formation in dark matter halos.

Recent observations of the Crab Nebula (Rudy et al 2015) have maintained its reputation for high energy astrophysical enlightenment and its use as a testbed for theories of the behaviour of magnetized, relativistic plasma. In particular, new observations of the inner knot located 0.65" SE from the pulsar confirm that it is compact, elongated transversely to the symmetry axis and curved concave towards the pulsar. 60 percent polarization has been measured along the symmetry axis (Moran et al 2013). The knot does not appear to be involved in the gamma ray flares. The new observations both reinforce the interpretation of the knot as dissipation of the pulsar wind at a strong shock and challenge the details of existing models of this process. In particular, it is argued that the compactness, high polarization and curvature are difficult to reconcile with simple relativistic shock models. Alternative possibilities include deflection of the outflow ahead of the shock and spatial variation in which the knot is interpreted as a caustic. Some future observations are proposed and new theoretical investigations are suggested.

In order to understand the conditions which lead a highly magnetized, relativistic plasma to become unstable, and in such cases how the plasma evolves, we study a prototypical class of magnetostatic equilibria where the magnetic field satisfies $\nabla \times\mathbf B = \alpha \mathbf B$, where $\alpha$ is spatially uniform, on a periodic domain. Using numerical solutions we show that generic examples of such equilibria are unstable to ideal modes (including incompressible ones) which are marked by exponential growth in the linear phase. We characterize the unstable mode, showing how it can be understood in terms of merging magnetic and current structures, and explicitly demonstrate its instability using the energy principle. Following the nonlinear evolution of these solutions, we find that they rapidly develop regions with relativistic velocities and electric fields of comparable magnitude to the magnetic field, liberating magnetic energy on dynamical timescales and eventually settling into a configuration with the largest allowable wavelength. These properties make such solutions a promising setting for exploring the mechanisms behind extreme cosmic sources of gamma rays.

Observations on the high-redshift galaxies at $z>6$ imply that their ionizing emissivity is unable to fully reionize the Universe at $z\sim 6$. Either a high escape fraction of ionizing photons from these galaxies or a large population of faint galaxies below the detection limit are required. However, these requirements are somewhat in tension with present observations. In this work, we explored the combined contribution of mini-quasars and stars to the reionization of cosmic hydrogen and helium. Our model is roughly consistent with: (1) the low escape fractions of ionizing photons from the observed galaxies, (2) the optical depth of Cosmic Microwave Background (CMB) measured by the WMAP-7, and (3) the redshift of the end of hydrogen and helium reionization at $z\approx 6$ and $z\approx 3$, respectively. Neither an extremely high escape fraction nor a large population of fainter galaxies is required in this scenario. In our most optimistic model, more than $\sim20\%$ of the cosmic helium is reionized by $z\sim6$, and the ionized fraction of cosmic helium rapidly climbs to more than $50\%$ by $z\sim5$. These results may imply that better measurements of helium reionization, especially at high redshifts, could be helpful in constraining the growth of intermediate-mass black holes (IMBHs) in the early Universe, which would shed some light on the puzzles concerning the formation of supermassive black holes (SMBHs).

One of the most intriguing results from the gamma-ray instruments in orbit has been the detection of powerful flares from the Crab Nebula. These flares challenge our understanding of pulsar wind nebulae and models for particle acceleration. We report on the portion of a multiwavelength campaign using Keck, HST, and Chandra concentrating on a small emitting region, the Crab's inner knot, located a fraction of an arcsecond from the pulsar. We find that the knot's radial size, tangential size, peak flux, and the ratio of the flux to that of the pulsar are correlated with the projected distance of the knot from the pulsar. A new approach, using singular value decomposition for analyzing time series of images, was introduced yielding results consistent with the more traditional methods while some uncertainties were substantially reduced. We exploit the characterization of the knot to discuss constraints on standard shock-model parameters that may be inferred from our observations assuming the inner knot lies near to the shocked surface. These include inferences as to wind magnetization, shock shape parameters such as incident angle and poloidal radius of curvature, as well as the IR/optical emitting particle enthalpy fraction. We find that while the standard shock model gives good agreement with observation in many respects, there remain two puzzles: (a) The observed angular size of the knot relative to the pulsar--knot separation is much smaller than expected; (b) The variable, yet high degree of polarization reported is difficult to reconcile with a highly relativistic downstream flow.

Massive star-forming regions with observed infall motions are good sites for studying the birth of massive stars. In this paper, 405 compact sources have been extracted from the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) compact sources that also have been observed in the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey during Years 1 and 2. These observations are complemented with Spitzer GLIMPSE/MIPSGAL mid-IR survey data to help classify the elected star-forming clumps into three evolutionary stages: pre-stellar, proto-stellar and UCHII regions. The results suggest that 0.05 g cm$^{-2}$ is a reliable empirical lower bound for the clump surface densities required for massive-star formation to occur. The optically thick HCO$^{+}$(1-0) and HNC(1-0) lines, as well as the optically thin N$_{2}$H$^{+}$(1-0) line were used to search for infall motions toward these sources. By analyzing the asymmetries of the optically thick HCO$^{+}$(1-0) and HNC(1-0) lines and the mapping observations of HCO$^{+}$(1-0), a total of 131 reliable infall candidates have been identified. The HCO$^{+}$(1-0) line shows the highest occurrence of obvious asymmetric features, suggesting that it may be a better infall motion tracer than other lines such as HNC(1-0). The detection rates of infall candidates toward pre-stellar, proto-stellar and UCHII clumps are 0.3452, 0.3861 and 0.2152, respectively. The relatively high detection rate of infall candidates toward UCHII clumps indicates that many UCHII regions are still accreting matter. The peak column densities and masses of the infall candidates, in general, display a increasing trend with progressing evolutionary stages. However, the rough estimates of the mass infall rate show no obvious variation with evolutionary stage.

We have completed two years of photometric and spectroscopic monitoring of a large number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, we report on the result of the second phase of the campaign, during 2013--2014, and the measurements of five new H$\beta$ time lags out of eight monitored AGNs. All five objects were identified as super-Eddington accreting massive black holes (SEAMBHs). The highest measured accretion rates for the objects in this campaign are $\dot{\mathscr{M}}\gtrsim 200$, where $\dot{\mathscr{M}}= \dot{M}_{\bullet}/L_{\rm Edd}c^{-2}$, $\dot{M}_{\bullet}$ is the mass accretion rates, $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. We find that the H$\beta$ time lags in SEAMBHs are significantly shorter than those measured in sub-Eddington AGNs, and the deviations increase with increasing accretion rates. Thus, the relationship between broad-line region size ($R_{_{\rm H\beta}}$) and optical luminosity at 5100\AA, $R_{_{\rm H\beta}}-L_{5100}$, requires accretion rate as an additional parameter. We propose that much of the effect may be due to the strong anisotropy of the emitted slim-disk radiation. Scaling $R_{_{\rm H\beta}}$ by the gravitational radius of the black hole, we define a new radius-mass parameter ($Y$) and show that it saturates at a critical accretion rate of $\dot{\mathscr{M}}_c=6\sim 30$, indicating a transition from thin to slim accretion disk and a saturated luminosity of the slim disks. The parameter $Y$ is a very useful probe for understanding the various types of accretion onto massive black holes. We briefly comment on implications to the general population of super-Eddington AGNs in the universe and applications to cosmology.

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

We derive photometric redshifts (\zp) for sources in the entire ($\sim0.4$ deg$^2$) Hawaii-Hubble Deep Field-North (\hdfn) field with the EAzY code, based on point spread function-matched photometry of 15 broad bands from the ultraviolet (\bandu~band) to mid-infrared (IRAC 4.5 $\mu$m). Our catalog consists of a total of 131,678 sources. We evaluate the \zp~quality by comparing \zp~with spectroscopic redshifts (\zs) when available, and find a value of normalized median absolute deviation \sigm$=$0.029 and an outlier fraction of 5.5\% (outliers are defined as sources having $\rm |\zp - \zs|/(1+\zs) > 0.15$) for non-X-ray sources. More specifically, we obtain \sigm$=0.024$ with 2.7\% outliers for sources brighter than $R=23$~mag, \sigm$=0.035$ with 7.4\% outliers for sources fainter than $R=23$~mag, \sigm$=$0.026 with 3.9\% outliers for sources having $z<1$, and \sigm$=$0.034 with 9.0\% outliers for sources having $z>1$. Our \zp\ quality shows an overall improvement over an earlier \zp\ work that focused only on the central \hdfn\ area. We also classify each object as star or galaxy through template spectral energy distribution fitting and complementary morphological parametrization, resulting in 4959 stars and 126,719 galaxies. Furthermore, we match our catalog with the 2~Ms {\it Chandra} Deep Field-North main \xray~catalog. For the 462 matched non-stellar \xray~sources (281 having \zs), we improve their \zp~quality by adding three additional AGN templates, achieving \sigm$=0.035$ and an outlier fraction of 12.5\%. We make our catalog publicly available presenting both photometry and \zp, and provide guidance on how to make use of our catalog.

We present new Jansky Very Large Array (VLA) radio images of the Crab Nebula at 5.5 GHz, taken at two epochs separated by 6 days about two months after a gamma-ray flare in 2012 July. We find no significant change in the Crab's radio emission localized to a region of <2 light-months in radius, either over the 6-day interval between our present observations or between the present observations and ones from 2001. Any radio counterpart to the flare has a radio luminosity of <~ $2 \times 10^{-4}$ times that of the nebula. Comparing our images to one from 2001, we do however find changes in radio brightness, up to 10% in amplitude, which occur on decade timescales throughout the nebula. The morphology of the changes is complex suggesting both filamentary and knotty structures. The variability is stronger, and the timescales likely somewhat shorter, nearer the centre of the nebula. We further find that even with the excellent uv~coverage and signal-to-noise of the VLA, deconvolution errors are much larger than the noise, being up to 1.2% of peak brightness of the nebula in this particular case.

Physicists have pondered the origin of cosmic rays for over a hundred years. However the last few years have seen an upsurge in the observation, progress in the theory and a genuine increase in the importance attached to the topic due to its intimate connection to the indirect detection of evidence for dark matter. The intent of this talk is to set the stage for the meeting by reviewing some of the basic features of the entire cosmic ray spectrum from GeV to ZeV energy and some of the models that have been developed. The connection will also be made to recent developments in understanding general astrophysical particle acceleration in pulsar wind nebulae, relativistic jets and gamma ray bursts. The prospects for future discoveries, which may elucidate the origin of cosmic rays, are bright.

It was suggested that the relativistic jets in gamma-ray bursts (GRBs) are powered via the Blandford-Znajek (BZ) mechanism or the annihilation of neutrinos and anti-neutrinos from a neutrino cooling-dominated accretion flow (NDAF). The advection and diffusion of the large-scale magnetic field of a NDAF is calculated, and the external magnetic field is found to be dragged inward efficiently by the accretion flow for a typical magnetic Prandtl number P_m=1. The maximal BZ jet power can be ~10^53-10^54 erg/sec for an extreme Kerr black hole, if an external magnetic field with 10^14 Gauss is advected by the NDAF. This is roughly consistent with the field strength of the disk formed after a tidal disrupted magnetar. The accretion flow near the black hole horizon is arrested by the magnetic field if the accretion rate is below than a critical value for a given external field. The arrested accretion flow fails to drag the field inward and the field strength decays, and then the accretion re-starts, which leads to oscillating accretion. The typical timescale of such episodic accretion is in an order of one second. This can qualitatively explain the observed oscillation in the soft extend emission of short-type GRBs.

We observed the H2CO(110-111) absorption lines and H110{\alpha} radio recombination lines (RRL) toward 180 NH3 sources using the Nanshan 25-m radio telescope. In our observation, 138 sources were found to have H2CO lines and 36 have H110{\alpha} RRLs. Among the 138 detected H2CO sources, 38 sources were first detected. The detection rates of H2CO have a better correlation with extinction than with background continuum radiation. Line center velocities of H2CO and NH3 agree well. The line width ratios of H2CO and NH3 are generally larger than 1 and are similar to that of 13CO. The correlation between column densities of H2CO and extinction is better than that between NH3 and extinction. These line width relation and column density relation indicate H2CO is distributed on a larger scale than that of NH3, being similar to the regions of 13CO. The abundance ratios between NH3 and H2CO were found to be different in local clouds and other clouds.

We report on observations of the supernova remnant Cassiopeia A in the energy range from 100 MeV to 100 GeV using 44 months of observations from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. We perform a detailed spectral analysis of this source and report on a low-energy break in the spectrum at $1.72^{+1.35}_{-0.89}$ GeV. By comparing the results with models for the gamma-ray emission, we find that hadronic emission is preferred for the GeV energy range.

Context. The progenitors of short gamma-ray bursts (SGRBs) have not yet been well identified. The most popular model is the merger of compact object binaries (NS-NS/NS-BH). However, other progenitor models cannot be ruled out. The delay-time distribution of SGRB progenitors, which is an important property to constrain progenitor models, is still poorly understood. Aims. We aim to better constrain the luminosity function of SGRBs and the delay-time distribution of their progenitors with newly discovered SGRBs. Methods. We present a low-contamination sample of 16 Swift SGRBs that is better defined by a duration shorter than 0.8 s. By using this robust sample and by combining a self-consistent star formation model with various models for the distribution of time delays, the redshift distribution of SGRBs is calculated and then compared to the observational data. Results. We find that the power-law delay distribution model is disfavored and that only the lognormal delay distribution model with the typical delay tau >= 3 Gyr is consistent with the data. Comparing Swift SGRBs with T90 > 0.8 s to our robust sample (T90 < 0.8 s), we find a significant difference in the time delays between these two samples. Conclusions. Our results show that the progenitors of SGRBs are dominated by relatively long-lived systems (tau >= 3 Gyr), which contrasts the results found for Type Ia supernovae. We therefore conclude that primordial NS-NS systems are not favored as the dominant SGRB progenitors. Alternatively, dynamically formed NS-NS/BH and primordial NS-BH systems with average delays longer than 5 Gyr may contribute a significant fraction to the overall SGRB progenitors.

In this paper, we study the interaction between the supermassive binary black holes in elliptical orbit and their surrounding disk with a gap. The gap in the disk is a low density region formed due to the tidal effects of the less massive black hole. The binary we have investigated has a sub-parsec separation and is coplanar with the disk. We find that the maximum variation of the surface density in the gap reaches 50% during an orbital period. However, in other regions of the disk, the density variation is much less than 1%. Furthermore, we calculate the corresponding variation of spectral energy distribution within a period, but little variation is found. The reason for these results is that the viscosity timescale of the disk at the binary radius is much longer than the orbital period of the binary.

Observations of long gamma-ray bursts (LGRBs) offer a unique opportunity for probing the cosmic star formation history, although whether or not LGRB rates are biased tracers of star formation rate history is highly debated. Based on an extensive sample of LGRBs compiled by Robertson & Ellis (2012), we analyze various models of star formation rate and the possible effect of the evolution of cosmic metallicity under the assumption that LGRBs tend to occur in low-metallicity galaxies. The models of star formation rate tested in this work include empirical fits from observational data as well as a self-consistent model calculated in the framework of the hierarchical structure formation. Comparing with the observational data, we find a relatively higher metallicity cut of $Z\gtrsim0.6Z_{\odot}$ for the empirical fits and no metallicity cut for the self-consistent model. These results imply that there is no strong metallicity preference for the host galaxy of LGRBs, in contrast to previous work which suggest a cut of $Z\sim0.1-0.3Z_{\odot}$, and that the inferred dependencies of LGRBs on their host galaxy properties are strongly related to the specific models of star formation rate. Furthermore, a significant fraction of LGRBs occur in small dark matter halos down to $3\times10^{8}\,\mathrm{M_{\odot}}$ can provide an alternative explanation for the discrepancy between the star formation rate history and LGRB rate history.

We performed a detailed study of the high-mass clump interacting with bubble N10 based on the spectral lines $^{12}CO(3-2)$, $HCO^+(4-3)$, $N_2H^+(4-3)$ and $CH_3OH(7(0,7)-6(0,6))$ and continuum emission data at 450 $\mu$m and 850 $\mu$m released on CADC and Spitzer data. Blue-shifted optically thick line $^{12}CO (3-2)$ seems to indicate that the outer envelope of the high-mass clump is still falling toward the center. Detection of $CH_3OH(7(0,7)-6(0,6))$ suggests that a hot core has formed around YSO N10-7. And position-velocity diagram of $N_2H^+ (4-3)$ indicates the cold dense core of the clump has not been destroyed by the star formation activities. The mass of N10-7 is about 27.44 $M_\odot$. The ratio $HCO^+(4-3)/N_2H^+ (4-3)$ in the outer part of the clump is larger than that in the inner part of it. The reason may be that the CO abundance relative to $N_2H^+ (4-3)$ increased in the outer part of the high-mass clump, more $N_2H^+ (4-3)$ were converted into $HCO^+(4-3)$.

We present results from our analysis of Chandra X-ray Observatory, W. M. Keck Observatory, and Karl G. Jansky Very Large Array (VLA) images of the Crab Nebula that were contemporaneous with the gamma-ray flare of 2011 April. Despite hints in the X-ray data, we find no evidence for statistically significant variations that pinpoint the specific location of the flares within the Nebula. The Keck observations extend this conclusion to the "inner knot", i.e., the feature within an arcsecond of the pulsar. The VLA observations support this conclusion. We also discuss theoretical implications of the gamma-ray flares and suggest that the most dramatic gamma-ray flares are due to radiation-reaction-limited synchrotron emission associated with sudden, dissipative changes in the current system sustained by the central pulsar.

We calculate the structure of a standard accretion disk with corona surrounding a massive Kerr black hole in general relativistic frame, in which the corona is assumed to be heated by the reconnection of the strongly buoyant magnetic fields generated in the cold accretion disk. The emergent spectra of the accretion disk-corona systems are calculated by using the relativistic ray-tracing method. We propose a new method to calculate the emergent Comptonized spectra from the coronae. The spectra of the disk-corona systems with a modified $\alpha$-magnetic stress show that both the hard X-ray spectral index and the hard X-ray bolometric correction factor $L_{\rm bol}/L_{\rm X,2-10keV}$ increase with the dimensionless mass accretion rate, which are qualitatively consistent with the observations of active galactic nuclei (AGNs). The fraction of the power dissipated in the corona decreases with increasing black hole spin parameter $a$, which leads to lower electron temperatures of the coronas for rapidly spinning black holes. The X-ray emission from the coronas surrounding rapidly spinning black holes becomes weak and soft. The ratio of the X-ray luminosity to the optical/UV luminosity increases with the viewing angle, while the spectral shape in the X-ray band is insensitive with the viewing angle. We find that the spectral index in the infrared waveband depends on the mass accretion rate and the black hole spin $a$, which deviates from $f_\nu\propto\nu^{1/3}$ expected by the standard thin disk model.

We investigate the vertical structure of neutrino dominated accretion disks by self-consistently considering the detailed microphysics, such as the neutrino transport, vertical hydrostatic equilibrium, the conservation of lepton number, as well as the balance between neutrino cooling, advection cooling and viscosity heating. After obtaining the emitting spectra of neutrinos and antineutrinos by solving the one dimensional Boltzmann equation of neutrino and antineutrino transport in the disk, we calculate the neutrino/antineutrino luminosity and their annihilation luminosity. We find that the total neutrino and antineutrino luminosity is about $10^{54}$ ergs/s and their annihilation luminosity is about $5\times10^{51}$ ergs/s with an extreme accretion rate $10 M_{\rm {sun}}$/s and an alpha viscosity $\alpha=0.1$. In addition, we find that the annihilation luminosity is sensitive to the accretion rate and will not exceed $10^{50}$ ergs/s which is not sufficient to power the most fireball of GRBs, if the accretion rate is lower than $1 M_{\rm {sun}}$/s. Therefore, the effects of the spin of black hole or/and the magnetic field in the accretion flow might be introduced to power the central engine of GRBs.

We numerically solve the one dimensional Boltzmann equation of the neutrino and anti-neutrino transport in accretion disks and obtain the fully energy dependent and direction dependent neutrino and anti-neutrino emitting spectra, under condition that the distribution of the mass density,temperature and chemical components are given. Then, we apply the resulting neutrino and anti-neutrino emitting spectra to calculate the corresponding annihilation rate of neutrino pairs above the neutrino dominated accretion disk and find that the released energy resulting from the annihilation of neutrino pairs can not provide sufficient energy for the most energetic short gamma ray bursts whose isotropic luminosity can be as high as $10^{52}$ ergs/s unless the high temperature zone where the temperature is beyond 10 MeV can stretch over 200 km in the disk. We also compare the resulting luminosity of neutrinos and anti-neutrinos with the results from the two commonly used approximate treatment of the neutrino and anti-neutrino luminosity: the Fermi-Dirac black body limit and a simplified model of neutrino transport, i.e., the gray body model, and find that both of them overestimate the neutrino/anti-neutrino luminosity and their annihilation rate greatly. Additionally, as did in Sawyer (2003), we also check the validity of the two stream approximation, and find that it is a good approximation to high accuracy.

M87 is the first detected non-blazar extragalactic Tera-Electron-Volt (TeV) source with rapid variation and very flat spectrum in the TeV band. To explain the two-peaks in the spectral energy distribution (SED) of the nucleus of M87 which is similar to those of blazars, the most commonly adopted models are the synchrotron self-Compton scattering (SSC) models and the external inverse Compton scattering (EIC) models. Considering that there is no correlated variation in the soft band (from radio to X-ray) matching the TeV variation, and the TeV sources should not suffer from the gamma-gamma absorption due to the flat TeV spectrum, the EIC models are advantageous in modeling the TeV emission from M87. In this paper, we propose a self-consistent EIC model to explain the flat TeV spectrum of M87 within the framework of fully general relativity, where the background soft photons are from the advection-dominated accretion flow (ADAF) around the central black hole, and the high energy electrons are from the mini-jets which are powered by the magnetic reconnection in the main jet (Giannios et al. 2010). In our model, both the TeV flares observed in the years of 2005 and 2008 could be well explained: the gamma-gamma absorption for TeV photons is very low, even inside the region very close to the black hole 20Rg~50Rg; at the same region, the average EIC cooling time (~ 10^2-10^3s) is short, which is consistent with the observed time scale of TeV variation. Furthermore, we also discuss the possibility that the accompanying X-ray flare in 2008 is due to the direct synchrotron radiation of the mini-jets.

Spitzer-IRS maps of H_2 pure rotational lines from S(0) to S(7) in 3 outflows from Class 0 sources - L1448, BHR71, and NGC2071 - are presented. These lines are used, in conjunction with available rovibrational, near-IR H_2 lines, to probe the physical conditions of the warm gas between hundreds and thousands of Kelvin. We have constructed maps of the molecular hydrogen column density, ortho-to-para ratio and volume density, together with the index beta of the power law describing the distribution of gas temperature. In all three outflows, the present ortho-to-para ratio significantly deviates from the high temperature equilibrium of 3, being on average between 2.0 and 2.3. These low values, that reflect the young age of these flows, are found also in regions of relatively high temperature (~ 1000 K), likely indicating that shocks are occurring in a time shorter than that needed for a complete para to ortho conversion. Density maps indicate upper limits close to LTE conditions, i.e. between 10^6-10^7 cm^-3; moreover we demonstrate, based on the detections of HD emission spots (R(3) and R(4) lines), that a density stratification does exist, with the low density components (10^4-10^5 cm^-3) associated with the coldest gas. The beta index is found in all flows to be above 3.8: this value is consistent with predictions of multiple C-type bow shocks with a range of velocities, some of which insufficient to achieve the temperature at which H_2 partially dissociates. The contribution of H_2 to the total cooling is quantitatively similar to that of other abundant molecules emitting in the far-infrared, such as water and CO; moreover, the luminosity radiated away is comparable with the estimated kinetic energy of the swept-out outflow: this supports the scenario in which the shock from which H_2 is emitted is also capable of accelerating the molecular outflow driven at the shock working surface.

We have studied the shock-excited molecular regions associated with four supernova remnants (SNRs) - IC443C, W28, W44 and 3C391 - and two Herbig-Haro objects, HH7 and HH54, using Spitzer's Infrared Spectrograph (IRS). The physical conditions within the observed areas are inferred from spectroscopic data obtained from IRS and from SWS and LWS onboard ISO, together with photometric data from Spitzer's Infrared Array Camera (IRAC). Adopting a power-law distribution for the gas temperature in the observed region, the H2 S(0) to S(7) spectral line maps obtained with IRS were used to constrain the gas density, yielding estimated n(H2) in the range 2-4*10^3 cm^-3. The excitation of H2 S(9) to S(12) and high-J CO pure rotational lines, however, require environments several times denser. The inconsistency among the best-fit densities estimated from different species can be explained by density fluctuations within the observed regions. The best-fit power-law index b is smaller than the value 3.8 predicted for a paraboloidal C-type bow shock, suggesting that the shock front has a "flatter" shape than that of a paraboloid. The best-fit parameters for SNRs and Herbig-Haro objects do not differ significantly between the two classes of sources, except that for the SNRs the ortho-to-para ratio (OPR) of hot gas (T> 1000 K) is close to the LTE value 3, while for HH7 and HH54 even the hottest gas exhibits an OPR smaller than 3. Finally, we mapped the physical parameters within the regions observed with IRS and found that the mid-lying H2 emissions - S(3) to S(5) - tend to trace the hot component of the gas, while the intensities of S(6) and S(7) are more sensitive to the density of the gas compared to S(3) to S(5).

Both types of long and short gamma ray bursts involve a stage of a hyper-Eddington accretion of hot and dense plasma torus onto a newly born black hole. The prompt gamma ray emission originates in jets at some distance from this 'central engine' and in most events is rapidly variable, having a form of spikes and subpulses. This indicates at the variable nature of the engine itself, for which a plausible mechanism is an internal instability in the accreting flow. We solve numerically the structure and evolution of the neutrino-cooled torus. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. The models are calculated for both Schwarzschild and Kerr black holes. We find that for sufficiently large accretion rates (> 10 Msun/s for non-rotating black hole, and >1 Msun/s for rotating black hole, depending on its spin), the inner regions of the disk become opaque, while the helium nuclei are being photodissociated. The sudden change of pressure in this region leads to the development of a viscous and thermal instability, and the neutrino pressure acts similarly to the radiation pressure in sub-Eddington disks. In the case of rapidly rotating black holes, the instability is enhanced and appears for much lower accretion rates. We also find the important and possibly further destabilizing role of the energy transfer from the rotating black hole to the torus via the magnetic coupling.

We present an analysis of Spitzer-IRS spectroscopic maps of the L1157 protostellar outflow in the H2 pure-rotational lines from S(0) to S(7). The aim of this work is to derive the physical conditions pertaining to the warm molecular gas and study their variations within the flow. The mid-IR H2 emission follows the morphology of the precessing flow, with peaks correlated with individual CO clumps and H2 2.12{\mu}m ro-vibrational emission. More diffuse emission delineating the CO cavities is detected only in the low-laying transitions, with J(lower) less or equal to 2. The H2 line images have been used to construct 2D maps of N(H2), H2 ortho-to-para ratio and temperature spectral index beta, in the assumption of a gas temperature stratification where the H2 column density varies as T^(beta). Variations of these parameters are observed along the flow. In particular, the ortho-to-para ratio ranges from 0.6 to 2.8, highlighting the presence of regions subject to recent shocks where the ortho-to-para ratio has not had time yet to reach the equilibrium value. Near-IR spectroscopic data on ro-vibrational H2 emission have been combined with the mid-IR data and used to derive additional shock parameters in the brightest blue- and red-shifted emission knots. A high abundance of atomic hydrogen (H/H2 about 0.1-0.3) is implied by the observed H2 column densities, assuming n(H2) values as derived by independent SiO observations. The presence of a high fraction of atomic hydrogen, indicates that a partially-dissociative shock component should be considered for the H2 excitation in these localized regions. However, planar shock models, either of C- or J-type, are not able to consistently reproduce all the physical parameters derived from our analysis of the H2 emission. Globally, H2 emission contributes to about 50% of the total shock radiated energy in the L1157 outflow.

The accelerating expansion of the Universe at the present stage is a process that will change the spectrum of relic gravitational waves. Here we present a numerical calculation for the power spectrum of relic gravitational waves in the accelerating Universe. The results show that although the overall features of the power spectrum is similar to that in the non-accelerating models, the amplitude is smaller by an order of $10^{-1}$. We also find that the spectrum is very sensitive to the index $\beta$ of the inflationary expansion with the scale factor $a(\tau) \propto |\tau|^{1+\beta}$. With increase of $\beta$, the resulting spectrum is tilted to be flatter with more power on high frequencies, and the sensitivity of the second science run of the LIGO detectors puts a restriction on the parameter $\beta < -1.8$. The influence of reheating following the inflation has been examined.

Based on the steady state equilibrium condition for neutron-proton-electron-positron gas in the neutrino-driven wind from protoneutron star, we estimate the initial electron fraction in the wind in a simple and effective way. We find that the condition in the wind might be propriate for the r-process nucleosynthesis.

Based on a new estimation of their thickness, the global properties of relativistic slim accretion disks are investigated in this work. The resulting emergent spectra are calculated using the relativistic ray-tracing method, in which we neglect the self-irradiation of the accretion disk. The angular dependence of the disk luminosity, the effects of the heat advection and the disk thickness on the estimation of the black hole spin are discussed. Compare to the previous works, our improvements are that we use the self-consistent disk equations and we consider the disk self-shadowing effect. We find that at the moderate accretion rate, the radiation trapped in the outer region of the accretion disks will escape in the inner region of the accretion disk and contribute to the emergent spectra. At the high accretion rate, for the large inclination and large black hole spin, both the disk thickness and the heat advection have significant influence on the emergent spectra. Consequently, these effects will influence the measurement of the black hole spin based on the spectra fitting and influence the angular dependence of the luminosity. For the disks around Kerr black holes with $a=0.98$, if the disk inclination is greater than $60^\circ$, and their luminosity is beyond 0.2 Eddington luminosity, the spectral model which is based on the relativistic standard accretion disk is no longer applicable for the spectra fitting. We also confirm that the effect of the self-shadowing is significantly enhanced by the light-bending, which implies that the non-relativistic treatment of the self-shadowing is inaccurate. According to our results, the observed luminosity dependence of the measured spin suggests that the disk self-shadowing significantly shapes the spectra of GRS 1915+105, which might lead to the underestimation of the black hole spin for the high luminosity states.

We report on the third phase of our study of the neutrino-cooled hyperaccreting torus around a black hole that powers the jet in Gamma Ray Bursts. We focus on the influence of the black hole spin on the properties of the torus. The structure of a stationary torus around the Kerr black hole is solved numerically. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. We find, that in the case of rapidly rotating black holes, the thermal instability discussed in our previous work is enhanced and develops for much lower accretion rates. We also find the important role of the energy transfer from the rotating black hole to the torus, via the magnetic coupling.

We report the results of spectroscopic mapping observations carried out toward protostellar outflows in the BHR71, L1157, L1448, NGC 2071, and VLA 1623 molecular regions using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, provide detailed maps of the 8 lowest pure rotational lines of molecular hydrogen and of the [SI] 25.25 micron and [FeII] 26.0 micron fine structure lines. The molecular hydrogen lines, believed to account for a large fraction of the radiative cooling from warm molecular gas that has been heated by a non-dissociative shock, allow the energetics of the outflows to be elucidated. Within the regions mapped towards these 5 outflow sources, total H2 luminosities ranging from 0.02 to 0.75 L(solar) were inferred for the sum of the 8 lowest pure rotational transitions. By contrast, the much weaker [FeII] 26.0 micron fine structure transition traces faster, dissociative shocks; here, only a small fraction of the fast shock luminosity emerges as line radiation that can be detected with Spitzer/IRS.

By using an average heavy nuclei model, the effects of the electron screening on electron capture (EC) in core-collapse supernovae are investigated. A one-dimension code based on the Ws15$M_{\sun}$ progenitor model is utilized to test the effects of electron screening during the collapsing process. The results show that, at high densities, the effects of EC on electron capture becomes significant. During the collapsing stage, the EC rate is decreased, the collapse timescale is prolonged and the leakage of the neutrino energy is increased. These effects lead to an appreciable decrease in the initial energy of the bounce shock wave. The effects of electron screeening in the other progenitor models are also discussed.

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.

The fast variability of energetic TeV photons from the center of M87 has been detected, offering a new clue to estimate spins of supermassive black holes (SMBHs). We extend the study of Wang et al. (2008) by including all of general relativistic effects. We numerically solve the full set of relativistic hydrodynamical equations of the radiatively inefficient accretion flows (RIAFs) and then obtain the radiation fields around the black hole. The optical depth of the radiation fields to TeV photons due to pair productions are calculated in the Kerr metric. We find that the optical depth strongly depends on: (1) accretion rates as $\tautev\propto \dot{M}^{2.5-5.0}$; (2) black hole spins; and (3) location of the TeV source. Jointly considering the optical depth and the spectral energy distribution radiated from the RIAFs, the strong degeneration of the spin with the other free parameters in the RIAF model can be largely relaxed. We apply the present model to M87, wherein the RIAFs are expected to be at work, and find that the minimum specific angular momentum of the hole is $a\sim0.8$. The present methodology is applicable to M87-like sources with future detection of TeV emissions to constrain the spins of SMBHs.

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 present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities and the predictions of detailed shock models indicates that the emission arises in both slow (12 - 24 km/s) and fast (36 - 53 km/s) C-type shocks with an initial ortho-to-para ratio of ~ 1. The present H2 ortho-to-para ratio exhibits a large degree of spatial variations. In the post-shocked gas, it is usually about 2, i.e. close to the equilibrium value (~ 3). However, around at least two outflows, we observe a region with a much lower (~ 0.5) ortho-to-para ratio. This region probably corresponds to gas which has been heated-up recently by the passage of a shock front, but whose ortho-to-para has not reached equilibrium yet. This, together with the low initial ortho-to-para ratio needed to reproduce the observed emission, provide strong evidence that H2 is mostly in para form in cold molecular clouds. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and thus can be used to ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for core disruption.

After considering the effects of negative feedback on the process of star formation, we explore the relationship between star formation process and the associated feedback, by investigating how the mechanical feedback from supernovae(SNe) and radiative feedback from luminous objects regulate the star formation rate and therefore affect the cosmic reionization.Based on our present knowledge of the negative feedback theory and some numerical simulations, we construct an analytic model in the framework of the Lambda cold dark matter model. In certain parameter regions, our model can explain some observational results properly. In large halos(T_vir>10000 K), both mechanical and radiative feedback have a similar behavior: the relative strength of negative feedback reduces as the redshift decreases. In contrast, in small halos (T_vir<10000 K$) that are thought to breed the first stars at early time, the radiative feedback gets stronger when the redshift decreases. And the star formation rate in these small halos depends very weakly on the star-formation efficiency. Our results show that the radiative feedback is important for the early generation stars. It can suppress the star formation rate considerably. But the mechanical feedback from the SNe explosions is not able to affect the early star formation significantly. The early star formation in small-halo objects is likely to be self-regulated. The radiative and mechanical feedback dominates the star formation rate of the PopII/I stars all along. The feedback from first generation stars is very strong and should not be neglected. However, their effects on the cosmic reionization are not significant, which results in a small contribution to the optical depth of Thomson scattering.

We report the results of spectroscopic mapping observations carried out toward a ~1 min x 1 min region within the northern lobe of the outflow from NGC 2071 using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations covered the 5.2-37 um spectral region and have led to the detection of a number of ionic, atomic, and molecular lines, including fine-structure emission of Si+, Fe+, S++, S, the S(0)-S(7) pure rotational lines of H2, the R(3) and R(4) transitions of HD, and at least 11 transitions of H2O. In addition, the 6.2, 7.4, 7.6, 7.9, 8.6 and 11.3 um PAH emission bands were also observed and several transitions of OH were tentatively detected. Most of the detected line transitions were strong enough to map including, for the first time, three transitions of hot H2O. We find that: (1) the water emission is extended; (2) the extended emission is aligned with the outflow; and, (3) the spatial distribution of the water emission generally follows that observed for H2. Based on the measured line intensities, we derive an HD abundance relative to H2 of 1.1-1.8 10^-5 and an H2O number density of 12-2 cm^3. The H2 density in the water-emitting region is not well constrained by our observations, but is likely between 3 10^4 and 10^6 cm^3, yielding an H2O abundance relative to H2 of between 2 10^-5 and 6 10^-4. Future observations planned for the Herschel Space Observatory should greatly improve the density estimate, and thus our knowledge of the H2O abundance, for the water-emitting regions reported here. Finally, we note a possible departure from the H2O ortho-to-para ratio of 3:1 expected for water formed in hot post-shocked gas, suggesting that a significant fraction of the water vapor we detect may arise from H2O sputtered from cold dust grains.

Photometric maps, obtained with Spitzer's Infrared Array Camera (IRAC), can provide a valuable probe of warm molecular hydrogen within the interstellar medium. IRAC maps of the supernova remnant IC443, extracted from the Spitzer archive, are strikingly similar to spectral line maps of the H2 pure rotational transitions that we obtained with the Infrared Spectrograph (IRS) instrument on Spitzer. IRS spectroscopy indicates that IRAC Bands 3 and 4 are indeed dominated by the H2 v=0-0 S(5) and S(7) transitions, respectively. Modeling of the H2 excitation suggests that Bands 1 and 2 are dominated by H2 v=1-0 O(5) and v=0-0 S(9). Large maps of the H2 emission in IC433, obtained with IRAC, show band ratios that are inconsistent with the presence of gas at a single temperature. The relative strengths of IRAC Bands 2, 3, and 4 are consistent with pure H2 emission from shocked material with a power-law distribution of gas temperatures. CO vibrational emissions do not contribute significantly to the observed Band 2 intensity. Assuming that the column density of H2 at temperatures T to T+dT is proportional to T raised to the power -b for temperatures up to 4000 K, we obtained a typical estimate of 4.5 for b. The power-law index, b, shows variations over the range 3 to 6 within the set of different sight-lines probed by the maps, with the majority of sight-lines showing b in the range 4 to 5. The observed power-law index is consistent with the predictions of simple models for paraboloidal bow shocks.

Optically thin two-temperature accretion flows may be thermally and viscously stable, but acoustically unstable. Here we propose that the O-mode instability of a cooling-dominated optically thin two-temperature inner disk may explain the 23-day quasi-periodic oscillation (QPO) period observed in the TeV and X-ray light curves of Mkn~501 during its 1997 high state. In our model the relativistic jet electrons Compton upscatter the disk soft X-ray photons to TeV energies, so that the instability-driven X-ray periodicity will lead to a corresponding quasi-periodicity in the TeV light curve and produce correlated variability. We analyse the dependence of the instability-driven quasi-periodicity on the mass (M) of the central black hole, the accretion rate ($\rm{\dot{M}}$) and the viscous parameter ($\alpha$) of the inner disk. We show that in the case of Mkn~501 the first two parameters are constrained by various observational results, so that for the instability occurring within a two-temperature disk where $\alpha=0.05-1.0$, the quasi-period is expected to lie within the range of 8 to 100 days, as indeed the case. In particular, for the observed 23-day QPO period our model implies a viscosity coefficient $\alpha \leq 0.28$, a sub-Eddington accretion rate $\dot{M} \simeq 0.02 \dot{M}_{\rm Edd}$ and a transition radius to the outer standard disk of $r_0 \sim 60 r_g$, and predicts a period variation $\delta P/P \sim 0.23$ due to the motion of the instability region.

We investigate the properties and evolution of accretion tori formed after the coalescence of two compact objects. At these extreme densities and temperatures, the accreting torus is cooled mainly by neutrino emission produced primarily by electron and positron capture on nucleons (beta reactions). We solve for the disc structure and its time evolution by introducing a detailed treatment of the equation of state which includes photodisintegration of helium, the condition of beta-equilibrium, and neutrino opacities. We self-consistently calculate the chemical equilibrium in the gas consisting of helium, free protons, neutrons and electron-positron pairs and compute the chemical potentials of the species, as well as the electron fraction throughout the disc. We find that, for sufficiently large accretion rates (> 10 solar masses per second), the inner regions of the disk become opaque and develop a viscous and thermal instability. The identification of this instability might be relevant for GRB observations.

In this paper, we used the database of the university of Michigan Radio Astronomy Observatory (UMRAO) at three (4.8 GHz, 8.0 GHZ, and 14.5 GHz) radio frequency to analyze the radio light curves by the power spectral analysis method in search of possible periodicity. The analysis results showed that the radio sources display astrophysically meaningful periodicity ranging from 2.2 to 20.8 years in their light curves at the three frequencies. We also calculated the variability parameters and investigated the correlations between the variability parameter and the flux density. For the variability parameters, we found that the parameters at higher frequency are higher than those in the lower frequency. In addition, the variability parameters of BL Lacertae objects are larger than those of flat-spectrum radio quasars. suggesting that they are more variable than flat spectrum radio quasars.

In this paper, using the database of the university of Michigan Radio Astronomy Observatory (UMRAO) at three (4.8 GHz, 8 GHZ, and 14.5 GHz) radio frequencies, we studied the polarization properties for 47 BL Lacertae objects(38 radio selected BL Lacertae objects, 7 X-ray selected BL Lacertae, and two inter-middle objects (Mkn 421 and Mkn 501), and found that (1) The polarizations at higher radio frequency is higher than those at lower frequency, (2) The variability of polarization at higher radio frequency is higher than those at lower frequency, (3) The polarization is correlated with the radio spectral index, and (4) The polarization is correlated with core-dominance parameter for those objects with known core-dominance parameters suggesting that the relativistic beaming could explain the polarization characteristic of BL Lacs.

We reported the result of long term optical variability of the blazar 4C 29.45 (QSO 1156+295, Ton 599), carried out optical photometric observations in Johnson V, Cousins RI passbands during April 1997 to March 2002 using the 1.56 meter telescope of the Shanghai Astronomical Observatory (SHAO) at Sheshan, China, compiled the post-1974 optical photometric data of the blazar by combining our new observations with the published optical data, and found maximum variations in different passbands: $\Delta$U = 4.41 mag, $\Delta$B = 5.55 mag, $\Delta$V = 4.53 mag, $\Delta$R = 5.80 mag, and $\Delta$I = 5.34 mag. The average color indices are: U$-$B = $-0.54\pm$0.18 mag, B$-$V = 0.56$\pm$0.21 mag, B$-$R = 0.93$\pm$0.18 mag, B$-$I = 1.51$\pm$0.24 mag, V$-$R = 0.44$\pm$0.15 mag and V$-$I = 1.03$\pm$0.23 mag. The post-1974 data give us an excellent opportunity to search for the existence of possible periodicity in the light curve. In search for periodicity in the R passband light curve, we performed Jurkevich test and power spectral (Fourier) analysis methods, and CLEANest algorithms to remove false signals. We found possible periods of 3.55 and 1.58 years. The possible mechanism for the periodic variability was discussed.

The reaction rates of the beta processes for all particles at arbitrary degeneracy are derived, and an {\it analytic} steady state equilibrium condition $\mu_n=\mu_p+2\mu_e$ which results from the equality of electron and positron capture rates in the hot electron-positron plasma with nucleons is also found, if the matter is transparent to neutrinos. This simple analytic formula is valid only if electrons are nondegenerate or mildly degenerate, which is generally satisfied in the hot electron-positron plasma. Therefore, it can be used to efficiently determine the steady state of the hot matter with plenty of positrons. Based on this analytic condition, given the baryon number density and the temperature, if the nucleons are nondegenerate, only one algebraic equation for determining the electron fraction is obtained, which shows the great advantage of the analytic equilibrium condition.

At the extreme densities and temperatures typical of the central engine of GRBs, the accreting torus is cooled mainly by advection and by neutrino emission. The latter process is dominated by electron and positron capture onto nucleons ($\beta$ reactions). We calculate the reaction rates and the nuclear composition of matter, assuming that the torus consists of helium, eletron-positron pairs, free neutrons and protons. After determining the equation of state and solving for the disk structure for a given initial accretion rate, we subsequently follow its time evolution. We find that, for accretion rates of the order of $10 M_{\odot}$/s, likely typical for the early stages of the accretion event, the disk becomes unstable, giving rise to variable energy output. This instability may play an important role for producing internal shocks.

Quasars represent a brief phase in the life-cycle of most massive galaxies, but the evolutionary connection between central black holes and their host galaxies remains unclear. While quasars are active and shining brighter than the Compton-limit luminosity, their radiation heats the surrounding medium to the Compton temperature, forming Compton spheres extending to the Str\"omgren radius of Fe$^{26+}$/He$^{2+}$. After the quasars shut off, their ``afterglow'' can be detected through three signatures: (1) an extended X-ray envelope, with a characteristic temperature of $\sim 3\times 10^7$ K; (2) Ly$\alpha$ and Ly$\beta$ lines and the $K-$edge of Fe$^{26+}$; and (3) nebulosity from hydrogen and helium recombination emission lines. We discuss the possibility of detecting these signatures using {\em Chandra}, the planned {\em XEUS} mission, and ground-based optical telescopes. The luminosity and size of quasar afterglows can be used to constrain the lifetime of quasars.

Observations have recently indicated that the Universe at the present stage is in an accelerating expansion, a process that has great implications. We evaluate the spectrum of relic gravitational waves in the current accelerating Universe and find that there are new features appearing in the resulting spectrum as compared to the decelerating models. In the low frequency range the peak of spectrum is now located at a frequency $\nu_E \simeq (\frac{\Omega_m}{\Omega_{\Lambda}})^{1/3} \nu_H$, where $\nu_H$ is the Hubble frequency, and there appears a new segment of spectrum between $\nu_E$ and $\nu_H$. In all other intervals of frequencies $\geq \nu_H$, the spectral amplitude acquires an extra factor $\frac{\Omega_m}{\Omega_{\Lambda}}$, due to the current acceleration, otherwise the shape of spectrum is similar to that in the decelerating models. The recent WMAP result of CMB anisotropies is used to normalize the amplitude for gravitational waves. The slope of the power spectrum depends sensitively on the scale factor $a(\tau) \propto |\tau|^{1+\beta}$ during the inflationary stage with $\beta=-2$ for the exact de Sitter space. With the increasing of $\beta$, the resulting spectrum is tilted to be flatter with more power on high frequencies, and the sensitivity of the second science run of the LIGO detectors puts a restriction on the parameter $\beta \leq -1.8$. We also give a numerical solution which confirms these features.

The existence of black holes remains open to doubt until other conceivable options are excluded. With this motivation, we consider a model of a compact star in which most of the mass consists of dark particles of some kind, and a small fraction of the mass is in the form of ordinary nucleonic gas. The gas does not interact with the dark matter other than via gravity, but collects at the center as a separate fermionic fluid component. Depending on whether the dark mass is made of fermions or bosons, the objects may be called fermion-fermion stars or boson-fermion stars, respectively. For appropriate choices of the mass of the dark matter particles, these objects are viable models of black hole candidates in X-ray binaries. We consider models with a dark mass of 10 solar masses and a range of gas mass from 10^{-6} to nearly one solar mass, and analyse the bursting properties of the models when they accrete gas. We show that all the models would experience thermonuclear Type I X-ray bursts at appropriate mass accretion rates. Since no Type I bursts have been reported from black hole candidates, the models are ruled out. The case for identifying black hole candidates in X-ray binaries as true black holes is thus strengthened.

The exact reaction rates of beta processes for all particles at arbitrary degeneracy are derived, and an analytic beta-equilibrium condition $\mu_n=\mu_p+2\mu_e$ for the hot electron-positron plasma with nucleons is found, if the matter is transparent to neutrinos. This simple analytic formula is valid only if electrons are nondegenerate, which is generally satisfied in the hot electron-positron plasma. Therefore, it can be used to efficiently determine the steady state of the hot matter with positrons.

We study the rotational evolution of a protoneutron star with hyperons and nucleons or solely nucleons in its core due to the escape of the trapped neutrinos. It is found that at the early stage of its evolution, the stellar crust contracts significantly, consequently the star spins up. At the late stage, a the protoneutron star with hyperons, it keeps shrinking and spinning up till all the trapped neutrinos escape. Consequently, the distribution of the stellar initial spin periods is skewed toward shorter periods. For a protoneutron star with only nucleons, the expansion of its core dominates, and the stellar rotation slows down. After the neutrinos escape, the range of the spin periods is narrower than the initial one, but the distribution is still nearly uniform. If the hyperonic star is metastable, its rotational frequency accelerates distinguishedly before it collapses to a black hole.

There is increasing evidence for the presence of an optically thin advection-dominated accretion flow (ADAF) in low luminosity active galactic nuclei and radio-loud quasars. The present paper is devoted to explore the fate of a blob ejected from an ADAF, and to discuss its observational consequences. It is inevitable for the ejected blob to drastically expand into its surroundings. Consequently, it is expected that a group of relativistic electrons should be accelerated, which may lead to nonthermal flares, since a strong shock will be formed by the interaction between the blob and its surroundings. Then the blob cools down efficiently, leading to the appearance of recombination lines about $10^5$s after its ejection from an ADAF. We apply this model to NGC 4258 for some observational prediction, and to PKS 2149--306 for the explanation of observational evidence. Future simultaneous observations of recombination X-ray lines and continuum emission are highly desired to test the present model.

Recent observations show strong evidence for the presence of Doppler-shifted emission lines in the spectrum of both black hole candidates and active galactic nuclei. These lines are likely to originate from relativistic outflows (or jets) in the vicinity of the central black hole. Consequently, the profile of such a line should be distorted by strong gravitational effects near the black hole, as well as special relativistic effects. In this paper, we present results from a detailed study on how each process affects the observed line profile. We found that the profile is sensitive to the intrinsic properties of the jets (Lorentz factor, velocity profile, and emissivity law), as well as to the spin of the black hole and the viewing angle (with respect to the axis of the jets). More specifically, in the case of approaching jets, an intrisically narrow line (blue-shifted) is seen as simply broadened at small viewing angles, but it shows a doubly peaked profile at large viewing angles for extreme Kerr black holes (due to the combination of gravitational focusing and Doppler effects); the profile is always singly peaked for Schwarzschild black holes. For receding jets, however, the line profile becomes quite complicated owing to complicated photon trajectories. To facilitate comparison with observations, we searched a large parameter space to derive representative line profiles. We show the results and discuss how to use emission lines as a potential tool for probing the inner region of a black hole jet system.

The difference between the Hipparcos open cluster distance scale and that obtained from main sequence fitting is examined. The two color main sequence fitting technique of Pinsonneault et al. (1998) is extended to NGC 2516, NGC 6475, and NGC 6633. The error sources for main sequence fitting are examined, and possible evidence for systematic errors in the Hipparcos parallaxes is discussed.

We have obtained high-resolution spectra of 89 M dwarf members of the Pleiades and Hyades and have derived radial velocities, H-alpha equivalent widths, and spectroscopic rotational velocities for these stars. Typical masses of the newly-observed Pleiades and Hyades stars are ~ 0.4 M_{\sun} and ~ 0.2 M_{\sun}, respectively. We combine our new observations with previously published data to explore the rotational evolution of young stars with M < 0.4 M_\sun. The average rotation rate in the Hyades (age 600 Myr) is about 0.4 that of the Pleiades (110 Myr), and the mean equivalent widths of H-alpha are also lower. As found in previous studies, the correlation between rotation and chromospheric activity is identical in both clusters, implying that the lower activity in the Hyades is a result of the lower rotation rates. We show that a simple scaling of the Pleiades rotational distribution for M \leq 0.4 M_{\sun}, corrected for the effects of structural evolution, matches that of the Hyades if the average angular momentum loss from the Pleiades to the Hyades age is factor of \approx 6. This suggests that the distribution of initial angular momenta and disk-locking lifetimes for the lowest mass stars was similar in both clusters. We argue that this result provides further evidence for a saturation of the angular momentum loss rate at high rotational velocities.