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

We entertain the exotic possibility that dark matter (DM) decays or annihilations taking place in our galaxy may produce a flux of relativistic very weakly-coupled bosons, axions or dark photons. We show that there exist several upper bounds for this flux on Earth assuming generic minimal requirements for DM, such as a lifetime longer than the age of the Universe or an annihilation rate that leaves unaffected the background evolution during matter domination. These bounds do not depend on the identity or the couplings of the bosons. We then show that this new flux cannot be large enough to explain the recent XENON1T excess, while assuming that the bosons' couplings to the Standard Model are consistent with all current experimental and observational constraints. We also discuss a possible caveat to these bounds and a route to explain the excess.

The advent of the Gaia era has led to potentially revolutionary understanding of dark matter (DM) dynamics in our galaxy, which has important consequences for direct detection (DD) experiments. In this paper, we study how the empirical DM velocity distribution inferred from Gaia-Sausage, a dominant substructure in the solar neighborhood, affects the interpretation of DD data. We survey different classes of operators in the non-relativistic effective field theory that could arise from several relativistic benchmark models and emphasize that the Gaia velocity distribution could modify both the total number of events as well as the shape of the differential recoil spectra, the two primary observables in DD experiments. Employing the euclideanized signal method, we investigate the effects of the Gaia distribution on reconstructing DM model parameters and identifying the correct DM model given a positive signal at future DD experiments. We find that for light DM with mass ~10 GeV, the Gaia distribution poses an additional challenge for characterizing DM interactions with ordinary matter, which may be addressed by combining complementary DD experiments with different targets and lowering the detection threshold. Meanwhile, for heavy DM with mass above ~30 GeV, depending on the type of DM model, there could be a (moderate) improvement in the sensitivity at future DD experiments.

Light fields with spatially varying backgrounds can modulate cosmic preheating, and imprint the nonlinear effects of preheating dynamics at tiny scales on large scale fluctuations. This provides us a unique probe into the preheating era which we dub the "cosmic microscope.'' We identify a distinctive effect of preheating on scalar perturbations that turns the Gaussian primordial fluctuations of a light scalar field into square waves, like a diode. The effect manifests itself as local non-Gaussianity. We present a model, "modulated partial preheating," where this nonlinear effect is consistent with current observations and can be reached by near future cosmic probes.

A scalar Higgs field can be repeatedly switched on and off when it couples to a classically oscillating scalar modulus field. The modulus flips the Higgs mass term between stable and tachyonic values. We study a cosmological scenario in which such repeated phase transitions occur during inflation. An irrelevant operator coupling the Higgs field to the inflaton can then imprint the pattern of phase transitions in the correlation functions of the inflaton. Using both numerical and analytic studies, we show that the inflaton 2-point function carries characteristic imprints of the modulus oscillation and its effect on the Higgs boson. We briefly remark on the potential observability of such patterns and how they might be distinguished from other dynamics in the early universe.

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.

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

In our previous work in Xiao et al. (2019), we suggested that 6 superluminal sources could be gamma-ray candidates, and in fact 5 of them have been confirmed in the fourth Fermi-LAT source catalogue (4FGL). In this work, based on the 4FGL, we report a sample of 229 Fermi detected superluminal sources (FDSs) including 40 new FDSs and 62 non-Fermi detected superluminal sources (non-FDSs). Thus, we believe that all superluminal sources should have $\gamma$-ray emissions, and superluminal motion could also be a clue to detect $\gamma$-ray emission from active galactic nuclei (AGN). We present a new approach of Doppler factor estimate through the study of the $\gamma$-ray luminosity ($L_{\rm \gamma}$) and of the viewing angle ($\phi$).

We continuously monitored the blazar S5 0716+714 in the optical $g$, $r$ and $i$ bands from Nov. 10, 2017 to Jun. 06, 2019. The total number of observations is 201 nights including 26973 data points. This is a very large quasi-simultaneous multicolor sample for the blazar. The average time spans and time resolutions are 3.4 hours and 2.9 minutes per night, respectively. During the period of observations, the target source in the $r$ band brightens from $14^{\rm m}.16$ to $12^{\rm m}.29$ together with five prominent sub-flares, and then first becomes fainter to $14^{\rm m}.76$ and again brightens to $12^{\rm m}.94$ with seven prominent sub-flares. For the long-term variations, we find a strong flatter when brighter (FWB) trend at a low flux state and then a weak FWB trend at a higher flux state. A weak FWB trend at a low flux state and then a strong FWB trend at a higher flux state are also reported. Most of sub-flares show the strong FWB trends, except for two flares with a weak FWB trend. The particle acceleration and cooling mechanisms together with the superposition of different FWB-slopes from sub-flares are likely to explain the optical color behaviours. A scenario of bent jet is discussed.

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

The recently published fourth Fermi Large Area Telescope source catalog (4FGL) reports 5065 gamma-ray sources in terms of direct observational gamma-ray properties. Among the sources, the largest population is the Active Galactic Nuclei (AGN), which consists of 3137 blazars, 42 radio galaxies, and 28 other AGNs. The blazar sample comprises 694 flat-spectrum radio quasars (FSRQs), 1131 BL Lac-type objects (BL Lacs), and 1312 blazar candidates of an unknown type (BCUs). The classification of blazars is difficult using optical spectroscopy given the limited knowledge with respect to their intrinsic properties, and the limited availability of astronomical observations. To overcome these challenges, machine learning algorithms are being investigated as alternative approaches. Using the 4FGL catalog, a sample of 3137 Fermi blazars with 23 parameters is systematically selected. Three established supervised machine learning algorithms (random forests (RFs), support vector machines (SVMs), artificial neural networks (ANNs)) are employed to general predictive models to classify the BCUs. We analyze the results for all of the different combinations of parameters. Interestingly, a previously reported trend the use of more parameters leading to higher accuracy is not found. Considering the least number of parameters used, combinations of eight, 12 or 10 parameters in the SVM, ANN, or RF generated models achieve the highest accuracy (Accuracy $\simeq$ 91.8\%, or $\simeq$ 92.9\%). Using the combined classification results from the optimal combinations of parameters, 724 BL Lac type candidates and 332 FSRQ type candidates are predicted; however, 256 remain without a clear prediction.

We present a detailed study of 11 narrow-line Seyfert 1 galaxies (NLS1s) from the Six-degree Field Galaxy Survey (6dFGS) that both have optical and X-ray spectroscopic observations. There are five complex NLS1s (C-NLS1s) and six simple NLS1s (S-NLS1s). We propose a possible correlation between [O III] line asymmetry and X-ray complexity. The outflow or wind from the inner accretion disk is commonly present in NLS1s and mostly directed along the system axis. In C-NLS1s only weak wind effects are measured, the X-ray spectral complexity might be caused by the presence of ionized material in the wind. On the contrary, the wind in S-NLS1s is fast, the ionized material could be swept by such a strong wind, thus the complex feature is missing which results in a simple X-ray spectrum. Furthermore, this outflow scenario seems to be an inclination effect. Since the speed of the wind is higher in a small inclination while lower in a large inclination, S-NLS1s might be sources viewed at small angles while C-NLS1s might be sources viewed at large angles.

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.

