Massive bosons, such as light scalars and vector bosons, can lead to instabilities of rotating black holes by the superradiance effect, which extracts energy and angular momentum from rapidly-rotating black holes effectively. This process results in spinning-down of black holes and the formation of boson clouds around them. In this work, we used the masses and spins of supermassive black holes measured from the ultraviolet/optical or X-ray observations to constrain the model parameters of the light bosons. We find that the mass range of light bosons from $10^{-22}$ eV to $10^{-17}$ eV can be largely excluded by a set of supermassive black holes (including also the extremely massive ones OJ 287, Ton 618 and SDSS J140821.67+025733.2), particularly for the vector boson scenario, which eliminates a good fraction of the so-called fuzzy dark matter parameter regions. For the scalar bosons with self-interaction, most part of the mass range from $\sim 3 \times 10^{-19}$ eV to $10^{-17}$ eV with a decay constant $f_a > 10^{15}$ GeV can be excluded, which convincingly eliminate the QCD axions at these masses.

Recently the AMS-02 reported the precise measurements of the energy spectra of medium-mass compositions (Neon, Magnesium, Silicon) of primary cosmic rays, which reveal different properties from those of light compositions (Helium, Carbon, Oxygen). Here we propose a nearby source scenario, together with the background source contribution, to explain the newly measured spectra of cosmic ray Ne, Mg, Si, and particularly their differences from that of He, C, O. Their differences at high energies can be naturally accounted for by the element abundance of the nearby source. Specifically, the abundance ratio of the nearby source to the background of the Ne, Mg, Si elements is lower by a factor of $\sim1.7$ than that of the He, C, O elements. Such a difference could be due to the abundance difference of the stellar evolution of the progenitor star or the acceleration process/environment, of the nearby source. This scenario can simultaneously explain the high-energy spectral softening features of cosmic ray spectra revealed recently by CREAM/NUCLEON/DAMPE, as well as the energy-dependent behaviors of the large-scale anisotropies. It is predicted that the dipole anisotropy amplitudes below PeV energies of the Ne, Mg, Si group are smaller than that of the He, C, O group, which can be tested with future measurements.

We investigate the differences in the stellar population properties, the structure, and the environment between massive compact star-forming galaxies (cSFGs) with or without active galactic nucleus (AGN) at $2<z<3$ in the five 3D-HST/CANDELS fields. In a sample of 221 massive cSFGs, we constitute the most complete AGN census so far, identifying 66 AGNs by the X-ray detection, the mid-infrared color criterion, and/or the SED fitting, while the rest (155) are non-AGNs. Further dividing these cSFGs into two redshift bins, i.e., $2<z<2.5$ and $2.5 \leq z<3$, we find that in each redshift bin the cSFGs with AGNs have similar distributions of the stellar mass, the specific star formation rate, and the ratio of $L_{\rm IR}$ to $L_{\rm UV}$ to those without AGNs. After having performed a two-dimensional surface brightness modeling for those cSFGs with X-ray-detected AGNs (37) to correct for the influence of the central point-like X-ray AGN on measuring the structural parameters of its host galaxy, we find that in each redshift bin the cSFGs with AGNs have comparable distributions of all concerned structural parameters, i.e., the Sersic index, the 20\%-light radius, the Gini coefficient, and the concentration index, to those without AGNs. With a gradual consumption of available gas and dust, the structure of cSFGs, indicated by the above structural parameters, seem to be slightly more concentrated with decreasing redshift. At $2<z<3$, the similar environment between cSFGs with and without AGNs suggests that their AGN activities are potentially triggered by internal secular processes, such as gravitational instabilities or/and dynamical friction.

Gamma Ray Burst (GRBs) are empirically classified as long-duration GRBs (LGRBs, $>$ 2s) and short-duration GRBs (SGRBs, $<$ 2s). Physically they can be grouped into two distinct progenitor categories: those originating from collapse of massive stars (also known as Type II) and those related to mergers of compact stars (also known as Type I). Even though most LGRBs are Type II and most SGRBs are Type I, the duration criterion is not always reliable to determine the physical category of a certain GRB. Based on our previous comprehensive study of the multi-wavelength properties of long and short GRBs, here we utilize the Naive Bayes method to physically classify GRBs as Type I and Type II GRBs based on multi-wavelength criteria. It results in 0.5\% training error rate and 1\% test error rate. Moreover, there is a gap [-1.2, -0.16] in the distribution of the posterior Odds, $\log O({\rm II:I})$, the Type II to Type I probability ratio. Therefore, we propose to use ${\cal O} = \log O({\rm II:I})+0.7$ as the parameter to classify GRBs into Type I ($<0$) or Type II ($>0$). The only confirmed Type I GRB, GRB 170817A, has log $O({\rm II:I})=-10$. According to this criterion, the supernova-less long GRBs 060614 and 060505 belong to Type I, and two controversial short GRBs 090426 and 060121 belong to Type II.

Measuring high energy cosmic ray electrons/positrons (CRE) provides important means for probing the nearby galactic sources and for the dark matter (DM) indirect detection. In this work, we perform a systematical analysis of the flavor structure of DM annihilations into charged leptons based on the cosmic ray CRE spectra measured by DAMPE, Fermi-LAT, AMS-02, and CALET experiments. We study the annihilations of possible TeV scale DM particles in a nearby subhalo, which is proposed to explain the possible peak-like structure of the DAMPE CRE data. We pay special attention to the possible non-resonant excess (besides the possible peak-like structure) and demonstrate that such non-resonant excess can mainly arise from the decay of muons produced by the DM annihilations in the subhalo. With these we study the flavor composition of the lepton final states of $\chi\chi \to e^+e^-, \mu^+\mu^-, \tau^+\tau^-$ by fitting the CRE data. We further analyze the constraints on the lepton flavor composition using the Fermi-LAT $\gamma$-ray measurements. We find that the flavor composition (favored by the CRE data) is consistent with the Fermi-LAT data at relatively low Galactic latitudes, while the fraction of the final state $\tau^\pm$ is severely bounded.

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

Significant flares of GeV $\gamma$-ray emission from the Crab Nebula have been found by AGILE and Fermi-LAT years ago, indicating that extreme particle acceleration and radiation occurs in young pulsar wind nebulae. To enlarge the flare sample and to investigate their statistical properties will be very useful in understanding the nature of the $\gamma$-ray flares. In this paper, we investigate the flaring emission from the Crab Nebula with eleven year observations of the Fermi-LAT. We identify 17 significant flares in the light curve of the low-energy (synchrotron) component of the $\gamma$-ray emission. The flare rate is about 1.5 per year, without any significant change or clustering during the 11 years of the observation.The spectral energy distributions of such flares can be fitted by a steady power-law background and a flare component with an expotentially cutoff power-law spectrum. We find an approximately linear correlation between the cutoff energies and the total energy fluxes above 100 MeV. We also detect a special flare with an extremely long duration of nearly one month, occurred in October, 2018. Again the energy fluxes at different time correlate linearly with the cutoff energies for this single flare. Our results do not apparently support the Doppler boosting scenario of the flare emission, and suggest more complicated acceleration and/or cooling processes instead.

