Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

The interstellar oxygen isotopic ratio of 18O/17O can reflect the relative amount of the secular enrichment by ejecta from high-mass versus intermediate-mass stars. Previous observations found a Galactic gradient of 18O/17O, i.e., low ratios in the Galactic center and large values in the Galactic disk, which supports the insideout formation scenario of our Galaxy. However, the observed objects are not many and, in particular, not so many at large galactocentric distances. Thus we started a systematic study on Galactic interstellar 18O/17O, through observations of C18O and C17O multi-transition lines toward a large sample of 286 sources (at least one order of magnitude larger than previous ones), from the Galactic center region to the far outer Galaxy (~22 kpc). Here we present our observations of J=1-0 lines of C18O and C17O, with the ARO12m and the IRAM 30m telescope. We detected successfully both C18O and C17O 1-0 lines for 34 sources among our IRAM30m sample of 50 targets and for 166 sources among our ARO12m sample of 260 targets. The C18O optical depth effect on our ratio results, evaluated by fitting results of C17O spectra with hyperfine components and our RADEX non-LTE model calculation for the strongest source, was found to be insignificant. Beam dilution does not seem to be a problem either, which was supported by the fact of no systematic variation between the isotopic ratio and the heliocentric distance, and consistent measured ratios from two telescopes for most of those detected sources. Our results, though there are still very few detections made for sources in the outer Galaxy, confirm the apparent 18O/17O gradient of 18O/17O = (0.10+-0.03)R_GC+(2.95+-0.30), with a Pearson's rank correlation coefficient R = 0.69. This is supported by the newest Galactic chemical evolution model including the impact of massive stellar rotators and novae.

We present spectropolarimetric and spectrophotometric observations of the peculiar active galactic nucleus (AGN) SDSS J120300.19+162443.7 (hereafter J1203+1624) at $z=0.1656$. Its optical total flux spectra clearly show broad emission lines (BELs) in H$\alpha$ and H$\beta$. After removal of narrow emission lines (NELs), the full width at half maximum of the lorentzian BEL is $FWHM_{\rm BEL} \approx 1,270$ km s$^{-1}$, fulfilling the conventional definition of a narrow line Seyfert 1 (NLS1) galaxy. However, its NELs are unprecedentedly strong when compared to type 1AGNs. This, together with its large MIR excess ($g-W_4 = 13.172$ mag), implies that the observer and the narrow emission line region (NELR) might see a different ionization continuum. Our optical spectropolarimetry confirms its type 2 nature by detecting a polarized blue continuum and Balmer BELs ($FWHM_{\rm Polarized~BEL} \approx 1,183$ km s$^{-1}$), with a high polarization degree of $> 20$% in the blue wing. The mass and Eddington rate of the central black hole is estimated based on both transmitted and scattered AGN radiation, which is $M_{\bullet} < 2.9 \times 10^7 M_\odot$ and $L_{\rm bol}/L_{\rm Edd} > 1.5$. Severe extinction of the AGN emission also enables a clear view of the compact host galaxy. Discovery of J1203+1624 suggests that NLS1 follows AGN unfication scheme, studying of its analogs could blaze a new trail for exploring the connection between black hole growth and star formation in the host galaxy. Interesting features of J1203+1624 like the peculiar NELs, inflowing scattering clouds within sublimation radius are worth detailed follow-ups in the future.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

A new, complete sample of 14,584 broad-line AGNs at $z<0.35$ is presented, which are uncovered homogeneously from the complete database of galaxies and quasars observed spectroscopically in the Sloan Digital Sky Survey Seventh Data Release. The stellar continuum is properly removed for each spectrum with significant host absorption line features, and careful analyses of the emission-line spectra, particularly in the H$\rm \alpha$ and H$\rm \beta$ wavebands, are carried out. The broad Balmer emission line, particularly, H$\rm \alpha$, is used to indicate the presence of an AGN. The broad H$\rm \alpha$ lines have luminosities in a range of $10^{38.5}$-$10^{44.3}$ erg s$^{-1}$, and line widths (FWHMs) of 500-34,000 km s$^{-1}$. The virial black hole masses, estimated from the broad line measurements, span a range of $10^{5.1}$-$10^{10.3}$ $M_\odot$, and the Eddington ratios vary from $-3.3$ to $1.3$ in logarithmic scale. Other quantities such as multi-wavelength photometric properties and flags denoting peculiar line profiles are also included in this catalog. We describe the construction of this catalog and briefly discuss its properties. The catalog is publicly available online. This homogeneously selected AGN catalog, along with the accurately measured spectral parameters, provide the most updated, largest AGN sample data, which will enable further comprehensive investigations of the properties of the AGN population in the low-redshift universe.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

Vortices play an unique role in heat and momentum transports in astro- and geo-physics, and it is also the origin of the Earth's dynamo. A question existing for a long time is whether the movement of vortices can be predicted or understood based on their historical data. Here we use both the experiments and numerical simulations to demonstrate some generic features of vortex motion and distribution. It can be found that the vortex movement can be described on the framework of Brownian particles where they move ballistically for the time shorter than some critical timescales, and then move diffusively. Traditionally, the inertia of vortex has often been neglected when one accounts for their motion, our results imply that vortices actually have inertial-induced memory such that their short term movement can be predicted. Extending to astro- and geo-physics, the critical timescales of transition are in the order of minutes for vortices in atmosphere and ocean, in which this inertial effect may often be neglected compared to other steering sources. However, the timescales for vortices are considerably larger which range from days to a year. It infers the new concept that not only the external sources alone, for example the solar wind, but also the internal source, which is the vortex inertia, can contribute to the short term Earth's magnetic field variation.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

The limiting rotational frequency of neutron stars may be determined by the strength of their crusts. As a star spins up from accretion, centrifugal forces will cause the crust to fail. If the crust breaks unevenly, a rotating mass quadrupole moment will radiate gravitational waves (GW). This radiation can prevent further spin up and may be a promising source for continuous GW searches. We calculate that for a breaking strain (strength) of neutron star crust that is consistent with molecular dynamics simulations, the crust may fail at rotational frequencies in agreement with observations.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