We use stellar kinematics from the latest Gaia data release (DR2) to measure the local dark matter density $\rho_{\rm DM}$ in a heliocentric cylinder of radius $R= 150 \ {\rm pc}$ and half-height $z= 200 \ {\rm pc}$. We also explore the prospect of using our analysis to estimate the DM density in local substructure by setting constraints on the surface density and scale height of a thin dark disk aligned with the baryonic disk and formed due to dark matter self-interaction. Performing the statistical analysis within a Bayesian framework for three types of tracers, we obtain ${\rho_{\rm DM}= 0.016 \pm 0.010}$ M$_\odot$/pc$^3$ for A stars; early G stars give a similar result, while F stars yield a significantly higher value. For a thin dark disk, A stars set the strongest constraint: excluding surface densities (5-12) M$_\odot$/pc$^2$ for scale heights below 100 pc with 95% confidence. Comparing our results with those derived using Tycho-Gaia Astrometric Solution (TGAS) data, we find that the uncertainty in our measurements of the local DM content is dominated by systematic errors that arise from assumptions of our kinematic analysis in the low $z$ region. Furthermore, there will only be a marginal reduction in these uncertainties with more data in the Gaia era. We comment on the robustness of our method and discuss potential improvements for future work.

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 study properties of boson stars consisting of ultra-light scalar dark matter with repulsive self-interactions. We investigate the origin of the maximum mass of spherically symmetric stable stars which emerges only when solving the full equations of motion in curved space-time, but not when solving the approximated Schr\"odinger-Newton equations. When the repulsion is weak, the backreaction of the curvature on the scalars acts as an additional source of attraction and can overcome the repulsion, resulting in a maximum star mass and compactness. We also point out that the potential in a UV completed particle physics model of light scalar dark matter is generally more complicated than the widely used $\phi^4$ interaction. Additional interactions beyond $\phi^4$ in the potential can dramatically change the properties of boson stars as well as modify the prospect of LIGO gravitational wave detection for binary mergers of boson stars.

In the third catalog of active galactic nuclei detected by the Fermi-LAT (3LAC) Clean Sample, there are 402 blazars candidates of uncertain type (BCU). Due to the limitations of astronomical observation or intrinsic properties, it is difficult to classify blazars using optical spectroscopy. The potential classification of BCUs using machine learning algorithms is essential. Based on the 3LAC Clean Sample, we collect 1420 Fermi blazars with 8 parameters of {\gamma}-ray photon spectral index, radio flux, flux density, curve significance, the integral photon flux in 100 to 300 MeV, 0.3 to 1 GeV, 10 to 100 GeV and variability index. Here, we apply 4 different supervised machine learning (SML) algorithms (\emph{Decision trees, Random forests, support vector machines and Mclust Gaussian finite mixture models}) to evaluate the classification of BCUs based on the direct observational properties. All the 4 methods can perform exceedingly well with a more accuracy and can effective forecast the classification of Fermi BCUs. The evaluating results show the results of these methods (SML) are valid and robust, where, about 1/4 sources are FSRQs and 3/4 are BL Lacs in 400 BCUs, which are consistent with some other recent results. Although a number of factors influence the accuracy of SML, the results are stable at a fixed ratio 1:3 between FSRQs and BL Lacs, which suggests that the SML can provides an effective method to evaluate the potential classification of BCUs. Among the 4 methods, Mclust Gaussian Mixture Modelling has the highest accuracy for our training sample (4/5, seed=123).

The Higgs potential appears to be fine-tuned, hence very sensitive to values of other scalar fields that couple to the Higgs. We show that this feature can lead to a new epoch in the early universe featuring violent dynamics coupling the Higgs to a scalar modulus. The oscillating modulus drives tachyonic Higgs particle production. We find a simple parametric understanding of when this process can lead to rapid modulus fragmentation, resulting in gravitational wave production. A nontrivial equation-of-state arising from the nonlinear dynamics also affects the time elapsed from inflation to the CMB, influencing fits of inflationary models. Supersymmetric theories automatically contain useful ingredients for this picture.

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.

Utilizing all the 16-year $RXTE$ observations, we analyze the X-ray spectra of 32 TeV blazars, and perform a systematic investigation of X-ray spectral variability for the 5 brightest sources during their major flares that lasted several days. We obtain photon spectral index ($\alpha$), flux and synchrotron radiation peak energy ($E_{\rm{p}}$) from empirical spectral fitting, and electron spectral index ($p$) from theoretical synchrotron radiation modeling. We find that both $\alpha$ and $p$ generally display a harder-when-brighter trend, confirming the results of many previous works. Furthermore, we confirm and strengthen the result that $p$ must vary in order to explain the observed X-ray spectral variability during flares, which would have useful implications for interpreting the associated higher-energy spectral variability. We see apparent electron spectral hysteresis in many but not all $p$-flux plots that takes a form of "loop" or oblique"8". We obtain a tight $p$-hardness ratio (HR) relation and a tighter $p$-$\alpha$ relation using spectra of flaring periods, both of which are also applicable to stacked data of quiescent periods. We demonstrate that these two empirical relations can be used efficiently to estimate $p$ from HR or $\alpha$ that is readily achieved. Finally, we find that, when considering TeV blazars as a whole, $\alpha$ and X-ray luminosity are positively correlated, $E_{\rm{p}}$ is negatively correlated with $p$ and $\alpha$, and $E_{\rm{p}}$ is positively correlated with HR; all these correlations are in line with the blazar sequence. However, after correcting for the Doppler boosting effect, $\alpha$ and intrinsic X-ray luminosity follow an anti-correlation.

We present a new accurate sample of narrow-line Seyfert 1 galaxies (NLS1s) in the southern hemisphere from the Six-degree Field Galaxy Survey (6dFGS). Based on the optical spectral features, 167 sources were classified as NLS1s. We derived flux-calibrated spectra in the sample for the first time. Strong luminosity correlations between the continuum and the emission lines were found. We estimated their central black hole masses and Eddington ratios, which are lying in a typical range of NLS1s. In the sample, 23 NLS1s were detected at radio frequencies and 12 of them are radio-loud. We publish the X-ray data analysis of a campaign of observations carried out by the Swift X-ray Telescope.

We present a new accurate catalog of narrow-line Seyfert 1 galaxies (NLS1s) in the southern hemisphere from the Six-degree Field Galaxy Survey (6dFGS) final data release, which is currently the most extensive spectroscopic survey available in the southern sky whose database has not yet been systematically explored. We classified 167 sources as NLS1s based on their optical spectral properties. We derived flux-calibrated spectra for the first time that the 6dFGS indeed does not provide. By analyzing these spectra, we obtained strong correlations between the monochromatic luminosity at 5100 A and the luminosity of Hbeta and [O III]5007 lines. The estimated central black hole mass and Eddington ratio have an average value of 8.6 x 10^6 M_Sun and 0.96 L_Edd respectively, which is a typical value for NLS1s. In the sample, 23 (13.8%) NLS1s were detected at radio frequencies, and 12 (7.0%) of them are radio-loud. Our results confirmed that radio-loud sources tend to have a higher redshift, more massive black hole, and higher radio and optical luminosity than radio-quiet sources.

Thanks to the execution of extensive spectroscopic surveys that have covered large fractions of the sky down to magnitudes as faint as $i \approx 19$, it has been possible to identify several narrow-line Seyfert 1 galaxies (NLS1s) and to investigate their properties over a large range of the electro-magnetic spectrum. The interpretation of their nature, however, is still hampered by the statistical uncertainties related to the way in which NLS1 candidates are selected. In this contribution, we present a study to detect and to model emission lines in optical spectra extracted from the Sloan Digital Sky Survey (SDSS), adopting the most proper strategy to identify the source of line excitation and to produce a detailed model with measurements of several emission-line parameters. We demonstrate the application of this technique to explore fundamental questions, such as the presence of gas and dust around the core of AGNs and the spectral energy distribution of their ionizing radiation. We compare the spectral properties of NLS1s with those of other type 1 active galaxies and we summarize the potential of this approach to identify NLS1s in present day and future spectroscopic surveys. We finally consider the implications of multi-frequency data analysis in the debate concerning the intrinsic nature of NLS1s.