We present a study on structures and physical properties of massive ($M_* >10^{10} M_{\sun} $) compact galaxies at $1.0<z<2.0$ in five 3D--{\it HST}/CANDELS fields. Compared with the extended star-forming galaxies (eSFGs), compact star-forming galaxies (cSFGs) are found to have the lower level of star formation, and mainly distribute in the quiescent region of the {\it UVJ} diagram. The distributions of dust attenuation and S{\'e}rsic index support that the progenitors of cQGs are cSFGs, and cSFGs are at a transitional phase between eSFGs and cQGs. The prevalence of X-ray selected AGNs ($\sim 28\%$) is confirmed in the cSFGs at $1<z<2$ which indicates that the violent gas-rich processes such as merger and disk instability could drive the structure to be more compact, and trigger both star formation and black hole growth in the central regions. Our results support the "two-step" scenario that the cSFGs at $1<z<2$ are the intermediate population after compaction but before a quick quenching. Our analysis of parametric and nonparametric morphologies shows that cQGs (eQGs) are more concentrated and have less substructures than cSFGs (eSFGs), and quenching and compactness should be associated with each other. The cSFGs at $1.5<z<2$ ($1<z<1.5$) prefer to be in higher (lower) density environment, similar as cQGs (eSFGs). It suggests that merger or strong interaction might be the main driving mechanism of compaction at higher redshifts, whereas the disk instability of individual galaxies might play a more important role on the formation of cSFGs at lower redshifts.

With high spatial resolution, polarimetric imaging of a supermassive black hole, like M87$^\star$ or Sgr A$^\star$, by the Event Horizon Telescope can be used to probe the existence of ultralight bosonic particles, such as axions. Such particles can accumulate around a rotating black hole through the superradiance mechanism, forming an axion cloud. When linearly polarized photons are emitted from an accretion disk near the horizon, their position angles oscillate due to the birefringent effect when traveling through the axion background. In particular, the observations of supermassive black holes M87$^\star$ (Sgr A$^\star$) can probe the dimensionless axion-photon coupling $c = 2 \pi g_{a \gamma} f_a$ for axions with mass around $O(10^{-20})$~eV ($O( 10^{-17}$)~eV) and decay constant $f_a < O(10^{16})$ GeV, which is complimentary to other axion measurements.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

The DArk Matter Particle Explorer (DAMPE) can measure $\gamma$-rays in the energy range from a few GeV to about 10 TeV. The direction of each $\gamma$-ray is reconstructed with respect to the reference system of the DAMPE payload. In this paper, we adopt a maximum likelihood method and use the $\gamma$-ray data centered around several bright point-like sources to measure and correct the angular deviation from the real celestial coordinate system, the so called ``boresight alignment'' of the DAMPE payload. As a check, we also estimate the boresight alignment for some sets of simulation data with artificial orientation and obtain consistent results. The time-dependent boresight alignment analysis does not show evidence for significant variation of the parameters.

The spectral hardenings of cosmic ray nuclei above $\sim 200$ GV followed by softenings around 10 TV, the knee of the all-particle spectrum around PeV energies, as well as the pattern change of the amplitude and phase of the large-scale anisotropies around 100 TeV indicate the complexities of the origin and transportation of Galactic cosmic rays. It has been shown that nearby source(s) are most likely to be the cause of such spectral features of both the spectra and the anisotropies. In this work, we study the anisotropy features of different mass composition (or mass groups) of cosmic rays in this nearby source model. We show that even if the spectral features from the nearby source component is less distinctive compared with the background component from e.g., the population of distant sources, the anisotropy features are more remarkable to be identified. Measurements of the anisotropies of each mass composition (group) of cosmic rays by the space experiments such as DAMPE and HERD and the ground-based experiments such as LHAASO in the near future are expected to be able to critically test this scenario.

Precise measurements of the energy spectra of cosmic rays (CRs) show various kinds of features deviating from single power-laws, which give very interesting and important implications on their origin and propagation. Previous measurements from a few balloon and space experiments indicate the existence of spectral softenings around 10 TV for protons (and probably also for Helium nuclei). Very recently, the DArk Matter Particle Explorer (DAMPE) measurement about the proton spectrum clearly reveals such a softening with a high significance. Here we study the implications of these new measurements, as well as the groundbased indirect measurements, on the origin of CRs. We find that a single component of CRs fails to fit the spectral softening and the air shower experiment data simultaneously. In the framework of multiple components, we discuss two possible scenarios, the multiple source population scenario and the background plus nearby source scenario. Both scenarios give reasonable fits to the wide-band data from TeV to 100 PeV energies. Considering the anisotropy observations, the nearby source model is favored.

The conversion between photons and axionlike particles (ALPs) in the Milky Way magnetic field could result in the detectable oscillation phenomena in $\gamma$-ray spectra of Galactic sources. In this work, the GeV (Fermi-LAT) and TeV (MAGIC/VERITAS/H.E.S.S.) data of three bright supernova remnants (SNRs, ie. IC443, W51C and W49B) have been adopted together to search such the oscillation effect. Different from our previous analysis of the sole Fermi-LAT data of IC443, we do not find any reliable signal for the photon-ALP oscillation in the joint broadband spectrum of each SNR. The reason for the inconsistence is that in this work we use the latest revision (P8R3) of Fermi-LAT data, updated diffuse emission templates and the new version of the source catalog (4FGL), which lead to some modification of the GeV spectrum of IC443. Then we set constraints on ALP parameters based on the combined analysis of all the three sources. Though these constraints are somewhat weaker than limits from the CAST experiment and globular clusters, they are supportive of and complementary to these other results.

The origin of Galactic cosmic rays (GCRs) remains a mystery after more than one century of their discovery. The diffusive propagation of charged particles in the turbulent Galactic magnetic field makes us unable to trace back to their acceleration sites. Nevertheless, nearby GCR source(s) may leave imprints on the locally measured energy spectra and the anisotropies of the arrival direction. In this work we propose a simple but natural description of the GCR production and propagation, within a two-zone disk-halo diffusion scenario together with a nearby source, to understand the up-to-date precise measurements of the energy spectra and anisotropies of GCRs. We find that a common energy scale of $\sim100$ TeV appears in both energy spectra of protons and Helium nuclei measured recently by CREAM and large-scale anisotropies detected by various experiments. These results indicate that one or more local sources are very likely important contributors to GCRs below $100$ TeV. This study provides a probe to identify source(s) of GCRs by means of joint efforts of spectral and anisotropy measurements.

The spectra of emission-line galaxies (ELGs) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digit Sky Survey (SDSS) are used to study the mass-metallicity relation (MZR) at $z\sim0.8$. The selected sample contains about 180,000 massive star-forming galaxies with $0.6 < z < 1.05$ and $9 < {\rm log}(M_{\star}/M_{\odot}) < 12$. The spectra are stacked in bins of different parameters including redshift, stellar mass, star formation rate (SFR), specific star formation rate (sSFR), half-light radius, mass density, and optical color. The average MZR at $z\sim0.83$ has a downward evolution in the MZR from local to high-redshift universe, which is consistent with previous works. At a specified stellar mass, galaxies with higher SFR/sSFR and larger half-light radius have systematically lower metallicity. This behavior is reversed for galaxies with larger mass density and optical color. Among the above physical parameters, the MZR has the most significant dependency on SFR. Our galaxy sample at $0.6<z<1.05$ approximately follows the fundamental metallicity relation (FMR) in the local universe, although the sample inhomogeneity and incompleteness might have effect on our MZR and FMR.

Using the integral field unit (IFU) data from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we study the kinematics of gas and stellar components in an edge-on Seyfert 2 galaxy, SDSS J171359.00+333625.5, with X-shaped bi-conical outflows. The gas and stars therein are found to be counter-rotating, indicating that the collision between the inner and external gas might be an effective way to dissipate the angular momentum, which leads to remarkable gas accretion into the galaxy center. Large [OIII]$\lambda$5007 equivalent width and AGN-like line ratio in the large bi-conical region suggest that the gas is ionized by the central AGN. The gas velocity in the bi-cone region shows that ionized gas is receding relative to the galaxy center, which could be the joint effect of inflows, outflows and disk rotation. We are probably witnessing the case where a great amount of gas in the disk is being efficiently accreted into the central black hole, and the AGN-driven galactic winds are blown out along the bi-cone. The kinematics of oxygen, including rotation velocity and velocity dispersion, is different from other elements, like hydrogen, nitrogen and sulfur. The rotation velocity estimated from oxygen is slower than from other elements. The velocity dispersion of other elements follows galactic gravitational potential, while the velocity dispersion of oxygen stays roughly constant along the galactic major-axis. The further advanced observations, e.g. of cold gas or with an IFU of higher spatial resolution, are required to better understand this object.

In this Letter, properties of black hole (BH) mass are well checked for the interesting object SDSS J0159, a changing-look AGN and also a host galaxy of a tidal disruption event (TDE). Through spectral absorption features, the stellar velocity dispersion of SDSS J0159 can be well measured as $\sigma\sim81~{\rm km/s}$, leading to SDSS J0159 being an apparent outlier in the \msig space, because of the BH mass estimated through the \msig relation about two magnitudes lower than the reported virial BH mass of about $10^8~{\rm M_\odot}$. After considerations of contributions of stellar debris from the central TDE to broad line emission clouds, the over-estimated virial BH mass could be well explained in SDSS J0159. Therefore, over-estimated virial BH masses through broad line properties in the \msig space could be treated as interesting clues to support central TDEs.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

To investigate the mass dependence of structural transformation and star formation quenching, we construct three galaxy samples using massive ($M_* > 10^{10} M_{\odot}$) red, green, and blue galaxy populations at $0.5<z<2.5$ in five 3D--{\it HST}/CANDELS fields. The structural parameters, including effective radius ($r_{\rm e}$), galaxy compactness ($\Sigma_{1.5}$), and second order moment of 20\% brightest pixels ($M_{20}$) are found to be correlated with stellar mass. S\'{e}rsic index ($n$), concentration ($C$), and Gini coefficient ($G$) seem to be insensitive to stellar mass. The morphological distinction between blue and red galaxies is found at a fixed mass bin, suggesting that quenching processes should be accompanied with transformations of galaxy structure and morphology. Except for $r_e$ and $\Sigma_{1.5}$ at high mass end, structural parameters of green galaxies are intermediate between red and blue galaxies in each stellar mass bin at $z < 2$, indicating green galaxies are at a transitional phase when blue galaxies are being quenched into quiescent statuses. The similar sizes and compactness for the blue and green galaxies at high-mass end implies that these galaxies will not appear to be significantly shrunk until they are completely quenched into red QGs. For the green galaxies at $0.5<z<1.5$, a morphological transformation sequence of bulge buildup can be seen as they are gradually shut down their star formation activities, while a faster morphological transformation is verified for the green galaxies at $1.5<z<2.5$.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

In order to search for the line-like signals in the Fermi-LAT data, we have analyzed totally 49152 regions of interest (ROIs) that cover the whole sky. No ROI displays a line signal with test statistic (TS) value above 25, while for 50 ROIs weak line-like excesses with ${\rm TS}>16$ are presented. The intrinsic significances of these potential signals are further reduced by the large trial factor introduced in such kind of analysis. For the largest TS value of 24.3 derived in our analysis, the corresponding global significance is only $0.54\sigma$. We thus do not find any significant line-like signal and set up constraints on the cross section of dark matter annihilating to gamma-ray lines, $\left<\sigma{v}\right>_{\gamma\gamma}$.

Based on a large sample of massive ($M_{*}\geq 10^{10} M_{\odot}$) compact galaxies at $1.0 < z < 3.0$ in five 3D-HST/CANDELS fields, we quantify the fractional abundance and comoving number density of massive compact galaxies as a function of redshift. The samples of compact quiescent galaxies (cQGs) and compact star-forming galaxies (cSFGs) are constructed by various selection criteria of compact galaxies in literatures, and the effect of compactness definition on abundance estimate is proved to be remarkable, particularly for the cQGs and cSFGs at high redshifts. Regardless of the compactness criteria adopted, their overall redshift evolutions of fractional abundance and number density are found to be rather similar. Large samples of the cQGs exhibit a sustaining increase in number density from $z \sim 3$ to 2 and a plateau at $1<z<2$. For massive cSFGs, a plateau in the number density at $2<z<3$ can be found, as well as a continuous drop from $z \sim 2$ to 1. The evolutionary trends of the cQG and cSFG abundances support the scenario that the cSFGs at $z \geq 2$ may have been rapidly quenched into quiescent phase via violent dissipational processes such as major merger and disk instabilities. Rarity of the cSFGs at lower redshifts ($z < 1$) can be interpreted by the decrease of gas reservoirs in dark matter halos and the consequent low efficiency of gas-rich dissipation.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

A class of dark photon dark matter models with ultralight masses would lead to oscillation of a test body through a coupling with baryons or $B-L$ charge. This periodical oscillation of an observer results in swing of a star's apparent position due to the effect of aberration of light, which could be probed with high-precision astrometry observations of stars in the Milky Way. We propose to use the observations of stellar positions of a number of stars by Gaia to search for this kind of dark photon dark matter. We show that this astrometry method is able to give promising sensitivities to search for the dark photon dark matter in the mass range of $10^{-23}\sim10^{-21}$ eV.

GeV gamma ray is an important observation target of DArk Matter Particle Explorer (DAMPE) for indirect dark matter searching and high energy astrophysics. We present in this work a set of accurate instrument response functions of DAMPE (DmpIRFs) including the effective area, point-spread function and energy dispersion that are crucial for the gamma-ray data analysis based on the high statistics simulation data. A dedicated software named DmpST is developed to facilitate the scientific analyses of DAMPE gamma-ray data. Considering the limited number of photons and the angular resolution of DAMPE, the maximum likelihood method is adopted in the DmpST to better disentangle different source components. The basic mathematics and the framework regarding this software are also introduced in this paper.

The $\gamma$-ray observation of interstellar gas provides a unique way to probe the cosmic rays (CRs) outside the solar system. In this work, we use an updated version of Fermi-LAT data and recent multi-wavelength tracers of interstellar gas to re-analyze a mid-latitude region in the third Galactic quadrant and estimate the local CR proton spectrum. Two $\gamma$-ray production cross section models for $pp$ interaction, the commonly used one from Kamae et al. (2006) and the up-to-date one from Kafexhiu et al. (2014), are adopted separately in the analysis. Both of them can well fit the emissivity and the derived proton spectra roughly resemble the direct measurements from AMS-02 and Voyager 1, but rather different spectral parameters are indicated. A break at $4\pm1~{\rm GeV}\;c^{-1}$ is shown if the cross section model by Kamae et al. (2006) is adopted. The resulting spectrum is $\lesssim 20\%$ larger than the AMS-02 observation above $15~\rm GeV$ and consistent with the de-modulated spectrum within $2\%$. The proton spectrum based on the cross section model of Kafexhiu et al. (2014) is about $1.4-1.8$ times that of AMS-02 at $2-100~\rm GeV$, however the difference decreases to $20\%$ below $10~\rm GeV$ with respect to the de-modulated spectrum. A spectral break at $20\pm11~{\rm GeV}\;c^{-1}$ is required in this model. An extrapolation down to $300~\rm MeV$ is performed to compare with the observation of Voyager 1, and we find a deviation of $\lesssim 2.5\sigma$ for both the models. In general, an approximately consistent CR spectrum can be obtained using $\gamma$-ray observation nowadays, but we still need a better $\gamma$-ray production cross section model to derive the parameters accurately.

Based on a large sample of 254 220 galaxies in 81 089 groups, which are selected from the spectroscopic galaxy sample of the SDSS DR12, we investigate the radial distribution of incidences, morphologies, environmental densities, and star formation properties of the active galactic nucleus (AGN) host galaxies and star-forming galaxies (SFGs) in the groups at z<0.2, as well as their changes with group richness ($N_{\rm rich}$). It is found that AGN fraction slightly declines with richness for the groups/clusters. The SFG fraction is on average about 2 times larger than the AGN fraction, with a significant declining trend with richness. The group AGNs are preferentially reside in spheroidal and bulge-dominated disc galaxies, whereas the majority of SFGs are late-type discs. Compared with the SFGs, the AGNs in poor groups ($5 \leqslant N_{\rm rich} \leqslant 10$) are closer to group center. The AGN fraction does not change with the distance to the group center, whereas the SFG fraction tends to be higher in the outskirts. The AGNs in groups have a higher incidence than the SFGs for the massive ($\log(M_*/M_{\odot}) > 10.7$) galaxies, and the mean SFG fraction is about 6 times as that of AGNs in the late-type galaxies with lower masses at larger radius. The distribution of environmental luminosity densities shows that the AGNs are likely to be reside in a denser environment relative to the SFGs. Compared with the SFGs in groups, the group AGNs are found to have a higher mean stellar mass, a lower mean star formation rate, and an older mean stellar age.

We study the propagation and injection models of cosmic rays using the latest measurements of the Boron-to-Carbon ratio and fluxes of protons, Helium, Carbon, and Oxygen nuclei by the Alpha Magnetic Spectrometer and the Advanced Composition Explorer at top of the Earth, and the Voyager spacecraft outside the heliosphere. The ACE data during the same time interval of the AMS-02 data are extracted to minimize the complexity of the solar modulation effect. We find that the cosmic ray nucleus data favor a modified version of the diffusion-reacceleration scenario of the propagation. The diffusion coefficient is, however, required to increase moderately with decreasing rigidity at low energies, which has interesting implications on the particle and plasma interaction in the Milky Way. We further find that the low rigidity ($<$ a few GV) injection spectra are different for different compositions. The injection spectra are softer for lighter nuclei. These results are expected to be helpful in understanding the acceleration process of cosmic rays.

One of the main purposes of the DArk Matter Particle Explorer (DAMPE) is to measure the cosmic ray nuclei up to several tens of TeV or beyond, whose origin and propagation remains a hot topic in astrophysics. The Plastic Scintillator Detector (PSD) on top of DAMPE is designed to measure the charges of cosmic ray nuclei from H to Fe and serves as a veto detector for discriminating gamma-rays from charged particles. We propose in this paper a charge reconstruction procedure to optimize the PSD performance in charge measurement. Essentials of our approach, including track finding, alignment of PSD, light attenuation correction, quenching and equalization correction are described detailedly in this paper after a brief description of the structure and operational principle of the PSD. Our results show that the PSD works very well and almost all the elements in cosmic rays from H to Fe are clearly identified in the charge spectrum.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

The stellar mass-metallicity relation ($M_* - Z$, MZR) indicates that the metallicities of galaxies increase with increasing stellar masses. The fundamental metallicity relation (FMR) suggests that the galaxies with higher star formation rates (SFRs) tend to have lower metallicities for a given stellar mass. To examine whether the MZR and FMR still hold at poorer metallicities and higher redshifts, we compile a sample of 35 star-forming galaxies (SFGs) at $0.6<z<0.9$ using the public spectral database ($\rm v5\_10\_0$) of emission-line galaxies from the extended Baryon Oscillation Spectroscopic Survey (eBOSS). These galaxies are identified for their significant auroral $\rm [OIII]\lambda4363$ emission line ($\rm S/N \geq 3$). With the electronic temperature metallicity calibration, we find nine SFGs are extremely metal-poor galaxies with $\rm 12 + log(O/H) \leq 7.69 \ (1/10 \ Z_\odot)$. The metallicity of the most metal-deficient galaxy is $7.35\pm 0.09$ (about 1/20 $Z_{\odot}$). Compared with the SFGs with normal metallicities in local and high redshift universe, our metal-poor SFGs have more than ten times higher SFRs at a fixed stellar mass. We create a new mass -- SFR relation for these metal-poor galaxies at $0.6<z<0.9$. Due to the higher SFRs and younger stellar ages, our metal-poor SFGs deviate from the MZR and FMR in the local universe toward lower metallicities, confirming the existence of FMR, as well as the cosmic evolution of MZR and FMR with redshift.

The precise observations of Galactic cosmic ray fluxes of the secondary family, such as Li, Be, B, are expected to have significant implications on our understanding of the cosmic ray origin and propagation. Here we employ the recent very precise measurements of those species by the Alpha Magnetic Spectrometer on the International Space Station, together with their parent species (C and O), as well as the data collected by the Voyager-1 spacecraft outside the heliosphere and the Advanced Composition Explorer, we investigate the propagation of cosmic rays in the Milky Way. We find that these new data strongly favor the existence of reacceleration of cosmic rays during the propagation rather than the convective transport. We further find that for the reacceleration model, the slope of rigidity dependence of the diffusion coefficient is close to the Kolmogorov-type of turbulence, and the spectral hardenings of both the primary and secondary particles can be well described by the injection hardening rather than the propagation hardening.

HESS J1640-465 is an extended TeV $\gamma$-ray source and its $\gamma$-ray emission whether from the shell of a supernova remnant (SNR) or a pulsar wind nebula (PWN) is still under debate. We reanalyze the GeV $\gamma$-ray data in the field of HESS J1640-465 using eight years of Pass 8 data recorded by the Fermi Large Area Telescope. An extended GeV $\gamma$-ray source positionally coincident with HESS J1640-465 is found. Its photon spectrum can be described by a power-law with an index of $1.42\pm0.19$ in the energy range of 10-500 GeV, and smoothly connects with the TeV spectrum of HESS J1640-465. The broadband spectrum of HESS J1640-465 can be well fit by a leptonic model with a broken power-law spectrum of electrons with an exponential cut-off at $\sim$ 300 TeV. The spectral properties of HESS J1640-465 are broadly consistent with the characteristics of other sources identified as PWNe, such as the correlations between high-energy luminosity ratios and the physical parameters of pulsar, including spin-down luminosity $\dot{E}$ and characteristic age $\tau_c$. All these pieces of evidence support that the $\gamma$-ray emission of HESS J1640-465 may originate from the PWN powered by PSR J1640-4631 rather than the shell of the SNR G338.3-0.0.

GRS 1915+105 has been active for more than 26 years since it was discovered in 1992. There are hundreds of RXTE pointed observations on this source, and the quasi-regular flares with a slow rise and a sharp decrease (i.e. the "heartbeat" state) were recorded in more than 200 observations. The connections among the disk/corona, jet, and the disk wind at the heartbeat state have been extensively studied. In this work, we firstly perform a statistical analysis of the light curves and the X-ray spectra to investigate this peculiar state. We calculate the parameters for heartbeat cycles, including the recurrence time, the maximum and the minimum count rate, the flare amplitude, and the cumulative radiation for each cycle. The recurrence time has a bimodal distribution ranging from $\sim 20$ to $\sim 200$ s. The minimum count rate increases with increasing recurrence time; while the maximum count rate remains nearly constant around 2 Crab. Fitting the averaged spectrum for each observation, we find the strong correlations among the recurrence time, the apparent inner radius of the accretion disk (or the color correction factor), and the (nonthermal) X-ray luminosity. We suggest that the true inner edge of the accretion disk might always extend to the marginally stable orbit, while the change in corona size should result in the observed correlations.

Recently the IceCube collaboration and 15 other collaborations reported the spatial and temporal coincidence between the neutrino event IceCube-170922A and the radio-TeV activity of the blazar TXS 0506+056. Their further analysis on 9.5 years of IceCube data discovered neutrino flare between September 2014 and March 2015, when TXS 0506+056 is however in "quiescent" state. We analyze the Fermi-LAT data in that direction, and find another bright GeV source PKS 0502+049, which is at an angle of $1.2^{\circ}$ from TXS 0506+056, with strong activties during the neutrino flare. No other bright GeV source was detected in the region of interest. Though PKS 0502+049 is $1.2^\circ$ separated from TXS 0506+056, it locates within the directional reconstruction uncertainties of 7 neutrinos, out of the 13 neutrino events during the neutrino flare. Together with the observed high flux of the $\gamma$-ray flare, it may be unreasonable to fully discard the (partial) contribution of PKS 0502+049 to the neutrino flare. The single source assumption used in the neutrino data analysis might need to be revisited.

The Plastic Scintillator Detector (PSD) of the DArk Matter Particle Explorer (DAMPE) is designed to measure cosmic ray charge (Z) and to act as a veto detector for gamma-ray identification. In order to fully exploit the charge identification potential of the PSD and to enhance its capability to identify the gamma ray events, we develop a PSD detector alignment method. The path length of a given track in the volume of a PSD bar is derived taking into account the shift and rotation alignment corrections. By examining energy spectra of corner-passing events and fully contained events, position shifts and rotations of all PSD bars are obtained, and found to be on average about 1mm and 0.0015 radian respectively. To validate the alignment method, we introduce the artificial shifts and rotations of PSD bars in the detector simulation. These shift and rotation parameters can be recovered successfully by the alignment procedure. As a result of the PSD alignment procedure, the charge resolution of the PSD is improved from $4\%$ to $8\%$ depending on the nuclei.

This paper presents the second data release (DR2) of the Beijing-Arizona Sky Survey (BASS). BASS is an imaging survey of about 5400 deg$^2$ in $g$ and $r$ bands using the 2.3 m Bok telescope. DR2 includes the observations as of July 2017 obtained by BASS and Mayall $z$-band Legacy Survey (MzLS). This is our first time to include the MzLS data covering the same area as BASS. BASS and MzLS have respectively completed about 72% and 76% of their observations. The two surveys will be served for the spectroscopic targeting of the upcoming Dark Energy Spectroscopic Instrument. Both BASS and MzLS data are reduced by the same pipeline. We have updated the basic data reduction and photometric methods in DR2. In particular, source detections are performed on stacked images, and photometric measurements are co-added from single-epoch images based on these sources. The median 5$\sigma$ depths with corrections of the Galactic extinction are 24.05, 23.61, and 23.10 mag for $g$, $r$, and $z$ bands, respectively. The DR2 data products include stacked images, co-added catalogs, and single-epoch images and catalogs. The BASS website (http://batc.bao.ac.cn/BASS/) provides detailed information and links to download the data.

The DArk Matter Particle Explorer (DAMPE) is a satellite-borne, high-energy particle and $\gamma$-ray detector, which is dedicated to indirectly detecting particle dark matter and studying high-energy astrophysics. The first results about precise measurement of the cosmic ray electron plus positron spectrum between 25 GeV and 4.6 TeV were published recently. The DAMPE spectrum reveals an interesting spectral softening around $0.9$ TeV and a tentative peak around $1.4$ TeV. These results have inspired extensive discussion. The detector of DAMPE, the data analysis, and the first results are introduced. In particular, the physical interpretations of the DAMPE data are reviewed.

The recent HAWC observations of very-high-energy $\gamma$-ray halo around Geminga and Monogem indicate a very slow diffusion of cosmic rays, which results in tiny contribution of positrons from these two pulsars to the local flux. This makes the cosmic positron excess anomaly observed by PAMELA and AMS-02 even more puzzling. However, from the Boron-to-Carbon ratio data one can infer that the average diffusion coefficient in the Galaxy should be much larger. In this work we propose a two-zone diffusion model that the diffusion is slow only in a small region around the source, outside of which the propagation is as fast as usual. We find that such a scenario can naturally explain the positron excess data with parameters even more reasonable than that in the conventional one-zone diffusion model. The reason is that during the life time of Geminga ($\sim 300$ kyr) the electrons/positrons have propagated too far away with a fast diffusion and lead to a low local flux. The slow diffusion region in the two-zone model helps to confine the electrons/positrons for a long time and lead to an enhancement of the local flux. So under the constraint of the HAWC observations, pulsars are still the probable origin of the cosmic-ray positron excess.

The acceleration of cosmic ray particles and their propagation in the Milky Way and the heliosphere tangle with each other, leading to complexity and degeneracy of the modeling of Galactic cosmic rays (GCRs). The recent measurements of the GCR spectra by Voyager-1 from outside of the heliosphere gave a first direct observation of GCRs in the local interstellar (LIS) environment. Together with the high-precision data near the Earth taken by ACE and AMS-02, we derive the LIS spectra of Helium, Lithium, Beryllium, Boron, Carbon, and Oxygen nuclei from a few MeV/n to TeV/n, using a non-parameterization method. These LIS spectra are helpful in further studying the injection and propagation parameters of GCRs. The nearly 20 years of data recorded by ACE are used to determine the solar modulation parameters over the solar cycles 23 and 24, based on the force-field approximation. We find general agreements of the modulation potential with the results inferred from neutron monitors and other cosmic ray data.

On 2017 September 22 the IceCube neutrino observatory detected a track-like, very-high-energy event (IceCube-170922A) that is spatially associated with TXS 0506+056, a quasar at a redshift of $z=0.3365$. This source is characterized by the increased acitivies in a very wide energy range (from radio to TeV) during these days. To investigate the possible connection of the PeV neutrino emission with the GeV activity of blazars, in this work we select 116 bright sources and analyze their lightcurves and spectra. We focus on the sources displaying GeV activities. Among these blazars, TXS 0506+056 seems to be typical in many aspects but is distinguished by the very strong GeV activties. We suggest to search for neutrino outburst in the historical data of IceCube, as recently done for TXS 0506+056, from the directions of these more energetic and harder blazars with strong GeV activities.

More than fifty years after the discovery of the knee in the cosmic ray spectra, its physical origin remains a mystery. This is partly due to the ambiguity of the energy spectrum of individual composition. Recently, direct measurements from several balloon/satellite-borne detectors found significant spectral hardenings of cosmic ray nuclei at a few hundred GV rigidities. A joint modeling of the direct and indirect measurements of the cosmic ray spectra may help to understand the experimental systematics and probably the physics of the knee. In this work, we update the phenomenological "poly-gonato" model of the knee to include the spectral hardenings, with a changing spectral index of $\gamma + \beta \cdot \log E$ with respect to energy. We find that there is some tension between the direct and indirect measurements. Specifically, the hard spectra of CREAM suggest a relatively low energy cutoff of the individual cosmic ray spectrum around the knee, which is inconsistent with the light component spectra of the air shower measurements. Furthermore, the light component data from AS$\gamma$/ARGO-YBJ/WFCTA/KASCADE are consistent with the all-particle spectra. We expect that more precise measurements of the cosmic ray spectrum of individual composition at TeV energies and beyond by the operating space detectors such as CALET and DAMPE, and the ground based experiments such as LHAASO, will eventually solve these discrepancies. Finally, as an illustration, we show that a spatial-dependent diffusion under two-halo model is able to reproduce the results of this "modified poly-gonato" model.

We present the photometric calibration of the Beijing-Arizona Sky Survey (BASS) and Mayall z-band Legacy Survey (MzLS), which are two of the three wide-field optical legacy imaging surveys to provide the baseline targeting data for the Dark Energy Spectroscopic Instrument (DESI) project. The method of our photometric calibration is subdivided into the external and internal processes. The former utilizes the point-source objects of Pan-STARRS1 survey (PS1) as the reference standards to achieve the zero points of the absolute flux for individual exposures. And then the latter revise the zero points to make them consistent across the survey based on multiple tilings and large offset overlaps. Our process achieves a homogeneous photometric calibration over most of the sky with precision better than 10 mmag for g and r bands, 15 mmag for z band. The accuracy of the calibration is better than 1% at the bright end (16-18 mag) over most of the survey area.

The potential antiproton excess in the AMS-02 data is of much interest and can probably come from dark matter annihilations. Based on the effective field theory approach, in this work we investigate the compatibility of the DM interpretation of the AMS-02 antiproton excess and the null results from direct detection experiments, LUX, PandaX-II, and XENON1T. We focus on dimension-five and -six operators with fermion DM. Only one of dimension-five and one of dimension-six operators can successfully account for the antiproton excess, while the rest either are excluded by direct detection or require very small cut-off scales which invalidate the effective field theory approach.

We study the cosmic ray antiprotons with updated constraints on the propagation, proton injection, and solar modulation parameters based on the newest AMS-02 data near the Earth and Voyager data in the local interstellar space, and on the cross section of antiproton production due to proton-proton collisions based on new collider data. We use a Bayesian approach to properly consider the uncertainties of the model predictions of both the background and the dark matter (DM) annihilation components of antiprotons. We find that including an extra component of antiprotons from the annihilation of DM particles into a pair of quarks can improve the fit to the AMS-02 antiproton data considerably. The favored mass of DM particles is about $60\sim100$ GeV, and the annihilation cross section is just at the level of the thermal production of DM ($\langle\sigma v\rangle \sim O(10^{-26})$ cm$^3$~s$^{-1}$).

In this work, we search for $\gamma$-ray emission from the directions of some nearby Milky Way dwarf spheroidal galaxies (dSphs) and candidates with the publicly-available Pass 8 data of Fermi-LAT. Our sample includes 12 sources with the distances $<50$ kpc. Very weak $\gamma$-ray excesses ($\sim 2\sigma$) are found in some dSphs/candidates, consistent with that reported in the previous literature. Intriguingly, the peak test statistic (TS) value of the weak emission from the direction of Reticulum II rises continually. If interpreted as dark matter (DM) annihilation, the peak TS value is 13.5 for the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ and the DM mass of $m_\chi \sim 16$ GeV. The combination of all these nearby sources yields a more significant (with local significance $> 4\sigma$) $\gamma$-ray signal.

Axion-like particles (ALPs) can oscillate to photons and vice versa in electromagnetic fields. The photon-ALP oscillation provides an attractive solution to the apparent transparency of the Universe to TeV photons. The allowed parameter regions for the ALP mass $m_{\rm a}\leq 10^{-7}$ eV have been tightly constrained by the Fermi-LAT and H.E.S.S observations of some extragalactic sources. In this work we show for the first time that the H.E.S.S observations of some TeV sources in the Galactic plane exclude the highest ALP mass region (i.e., $m_{\rm a}\sim {\rm a ~few\times 10^{-7}}$ eV) that accounts for the TeV transparency of the Universe. The upcoming CTA observations of the Galactic TeV sources are shown to be able to improve the constraints significantly.

To explore the evolutionary connection among red, green, and blue galaxy populations, based on a sample of massive ($M_* > 10^{10} M_{\odot} $) galaxies at 0.5<z<2.5 in five 3D-HST/CANDELS fields, we investigate the dust content, morphologies, structures, AGN fractions, and environments of these three galaxy populations. Green valley galaxies are found to have intermediate dust attenuation, and reside in the middle of the regions occupied by quiescent and star-forming galaxies in the UVJ diagram. Compared with blue and red galaxy populations at z<2, green galaxies have intermediate compactness and morphological parameters such as Sersic index, concentration, Gini coefficient, and the second order moment of the 20% brightest pixels of a galaxy. Above findings seem to favor the scenario that green galaxies are at transitional phase when star-forming galaxies are being quenched into quiescent status. The green galaxies at z<2 show the highest AGN fraction, suggesting that AGN feedback may have played an important role in star formation quenching. For the massive galaxies at 2<z<2.5, both red and green galaxies are found to have a similarly higher AGN fraction than the blue ones, which implies that AGN feedback may help to keep quiescence of red galaxies at z>2. A significant environmental difference is found between green and red galaxies at z<1.5. Green and blue galaxies at z>0.5 seem to have similar local density distributions, suggesting that environment quenching is not the major mechanism to cease star formation at z>0.5. The fractions of three populations as functions of mass support a "downsizing" quenching picture that the bulk of star formation in more massive galaxies is completed earlier than that of lower mass galaxies.

We propose a flavored $U(1)_{e\mu}$ neutrino mass and dark matter~(DM) model to explain the recent DArk Matter Particle Explorer (DAMPE) data, which feature an excess on the cosmic ray electron plus positron flux around 1.4 TeV. Only the first two lepton generations of the Standard Model are charged under the new $U(1)_{e\mu}$ gauge symmetry. A vector-like fermion $\psi$, which is our DM candidate, annihilates into $e^{\pm}$ and $\mu^{\pm}$ via the new gauge boson $Z'$ exchange and accounts for the DAMPE excess. We have found that the data favors a $\psi$ mass around 1.5~TeV and a $Z'$ mass around 2.6~TeV, which can potentially be probed by the next generation lepton colliders and DM direct detection experiments.

The jet breaks in the afterglow lightcurves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta ($\theta_{\rm j}$) and hence the neutron star merger rate. In this work we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B and find $\theta_{\rm j}\sim 0.1$ rad. Together with five events reported before 2015 and other three "identified" recently (GRB 050709, GRB 060614 and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated $\theta_{\rm j}$. In particular, three {\it Swift} bursts in the sample have redshifts $z\leq 0.2$, with which we estimate the local neutron star merger rate density {to be $\sim 1109^{+1432}_{-657}~{\rm Gpc^{-3}~yr^{-1}}$ or $162^{+140}_{-83} {\rm Gpc^{-3}yr^{-1}}$ if the narrowly-beamed GRB 061201 is excluded}. Inspired by the typical $\theta_{\rm j}\sim 0.1$ rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association or not. For the former the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of $\sim 10\%$ is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of $\sim 3100$ K at $\sim 3.6$ days after the burst and more data are needed to ultimately clarify.

Axionlike-particles (ALPs) are one promising type of dark matter candidate particle that may generate detectable effects on $\gamma$-ray spectra other than the canonical weakly interacting massive particles. In this work we search for such oscillation effects in the spectra of supernova remnants caused by the photon-ALP conversion, using the Fermi Large Area Telescope data. Three bright supernova remnants, IC443, W44, and W51C, are analyzed. The inclusion of photon-ALP oscillations yields an improved fit to the $\gamma$-ray spectrum of IC443, which gives a statistical significance of $4.2\sigma$ in favor of such spectral oscillation. However, the best-fit parameters of ALPs ($m_{a}=6.6\,{\rm neV}$, $g_{a\gamma}=13.4 \times 10^{-11}\,{\rm GeV}^{-1}$) are in tension with the upper bound ($g_{a\gamma}< 6.6 \times 10^{-11}\,{\rm GeV}^{-1}$) set by the CAST experiment. It is difficult to explain the results using the systematic uncertainties of the flux measurements. We speculate that the "irregularity" displayed in the spectrum of IC443 may be due to the superposition of the emission from different parts of the remnant.

Recent direct measurements of Galactic cosmic ray spectra by balloon/space-borne detectors reveal spectral hardenings of all major nucleus species at rigidities of a few hundred GV. The all-sky diffuse gamma-ray emissions measured by the Fermi Large Area Telescope also show spatial variations of the intensities and spectral indices of cosmic rays. These new observations challenge the traditional simple acceleration and/or propagation scenario of Galactic cosmic rays. In this work we propose a spatially-dependent diffusion scenario to explain all these phenomena. The diffusion coefficient is assumed to be anti-correlated with the source distribution, which is a natural expectation from the charged particle transportation in turbulent magnetic field. The spatially-dependent diffusion model also gives a lower level of anisotropies of cosmic rays, which are consistent with observations by underground muons and air shower experiments. The spectral variations of cosmic rays across the Galaxy can be properly reproduced by this model.

The axion-like particle (ALP)-photon mixing in the magnetic field around $\gamma$-ray sources or along the line-of-sight could induce oscillation between photons and ALPs, which then causes irregularities in the $\gamma$-ray spectra. In this work we try to search for such spectral irregularities in the spectrum of PKS $2155-304$ using 8.6 years of the Fermi Large Area Telescope (Fermi-LAT) data. No significant evidence for the presence of ALP-photon oscillation is obtained, and the parameter space of ALPs is constrained. The exclusion region sensitively depends on the poorly known magnetic field of host galaxy cluster of PKS $2155-304$. If the magnetic field is as high as $\sim 10~{\rm \mu G}$, the "hole"-like parameter region allowed in Ref.~\cite{Fermi} can be ruled out.

The remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions.

The $\gamma$-ray emission from molecular clouds is widely believed to have a hadronic origin, while unequivocal evidence is still lacking. In this work, we analyze the Fermi-LAT Pass 8 publicly available data accumulated from 2008 August 4 to 2017 June 29 and report the significant detection of a characteristic $\pi^0$-decay feature in the $\gamma$-ray spectra of the molecular clouds Orion A and Orion B. This detection provides a direct evidence for the hadronic origin of their $\gamma$-ray emission.

The Be X-ray pulsar SMC X-3 underwent an extra long and ultraluminous giant outburst from 2016 August to 2017 March. The peak X-ray luminosity is up to $\sim 10^{39}$ erg/s, suggesting a mildly super-Eddington accretion onto the strongly magnetized neutron star. It therefore bridges the gap between the Galactic Be/X-ray binaries ($L_{\rm X}^{\rm peak} \leq 10^{38}$ erg/s) and the ultraluminous X-ray pulsars ($L_{\rm X}^{\rm peak} \geq 10^{40}$ erg/s) found in nearby galaxies. A number of observations were carried out to observe the outburst. In this paper, we perform a comprehensive phase-resolved analysis on the high quality data obtained with the Nustar and XMM-Newton, which were observed at a high and intermediate luminosity levels. In order to get a better understanding on the evolution of the whole extreme burst, we take the Swift results at the low luminosity state into account as well. At the early stage of outburst, the source shows a double-peak pulse profile, the second main peak approaches the first one and merges into the single peak at the low luminosity. The second main peak vanishes beyond 20 keV, and its radiation becomes much softer than that of the first main peak. The line widths of fluorescent iron line vary dramatically with phases, indicating a complicated geometry of accretion flows. In contrast to the case at low luminosity, the pulse fraction increases with the photon energy. The significant small pulse fraction detected below 1 keV can be interpreted as the existence of an additional thermal component located at far away from the central neutron star.

We report the detection of GeV $\gamma$-ray emission from supernova remnant HESS J1731-347 using 9 years of {\it Fermi} Large Area Telescope data. We find a slightly extended GeV source in the direction of HESS J1731-347. The spectrum above 1 GeV can be fitted by a power-law with an index of $\Gamma = 1.77\pm0.14$, and the GeV spectrum connects smoothly with the TeV spectrum of HESS J1731-347. Either a hadronic-leptonic or a pure leptonic model can fit the multi-wavelength spectral energy distribution of the source. However, the hard GeV $\gamma$-ray spectrum is more naturally produced in a leptonic (inverse Compton scattering) scenario, under the framework of diffusive shock acceleration. We also searched for the GeV $\gamma$-ray emission from the nearby TeV source HESS J1729-345. No significant GeV $\gamma$-ray emission is found, and upper limits are derived.

The DArk Matter Particle Explorer (DAMPE), also known as Wukong in China, launched on December 17, 2015, is a new high energy cosmic ray and {\gamma}-ray satellite-borne observatory in space. One of the main scientific goals of DAMPE is to observe GeV-TeV high energy {\gamma}-rays with accurate energy, angular, and time resolution, to indirectly search for dark matter particles and for the study of high energy astrophysics. Due to the comparatively higher fluxes of charged cosmic rays with respect to {\gamma}-rays, it is challenging to identify {\gamma}-rays with sufficiently high efficiency minimizing the amount of charged cosmic ray contamination. In this work we present a method to identify {\gamma}-rays in DAMPE data based on Monte Carlo simulations, using the powerful electromagnetic/hadronic shower discrimination provided by the calorimeter and the veto detection of charged particles provided by the plastic scintillation detector. Monte Carlo simulations show that after this selection the number of electrons and protons that contaminate the selected {\gamma}-ray events at $\sim10$ GeV amounts to less than 1% of the selected sample. Finally, we use flight data to verify the effectiveness of the method by highlighting known {\gamma}-ray sources in the sky and by reconstructing preliminary light curves of the Geminga pulsar.

On 17 August 2017, a gravitational wave event (GW170817) and an associated short gamma-ray burst (GRB 170817A) from a binary neutron star merger had been detected. The followup optical/infrared observations also identified the macronova/kilonova emission (AT2017gfo). In this work we discuss some implications of the remarkable GW170817/GRB 170817A/AT2017gfo association. We show that the $\sim 1.7$s time delay between the gravitational wave (GW) and GRB signals imposes very tight constraint on the superluminal movement of gravitational waves (i.e., the relative departure of GW velocity from the speed of light is $\leq 4.3\times 10^{-16}$) or the possible violation of weak equivalence principle (i.e., the difference of the gamma-ray and GW trajectories in the gravitational field of the galaxy and the local universe should be within a factor of $\sim 3.4\times 10^{-9}$). The so-called Dark Matter Emulators and a class of contender models for cosmic acceleration ("Covariant Galileon") are ruled out, too. The successful identification of Lanthanide elements in the macronova/kilonova spectrum also excludes the possibility that the progenitors of GRB 170817A are a binary strange star system. The high neutron star merger rate (inferred from both the local sGRB data and the gravitational wave data) together with the significant ejected mass strongly suggest that such mergers are the prime sites of heavy r-process nucleosynthesis.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

Weakly interacting massive particles are a widely well-probed dark matter candidate by the dark matter direct detection experiments. Theoretically, there are a large number of ultraviolet completed models that consist of a weakly interacting massive particle dark matter. The variety of models makes the comparison with the direct detection data complicated and often non-trivial. To overcome this, in the non-relativistic limit, the effective theory was developed in the literature which works very well to significantly reduce the complexity of dark matter-nucleon interactions and to better study the nuclear response functions. In the effective theory framework for a spin-1/2 dark matter, we combine three independent likelihood functions from the latest PandaX, LUX, and XENON1T data, and give a joint limit on each effective coupling. The astrophysical uncertainties of the dark matter distribution are also included in the likelihood. We further discuss the isospin violating cases of the interactions. Finally, for both dimension-five and dimension-six effective theories above the electroweak scale, we give updated limits of the new physics mass scales.

The routinely flaring events from \sgras\ trace dynamic, high-energy processes in the immediate vicinity of the supermassive black hole. We statistically study temporal and spectral properties, as well as fluence and duration distributions, of the flares detected by the \chandra\ X-ray Observatory from 1999 to 2012. The detection incompleteness and bias are carefully accounted for in determining these distributions. We find that the fluence distribution can be well characterized by a power-law with a slope of $1.73^{+0.20}_{-0.19}$, while the durations ($\tau$ in seconds) by a log-normal function with a mean $\log(\tau)=3.39^{+0.27}_{-0.24}$ and an intrinsic dispersion $\sigma=0.28^{+0.08}_{-0.06}$. No significant correlation between the fluence and duration is detected. The apparent positive correlation, as reported previously, is mainly due to the detection bias (i.e., weak flares can be detected only when their durations are short). These results indicate that the simple self-organized criticality model has difficulties in explaining these flares. We further find that bright flares usually have asymmetric lightcurves with no statistically evident difference/preference between the rising and decaying phases in terms of their spectral/timing properties. Our spectral analysis shows that although a power-law model with a photon index of $2.0\pm0.4$ gives a satisfactory fit to the joint spectra of strong and weak flares, there is weak evidence for a softer spectrum of weaker flares. This work demonstrates the potential to use statistical properties of X-ray flares to probe their trigger and emission mechanisms, as well as the radiation propagation around the black hole.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Recent precise measurements of cosmic ray (CR) spectra show that the energy distribution of protons is softer than those of heavier nuclei, and there are spectral hardenings for all nuclear compositions above $\sim$200 GV. Models proposed for these anomalies generally assume steady-state solutions of the particle acceleration process. We show that, if the diffusion coefficient has a weak dependence on the particle rigidity near shock fronts of supernova remnants (SNRs), time-dependent solutions of the linear diffusive shock acceleration at two stages of SNR evolution can naturally account for these anomalies. The high-energy component of CRs is dominated by acceleration in the free expansion and adiabatic phases with enriched heavy elements and a high shock speed. The low energy component may be attributed to acceleration by slow shocks propagating in dense molecular clouds with low metallicity in the radiative phase. Instead of a single power-law distribution, the spectra of time-dependent solutions soften gradually with the increase of energy, which may beresponsible for the "knee" of CRs.

The Be X-ray pulsar, SMC X-3 underwent a giant outburst from 2016 August to 2017 March, which was monitored with the Swift satellite. During the outburst, its broadband flux increased dramatically, and the unabsorbed X-ray luminosity reached an extreme value of $\sim 10^{39}$ erg/s around August 24. Using the Swift/XRT data, we measure the observed pulse frequency of the neutron star to compute the orbital parameters of the binary system. After applying the orbital corrections to Swift observations, we find that the spin frequency increases steadily from 128.02 mHz on August 10 and approach to spin equilibrium $\sim 128.74$ mHz in 2017 January with an unabsorbed luminosity of $L_{\rm X} \sim 2\times10^{37}$ erg/s, indicating a strong dipolar magnetic field $B \sim 6.8\times10^{12}$ G at the neutron star surface. The spin-up rate is tightly correlated with its X-ray luminosity during the super-Eddington outburst. The pulse profile in the Swift/XRT data is variable, showing double peaks at the early stage of outburst and then merging into a single peak at low luminosity. Additionally, we report that a low temperature ($kT \sim 0.2$ keV) thermal component emerges in the phase-averaged spectra as the flux decays, and it may be produced from the outer truncated disk or the boundary layer between the exterior flow and the magnetosphere.

GC-1 and GC-2 are two globular clusters (GCs) in the remote halo of M81 and M82 in the M81 group discovered by Jang et al. using the {\it Hubble Space Telescope} ({\it HST}) images. These two GCs were observed as part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey, using 14 intermediate-band filters covering a wavelength range of 4000--10000 \AA. We accurately determine these two clusters' ages and masses by comparing their spectral energy distributions (from 2267 to 20000~{\AA}, comprising photometric data in the near-ultraviolet of the {\it Galaxy Evolution Explorer}, 14 BATC intermediate-band, and Two Micron All Sky Survey near-infrared $JHK_{\rm s}$ filters) with theoretical stellar population-synthesis models, resulting in ages of $15.50\pm3.20$ for GC-1 and $15.10\pm2.70$ Gyr for GC-2. The masses of GC-1 and GC-2 obtained here are $1.77-2.04\times 10^6$ and $5.20-7.11\times 10^6 \rm~M_\odot$, respectively. In addition, the deep observations with the Advanced Camera for Surveys and Wide Field Camera 3 on the {\it HST} are used to provide the surface brightness profiles of GC-1 and GC-2. The structural and dynamical parameters are derived from fitting the profiles to three different models; in particular, the internal velocity dispersions of GC-1 and GC-2 are derived, which can be compared with ones obtained based on spectral observations in the future. For the first time, in this paper, the $r_h$ versus $M_V$ diagram shows that GC-2 is an ultra-compact dwarf in the M81 group.

Puppis A is a very famous and extensively studied supernova remnant (SNR) that shows strong evidence of shock-cloud interaction. We re-analyze the GeV $\gamma$-ray emission of it using seven years Pass 8 data recorded by the Fermi Large Area Telescope (Fermi-LAT). The morphology of the $\gamma$-ray emission is more compatible with that of the thermal X-ray and IR emissions than the radio image, which suggests a possible correlation between the gamma-ray emitting region and dense clouds. The $\gamma$-ray spectrum in the energy range of 1-500 GeV shows a break at $7.92\pm1.91$ GeV with the photon indices of $1.81\pm0.08$ below the break and $2.53\pm0.12$ above the break, which can naturally explain the lack of TeV $\gamma$-ray emission from Puppis A. The multi-wavelength observations favor a hadronic origin for the $\gamma$-ray emission.

Compared with numerous X-ray dominant active galactic nuclei (AGNs) without emission-line signatures in their optical spectra, the X-ray selected AGNs with optical emission lines are probably still in the high-accretion phase of black hole growth. This paper presents an investigation on the fraction of these X-ray detected AGNs with optical emission-line spectra in 198 galaxy groups at $z<1$ in a rest frame 0.1-2.4 keV luminosity range 41.3 <log(L_X/erg s-1) < 44.1 within the COSMOS field, as well as its variations with redshift and group richness. For various selection criteria of member galaxies, the numbers of galaxies and the AGNs with optical emission lines in each galaxy group are obtained. It is found that, in total 198 X-ray groups, there are 27 AGNs detected in 26 groups. AGN fraction is on everage less than $4.6 (\pm 1.2)\%$ for individual groups hosting at least one AGN. The corrected overall AGN fraction for whole group sample is less than $0.98 (\pm 0.11) \%$. The normalized locations of group AGNs show that 15 AGNs are found to be located in group centers, including all 6 low-luminosity group AGNs. A week rising tendency with $z$ are found: overall AGN fraction is 0.30-0.43% for the groups at $z<0.5$, and 0.55-0.64% at 0.5 < z < 1.0. For the X-ray groups at $z>0.5$, most member AGNs are X-ray bright, optically dull, which results in a lower AGN fractions at higher redshifts. The AGN fraction in isolated fields also exhibits a rising trend with redshift, and the slope is consistent with that in groups. The environment of galaxy groups seems to make no difference in detection probability of the AGNs with emission lines. Additionally, a larger AGN fractions are found in poorer groups, which implies that the AGNs in poorer groups might still be in the high-accretion phase, whereas the AGN population in rich clusters is mostly in the low-accretion, X-ray dominant phase.

Using the latest AMS-02 cosmic ray antiproton flux data, we search for potential dark matter annihilation signal. The background parameters about the propagation, source injection, and solar modulation are not assumed {\it a priori}, but based on the results inferred from the recent B/C ratio and proton data measurements instead. The possible dark matter signal is incorporated into the model self-consistently under a Bayesian framework. Compared with the astrophysical background only hypothesis, we find that a dark matter signal is favored. The rest mass of the dark matter particles is $\sim 20-80$ GeV and the velocity-averaged hadronic annihilation cross section is about $(0.2-5)\times 10^{-26}$ cm$^{3}$s$^{-1}$, in agreement with that needed to account for the Galactic center GeV excess and/or the weak GeV emission from dwarf spheroidal galaxies Reticulum 2 and Tucana III. Tight constraints on the dark matter annihilation models are also set in a wide mass region.

Tidal disruption events (TDEs) occur when a star or sub-stellar object passes close enough to a galaxy's supermassive black hole to be disrupted by tidal forces. NGC 4845 (d=17 Mpc) was host to a TDE, IGR J12580+0134, detected in November 2010. Its proximity offers us a unique close-up of the TDE and its aftermath. We discuss new Very Long Baseline Array (VLBA) and Karl G. Jansky Very Large Array (JVLA) observations, which show that the radio flux from the active nucleus created by the TDE has decayed in a manner consistent with predictions from a jet-circumnuclear medium interaction model. This model explains the source's broadband spectral evolution, which shows a spectral peak that has moved from the submm (at the end of 2010) to GHz radio frequencies (in 2011-2013) to <1 GHz in 2015. The milliarcsecond-scale core is circularly polarized at 1.5 GHz but not at 5 GHz, consistent with the model. The VLBA images show a complex structure at 1.5 GHz that includes an east west extension ~40 milliarcsec (3 pc) long as well as a resolved component 52 milliarcsec (4.1 pc) northwest of the flat-spectrum core, which is all that can be seen at 5 GHz. If ejected in 2010, the NW component must have had v=0.96 c over five years. However, this is unlikely, as our model suggests strong deceleration to speeds < 0.5c within months and a much smaller, sub-parsec size. In this interpretation, the northwest component could have either a non-nuclear origin or be from an earlier event.

In this work we use the newly reported Boron-to-Carbon ratio (B/C) from AMS-02 and the time-dependent proton fluxes from PAMELA and AMS-02 to constrain the source and propagation parameters of cosmic rays in the Milky Way. A linear correlation of the solar modulation parameter with solar activities is assumed to account for the time-varying cosmic ray fluxes. A comprehensive set of propagation models, with/without reacceleration or convection, have been discussed and compared. We find that only the models with reacceleration can self-consistently fit both the proton and B/C data. The rigidity dependence slope of the diffusion coefficient, $\delta$, is found to be about $0.38-0.50$ for the diffusion-reacceleration models. The plain diffusion and diffusion-convection models fit the data poorly. We compare different model predictions of the positron and antiproton fluxes with the data. We find that the diffusion-reacceleration models over-produce low energy positrons, while non-reacceleration models give better fit to the data. As for antiprotons, reacceleration models tend to under-predict low energy antiproton fluxes, unless a phenomenological modification of the velocity-dependence of the diffusion coefficient is applied. Our results suggest that there could be important differences of the propagation for nuclei and leptons, in either the Milky Way or the solar heliosphere.

Milky Way-like galaxies are predicted to host a very large number of dark matter subhalos. Some massive and nearby subhalos could generate detectable gamma-rays, appearing as unidentified, spatially-extended and stable gamma-ray sources. We search for such sources in the third Fermi Large Area Telescope source List (3FGL) and report the identification of a new candidate, 3FGL J1924.8-1034. With the Fermi-LAT Pass 8 data, we find that 3FGL J1924.8-1034 is spatially-extended at a high confidence level of $5.4\sigma$, with a best-fit extension radius of $\sim0.15^{\circ}$. No significant variability has been found and its gamma-ray spectrum is well fitted by the dark matter annihilation into $b\bar{b}$ with a mass of $\sim 43$ GeV. All these facts make 3FGL J1924.8-1034 a possible dark matter subhalo candidate. However, due to the limited angular resolution, the possibility of that the spatial extension of 3FGL J1924.8-1034 is caused by the contamination from the other un-resolved point source can not be ruled out.

We consider indirect detection of meta-stable dark matter particles decaying into a stable neutral particle and a pair of standard model fermions. Due to the softer energy spectra from the three-body decay, such models could potentially explain the AMS-02 positron excess without being constrained by the Fermi-LAT gamma-ray data and the cosmic ray anti-proton measurements. We scrutinize over different final state fermions, paying special attention to handling of the cosmic ray background and including various contributions from cosmic ray propagation with the help of the \textsc{LikeDM} package. It is found that primary decays into an electron-positron pair and a stable neutral particle could give rise to the AMS-02 positron excess and, at the same time, stay unscathed against the gamma-ray and anti-proton constraints. Decays to a muon pair or a mixed flavor electron-muon pair may also be viable depending on the propagation models. Decays to all other standard model fermions are severely disfavored.

With the large progress in searches for dark matter (DM) particles with indirect and direct methods, we develop a numerical tool that enables fast calculations of the likelihoods of specified DM particle models given a number of observational data, such as charged cosmic rays from space-borne experiments (e.g., PAMELA, AMS-02), gamma-rays from the Fermi space telescope, and underground direct detection experiments. The purpose of this tool --- LikeDM, likelihood calculator for dark matter detection --- is to bridge the gap between a particle model of DM and the observational data. The intermediate steps between these two, including the astrophysical backgrounds, the propagation of charged particles, the analysis of Fermi gamma-ray data, as well as the DM velocity distribution and the nuclear form factor, have been dealt with in the code. We release the first version (v1.0) focusing on the constraints from indirect detection of DM with charged cosmic and gamma rays. Direct detection will be implemented in the next version. This manual describes the framework, usage, and related physics of the code. The code LikeDM can be download from https://likedm.hepforge.org/

The spectral and temporal evolution during X-ray outbursts give important clues on the accretion process and radiation mechanism in black-hole X-ray binaries (BH XRBs). A set of Swift and RXTE observations were executed to monitor the 2008 outburst of the black-hole candidate Swift J1842.5-1124. We investigate these data to explore the accretion physics in BH XRBs. We carry out a comprehensive spectral and timing analysis on all the available pointing observations, including fitting both X-ray spectra and power density spectra, measuring the optical and near-ultraviolet flux density. We also search for correlations among the spectral and timing parameters. The observed properties of Swift J1842.5-1124 are similar to other BH XRBs in many respects, for example the hardness-intensity diagram and hardness-rms diagram. The type-C quasi-periodic oscillations (QPOs) were observed as the source started to transit from the low-hard state to the high-soft state. The frequency of QPOs correlate with intensity and the hard component index, and anti-correlate with the hardness and the total fractional rms. These relations are consistent with the Lense-Thirring precession model. The estimated U-band flux changed with the X-ray flux, while the flux density at the V band remained 0.26 mJy. These results imply that the X-ray reprocessing or the tail of thermal emission from the outer disk contributes a significant fraction of the U-band radiation; alternatively, the companion star or the jet dominates the flux at longer wavelengths.

We report the detection of a GeV $\gamma$-ray source which is spatially overlapping and thus very likely associated with the unidentified very-high-energy (VHE) $\gamma$-ray source HESS J1427-608 with the Pass 8 data recorded by the Fermi Large Area Telescope. The photon spectrum of this source is best described by a power-law with an index of $1.85\pm0.17$ in the energy range of $3-500$ GeV, and the measured flux connects smoothly with that of HESS J1427-608 at a few hundred GeV. This source shows no significant extension and time variation. The broadband GeV-TeV emission over four decades of energies can be well fitted by a single power-law function with an index of 2.0, without obvious indication of spectral cutoff toward high energies. Such a result implies that HESS J1427-608 may be a PeV particle accelerator. We discuss possible nature of HESS J1427-608 according to the multi-wavelength spectral fittings. Given the relatively large errors, either a leptonic or a hadronic model can explain the multi-wavelength data from radio to VHE $\gamma$-rays. The inferred magnetic field strength is a few $\mu$G, which is smaller than typical values of supernova remnants (SNRs), and is consistent with some pulsar wind nebulae (PWNe). On the other hand, the flat $\gamma$-ray spectrum is slightly different from typical PWNe but similar to that of some known SNRs.

Gamma-ray is a good probe of dark matter (DM) particles in the Universe. We search for the DM annihilation signals in the direction of the Andromeda galaxy (M31) using 7.5 year Fermi-LAT pass 8 data. Similar to Pshirkov et al. (2016), we find that there is residual excess emission from the direction of M31 if only the galactic disk as traced by the far infrared emission is considered. Adding a point-like source will improve the fitting effectively, although additional slight improvements can be found if an extended component such as a uniform disk or two bubbles is added instead. Taking the far infrared disk plus a point source as the background model, we search for the DM annihilation signals in the data. We find that there is strong degeneracy between the emission from the galaxy and that from 10s GeV mass DM annihilation in the main halo with quark final state. However, the required DM annihilation cross section is about $10^{-25}-10^{-24}$ cm$^3$s$^{-1}$, orders of magnitude larger than the constraints from observations of dwarf spheroidal galaxies, indicating a non-DM origin of the emission. If DM subhalos are taken into account, the degeneracy is broken. When considering the enhancement from DM subhalos, the constraints on DM model parameters are comparable to (or slightly weaker than) those from the population of dwarf spheroidal galaxies. We also discuss the inverse Compton scattering component from DM annihilation induced electrons/positrons. For the first time we include an energy dependent template of the inverse Compton emission (i.e., a template cube) in the data analysis to take into account the effect of diffusion of charged particles. We find a significant improvement of the constraints in the high mass range of DM particles after considering the inverse Compton emission.

N-body simulations predict that galaxies at the Milky Way scale host a large number of dark matter (DM) subhalos. Some of these subhalos, if they are massive enough or close enough to the Earth, might be detectable in $\gamma$ rays due to the DM annihilation. 3FGL J2212.5+0703, an unidentified gamma-ray source, has been suggested to be the counterpart candidate of a DM subhalo by Bertoni et al. (2015, 2016). In this work we analyze the Fermi-LAT Pass 8 data of 3FGL J2212.5+0703 to independently test the DM subhalo hypothesis of this source. In order to suppress the possible contamination from two nearby very-bright blazars, we just take into account the front-converting gamma-rays which have better angular resolutions than that of the back-converting photons. In addition to the spatial distribution analysis, we have extended the spectrum analysis down to the energies of $\sim 100$ MeV, and thoroughly examined the variability of the emission during the past 8 years. We confirm that 3FGL J2212.5+0703 is a steady and spatially-extended gamma-ray emitter at a high confidence level. The spectrum is well consistent with that expected from DM annihilation into $b\bar{b}$. The introduction of a phenomenological LogParabola spectrum just improves the fit slightly. All these results suggest that 3FGL J2212.5+0703 could be indicative of a DM subhalo.

With 91 months of the publicly available Fermi-LAT Pass 8 data, we analyze the gamma-ray emission from the Milky Way satellites to search for potential line signals due to the annihilation of dark matter particles into double photons. The searched targets include a sample of dwarf spheroidal galaxies, the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). No significant line emission has been found neither in the stacked dwarf galaxy sample nor in the direction of LMC/SMC. The corresponding upper limits on the cross section of DM annihilation into two photons are derived. Compared with results of previous gamma-ray line searches with the Pass 7 data, the current constraints on the line emission from dwarf spheroidal galaxies has been significantly improved in a wide energy range. With the rapid increase of the sample of dwarf spheroidal galaxies (candidates), we expect that the sensitivity of gamma ray line searches will be significantly improved in the near future.

The merger-driven Gamma-ray Bursts (GRBs) and their associated gravitational wave (GW) radiation, if both successfully detected, have some far-reaching implications, including for instance: (i) The statistical comparison of the physical properties of the short/long-short GRBs with and without GW detection can test the general origin model; (ii) Revealing the physical processes taking place at the central engine; (iii) Measuring the velocity of the Gravitational wave directly/accurately. In this work we discuss these implications in the case of possible association of GW150914/ GBM transient 150914. We compared GBM transient 150914 with other SGRBs and found that such an event {may be} a distinct outlier in some statistical diagrams, possibly due to its specific binary-black-hole merger origin. However, the presence of a "new" group of SGRBs with "unusual" physical parameters is also possible. If the outflow of GBM transient 150914 was launched by the accretion onto the nascent black hole, the magnetic activity rather than the neutrino process is likely responsible for the energy extraction and the accretion disk mass is estimated to be $\sim 10^{-5}~M_\odot$. The GW150914/GBM transient 150914 association, {if confirmed, would} provide the first opportunity to directly measure the GW velocity and its departure from the speed of the light {should be within} a factor of $\sim 10^{-17}$.

Relativistic jets can form from at least some tidal disruption events (TDEs) of (sub-)stellar objects around supermassive black holes. We detect the millimeter (MM) emission of IGR J12580+0134 --- the nearest TDE known in the galaxy NGC 4845 at the distance of only 17 Mpc, based on Planck all-sky survey data. The data show significant flux jumps after the event, followed by substantial declines, in all six high frequency Planck bands from 100 GHz to 857 GHz. We further show that the evolution of the MM flux densities are well consistent with our model prediction from an off-axis jet, as was initially suggested from radio and X-ray observations. This detection represents the second TDE with MM detections; the other is Sw J1644+57, an on-axis jetted TDE at redshift of 0.35. Using the on- and off-axis jet models developed for these two TDEs as templates, we estimate the detection potential of similar events with the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/submillimeter Array (ALMA). Assuming an exposure of one hour, we find that the LMT (ALMA) can detect jetted TDEs up to redshifts $z\sim1$ (2), for a typical disrupted star mass of $\sim1$ M$_\odot$. The detection rates of on- and off-axis TDEs can be as high as $\sim0.6$ (13) and 10 (220) per year, respectively, for the LMT (ALMA). We briefly discuss how such observations, together with follow-up radio monitoring, may lead to major advances in understanding the jetted TDEs themselves and the ambient environment of the CNM.

Supernova remnants (SNRs) are believed to be the most important acceleration sites for cosmic rays (CRs) below $\sim10^{15}$ eV in the Galaxy. High energy photons, either directly from the shocks of the SNRs or indirectly from the interaction between SNRs and the nearby clouds, are crucial probes for the CR acceleration. Big progresses on observations of SNRs have been achieved by space- and ground-based $\gamma$-ray facilities. However, whether $\gamma$-rays come from accelerated hadrons or not, as well as their connection with the CRs observed at Earth, remains in debate. Large High Altitude Air Shower Observatory (LHAASO), the next generation experiment, is designed to survey the northern part of the very high energy $\gamma$-ray sky from $\sim 0.3$ TeV to PeV with the sensitivity of $\lesssim1\%$ of the Crab nebula flux. In this paper, we indicate that LHAASO will be dedicated to enlarging the $\gamma$-ray SNR samples and improving the spectral and morphological measurements. These measurements, especially at energies above 30 TeV, will be important for us to finally understand the CR acceleration in SNRs.

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $\gamma$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $\gamma$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

Supernova remnants (SNRs) are the most attractive candidates for the acceleration sites of Galactic cosmic rays. We report the detection of GeV $\gamma$-ray emission with the Pass 8 events recorded by Fermi Large Area Telescope (Fermi-LAT) in the vicinity of the shell type SNR CTB 37B that is likely associated with the TeV $\gamma-$ray source HESS J1713-381. The photon spectrum of CTB 37B is consistent with a power-law with an index of $1.89\pm0.08$ in the energy range of $0.5-500$ GeV, and the measured flux connects smoothly with that of HESS J1713-381 at a few hundred GeV. No significant spatial extension and time variation are detected. The multi-wavelength data can be well fitted with either a leptonic model or a hadronic one. However, parameters of both models suggest more efficient particle acceleration than typical SNRs. Meanwhile, the X-ray and $\gamma$-ray spectral properties of CTB 37B show that it is an interesting source bridging young SNRs dominated by non-thermal emission and old SNRs interacting with molecular clouds.

We have monitored the BL Lacertae object S5 0716+714 simultaneously in the B, R and I bands on three nights in November 2014. The average time resolution is quite high (73s, 34s, 58s for the filters B, R and I), which can help us trace the profile of the variation and search for the short inter-band time delay. Intra-day variability was about 0.1 mag on the first two nights and more than 0.3 mag on the third. A bluer-when-brighter color behavior was found. An clear loop path can be seen on the color-magnitude diagram of the third night, revealing possible time delays between variations at high and low energies. It is the first time that the intra-day spectral hysteresis loop has been found so obviously in the optical band. We used the interpolated cross-correlation function method to further confirm the time delay and calculated the values of lag between light curves at different wavelengths on each night. On the third night, variations in the R and B bands is approximately 1.5 minutes lagging behind the I band. Such optical time delay is probably due to the interplay of different processes of electrons in the jet of the blazar.

Daily X-ray flaring represents an enigmatic phenomenon of Sgr A$^{\star}$ --- the supermassive black hole at the center of our Galaxy. We report initial results from a systematic X-ray study of this phenomenon, based on extensive {\it Chandra} observations obtained from 1999 to 2012, totaling about 4.5 Ms. We detect flares, using a combination of the maximum likelihood and Markov Chain Monte Carlo methods, which allow for a direct accounting for the pile-up effect in the modeling of the flare lightcurves and an optimal use of the data, as well as the measurements of flare parameters, including their uncertainties. A total of 82 flares are detected. About one third of them are relatively faint, which were not detected previously. The observation-to-observation variation of the quiescent emission has an average root-mean-square of $6\%-14\%$, including the Poisson statistical fluctuation of faint flares below our detection limits. We find no significant long-term variation in the quiescent emission and the flare rate over the 14 years. In particular, we see no evidence of changing quiescent emission and flare rate around the pericenter passage of the S2 star around 2002. We show clear evidence of a short-term clustering for the ACIS-S/HETG 0th-order flares on time scale of $20-70$ ks. We further conduct detailed simulations to characterize the detection incompleteness and bias, which is critical to a comprehensive follow-up statistical analysis of flare properties. These studies together will help to establish Sgr A$^{\star}$ as a unique laboratory to understand the astrophysics of prevailing low-luminosity black holes in the Universe.

Supermassive black holes (SMBHs) can capture and tidally disrupt stars or sub-stellar objects orbiting nearby. The detections of Sw J1644+57-like events suggest that at least some TDEs can launch a relativistic jet beaming towards Earth. A natural expectation would be the existence of TDEs with a relativistic jet beaming away from Earth. The nearby TDE candidate IGR J12580+0134 provides new insights into the jet phenomenon. Combining several constraints, we find that the event invokes a $8-40$ Jupiter mass object tidally disrupted by a $3 \times 10^5 - 1.8 \times 10^7 M_\sun$ SMBH. Recently, a bright radio transient was discovered by Irwin et al. in association with IGR J12580+0134. We perform detailed modeling of the event based on a numerical jet model previously developed for the radio emission of Sw J1644+57. We find that the radio data of IGR J12580+0134 can be interpreted within an external forward shock model in the Newtonian regime. Using Sw J1644+57 as a template and properly correcting for its luminosity, we argue that the observed X-ray flux in early times is too faint to allow an on-beam relativistic jet unless the Lorentz factor is very small. Rather, the X-ray emission is likely from the disk or corona near the black hole. From various constraints, we find that the data are consistent with an off-beam relativistic jet with a viewing angle $\theta_{\rm obs} \gtrsim 30^{\rm o}$, and an initial Lorentz factor $\Gamma_j \gtrsim $ a few.This scenario can readily be tested in the upcoming VLBI observations.

X-ray flares have routinely been observed from the supermassive black hole, Sagittarius A$^\star$ (Sgr A$^\star$), at our Galactic center. The nature of these flares remains largely unclear, despite of many theoretical models. In this paper, we study the statistical properties of the Sgr A$^\star$ X-ray flares, by fitting the count rate (CR) distribution and the structure function (SF) of the light curve with a Markov Chain Monte Carlo (MCMC) method. With the 3 million second \textit{Chandra} observations accumulated in the Sgr A$^\star$ X-ray Visionary Project, we construct the theoretical light curves through Monte Carlo simulations. We find that the $2-8$ keV X-ray light curve can be decomposed into a quiescent component with a constant count rate of $\sim6\times10^{-3}~$count s$^{-1}$ and a flare component with a power-law fluence distribution $dN/dE\propto E^{-\alpha_{\rm E}}$ with $\alpha_{\rm E}=1.65\pm0.17$. The duration-fluence correlation can also be modelled as a power-law $T\propto E^{\alpha_{\rm ET}}$ with $\alpha_{\rm ET} < 0.55$ ($95\%$ confidence). These statistical properties are consistent with the theoretical prediction of the self-organized criticality (SOC) system with the spatial dimension $S = 3$. We suggest that the X-ray flares represent plasmoid ejections driven by magnetic reconnection (similar to solar flares) in the accretion flow onto the black hole.

Abell 1767 is a dynamically relaxed, cD cluster of galaxies with a redshift of 0.0703. Among 250 spectroscopically confirmed member galaxies within a projected radius of 2.5r_{200}, 243 galaxies (~ 97%) are spectroscopically covered by the Sloan Digital Sky Survey (SDSS). Based on this homogeneous spectral sample, the stellar evolutionary synthesis code, STARLIGHT, is applied to investigate the stellar populations and star formation histories (SFHs) of cluster galaxies. The star formation properties of galaxies, such as mean stellar ages, metallicities, stellar masses, and star formation rates (SFRs), are presented as the functions of local galaxy density. Strong environmental effect is found in the manner that massive galaxies in the high-density core region of cluster tend to have higher metallicities, longer mean stellar ages, and lower specific star formation rates (SSFRs), and their recent star formation activities have been remarkably suppressed. In addition, the correlations of the metallicity and SSFR with stellar mass are confirmed.

Asymmetric X-ray emission and powerful cluster-scale radio halo indicate that A2319 is a merging cluster of galaxies. This paper presents our multicolor photometry for A2319 with 15 optical intermediate filters in the Beijing-Arizona-Taiwan-Connecticut (BATC) system. There are 142 galaxies with known spectroscopic redshifts within the viewing field, including 128 member galaxies (called sample I).A large velocity dispersion in the rest frame suggests a merger dynamics in A2319. The contour map of projected density and localized velocity structure confirm the so-called A2319B substructure, at ~ 10' NW to the main concentration A2319A. The spectral energy distributions (SEDs) of more than 30,000 sources are obtained in our BATC photometry down to V ~ 20 mag. With color-color diagrams and photometric redshift technique, 233 galaxies brighter than h=19.0 are newly selected as member candidates. The early-type galaxies are found to follow a tight color-magnitude correlation. Based on sample I and the enlarged sample of member galaxies (called sample II), subcluster A2319B is confirmed. A strong environmental effect on star formation histories is found in the manner that galaxies in the sparse regions have various star formation histories, while galaxies in the dense regions are found to have shorter SFR time scales, older stellar ages, and higher ISM metallicities. For the merging cluster A2319, local surface density is a better environmental indicator rather than the clustercentric distance. Compared with the well-relaxed cluster A2589, a higher fraction of star-forming galaxies is found in A2319, indicating that the galaxy-scale turbulence stimulated by the subcluster merger might have played a role in triggering the star formation activity.

The operation of AMS-02 opens a new era for the study of cosmic ray physics with unprecedentedly precise data which are comparable with the laboratory measurements. The high precision data allow a quantitative study on the cosmic ray physics and give strict constraints on the nature of cosmic ray sources. However, the intrinsic errors from the theoretical models to interpret the data become dominant over the errors in the data. In the present work we try to give a systematic study on the uncertainties of the models to explain the AMS-02 positron fraction data, which shows the cosmic ray $e^+e^-$ excesses together with the PAMELA and Fermi-LAT measurements. The excesses can be attributed to contributions from the extra $e^+e^-$ sources, such as pulsars or the dark matter annihilation. The possible systematic uncertainties of the theoretical models considered include the cosmic ray propagation, the treatment of the low energy data, the solar modulation, the $pp$ interaction models, the nuclei injection spectrum and so on. We find that in general a spectral hardening of the primary electron injection spectrum above $\sim50-100$ GeV is favored by the data. Furthermore, the present model uncertainties may lead to a factor of $\sim2$ enlargement in the determination of the parameter regions of the extra source, such as the dark matter mass, annihilation rate and so on.

The AMS-02 has just published the unprecedentedly precise measurement of the cosmic electron and positron spectra. In this paper we try to give a quantitative study on the AMS-02 results by a global fitting to the electron and positron spectra, together with the updated positron fraction data. The Markov Chain Monte Carlo algorithm is adopted to do the fitting. The primary electron spectrum and the parameters for pulsars or dark matter which contribute extra positrons are determined simultaneously. We find that there is a hardening of the primary electron spectrum at $\sim 60$ GeV. With such a new feature at the background spectrum both the pulsars and dark matter can explain the AMS-02 results very well. The dark matter scenario shows a drop at positron fraction at $\sim 300$ GeV, however, suffers very strong constraints from Fermi $\gamma$-ray observations. The fitting results also suggest that the propagation model with convection may be more favored by the lepton data than the reacceleration model.

The recent analyses of the Fermi Large Area Telescope data show an extended GeV $\gamma$-ray excess on top of the expected diffuse background in the Galactic center region, which can be explained with annihilating dark matter or a population of millisecond pulsars (MSPs). We propose to observe the very high energy $\gamma$-rays for distinguishing the MSP scenario from the dark matter scenario. The GeV $\gamma$-ray MSPs should release most energy to the relativistic $e^{\pm}$ wind, which will diffuse in the Galaxy and radiate TeV $\gamma$-rays through inverse Compton scattering and bremsstrahlung processes. By calculating the spectrum and spatial distribution, we show that such emission is detectable with the next generation very high energy $\gamma$-ray observatory, the Cherenkov Telescope Array (CTA), under reasonable model parameters. It is essential to search for the multi-wavelength counterparts to the GeV $\gamma$-ray excess for solving this mystery in the high energy universe.

This paper presents CCD multicolor photometry for 304 old star clusters in the nearby spiral galaxy M31. Of which photometry of 55 star clusters is first obtained. The observations were carried out as a part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey from 1995 February to 2008 March, using 15 intermediate-band filters covering 3000--10000 \AA. Detailed comparisons show that our photometry is in agreement with previous measurements. Based on the ages and metallicities from Caldwell et al. and the photometric measurements here, we estimated the clusters' masses by comparing their multicolor photometry with stellar population synthesis models. The results show that the sample clusters have masses between $\sim 3\times10^4 M_\odot$ and $\sim 10^7 M_\odot$ with the peak of $\sim 4\times10^5 M_\odot$. The masses here are in good agreement with those in previous studies. Combined with the masses of young star clusters of M31 from Wang et al., we find that the peak of mass of old clusters is ten times that of young clusters.

The origin of the knee in cosmic ray spectrum remains to be an unsolved fundamental problem. There are various kinds of models which predict different break positions and the compositions of the knee. In this work, we suggest to use diffuse $\gamma$-rays and neutrinos as probes to test these models. Based on several typical types of the composition models, the diffuse $\gamma$-ray and neutrino spectra are calculated, which show distinctive cutoff behaviours at energies from tens of TeV to multi-PeV. The expected flux will be observable by the newly upgraded Tibet-AS$\gamma$+MD (muon detector) experiment as well as more sensitive future projects, such as LHAASO and HiSCORE. By comparing the neutrino spectrum with the recent observations by IceCube experiment, we find that the diffuse neutrinos from interactions between the cosmic rays and the interstellar medium may not be responsible to the majority of the IceCube events. Future measurements of the neutrinos may be able to identify the Galactic diffuse component and further shed light on the problem of the knee of cosmic rays.

The observations of high energy $\gamma$-ray emission from the Galactic center (GC) by HESS, and recently by Fermi, suggest the cosmic ray acceleration in the GC and possibly around the supermassive black hole. In this work we propose a lepton-hadron hybrid model to explain simultaneously the GeV-TeV $\gamma$-ray emission. Both electrons and hadronic cosmic rays were accelerated during the past activity of the GC. Then these particles would diffuse outwards and interact with the interstellar gas and background radiation field. The collisions between hadronic cosmic rays with gas is responsible to the TeV $\gamma$-ray emission detected by HESS. With fast cooling in the strong radiation field, the electrons would cool down and radiate GeV photons through inverse Compton scattering off the soft background photons. This scenario provides a natural explanation of the observed GeV-TeV spectral shape of $\gamma$-rays.

We present the results of the first systematic long-term multi-color optical monitoring of the BL Lacertae object 1ES 0806+52.4. The monitoring was performed in multiple passbands with a 60/90 cm Schmidt telescope from December 2005 to February 2011. The overall brightness of this object decreased from 2005 December to 2008 December, and regained after that. A sharp outburst probably occurred around the end of our monitoring program. Overlapped on the long-term trend are some short-term small-amplitude oscillations. No intra-night variability was found in the object, which is in accord with the historical observations before 2005. By investigating the color behavior, we found strong bluer-when-brighter chromatism for the long-term variability of 1ES 0806+52.4. The total amplitudes at the c, i and o bands are 1.18, 1.12, and 1.02 mags, respectively. The amplitudes tend to increase toward shorter wavelength, which may be the major cause of bluer-when-brighter. Such bluer-when-brighter is also found in other blazars like S5 0716+714, OJ 287, etc. The hard X-ray data collected from the Swift/BAT archive was correlated with our optical data. No positive result was found, the reason of which may be that the hard X-ray flux is a combination of the synchrotron and inverse Compton emission but with different timescales and cadences under the leptonic Synchrotron-Self-Compton (SSC) model.

It was found in the Fermi-LAT data that there is an extended $\gamma$-ray excess in the Galactic center region. The proposed sources to be responsible for the excess include the dark matter annihilation or an astrophysical alternative from a population of millisecond pulsars (MSPs). Whether or not the MSP scenario can explain the data self-consistently has very important implications for the detection of particle dark matter, which is however, subject to debate in the literature. In this work we study the MSP scenario in detail, based on the detected properties of the MSPs by Fermi-LAT. We build a model of the Milky Way MSPs which can reproduce the $\gamma$-ray properties of the Fermi-LAT MSPs, and derive the intrinsic luminosity function of the MSPs. The model is then applied to a bulge population of MSPs. We find that the extended $\gamma$-ray excess can be well explained by the bulge MSPs without violating the detectable flux distribution of MSPs by Fermi-LAT. The spatial distribution of the bulge MSPs as implied by the distribution of low mass X-ray binaries follows a $r^{-2.4}$ profile, which is also consistent with the $\gamma$-ray excess data. We conclude that the MSP model can explain the Galactic center $\gamma$-ray excess self-consistently, satisfying all the current observational constraints.

Using $5.4$ year Fermi-LAT data, we report the detection of GeV $\gamma$-ray emission from the shell-type supernova remnant RCW 86 (G315.4-2.3) with a significance of $\sim5.1\sigma$. The data slightly favors an extended emission of this supernova remnant. The spectral index of RCW 86 is found to be very hard, $\Gamma\sim1.4$, in the $0.4$ to $300$ GeV range. A one zone leptonic model can well fit the multi-wavelength data from radio to very high energy $\gamma$-rays. The very hard GeV $\gamma$-ray spectrum and the inferred low gas density seem to disfavor the hadronic origin of the $\gamma$-rays. The $\gamma$-ray behavior of RCW 86 is very similar to several other TeV shell-type supernova remnants, e.g., RX J1713.7-3946, RX J0852.0-4622, SN 1006 and HESS J1731-347.

In light of the recent discovery by the ATLAS and CMS experiments at the Large Hadron Collider (LHC) of a Higgs-like particle with a narrow mass range of 125-126 GeV, we perform an updated analysis on one of the popular scalar dark matter models, the Inert Higgs Doublet Model (IHDM). We take into account in our likelihood analysis of various experimental constraints, including recent relic density measurement, dark matter direct and indirect detection constraints as well as the latest collider constraints on the invisible decay width of the Higgs boson and monojet search at the LHC. It is shown that if the invisible decay of the standard model Higgs boson is open, LHC as well as direct detection experiments like LUX and XENON100 could put stringent limits on the Higgs boson couplings to dark matter. We find that the most favoured parameter space for IHDM corresponds to dark matter with a mass less than 100 GeV or so. In particular, the best-fit points are at the dark matter mass around 70 GeV where the invisible Higgs decay to dark matter is closed. Scalar dark matter in the higher mass range of 0.5-4 TeV is also explored in our study. Projected sensitivities for the future experiments of monojet at LHC-14, XENON1T and AMS-02 one year antiproton flux are shown to put further constraints on the existing parameter space of IHDM.

Context: HESS J1731-347 has been identified as one of the few TeV-bright shell-type supernova remnants (SNRs). These remnants are dominated by nonthermal emission, and the nature of TeV emission has been continuously debated for nearly a decade. Aims: We carry out the detailed modeling of the radio to gamma-ray spectrum of HESS J1731-347 to constrain the magnetic field and energetic particles sources, which we compare with those of the other TeV-bright shell-type SNRs explored before. Methods: Four years of data from Fermi Large Area Telescope (LAT) observations for regions around this remnant are analyzed, leading to no detection correlated with the source discovered in the TeV band. The Markov Chain Monte Carlo method is used to constrain parameters of one-zone models for the overall emission spectrum. Results: Based on the 99.9% upper limits of fluxes in the GeV range, one-zone hadronic models with an energetic proton spectral slope greater than 1.8 can be ruled out, which favors a leptonic origin for the gamma-ray emission, making this remnant a sibling of the brightest TeV SNR RX J1713.7-3946, the Vela Junior SNR RX J0852.0-4622, and RCW 86. The best-fit leptonic model has an electron spectral slope of 1.8 and a magnetic field of about 30 muG, which is at least a factor of 2 higher than those of RX J1713.7-3946 and RX J0852.0-4622, posing a challenge to the distance estimate and/or the energy equipartition between energetic electrons and the magnetic field of this source. A measurement of the shock speed will address this challenge and has implications on the magnetic field evolution and electron acceleration driven by shocks of SNRs.

The AMS-02 collaboration has just released its first result of the cosmic positron fraction $e^+/(e^-+e^+)$ with high precision up to $\sim 350$ GeV. The AMS-02 result shows the same trend with the previous PAMELA result, which requires extra electron/positron sources on top of the conventional cosmic ray background, either from astrophysical sources or from dark matter annihilation/decay. In this paper we try to figure out the nature of the extra sources by fitting to the AMS-02 $e^+/(e^-+e^+)$ data, as well as the electron and proton spectra by PAMELA and the $(e^-+e^+)$ spectrum by Fermi and HESS. We adopt the GALPROP package to calculate the propagation of the Galactic cosmic rays and the Markov Chain Monte Carlo sampler to do the fit. We find that the AMS-02 data have implied essential difference from the PAMELA data. There is {\rm tension} between the AMS-02 $e^+/(e^-+e^+)$ data and the Fermi/HESS $(e^-+e^+)$ spectrum, that the AMS-02 data requires less contribution from the extra sources than Fermi/HESS. Then we redo the fit without including the Fermi/HESS data. In this case both the pulsars and dark matter annihilation/decay can explain the AMS-02 data. The pulsar scenario has a soft inject spectrum with the power-law index $\sim 2$, while the dark matter scenario needs $\tau^+\tau^-$ final state with mass $\sim 600$ GeV and a boost factor $\sim 200$.

The recently reported positron fraction up to $\sim 350$ GeV by AMS-02 seems to have tension with the total electron/positron spectra detected by Fermi and HESS, for either pulsar or dark matter annihilation/decay scenario as the primary positron sources. In this work we will show that the tension will be removed by an adjustment of the primary electron spectrum. If the primary electron spectrum becomes harder above $\sim50$ GeV, similar as the cosmic ray nuclei spectrum, the AMS-02 positron fraction and Fermi/HESS data can be well fitted by both the pulsar and dark matter models. This result indicates that there should be a common origin of the cosmic ray nuclei and the primary electrons. Furthermore, this study also implies that the properties of the extra sources derived from the fitting to the AMS-02 data should depend on the form of background.

MGRO J2019+37 within the Cygnus region is a bright and extended source revealed by Milagro at 12-35 TeV. This source is almost as bright as Crab Nebula in northern sky, while it is not confirmed by ARGO-YBJ around TeV. Up to now, no obvious counterpart at low energy wavelengths has been found. Hence, MGRO J2019+37 becomes mysterious and its VHE gamma-ray emission mechanism is attractive. In this paper, a brief summary of the multi-wavelength observations from Radio to gamma-ray is presented. All the available data from the XMM-Newton and INTEGRAL at X-ray, and the Fermi-LAT at gamma-ray bands were used to get constraint on its emission flux at low energy wavelengths. Then, its possible counterparts and the VHE emission mechanism are discussed.

The AMS-02 collaboration has just published a high precision measurement of the cosmic positron fraction $e^+/(e^- + e^+)$, which rises with energy from $\sim 5$ GeV to $\sim 350$ GeV. The result indicates the existence of primary electron/positron sources to account for the positron excess. In this work, we investigate the possibility that the nearby mature pulsars are the primary positron sources. By fitting the data we find that the positrons from a single nearby pulsar, such as Geminga or Monogem, with the spectral index $\alpha \sim 2$ can interpret the AMS-02 result. We also investigate the possibility that high energy positrons are generated by multiple known pulsars in the ATNF catalogue. Such a scenario can also fit the AMS-02 data well. Future precise measurements of fine structures in the positron spectrum would be a support to the pulsar scenario.

Observation of gamma-rays from dwarf galaxies is an effective way to search for particle dark matter. Using 4-year data of Fermi-LAT observations on a series of Milky Way satellites, we develop a general way to search for the signals from dark matter annihilation in such objects. Instead of giving prior information about the energy spectrum of dark matter annihilation, we bin the Fermi-LAT data into several energy bins and build a likelihood map in the "energy bin - flux" plane. The final likelihood of any spectrum can be easily derived through combining the likelihood of all the energy bins. It gives consistent result with that directly calculated using the Fermi Scientific Tool. This method is very efficient for the study of any specific dark matter models with gamma-rays. We use the new likelihood map with Fermi-LAT 4 year data to fit the parameter space in three representative dark matter models: i) toy dark matter model, ii) effective dark matter operators, and iii) supersymmetric neutralino dark matter.