We present the cold neutral content (H I and C I gas) of 13 quasar 2175 \AA$ $ dust absorbers (2DAs) at $z$ = 1.6 - 2.5 to investigate the correlation between the presence of the UV extinction bump with other physical characteristics. These 2DAs were initially selected from the Sloan Digital Sky Surveys I - III and followed up with the Keck-II telescope and the Multiple Mirror Telescope as detailed in our Paper I. We perform a correlation analysis between metallicity, redshift, depletion level, velocity width, and explore relationships between 2DAs and other absorption line systems. The 2DAs on average have higher metallicity, higher depletion levels, and larger velocity widths than Damped Lyman-$\alpha$ absorbers (DLAs) or subDLAs. The correlation between [Zn/H] and [Fe/Zn] or [Zn/H] and log$\Delta$V$_{90}$ can be used as alternative stellar mass estimators based on the well-established mass-metallicity relation. The estimated stellar masses of the 2DAs in this sample are in the range of $\sim$ 10$^9$ to $\sim$2 $\times$ 10$^{11}$ $M_{\odot}$ with a median value of $\sim$2 $\times$ 10$^{10}$ $M_{\odot}$. The relationship with other quasar absorption line systems can be described as (1) 2DAs are a subset of Mg II and Fe II absorbers, (2) 2DAs are preferentially metal-strong DLAs/subDLAs, (3) More importantly, all of the 2DAs show C I detections with logN(C I) $>$ 14.0 cm$^{-2}$, (4) 2DAs can be used as molecular gas tracers. Their host galaxies are likely to be chemically enriched, evolved, massive (more massive than typical DLA/subDLA galaxies), and presumably star-forming galaxies.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

By the end of 2014, a cosmic ray muon telescope was installed at Zhongshan Station in Antarctic and has been continuously collecting data since then. It is the first surface muon telescope to be built in Antarctic. In June 2015, five CMEs were ejected towards the Earth initiating a big large Forbush decrease (FD) event. We conduct a comprehensive study of the galactic cosmic ray intensity fluctuations during the FD using the data from cosmic ray detectors of multiple stations (Zhongshan, McMurdo, South Polar and Nagoya) and he solar wind measurements from ACE and WIND. A pre-increase before the shock arrival was observed. Distinct differences exist in the timelines of the galactic cosmic ray recorded by the neutron monitors and the muon telescopes. FD onset for Zhongshan muon telescope is delayed (2.5h) with respect to SSC onset. This FD had a profile of four-step decrease. The traditional one- or two-step classification of FDs was inadequate to explain this FD.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

We obtained new high-resolution spectra using the Lijiang 1.8-m and 2.4-m telescopes to investigate the chromospheric activities of V1355 Ori as indicated in the behaviors of Ca ii H&K, Hdelta, Hgamma, Hbeta, Na i D1, D2, Halpha and Ca ii infrared triplet (IRT) lines. The observed spectra show obvious emissions above the continuum in Ca ii H&K lines, absorptions in the Hdelta, Hgamma, Hbeta and Na i D1, D2 lines, variable behavior (filled-in absorption, partial emission with a core absorption component or emission above the continuum) in the Halpha line, and weak self-reversal emissions in the strong filled-in absorptions of the Ca ii IRT lines. We used a spectral subtraction technique to analyze our data. The results show no excess emission in the Hdelta and Hgamma lines, very weak excess emissions in the Na i D1, D2 lines, excess emission in the Hbeta line, clear excess emission in the Halpha line, and excess emissions in the Ca ii IRT lines. The value of the ratio of EW8542/EW8498 is in the range 0.9 to 1.7, which implies that chromospheric activity might have been caused by plage events. The value of the ratio E Halpha/E Hbeta is above 3, indicating that the Balmer lines would arise from prominence-like material. We also found time variations in light curves associated with equivalent widths of chromospheric activity lines in the Na i D1, D2, Ca ii IRT and Halpha lines in particular. These phenomena can be explained by plage events, which are consistent with the behavior of chromospheric activity indicators.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

We report the identification of an unusual absorption line system in the quasar SDSS J080248.18$+$551328.9 and present a detailed study of the system, incorporating follow-up optical and NIR spectroscopy. A few tens of absorption lines are detected, including He I*, Fe II* and Ni II* that arise from metastable or excited levels, as well as resonant lines in Mg I, Mg II, Fe II, Mn II, and Ca II. All of the isolated absorption lines show the same profile of width $\Delta v\sim 1,500$km s$^{-1}$ centered at a common redshift as that of the quasar emission lines, such as [O II], [S II], and hydrogen Paschen and Balmer series. With narrow Balmer lines, strong optical Fe II multiplets, and weak [O III] doublets, its emission line spectrum is typical for that of a narrow-line Seyfert 1 galaxy (NLS1). We have derived reliable measurements of the gas-phase column densities of the absorbing ions/levels. Photoionization modeling indicates that the absorber has a density of $n_{\rm H} \sim (1.0-2.5)\times 10^5~ {\rm cm}^{-3}$ and a column density of $N_{\rm H} \sim (1.0-3.2)\times 10^{21} \sim {\rm cm}^{-2}$, and is located at $R\sim100-250$ pc from the central super-massive black hole. The location of the absorber, the symmetric profile of the absorption lines, and the coincidence of the absorption and emission line centroid jointly suggest that the absorption gas is originated from the host galaxy and is plausibly accelerated by stellar processes, such as stellar winds \zhy{and/or} supernova explosions. The implications for the detection of such a peculiar absorption line system in an NLS1 are discussed in the context of co-evolution between super-massive black hole growth and host galaxy build-up.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

We present a three dimensional extinction map in $r$ band. The map has a spatial angular resolution, depending on latitude, between 3 -- 9\,arcmin and covers the entire XSTPS-GAC survey area of over 6,000\,$\rm deg^2$ for Galactic longitude $\rm 140 \leq$ $l$ $ \leq 220\deg$ and latitude $\rm -40\leq$ $b$ $ \leq 40\deg$. By cross-matching the photometric catalog of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) with those of 2MASS and WISE, we have built a multi-band photometric stellar sample of about 30 million stars and applied spectral energy distribution (SED) fitting to the sample. By combining photometric data from the optical to the near-infrared, we are able to break the degeneracy between the intrinsic stellar colours and the amounts of extinction by dust grains for stars with high photometric accuracy, and trace the extinction as a function of distance for low Galactic latitude and thus highly extincted regions. This has allowed us to derive the best-fit extinction and distance information of more than 13 million stars, which are used to construct the three dimensional extinction map. We have also applied a Rayleigh-Jeans colour excess (RJCE) method to the data using the 2MASS and WISE colour $(H-W2)$. The resulting RJCE extinction map is consistent with the integrated two dimensional map deduced using the best-fit SED algorithm. However for individual stars, the amounts of extinction yielded by the RJCE method suffer from larger errors than those given by the best-fit SED algorithm.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