Many cosmological models rely on large couplings of axions to gauge fields. Examples include theories of magnetogenesis, inflation on a steep potential, chiral gravitational waves, and chromonatural inflation. Such theories require a mismatch between the axion field range and the mass scale appearing in the $a F \widetilde{F}$ coupling. This mismatch suggests an underlying monodromy, with the axion winding around its fundamental period a large number of times. We investigate the extent to which this integer can be explained as a product of smaller integers in a UV completion: in the parlance of our times, can the theory be "clockworked"? We argue that a clockwork construction producing a potential $\mu^4 \cos(\frac{a}{j F_a})$ for an axion of fundamental period $F_a$ will obey the constraint $\mu < F_a$. For some applications, including chromonatural inflation with sub-Planckian field range, this constraint obstructs a clockwork UV completion. Alternative routes to a large coupling include fields of large charge (an approach limited by strong coupling) or kinetic mixing (requiring a lighter axion). Our results suggest that completions of axion cosmologies that explain the large parameter in the theory potentially alter the phenomenological predictions of the model.

In this work, we collected a sample of 69 TeV blazars from TeVCat, obtained their multi-wavelength observations, and fitted their spectral energy distributions by using the second degree polynomial function. The structure parameters of the synchrotron bumps for 68 blazars and those of the inverse-Compton bumps for 56 blzars are obtained. Then, we adopted statistical analysis to the parameters (curvature, peak frequency, peak luminosity, bolometric luminosity, and X/$\gamma$-ray spectral indexes). From our analysis and discussions, we can get following conclusions: 1. There is a clear positive correlation between the synchrotron peak frequency, $\log \nu_p^s$, and the inverse-Compton peak frequency $\log \nu_p^{IC}$, and that between the synchrotron peak luminosity, $\log \nu_p^s L_p^s$, and the inverse-Compton peak luminosity, $\log \nu_p^{IC}L_p^{IC}$. 2. The correlation between the peak frequency and the curvature of synchrotron bump is clearly different from that of the inverse-Compton bump, which further indicates that there are different emission mechanisms between them. 3. There is a correlation between $\log \nu_p^{IC}$ and $\gamma$-ray spectral index, $\alpha_\gamma$, for the TeV blazars: $\log \nu^{IC}_p = -(4.59 \pm 0.30) \alpha_\gamma + (32.67 \pm 0.59)$, which is consistent with previous work of Abdo et al.(2010). 4. An "L-shape" relation is found between $\log \nu^s_p$ and $\alpha_X$ for both TeV blazars and Fermi blazars. A significant correlation between $\log \nu^s_p$ and X-ray photon index ($\alpha_X$) is found for the TeV blazars with high synchrotron peak frequency: $\log \nu^{s}_{p} = -(3.20 \pm 0.34) \alpha_X + (24.33 \pm 0.79)$, while the correlation is positive for low synchrotron peaked TeV sources. 5. In the $\alpha_X-\alpha_\gamma$ diagram, there is also an "L-shape", the anti-correlation is consistent with the available results in the literature, (...)

The QCD axion's coupling to photons is often assumed to lie in a narrow band as a function of the axion mass. We demonstrate that several simple mechanisms, in addition to the photophilic clockwork axion already in the literature, can significantly extend the allowed range of couplings. Some mechanisms we present generalize the KNP alignment scenario, widely studied as a model of inflation, to the phenomenology of a QCD axion. In particular we present KSVZ-like realizations of two-axion KNP alignment and of the clockwork mechanism. Such a "confinement tower" realization of clockwork may prove useful in a variety of model-building contexts. We also show that kinetic mixing of the QCD axion with a lighter axion-like particle can dramatically alter the QCD axion's coupling to photons, differing from the other models we present by allowing non-quantized couplings. The simple models that we present fully cover the range of axion-photon couplings that could be probed by experiments. They motivate growing axion detection efforts over a wide space of masses and couplings.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

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.

The possibility of a subdominant component of dark matter dissipating energy could lead to dramatic new phenomenology such as the formation of a dark disk. One rigorous way to assess this possibility and settle the debate on its feasibility is to include the dissipative dark component in a numerical hydrodynamical simulation. A necessary input to such a simulation is a prescription including energy dissipation rates of different processes and rates of processes that change the number densities of dark ions and atoms. In this article, we study the simplest dissipative dark sector which consists of a dark electron and proton, both charged under a dark gauged U(1). We present approximate analytic formulas for energy loss rates due to Compton scattering, bremsstrahlung, recombination, collisional ionization and collisional excitation as well as the rates of number density change. We also include the heating rate due to photoionization. The work serves as the first step to realize a numerical simulation including a dissipative dark sector, which hopefully can shed more light on the formation and properties of a dark disk originating from dark matter self-interactions.

Beaming effect is important for the observational properties of blazars. In this work, we collect 91 $Fermi$ blazars with available radio Doppler factors. $\gamma$-ray Doppler factors are estimated and compared with radio Doppler factors for some sources. The intrinsic (de-beamed) $\gamma$-ray flux density ($f^{\rm in}_{\gamma}$), intrinsic $\gamma$-ray luminosity ($L^{\rm in}_{\gamma}$), and intrinsic synchrotron peak frequency ($\nu_{\rm p}^{\rm in}$) are calculated. Then we study the correlations between $f^{\rm in}_{\gamma}$ and redshift and find that they follow the theoretical relation: $\log f = -2.0 \log z + {\rm const}$. When the subclasses are considered, we find that stationary jets are perhaps dominant in low synchrotron peaked blazars. 63 $Fermi$ blazars with both available short variability time scales ($\Delta T$) and Doppler factors are also collected. We find that the intrinsic relationship between $L ^{\rm in}_{\gamma}$ and $\Delta T^{\rm in}$ obeys the Elliot \& Shapiro and the Abramowicz \& Nobili relations. Strong positive correlation between $f_{\gamma}^{\rm in}$ and $\nu_{\rm p}^{\rm in}$ is found, suggesting that synchrotron emissions are highly correlated with $\gamma$-ray emissions.

In a two-component jet model, the emissions are the sum of the core and extended emissions: $S^{\rm ob}=S_{\rm core}^{\rm ob}+S_{\rm ext}^{\rm ob}$, with the core emissions, $S_{\rm core}^{\rm ob}= f S_{\rm ext}^{\rm ob}\delta^{q}$, being a function of the Doppler factor, $\delta$, the extended emission, $S_{\rm ext}^{\rm ob}$, jet type dependent factor, $q$, and the ratio of the core to the extended emissions in the comoving frame, $f$. The $f$ is an unobservable but important parameter. Following our previous work, we collect 65 blazars with available Doppler factor, $\delta$, superluminal velocity, $\beta_{app}$, and core-dominance parameter, $R$, calculate the ratio, $f$, and peform statistical analyses. We find that the ratio, $f$, in BL Lacs is on average larger than that in FSRQs. We suggest that the difference of the ratio $f$ between FSRQs and BL Lacs is one of the possible reasons that cause the difference of other observed properties between them. We also find some significant correlations between $\log f$ and other parameters, including intrinsic (de-beamed) peak frequency, $\log \nu _{\rm p}^{\rm in}$, intrinsic polarization, $\log P^{\rm in}$, and core-dominance parameter, $\log R$, for the whole sample. In addition, we show that the ratio, $f$, can be estimated by $R$.

