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}$).

In this work, Machine Learning (ML) methods are used to efficiently identify the unassociated sources and the Blazar Candidate of Uncertain types (BCUs) in the Fermi-LAT Third Source Catalog (3FGL). The aims are twofold: 1) to distinguish the Active Galactic Nuclei (AGNs) from others (non-AGNs) in the unassociated sources; 2) to identify BCUs into BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs). Two dimensional reduction methods are presented to decrease computational complexity, where Random Forest (RF), Multilayer Perceptron (MLP) and Generative Adversarial Nets (GAN) are trained as individual models. In order to achieve better performance, the ensemble technique is further explored. It is also demonstrated that grid search method is of help to choose the hyper-parameters of models and decide the final predictor, by which we have identified 748 AGNs out of 1010 unassociated sources, with an accuracy of 97.04%. Within the 573 BCUs, 326 have been identified as BL Lacs and 247 as FSRQs, with an accuracy of 92.13%.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

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

Active galactic nuclei (AGNs) can be divided into two major classes, namely radio-loud and radio-quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff-Riley type I and type II (FRI&II) radio galaxies. Following our previous work, we present a latest sample of 966 sources with measured radio flux densities of the core and extended components. The sample includes 83 BL Lacs, 473 FSRQs, 101 Seyferts, 245 galaxies, 52 FRIs&IIs and 12 unidentified sources. We then calculate the radio core-dominance parameters and spectral indices and study their relationship. Our analysis shows that the core-dominance parameters and spectral indices are quite different for different types of sources. We also corroborate that the correlation between core-dominance parameter and radio spectral index extends over all the sources in a large sample presented.

Aims. We present a sample of 4388 AGNs with available radio core-dominance parameters defined as the ratio of the core flux densities to extended ones, namely, R = Score/Sext., which includes 630 Fermi-detected AGNs respect to the catalog of 4FGL, the fourth Fermi Large Area Telescope source catalog, and the rest of them are non-Fermi-detected AGNs. In our sample, 584 blazars are Fermi-detected and 1310 are not, and also consists of other subclasses such as Seyfert, Fanaroff-Riley I/II and normal galaxies. We investigate various different properties between Fermi-detected AGNs and non-Fermi-detected ones by using the core-dominance parameters as the previous study has shown that R is a good indication of beaming effect. Methods. We then calculate the radio spectral indices for whole sample and adopt {\gamma}-ray photon indices for Fermi AGNs from 4FGL catalog for discussing their different performances on different subclasses, and obtain the relation between core-dominance parameters and radio spectral indices for both Fermi and non-Fermi sources according to the two components model on radio band, which are consistent with our previous study. Results. We found that the core-dominance parameters and radio spectral indices are quite different for different subclasses of AGNs, not only for Fermi sources but also non-Fermi sources, particularly, R for the former ones is averagely higher than later ones. We also adopt the same relation on core-dominance parameters and {\gamma}-ray photon indices for Fermi sources by taking the same assumption with two components model on {\gamma}-ray band, and obtain the fitting results indicating that the {\gamma}-ray emissions of Fermi blazars are mainly from the core component, which is perhaps associated with the beaming effect. Therefore, Fermi blazars are beamed.

Active galactic nuclei (AGNs) can be divided into two major classes, namely radio loud and radio quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff & Riley type I/II (FRI/II) radio galaxies. Following our previous work (Fan et al. 2011), we present a sample of 2400 sources with measured radio flux densities of the core and extended components. The sample contains 250 BL Lacs, 520 quasars, 175 Seyferts, 1178 galaxies, 153 Fanaroff & Riley type I and type II galaxies and 104 unidentified sources. We then calculate the radio core-dominance parameters and spectral indices and study their relationship. Our analysis shows that the core-dominance parameters and spectral indices are quite different for different types of sources. We also confirm that the correlation between core-dominance parameter and spectral index exists for a large sample presented in this work.

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.

Blazars are an extreme subclass of active galactic nuclei. Their rapid variability, luminous brightness, superluminal motion, and high and variable polarization are probably due to a beaming effect. However, this beaming factor (or Doppler factor) is very difficult to measure. Currently, a good way to estimate it is to use the time scale of their radio flares. In this $Letter$, we use multiwavelength data and Doppler factors reported in the literatures for a sample of 86 flaring blazars detected by Fermi to compute their intrinsic multiwavelength data and intrinsic spectral energy distributions, and investigate the correlations among observed and intrinsic data. Quite interestingly, intrinsic data show a positive correlation between luminosity and peak frequency, in contrast with the behavior of observed data, and a tighter correlation between $\gamma$-ray luminosity and the lower energy ones. For flaring blazars detected by Fermi, we conclude that (1) Observed emissions are strongly beamed; (2) The anti-correlation between luminosity and peak frequency from the observed data is an apparent result, the correlation between intrinsic data being positive; and (3) Intrinsic $\gamma$-ray luminosity is strongly correlated with other intrinsic luminosities.

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