A new multi-wavelength solar telescope, Optical and Near-infrared Solar Eruption Tracer (ONSET) of Nanjing University, was constructed, being fabricated by Nanjing Institute of Astronomical Optics & Technology and run in cooperation with Yunnan Astronomical Observatory. ONSET is able to observe the Sun in three wavelength windows: He {\small I} 10830 \AA, H$\alpha$, and white-light at 3600 \AA and 4250 \AA, which are selected in order to obtain the dynamics in the corona, chromosphere, and the photosphere simultaneously. Full-disk or partial-disk solar images with a field of 10 arcmin at three wavelengths can be obtained nearly simultaneously. It is designed to trace solar eruptions with high spatial and temporal resolutions. This telescope was installed at a new solar observing site near Fuxian Lake in Yunnan Province, southwest China. The site is located at E102N24, with an altitude of 1722 m. The seeing is stable and of high quality. We give a brief description of the scientific objectives and the basic structure of the ONSET. Some preliminary results are also presented.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

The ensemble optical/ultraviolet variability of narrow-line Seyfert 1 (NLS1) type active galactic nuclei (AGNs) is investigated, based on a sample selected from the Sloan Digital Sky Survey (SDSS) Stripe-82 region with multi-epoch photometric scanning data. As a comparison a control sample of broad-line Seyfert 1 (BLS1) type AGNs is also incorporated. To quantify properly the intrinsic variation amplitudes and their uncertainties, a novel method of parametric maximum-likelihood is introduced, that has, as we argued, certain virtues over previously used methods. The majority of NLS1-type AGNs exhibit significant variability on timescales from about ten days to a few years with, however, on average smaller amplitudes compared to BLS1-type AGNs. About 20 NLS1- type AGNs showing relatively large variations are presented, that may deserve future monitoring observations, for instance, reverberation mapping. The averaged structure functions of variability, constructed using the same maximumlikelihood method, show remarkable similarity in shape for the two types of AGNs on timescales longer than about 10 days, which can be approximated by a power-law or an exponential function. This, along with other similar properties, such as the wavelength-dependent variability, are indicative of a common dominant mechanism responsible for the long-term optical/UV variability of both NLS1- and BLS1-type AGNs. Towards the short timescales, however, there is tentative evidence that the structure function of NLS1-type AGNs continues declining, whereas that of BLS1-type AGNs flattens with some residual variability on timescales of days. If this can be confirmed, it may suggest that an alternative mechanism, such as X-ray reprocessing, starts to become dominating in BLS1-type AGNs, but not in NLS1-, on such timescales.

The durations (T90) of 315 GRBs detected with Fermi/GBM (8-1000 keV) by 2011 September are calculated using the Bayesian Block method. We compare the T90 distributions between this sample and those derived from previous/current GRB missions. We show that the T90 distribution of this GRB sample is bimodal, with a statistical significance level being comparable to those derived from the BeppoSAX/GRBM sample and the Swift/BAT sample, but lower than that derived from the CGRO/BATSE sample. The short-to-long GRB number ratio is also much lower than that derived from the BATSE sample, i.e., 1:6.5 vs 1:3. We measure T90 in several bands, i.e., 8-15, 15-25, 25-50, 50-100, 100-350, and 350-1000 keV, to investigate the energy-dependence effect of the bimodal T90 distribution. It is found that the bimodal feature is well observed in the 50-100 and 100-350 keV bands, but is only marginally acceptable in the 25-50 keV and 350-1000 keV bands. The hypothesis of the bimodality is confidently rejected in the 8-15 and 15-25 keV bands. The T90 distributions in these bands are roughly consistent with those observed by missions with similar energy bands. The parameter T90 as a function of energy follows \bar T90 \propto E^{-0.20\pm 0.02} for long GRBs. Considering the erratic X-ray and optical flares, the duration of a burst would be even much longer for most GRBs. Our results, together with the observed extended emission of some short GRBs, indicate that the central engine activity time scale would be much longer than T90} for both long and short GRBs and the observed bimodal T90 distribution may be due to an instrumental selection effect.

(Abridged for arXiv) We report the discovery of an unusual, extremely dust-rich and metal-strong damped \lya absorption system (DLA) at a redshift $z_{a}=2.4596$ toward the quasar SDSS J115705.52+615521.7 (hereafter J1157+6155) with an emission-line redshift $z_{e}=2.5125$. Its neutral hydrogen column density $N_{\hi} = 10^{21.8\pm0.2}$ cm$^{-2}$ is among the highest values measured in quasar DLAs. The measured metal column density is $N_{ZnII}\approx 10^{13.8}$ cm$^{-2}$, which is about 1.5 times larger than the largest value in any previously observed quasar DLAs. The best-fit curve is a MW-like law with a significant broad feature centered around 2175 {\AA} in the rest frame of the absorber. The measured extinction $A_V \approx 0.92$ mag is unprecedentedly high in quasar DLAs. After applying an extinction correction, the $i$ band absolute magnitude of the quasar is as high as $M_{i} \approx -29.4$ mag, placing it one of the most luminous quasars ever known. This discovery is suggestive of the existence of a rare yet important population of dust-rich DLAs with both high metallicities and high column densities, which may have significant impact on the measurement of the cosmic evolution of neutral gas mass density and metallicity.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.

We present a statistical study of a large, homogeneously analyzed sample of narrow-line Seyfert 1 (NLS1) galaxies, accompanied by a comparison sample of broad-line Seyfert 1 (BLS1) galaxies. Optical emission-line and continuum properties are subjected to correlation analyses, in order to identify the main drivers of active galactic nuclei (AGN) correlation space, and of NLS1 galaxies in particular. For the first time, we have established the density of the narrow-line region as a key parameter in Eigenvector 1 space, as important as the Eddington ratio L/Ledd. This is important because it links the properties of the central engine with the properties of the host galaxy; i.e., the interstellar medium (ISM). We also confirm previously found correlations involving the line width of Hbeta, and the strength of the FeII and [OIII]5007 emission lines, and we confirm the important role played by L/Ledd in driving the properties of NLS1 galaxies. A spatial correlation analysis shows that large-scale environments of the BLS1 and NLS1 galaxies of our sample are similar. If mergers are rare in our sample, accretion-driven winds on the one hand, or bar-driven inflows on the other hand, may account for the strong dependence of Eigenvector 1 on ISM density.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