We have monitored the quasar 3C 273 in optical $V$, $R$ and $I$ bands from 2005 to 2016. Intraday variability (IDV) is detected on seven nights. The variability amplitudes for most of nights are less than 10\% and four nights more than 20\%. When considering the nights with time spans $>4$ hours, the value of duty cycle (DC) is 14.17 per cent. Over the twelve years, the overall magnitude and color index variabilities are $\bigtriangleup I=0^{\rm m}.67$, $\bigtriangleup R=0^{\rm m}.72$, $\bigtriangleup V=0^{\rm m}.68$, and $\bigtriangleup (V-R)=0^{\rm m}.25$ respectively. The largest clear IDV has an amplitude of 42% over just 5.8 minutes and the weakest detected IDV is 5.4% over 175 minutes. The BWB (bluer when brighter) chromatic trend is dominant for 3C 273 and appears at different flux levels on intraday timescales. The BWB trend exists for short-term timescales and intermediate-term timescales but different timescales have different correlations. There is no BWB trend for our whole time-series data sets. A significant anti-correlation between BWB trend and length of timescales is found. Combining with $V$-band data from previous works, we find a possible quasi-periodicity of $P=3918\pm1112$ days. The possible explanations for the observed variability, BWB chromatic trend and periodicity are discussed.

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.

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.

We present spectral line mapping observations toward four massive star-forming regions (Cepheus A, DR21S, S76E and G34.26+0.15), with the IRAM 30 meter telescope at 2 mm and 3 mm bands. Totally 396 spectral lines from 51 molecules, one helium recombination line, ten hydrogen recombination lines, and 16 unidentified lines were detected in these four sources. An emission line of nitrosyl cyanide (ONCN, 14$_{0,14}$-13$_{0,13}$) was detected in G34.26+0.15, as first detection in massive star-forming regions. We found that the $c$-C$_{3}$H$_{2}$ and NH$_{2}$D show enhancement in shocked regions as suggested by evidences of SiO and/or SO emission. Column density and rotational temperature of CH$_{3}$CN were estimated with the rotational diagram method for all four sources. Isotope abundance ratios of $^{12}$C/$^{13}$C were derived using HC$_{3}$N and its $^{13}$C isotopologue, which were around 40 in all four massive star-forming regions and slightly lower than the local interstellar value ($\sim$65). $^{14}$N/$^{15}$N and $^{16}$O/$^{18}$O abundance ratios in these sources were also derived using double isotopic method, which were slightly lower than that in local interstellar medium. Except for Cep A, $^{33}$S/$^{34}$S ratio in the other three targets were derived, which were similar to that in the local interstellar medium. The column density ratios of N(DCN)/N(HCN) and N(DCO$^{+}$)/N(HCO$^{+}$) in these sources were more than two orders of magnitude higher than the elemental [D]/[H] ratio, which is 1.5$\times$10$ ^{-5}$. Our results show the later stage sources, G34.26+0.15 in particular, present more molecular species than earlier stage ones. Evidence of shock activity is seen in all stages studied.

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability.

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

BL Lacertae objects (BL Lacs) is the dominant population of TeV emitting blazars. In this work, we investigate whether there is any special observational properties for TeV sources. To do so, we will compare the observational properties of TeV detected BL Lacs (TeV BLs) and non-TeV detected BL Lac objects (non-TeV BLs). From the 3rd $Fermi$/LAT catalog (3FGL), we can get 662 BL Lacs, out of which, 47 are TeV BLs and 615 are non-TeV BLs. Their multi-wavelength flux densities ($F_{\rm R}$, $F_{\rm O}$, $F_{\rm X}$, $F_{\gamma}$), photon spectral indexes ($\alpha_{\rm X}^{\rm ph}$, $\alpha_{\gamma}^{\rm ph}$), and effective spectral indexes ($\alpha_{\rm RO}$ and $\alpha_{\rm OX}$) are compiled from the available literatures. Then the luminosities ($\log\,{\nu}L_{\rm R}$, $\log\,{\nu}L_{\rm O}$, $\log\,{\nu}L_{\rm X}$, $\log\,{\nu}L_{\gamma}$) are calculated. From comparisons, we found that TeV BLs are different from low-synchrotron-peaked BLs (LSP) and intermediate-synchrotron-peaked BLs (ISP), but TeV BLs show similar properties as high-synchrotron-peaked BLs (HSP). Therefore, we concentrated on comparison between TeV HSP BLs and non-TeV HSP BLs. Analysis results suggest that TeV HSP BLs and non-TeV HSP BLs show some differences in their $\alpha_{\rm RO}$ and $\alpha_{\gamma}^{\rm ph}$, while their other properties are quite similar.

Ultralight scalar dark matter with mass at or below the eV scale and pressure from repulsive self-interaction could form a Bose-Einstein condensate in the early Universe and maybe in galaxies as well. It has been suggested to be a possible solution to the cusp/core problem or even to explain MOND phenomenology. In this paper, I initiate a study of possible self-interactions of ultralight scalar dark matter from the particle physics point of view. To protect its mass, the scalar dark matter is identified as a pseudo Nambu-Goldstone boson (pNGB). Quite a few pNGB models with different potentials such as the QCD axion and the dilaton lead to attractive self-interactions. Yet if an axion is a remnant of a 5D gauged U(1) symmetry, its self-interactions could be repulsive provided the masses and charges of the 5D matter contributing to its potential satisfy certain constraints. Collective symmetry breaking could also lead to a repulsive self-interaction yet with too large a strength that is ruled out by Bullet Cluster constraints. I also discuss cosmological and astrophysical constraints on ultralight repulsive dark matter in terms of a parametrization motivated by particle physics considerations.

We analyzed the multi-band optical behaviour of the BL Lacertae object, S5 0716+714, during its outburst state from 2014 November - 2015 March. We took data on 23 nights at three observatories, one in India and two in Bulgaria, making quasi-simultaneous observations in B, V, R, and I bands. We measured multi-band optical fluxes, colour and spectral variations for this blazar on intraday and short timescales. The source was in a flaring state during the period analyzed and displayed intense variability in all wavelengths. R band magnitude of 11.6 was attained by the target on 18 Jan 2015, which is the brightest value ever recorded for S5 0716+714. The discrete correlation function method yielded good correlation between the bands with no measurable time lags, implying that radiation in these bands originate from the same region and by the same mechanism. We also used the structure function technique to look for characteristic timescales in the light curves. During the times of rapid variability, no evidence for the source to display spectral changes with magnitude was found on either of the timescales. The amplitude of variations tends to increase with increasing frequency with a maximum of $\sim$ 22% seen during flaring states in B band. A mild trend of larger variability amplitude as the source brightens was also found. We found the duty cycle of our source during the analyzed period to be $\sim$ 90%. We also investigated the optical spectral energy distribution of S5 0716+714 using B, V, R, and I data points for 21 nights. We briefly discuss physical mechanisms most likely responsible for its flux and spectral variations.

High-scale supersymmetry (SUSY) with a split spectrum has become increasingly interesting given the current experimental results. A SUSY scale above the weak scale could be naturally associated with a heavy unstable gravitino, whose decays populate the dark matter (DM) particles. In the mini-split scenario with gravitino at about the PeV scale and the lightest TeV scale neutralino being (a component of) DM, the requirement that the DM relic abundance resulting from gravitino decays does not overclose the Universe and satisfies the indirect detection constraints demand the reheating temperature to be below 10^9 - 10^{10} GeV. On the other hand, the BICEP2 result prefers a heavy inflaton with mass at around 10^{13} GeV and a reheating temperature at or above 10^9 GeV with some general assumptions. The mild tension could be alleviated if SUSY scale is even higher with the gravitino mass above the PeV scale. Intriguingly, in no-scale supergravity, gravitinos could be very heavy at about 10^{13} GeV, the inflaton mass scale, while gauginos could still be light at the TeV scale.