Recently, a tentative 130 GeV $\gamma$-ray line signal was identified by quite a few groups. If correct it would constitute a ``smoking gun'' for dark matter annihilations. Interestingly, the spectra of the cosmic ray electrons detected by PAMELA and Fermi-LAT both show tiny wiggle-like structure at $\sim 100$ GeV, which might indicate a weak signal of the annihilation of $\sim 130$ GeV dark matter particles into electrons/positrons with a velocity-weighted cross section $\langle\sigma v\rangle_{\rm \chi\chi\rightarrow e^{+}e^{-}} \sim 4\times10^{-26}~{\rm cm^{3}~s^{-1}}$. The prospect of identifying such a potential weak dark-matter-annihilation electron and/or positron component by AMS-02, a mission to measure the high energy cosmic ray spectra with unprecedented accuracy, is investigated.

We investigate the electron energy distributions (EEDs) and the acceleration processes in the jet of Mrk 421 through fitting the spectral energy distributions (SEDs) in different active states in the frame of a one-zone synchrotron self-Compton (SSC) model. After assuming two possible EEDs formed in different acceleration models: the shock accelerated power-law with exponential cut-off (PLC) EED and the stochastic turbulence accelerated log-parabolic (LP) EED, we fit the observed SEDs of Mrk 421 in both low and giant flare states by using the Markov Chain Monte Carlo (MCMC) method which constrains the model parameters in a more efficient way. Our calculating results indicate that (1) the PLC and LP models give comparably good fits for the SED in low state, but the variations of model parameters from low state to flaring can be reasonably explained only in the case of the PLC in low state; and (2) the LP model gives better fits compared to the PLC model for the SED in flare state, and the intra-day/night variability observed at GeV-TeV bands can be accommodated only in the LP model. The giant flare may be attributed to the stochastic turbulence re-acceleration of the shock accelerated electrons in low state. Therefore, we may conclude that shock acceleration is dominant in low state, while stochastic turbulence acceleration is dominant in flare state. Moreover, our result shows that the extrapolated TeV spectra from the best-fit SEDs from optical through GeV with the two EEDs are different. It should be considered in caution when such extrapolated TeV spectra are used to constrain extragalactic background light (EBL) models.