We present empirical relations between the global dust reddening and other physical galaxy properties including the Halpha luminosity, Halpha surface brightness, metallicity and axial ratio for star-forming disc galaxies. The study is based on a large sample of ~22 000 well-defined star-forming galaxies selected from the Sloan Digital Sky Survey (SDSS). The reddening parameterized by color excess E(B-V) is derived from the Balmer decrement. Besides the dependency of reddening on Halpha luminosity / surface brightness and gas phase metallicity, it is also correlated with the galaxy inclination, in the sense that edge-on galaxies are more attenuated than face-on galaxies at a give intrinsic luminosity. In light of these correlations, we present the empirical formulae of E(B-V) as a function of these galaxy properties, with a scatter of only 0.07 mag. The empirical relation can be reproduced if most dust attenuation to the HII region is due to diffuse interstellar dust distributing in a disc thicker than that of HII regions. The empirical formulae can be incorporated into semi-analytical models of galaxy formation and evolution to estimate the dust reddening and enable comparison with observations more practically.

Variable super-strong coronal emission lines were observed in one galaxy, SDSS J095209.56+214313.3, and their origin remains controversy. In this paper, we report the detection of variable broad spectral bumps, reminiscent of supernova (SN) II-Plateau (II-P) spectra taken a few days after the shock breakout, in the second galaxy with variable super-strong coronal lines, SDSS J074820.67+471214.3. The coronal line spectrum shows unprecedented high ionization with strong [Fe X], [Fe XI], [Fe XIV], [S XII] and [Ar XIV], but without detectable optical [Fe VII] lines. The coronal line luminosities are similar to that observed in bright Seyfert galaxies, and 20 times more luminous than that reported in the hottest Type IIn SN 2005ip. The coronal lines ($\sigma ~120-240$ km s-1) are much broader than the narrow lines ($\sigma \sim 40$ km/s) from the star forming regions in the galaxy, but at nearly the same systematic redshift. We also detected a variable non-stellar continuum in optical and UV. In the follow-up spectra taken 4-5 years later, the coronal lines, SN-like feature, and non-stellar continuum disappeared, while the [O III] intensity increased by about a factor of ten. Our analysis suggests that the coronal line region should be at least ten light days in size, and be powered either by a quasi-steady ionizing source with a soft X-ray luminosity at least a few 10^{42} erg s-1 or by a very luminous soft X-ray outburst. These findings can be more naturally explained by a star tidally disrupted by the central black hole than by an SN explosion.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

We report on a dusty Mg~II absorber associated with the quasar SDSSJ003545.13+011441.2 (hereafter J0035+0114) at $z$=1.5501, which is the strongest one among the three Mg~II absorbers along the sight line of quasar. The two low redshift intervening absorbers are at $z$=0.7436, 0.5436, respectively. Based on the photometric and spectroscopic data of Sloan Digital Sky Survey (hereafter SDSS), we infer the rest frame color excess E(\bv) due to the associated dust is more than 0.07 by assuming a Small Magellanic Cloud (hereafter SMC) type extinction curve. Our follow-up moderate resolution spectroscopic observation at the 10-m Keck telescope with the ESI spectrometer enable us to reliably identify most of the important metal elements, such as Zn, Fe, Mn, Mg, Al, Si, Cr, and Ni in the associated system. We measure the column density of each species and detect significant dust depletion. In addition, we develop a simulation technique to gauge the significance of 2175-{\AA} dust absorption bump on the SDSS quasar spectra. By using it, we analyze the SDSS spectrum of J0035+0114 for the presence of a associated 2175-{\AA} extinction feature and report a tentative detection at $\sim$2$\sigma$ significant level.

We investigate the cosmological evolution of the tachyon and phantom-tachyon scalar field by considering the potential parameter $\Gamma$($=\frac{V V"}{V'^2}$) as a function of another potential parameter $\lambda$($=\frac{V'}{\kappa V^{3/2}}$), which correspondingly extends the analysis of the evolution of our universe from two-dimensional autonomous dynamical system to the three-dimension. It allows us to investigate the more general situation where the potential is not restricted to inverse square potential and .One result is that, apart from the inverse square potential, there are a large number of potentials which can give the scaling and dominant solution when the function $\Gamma(\lambda)$ equals $3/2$ for one or some values of $\lambda_{*}$ as well as the parameter $\lambda_{*}$ satisfies condition Eq.(18) or Eq.(19). We also find that for a class of different potentials the dynamics evolution of the universe are actually the same and therefore undistinguishable.

The dependence of the long-term optical/UV variability on the spectral and the fundamental physical parameters for radio-quiet active galactic nuclei (AGNs) is investigated. The multi-epoch repeated photometric scanning data in the Stripe-82 region of the Sloan Digital Sky Survey (SDSS) are exploited for two comparative AGN samples (mostly quasars) selected therein, a broad-line Seyfert\,1 (BLS1) type sample and a narrow-line Seyfert\,1 (NLS1) type AGN sample within redshifts 0.3--0.8. Their spectral parameters are derived from the SDSS spectroscopic data. It is found that on rest-frame timescales of several years the NLS1-type AGNs show systematically smaller variability compared to the BLS1-type. In fact, the variability amplitude is found to correlate, though only moderately, with the Eigenvector\,1 parameters, i.e., the smaller the \hb\ linewidth, the weaker the [O\,III] and the stronger the \feii\ emission, the smaller the variability amplitude is. Moreover, an interesting inverse correlation is found between the variability and the Eddington ratio, which is perhaps more fundamental. The previously known dependence of the variability on luminosity is not significant, and that on black hole mass---as claimed in recent papers and also present in our data---fades out when controlling for the Eddington ratio in the correlation analysis, though these may be partly due to the limited ranges of luminosity and black hole mass of our samples. Our result strongly supports that an accretion disk is likely to play a major role in producing the opitcal/UV variability.

$Om$ diagnostic can differentiate between different models of dark energy without the accurate current value of matter density. We apply this geometric diagnostic to dilaton dark energy(DDE) model and differentiate DDE model from LCDM. We also investigate the influence of coupled parameter $\alpha$ on the evolutive behavior of $Om$ with respect to redshift $z$. According to the numerical result of $Om$, we get the current value of equation of state $\omega_{\sigma0}$=-0.952 which fits the WMAP5+BAO+SN very well.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

Quasar outflows may play important role in the evolution of its host galaxy and central black hole. In this paper, we present a detailed study of multiple outflows in the obscured ultra-luminous infrared quasar Q1321+058. The outflows reveal themselves in the complex optical and UV emission line spectrum, with a broad component blueshifted by 1650 km/s and a narrow component by 360 km/s, respectively.The higher velocity component shows ever strong N III] and strong Si III], in addition to strong [O III]5007 and [Ne III]3869 emission lines, suggesting an overabundance of N and Si relative to C. The abundance pattern is consistent a fast chemical enriching process associated with a recent starburst. The outflow extends to several tens to hundred parsecs from the quasar, and covers only a very small sky. We find that the outflow with line emitting gas is energetically insufficient to remove the ISM of the host galaxy. The velocity range and the column density suggest that the outflow might be part of the low ionization broad absorption line region as seen in a small class of quasars. The optical and UV continuum is starlight-dominated and can be modelled with a young-aged (1 Myr) plus an intermediate-aged (~0.5-1 Gyr) stellar populations, suggesting a fast building of the stellar mass in the host galaxy, consistent with the starburst-type metal abundances inferred from the high velocity outflow spectrum. The broad band spectral energy distribution shows that it is an obscured quasar with its bulk emission in the middle infrared. The star formation rate, independently estimated from UV, far-infrared, and emission line luminosity, is much lower than that is required for the co-evolution of the black hole and its host spheroid.