We study the constraints from direct detection and solar capture on dark matter scenarios with a subdominant dissipative component. This dissipative dark matter component in general has both a symmetric and asymmetric relic abundance. Dissipative dynamics allow this subdominant dark matter component to cool, resulting in its partial or total collapse into a smaller volume inside the halo (e.g., a dark disk) as well as a reduced thermal velocity dispersion compared to that of normal cold dark matter. We first show that these features considerably relax the limits from direct detection experiments on the couplings between standard model (SM) particles and dissipative dark matter. On the other hand, indirect detection of the annihilation of the symmetric dissipative dark matter component inside the Sun sets stringent and robust constraints on the properties of the dissipative dark matter. In particular, IceCube observations force dissipative dark matter particles with mass above 50 GeV to either have a small coupling to the SM or a low local density in the solar system, or to have a nearly asymmetric relic abundance. Possible helioseismology signals associated with purely asymmetric dissipative dark matter are discussed, with no present constraints.

We examine the evolution of cosmological perturbations in a non-thermal post-inflationary history with a late-time matter domination period prior to BBN. Such a cosmology could arise naturally in the well-motivated moduli scenario in the context of supersymmetry (SUSY) -- in particular in models of Split-SUSY. Sub-horizon dark matter perturbations grow linearly during the matter dominated phase before reheating and can lead to an enhancement in the growth of substructure on small scales, even in the presence of dark matter annihilations. This suggests that a new scale (the horizon size at reheating) could be important for determining the primordial matter power spectrum. However, we find that in many non-thermal models free-streaming effects or kinetic decoupling after reheating can completely erase the enhancement leading to small-scale structures. In particular, in the moduli scenario with wino or higgsino dark matter we find that the dark matter particles produced from moduli decays would thermalize with radiation and kinetically decouple below the reheating temperature. Thus, the growth of dark matter perturbations is not sustained, and the predictions for the matter power spectrum are similar to a standard thermal history. We comment on possible exceptions, but these appear difficult to realize within standard moduli scenarios. We conclude that although enhanced structure does not provide a new probe for investigating the cosmic dark ages within these models, it does suggest that non-thermal histories offer a robust alternative to a strictly thermal post-inflationary history.

We point out that current constraints on dark matter imply only that the majority of dark matter is cold and collisionless. A subdominant fraction of dark matter could have much stronger interactions. In particular, it could interact in a manner that dissipates energy, thereby cooling into a rotationally-supported disk, much as baryons do. We call this proposed new dark matter component Double-Disk Dark Matter (DDDM). We argue that DDDM could constitute a fraction of all matter roughly as large as the fraction in baryons, and that it could be detected through its gravitational effects on the motion of stars in galaxies, for example. Furthermore, if DDDM can annihilate to gamma rays, it would give rise to an indirect detection signal distributed across the sky that differs dramatically from that predicted for ordinary dark matter. DDDM and more general partially interacting dark matter scenarios provide a large unexplored space of testable new physics ideas.

Recently evidence for gamma ray lines at energies of approximately 111 and 128 GeV has been found in Fermi-LAT data from the center of the galaxy and from unassociated point sources. Many explanations in terms of dark matter particle pairs annihilating to gamma gamma and gamma Z have been suggested, but these typically require very large couplings or mysterious coincidences in the masses of several new particles to fit the signal strength. We propose a simple novel explanation in which dark matter is part of a multiplet of new states which all have mass near 260 GeV as a result of symmetry. Two dark matter particles annihilate to a pair of neutral particles in this multiplet which subsequently decay to gamma gamma and gamma Z. For example, one may have a triplet of pseudo-Nambu-Goldstone bosons, charged pions and a neutral pion, where the charged pions are stabilized by their charge under a new U(1) symmetry and the slightly lighter neutral pion state decays to gamma gamma and gamma Z. The symmetry structure of such a model explains the near degeneracy in masses needed for the resulting photons to have a line-like shape and the large observed flux. The tunable lifetime of the neutral state allows such models to go unseen at direct detection or collider experiments that can constrain most other explanations. However, nucleosynthesis constraints on the neutral pion lifetime fix a minimum necessary coupling between the new multiplet and the Standard Model. The spectrum is predicted to be not a line but a box with a width of order a few GeV, smaller than but on the order of the Fermi-LAT resolution.

Based on observational constraints on large scale structure and halo structure, dark matter is generally taken to be cold and essentially collisionless. On the other hand, given the large number of particles and forces in the visible world, a more complex dark sector could be a reasonable or even likely possibility. This hypothesis leads to testable consequences, perhaps portending the discovery of a rich hidden world neighboring our own. We consider a scenario that readily satisfies current bounds that we call Partially Interacting Dark Matter (PIDM). This scenario contains self-interacting dark matter, but it is not the dominant component. Even if PIDM contains only a fraction of the net dark matter density, comparable to the baryonic fraction, the subdominant component's interactions can lead to interesting and potentially observable consequences. Our primary focus will be the special case of Double-Disk Dark Matter (DDDM), in which self-interactions allow the dark matter to lose enough energy to lead to dynamics similar to those in the baryonic sector. We explore a simple model in which DDDM can cool efficiently and form a disk within galaxies, and we evaluate some of the possible observational signatures. The most prominent signal of such a scenario could be an enhanced indirect detection signature with a distinctive spatial distribution. Even though subdominant, the enhanced density at the center of the galaxy and possibly throughout the plane of the galaxy can lead to large boost factors, and could even explain a signature as large as the 130 GeV Fermi line. Such scenarios also predict additional dark radiation degrees of freedom that could soon be detectable and would influence the interpretation of future data, such as that from Planck and from the Gaia satellite. We consider this to be the first step toward exploring a rich array of new possibilities for dark matter dynamics.

Indirect detection constraints on gamma rays (both continuum and lines) have set strong constraints on wino dark matter. By combining results from Fermi-LAT and HESS, we show that: light nonthermal wino dark matter is strongly excluded; thermal wino dark matter is allowed only if the Milky Way dark matter distribution has a significant (>~0.4 kpc) core; and for plausible NFW and Einasto distributions the entire range of wino masses from 100 GeV up to 3 TeV can be excluded. The case of light, nonthermal wino dark matter is particularly interesting in scenarios with decaying moduli that reheat the universe to a low temperature. Typically such models have been discussed for low reheating temperatures, not far above the BBN bound of a few MeV. We show that constraints on the allowed wino relic density push such models to higher reheating temperatures and hence heavier moduli. Even for a flattened halo model consisting of an NFW profile with constant-density core inside 1 kpc and a density near the sun of 0.3 GeV/cm^3, for 150 GeV winos current data constrains the reheat temperature to be above 1.4 GeV. As a result, for models in which the wino mass is a loop factor below the gravitino mass, the data favor moduli that are more than an order of magnitude heavier than the gravitino. We discuss some of the sobering implications of this result for the status of supersymmetry. We also comment on other neutralino dark matter scenarios, in particular the case of mixed bino/higgsino dark matter. We show that in this case, direct and indirect searches are complementary to each other and could potentially cover most of the parameter space.

The \gamma-ray loud blazars (flat spectrum radio quasars--FSRQs and BL Lacertae objects-BLs) are very bright in the \gamma-ray bands, which is perhaps associated with a beaming effect. Therefore, one can expect that the \gamma-ray luminosity is correlated with the beaming factor. In this paper, we investigated the relation between the radio Doppler factors and the gamma-ray luminosities. Our analysis suggests that the \gamma-ray luminosity be strongly correlated with the factor of \delta_R for the whole sample, FSRQs, and BLs. When the effect of a common redshift is excluded, the correlation still exists for the FSRQs sub-sample suggesting that the \gamma-rays are strongly beamed. However, the partial correlation analysis does not show a correlation for the small BL Lac sample.