We revisit the idea that the Galactic center (GC) is the dominant source of Galactic cosmic rays (GCRs), based on a series of new observational evidence. A unified model is proposed to explain the new phenomena of GCRs and $\gamma$-rays simultaneously. The GCRs are thought to be accelerated during past activities of the GC. The pair production process of GCRs in the strong radiation field due to the GC activity is responsible for the knee structure of the cosmic ray spectra. A f raction of $e^+e^-$ produced by pair production interactions, can be reaccelerated in the induced bipolar jets and be transported into the halo, leaving the Fermi $\gamma$-ray bubbles and WMAP microwave haze as the remnant signal. Finally, the CRs diffuse in the bulge could further interact with the interstellar medium to produce low energy $e^+e^-$. After cooling down, these positrons may annihilate to produce the 511 keV line emission as discovered by INTEGRAL.

The extragalactic background light (EBL) contains important information about stellar and galaxy evolution. It leaves imprint on the very high energy $\gamma$-ray spectra from sources at cosmological distances due to the process of pair production. In this work we propose to {\em measure} the EBL directly by extracting the collective attenuation effects in a number of $\gamma$-ray sources at different redshifts. Using a Markov Chain Monte Carlo fitting method, the EBL intensities and the intrinsic spectral parameters of $\gamma$-ray sources are derived simultaneously. No prior shape of EBL is assumed in the fit. With this method, we can for the first time to derive the spectral shape of the EBL model-independently. Our result shows the expected features predicted by the present EBL models and thus support the understanding of the EBL origin.