The galaxy SDSSJ0952+2143 showed remarkable emission-line properties first reported in 2008 (paper I), which are the consequence of a powerful high-energy flare. Here we report follow-up observations of SDSSJ0952+2143, and discuss outburst scenarios in terms of stellar tidal disruption by a SMBH, peculiar variability of an AGN, and a supernova explosion. The optical spectrum of SDSSJ0952+2143 exhibits several peculiarities: an exceptional ratio of [FeVII] transitions over [OIII], a dramatic decrease by a factor of 10 of the highest-ionization lines, a very unusual and variable Balmer line profile including a triple-peaked narrow component with two unresolved horns, and a large Balmer decrement. The MIR emission measured with the Spitzer IRS in the narrow 10-20mu band is extraordinarily luminous (3.5 x 10^{43} erg\s). The IRS spectrum shows a bump around ~11mu and an increase towards longer wavelengths, reminiscent of silicate emission. The strong MIR excess over the NIR implies the dominance of relatively cold dust. The X-ray luminosity of 10^{41} erg\s measured with Chandra is below that typically observed in AGN. Similarities of SDSSJ0952+2143 with some extreme supernovae suggest the explosion of a supernova of Type IIn. However, an extreme accretion event in a low-luminosity AGN or inactive galaxy, especially stellar tidal disruption, remain possibilities, which could potentially produce a very similar emission-line response. If indeed a supernova, SDSSJ0952+2143 is one of the most distant X-ray and MIR detected SNe known so far, the most MIR luminous, and one of the most X-ray luminous. It is also by far the most luminous (>10^{40} erg\s) in high-ionization coronal lines, exceeding previous SNe by at least a factor of 100 [abridged].

By considering the potential parameter $\Gamma$ as a function of another potential parameter $\lambda$[47], We successfully extend the analysis of two-dimensional autonomous dynamical system of quintessence scalar field model to the analysis of three-dimension, which makes us be able to research the critical points of a large number of potentials beyond the exponential potential exactly. We find that there are ten critical points in all, three points $P_{3, 5, 6}$} are general points which are possessed by all quintessence models regardless of the form of potentials and the rest points are closely connected to the concrete potentials. It is quite surprising that, apart from the exponential potential, there are a large number of potentials which can give the scaling solution when the function $f(\lambda)(=\Gamma(\lambda)-1)$ equals zero for one or some values of $\lambda_{*}$ and if the parameter $\lambda_{*}$ also satisfies the condition Eq.(16) or Eq.(17) at the same time. We give the differential equations to derive these potentials $V(\phi)$ from $f(\lambda)$. We also find that, if some conditions are satisfied, the de-Sitter-like dominant point $P_4$ and the scaling solution point $P_9$(or $P_{10}$) can be stable simultaneously but $P_9$ and $P_{10}$ can not be stable simultaneity. Although we survey scaling solutions beyond the exponential potential for ordinary quintessence models in standard general relativity, this method can be applied to other extensively scaling solution models studied in literature[46] including coupled quintessence, (coupled-)phantom scalar field, k-essence and even beyond the general relativity case $H^2 \propto\rho_T^n$. we also discuss the disadvantage of our approach.

We have conducted a series of investigations on the geometry of the reddening material in AGNs, which have important implications to the AGN unification and SMBH demography. According to our statistics of partially obscured quasars, we found that SMBHs in partially obscured type/phase (i.e., intermediate type) are at least as abundant as normal quasars in the local Universe; the reddening material in most objects are located in between the BLR and NLR. According to our comparison of narrow lines in type 1 and 2 AGNs, we found that for high-ionization or high-critical density narrow lines (e.g. OIII 5007, Balmer lines and FeII) a significant fraction of the emission arises from the inner dense part of the NLR; this inner NLR is located very close to the central engine and thus can be covered by the torus.

We present a detailed study of the X-ray and optical spectra of the luminous Seyfert galaxy Mrk 1393, which revealed variable partial obscuration of the active nucleus. The X-ray spectra obtained by XMM-Newton and Swift show moderate absorption with a column density around 3x10^21 cm^-2, consistent with a dust-reddening interpretation of the steep Balmer decrement seen in recent optical spectra. The X-ray flux in the 0.5 to 2 keV band during the XMM-Newton observation in 2005 and Swift observation in 2006 was a factor 6 brighter than that of the ROSAT All Sky Survey in 1991. In the past 4 years, the broad H\alpha line brightened by a factor of 4 accompanied by a decrease in the Balmer decrement. A comparison with literature spectra reveals variations in the dust extinction on time scales of several years, suggesting that the obscuring material is very close to the active nucleus. These observations indicate that a dust-to-gas ratio as high as the Galactic value can be present in moderately thick gas in the vicinity of the central engine within a few parsecs. We suggest that the obscuring material may be debris disrupted from the dusty torus.