Gamma-ray bursts (GRBs) are divided into two classes according to their durations. We investigate if the softness of bursts plays a role in the conventional classification of the objects. We employ the BATSE (Burst and Transient Source Experiment) catalog and analyze the duration distributions of different groups of GRBs associated with distinct softness. Our analysis reveals that the conventional classification of GRBs with the duration of bursts is influenced by the softness of the objects. There exits a bimodality in the duration distribution of GRBs for each group of bursts and the time position of the dip in the bimodality histogram shifts with the softness parameter. Our findings suggest that the conventional classification scheme should be modified by separating the two well-known populations in different softness groups, which would be more reasonable than doing so with a single sample. According to the relation between the dip position and the softness parameter, we get an empirical function that can roughly set apart the short-hard and long-soft bursts: $SP = (0.100 \pm 0.028) T_{90}^{-(0.85 \pm 0.18)}$, where $SP$ is the softness parameter adopted in this paper.

To shed light onto the circumnuclear environment of 22 GHz H2O maser galaxies, we have analyzed some of their multi-wavelength properties, including the far infrared luminosity (FIR), the luminosity of the [O III]\lambda5007 emission line, the nuclear X-ray luminosity, and the equivalent width of the neutral iron Ka emission line (EW (Ka)). Our statistical analysis includes a total of 85 sources, most of them harboring an active galactic nucleus (AGN). There are strong anti-correlations between EW (Ka) and two "optical thickness parameters", i.e. the ratios of the X-ray luminosity versus the presumably more isotropically radiated [O III] and far infrared (FIR) luminosities. Based on these anti-correlations, a set of quantitative criteria, EW (Ka) > 300eV, L_{2-10keV} < 2L_[O III] and L_{FIR} > 600L_{2-10keV} can be established for Compton-thick nuclear regions. 18 H2O maser galaxies belong to this category. There are no obvious correlations between the EW (Ka), the [O III] luminosity and the isotropic H2O maser luminosity. When comparing samples of Seyfert 2s with and without detected H2O maser lines, there seem to exist differences in EW (Ka) and the fraction of Compton-thick nuclei. This should be studied further. For AGN masers alone, there is no obvious correlation between FIR and H2O maser luminosities. However, including masers associated with star forming regions, a linear correlation is revealed. Overall, the extragalactic FIR-H2O data agree with the corresponding relation for Galactic maser sources, extrapolated by several orders of magnitude to higher luminosities.

We monitored the BL Lac object S5 0716+714 in the optical band during October 2008, December 2008 and February 2009 with a best temporal resolution of about 5 minutes in the BVRI bands. Four fast flares were observed with amplitudes ranging from 0.3 to 0.75 mag. The source remained active during the whole monitoring campaign, showing microvariability in all days except for one. The overall variability amplitudes are delta B ~ 0.89 mag, delta V ~ 0.80 mag, delta R ~ 0.73 mag and delta I ~0.51 mag. Typical timescales of microvariability range from 2 to 8 hours. The overall V - R color index ranges from 0.37 to 0.59. Strong bluer- when-brighter chromatism was found on internight timescales. However, different spectral behavior was found on intranight timescales. A possible time lag of ~ 11 mins between B and I bands was found on one night. The shock-in-jet model and geometric effects can be applied to explain the source's intranight behavior.

The softening process observed in the steep decay phase of early X-ray afterglows of Swift bursts has remained a puzzle since its discovery. The softening process can also be observed in the later phase of the bursts and its cause has also been unknown. Recently, it was suggested that, influenced by the curvature effect, emission from high latitudes would shift the Band function spectrum from higher energy band to lower band, and this would give rise to the observed softening process accompanied by a steep decay of the flux density. The curvature effect scenario predicts that the terminating time of the softening process would be correlated with the duration of the process. In this paper, based on the data from the UNLV GRB group web-site, we found an obvious correlation between the two quantities. In addition, we found that the softening process can be divided into two classes: the early type softening ($t_{s,max}\leq "4000"s$) and the late type softening ($t_{s,max} > "4000"s$). The two types of softening show different behaviors in the duration vs. terminating time plot. In the relation between the variation rates of the flux density and spectral index during the softening process, a discrepancy between the two types of softening is also observed. According to their time scales and the discrepancy between them, we propose that the two types are of different origins: the early type is of internal shock origin and the late type is of external shock origin. The early softening is referred to the steep decay just following the prompt emission, whereas the late decay typically conceives the transition from flat decay to late afterglow decay. We suspect that there might be a great difference of the Lorentz factor in two classes which is responsible for the observed discrepancy.

In the present paper, we have reported the result of simultaneous multi-wavelength observations of the TeV blazar Mrk 421 during February $-$ March 2003. In this period, we have observed Mrk 421 using Pachmarhi Array of \v{C}erenkov Telescopes (PACT) of Tata Institute of Fundamental Research at Pachmarhi, India. Other simultaneous data were taken from the published literature and public data archives. We have analyzed the high quality X-ray (2-20 keV) observations from the NASA Rossi X-Ray Timing Explorer (RXTE). We have seen a possible correlated variability between X-ray and J band (1.25 $\mu$) near infrared (NIR) wavelength. This is the first case of X-ray and NIR correlated variability in Mrk 421 or any high energy peaked (HBL) blazar. The correlated variability reported here is indicating a similar origin for NIR and X-ray emission. The emission is not affected much by the environment of the surrounding medium around the central engine of the Mrk 421. The observations are consistent with the shock-in-jet model for the emission of radiations.

We selected a sample of a dozen blazars which are the prime candidates for simultaneous multi-wavelength observing campaigns in their outburst phase. We searched for optical outbursts, intra-day variability and short term variability in these blazars. We carried out optical photometric monitoring of nine of these blazars in 13 observing nights during our observing run October 27, 2006 - March 20, 2007 by using the 1.02 meter optical telescope. From our observations, our data favor the hypothesis that three blazars were in the outburst state; one blazar was in the post outburst state; three blazars were in the pre/post outburst state; one blazar was in the low-state; and the state of one blazar was not known because there is not much optical data available for the blazar to compare with our observations. Out of three nights of observations of AO 0235+164, intra-day variability was detected in two nights. Out of five nights of observations of S5 0716+714, intra-day variability was detected in two nights. In one night of observations of PKS 0735+178, intra-day variability was detected. Out of six nights of observations of 3C 454.3, intra-day variability was detected in three nights. No intra-day variability was detected in S2 0109+224, OJ 287, ON 231, 3C 279 and 1ES 2344+514 in their 1, 4, 1, 2 and 1 nights of observations respectively. AO 0235+164, S5 0716+714, OJ 287, 3C 279 and 3C 454.3 were observed in more than one night and short term variations in all these blazars were also noticed. From our observations and the available data, we found that the predicted optical outburst with the time interval of ~ 8 years in AO 0235+164 and ~ 3 years in S5 0716+714 have possibly occurred.

We report V, R, and I band CCD photometry of the radio galaxy 3C 390.3 obtained with the 1.56-m telescope of the Shanghai Astronomical Observatory from March 1995 to August 2004. Combining these data with data from the literature, we have constructed a historical light curve from 1894 to 2004 and searched for periodicities using the CLEANest program. We find possible periods of 8.30+-1.17, 5.37+-0.49, 3.51+-0.21, and 2.13+-0.08 years.

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.

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.