We have monitored the BL Lac object S5 0716+714 in five intermediate optical wavebands from 2004 September to 2011 April. Here we present the data that include 8661 measurements which represents one of the largest databases obtained for an object at optical domain. A simple analysis of the data indicates that the object was active in most time, and intraday variability was frequently observed. In total, the object varied by 2.614 magnitudes in the $i$ band. Strong bluer-when-brighter chromatism was observed on long, intermediate, and short timescales.

In this paper, we presented metal abundance properties of 144 M81 globular clusters. These globulars consist of the largest globular cluster sample in M81 till now. Our main results are: the distribution of metallicities are bimodal, with metallicity peaks at [Fe/H]\sim-1.51 and -0.58, and the metal-poor globular clusters tend to be less spatially concentrated than the metal-rich ones; the metal-rich globular clusters in M81 do not demonstrate a centrally concentrated spatial distribution as the metal-rich ones in M31 do; like our Galaxy and M31, the globular clusters in M81 have a small radial metallicity gradient. These results are consistent with those obtained based on a small sample of M81 globular clusters. In addition, this paper showed that there is evidence that a strong rotation of the M81 globular cluster system around the minor axis exists, and that rotation is present in the metal-rich globular cluster subsample, while the metal-poor globular cluster subsample shows no evidence for rotation. The most significant difference between the rotation of the metal-rich and metal-poor globular clusters occurs at intermediate projected galactocentric radii. The results of this paper confirm the conclusion of Schroder et al. that M81's metal-rich globular clusters at intermediate projected radii were associated with a thick disk of M81.