In the unification scheme of Seyfert galaxies, a dusty torus blocks the continuum source and broad line region in Seyfert 2 galaxies. However it is not clear whether or not and to what extent the torus affects the narrow line spectra. In this paper, we show that Seyfert 1 and Seyfert 2 galaxies have different distributions on the [OIII]/H$\beta $ vs [NII]/H$\alpha$ diagram (BPT diagram) for narrow lines. Seyfert 2 galaxies display a clear left boundary on the BPT diagram and only 7.3% of them lie on the left. By contrast, Seyfert 1 galaxies do not show such a cutoff and 33.0% of them stand on the left side of the boundary. Among Seyfert 1 galaxies, the distribution varies with the extinction to broad lines. As the extinction increases, the distribution on BPT diagram moves to larger [NII]/H$\alpha$ value. We interpret this as an evidence for the obscuration of inner dense narrow line region by the dusty torus. We also demonstrate that the [OIII] and broad line luminosity correlation depends on the extinction of broad lines in the way that high extinction objects have lower uncorrected [OIII] luminosities, suggesting that [OIII] is partially obscured in these objects. Therefore, using [OIII] as an indicator for the nuclear luminosity will systematically under-estimate the nuclear luminosity of Seyfert 2 galaxies.

It has been demonstrated that active galactic nuclei are powered by gas accretion onto supermassive black holes located at their centres. The paradigm that the nuclei of inactive galaxies are also occupied by black holes was predicted long ago by theory. In the last decade, this conjecture was confirmed by the discovery of giant-amplitude, non-recurrent X-ray flares from such inactive galaxies and explained in terms of outburst radiation from stars tidally disrupted by a dormant supermassive black hole at the nuclei of those galaxies. Due to the scarcity of detected tidal disruption events, the confirmation and follow-up of each new candidate is needed to strengthen the theory through observational data, as well as to shed new light on the characteristics of this type of events. Two tidal disruption candidates have been detected with XMM-Newton during slew observations. Optical and X-ray follow-up, post-outburst observations were performed on these highly variable objects in order to further study their classification and temporal evolution. We show that the detected low-state X-ray emission for these two candidates has properties such that it must still be related to the flare. The X-ray luminosity of the objects decreases according to theoretical predictions for tidal disruption events. At present, optical spectra of the sources do not present any evident signature of the disruption event. In addition, the tidal disruption rate as derived from the XMM-Newton slew survey has been computed and agrees with previous studies.

(abridged) We present a comprehensive study of a sample of 23 genuine radio-loud NLS1 galaxies which have the radio-loudness parameters greater than 100. The radio sources of the sample are ubiquitously compact. A significant fraction of these objects show interesting radio to X-ray properties that are unusual to most of the previously known radio-loud NLS1 AGN, but are reminiscent of blazars. These include flat radio spectra, large amplitude flux and spectral variability, compact VLBI cores, very high brightness temperatures derived from variability, enhanced optical emission in excess of the normal ionising continuum, flat X-ray spectra, and blazar-like SEDs. We interpret them as evidence for the postulated blazar nature of these very radio-loud NLS1 AGN, which might possess at least moderately relativistic jets. Intrinsically, some of the objects have relatively low radio power and would have been classified as radio-intermediate AGN. The black hole masses are estimated to be within 10^{6-8}Msun, and the inferred Eddington ratios are around unity. The results imply that radio-loud AGN may be powered by black holes with moderate masses (10^{6-7}Msun) accreting at high rates. We find that a significant fraction of the objects, despite having strong emission lines, resemble high-energy peaked BL Lacs (HBL) in their SED. Given the peculiarities of blazar-like NLS1 galaxies, questions arise as to whether they are plain downsizing extensions of normal radio-loud AGN, or whether they form a previously unrecognised population.

We present SDSSJ092712.65+294344.0 as the best candidate to date for a recoiling supermassive black hole (SMBH). SDSSJ0927+2943 shows an exceptional optical emission-line spectrum with two sets of emission lines: one set of very narrow emission lines, and a second set of broad Balmer and broad high-ionization forbidden lines which are blueshifted by 2650 km\s relative to the set of narrow emission lines. This observation is most naturally explained if the SMBH was ejected from the core of the galaxy, carrying with it the broad-line gas while leaving behind the bulk of the narrow-line gas. We show that the observed properties of SDSSJ0927+2943 are consistent with predictions and expectations from recent numerical relativity simulations which demonstrate that SMBHs can receive kicks up to several thousand km\s due to anisotropic emission of gravitational waves during the coalescence of a binary. Our detection of a strong candidate for a rapidly recoiling SMBH implies that kicks large enough to remove SMBHs completely from their host galaxies do occur, with important implications for models of black hole and galaxy assembly at the epoch of structure formation, and for recoil models.

We report the discovery of superstrong, fading, high-ionization iron line emission in the galaxy SDSSJ095209.56+214313.3 (SDSSJ0952+2143 hereafter), which must have been caused by an X-ray outburst of large amplitude. SDSSJ0952+2143 is unique in its strong multiwavelength variability; such a broadband emission-line and continuum response has not been observed before. The strong iron line emission is accompanied by unusual Balmer line emission with a broad base, narrow core and double-peaked narrow horns, and strong HeII emission. These lines, while strong in the SDSS spectrum taken in 2005, have faded away significantly in new spectra taken in December 2007. Comparison of SDSS, 2MASS, GALEX and follow-up GROND photometry reveals variability in the NUV, optical and NIR band. Taken together, these unusual observations can be explained by a giant outburst in the EUV--X-ray band, detected even in the optical and NIR. The intense and variable iron, Helium and Balmer lines represent the ``light echo'' of the flare, as it traveled through circumnuclear material. The outburst may have been caused by the tidal disruption of a star by a supermassive black hole. Spectroscopic surveys such as SDSS are well suited to detect emission-line light echoes of such rare flare events. Reverberation-mapping of these light echoes can then be used as a new and efficient probe of the physical conditions in the circumnuclear material in non-active or active galaxies.

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.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

We obtain an exact solution of the supergravity equations of motion in which the four-dimensional observed universe is one of a number of colliding D3-branes in a Calabi-Yau background. The collision results in the ten-dimensional spacetime splitting into disconnected regions, bounded by curvature singularities. However, near the D3-branes the metric remains static during and after the collision. We also obtain a general class of solutions representing $p$-brane collisions in arbitrary dimensions, including one in which the universe ends with the mutual annihilation of a positive-tension and negative-tension 3-brane.

The cosmological evolution in Nonlinear Born-Infeld(hereafter NLBI) scalar field theory with negative potentials was investigated. The cosmological solutions in some important evolutive epoches were obtained. The different evolutional behaviors between NLBI and linear(canonical) scalar field theory have been presented. A notable characteristic is that NLBI scalar field behaves as ordinary matter nearly the singularity while the linear scalar field behaves as "stiff" matter. We find that in order to accommodate current observational accelerating expanding universe the value of potential parameters $|m|$ and $|V_0|$ must have an {\it upper bound}. We compare different cosmological evolutions for different potential parameters $m, V_0$.