Magnetohydronamic turbulence is believed to play a crucial role in heating the laboratorial, space, and astrophysical plasmas. However, the precise connection between the turbulent fluctuations and the particle kinetics has not yet been established. Here we present clear evidence of plasma turbulence heating based on diagnosed wave features and proton velocity distributions from solar wind measurements by the Wind spacecraft. For the first time, we can report the simultaneous observation of counter-propagating magnetohydrodynamic waves in the solar wind turbulence. Different from the traditional paradigm with counter-propagating Alfv\'en waves, anti-sunward Alfv\'en waves (AWs) are encountered by sunward slow magnetosonic waves (SMWs) in this new type of solar wind compressible turbulence. The counter-propagating AWs and SWs correspond respectively to the dominant and sub-dominant populations of the imbalanced Els\"asser variables. Nonlinear interactions between the AWs and SMWs are inferred from the non-orthogonality between the possible oscillation direction of one wave and the possible propagation direction of the other. The associated protons are revealed to exhibit bi-directional asymmetric beams in their velocity distributions: sunward beams appearing in short and narrow patterns and anti-sunward broad extended tails. It is suggested that multiple types of wave-particle interactions, i.e., cyclotron and Landau resonances with AWs and SMWs at kinetic scales, are taking place to jointly heat the protons perpendicularly and parallel.

We present observations of the power spectral anisotropy in wave-vector space of solar wind turbulence, and study how it evolves in interplanetary space with increasing heliocentric distance. For this purpose we use magnetic field measurements made by the Helios-2 spacecraft at three positions between 0.29 and 0.9 AU. To derive the power spectral density (PSD) in $(k_\parallel, k_\bot)$-space based on single-satellite measurements is a challenging task not yet accomplished previously. Here we derive the spectrum $\rm{PSD}_{\rm{2D}}$($\rm{k}_\parallel$, $\rm{k}_\bot$) from the spatial correlation function $\rm{CF}_{\rm{2D}}(r_\parallel, r_\bot)$ by a transformation according to the projection-slice theorem. We find the so constructed PSDs to be distributed in k-space mainly along a ridge that is more inclined toward the $\rm{k}_\bot$ than $\rm{k}_\parallel$ axis, a new result which probably indicates preferential cascading of turbulent energy along the $\rm{k}_\bot$ direction. Furthermore, this ridge of the distribution is found to gradually get closer to the $\rm{k}_\bot$ axis, as the outer scale length of the turbulence becomes larger while the solar wind flows further away from the Sun. In the vicinity of the $\rm{k}_\parallel$ axis, there appears a minor spectral component that probably corresponds to quasi-parallel Alfv\'enic fluctuations. Their relative contribution to the total spectral density tends to decrease with radial distance. These findings suggest that solar wind turbulence undergoes an anisotropic cascade transporting most of its magnetic energy towards larger $\rm{k}_\bot$, and that the anisotropy in the inertial range is radially developing further at scales that are relatively far from the ever increasing outer scale.

We present results from multi-epoch neutral hydrogen (HI) absorption observations of six bright pulsars with the Arecibo telescope. Moving through the interstellar medium (ISM) with transverse velocities of 10--150 AU/yr, these pulsars have swept across 1--200 AU over the course of our experiment, allowing us to probe the existence and properties of the tiny scale atomic structure (TSAS) in the cold neutral medium (CNM). While most of the observed pulsars show no significant change in their HI absorption spectra, we have identified at least two clear TSAS-induced opacity variations in the direction of B1929+10. These observations require strong spatial inhomogeneities in either the TSAS clouds' physical properties themselves or else in the clouds' galactic distribution. While TSAS is occasionally detected on spatial scales down to 10 AU, it is too rare to be characterized by a spectrum of turbulent CNM fluctuations on scales of 10-1000 AU, as previously suggested by some work. In the direction of B1929+10, an apparent correlation between TSAS and interstellar clouds inside the warm Local Bubble (LB) indicates that TSAS may be tracing the fragmentation of the LB wall via hydrodynamic instabilities. While similar fragmentation events occur frequently throughout the ISM, the warm medium surrounding these cold cloudlets induces a natural selection effect wherein small TSAS clouds evaporate quickly and are rare, while large clouds survive longer and become a general property of the ISM.