The weakly interacting massive particle (WIMP) often serves as a candidate for the cold dark matter, however when produced non-thermally it could behave like warm dark matter. In this paper we study the properties of the $\gamma$-ray emission from annihilation of WIMP dark matter in the halo of our own Milky-Way Galaxy with high resolution $N$-body simulations of a Milky-Way like dark matter halo, assuming different nature of WIMPs. Due to the large free-streaming length in the scenario of warm WIMPs, the substructure contend of the dark matter halo is significantly different from that of the cold WIMP counterpart, resulting in distinct predictions of the $\gamma$-ray signals from the dark matter annihilation. We illustrate these by comparing the predicted $\gamma$-ray signals from the warm WIMP annihilation to that of cold WIMPs. Pronounced differences from the subhalo skymap and statistical properties between two WIMP models are demonstrated. Due to the potentially enhanced cross section of the non-thermal production mechanism in warm WIMP scenario, the Galactic center might be prior for the indirect detection of warm WIMPs to dwarf galaxies, which might be different from the cold dark matter scenario. As a specific example we consider the non-thermally produced neutralino of supersymmetric model and discuss the detectability of warm WIMPs with Fermi $\gamma$-ray telescope.

It was recently reported that there may exist monochromatic $\gamma$-ray emission at $\sim 130$ GeV from the Galactic center in the Fermi Large Area Telescope data, which might be related with dark matter (DM) annihilation. In this work we carry out a comprehensive check of consistency of the results with the DM annihilation scenario, using the 3.7 yrs Fermi observation of the inner Galaxy, Galactic halo, clusters of galaxies and dwarf galaxies. The results found are as follows. 1) Very strong constraints on the DM annihilation into continuous $\gamma$-rays from the Galactic center are set, which are as stringent as the "natural" scale assuming thermal freeze-out of DM. Such limit sets strong constraint on the DM models to explain the line emission. 2) No line emission from the Galactic halo is found in the Fermi data, and the constraints on line emission is marginally consistent with the DM annihilation interpretation of the $\sim 130$ GeV line emission from the inner Galaxy. 3) No line emission from galaxy clusters and dwarf galaxies is detected, although possible concentration of photons from clusters in 120-140 GeV is revealed. The constraints from clusters and dwarf galaxies are weak and consistent with the DM annihilation scenario to explain the $\sim 130$ GeV line emission.

Shocks of supernova remnants (SNRs) are important (and perhaps the dominant) agents for production of the Galactic cosmic rays. Recent $\gamma$-ray observations of several SNRs have made this case more compelling. However, these broadband high-energy measurements also reveal a variety of spectral shape demanding more comprehensive modeling of emissions from SNRs. According to the locally observed fluxes of cosmic ray protons and electrons, the electron-to-proton number ratio is known to be about 1%. Assuming such a ratio is universal for all SNRs and identical spectral shape for all kinds of accelerated particles, we propose a unified model that ascribes the distinct $\gamma$-ray spectra of different SNRs to variations of the medium density and the spectral difference between cosmic ray electrons and protons observed at Earth to transport effects. For low density environments, the $\gamma$-ray emission is inverse-Compton dominated. For high density environments like systems of high-energy particles interacting with molecular clouds, the $\gamma$-ray emission is $\pi^0$-decay dominated. The model predicts a hadronic origin of $\gamma$-ray emission from very old remnants interacting mostly with molecular clouds and a leptonic origin for intermediate age remnants whose shocks propagate in a low density environment created by their progenitors via e.g., strong stellar winds. These results can be regarded as evidence in support of the SNR-origin of the Galactic cosmic rays.

Recently some hints of the existence of $\gamma$-ray line around 130 GeV are reported according to the analysis of Fermi-LAT data. If confirmed it would be the first direct evidence to show the existence of new physics beyond the standard model. Here we suggest that using the forthcoming high energy resolution $\gamma$-ray detectors, such as CALET and DAMPE, we may test whether it is real line structure or just the background effect. For DAMPE like detector with designed energy resolution $\sim1.5%$, a line significance will reach $11\sigma$ for the same statistics as Fermi-LAT. For about 1.4 yr survey observation, DAMPE may detect a $5\sigma$ signal of such a $\gamma$-ray line.

Recently, several groups identified a tentative $\gamma$-ray line signal with energy $\sim 130$ GeV in the central Galaxy from the Fermi-LAT data. %The morphology study shows that the signal is consistent with dark matter %annihilation, but with an offset $\sim 220$ pc ($1.5^{\circ}$) of the %center from the Galactic center Sgr A$^{\star}$, Such a $\gamma-$ray line can be interpreted as the signal of dark matter annihilation. However, the offset $\sim 220$ pc ($1.5^{\circ}$) of the center of the most prominent signal region from the Galactic center Sgr A$^{\star}$ has been thought to challenge the dark matter annihilation interpretation. Considering the fact that such a 130 GeV $\gamma$-ray line signal consists of only $\sim14$ photons, we suggest that the "imperfect" consistency of these photons with the expected dark matter distribution is due to the limited statistics. The offset will be smaller as more signal photons have been collected in the near future. Our Monte Carlo simulation supports the above speculation.

This paper presents multicolor optical photometry of the nearby galaxy cluster Abell 119 (z = 0:0442) with the Beijing-Arizona-Taiwan-Connecticut (BATC) system of 15 intermediate bands. Within the BATC viewing field of 58'* 58', there are 368 galaxies with known spectroscopic redshifts, including 238 member galaxies (called sample I). Based on the spectral energy distributions (SEDs) of 1376 galaxies brighter than iBATC = 19:5, photometric redshift technique and the color-magnitude relation of earlytype galaxies are applied to select faint member galaxies. As a result, 117 faint galaxies were selected as new member galaxies. Combined with sample I, an enlarged sample (called sample II) of 355 member galaxies is obtained. Spatial distribution and localized velocity structure for two samples demonstrate that A119 is a dynamically complex cluster with at least three prominent substructures in the central region within 1 Mpc. A large velocity dispersion for the central clump indicates a merging along the line of sight. No significant evidences for morphology and luminosity segregations are found in both samples. With the evolutionary synthesis model PEGASE, environmental effect on the star formation properties is confirmed. Faint galaxies in low-density region tend to have longer time scales of star formation, smaller mean stellar ages, and lower metallicities of interstellar medium, which is in agreement with the context of hierarchical cosmological scenario.

The globular clusters are probably good targets for dark matter (DM) searches in $\gamma$-rays due to the possible adiabatic contraction of DM by baryons. In this work we analyse the three-year data collected by {\it Fermi} Large Area Telescope of globular clusters NGC 6388 and M 15 to search for possible DM signals. For NGC 6388 the detection of $\gamma$-ray emission was reported by {\it Fermi} collaboration, which is consistent with the emission of a population of millisecond pulsars. The spectral shape of NGC 6388 is also shown to be consistent with a DM contribution if assuming the annihilation final state is $b\bar{b}$. No significant $\gamma$-ray emission from M 15 is observed. We give the upper limits of DM contribution to $\gamma$-ray emission in both NGC 6388 and M 15, for annihilation final states $b\bar{b}$, $W^+W^-$, $\mu^+\mu^-$, $\tau^+\tau^-$ and monochromatic line. The constraints are stronger than that derived from observation of dwarf galaxies by {\it Fermi}.

A self-consistent global fitting method based on the Markov Chain Monte Carlo technique to study the dark matter (DM) property associated with the cosmic ray electron/positron excesses was developed in our previous work. In this work we further improve the previous study to include the hadronic branching ratio of DM annihilation/decay. The PAMELA $\bar{p}/p$ data are employed to constrain the hadronic branching ratio. We find that the 95% ($2\sigma$) upper limits of the quark branching ratio allowed by the PAMELA $\bar{p}/p$ data is $\sim 0.032$ for DM annihilation and $\sim 0.044$ for DM decay respectively. This result shows that the DM coupling to pure leptons is indeed favored by the current data. Based on the global fitting results, we further study the neutrino emission from DM in the Galactic center. Our predicted neutrino flux is some smaller than previous works since the constraint from $\gamma$-rays is involved. However, it is still capable to be detected by the forth-coming neutrino detector such as IceCube. The improved points of the present study compared with previous works include: 1) the DM parameters, both the particle physical ones and astrophysical ones, are derived in a global fitting way, 2) constraints from various species of data sets, including $\gamma$-rays and antiprotons are included, and 3) the expectation of neutrino emission is fully self-consistent.

Recently PAMELA collaboration published the cosmic nuclei and electron spectra with high precision, together with the cosmic antiproton data updated, and the Fermi-LAT collaboration also updated the measurement of the total $e^+e^-$ spectrum to lower energies. In this paper we develop a Markov Chain Monte Carlo (MCMC) package {\it CosRayMC}, based on the GALPROP cosmic ray propagation model to study the implications of these new data. It is found that if only the background electrons and secondary positrons are considered, the fit is very bad with $\chi_{\rm red}^2 \approx 3.68$. Taking into account the extra $e^+e^-$ sources of pulsars or dark matter annihilation we can give much better fit to these data, with the minimum $\chi_{\rm red}^2 \approx 0.83$. This means the extra sources are necessary with a very high significance in order to fit the data. However, the data show little difference between pulsar and dark matter scenarios. Both the background and extra source parameters are well constrained with this MCMC method. Including the antiproton data, we further constrain the branching ratio of dark matter annihilation into quarks $B_q<0.5%$ at $2\sigma$ confidence level. The possible systematical uncertainties of the present study are discussed.

The superluminal neutrinos detected by OPERA indicates Lorentz invariance violation (LIV) of the neutrino sector at the order of $10^{-5}$. We study the implications of the result in this work. We find that such a large LIV implied by OPERA data will make the neutrino production process $\pi \to \mu + \nu_\mu$ kinematically forbidden for neutrino energy greater than about 5 GeV. The OPERA detection of neutrinos at 40 GeV can constrain the LIV parameter to be smaller than $3\times 10^{-7}$. Furthermore the neutrino decay in the LIV framework will modify the neutrino spectrum greatly. The atmospheric neutrino spectrum measured by IceCube can constrain the LIV parameter to the level of $10^{-12}$. The future detection of astrophysical neutrinos of Galactic sources is expected to be able to give even stronger constraint on the LIV parameter of neutrinos.

In this paper we present a photometric and spectroscopic study of the nearby galaxy cluster Abell 0671 (A671) with 15 intermediate-band filters in the Beijing-Arizona-Taiwan-Connecticut (BATC) system and the Sloan Digital Sky Survey (SDSS) data. The photometric redshift technique is applied to the galaxy sample for further membership determination. After the color-magnitude relation is taken into account, 97 galaxies brighter than h_batc=19.5 mag are selected as new member galaxies. Based on the enlarged sample of cluster galaxies, spatial distribution, dynamics of A671 are investigated. The substructures of A671 are well shown by the sample of bright members, but it appears less significant based on the enlarged sample, which is mainly due to larger uncertainties in the light-of-sight velocities of the newly-selected faint members. The SDSS r-band luminosity function of A671 is flat at faint magnitudes, with the faint end slope parameter alpha=-1.12. The SDSS spectra allow us to investigate the star formation history of bright cluster galaxies, and the galaxies in the core region are found to be older than those in the outskirts. No environmental effect is found for metallicities of the early-type galaxies (ETGs). The positive correlation between age and stellar mass supports the downsizing scenario. By comparing ETG absorption-line indices with the state-of-art stellar population models, we derive the relevant parameters of simple stellar population (such as age, [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe]). The ETGs at cluster center tend to have smaller hb indices, indicating that central ETGs are likely to be older. The distribution of total metallicity indicator, [MgFe]', does not show any environmental effects. The relations between the simple stellar population parameters and velocity dispersion in A671 are in good agreement with previous studies.

The cosmic ray (CR) energy spectra measured with ATIC, CREAM and PAMELA showed that there is remarkable hardening for rigidity of several hundred GV. We propose that this hardening is due to the superposition of spectra from a population of sources, e.g., supernova remnants (SNRs), whose injection spectral indices have a dispersion. Adopting proper model parameters the observational data can be well explained. It is interesting that the injection source parameters are similar with that derived from gamma-ray observations of SNRs, which may support the SNR-origin of CRs. Furthermore this mechanism provides an alternative explanation of the "ankle-cutoff" structure of the ultra high energy CR spectra.

Recently the AGILE and Fermi/LAT detectors uncovered giant $\gamma$-ray flares from the Crab nebula. The duration of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, we propose that the $\gamma$-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to $\sim$PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a picture can reproduce the observational data quite well given proper model parameters.

Recent cosmic ray (CR) experiments discovered that the CR spectra experience a remarkable hardening for rigidity above several hundred GV. We propose that this is caused by the superposition of the CR energy spectra of many sources that have a dispersion in the injection spectral indices. Adopting similar parameters as those of supernova remnants derived from the Fermi $\gamma$-ray observations, we can reproduce the observational CR spectra of different species well. This may be interpreted as evidence to support the supernova remnant origin of CRs below the knee. We further propose that the same mechanism may explain the "ankle" of the ultra high energy CR spectrum.

In the paper we study the Breit-Wigner enhancement of dark matter (DM) annihilation considering the kinetic decoupling in the evolution of DM freeze-out at the early universe. Since the DM temperature decreases much faster (as $1/R^2$) after kinetic decoupling than that in kinetic equilibrium (as 1/R) we find the Breit-Wigner enhancement of DM annihilation rate after the kinetic decoupling will affect the DM relic density significantly. Focusing on the model parameters that trying to explain the anomalous cosmic positron/electron excesses observed by PAMELA/Fermi/ATIC we find the elastic scattering $Xf\to Xf$ is not efficient to keep dark matter in kinetic equilibrium, and the kinetic decoupling temperature $T_{kd}$ is comparable to the chemical decoupling temperature $T_f\sim O(10) GeV$. The reduction of the relic density after $T_{kd}$ is significant and leads to a limited enhancement factor $\sim O(10^2)$. Therefore it is difficult to explain the anomalous positron/electron excesses in cosmic rays by DM annihilation and give the correct DM relic density simultaneously in the minimal Breit-Wigner enhancement model.

An optical photometric observation in 15 bands was carried out for nearby galaxy cluster Abell 1775 by the Beijing-Arizona-Taiwan-Connecticut (BATC) multi-color system. Over 5000 sources' spectral energy distributions (SEDs) were obtained. Since this cluster has also been observed by the Sloan Digital Sky Survey (SDSS), the BATC SEDs were combined with the SDSS five-band photometric data. Using the combined SEDs, 146 faint galaxies were selected as new member galaxies by the photometric redshift technique. Based on the positions, redshifts and 20-band SEDs of member galaxies, dynamical substructures and luminosity functions (LFs) of A1775 were investigated. The previous reported bimodal structure of A1775 has been confirmed: a poor subcluster with lower redshift, A1775A, is located ~14' southeast to the main concentration A1775B. After taking into account the new supplemented member galaxies, a new subcluster A1775C was found along the aligned direction of A1775A and A1775B. The different LF faint ends of the two subclusters indicate that A1775B is a more dynamically evolved system, while A1775A is still dynamically young. By the STARLIGHT spectral synthesis code, the star- formation histories of the member galaxies were studied. The dependence of the mean stellar ages upon the Hubble type was confirmed, and the environmental effect on star-formation activities for galaxies in A1775B has been explored.

Emission mechanisms of the shell-type supernova remnant (SNR) RX J1713.7-3946 are studied with multi-wavelength observational data from radio, X-ray, GeV $\gamma$-ray to TeV $\gamma$-ray band. A Markov Chain Monte Carlo method is employed to explore the high-dimensional model parameter space systematically. Three scenarios for the $\gamma$-ray emission are investigated: the leptonic, the hadronic and a hybrid one. Thermal emission from the background plasma is also included to constrain the gas density, assuming ionization equilibrium, and a 2$\sigma$ upper limit of about 0.03 cm$^{-3}$ is obtained as far as thermal energies account for a significant fraction of the dissipated kinetic energy of the SNR shock. Although systematic errors dominate the $\chi^2$ of the spectral fit of all models, we find that 1) the leptonic model has the best constrained model parameters, whose values can be easily accommodated with a typical supernova, but gives relatively poor fit to the TeV $\gamma$-ray data; 2) The hybrid scenario has one more parameter than the leptonic one and improves the overall spectral fit significantly; 3) The hadronic one, which has three more parameters than the leptonic model, gives the best fit to the overall spectrum with relatively not-well-constrained model parameters and very hard spectra of accelerated particles. The uncertainties of the model parameters decrease significantly if the spectral indices of accelerated electrons and protons are the same. The hybrid and hadronic models also require an energy input into high-energy protons, which seems to be too high compared with typical values of a supernova explosion. Further investigations are required to reconcile these observations with SNR theories.

It has been proposed that during the formation of the first generation stars there might be a "dark star" phase in which the power of the star comes from dark matter annihilation. The adiabatic contraction process to form the dark star would result in a highly concentrated density profile of the host halo at the same time, which may give enhanced indirect detection signals of dark matter. In this work we investigate the extragalactic $\gamma$-ray background from dark matter annihilation with such a dark star formation scenario, and employ the isotropic $\gamma$-ray data from Fermi-LAT to constrain the model parameters of dark matter. The results suffer from large uncertainties of both the formation rate of the first generation stars and the subsequent evolution effects of the host halos of the dark stars. We find, in the most optimistic case for $\gamma$-ray production via dark matter annihilation, the expected extragalactic $\gamma$-ray flux will be enhanced by 1-2 orders of magnitude. In such a case, the annihilation cross section of the supersymmetric dark matter can be constrained to the thermal production level, and the leptonic dark matter model which is proposed to explain the positron/electron excesses can be well excluded. Conversely, if the positron/electron excesses are of a dark matter annihilation origin, then the early Universe environment is such that no dark star can form.

The Fermi $\gamma$-ray space telescope reported the observation of several Galactic supernova remnants recently, with the $\gamma$-ray spectra well described by hadronic $pp$ collisions. The possible neutrino emissions from these Fermi detected supernova remnants are discussed in this work, assuming the hadronic origin of the $\gamma$-ray emission. The muon event rates induced by the neutrinos from these supernova remnants on typical km$^3$ neutrino telescopes, such as the IceCube and the KM3NeT, are calculated. The results show that for most of these supernova remnants the neutrino signals are too weak to be detected by the on-going or up-coming neutrino experiment. Only for the TeV bright sources RX J1713.7-3946 and possibly W28 the neutrino signals can be comparable with the atmospheric background in the TeV region, if the protons can be accelerated to very high energies. The northern hemisphere based neutrino telescope might detect the neutrinos from these two sources.