We present the scalar-tensor gravitational theory with an exponential potential in which pauli metric is regarded as the physical space-time metric. We show that it is essentially equivalent to coupled quintessence(CQ) model. However for baryotropic fluid being radiation there are in fact no coupling between dilatonic scalar field and radiation. We present the critical points for baryotropic fluid and investigate the properties of critical points when the baryotropic matter is specified to ordinary matter. It is possible for all the critical points to be attractors as long as the parameters $\lambda$ and $\beta$ satisfy certain conditions. To demonstrate the attractor behaviors of these critical points, We numerically plot the phase plane for each critical point. Finally with the bound on $\beta$ from the observation and the fact that our universe is undergoing an accelerating expansion, we conclude that present accelerating expansion is not the eventual stage of universe. Moreover, we numerically describe the evolution of the density parameters $\Omega$ and the decelerating factor $q$, and computer the present values of some cosmological parameters, which are consistent with current observational data.

The purpose of this work is to search the Yangbajing neutron monitor data obtained between 1998 October and 2000 June for solar neutrons associated with solar flares. Using the onset times of 166 BATSE-detected flares with the GOES peak flux (1 -- 8 \AA) higher than $1.0 \times 10^{-5}$ $\mathrm{Wm^{-2}}$, we prepare for each flare a light curve of the Yangbajing neutron monitor, spanning $\pm$ 1.5 hours from the BATSE onset time. Based on the light curves, a systematic search for solar neutrons in energies above 100 MeV from the 166 flares was performed. No statistically significant signals due to solar neutrons were found in the present work. Therefore, we put upper limits on the $>$ 100 MeV solar-neutron flux for 18 events consisting of 2 X and 16 M class flares. The calculation assumed a power-law shaped neutron energy spectrum and three types of neutron emission profiles at the Sun. Compared with the other positive neutron detections associated with X-class flares, typical 95% confidence level upper limits for the two X-class flares are found to be comparable to the lowest and second lowest neutron fluxes at the top of the atmosphere.In addition, the upper limits for M-class flares scatter in the range of $10^{-2}$ to 1 neutrons $\mathrm{cm^{-2}s^{-1}}$. This provides the first upper limits on the solar-neutron flux from M-class solar flares, using space observatories as well as ground-based neutron monitors.

We present an analysis of the emission line property and the broad band spectral energy distribution of the ultra-luminous infrared Type II quasar Q132+058. The optical and ultraviolet emission lines show four distinct components: a LINER-like component at the systematic velocity, a heavily reddened HII-like component blueshifted by about 400 km/s, and two broad components blueshifted by about 400 km/s and 1900 km/s, respectively. The emission line ratios suggest that broad components are produced in dense and alpha-enriched outflows with a metalicity of 3-5 Zsun and a density of n_H ~ 10^7 cm^-3. The optical--UV continuum is dominated by stellar light and can be modeled with a young (<1 Myr) plus an intermediate age (0.5-0.8Gyr) stellar populations. The near to mid-infrared light is dominated by hot and warm dust heated by the hidden quasar. We derive a star formation rate (SFR) of 300-500 M\sun yr$^{-1}$ from the UV spectrum and far-infrared luminosity, which is two orders of magnitude larger than that indicated by reddening uncorrected [OII] luminosity.

We report the identification of 2MASX J032441.19+341045.9 (hereafter 2MASX J0324+3410) with an appealing object which shows the dual properties of both a narrow line Seyfert 1 galaxy (NLS1) and a blazar. Its optical spectrum, which has a H\beta line width about 1600 km s^-1 (FWHM), an [OIII] to H\beta line ratio ~0.12, and strong FeII emission, clearly fulfills the conventional definition of NLS1s. On the other hand, 2MASX J0324+3410 also exhibits some behavior which is characteristic of blazars, including a flat radio spectrum above 1 GHz, a compact core plus a one-sided jet structure on mas-scale at 8.4 GHz, highly variable fluxes in the radio, optical, and X-ray bands, and a possible detection of TeV gamma-ray emission. On its optical image, obtained with the HST WFPC2, the active nucleus is displaced from the center of the host galaxy, which exhibits an apparent one-armed spiral structure extended to 16 kpc. The remarkable hybrid behavior of this object presents a challenge to current models of NLS1 galaxies and $\gamma$-ray blazars.

We generally investigate the scalar field model with the lagrangian $L=F(X)-V(\phi)$, which we call it {\it General Non-Canonical Scalar Field Model}. We find that it is a special square potential(with a negative minimum) that drives the linear field solution($\phi=\phi_0t$) while in K-essence model(with the lagrangian $L=-V(\phi)F(X)$) the potential should be taken as an inverse square form. Hence their cosmological evolution are totally different. We further find that this linear field solutions are highly degenerate, and their cosmological evolutions are actually equivalent to the divergent model where its sound speed diverges. We also study the stability of the linear field solution. With a simple form of $F(X)=1-\sqrt{1-2X}$ we indicate that our model may be considered as a unified model of dark matter and dark energy. Finally we study the case when the baryotropic index $\gamma$ is constant. It shows that, unlike the K-essence, the detailed form of F(X) depends on the potential $V(\phi)$. We analyze the stability of this constant $\gamma_0$ solution and find that they are stable for $\gamma_0\leq1$. Finally we simply consider the constant c_s^2 case and get an exact solution for F(X)

We report the discovery of a dwarf Seyfert 1 active galactic nucleus (AGN) with a candidate intermediate-mass black hole hosted by the dwarf galaxy SDSS J160531.84+174826.1 at z=0.032. A broad component of the H-alpha line with FWHM=781 km/s is detected in its optical spectrum, and a bright, point-like nucleus is evident from a HST imaging observation. Non-thermal X-ray emission is also detected from the nucleus. The black hole mass, as estimated from the luminosity and width of the broad H-alpha component, is about 7x10^4 \msun. The host galaxy appears to be a disk galaxy with a boxy bulge or nuclear bar; with an absolute magnitude of M_R = -17.8 (M_B = -16.4), it is among the least luminous host galaxies ever identified for a Seyfert 1.

We have conducted a systematic search of AGNs with IMBHs from the SDSS DR4. As results we found 245 candidates of broad-line AGN with M_{BH}<10^6 \Msun estimated from the luminosity and width of the broad Halpha component. Compared to the pioneer Greene & Ho (2004) sample of 19 IMBH AGNs, our sample has improved in covering a larger range of the Eddington ratio, as well as black hole mass and redshift, taking the advantage of our AGN-galaxy spectral decomposition algorithm. Among these, thirty-six have L_{bol}/L_{Edd} < 0.1, hinting that a significant fraction of IMBHs might exist with weak or no nuclear activity.