We have investigated the gamma-ray bursts pulses assumed to arise from relativistically expending fireballs. We find, due to the curvature effect, the evolutionary curve of the raw hardness ratio (when the background count is included) varies significantly, depending on the hardness of a burst. For a soft burst, an upside-down pulse of its raw hardness ratio curve (HRC) would be observed; for a hard burst, its raw HRC shows a pulse-like profile with a sinkage in its decaying phase; for a very hard burst, the raw HRC possesses a pulse-like profile without a sinkage in its decaying phase. For a pulse-like raw HRC as shown in the case of the hard and very hard bursts, its peak would appear in advance of that of the corresponding light curve. We have studied here the HRC of GRB 920216, GRB 920830 and GRB 990816 in detail. The features of the raw HRC expected in the hard burst are observed in these bursts. A fit to the three bursts shows that the curvature effect alone could indeed account for the predicted characteristics of HRCs. In addition, we find that the observed hardness ratio tends to be harder at the beginning of the pulses than what the curvature effect could predict and be softer at the late time of the pulses. We believe this is an evidence showing the existence of intrinsic hard-to-soft radiation.

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.

We determined the basic parameters, such as the central black hole mass ($M$), the boosting factor (or Doppler factor) ($\delta$), the propagation angle ($\Phi$) and the distance along the axis to the site of $\gamma$-ray production ($d$) for 23 $\gamma$-ray-loud blazars using their available variability timescales. In this method, the absorption effect depends on the $\gamma$-ray energy, emission size and property of the accretion disk. Using the intrinsic $\gamma$-ray luminosity as a fraction $\lambda$ of the Eddington luminosity, $L^{in}_{\gamma}=\lambda L_{Ledd.}$ and the optical depth equal to unity, we can determine the upper limit of the central black hole masses. We found that the black hole masses range between $10^{7}M_{\odot}$ and $10^{9}M_{\odot}$ when $\lambda$ = 0.1 and 1.0 are adopted. Since this method is based on gamma-ray emissions and the short time-scale of the sources, it can also be used for central black hole mass determination of high redshift gamma-ray sources. In the case of the upper limit of black hole mass there is no clear difference between BLs and FSRQs, which suggests that the central black hole masses do not play an important role in the evolutionary sequence of blazars.

Recent analysis of the long term radio lightcurve of the extremely variable BL Lacertae object AO 0235+16 by Raiteri et al. (2001) have revealed the presence of recurrent outbursts in this source with a period of ~ 5.7 +/- 0.5 yr. Periodicity analysis of the optical lightcurve also show evidence for a shorter period. Here we discuss whether such a behavior can be explained by a binary black hole model where the accretion disk of one of the supermassive black holes is precessing due to the tidal effects of the companion. We estimate the mass of the accreting hole and analyze under what constraints onto the secondary mass and the orbital parameters of the system it is possible to provide a viable interpretation of the available multiwavelength data.

BL Lacertae objects are an extreme subclass of AGNs showing rapid and large-amplitude variability, high and variable polarization, and core-dominated radio emissions. If a strong beaming effect is the cause of the extreme observation properties, one would expect that these properties would be correlated with each other. Based on the relativistic beaming model, relationships between the polarization and the magnitude variation in brightness, as well as the core- dominance parameter are derived and used statistically to compare with the observational data of a BL Lacertae object sample. The statistical results are consistent with these correlations, which suggests that the polarization, the variation, and the core-dominance parameter are possible indications of the beaming effect.

Based on the $\gamma$-ray data of blazars in the third EGRET catalog and radio data at 5 GHz, we studied the correlation between the radio and $\gamma$-ray luminosities using two statistical methods. The first method was the partial correlation analysis method, which indicates that there exist correlations between the radio and $\gamma$-ray luminosities in both high and low states as well as in the average case. The second method involved a comparison of expected $\gamma$-ray luminosity distribution with the observed data using the Kolmogorov-- Smirnov (KS) test. In the second method, we assumed that there is a correlation between the radio and $\gamma$-ray luminosities and that the $\gamma$-ray luminosity function is proportional to the radio luminosity function. The KS test indicates that the expected gamma-ray luminosity distributions are consistent with the observed data in a reasonable parameter range. Finally, we used different $\gamma$-ray luminosity functions to estimate the possible 'observed' $\gamma$-ray luminosity distributions by GLAST.

We present the optical (V, R, and I) photometry for BL Lacertae, which was observed from 1997 through 1999, with the 1.56-m telescope at the Shanghai astronomical observatory (SHAO). After the 1997 outburst, it dimmed to a low state and then brightened again. During the period JD 2450701 to JD 2450701.5, variations of 0.40mag, 0.27mag, and 0.21mag over a time scale of 100 minutes were found for V, R, and I bands, suggesting that the variations were decreasing with wavelength. The correlation between V, R, and I is also analyzed using the DCF (Discrete Correlation Function) method. This shows that the variability in the V, R, and I bands are correlated with no time delay longer than 0.2 day.

The relation between the $\gamma$-ray and the emission line luminosities for a sample of 36 $\gamma$-ray loud blazars is investigated; an apparent correlation between them, $L_{\gamma} \propto L_{Line}^{0.69\pm0.11}$, with a correlation coefficient $r=0.741$ and a chance probability of $p = 1.9\times10^{-6}$, is found. It is found, however, that there is no intrinsic correlation between them: the apparent correlation is due to the redshift dependence in a flux-limited sample. Thus no evidence is found to support the argument that the up-scattered soft photons are from the broad emission lines. Our analysis does not conflict with the SSC model. The disk-jet symbiosis and radio/$\gamma$-ray correlation found in the literature are also discussed. The radio/$\gamma$-ray correlation may be an apparent correlation caused by the boosting effect since both bands are strongly beamed.

In this paper, the post-1977 photometric observations of PKS 2155-304 are compiled and used to discuss the variation periodicity. Largest amplitude variations ($\Delta U = 1^{m}.5$; $\Delta B = 1^{m}.65$; $\Delta V = 1^{m}.85 $; $\Delta R = 1^{m}.25$; $\Delta I = 1^{m}.14 $) and color indices ($(B-V) = 0.30\pm 0.06$; $(U-B) = -0.72\pm 0.08$; $(B-R) = 0.62\pm 0.07$; $(V-R) = 0.32\pm 0.04$) are found. The Jurkevich's method and DCF (Discrete Correlation Function) method indicate possible periods of 4.16-year and 7.0-year in the V light curve.

In this paper, the compiled long-term optical and infrared measurements of some blazars are used to analyze the variation properties and the optical data are used to search for periodicity evidence in the lightcurve by means of the Jurkevich technique and the discrete correlation function (DCF) method. Following periods are found: 4.52-year for 3C 66A; 1.56 and 2.95 years for AO 0235+164; 14.4, 18.6 years for PKS 0735+178; 17.85 and 24.7 years for PKS 0754+100; 5.53 and 11.75 for OJ 287. 4.45, and 6.89 years for PKS 1215; 9 and 14.84 years for PKS 1219+285; 2.0, 13.5 and 22.5 for 3C273; 7.1 year for 3C279; 6.07 for PKS 1308+326; 3.0 and 16.5 years for PKS 1418+546; 2.0 and 9.35 years for PKS 1514-241; 18.18 for PKS 1807+698; 4.16 and 7.0 for 2155-304; 14 and 20 years for BL Lacertae. Some explanations have been discussed.

Both observations and theoretical models suggest that the emissions in gamma-ray bursts (GRBs) and the afterglows are beamed. We argue that the recent polarization measured in the afterglows gives further evidence of beaming in GRBs. In this approach, we adopted the polarization-magnitude relation of BL Lacertae objects to 4 GRBs with available polarization measurements and found that the data of the 4 GRBs are consistent with the relation of BL Lacertae objects. This result suggests that the emissions in GRBs are beamed and there are similarities between GRBs and BL Lacertae objects (or generally blazars). We suggest that the weak polarization of GRBs may result from low percentage of orderly magnetic field in the shock front.