A statistical scenario is proposed to explain the $\gamma$-ray variability and flares of the Crab nebula, which were observed recently by the Fermi/LAT. In this scenario electrons are accelerated in a series of knots, whose sizes follow a power-law distribution. These knots presumably move outwards from the pulsar and have a distribution in the Doppler boost factor. The maximal electron energy is assumed to be proportional to the size of the knot. Fluctuations at the highest energy end of the overall electron distribution will result in variable $\gamma$-ray emission via the synchrotron process in the $\sim 100$ MeV range. Since highly boosted larger knots are rarer than smaller knots, the model predicts that the variability of the synchrotron emission increases with the photon energy. We realize such a scenario with a Monte-Carlo simulation and find that the model can reproduce both the two $\gamma$-ray flares over a period of $\sim$ year and the monthly scale $\gamma$-ray flux fluctuations as observed by the Fermi/LAT. The observed $\gamma$-ray spectra in both the steady and flaring states are also well reproduced.

In this paper we present a multicolor photometry for A2589 ($z=0.0414$) with 15 intermediate bands in the Beijing-Arizona-Taiwan-Connecticut (BATC) system which covers an optical wavelength range from 3000 \AA\ to 10000 \AA. The spectral energy distributions (SEDs) for more than 5000 sources are achieved down to {\it V} $\sim$ 20 mag in about 1 deg$^{2}$ field. A2589 has been also covered by the Sloan Digital Sky Survey (SDSS) in photometric mode only. A cross-identification of the BATC-detected galaxies with the SDSS photometric catalog achieves 1199 galaxies brighter than $i=19.5$ mag, among which 68 member galaxies with known spectroscopic redshifts are found. After combining the SDSS five-band photometric data and the BATC SEDs, the technique of photometric redshift is applied to these galaxies for selecting faint member galaxies. The color-magnitude relation is taken as a further restriction of early-type cluster galaxies. As a result, 106 galaxies are newly selected as member galaxies. Spatial distribution of member galaxies shows a north-south elongation which agrees with the X-ray brightness profile and the orientation of central cD galaxy, NGC 7647. No substructures are detected on the basis of positions and radial velocities of cluster galaxies, indicating that A2589 is a well-relaxed system. The luminosity function of A2589 exhibits a peak at $M_{R} \sim -20$ mag and a dip at $M_{R} \sim -19 $ mag. The low-density outer regions are the preferred habitat of faint galaxies. With the evolutionary population synthesis model, PEGASE, the environmental effect on the star formation properties for 68 spectroscopically confirmed member galaxies is studied. The outlier faint galaxies tend to have longer time scales of star formation, shorter mean stellar ages, and lower metallicities of interstellar medium, which can be interpreted in the context of hierarchical cosmological scenario.