This proceedings contribution summarizes first results obtained from a systematic study of 2011 NLS1 galaxies, as presented in Zhou et al. (2006, ApJS 166, 128). The sample was compiled by examining the spectral parameters of galaxies and QSOs derived from SDSS DR3 data. We discuss some preliminary results on the statistic properties of the sample, such as the fraction of NLS1, the properties of broad and narrow emission lines, and emission in other wavebands. The black hole mass--velocity dispersion relation for NLS1s was re-examined using the velocity dispersion values estimated from the stellar absorption spectra of the host galaxies. Preliminary result from an X-ray study for a small subset using data obtained by XMM-Newton is briefly discussed.

In this paper, we regard dilaton in Weyl-scaled induced gravitational theory as a coupled quintessence. Based on this consideration, we investigate the dilaton coupled quintessence(DCQ) model in $\omega-\omega'$ plane, which is defined by the equation of state parameter for the dark energy and its derivative with respect to $N$(the logarithm of the scale factor $a$). We find the scalar field equation of motion in $\omega-\omega'$ plane, and show mathematically the property of attractor solutions which correspond to $\omega_\sigma\sim-1$, $\Omega_\sigma=1$. Finally, we find that our model is a tracking one which belongs to "freezing" type model classified in $\omega-\omega'$ plane.

Based on dilatonic dark energy model, we consider two cases: dilaton field with positive kinetic energy(coupled quintessence) and with negative kinetic energy(phantom). In the two cases, we investigate the existence of attractor solutions which correspond to an equation of state parameter $\omega=-1$ and a cosmic density parameter $\Omega_\sigma=1$. We find that the coupled term between matter and dilaton can't affect the existence of attractor solutions. In the Mexican hat potential, the attractor behaviors, the evolution of state parameter $\omega$ and cosmic density parameter $\Omega$, are shown mathematically. Finally, we show the effect of coupling term on the evolution of $X(\frac{\sigma}{\sigma_0})$ and $Y(\frac{\dot{\sigma}}{\sigma^2_0})$ with respect to $N(lna)$ numerically.

This is the first paper in a series dedicated to the study of the emission-line and continuum properties of narrow line Seyfert 1 galaxies (NLS1s). We carried out a systematic search for NLS1s from objects assigned as ``QSOs'' or ``galaxies'' in the spectroscopic sample of the Sloan Digital Sky Survey Data Release 3 (SDSS DR3) by a careful modeling of their emission lines and continua. The result is a uniform sample comprising $\sim 2~000$ NLS1. This sample dramatically increases the number of known NLS1 by a factor of $\sim 10$ over previous compilations. This paper presents the parameters of the prominent emission lines and continua, which were measured accurately with typical uncertainties $<10%$. Taking advantage of such an unprecedented large and uniform sample with accurately measured spectral parameters, we carried out various statistical analysis, some of which were only possible for the first time (abbreviated).

We study the Non-Linear Born-Infeld(NLBI) scalar field model and quintessence model with two different potentials($V(\phi)=-s\phi$ and ${1/2}m^2\phi^2$). We investigate the differences between those two models. We explore the equation of state parameter w and the evolution of scale factor $a(t)$ in both NLBI scalar field and quintessence model. The present age of universe and the transition redshift are also obtained. We use the Gold dataset of 157 SN-Ia to constrain the parameters of the two models. All the results show that NLBI model is slightly superior to quintessence model.

In this paper, we employe a new statistical analysis technique, Ensemble Learning for Independent Component Analysis (EL-ICA), on the synthetic galaxy spectra from a newly released high resolution evolutionary model by Bruzual & Charlot. We find that EL-ICA can sufficiently compress the synthetic galaxy spectral library to 6 non-negative Independent Components (ICs), which are good templates to model huge amount of normal galaxy spectra, such as the galaxy spectra in the Sloan Digital Sky Survey (SDSS). Important spectral parameters, such as starlight reddening, stellar velocity dispersion, stellar mass and star formation histories, can be given simultaneously by the fit. Extensive tests show that the fit and the derived parameters are reliable for galaxy spectra with the typical quality of the SDSS.

We considered the phantom cosmology with a lagrangian $\displaystyle L=\frac{1}{\eta}[1-\sqrt{1+\eta g^{\mu\nu}\phi_{, \mu}\phi_{, \nu}}]-u(\phi)$, which is original from the nonlinear Born-Infeld type scalar field with the lagrangian $\displaystyle L=\frac{1}{\eta}[1-\sqrt{1-\eta g^{\mu\nu}\phi_{, \mu}\phi_{, \nu}}]-u(\phi)$. This cosmological model can explain the accelerated expansion of the universe with the equation of state parameter $w\leq-1$. We get a sufficient condition for a arbitrary potential to admit a late time attractor solution: the value of potential $u(X_c)$ at the critical point $(X_c,0)$ should be maximum and large than zero. We study a specific potential with the form of $u(\phi)=V_0(1+\frac{\phi}{\phi_0})e^{(-\frac{\phi}{\phi_0})}$ via phase plane analysis and compute the cosmological evolution by numerical analysis in detail. The result shows that the phantom field survive till today (to account for the observed late time accelerated expansion) without interfering with the nucleosynthesis of the standard model(the density parameter $\Omega_{\phi}\simeq10^{-12}$ at the equipartition epoch), and also avoid the future collapse of the universe.

We construct exact time-dependent solutions of the supergravity equations of motion in which two initially non-singular branes, one with positive and the other with negative tension, move together and annihilate each other in an all-enveloping spacetime singularity. Among our solutions are the Horava-Witten solution of heterotic M-theory and a Randall-Sundrum I type solution, both of which are supersymmetric, i.e. BPS, in the time-independent case. In the absence of branes our solutions are of Kasner type, and the source of instability may ascribed to a failure to stabilise some of the modulus fields of the compactification. It also raises questions about the viability of models based on some sorts of negative tension brane.

A recent paper on the large-scale structure of the Universe presented evidence for a rectangular three-dimensional lattice of galaxy superclusters and voids, with lattice spacing ~120 Mpc and called for some ``hitherto unknown process'' to explain it. Here we report that a rectangular three-dimensional lattice of intersecting domain walls, with arbitrary spacing, emerges naturally as a classical solution of M-theory.