In this paper, a method is proposed to determine the basic properties of $\gamma$-ray loud blazars, among them the central black hole mass, M, the Doppler factor, $\delta$, the propagation angle of the $\gamma$-rays with respect to the symmetric axis of a two-temperature accretion disk, $\Phi$, and the distance (i.e. the height above the accretion disk), d at which the $\gamma$-rays are created, for seven $\gamma$-ray loud blazars with available GeV variability timescales and in which the absorption effect of a $\gamma$-ray and the beaming effect have been taken into account. Our results indicate that, if we take the intrinsic $\gamma$-ray luminosity to be $\lambda$ times the Eddington luminosity, $L_{\gamma}^{in} = \lambda L_{Edd.}$, the masses of the blazars are in the range of $(4 \sim 131)\times 10^{7}M_{\odot}$, the Doppler factors ($\delta$) lie in the range of 0.57 to 5.33 the angle ($\Phi$) is in the range of $13^{\circ}$ to 43$^{\circ}$ and the distance (d) is in the range of 26R_{g} to 411R_{g}. Our model results are independent of $\gamma$-ray emission mechanisms but they do depend on the X-ray emission mechanism of the accretion disk.

In this paper, The central black hole masses and the Doppler factors are derived for PKS 0528+134, PKS 0537-441, 3C279, PKS 1406-074, PKS 1622-297, Q1633+382, Mkn 501, and BL Lacertae. The masses obtained are in the range of (1 -7)$\times 10^{7}M_{\odot}$ and compared with that obtained with the Klein-Nishina cross section considered (Dermer & Gehrels 1995). If we considered only the Thomson cross section, the masses are in the range of 2.6$\times 10^{6}M_{\odot}$ - 2$\times 10^{11}M_{\odot}$. The masses obtained from our method are less sensitive to the flux than those obtained from Dermer & Gehrels (1995) method. The masses obtained from two flares (1991 and 1996 flares) of 3C279 are almost the same. For 3C279 and BL Lacertae, viewing angle, $\theta$, and Lorentz factor, $\Gamma$, are estimated from the derived Doppler factor and the measured superluminal velocity. For 3C279, $\theta = 10^{\circ}.9-15^{\circ}.6$, $\Gamma$ = 2.4-14.4 for $\delta$ = 3.37; $\theta = 8^{\circ}.45-9^{\circ}.7$, $\Gamma$ = 2.95-11.20 for $\delta$ = 4.89; For BL Lacertae, $\theta = 25^{\circ}-29^{\circ}.4$, $\Gamma$ = 2.0-4.0.

In this paper, the historical infrared (JHK) data compiled from the published literature are presented here for 30 blazars. Maximum near-IR variations are found and compared with the optical ones. For the color-magnitude relation, some objects show that the color index increases with magnitude while 1253-055 shows complex behaviour, which perhaps suggests that the emission mechanism consists of, at least, two parts. The spectral indexes are in the range of $\alpha_{IR}$ = 0.77 to 2.37.

In this paper, the historical infrared (JHK) data compiled from the published literature are presented in electronic form for 40 radio selected BL Lacertae objects (RBLs) for the first time. Largest variations are found and compared with the largest optical variation. Relations between color index and magnitude and between color-color indices are discussed individually. For the color-magnitude relation, some objects (0048-097, 0735+178, 0851+202, 1215+303, 1219+285, 1749+096, and perhaps 0219+428, 0537-441, 1514-241) show that color index increases with magnitude indicating that the spectrum flattens when the source brightens while some other objects (0754+100, 1147+245, 1418+546, and 1727+502) perhaps show an opposite behaviour, while remaining objects do not show any clear tendency; For color-color relation, it is found common that (J-K) is closely correlated with (J-H) while (J-H) is not correlated with (H-K).

In this paper, the historical infrared (JHK) data compiled from 22 publications are presented for 3C279 and 4C29.45. Relation between color index and magnitude and between color-color indices are discussed. (J-K) is found closely correlated with (J-H) while (J-H) is almost not correlated with (H-K), which is perhaps due to two facts of (1) (J-K) shows wider distribution than (J-H) and (H-K) and (2) the spectrum deviates from the power law. For color-magnitude relation, there is a correlation of color index increasing with magnitude for 3C279 suggesting that the spectrum flattens when the source brightens. But an opposite tendency is also noted between (H-K) and J when I is fainter than 14 magnitude, which suggests that the near-IR spectrum of 3C279 consists of, at least, two components. No similar result is found for 4C29.45. The largest amplitude infrared variation is smaller than but comparable with the largest amplitude variation for the two objects.

In this paper, historical optical(UBVRI) data and newly observed data from the Yunnan Observatory of China(about100 years) are presented for BL Lacertae. Maximum variations in UBVRI: 5.12, 5.31, 4.73, 2.59, and 2.54 and color indices of U-B = -0.11 +/- 0.20, B-V= 1.0 +/- 0.11, V-R= 0.73 +/- 0.19, V-I= 1.42 +/- 0.25, R-I= 0.82 +/- 0.11, and B-I= 2.44 +/- 0.29 have been obtained from the literature; The Jurkevich method is used to investigate the existence of periods in the B band light curve, and a long-term period of 14 years is found. The 0.6 and 0.88 year periods reported by Webb et al.(1988) are confirmed. In addition, a close relation between B-I and B is found, suggesting that the spectra flattens when the source brightens.

In this paper, the correlation between the gamma-ray and the radio bands is investigated. The results show that there is a closer correlation between the gamma-ray emission and the high frequency (1.3mm, 230GHz) radio emission for maximum data than between the gamma-ray and the lower frequency (5GHz) radio emissions, which means that the gamma-ray is associated with the radio emission from the jet.

In this paper: 1. The assumption of the dependence of Doppler factor on the emission frequency (Fan et al. 1993) is used to explain the observational differences between the radio-selected BL Lac objects(RBLs) and the X-Ray-Selected BL Lac objects (XBLs): a) Hubble relation; b) different multiwavelength correlations; c) different regions in the effective spectral index diagram; d) different polarization. The results suggest that RBLs and XBLs are the same. 2. From the analysis of the relation between infrared magnitude and redshift, it is proposed that the parent population of BL Lac objects should be FRI radio galaxies and FRII(G) radio galaxies showing the optical spectra of a galaxy. 3. From the superluminal motion, the assumption (Fan et al. 1993) is confirmed. 4. Based on the relation between polarization and Doppler factor (Fan et al. 1997, A&A, 327, 947), it is proposed that the ratio of the beamed luminosity to the unbeamed luminosity in the source frame of OVVs/HPQs is smaller than that of BL Lac objects.

In this paper, we will discuss the long-term variations in the optical and the infrared bands of some AGNs. It is interesting to note that the reported periods of AGNs are of the similar values (about 10 years, see Fan et al. 1998a, ApJ, 507 in press, and references therein). DCF method shows the optical and infrared bands are strongly correlated suggesting that the emission mechanisms in the two bands have a common origin.

In this paper, the long-term historical optical (UBVRI) and near-infrared (JHK) data are presented with some new observations in the optical (February 1994-January 1995) and near-infrared (November 1995) bands included for BL Lac object OJ287. The new optical data in V-band are in agreement with the results reported by other authors (Sillanpaa et al. 1996a; Arimoto et al. 1997), a close correlation between the color index of B-V and the magnitude V has been obtained from our new observations; The new infrared observations presented here indicate that the source was at a high level in the infrared band during the observation period; From the available literature, we have got that the largest variations for UBVRIJHK bands are respectively: 4.72mag, 5.93mag, 5.18mag, 4.45mag, 4.07mag, 3.87mag, 3.78mag, and 3.54mag. A strong correlation is found between the optical and near-infrared bands when the DCF method is used, which suggests that these two bands have the same emission mechanism.