We calculate the CMB temperature distortion due to the energetic electrons and positrons produced by dark matter annihilation (Sunyaev-Zel'dovich effect), in dwarf spheroidal galaxies (dSphs). In the calculation we have included two important effects which were previously ignored. First we show that the electron-positron pairs with energy less than GeV, which were neglected in previous calculation, could contribute a significant fraction of the total signal. Secondly we also consider the full effects of diffusion loss, which could significantly reduce the density of electron-positron pairs at the center of cuspy halos. For neutralinos, we confirm that detecting such kind of SZ effect is beyond the capability of the current or even the next generation experiments. In the case of light dark matter (LDM) the signal is much larger, but even in this case it is only marginally detectable with the next generation of experiment such as ALMA. We conclude that similar to the case of galaxy clusters, in the dwarf galaxies the $SZ_2DM}$ effect is not a strong probe of DM annihilations.

The cosmic electron and positron excesses have been explained as possible dark matter (DM) annihilation products. In this work we investigate the possible effects of such a DM annihilation scenario during the evolution history of the Universe. We first calculate the extragalactic $\gamma$-ray background (EGRB), which is produced through the final state radiation of DM annihilation to charged leptons and the inverse Compton scattering between electrons/positrons and the cosmic microwave background. The DM halo profile and the minimal halo mass, which are not yet well determined from the current N-body simulations, are constrained by the EGRB data from EGRET and Fermi telescopes. Then we discuss the impact of such leptonic DM models on cosmic evolution, such as the reionization and heating of intergalactic medium, neutral Hydrogen 21 cm signal and suppression of structure formation. We show that the impact on the Hydrogen 21 cm signal might show interesting signatures of DM annihilation, but the influence on star formation is not remarkable. Future observations of the 21 cm signals could be used to place new constraints on the properties of DM.

Monochromatic gamma-rays are thought to be the smoking gun signal for identifying the dark matter annihilation. However, the flux of monochromatic gamma-rays is usually suppressed by the virtual quantum effects since dark matter should be neutral and does not couple with gamma-rays directly. In the work we study the detection strategy of the monochromatic gamma-rays in a future space-based detector. The monochromatic gamma-ray flux is calculated by assuming supersymmetric neutralino as a typical dark matter candidate. We discuss both the detection focusing on the Galactic center and in a scan mode which detects gamma-rays from the whole Galactic halo are compared. The detector performance for the purpose of monochromatic gamma-rays detection, with different energy and angular resolution, field of view, background rejection efficiencies, is carefully studied with both analytical and fast Monte-Carlo method.

(Aims.) SNR RX J1713.7-3946 is perhaps one of the best observed shell-type supernova remnants with emissions dominated by energetic particles accelerated near the shock front. The nature of the TeV emission, however, is an issue still open to investigation. (Methods.) We carry out a systematic study of four lepton models for the TeV emission with the Markov chain Monte Carlo method. (Results.) It is shown that current data already give good constraints on the model parameters. Two commonly used parametric models do not appear to fit the observed radio, X-ray, and gamma-ray spectra. Models motivated by diffusive shock acceleration and by stochastic acceleration by compressive waves in the shock downstream give comparably good fits. The former has a sharper spectral cutoff in the hard X-ray band than the latter. Future observations with the HXMT and NuSTAR may distinguish these two models.

In this work, we revisit the all-sky Galactic diffuse $\gamma$-ray emission taking into account the new measurements of cosmic ray electron/positron spectrum by PAMELA, ATIC and Fermi, which show excesses of cosmic electrons/positrons beyond the expected fluxes in the conventional model. Since the origins of the extra electrons/positrons are not clear, we consider three different scenarios to account for the excesses: the astrophysical sources such as the Galactic pulsars, dark matter decay and annihilation. Further, new results from Fermi-LAT of the (extra-)Galactic diffuse $\gamma$-ray are adopted. The background cosmic rays without the new sources give lower diffuse $\gamma$ rays compared to Fermi-LAT observation, which is consistent with previous analysis. The scenario with astrophysical sources predicts diffuse $\gamma$-rays with little difference with the background. The dark matter annihilation models with $\tau^{\pm}$ final state are disfavored by the Fermi diffuse $\gamma$-ray data, while there are only few constraints on the decaying dark matter scenario. Furthermore, these is always a bump at higher energies ($\sim$ TeV) of the diffuse $\gamma$-ray spectra for the dark matter scenarios due to final state radiation. Finally we find that the Fermi-LAT diffuse $\gamma$-ray data can be explained by simply enlarging the normalization of the electron spectrum without introduce any new sources, which may indicate that the current constraints on the dark matter models can be much stronger given a precise background estimate.

The $\gamma$-ray and neutrino emissions from dark matter (DM) annihilation in galaxy clusters are studied. After about one year operation of Fermi-LAT, several nearby clusters are reported with stringent upper limits of GeV $\gamma$-ray emission. We use the Fermi-LAT upper limits of these clusters to constrain the DM model parameters. We find that the DM model distributed with substructures predicted in cold DM (CDM) scenario is strongly constrained by Fermi-LAT $\gamma$-ray data. Especially for the leptonic annihilation scenario which may account for the $e^{\pm}$ excesses discovered by PAMELA/Fermi-LAT/HESS, the constraint on the minimum mass of substructures is of the level $10^2-10^3$ M$_{\odot}$, which is much larger than that expected in CDM picture, but is consistent with a warm DM scenario. We further investigate the sensitivity of neutrino detections of the clusters by IceCube. It is found that neutrino detection is much more difficult than $\gamma$-rays. Only for very heavy DM ($\sim 10$ TeV) together with a considerable branching ratio to line neutrinos the neutrino sensitivity is comparable with that of $\gamma$-rays.

It is found that \feii emission contributes significantly to the optical and ultraviolet spectra of most active galactic nuclei. The origin of the optical/UV \feii emission is still a question open to debate. The variability of \feii would give clues to this origin. Using 7.5 yr spectroscopic monitoring data of one Palomer-Green (PG) quasi-stellar object (QSO), PG 1700+518, with strong optical \feii emission, we obtain the light curves of the continuum \lv, \feii, the broad component of \hb, and the narrow component of \hb by the spectral decomposition. Through the interpolation cross-correlation method, we calculate the time lags for light curves of \feii, the total \hb, the broad component of \hb, and the narrow component of \hb with respect to the continuum light curve. We find that the \feii time lag in PG1700+518 is $209^{+100}_{-147}$ days, and the \hb time lag cannot be determined. Assuming that \feii and \hb emission regions follow the virial relation between the time lag and the FWHM for the \hb and \feii emission lines, we can derive that the \hb time lag is $148^{+72}_{-104}$ days. The \hb time lag calculated from the empirical luminosity--size relation is 222 days, which is consistent with our measured \feii time lag. Considering the optical \feii contribution, PG 1700+518 shares the same characteristic on the spectral slope variability as other 15 PG QSOs in our previous work, i.e., harder spectrum during brighter phase.

The paper investigates the overall and detailed features of cosmic ray (CR) spectra in the knee region using the scenario of nuclei-photon interactions around the acceleration sources. Young supernova remnants can be the physical realities of such kind of CR acceleration sites. The results show that the model can well explain the following problems simultaneously with one set of source parameters: the knee of CR spectra and the sharpness of the knee, the detailed irregular structures of CR spectra, the so-called "component B" of Galactic CRs, and the electron/positron excesses reported by recent observations. The coherent explanation serves as evidence that at least a portion of CRs might be accelerated at the sources similar to young supernova remnants, and one set of source parameters indicates that this portion mainly comes from standard sources or from a single source.

In this paper we develop a Markov Chain Monte Carlo code to study the dark matter properties in interpreting the recent observations of cosmic ray electron/positron excesses. We assume that the dark matter particles couple dominantly to leptons and consider two cases separately with annihilating or decaying into lepton pairs. The constraint on the central density profile from diffuse $\gamma$-rays around the Galactic center is also included in the Markov Chain Monte Carlo code self-consistently. In the numerical study, 7 parameters are introduced. We fit two data sets independently and find that for the Data set I (PAMELA+ATIC), dark matter with $m_{\chi}\approx0.7$ TeV for annihilation (or 1.4 TeV for decay) and pure $e^+e^-$ channel is favored, while for the Data set II (PAMELA+Fermi-LAT+H.E.S.S.) $m_{\chi}\approx 2$ TeV for annihilation (or 4 TeV for decay) and the combination of $\mu^+\mu^-$ and $\tau^+\tau^-$ final states can best fit the data. The H.E.S.S. observation of the Galactic center $\gamma$-rays puts a strong constraint on the central density profile of the dark matter halo for the annihilation dark matter scenario. In this case the NFW profile which is regarded as the typical predication from the cold dark matter scenario, is excluded with a high significance ($>5\sigma$). For the decaying dark matter scenario, the constraint is much weaker.

Boost factors of dark matter annihilation into antiprotons and electrons/positrons due to the clumpiness of dark matter distribution are studied in detail in this work, taking the Sommerfeld effect into account. It has been thought that the Sommerfeld effect, if exists, will be more remarkable in substructures because they are colder than the host halo, and may result in a larger boost factor. We give a full calculation of the boost factors based on the recent N-body simulations. Three typical cases of Sommerfeld effects, the non-resonant, moderately resonant and strongly resonant cases are considered. We find that for the non-resonant and moderately resonant cases the enhancement effects of substructures due to the Sommerfeld effect are very small ($\sim \mathcal{O}(1)$) because of the saturation behavior of the Sommerfeld effect. For the strongly resonant case the boost factor is typically smaller than $\sim \mathcal{O}(10)$. However, it is possible in some very extreme cases that DM distribution is adopted to give the maximal annihilation the boost factor can reach up to $\sim 1000$. The variances of the boost factors due to different realizations of substructures distribution are also discussed in the work.

In this work we study how the cosmological parameter, the Hubble constant $H_0$, can be constrained by observation of very high energy (VHE) $\gamma$-rays at the TeV scale. The VHE $\gamma$-rays experience attenuation by background radiation field through $e^+e^-$ pair production during the propagation in the intergalactic space. This effect is proportional to the distance that the VHE $\gamma$-rays go through. Therefore the absorption of TeV $\gamma$-rays can be taken as cosmological distance indicator to constrain the cosmological parameters. Two blazars Mrk 501 and 1ES 1101-232, which have relatively good spectra measurements by the atmospheric Cerenkov telescope, are studied to constrain $H_0$. The mechanism constraining the Hubble constant adopted here is very different from the previous methods such as the observations of type Ia supernovae and the cosmic microwave background. However, at $2\sigma$ level, our result is consistent with other methods.

In this work we calculate the Sunyaev-Zel'dovich (SZ) effect due to the $e^+e^-$ from dark matter (DM) annihilation in galaxy clusters. Two candidates of DM particle, (1) the weakly-interacting massive particle (WIMP) and (2) the light dark matter (LDM) are investigated. For each case, we also consider several DM profiles with and without central cusp. We generally find smaller signals than previously reported. Moreover, the diffusion of electrons and positrons in the galaxy clusters, which was generally thought to be negligible, is considered and found to have significant effect on the central electron/positron distribution for DM profile with large spatial gradient. We find that the SZ effect from WIMP is almost always non-observable, even for the highly cuspy DM profile, and using the next generation SZ interferometer such as ALMA. Although the signal of the LDM is much larger than that of the WIMP, the final SZ effect is still very small due to the smoothing effect of diffusion. Only for the configuration with large central cusp and extremely small diffusion effect, the LDM induced SZ effect might have a bit chance of being detected.

An optical photometric observation with the Beijing-Arizona-Taiwan-Connecticut (BATC) multicolor system is carried out for A98 (z=0.104), a galaxy cluster with two large enhancements in X-ray surface brightness. The spectral energy distributions (SEDs) covering 15 intermediate bands are obtained for all sources detected down to V ~ 20 mag in a field of $58' \times 58'$. After the star-galaxy separation by the color-color diagrams, a photometric redshift technique is applied to the galaxy sample for further membership determination. The color-magnitude relation is taken as a further restriction of the early-type cluster galaxies. As a result, a list of 198 faint member galaxies is achieved. Based on newly generated sample of member galaxies, the dynamical substructures, A98N, A98S, and A98W, are investigated in detail. A separate galaxy group, A98X, is also found to the south of main concentration of A98, which is gravitationally unbound to A98. For 74 spectroscopically confirmed member galaxies, the environmental effect on the star formation histories is found. The bright galaxies in the core region are found to have shorter time scales of star formation, longer mean stellar ages, and higher metallicities of interstellar medium, which can be interpreted in the context of hierarchical cosmological scenario.

The excesses of the cosmic positron fraction recently measured by PAMELA and the electron spectra by ATIC, PPB-BETS, Fermi and H.E.S.S. indicate the existence of primary electron and positron sources. The possible explanations include dark matter annihilation, decay, and astrophysical origin, like pulsars. In this work we show that these three scenarios can all explain the experimental results of the cosmic $e^\pm$ excess. However, it may be difficult to discriminate these different scenarios by the local measurements of electrons and positrons. We propose possible discriminations among these scenarios through the synchrotron and inverse Compton radiation of the primary electrons/positrons from the region close to the Galactic center. Taking typical configurations, we find the three scenarios predict quite different spectra and skymaps of the synchrotron and inverse Compton radiation, though there are relatively large uncertainties. The most prominent differences come from the energy band $10^4\sim 10^9$ MHz for synchrotron emission and $\gtrsim 10$ GeV for inverse Compton emission. It might be able to discriminate at least the annihilating dark matter scenario from the other two given the high precision synchrotron and diffuse $\gamma$-ray skymaps in the future.

Supernova remnants have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants.

The perspective of the detectability of Galactic dark matter subhaloes on the Fermi satellite is investigated in this work. Under the assumptions that dark matter annihilation accounts for the "GeV excess" of the Galactic diffuse $\gamma$-rays discovered by EGRET and the $\gamma$-ray flux is dominated by the contribution from subhaloes of dark matter, we calculate the expected number of dark matter subhaloes that Fermi may detect. We show that Fermi may detect a few tens to several hundred subhaloes in 1-year all sky survey. Since EGRET observation is taken as a normalization, this prediction is independent of the particle physics property of dark matter. The uncertainties of the prediction are discussed in detail. We find that the major uncertainty comes from the mass function of subhaloes, i.e., whether the subhaloes are "point like" (high-mass rich) or "diffuse like" (low-mass rich). Other uncertainties like the background estimation and the observational errors will contribute a factor of $2\sim 3$.

In this paper we propose a dark matter model and study aspects of its phenomenology. Our model is based on a new dark matter sector with a U(1)' gauge symmetry plus a discrete symmetry added to the Standard Model of particle physics. The new fields of the dark matter sector have no hadronic charges and couple only to leptons. Our model can not only give rise to the observed neutrino mass hierarchy, but can also generate the baryon number asymmetry via non-thermal leptogenesis. The breaking of the new U(1)' symmetry produces cosmic strings. The dark matter particles are produced non-thermally from cosmic string loop decay which allows one to obtain sufficiently large annihilation cross sections to explain the observed cosmic ray positron and electron fluxes recently measured by the PAMELA, ATIC, PPB-BETS, Fermi-LAT, and HESS experiments while maintaining the required overall dark matter energy density. The high velocity of the dark matter particles from cosmic string loop decay leads to a low phase space density and thus to a dark matter profile with a constant density core in contrast to what happens in a scenario with thermally produced cold dark matter where the density keeps rising towards the center. As a result, the flux of gamma rays radiated from the final leptonic states of dark matter annihilation from the Galactic center is suppressed and satisfies the constraints from the HESS gamma-ray observations.

The GZK cutoff predicted at the Ultra High Energy Cosmic Ray (UHECR) spectrum as been observed by the HiRes and Auger experiments. The results put severe constraints on the effect of Lorentz Invariance Violation(LIV) which has been introduced to explain the absence of GZK cutoff indicated in the AGASA data. Assuming homogeneous source distribution with a single power law spectrum, we calculate the spectrum of UHECRs observed on Earth by taking the processes of photopion production, $e^+e^-$ pair production and adiabatic energy loss into account. The effect of LIV is also taken into account in the calculation. By fitting the HiRes monocular spectra and the Auger combined spectra, we show that the LIV parameter is constrained to $\xi=-0.8^{+3.2}_{-0.5}\times10^{-23}$ and $0.0^{+1.0}_{-0.4}\times10^{-23}$ respectively, which is well consistent with strict Lorentz Invariance up to the highest energy.

Recently the Milagro experiment observed diffuse multi-TeV gamma-ray emission in the Cygnus region, which is significantly stronger than what predicted by the Galactic cosmic ray model. However, the sub-GeV observation by EGRET shows no excess to the prediction based on the same model. This TeV excess implies possible high energy cosmic rays populated in the region with harder spectrum than that observed on the Earth. In the work we studied this theoretical speculation in detail. We find that, a diffuse proton source with power index $\alpha_p\lesssim 2.3$, or a diffuse electron source with power index $\alpha_e\lesssim2.6$ can reproduce the Milagro's observation without conflicting with the EGRET data. Further detections on neutrinos, a diagnostic of the hadronic model, and hard X-ray synchrontron radiation, a diagnostic of the lepton model, help to break this degeneracy. In combination with the gamma ray observations to several hundred GeV by Fermi, we will be able to understand the diffuse emission mechanisms in the Cygnus region better.

We employ an analytical approach to investigate the signatures of Baryon Acoustic Oscillations (BAOs) on the convergence power spectrum of weak lensing by large scale structure. It is shown that the BAOs wiggles can be found in both of the linear and nonlinear convergence power spectra of weak lensing at about $40\le l\le600$, but they are weaker than that of matter power spectrum. Although the statistical error for LSST are greatly smaller than that of CFHT and SNAP survey especially at about $30<l<300$, they are still larger than the their maximum variations of BAOs wiggles. Thus, the detection of BAOs with the ongoing and upcoming surveys such as LSST, CFHT and SNAP survey confront a technical challenge.

We use observations of Sunyaev-Zel'dovich effect and X-ray surface brightness to reconstruct the radial profiles of gas temperature and density under the assumption of a spherically symmetric distribution of the gas. The method of reconstruction, first raised by Silk & White, depends directly on the observations of the Sunyaev-Zel'dovich effect and the X-ray surface brightness, without involving additional assumptions such as the equation of state of the gas or the conditions of hydrostatic equilibrium. We applied this method to the cluster RX J1347.5-1145, which has both the Sunyaev-Zel'dovich effect and X-ray observations with relative high precision. It is shown that it will be an effective method to obtain the gas distribution in galaxy clusters. Statistical errors of the derived temperature and density profiles of gas were estimated according to the observational uncertainties.

In this paper, we use a set of observational $H(z)$ data (OHD) to constrain the $\Lambda$CDM cosmology. This data set can be derived from the differential ages of the passively evolving galaxies. Meanwhile, the $\mathcal {A}$-parameter, which describes the Baryonic Acoustic Oscillation (BAO) peak, and the newly measured value of the Cosmic Microwave Background (CMB) shift parameter $\mathcal {R}$ are used to present combinational constraints on the same cosmology. The combinational constraints favor an accelerating flat universe while the flat $\Lambda$CDM cosmology is also analyzed in the same way. We obtain a result compatible with that by many other independent cosmological observations. We find that the observational $H(z)$ data set is a complementarity to other cosmological probes.

Recently PAMELA released their first results on the positron and antiproton ratios. Stimulated by the new data, we studied the cosmic ray propagation models and calculated the secondary positron and antiproton spectra. The low energy positron ratio can be consistent with data in the convection propagation model. Above $\sim 10$ GeV PAMELA data shows a clear excess on the positron ratio. However, the secondary antiproton is roughly consistent with data. The positron excess may be a direct evidence of dark matter annihilation or decay. We compare the positron and anti-proton spectra with data by assuming dark matter annihilates or decays into different final states. The PAMELA data actually excludes quark pairs being the main final states, disfavors gauge boson final states. Only in the case of leptonic final states the positron and anti-proton spectra can be explained simultaneously. We also compare the decaying and annihilating dark matter scenarios to account for the PAMELA results and prefer to the decaying dark matter. Finally we consider a decaying neutralino dark matter model in the frame of supersymmetry with R-parity violation. The PAMELA data is well fitted with neutralino mass $600\sim 2000$ GeV and life time $\sim 10^{26}$ seconds. We also demonstrate that neutralino with mass around 2TeV can fit PAMELA and ATIC data simultaneously.

In this work we develop a new propagation model for the Galactic cosmic rays based on the GALPROP code, including contributions from dark matter annihilation. The model predicts compatible Galactic diffuse $\gamma$ ray spectra with EGRET data in all sky regions. It also gives consistent results of the diffuse $\gamma$ ray longitude and latitude distributions. Further the results for B/C, $^{10}$Be/$^9$Be, proton, electron and antiproton spectra are also consistent with cosmic ray measurements. In the model we have taken a universal proton spectrum throughout the Galaxy without introducing large fluctuation for the proton energy loss is negligible. The dark matter annihilation signals are `boosted' after taking the contribution from subhalos into account. Another interesting feature of the model is that it gives better description of the diffuse $\gamma$ rays when taking the source distribution compatible with supernova remnants data, which is different from previous studies.

The diffuse Galactic $\gamma$-rays from EGRET observation shows excesses above 1 GeV in comparison with the expectations from conventional Galactic cosmic ray (CR) propagation model. In the work we try to solve the "GeV excess" problem by dark matter (DM) annihilation in the frame of supersymmetry (SUSY). Compared with previous works, there are three aspects improved in this work: first, the direction-independent "boost factor" for diffuse $\gamma$-rays from dark matter annihilation (DMA) is naturally reproduced by taking the DM substructures into account; second, there is no need for renormalization of the diffuse $\gamma$-ray background produced by CRs; last but not the least, in this work our new propagation model can give consistent results of both diffuse $\gamma$-rays and antiprotons, by directly adding the signals from DMA to the diffuse $\gamma$-ray background. This is a self-consistent model among several possible scenarios at present, and can be tested or optimized by the forthcoming experiments such as GLAST, PAMELA and AMS02.

Dark matter (DM) annihilations in the Galaxy may produce high energy neutrinos, which can be detected by the neutrino telescopes, for example IceCube, ANTARES and Super-Kamiokande. The neutrinos can also arise from hadronic interaction between cosmic ray and atmosphere around the Earth, known as atmospheric neutrino. Current measurements on neutrino flux is consistent with theoretical prediction of atmospheric neutrino within the uncertainties. In this paper, by requiring that the DM annihilation neutrino flux is less than the current measurements, we obtain an upper bound on the cross section of dark matter annihilation $ < {\sigma v} >$. Compared with previous investigations, we improve the bound by including DM substructure contributions. In our paper, two kinds of substructure effects are scrutinized. One is the substructure average contribution over all directions. The other is point source effect by single massive sub-halo. We found that the former can improve the bound by several times, while the latter can improve the bound by $ 10^1 \sim 10^4$ utilizing the excellent angular resolution of neutrino telescope IceCube. The exact improvement depends on the DM profile and the sub-halo concentration model. In some model, IceCube can achieve the sensitivity of $ < {\sigma v} > \sim 10^{- 26} cm^3 s^{- 1} $.

It has been revealed recently that, in the scale free range, i.e. from the scale of the onset of nonlinear evolution to the scale of dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the self-similar log-Poisson hierarchy. As a consequence of this evolution, the relations among various physical quantities of cosmic baryon fluid should be scale invariant, if the physical quantities are measured in cells on scales larger than the dissipation scale, regardless the baryon fluid is in virialized dark halo, or in pre-virialized state. We examine this property with the relation between the Compton parameter of the thermal Sunyaev-Zel'dovich effect, $y(r)$, and X-ray luminosity, $L_{\rm x}(r)$, where $r$ being the scale of regions in which $y$ and $L_{\rm x}$ are measured. According to the self-similar hierarchical scenario of nonlinear evolution, one should expect that 1.) in the $y(r)$-$L_x(r)$ relation, $y(r)=10^{A(r)}[L_{\rm x}(r)]^{\alpha(r)}$, the coefficients $A(r)$ and $\alpha(r)$ are scale-invariant; 2.) The relation $y(r)=10^{A(r)}[L_{\rm x}(r)]^{\alpha(r)}$ given by cells containing collapsed objects is also available for cells without collapsed objects, only if $r$ is larger than the dissipation scale. These two predictions are well established with a scale decomposition analysis of observed data, and a comparison of observed $y(r)$-$L_x(r)$ relation with hydrodynamic simulation samples. The implication of this result on the characteristic scales of non-gravitational heating is also addressed.

Using a large sample of 26623 quasars with redshifts in the range $1.5 \le z\le 5.1$ with \civ $\lambda$1549 \AA emission line in Fifth Data Release of the Sloan Digital Sky Survey (SDSS), we investigate the cosmological evolution of the Baldwin Effect, i.e. the relation between the equivalent width (EW) of the \civ emission line and continuum luminosity. We confirm the earlier result that there exists a strong correlation between the \civ EW and the continuum luminosity, and we find that, up to $z\approx 5$, the slope of the Baldwin Effect seems to have no effect of cosmological evolution. A sub-sample of 13960 quasars with broad \civ $\lambda$1549 \AA emission line from SDSS is used to explore the origin of the Baldwin Effect. We find that \civ EW have a strong correlation with the mass of supermassive black hole (SMBH), and a weak correlation with the Eddington ratio, $\lb/\ledd$. This suggests that the SMBH mass is probably the primary drive for the Baldwin Effect.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

Anti-proton and positron Galactic cosmic ray spectra are among the key targets for indirect detection of dark matter. The boost factors, corresponding to an enhancement of the signal and linked to the clumpiness properties of the dark matter distribution, have been taken as high as thousands in the past. The dramatic impact of these boost factors for indirect detection of antiparticles, for instance with the PAMELA satellite or the coming AMS-02 experiment, asks for their detailed calculation. We take into account the results of high resolution N-body dark matter simulations to calculate the most likely energy dependent boost factors linked to the cosmic ray propagation properties, for anti-protons and positrons. Starting from the mass and space distributions of sub-halos, the anti-proton and positron propagators are used to calculate the mean value and the variance of the boost factor for the primary fluxes. We take advantage of the statistical method introduced in Lavalle et al. (2007) and cross-check the results with Monte Carlo computations. By spanning some extreme configurations of sub-halo and propagation properties, we find that the average contribution of the clumps is negligible compared to that of the smooth dark matter component. Sub-halos do not lead to enhancement of the signals, unless they are taken with some extreme (unexpected) properties. This result is independent of the nature of the self-annihilating dark matter candidate considered, and provides precise estimates of the theoretical and the statistical uncertainties of the antimatter flux from dark matter substructures. Spectral distortions can still be expected in antimatter flux measurements, but scenarios invoking large and even mild clumpiness boost factors are strongly disfavoured by our analysis.

In this paper we calculate the Galactic positron flux from dark matter annihilation in the frame of supersymmetry, taking the enhancement of the flux by existence of dark matter substructures into account. The propagation of positrons in the Galactic magnetic field is solved in a realistic numerical model GALPROP. The secondary positron flux is recalculated in the GLAPROP model. The total positron flux from secondary products and dark matter annihilation can fit the HEAT data well when taking a cuspy density profile of the substructures.

Three sets of two-component profiles are used to simultaneously model the [O III]$\lambda\lambda$4959, 5007 and H$\beta$ lines for the Fe II-subtracted spectra of 149 narrow line Seyfert 1 galaxies (NLSls) from the Sloan Digital Sky Survey (SDSS). Using the linewidth of the narrow/core component of [O III]$\lambda$5007 to trace the stellar velocity dispersion instead of using the total linewidth of [O III]$\lambda$5007, we found that the SDSS NLSls are still deviated from the $M_{bh}-\sigma$ relation found in the nearby inactive galaxies, suggesting that the linewidth of the [O III] narrow/core component is likely not a good tracer of bulge velocity dispersion in NLSls, since some other studies indicate that NLS1s, like other AGN, should follow the $M_{bh}-\sigma$ relation. If we assume that the [O III]5007/H$\beta_{n}$ line ratio emitted in narrow line region ranges from one to ten, 63 objects are found to satisfied with this criterion and their H$\beta$ broad components should be used to calculate their virial black hole masses. These 63 objects are statically consistent with the $M_{bh}-\sigma_{[O III]}$ relation. With the Chandra observation of some SDSS NLSls, for one object in these 63 objects, J143030.22-001115.1, we found that it can't be classified as a genuine NLS1. Its narrow component of H$\beta$ is coming from narrow line regions (NLRs). This is consistent with its very flat hard X-ray spectrum found by williams et al. (2004).

The progress of the studies on the orientation of galaxies in the Local Supercluster (LSC) is reviewed and a summary of recent results is given. Following a brief introduction of the LSC, we describe the results of early studies based on two-dimensional analysis, which were mostly not conclusive. We describe next the three-dimensional analysis, which is used widely today. Difficulties and systematic effects are explained and the importance of selection effects is described. Then, results based on the new method and modern databases are given, which are summarized as follows. When the LSC is seen as a whole, galaxy planes tend to align perpendicular to the LSC plane with lenticulars showing the most pronounced tendency. Projections onto the LSC plane of the spin vectors of Virgo cluster member galaxies, and to some extent, those of the total LSC galaxies, tend to point to the Virgo cluster center. This tendency is more pronounced for lenticulars than for spirals. It is suggested that 'field' galaxies, i.e., those which do not belong to groups with more than three members, may be better objects than other galaxies to probe the information at the early epoch of the LSC formation through the analysis of galaxy orientations. Field lenticulars show a pronounced anisotropic distribution of spin vectors in the sense that they lay their spin vectors parallel to the LSC plane while field spirals show an isotropic spin-vector distribution.

We investigate the radial velocity difference between the narrow emission-line components of [O III] $\lambda$ 5007 and H$\beta$ in a sample of 150 SDSS NLS1s. Seven "blue outliers" with [O III] blueshifted by more than 250 \kms are found. A strong correlation between the [O III] blueshift and the Eddington ratio is found for these seven "blue outliers". For the entire sample, we found a modest correlation between the blueshift and the linewidth of the narrow component of the [O III] line. The reflected profile of [O III] indicates two kinematically and physically distinct regions. The [O III] linewidth depends not only on the bulge stellar gravitational potential, but also on the central black hole potential.

The merging cluster of galaxies A2255 is covered by the Sloan Digital Sky Survey (SDSS) survey. In this paper we perform a morphological classification on the basis of the SDSS imaging and spectral data, and investigate the morphological dependence of the star formation rates (SFRs) for these member galaxies. As we expect, a tight correlation between the normalized SFR by stellar mass (SFR/M$_*$) and the H$\alpha$ equivalent width is found for the late-type galaxies in A2255. The correlation of SFR/M$_*$ with the continuum break strength at 4000 \AA is also confirmed. The SFR/M$_*$ - M$_*$ correlation is found for both the early- and late-type galaxies, indicating that the star formation activity tends to be suppressed when the assembled stellar mass M$_*$) increases, and this correlation is tighter and steeper for the late-type cluster galaxies. Compared with the mass range of field spiral galaxies, only two massive late-type galaxies with M$_*>10^{11}$ M$_{\odot}$ are survived in A2255, suggesting that the gas disks of massive spiral galaxies could have been tidally stripped during cluster formation. Additionally, the SFR variation with the projected radial distance are found to be heavily dependent upon galaxy morphology: the early-type galaxies have a very weak inner decrease in SFR/M$_*$, while the inner late-type galaxies tend to have higher SFR/M$_*$ values than the outer late-types. This may suggest that the galaxy-scale turbulence stimulated by the merging of subclusters might have played different roles on early- and late-type galaxies, which leads to a suppression of the star formation activity for E/S0 galaxies and a SFR enhancement for spiral and irregular galaxies.

Convincing evidence of a past interaction between two rich clusters A399 and A401 was recently found by the X-ray imaging observations. In this paper we examine the structure and dynamics of this pair of galaxy clusters. A mixture-modeling algorithm has been applied to obtain a robust partition into two clusters, which allows us to discuss the virial mass and velocity distribution for each cluster. Assuming that these two clusters follow a linear orbit and they have once experienced a close encounter, we model the binary cluster as a two-body system. As a result, four gravitationally bound solutions are obtained. The recent X-ray observations seem to favor a scenario in which the two clusters with a true separation of $5.4h^{-1}$ Mpc are currently expanding at 583 km/s along the direction with a projection angle of 67.5 degree, and they will reach a maximum extent of $5.65h^{-1}$ Mpc in about $1.0h^{-1}$ Gyr.

This paper presents a multicolor study of the nearby cluster of galaxies Abell 168 (z=0.045) with 13 intermediate-band filters in the Beijing-Arizona-Taiwan-Connecticut (BATC) filter system. After a cross-identification between the photometric data obtained from the BATC and the Sloan Digital Sky Survey (SDSS), a catalog containing 1553 galaxies down to r'<20.0 mag is achieved, which includes 121 spectroscopically confirmed member galaxies. The technique of photometric redshift has been applied to all these galaxies with combined 18-band (13 from BATC and 5 from SDSS) spectral energy distributions (SEDs), in order to perform a faint membership selection in Abell 168. As a result, 255 galaxies are newly selected as the member candidates. Based on the enlarged sample of cluster galaxies, the spatial distribution and dynamics of Abell 168 are investigated. In the light of the spatial distribution of the member galaxies and the 0.2-3.5 keV X-ray image by Einstein observatory, it seems that Abell 168 consists of two merging subclusters with a relative radial velocity of 264 +/- 142 km/s. With the help of ROSTAT software, a detailed investigation of the dynamics shows the intrinsic difference in the velocity distributions for these two subclusters. Under a linear two-body model, they are found to be a gravitationally bound system with 92% probability. The slight deviation of the local velocity distribution from the overall distribution in the central region may suggest a picture of head-on collision in the projection plane. By using a double Schechter function, the close examination of the luminosity functions reveals the different galaxy content of these two subclusters, which implies the merging process is possibly at the early stage.

We present our optical multicolor photometry for the nearby cluster of galaxies Abell 2255 with 13 intermediate filters in the BATC system. The spectral energy distributions (SEDs) in the optical band for more than 7000 sources are achieved down to V ~ 20 in a field of 58' x 58' centered on this rich cluster. Abell 2255 has been recently observed by the SDSS. A method of combining the SDSS photometric data in five broad bands and the BATC SEDs is then explored. The application of the technique of photometric redshift shows that the combined SEDs with higher resolution could lead to a more accurate estimate of photometric redshift. Based on 214 spectroscopically confirmed member galaxies, the bimodality and large dispersion in the velocity distribution indicate that Abell 2255 is an unrelaxed system. A tight color-magnitude correlation for 188 known early-type cluster galaxies is found. After an exclusion of 254 extragalactic sources with known redshifts, the combined SEDs for 2522 galaxies allow a further membership selection by the photometric redshift technique and color-magnitude correlation. As a result, 313 galaxies with their photometric redshifts between 0.068 and 0.090 are selected as new cluster members. The combined SEDs and estimated redshifts for previously known and newly-selected member galaxies are cataloged. On the basis of the enlarged sample of member galaxies, the spatial distribution, localized velocity structure, color-magnitude relation, and luminosity function of cluster galaxies are discussed. The reverse peculiar velocities are found for two satellite subclusters located on opposite sides of the central concentration, which supports an ongoing merger in Abell 2255.

As one of the most luminous radio-loud quasars showing intrinsic ultraviolet (UV) and X-ray absorption, 3C 351 provides a laboratory for studying the kinematics and physical conditions of such ionized absorbers. We present an analysis of the intrinsic absorption lines in the high-resolution ($\sim$ 7 km/s) far-UV spectrum which was obtained from observations with the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST). The spectrum spans wavelengths from 1150 \AA to 1710 \AA, and shows strong emission lines from O VI and Ly$\alpha$. Associated absorption lines are present on the blue wings of the high-ionization emission doublets O VI $\lambda\lambda$ 1032,1038 and N V $\lambda\lambda$ 1238,1242, as well as the Lyman lines through Ly$\epsilon$. These intrinsic absorption features are resolved into several distinct kinematic components, covering rest-frame velocities from -40 to -2800 km/s, with respect to the systemic redshift of $z_{em}=0.3721$. For the majority of these absorption line regions, strong evidence of partial covering of both the background continuum source and the BELR is found, which supports the intrinsic absorption origin and rules out the possibility that the absorption arises in some associated cluster of galaxies. The relationship between the far-UV absorbers and X-ray `warm' absorbers are studied with the assistance of photoionization models. Most of the UV associated absorption components have low values of the ionization parameter and total hydrogen column densities, which is inconsistent with previous claims that the UV and X-ray absorption arises in the same material. Analysis of these components supports a picture with a wide range of ionization parameters, temperatures, and column densities in AGN outflows.

An optical photometric observation with the BATC multicolor system is carried out for the central 56' x 56' region of the nearby cluster of galaxies, Abell 2634. We achieved the spectral energy distributions (SEDs) of 5572 sources detected down to $V \sim 20$ mag, including 178 previously known galaxies, with fourteen filters covering a wavelength range from 3600 \AA\ to 10000 \AA. After excluding the foreground and background galaxies, a sample of 124 known members is formed for an investigation of the SED properties. Based on the knowledge of SED properties of member galaxies, we performed the selection of faint galaxies belonging to Abell 2634. The color-color diagrams are powerful in the star/galaxy separation, and 359 faint galaxies are selected by their color features. The technique of photometric redshift and color-magnitude correlation for the early-type galaxies are applied for these faint galaxies, and a list of 74 faint member galaxies is achieved. Basis on the enlarged sample of member galaxies, the spatial distribution and color-magnitude relation of the galaxies in core region of Abell 2634 are discussed. We find a tendency that the color index dispersion of the early-type members is larger for the outer region, which might reflect some clues about the environmental effect on the evolution of galaxies in a cluster.