The integrated Sachs-Wolfe (ISW) effect is caused by the decay of cosmological gravitational potential, and is therefore a unique probe of dark energy. However, its robust detection is still problematic. Various tensions between different data sets, different large scale structure (LSS) tracers, and between data and the $\Lambda$CDM theory prediction, exist. We propose a novel method of ISW measurement by cross correlating CMB and the LSS traced by "low-density-position" (LDP, \citet{2019ApJ...874....7D}). It isolates the ISW effect generated by low-density regions of the universe, but insensitive to selection effects associated with voids. We apply it to the DR8 galaxy catalogue of the DESI Legacy imaging surveys, and obtain the LDPs at $z\leq 0.6$ over $\sim$ 20000 $deg^2$ sky coverage. We then cross correlate with the Planck temperature map, and detect the ISW effect at $3.2\sigma$. We further compare the measurement with numerical simulations of the concordance $\Lambda$CDM cosmology, and find the ISW amplitude parameter $A_{ISW}=1.14\pm0.38$ when we adopt a LDP definition radius $R_s=3^{'}$, fully consistent with the prediction of the standard $\Lambda$CDM cosmology ($A_{ISW}=1$). This agreement with $\Lambda$CDM cosmology holds for all the galaxy samples and $R_s$ that we have investigated. Furthermore, the S/N is comparable to that of galaxy ISW measurement. These results demonstrate the LDP method as a competitive alternative to existing ISW measurement methods, and provide independent checks to existing tensions.

We perform two-dimensional hydrodynamic simulations to study the non-rotating gas flow in the region from 1-1000 parsec. This region connects the galaxy scale and the accretion disk scale. The gas is irradiated by the photons generated by the accretion disk system near the central black hole with $10^8$ solar mass. We assume that the luminosity of the central accretion system is constant and equals to $0.3$ Eddington luminosity. Gas with density higher than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally unstable. Two phase, cold and hot, gas flow will form due to thermal instability. We calculated the ratio of cold gas accretion rate to hot gas accretion rate. This ratio is a function of gas density and generally increases with the increase of gas density. The maximum value of this ratio is $16$, when gas density is $10^{-21} {\rm g \ cm^{-3}} $. Gas with density lower than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally stable and the gas flow is one-phase. The applications of the results are briefly discussed.

Filaments are very common physical phenomena on the Sun and are often taken as important proxies of solar magnetic activities. The study of filaments has become a hot topic in the space weather research. For a more comprehensive understanding of filaments, especially for an understanding of solar activities of multiple solar cycles, it is necessary to perform a combined multifeature analysis by constructing a data set of multiple solar cycle data. To achieve this goal, we constructed a centennial data set that covers the H$\alpha$ data from five observatories around the world. During the data set construction, we encountered varieties of problems, such as data fusion, accurate determination of the solar edge, classifying data by quality, dynamic threshold, and so on, which arose mainly due to multiple sources and a large time span of data. But fortunately, these problems were well solved. The data set includes seven types of data products and eight types of feature parameters with which we can implement the functions of data searching and statistical analyses. It has the characteristics of better continuity and highly complementary to space observation data, especially in the wavelengths not covered by space observations, and covers many solar cycles (including more than 60 yr of high-cadence data). We expect that this new comprehensive data set as well as the tools will help researchers to significantly speed up their search for features or events of interest, for either statistical or case study purposes, and possibly help them get a better and more comprehensive understanding of solar filament mechanisms.

The 0.5\,m Antarctic Survey Telescopes (AST3) were designed for time-domain optical/infrared astronomy. They are located in Dome~A, Antarctica, where they can take advantage of the continuous dark time during winter. Since the site is unattended in winter, everything for the operation, from observing to data reduction, had to be fully automated. Here, we present a brief overview of the AST3 project and some of its unique characteristics due to its location in Antarctica. We summarise the various components of the survey, including the customized hardware and software, that make complete automation possible.

The star formation rate function (SFRF) and specific star formation rate function (sSFRF) from the observation are impacted by the Eddington bias, due to the uncertainties on the estimated SFR. We develop a novel method to correct the Eddington bias and obtained the intrinsic SFRF and sSFRF from the Sloan Digital Sky Survey Data Release 7. The intrinsic SFRF is in good agreement with measurements from previous data in the literature that relied on UV SFRs but its high star-forming end is slightly lower than those IR and radio tracers. We demonstrate that the intrinsic sSFRF from SDSS has a bi-modal form with the one peak found at ${\rm sSFR \sim 10^{-9.7} yr^{-1}}$ representing the star-forming objects while the other peak is found at ${\rm sSFR \sim 10^{-12} yr^{-1}}$ representing the quenched population. Furthermore, we compare our observations with the predictions from the IllustrisTNG and Illustris simulations and show that the ``TNG'' model performs much better than its predecessor. However, we show that the simulated SFRF and cosmic star formation density (CSFRD) of TNG simulations are highly dependent on resolution, reflecting the limitations of the model and today state-of-the-art simulations. We demonstrate that the bi-modal, two peaked sSFRF implied by the SDSS observations does not appear in TNG regardless of the adopted box-size or resolution. This tension reflects the need for inclusion of an additional efficient quenching mechanism to the TNG model.

We investigate a new mechanism to create large curvature perturbations on small scales due to parameter resonance in a single-field inflationary model with a small periodic structure upon the potential. After reentering the horizon, the amplified curvature perturbations can lead to observable primordial black holes as well as stochastic gravitational waves. The mass of primordial black holes and frequency of the induced gravitational waves depend on the model parameters. The resulted primordial black hole could constitute all dark matter or a fraction of dark matter in the universe, and corresponding stochastic gravitational waves fall in the frequency band measurable for the pulsar timing array and the space-based gravitational wave detectors.

Using the galaxy catalog built from ELUCID N-body simulation and the semi-analytical galaxy formation model, we have built a mock HI intensity mapping map. We have implemented the Finger-of-God (FoG) effect in the map by considering the galaxy HI gas velocity dispersion. By comparing the HI power spectrum in the redshift space with the measurement from IllustrisTNG simulation, we have found that such FoG effect can explain the discrepancy between current mock map built from N-body simulation and Illustris TNG simulation. Then we built a parameter-free FoG model and a shot-noise model to calculate the HI power spectrum. We found that our model can accurately fit both the monopole and quadrupole moments of the HI matter power spectrum. Our method of building the mock HI intensity map and the parameter-free FoG model will be widely useful for the up-coming 21cm intensity mapping experiments, such as CHIME, Tianlai, BINGO, FAST and SKA. It is also crucial for us to study the non-linear effects in 21cm intensity mapping.

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

We analyze ground-based chromospheric data acquired at a high temporal cadence of 2 s in wings of the H$\alpha$ spectral line using Goode Solar Telescope (GST) operating at the Big Bear Solar Observatory. We inspected a 30 minute long H$\alpha$-0.08~nm data set to find that rapid blue-shifted H$\alpha$ excursions (RBEs), which are a cool component of type II spicules, experience very rapid morphological changes on the time scales of the order of 1 second. Unlike typical reconnection jets, RBEs very frequently appear \textit{in situ} without any clear evidence of H$\alpha$ material being injected from below. Their evolution includes inverted "Y", "V", "N", and parallel splitting (doubling) patterns as well as sudden formation of a diffuse region followed by branching. We also find that the same feature may undergo several splitting episodes within about 1 min time interval.

We constrain the rest-frame FUV (1546\AA), NUV (2345\AA) and U-band (3690\AA) luminosity functions (LFs) and luminosity densities (LDs) with unprecedented precision from $z\sim0.2$ to $z\sim3$ (FUV, NUV) and $z\sim2$ (U-band). Our sample of over 4.3 million galaxies, selected from the CFHT Large Area $U$-band Deep Survey (CLAUDS) and HyperSuprime-Cam Subaru Strategic Program (HSC-SSP) data lets us probe the very faint regime (down to $M_\mathrm{FUV},M_\mathrm{NUV},M_\mathrm{U} \simeq -15$ at low redshift) while simultaneously detecting very rare galaxies at the bright end down to comoving densities $<10^{-5}$ Mpc$^{-3}$. Our FUV and NUV LFs are well fitted by single Schechter functions, with faint-end slopes that are very stable up to $z\sim2$. We confirm, but self-consistently and with much better precision than previous studies, that the LDs at all three wavelengths increase rapidly with lookback time to $z\sim1$, and then much more slowly at $1<z<2$--$3$. Evolution of the FUV and NUV LFs and LDs at $z<1$ is driven almost entirely by the fading of the characteristic magnitude, $M^\star_{UV}$, while at $z>1$ it is due to the evolution of both $M^\star_{UV}$ and the characteristic number density $\phi^\star_{UV}$. In contrast, the U-band LF has an excess of faint galaxies and is fitted with a double-Schechter form; $M^\star_\mathrm{U}$, both $\phi^\star_\mathrm{U}$ components, and the bright-end slope evolve throughout $0.2<z<2$, while the faint-end slope is constant over at least the measurable $0.05<z<0.6$. We present tables of our Schechter parameters and LD measurements that can be used for testing theoretical galaxy evolution models and forecasting future observations.

The 21~cm intensity mapping experiments promise to measure the baryon acoustic oscillations (BAO) through the mapping of H{\small I} gas at the post-reionization epoch. Its success depends on the assumption that H{\small I} gas traces the matter density distribution {\it linearly} on large scales. However, both nonlinear halo clustering and nonlinear effects modulating H{\small I} gas in halos may spoil this. We employ three approaches to generate the mock H{\small I} density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of H{\small I} linearity is valid at the scale corresponding to the first BAO peak, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and H{\small I} content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the H{\small I} bias is scale-independent down to small scales. We also find that the linear H{\small I} bias scales approximately linearly with redshift for $z\le 3$.

We have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates from the helical motion jet in a supermassive binary black hole (SMBBH) system undergoing major merger, we estimate the primary black hole mass M~7.7*10^{9} M_sun. To explore the origin of the gamma-ray, we investigated the optical-gamma-ray correlations using discrete correlation function (DCF) method, and found that the correlation between the two bands is very significant. This strong correlation tends to imply lepton self-synchro-Compton (LSSC) model to produce the gamma-ray.

We test the Emergent Gravity(EG) theory using the galaxy-galaxy lensing technique based on SDSS DR7 data. In the EG scenario, we do not expect color dependence of the galaxy sample in the 'apparent dark matter' predicted by EG, which is exerted only by the baryonic mass. If the baryonic mass is similar, then the predicted lensing profiles from the baryonic mass should be similar according to the EG, regardless of the color of the galaxy sample. We use the stellar mass of the galaxy as a proxy of its baryonic mass. We divide our galaxy sample into 5 stellar mass bins, and further classify them as red and blue subsamples in each stellar mass bin. If we set halo mass and concentration as free parameters, $\Lambda$CDM is favored by our data in terms of the reduced $\chi^2$ while EG fails to explain the color dependence of ESDs from the galaxy-galaxy lensing measurement.

We continue our empirical study of the emission line flux originating in the cool ($T\sim10^4$ K) gas that populates the halos of galaxies and their environments. Specifically, we present results obtained for a sample of galaxy pairs with a range of projected separations, {\bf $10 < {S_p/\rm kpc} < 200$}, and mass ratios $<$ 1:5, intersected by 5,443 SDSS lines of sight at projected radii of 10 to 50 kpc from either or both of the two galaxies. We find significant enhancement in H$\alpha$ emission and a moderate enhancement in [N {\small II}]6583 emission for low mass pairs (mean stellar mass per galaxy, $\overline{\rm M}_*, <10^{10.4} {\rm M}_\odot$) relative to the results from a control sample. This enhanced H$\alpha$ emission comes almost entirely from sight lines located between the galaxies, consistent with a short-term, interaction-driven origin for the enhancement. We find no enhancement in H$\alpha$ emission, but significant enhancement in [N {\small II}]6583 emission for high mass ($\overline{\rm M}_* >10^{10.4}{\rm M}_\odot$) pairs. Furthermore, we find a dependence of the emission line properties on the galaxy pair mass ratio such that those with a mass ratio below 1:2.5 have enhanced [N {\small II}]6583 and those with a mass ratio between 1:2.5 and 1:5 do not. In all cases, departures from the control sample are only detected for close pairs ($S_p <$ 100 kpc). Attributing an elevated [N {\small II}]6583/H$\alpha$ ratio to shocks, we infer that shocks play a role in determining the CGM properties for close pairs that are among the more massive and have mass ratios closer to 1:1.

Here we confirm the short periodic variability of a subparsec supermassive binary black hole (SMBBH) candidate Mrk 231 in the extended optical photometric data set collected by the Catalina Real-Time Transient Survey (CRTS) and All-Sky Automated Survey for Supernovae (ASAS-SN). Using the Lomb-Scargle periodogram and 2DHybrid method, we detected the significant periodicity of ~ 1.1 yr beyond a damped random walk model in the CRTS+ASAS-SN optical data set. Mrk 231 has been previously proposed as a SMBBH candidate with a highly unequal mass ratio (q~ 0.03), very tight mutual separation of ~590 AU, and an orbital period of ~1.2 yr. Hence, our result further supports, even though not prove, the intriguing hypothesis that SMBBHs with low mass ratios may be more common than close-equal mass SMBBHs. This result, however, was obtained from the contribution of CRTS data with limited sampling cadence and photometric accuracy, and further monitoring of Mrk 231 is crucial to confirm the periodicity.

We use the ELUCID N-Body simulation to study the relation of structural properties of dark matter halos to their assembly history and environment. We find that the complexity of halo assembly histories can be well described by a small number of principal components. Using decision trees built with the random ensemble method, we find that the correlation between halo concentration and assembly history is tight: more than 60% of the variance in halo concentration can be explained by assembly history alone. Halo shape and spin parameters are found to depend not only on assembly history, but also on environment and initial condition, and these dependencies are not orthogonal. The structures of more massive halos depend more strongly on higher-order components of the assembly history. The relations between halo intrinsic properties and environment are weak compared to their variances, with the anisotropy of local tidal fields having the strongest correlation with halo properties. Our method of dimension reduction and regression can help simplify the characterization of the halo population and clarify the degeneracy among halo properties.

Chiral symmetry is maximally violated in weak interactions, and such microscopic asymmetries in the early Universe might leave observable imprints on astrophysical scales without violating the cosmological principle. In this Letter, we propose a helicity measurement to detect primordial chiral violation. We point out that observations of halo-galaxy angular momentum directions (spins), which are frozen in during the galaxy formation process, provide a fossil chiral observable. From the clustering mode of large scale structure of the Universe, we construct a spin mode in Lagrangian space and show in simulations that it is a good probe of halo-galaxy spins. In standard model, a strong symmetric correlation between the left and right helical components of this spin mode and galaxy spins is expected. Measurements of these correlations will be sensitive to chiral breaking, providing a direct test of chiral symmetry breaking in the early Universe.

It may be determined by non-parametric method if the dark energy evolves with time. We propose a method of combining PCA and biased estimation on the basis of ridge regression analysis to reconstruct parameters, meanwhile we present an interesting principal component selection criterion to avoid the arbitrariness of principal component selections, and use numerical integral by Lagrange interpolation to linearize the luminosity distance integral formula in nearly flat space to avoid instability of derivative for functional data. We get the preliminary test results that shows if $\Delta \overline w (z) = \overline {\left| {1 + w(z)} \right|} <=0.05$ included $w(z) = -1$, the probability of making a type I error for ${w_{recon}} \ne -1$ is almost zero ($1\% $) in the test; otherwise, if $\Delta \overline w (z) = \overline {\left| {1 + w(z)} \right|} > 0.05$, the probability of making a type I error for ${w_{recon}} = -1$ is not more than $10\% $. Finally, we use JLA sample to reconstruct $w(z)$, and the results reject ${\rm{w(z) }} \ne {\rm{ - 1}}$, which is agreement with $\Lambda CDM$ model.

A good compromise between the resistive model and the PIC model is Aristotelian electrodynamics, which can include the back-reaction of the radiative photons onto particle motion and allow for a local dissipation where the force-free condition is violated. We study the dissipative pulsar magnetosphere with Aristotelian electrodynamics where particle acceleration is fully balanced by radiation. The expression for the current density is defined by introducing a pair multiplicity. The 3D structure of the pulsar magnetosphere is then presented by solving the time-dependent Maxwell equations using a pseudo-spectral algorithm. It is found that the dissipative magnetosphere approaches the force-free solution and the dissipative region is more restricted to the current sheet outside the light-cylinder (LC) as the pair multiplicity increases. The spatial extension of the dissipative region is self-consistently controlled by the pair multiplicity. Our simulations show the high magnetospheric dissipation outside the LC for the low pair multiplicity.

In order to measure the Universal Time (UT1) in real time, National Time Service Center (NTSC) has built a VGOS-like (VLBI Global Observing System) broadband VLBI network, which includes three 13-m radio telescopes located in Jilin, Sanya and Kashi, and a data analysis center in Xi'an. Each station is equipped with a highly stable hydrogen atomic clock and a self-developed VLBI backend, and is co-located with two GPS receivers. This VGOS-like VLBI network may play an important role in improving the Chinese broadband VLBI technology and making valuable contributions to domestic VLBI measurements of UT1. In this paper, we introduce the specifications of this VLBI network, and present the UT1 measurements at C-band conducted in 2018 using the Jilin-Kashi baseline of this network. The comparisons between our UT1 estimates and those provided by IERS suggest that the NTSC VLBI network is capable to determine UT1 accurate at the level of 58.8 microseconds.

In this paper, we present a novel method to test the Einstein's Equivalence Principle (EEP) using (simultaneous) multi-wavelength radio observations of polarized blazars. We analyze simultaneous multi-wavelength polarization observations of 3C 279 at 22, 43, and 86 GHz obtained by two antennas of the Korean VLBI Network. We obtained 15 groups of polarization data, and applied the Metropolis-Hastings Markov Chain (MHMC) to simulate the parameters when considering the EEP effect and the simplest form of Faraday rotation (single external Faraday screen). The final results show the constraint of the parameterized post-Newtonian (PPN) parameter $\gamma$ discrepancy as $\Delta \gamma_{p} = (1.91\pm0.34)\times10^{-20}$. However, the single external Faraday screen is an oversimplification for blazars because there are numerous observations show complex Faraday rotation behavior for blazars due to internal/external Faraday dispersion, beam depolarization, etc. The value $\Delta \gamma_{p}$ results of this paper can only be considered as upper limits. Only if all other effects are revealed and considered, should the result be taken as a direct measurement of the violation of the EEP.

A sunspot is an ideal waveguide for a variety of magnetohydrodynamic waves, which carry a significant amount of energy to the upper atmosphere and could be used as a tool to probe magnetic and thermal structure of a sunspot. In this study, we used the New Vacuum Solar Telescope and took high-resolution image sequences simultaneously in both TiO (7058$\pm$10 \AA) and H$_\alpha$ (6562$\pm$2.5 \AA) bandpasses. We extracted the area and total emission intensity variations of sunspot umbra and analyzed the signals with synchrosqueezing transform. We found that the area and emission intensity varied with both three and five minute periodicity. Moreover, the area and intensity oscillated in phase with each other, this fact hold in both TiO and H$_\alpha$ data. We interpret this oscillatory signal as propagating slow sausage wave. The propagation speed is estimated at about 8 km$\cdot$s$^{-1}$. We infer that this sunspot's umbra could have temperature as low as 2800--3500 K.

For filter-based magnetographs, the linear calibration method under the weak-field assumption is usually adopted; this leads to magnetic saturation effect in the regions with strong magnetic field. This article explores a new method to overcome the above disadvantage using a multilayer perceptron network, which we call MagMLP, based on a back-propagation algorithm with one input layer, five hidden layers, and one output layer. We use the data from the \textit{Spectropolarimeter} (SP) on board \textit{Hinode} to simulate single-wavelength observations for the model training, and take into account the influence of the Doppler velocity field and the filling factor. The training results show that the linear fitting coefficient (LFC) of the transverse field reaches above 0.91, and that of the longitudinal field is above 0.98. The generalization of the models is good because the corresponding LFCs are above 0.9 for the test subsets. Compared with the linear calibration method, the MagMLP is much more effective on dealing with the magnetic saturation effect. Analyzing an active region, the results of the linear calibration present an evident magnetic saturation effect in the umbra regions; the corresponding systematic error reaches values greater than 1000 G in most areas, or even exceeds 2000 G at some pixels. However, the results of MagMLP at these locations are very close to the inversion results, and the systematic errors are basically within 300 G. In addition, we find that there are many "bright spots" and "dark spots" on the inclination angle images from the inversion results of \textit{Hinode}/SP with values of 180 and 0 degrees, respectively, where the inversion is not reliable and does not produce a good result; the MagMLP handles these points well.

There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \mu m}$, MIPS - 24, 70 and 160 ${\rm \mu m}$ or SPIRE - 250, 350, 500 ${\rm \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.

Numerical simulations of radiative two-temperature hot accretion flows (HAFs) around Neutron stars (NSs) are performed. We assume that all of the energy carried by the HAF around a NS will be thermalized and radiated out at the surface of the NS. The thermal photons will propagate outwards radially and cool the HAF vis Comptonization. We define $\dot m$ as the mass accretion rate at the surface of the central object in unit of Eddington accretion rate ($\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ being Eddington luminosity and speed of light, respectively). When $\dot m$ is lower than $\sim 10^{-4}$, the cooling of the HAF is not important and outflows are very strong. When $\dot m > \sim 10^{-3}$, cooling becomes important and outflows are significantly weak. In the range $10^{-4} < \dot m < 10^{-3}$, the HAFs transients from a strong outflow phase to a very weak outflow phase with increase of $\dot m$. The properties of the HAF around a NS are also compared to those of the HAF around a BH. We find that with a similar $\dot m$, the dynamical properties of the HAF around a NS are quite similar as those of the HAF around a BH. However, the emitted spectrum of a HAF around a NS can be quite different from that of a HAF around a BH due to the presence of a thermal soft X-ray component coming from the surface of the NS.

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

Superhydrogenated polycyclic aromatic hydrocarbons (PAHs) may be present in H-rich and ultraviolet-poor benign regions. The addition of excess H atoms to PAHs converts the aromatic bonds into aliphatic bonds, the strongest of which falls near 3.4 $\mu$m. Therefore, superhydrogenated PAHs are often hypothesized as a carrier of the 3.4 $\mu$m emission feature which typically accompanies the stronger 3.3 $\mu$m aromatic C--H stretching feature. To assess this hypothesis, we use density function theory to compute the IR vibrational spectra of superhydrogenated PAHs and their ions of various sizes (ranging from benzene, naphthalene to perylene and coronene) and of various degrees of hydrogenation (ranging from minimal hydrogenation to heavy hydrogenation). For each molecule, we derive the intrinsic oscillator strengths of the 3.3 $\mu$m aromatic C--H stretch ($A_{3.3}$) and the 3.4 $\mu$m aliphatic C--H stretch ($A_{3.4}$). By comparing the computationally-derived mean ratio of $\langle A_{3.4}/A_{3.3}\rangle\sim1.98$ with the mean ratio of the observed intensities $\langle I_{3.4}/I_{3.3}\rangle\sim0.12$, we find that the degree of superhydrogenation --- the fraction of C atoms attached with excess H atoms --- is only $\sim2.2\%$ for neutral PAHs which predominantly emit the 3.3 and 3.4 $\mu$m features. We also determine for each molecule the intrinsic band strengths of the 6.2 $\mu$m aromatic C--C stretch ($A_{6.2}$) and the 6.85 $\mu$m aliphatic C--H deformation ($A_{6.85}$). We derive the degree of superhydrogenation from the mean ratio of the observed intensities $\langle I_{6.85}/I_{6.2}\rangle\sim0.10$ and $\langle A_{6.85}/A_{6.2}\rangle\sim1.53$ for neutrals and $\langle A_{6.85}/A_{6.2}\rangle\sim1.23$ for cations to be $\lesssim 3.1\%$ for neutrals and $\lesssim 8.6\%$ for cations. We conclude that astrophysical PAHs are primarily aromatic and are only marginally superhydrogenated.

We propose a radically different picture for X-shaped radio galaxies compared to existing models as transition objects between cold mode accreting low spinning retrograde black holes and low spinning prograde black holes. The model explains their smaller average black hole masses, their general aversion for cluster compared to isolated environments, the negligible difference in spectral index between primary and secondary jets despite a time difference in their formation, their absence among the most powerful radio quasars and radio galaxies, and their connection to the elusive FRI quasar class, among others. The key to their formation is cold gas accretion onto Schwarzschild black holes.

As demonstrated in Paper I, the quenching properties of central and satellite galaxies are quite similar as long as both stellar mass and halo mass are controlled. Here we extend the analysis to the size and bulge-to-total light ratio (B/T) of galaxies. In general central galaxies have size-stellar mass and B/T-stellar mass relations different from satellites. However, the differences are eliminated when halo mass is controlled. We also study the dependence of size and B/T on halo-centric distance and find a transitional stellar mass (M$_{*,t}$) at given halo mass (M$_h$), which is about one fifth of the mass of the central galaxies in halos of mass M$_h$. The transitional stellar masses for size, B/T and quenched fraction are similar over the whole halo mass range, suggesting a connection between the quenching of star formation and the structural evolution of galaxies. Our analysis further suggests that the classification based on the transitional stellar mass is more fundamental than the central-satellite dichotomy, and provide a more reliable way to understand the environmental effects on galaxy properties. We compare the observational results with the hydro-dynamical simulation, EAGLE and the semi-analytic model, L-GALAXIES. The EAGLE simulation successfully reproduces the similarities of size for centrals and satellites and even M$_{*,t}$, while L-GALAXIES fails to recover the observational results.

Galaxy clustering provides insightful clues to our understanding of galaxy formation and evolution, as well as the universe. The redshift assignment for the random sample is one of the key steps to measure the galaxy clustering accurately. In this paper, by virtue of the mock galaxy catalogs, we investigate the effect of two redshift assignment methods on the measurement of galaxy two-point correlation functions (hereafter 2PCFs), the Vmax method and the "shuffled" method. We found that the shuffled method significantly underestimates both of the projected 2PCFs and the two-dimensional 2PCFs in redshift space. While the Vmax method does not show any notable bias on the 2PCFs for volume-limited samples. For flux-limited samples, the bias produced by the Vmax method is less than half of the shuffled method on large scales. Therefore, we strongly recommend the Vmax method to assign redshifts to random samples in the future galaxy clustering analysis.

We continue our empirical study of the emission line flux originating in the cool ($T\sim10^4$ K) gas that populates the halos of galaxies and their environments. Specifically, we present results obtained for a sample of nearly half a million individual galaxies, groups, and clusters of galaxies, intersected by more than two million SDSS lines of sight at projected separations of up to a quarter of the virial radius. Adopting simple power law relationships between the circumgalactic (CGM) cool gas fraction and either the halo or stellar mass, we present expressions for the CGM cool gas fraction as a function of either halo or stellar mass, $f_{\rm cool}(M_h) = (0.25^{+0.07}_{-0.06}) \times (M_h/10^{12}M_\odot)^{(-0.39^{+0.06}_{-0.07})}$ or $f_{\rm cool}(M_{*}) = (0.28^{+0.07}_{-0.06}) \times (M_{\rm *}/10^{10.0}M_\odot)^{(-0.33\pm 0.06)}$. Where we can compare, our results are consistent with previous constraints from absorption line studies, our own previous emission line work, and simulations. The cool gas can be the dominant baryonic CGM component, comprising a fraction as high as $> 90\%$ of halo gaseous baryons, in low mass halos, $M_h\sim$ $10^{10.5} M_\odot$, and a minor fraction, $<$ 5\%, in groups and clusters, $M_h > 10^{14} M_\odot$.

In order to find clues to the origin of the "winged" or "X-shaped" radio galaxies (XRGs) we investigate here the parent galaxies of a large sample of 106 XRGs for optical-radio axes alignment, interstellar medium, black hole mass, and large-scale environment. For 41 of the XRGs it was possible to determine the optical major axis and the primary radio axis and the strong tendency for the two axes to be fairly close is confirmed. However, several counter-examples were also found and these could challenge the widely discussed backflow diversion model for the origin of the radio wings. Comparison with a well-defined large sample of normal FR II radio galaxies has revealed that: (i) XRGs possess slightly less massive central black holes than the normal radio galaxies (average masses being log$M_{\rm BH} \sim$ 8.81 $M_{\odot}$ and 9.07 $M_{\odot}$, respectively); (ii) a much higher fraction of XRGs ($\sim$ 80%) exhibits red mid-IR colors ($W2 - W3 > 1.5$), indicating a population of young stars and/or an enhanced dust mass, probably due to relatively recent galaxy merger(s). A comparison of the large-scale environment (i.e., within $\sim$ 1 Mpc) shows that both XRGs and FRII radio galaxies inhabit similarly poor galaxy clustering environments (medium richness being 8.94 and 11.87, respectively). Overall, the origin of XRGs seems difficult to reconcile with a single dominant physical mechanism and competing mechanisms seem prevalent.

If the black holes detected by LIGO/VIRGO are primordial black holes (PBHs) sourcing from a large primordial curvature perturbation on small scales, the corresponding induced gravitational waves (GWs) would peak at nanohertz that is detectable by the current and future observations of pulsar timing array (PTA). In this paper we show that with the mass function estimated from the merger rate of LIGO O1 and O2 events, the induced GWs from such a curvature perturbation with a Gaussian narrow peak at some small scale would be in a seemingly mild tension with current constraints from PTA. However, if the curvature perturbation is of local-type non-Gaussianity with a non-linear parameter $f_\text{NL}\gtrsim\mathcal{O}(10)$, the tension could be relieved. Nevertheless, such an induced GWs must be detectable by the Square Kilometer Array in a decade or less.

The paper presents a reliable method using deep learning to recognize solar filaments in H-alpha full-disk solar images automatically. This method cannot only identify filaments accurately but also minimize the effects of noise points of the solar images. Firstly, a raw filament dataset is set up, consisting of tens of thousands of images required for deep learning. Secondly, an automated method for solar filament identification is developed using the U-Net deep convolutional network. To test the performance of the method, a dataset with 60 pairs of manually corrected H-alpha images is employed. These images are obtained from the Big Bear Solar Observatory/Full-Disk H-alpha Patrol Telescope (BBSO/FDHA) in 2013. Cross-validation indicates that the method can efficiently identify filaments in full-disk H-alpha images.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

Prior to the availability of digital cameras, the solar observational images are typically recorded on films, and the information such as date and time were stamped in the same frames on film. It is significant to extract the time stamp information on the film so that the researchers can efficiently use the image data. This paper introduces an intelligent method for extracting time stamp information, namely, the Convolutional Neural Network (CNN), which is an algorithm in deep learning of multilayer neural network structures and can identify time stamp character in the scanned solar images. We carry out the time stamp decoding for the digitized data from the National Solar Observatory from 1963 to 2003. The experimental results show that the method is accurate and quick for this application. We finish the time stamp information extraction for more than 7 million images with the accuracy of 98\%.

In Cassini ISS (Imaging Science Subsystem) images, contour detection is often performed on disk-resolved object to accurately locate their center. Thus, the contour detection is a key problem. Traditional edge detection methods, such as Canny and Roberts, often extract the contour with too much interior details and noise. Although the deep convolutional neural network has been applied successfully in many image tasks, such as classification and object detection, it needs more time and computer resources. In the paper, a contour detection algorithm based on H-ELM (Hierarchical Extreme Learning Machine) and DenseCRF (Dense Conditional Random Field) is proposed for Cassini ISS images. The experimental results show that this algorithm's performance is better than both traditional machine learning methods such as SVM, ELM and even deep convolutional neural network. And the extracted contour is closer to the actual contour. Moreover, it can be trained and tested quickly on the general configuration of PC, so can be applied to contour detection for Cassini ISS images.

A hot corona is suggested to be above the standard thin disk. The anisotropy of hard X-ray emission in radio-quiet active galactic nuclei implies that the corona is not static and probably moves outwards like winds. We perform two-dimensional magnetohydrodynamical simulations to study the outflowing corona driven by magnetic field and radiation force. In our simulations, as the initial state and the boundary condition at the disk surface, the corona temperature is set to $10^9$ K inside 10 Schwarzschild radius ($r_{\rm s}$) while the corona temperature is set to $10^7$ K a 10 $r_{\rm s}$. We employ a weak poloidal magnetic field as the initial magnetic field. A collimated outflow and a wide-angle ordered outflow are observed in our simulations. The collimated outflow is around the rotational axis and has a bulk velocity of $\sim$0.03--0.3 $c$ ( $c$ is speed of light) at 90 $r_{\rm s}$, while their mass outflow rate is very low. The collimated outflow is a weak jet. The wide-angle ordered outflow is distributed at the middle and high latitudes and moves outwards with a velocity of $10^2$--$10^4$ km s$^{-1}$. The outflow velocity depends on the disk luminosity. The gas around the disk surface is turbulent, especially outside of 10 $r_{\rm s}$. The other properties of outflows are discussed in detail.

We perform simulations to study the effects of active galactic nuclei (AGNs) radiation and wind feedback on the properties of slowly rotating accretion flow at parsec scale. We find that when only radiative feedback is considered, outflows can be produced by the radiation pressure due to Thomson scattering. The mass flux of outflow is comparable to that of inflow. Although strong outflow is present, the luminosity of AGN can be easily super-Eddington. When wind feedback is also taken into account, the mass flux of outflow does not change much. Consequently, the luminosity of the central AGN can still be super-Eddington. However, observations show that the luminosity of most AGNs is sub-Eddington. Some other mechanisms are needed to reduce the AGNs luminosity. Although the mass outflow rate does not change much by wind feedback, other properties of outflow (density, temperature, velocity and kinetic power) can be significantly changed by wind feedback. In the presence of wind feedback, the density of outflow becomes significantly lower, temperature of outflow becomes significantly higher, velocity of outflow is increased by one order of magnitude, the kinetic power of outflow is increased by a factor of $40-100$.

Quantum electrodynamics (QED) effects may be included in physical processes of magnetar and pulsar magnetospheres with strong magnetic fields. Involving the quantum corrections, the Maxwell electrodynamics is modified to non-linear electrodynamics. In this work, we study the force-free magnetosphere in non-linear electrodynamics in a general framework. The pulsar equation describing a steady and axisymmetric magnetosphere is derived, which now admits solutions with corrections. We derive the first-order non-linear corrections to the near-zone dipole magnetosphere in some popular non-linear effective theories. The field lines of the corrected dipole tend to converge on the rotational axis so that the fields in the polar region are stronger compared to the pure dipole case.

We investigate the abundance, small-scale clustering and galaxy-galaxy lensing signal of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). To this end, we present new measurements of the redshift and stellar mass dependence of the lensing properties of the galaxy sample. We analyse to what extent models assuming the Planck18 cosmology fit to the number density and clustering can accurately predict the small-scale lensing signal. In qualitative agreement with previous BOSS studies at redshift $z \sim 0.5$ and with results from the Sloan Digital Sky Survey, we find that the expected signal at small scales ($0.1 < r_{\rm p} < 3 \, h^{-1} \mathrm{Mpc}$) is higher by $\sim 25\%$ than what is measured. Here, we show that this result is persistent over the redshift range $0.1 < z < 0.7$ and for galaxies of different stellar masses. If interpreted as evidence for cosmological parameters different from the Planck CMB findings, our results imply $S_8 = \sigma_8 \sqrt{\Omega_{\rm m} / 0.3} = 0.744 \pm 0.015$, whereas $S_8 = 0.832 \pm 0.013$ for Planck18. However, in addition to being in tension with CMB results, such a change in cosmology alone does not accurately predict the lensing amplitude at larger scales. Instead, other often neglected systematics like baryonic feedback or assembly bias are likely contributing to the small-scale lensing discrepancy. We show that either effect alone, though, is unlikely to completely resolve the tension. Ultimately, a combination of the two effects in combination with a moderate change in cosmological parameters might be needed.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

In the CDM paradigm, the halo mass function is a sensitive probe of the cosmic structure. In observations, halo mass is typically estimated from its relation with other observables. The resulting halo mass function is subject to systematic bias, such as the Eddington bias, due to the scatter or uncertainty in the observable - mass relation. Exact correction for the bias is not easy, as predictions for the observables are typically model-dependent in simulations. In this paper, we point out an interesting feature in the halo mass function of the concordence $\Lambda$CDM model: the total halo mass within each evenly-spaced logarithmic mass bin is approximately the same over a large mass range. We show that this property allows us to construct an almost bias-free halo mass function using only an observable (as a halo mass estimator) and stacked weak lensing measurements as long as the scatter between the true halo mass and the observable-inferred mass has a stable form in logarithmic units. The method is not sensitive to the form of the mass-observable relation. We test the idea using cosmological simulations, and show that the method performs very well for realistic observables.

We study the morphological transformation from late types to early types and the quenching of galaxies with the seventh Data Release (DR7) of the Sloan Digital Sky Survey (SDSS). Both early type galaxies and late type galaxies are found to have bimodal distributions on the star formation rate versus stellar mass diagram ($\lg SFR - \lg M_*$). We therefore classify them into four types: the star-forming early types (sEs), the quenched early types (qEs), the star-forming late types (sLs) and the quenched late types (qLs). We checked many parameters on various environmental scales for their potential effects on the quenching rates of late types and early types, as well as the early type fractions among star-forming galaxies and those among quenched galaxies. These parameters include: the stellar mass $M_*$, and the halo mass $M_{halo}$; the small-scale environmental parameters, such as the halo centric radius $R_p/r_{180}$ and the third nearest neighbor distances ($d_{3nn}$); the large-scale environmental parameters, specifically whether they are located in clusters, filaments, sheets, or voids. We found that the morphological transformation is mainly regulated by the stellar mass. Quenching is mainly driven by the stellar mass for more massive galaxies and by the halo mass for galaxies with smaller stellar masses. In addition, we see an overall stronger halo quenching effect in early type galaxies, which might be attributed to their lacking of cold gas or earlier accretion into the massive host halos.

We present a catalog of 290 `winged' or X-shaped radio galaxies (XRGs) extracted from the latest (17 Dec. 2014) data release of the `Very Large Array Faint Images of the Radio Sky at Twenty centimeter' (VLA FIRST survey). We have combined these radio images with their counterparts in the TIFR GMRT sky survey (TGSS) at 150 MHz (TGSS alternative data release, TGSS_ADR1), in an attempt to identify any low surface-brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 `strong' XRG candidates and 184 `probable' XRG candidates whose XRG designation must be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs presently known. Twenty-five of our 290 XRG candidates (9 `strong' and 16 `probable') are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (2 `strong' and 1 `probable'). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index $\alpha$ flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 GHz and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

We present a catalog of 290 `winged' or X-shaped radio galaxies (XRGs) extracted from the latest (17 Dec. 2014) data release of the `Very Large Array Faint Images of the Radio Sky at Twenty centimeter' (VLA FIRST survey). We have combined these radio images with their counterparts in the TIFR GMRT sky survey (TGSS) at 150 MHz (TGSS alternative data release, TGSS_ADR1), in an attempt to identify any low surface-brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 `strong' XRG candidates and 184 `probable' XRG candidates whose XRG designation must be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs presently known. Twenty-five of our 290 XRG candidates (9 `strong' and 16 `probable') are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (2 `strong' and 1 `probable'). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index $\alpha$ flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 GHz and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

Interstellar graphene could be present in the interstellar medium (ISM), resulting from the photochemical processing of polycyclic aromatic hydrocarbon (PAH) molecules and/or collisional fragmentation of graphitic particles. Indeed, by comparing the observed ultraviolet (UV) extinction and infrared (IR) emission of the diffuse ISM with that predicted for graphene, as much as $\sim$ 2% of the total interstellar carbon could have been locked up in graphene without violating the observational constraints. While the possible detection of planar C$_{24}$, a small piece of a graphene sheet, has been reported towards several Galactic and extragalactic planetary nebulae, graphene has not yet been detected in interstellar environments. In this work, we search for the characteristic IR features of C$_{24}$ at $\sim$ 6.6, 9.8, 20 $\mu m$ toward Sgr B2, a high-mass star formation region, and find a candidate target toward R.A. (J2000) = $267^{\circ}.05855$ and Decl. (J2000) = $-28^{\circ}.01479$ in Sgr B2 whose Spitzer/IRS spectra exhibit three bands peaking at $\sim$ 6.637, 9.853 and 20.050 $\mu m$ which appear to be coincident with that of C$_{24}$. Possible features of C$_{60}$ are also seen in this region. The candidate region is a warm dust environment heated by massive stars or star clusters, associated with a WISE spot (a tracer of star-formation activity), close to the HII region candidate IRAS 17450-2759, and is surrounded by seven young stellar object candidates within $\sim$ 5$^{\prime}$, suggesting that the creation and/or excitation of C$_{24}$ could be related to star formation activities.

We identify analytically a multiple-peak structure in the energy-density spectrum of induced gravitational waves (GWs) generated at second-order from a primordial scalar perturbations also with multiple($n$) peaks at small scales $k_{*i}$. The energy-density spectrum of induced GWs exhibits at most $C_{n+1}^2$ and at least $n$ peaks at wave-vectors $k_{ij}\equiv(k_{*i}+k_{*j})/\sqrt{3}$ due to resonant amplification, and, under the narrow-width approximation, it contains an universal factor that can be interpreted as a result of momentum conservation. We also extend these discussions to the case of non-Gaussian perturbations.

Based on the Chinese historical sunspots drawings, a data set consisting of the scanned images and all their digitized parameters from 1925 to 2015 have been constructed. In this paper, we briefly describe the developmental history of sunspots drawings in China. This paper describes the preliminary processing processes that strat from the initial data (inputing to the scanning equipment) to the parameters extraction, and finally summarizes the general features of this dataset. It is the first systematic project in Chinese solar-physics community that the historical observation of sunspots drawings were digitized. Our data set fills in an almost ninety years historical gap, which span 60 degrees from east to west and 50 degrees from north to south and have no continuous and detailed digital sunspot observation information. As a complementary to other sunspots observation in the world, our dataset provided abundant information to the long term solar cycles solar activity research.

Low density regions are less affected by the nonlinear structure formation and baryonic physics. They are ideal places for probing the nature of dark energy, a possible explanation for the cosmic acceleration. Unlike void lensing, which requires identifications of individual voids, we study the stacked lensing signals around the low-density-positions (LDP), defined as places that are devoid of foreground bright galaxies in projection. The method allows a direct comparison with numerical results by drawing correspondence between the bright galaxies with halos. It leads to lensing signals that are significant enough for differentiating several dark energy models. In this work, we use the CFHTLenS catalogue to define LDPs, as well as measuring their background lensing signals. We consider several different definitions of the foreground bright galaxies (redshift range \& magnitude cut). Regarding the cosmological model, we run six simulations: the first set of simulations have the same initial conditions, with $\rm{w_{de}=-1,-0.5,-0.8,-1.2}$; the second set of simulations include a slightly different $\Lambda$CDM model and a w(z) model from \cite{2017NatAs...1..627Z}. The lensing results indicate that the models with $\rm{w_{de}=-0.5,-0.8}$ are not favored, and the other four models all achieve comparable agreement with the data.

In this study we report detailed observations of magnetic environment at four footpoints of two warm coronal loops observed on 5 May 2016 in NOAA AR 12542 (Loop I) and 17 Dec 2015 in NOAA AR 12470 (Loop II). These loops were connecting a plage region with sunspot periphery (Loop I) and a sunspot umbra (Loop II). We used Solar Dynamics Observatory (SDO) and Goode Solar Telescope (GST) data to describe the phenomenon and understand its causes. The study indicates loop brightening episodes were associated with magnetic flux emergence and cancellation processes observed in SDO's Helioseismic and Magnetic Imager (HMI) and GST's Near InfraRed Imaging Spectrapolarimeter (NIRIS) data. The observed activity was driven by magnetic reconnection between small-scale emerging dipoles and large-scale pre-existing fields, suggesting that the reconnection occurred in the lower chromosphere at the edge of an extended plage region, where the loops were rooted. We suggest that plasma, evaporated during these reconnection events, gradually filled the loops and as it cooled the visible density front propagated from one footpoint of the loop to another at a rate of 90-110 km s$^{-1}$. This study also indicates that at least some of the bright loops seen in SDO Atmospheric Imaging Assembly images rooted in sunspot umbra may be heated due to magnetic activity taking place at the remote (non-sunspot) footpoint.

We study the gamma-ray emission patterns and light curves in dissipative pulsar magnetospheres. We produce the gamma-ray light curves by using the geometric method and the particle trajectory method. For the geometric method, assuming the gamma-ray emission originates in a finite-width layer along the last closed lines, we generate the gamma-ray light curves based on the uniform emissivity along the magnetic field lines in the comoving frame (CF). For the particle trajectory method, we consider the spatial distribution of conductivity $\sigma$ by assuming a very high conductivity within the light cylinder and a finite conductivity outside the light cylinder (LC). Assuming that all the $\gamma$-ray emission originates in the outer magnetosphere outside the LC, we generate the gamma-ray light curves by computing realistic particle trajectories and Lorentz factors, taking into account both the accelerating electric field and curvature radiation loss. Further, we compare the modeling light curves to the observed light curves at $>0.1\, \rm GeV$ energies for Vela pulsar. Our results show that the magnetosphere with the low $\sigma$ value is more preferred for the geometric method. However, the magnetosphere with a near force-free regime within the LC and a high $\sigma$ value outside the LC is more favored for the particle trajectory method. It is noted that the particle trajectory method uses the parallel electric fields that are self-consistent with the magnetic fields of the magnetosphere. We suggest that the results from the particle trajectory method are more supported on the physical ground.

We perform as the first time hydrodynamic simulations to study the properties of hot accretion flow (HAF) around a neutron star (NS). The energy carried by the HAF will eventually be radiated out at the surface of the NS. The emitted photons can propagate inside the HAF and cool the HAF via Comptonization. We find that the Compton cooling can affect the properties of HAF around a NS significantly. We define the Eddington accretion rate as $\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ being the Eddington luminosity and the speed of light, respectively. We define $\dot m$ as the mass accretion rate at the NS surface in unit of $\dot M_{\rm Edd}$. When $\dot m > 10^{-4}$, Compton cooling can effectively cool the HAF and suppress wind. Therefore, the mass accretion rate is almost a constant with radius. The density profile is $\rho \propto r^{-1.4}$. When $\dot m < 10^{-4}$, the Compton cooling effects become weaker, wind becomes stronger, accretion rate is proportional to $r^{0.3-0.5}$. Consequently, the density profile becomes flatter, $\rho \propto r^{-1 \sim -0.8}$. When $\dot m < 10^{-6}$, the Compton cooling effects can be neglected. We find that with a same accretion rate, the temperature of HAF around a NS is significantly lower than that of HAF around a black hole (BH). Also, the Compton $y-$parameter of HAF around a NS is significantly smaller than that of HAF around a BH. This result predicts that HAF around a NS will produce a softer spectrum compared to HAF around a BH, which is consistent with observations.

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

Fiber collision is a persistent problem faced by modern spectroscopic galaxy surveys. In this work, we propose a new method to correct for this undesired effect, focusing on the clustering from the fiber-collision scale up to $\lesssim 10 \rm Mpc\,h^{-1}$. We assume that the fiber-collided galaxies are in association with their nearest three angular neighbors. Compared with the conventional nearest-neighbor method, we have properly accounted for the foreground (background) galaxies that are associated with the foreground (background) cosmic webs relative to the nearest neighbor. We have tested the new method with mock catalogs of the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). The test demonstrates that our new method can recover the projected two-point correlation functions at an accuracy better than 1\% on small (below the fiber-collision scale) to intermediate (i.e., $10 \rm Mpc\,h^{-1}$) scales, where the fiber collision takes effect and the SDSS main sample can probe. The new method also gives a better recovery of the redshift-space correlation functions almost on all scales that we are interested in.

We study the evolution of star-forming galaxies with $10^{10} M_\odot<M_*<10^{11.6} M_\odot$ over the redshift range of 0.7<z<1.2 using the emission line galaxies (ELGs) in the extended Baryon Oscillation Spectroscopic Survey (eBOSS). By applying the incomplete conditional stellar mass function (ICSMF) model proposed in Guo et al., we simultaneously constrain the sample completeness, the stellar--halo mass relation (SHMR) and the quenched galaxy fraction. We obtain the intrinsic stellar mass functions for star-forming galaxies in the redshift bins of 0.7<z<0.8, 0.8<z<0.9, 0.9<z<1.0 and 1.0<z<1.2, as well as the stellar mass function for all galaxies in the redshift bin of 0.7<z<0.8. We find that the eBOSS ELG sample only selects about 1%-10% of the star-forming galaxy population at the different redshifts, with the lower redshift samples to be more complete. There is only weak evolution in the SHMR of the ELGs from z=1.2 to z=0.7, as well as the intrinsic galaxy stellar mass functions for lower-mass galaxies of $M_*<10^{11} M_\odot$. There is significant decrease of the stellar mass function for star-forming galaxies with redshift at the massive end. Our best-fitting models show that the central ELGs at these redshifts live in halos of mass $M\sim10^{12} M_\odot$ while the satellite ELGs occupy slightly more massive halos of $M\sim10^{12.6} M_\odot$. The average satellite fraction of the observed ELGs varies from 13% to 17%, with the galaxy bias increasing from 1.1 to 1.4 from z=0.7 to 1.2.

The productions of data with internationally agreed standards are the common pursuits for all astronomical observation instruments both ground and spacebased undoubtedly, Since standard data to facilitate the use of scientific research. The routine solar observations are available at Huairou Solar Observing Station (HSOS) since 1987, which should be regarded as one main solar observations in the world. The data storage medium and format at HSOS experienced several changes, so there exist some inconveniences for solar physicist. This paper shows that the observations data of HSOS are further processed both for storage medium and format toward international standards, in order to explore HSOS observations data for scientific research.

The regular solar observations are operated at Huairou Solar Observing Station (HSOS) since 1987, which make the construction of long-term magnetic field datasets available to understand solar magnetic field and cycles. There exist some inconveniences for solar physicist to use these data, because the data storage medium and format at HSOS experienced some changes. Additionally, the processes of magnetic field calibration are not easy to deal with for who are not familiar with these data. Here shows that the magnetic field of HSOS are further processed toward international standards, in order to explore HSOS observations data for scientific research.

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.

Mass accretion rate is a key parameter in accretion disk theory. It determines black hole accretion mode. In large scale cosmological simulations studying galaxy formation and evolution, Bondi radius can at most be marginally resolved. In those simulations, Bondi accretion formula is always used to estimate black hole accretion rate. Bondi solution is too simple, which cannot represent the real accretion process. We perform simulations of hot accretion flow at parsec scale irradiated by a low luminosity active galactic nucleus (AGN). We perform 77 simulations with varying density and temperature at outer boundary (10 parsec). Our purpose is to find a formula to calculate real black hole accretion rate based on the gas density and temperature at parsec scale. We define Eddington accretion rate to be $\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ been Eddington luminosity and speed of light respectively. We set black hole mass to be $10^8M_{\odot}$, $M_{\odot}$ is solar mass. We find that black hole accretion rate can be expressed as $\frac{\dot M_{\rm BH} }{\dot M_{\rm Edd}}=10^{-3.11} (\rho_0/10^{-22}{\rm g cm}^{-3})^{1.36} (T_0/10^7 {\rm K})^{-1.9}$, with $\rho_0$ and $T_0$ being density and temperature at parsec scale, respectively. We find the formula can accurately predict the luminosity of observed low-luminosity AGNs (with black hole mass $\sim 10^8M_{\odot}$). This formula can be used in the sub-grid models in large scale cosmological simulations with a black hole mass of $\sim 10^8M_{\odot}$.

Observations of many dim galactic nuclei in local universe give good estimation of gas density and temperature at Bondi radius. If we assume the black hole accretes at Bondi accretion rate and radiates at efficiency of low-luminosity hot accretion flow, the predicted nuclei luminosity can be significantly higher than observations. Therefore, the real black hole mass accretion rate in these sources may be significantly smaller than the Bondi value. Active galactic nuclei (AGN) feedback may be responsible for decreasing black hole accretion rate to values much smaller than Bondi rate. We perform two-dimensional simulations of low angular momentum accretion flow at parsec and sub-parsec scales around low-luminosity active galactic nuclei (LLAGNs). We take into account the radiation and wind feedbacks of the LLAGN. The cross section of particle-particle interaction can be several orders of magnitude larger than that of photon-particle interaction. Therefore, we find that for the LLAGNs, effects of radiation feedback in decreasing black hole accretion rate are small. However, wind feedback can effectively decrease the black hole mass accretion rate. Due to the decrease of accretion rate, the black hole luminosity can be decreased by a factor of $\sim 33-400$. The results may be useful to explain why many galactic nuclei in the local universe are so dim.

Observations show that winds can be driven from the innermost region (inside a 50 Schwarschild radius) of a thin disk. It is interesting to study the winds launched from the innermost region. A hot corona above the black hole (BH) thin disk is irradiated by the disk. We perform two-dimensional hydrodynamical simulations to study the winds driven by radiation force from the corona in the innermost regions. The hard X-ray spectrum from active galactic nuclei (AGNs) suggests that the corona temperature is about $10^9$ K, so that we mainly analyze the properties of winds (or outflows) from the $10^9$ K corona. The disk luminosity plays an important role in driving the outflows. The more luminous the disk, the stronger the outflows. Mass outflow rate ($\dot{M}_{\rm out}$) at a 90 Schwarschild radius depends on disk luminosity, which can be described as $\dot{M}_{\rm out}\propto 10^{3.3 \Gamma}$ ($\Gamma$ is the ratio of the disk luminosity to the Eddington luminosity). In the case of high luminosity (e.g. $\Gamma=0.75$), the supersonic outflows with maximum speed $1.0 \times 10^4$ Km s$^{-1}$ are launched at $\sim17^{o}$ --$30^{o}$ and $\sim50^{o}$ --$80^{o}$ away from the pole axis. The Bernoulli parameter keeps increasing with the outward propagation of outflows. The radiation force keeps accelerating the outflows when outflows move outward. Therefore, we can expect the outflows to escape from the BH gravity and go to the galactic scale. The interaction between outflows and interstellar medium may be an important AGN feedback process.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

Halo merger trees are constructed from ELUCID, a constrained $N$-body simulation in the Sloan Digital Sky Survey (SDSS) volume. These merger trees are used to populate dark matter halos with galaxies according to an empirical model of galaxy formation. Mock catalogs in the SDSS sky coverage are constructed, which can be used to study the spatial distribution of galaxies in the low-$z$ Universe. These mock catalogs are used to quantify the cosmic variance in the galaxy stellar mass function (GSMF) measured from the SDSS survey. The GSMF estimated from the SDSS magnitude-limited sample can be affected significantly by the presence of the under-dense region at $z<0.03$, so that the low-mass end of the function can be underestimated significantly. Several existing methods designed to deal with the effects of the cosmic variance in the estimate of GSMF are tested, and none is found to be able to fully account for the cosmic variance. We propose a method based on the conditional stellar mass functions in dark matter halos, which can provide an unbiased estimate of the global GSMF. The application of the method to the SDSS data shows that the GSMF has a significant upturn at $M_*< 10^{9.5} h^{-1}{\rm M}_\odot$, which has been missed in many earlier measurements of the local GSMF.

Using a sample of galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and a catalog of bulge-disk decompositions, we study how the size distribution of galaxies depends on the intrinsic properties of galaxies, such as concentration, morphology, specific star formation rate (sSFR), and bulge fraction, and on the large-scale environments in the context of central/satellite decomposition, halo environment, the cosmic web: \cluster, \filament, \sheet ~and \void, as well as galaxy number density. We find that there is a strong dependence of the luminosity- or mass-size relation on the galaxy concentration, morphology, sSFR, and bulge fraction. Compared with late-type (spiral) galaxies, there is a clear trend of smaller sizes and steeper slope for early-type (elliptical) galaxies. Similarly, galaxies with high bulge fraction have smaller sizes and steeper slope than those with low bulge fraction. Fitting formula of the average luminosity- and mass-size relations are provided for galaxies of these different intrinsic properties. Examining galaxies in terms of their large scale environments, we find that the mass-size relation has some weak dependence on the halo mass and central/satellite segregation for galaxies within mass range $9.0\le \log M_{\ast} \le 10.5$, where satellites or galaxies in more massive halos have slightly smaller sizes than their counterparts. While the cosmic web and local number density dependence of the mass-size relation is almost negligible.

As we demonstrated in Paper I, the quenched fractions of central and satellite galaxies as function of halo mass are extremely similar, as long as one controls for stellar mass. The same holds for the quenched fractions as a function of central velocity dispersion, which is tightly correlated with black hole mass, as long as one controls for both stellar and halo mass. Here we use mock galaxy catalogs constructed from the latest semi-analytic model, L-GALAXIES, and the state-of-the-art hydrodynamical simulation, EAGLE, to investigate whether these models can reproduce the trends seen in the data. We also check how the group finder used to identify centrals and satellites impacts our results. We find that L-GALAXIES fails to reproduce the trends. The predicted quenched fraction of central galaxies increases sharply with halo mass around $10^{12.5}h^{-1}M_{\odot}$ and with black hole mass around $\sim10^{6.5}M_{\odot}$, while the predicted quenched fraction of satellites increases with both halo and black hole masses gradually. In contrast, centrals and satellites in EAGLE follow almost the same trend as seen in the data. We discuss the implications of our results for how feedback processes regulate galaxy quenching.

In paper I of this series (Yang et al. 2017, ApJ), we proposed an analytical method of blind separation ({\bf ABS}) to extract the cosmic magnification signal in galaxy number distribution and reconstruct the weak lensing power spectrum. Here we report a new version of the ABS method, with significantly improved performance. This version is characterized by a shift parameter $\mathcal{S}$, with the special case of $\mathcal{S}=0$ corresponding to the original ABS method. We have tested this new version, compared to the previous one, and confirmed its supreme performance in all investigated situations. Therefore it supersedes the previous version. The proof of concept studies presented in this paper demonstrate that it may enable surveys such as LSST and SKA to reconstruct the lensing power spectrum at $z\simeq 1$ with $1\%$ accuracy. We will test with more realistic simulations to verify its applicability in real data.

We extend the real-space mapping method developed in Shi et at. (2016) so that it can be applied to flux-limited galaxy samples. We use an ensemble of mock catalogs to demonstrate the reliability of this extension, showing that it allows for an accurate recovery of the real-space correlation functions and galaxy biases. We also demonstrate that, using an iterative method applied to intermediate-scale clustering data, we can obtain an unbiased estimate of the growth rate of structure $f\sigma_8$, which is related to the clustering amplitude of matter, to an accuracy of $\sim 10\%$. Applying this method to the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we construct a real-space galaxy catalog spanning the redshift range $0.01 \leq z \leq 0.2$, which contains 584,473 galaxies in the north Galactic cap (NGC). Using this data, we infer $\fss$ at a median redshift $z=0.1$, which is consistent with the WMAP9 cosmology at the $1\sigma$ level. By combining this measurement with the real-space clustering of galaxies and with galaxy-galaxy weak lensing measurements for the same sets of galaxies, we are able to break the degeneracy between $f$, $\sigma_8$, and $b$. From the SDSS DR7 data alone, we obtain the following cosmological constraints at redshift $z=0.1$: $f=$$0.464^{+0.040}_{-0.040}$, $\sigma_8=0.769^{+0.121}_{-0.089}$, and $b=1.910^{+0.234}_{-0.268}$, $1.449^{+0.194}_{-0.196}$, $1.301^{+0.170}_{-0.177}$, and $1.196^{+0.159}_{-0.161}~$ for galaxies within different absolute magnitude bins $^{0.1}{\rm M}_r-5\log h=[-23,0, -22.0], [-22,0, -21.0], [-21.0, -20.0]$ and $[-20.0, -19.0]$, respectively.

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 intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

It is often suggested that a super massive black hole is embedded in a nuclear bulge of size of a few $10^2$ $\text{ parsec}$ . The nuclear stars gravity is not negligible near $\sim 10 \text{ parsec}$. In order to study the effect of nuclear stars gravity on quasar radiation feedback on the parsec scale, we have simulated the parsec scale flows irradiated by a quasar by taking into account the gravitational potential of both the black hole and the nuclear star cluster. We find that the effect of nuclear stars gravity on the parsec-scale flows is related to the fraction of X-ray photons in quasar radiation. For the models in which the fraction of X-ray photons is not small (e.g. the X-ray photons contribute to $20\%$ of the quasar radiation), the nuclear stars gravity is very helpful to collimate the outflows driven by UV photons, significantly weakens the outflow power at the outer boundary and significantly enhances the net accretion rate onto the black hole. For the models in which X-ray photons are significantly decreased (e.g. the X-ray photons contribute to $5\%$ of the quasar radiation), the nuclear stars gravity can just slightly change properties of outflow and slightly enhance the net accretion rate onto the black hole.

We investigate the quenching properties of central and satellite galaxies, utilizing the halo masses and central-satellite identifications from the SDSS galaxy group catalog of Yang et al. We find that the quenched fractions of centrals and satellites of similar stellar masses have similar dependence on host halo mass. The similarity of the two populations is also found in terms of specific star formation rate and 4000 \AA\ break. The quenched fractions of centrals and satellites of similar masses show similar dependencies on bulge-to-total light ratio, central velocity dispersion and halo-centric distance in halos of given halo masses. The prevalence of optical/radio-loud AGNs is found to be similar for centrals and satellites at given stellar masses. All these findings strongly suggest that centrals and satellites of similar masses experience similar quenching processes in their host halos. We discuss implications of our results for the understanding of galaxy quenching.

We propose a novel method to constrain the missing fraction of galaxies using galaxy clustering measurements in the galaxy conditional stellar mass function (CSMF) framework, which is applicable to surveys that suffer significantly from sample selection effects. The clustering measurements, which are not sensitive to the random sampling (missing fraction) of galaxies, are widely used to constrain the stellar-halo mass relation (SHMR). By incorporating a missing fraction (incompleteness) component into the CSMF model (ICSMF), we use the incomplete stellar mass function and galaxy clustering to simultaneously constrain the missing fractions and the SHMRs. Tests based on mock galaxy catalogs with a few typical missing fraction models show that this method can accurately recover the missing fraction and the galaxy SHMR, and hence provides us reliable measurements of the galaxy stellar mass functions. We then apply it to the Baryon Oscillation Spectroscopic Survey (BOSS) over the redshift range of 0.1<z<0.8 for galaxies of $M_*>10^{11}M_\odot$. We find the sample completeness for BOSS is over 80% at z<0.6, but decreases at higher redshifts to about 30%. After taking these completeness factors into account, we provide accurate measurements of the stellar mass functions for galaxies with $10^{11}M_\odot<M_*<10^{12}M_\odot$, as well as the SHMRs, over the redshift range 0.1<z<0.8 in this largest galaxy redshift survey.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

We examine the thermal energy contents of the intergalactic medium (IGM) over three orders of magnitude in both mass density and gas temperature using thermal Sunyaev-Zel'dovich effect (tSZE). The analysis is based on {\it Planck} tSZE map and the cosmic density field, reconstructed for the SDSS DR7 volume and sampled on a grid of cubic cells of $(1h^{-1}{\rm Mpc})^3$, together with a matched filter technique employed to maximize the signal-to-noise. Our results show that the pressure - density relation of the IGM is roughly a power law given by an adiabatic equation of state, with an indication of steepening at densities higher than about $10$ times the mean density of the universe. The implied average gas temperature is $\sim 10^4\,{\rm K}$ in regions of mean density, $\rho_{\rm m} \sim {\overline\rho}_{\rm m}$, increasing to about $10^5\,{\rm K}$ for $\rho_{\rm m} \sim 10\,{\overline\rho}_{\rm m}$, and to $>10^{6}\,{\rm K}$ for $\rho_{\rm m} \sim 100\,{\overline\rho}_{\rm m}$. At a given density, the thermal energy content of the IGM is also found to be higher in regions of stronger tidal fields, likely due to shock heating by the formation of large scale structure and/or feedback from galaxies and AGNs. A comparison of the results with hydrodynamic simulations suggests that the current data can already provide interesting constraints on galaxy formation.

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 study the effects of accretion environment (gas density, temperature and angular momentum) at large radii ($\sim 10$pc) on luminosity of hot accretion flows. The radiative feedback effects from the accretion flow on the accretion environment are also self-consistently taken into account. We find that the slowly rotating flows at large radii can significantly deviate from Bondi accretion when radiation heating and cooling are considered. We further find that when the temperature of environment gas is low (e.g. $T=2\times 10^7$K), the luminosity of hot accretion flows is high. When the temperature of gas is high (e.g. $T\geq4\times 10^7$K), the luminosity of hot accretion flow significantly deceases. The environment gas density can also significantly influence the luminosity of accretion flows. When density is higher than $\sim 4\times 10^{-22}\text{g} \text{cm}^{-3}$ and temperature is lower than $2\times 10^7$K, hot accretion flow with luminosity lower than $2\%L_{\text{Edd}}$ is not present. Therefore, the pc-scale environment density and temperature are two important parameters to determine the luminosity. The results are also useful for the sub-grid models adopted by the cosmological simulations.

We study slowly rotating accretion flow at parsec and sub-parsec scale irradiated by a low luminosity active galactic nuclei. We take into account the Compton heating, photoionization heating by the central X-rays. The bremsstrahlung cooling, recombination and line cooling are also included. We find that due to the Compton heating, wind can be thermally driven. The power of wind is in the range $(10^{-6}-10^{-3}) L_{\rm Edd}$, with $L_{\rm Edd}$ being the Eddington luminosity. The mass flux of wind is in the range $(0.01-1) \dot M_{\rm Edd}$ ($\dot M_{\rm Edd}= L_{\rm Edd}/0.1c^2$ is the Eddington accretion rate, $c$ is speed of light). We define the wind generation efficiency as $\epsilon = P_W/\dot {M}_{\rm BH}c^2$, with $P_W$ being wind power, $\dot M_{\rm BH}$ being the mass accretion rate onto the black hole. $\epsilon$ lies in the rage $10^{-4}-1.18$. Wind production efficiency decreases with increasing mass accretion rate. The possible role of the thermally driven wind in the active galactic feedback is briefly discussed.

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general

Using a series of cosmological simulations that includes one dark-matter-only (DM-only) run, one gas cooling-star formation-supernovae feedback (CSF) run and one that additionally includes feedback from active galactic nuclei (AGNs), we classify the large-scale structures with both a velocity-shear-tensor code (Vweb) and a tidal-tensor code (Pweb). We find that the baryonic processes have almost no impact on large-scale structures -- at least not when classified using aforementioned techniques. More importantly, our results confirm that the gas component alone can be used to infer the filamentary structure of the Universe practically un-biased, which could be applied to cosmology constrains. In addition, the gas filaments are classified with its velocity (Vweb) and density (Pweb) fields, which can theoretically connect to the radio observations, such as HI surveys. This will help us to bias-freely link the radio observations with DM distributions at large scale.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

We apply a halo-based group finder to four large redshift surveys, the 2MRS, 6dFGS, SDSS and 2dFGRS, to construct group catalogs in the low-redshift Universe. The group finder is based on that of Yang et al. but with an improved halo mass assignment so that it can be applied uniformly to various redshift surveys of galaxies. Halo masses are assigned to groups according to proxies based on the stellar mass/luminosity of member galaxies. The performances of the group finder in grouping galaxies according to common halos and in halo mass assignments are tested using realistic mock samples constructed from hydrodynamical simulations and empirical models of galaxy occupation in dark matter halos. Our group finder finds $\sim 94\%$ of the correct true member galaxies for $90-95\%$ of the groups in the mock samples; the halo masses assigned by the group finder are un-biased with respect to the true halo masses, and have a typical uncertainty of $\sim0.2\,{\rm dex}$. The properties of group catalogs constructed from the observational samples are described and compared with other similar catalogs in the literature.

A matched filter technique is applied to the Planck all-sky Compton y-parameter map to measure the thermal Sunyaev-Zel'dovich (tSZ) effect produced by galaxy groups of different halo masses selected from large redshift surveys in the low-z Universe. Reliable halo mass estimates are available for all the groups, which allows us to bin groups of similar halo masses to investigate how the tSZ effect depends on halo mass over a large mass range. Filters are simultaneously matched for all groups to minimize projection effects. We find that the integrated y-parameter and the hot gas content it implies are consistent with the predictions of the universal pressure profile model only for massive groups above $10^{14}\,{\rm M}_\odot$, but much lower than the model prediction for low-mass groups. The halo mass dependence found is in good agreement with the predictions of a set of simulations that include strong AGN feedback, but simulations including only supernova feedback significantly over predict the hot gas contents in galaxy groups. Our results suggest that hot gas in galaxy groups is either effectively ejected or in phases much below the virial temperatures of the host halos.

In this paper we use high-resolution cosmological simulations to study halo intrinsic alignment and its dependence on mass, formation time and large-scale environment. In agreement with previous studies using N-body simulations, it is found that massive halos have stronger alignment. For given mass, older halos have stronger alignment than younger ones. By identifying the cosmic environment of halo using Hessian matrix, we find that for given mass, halos in cluster regions also have stronger alignment than those in filament. The existing theory has not addressed these dependencies explicitly. In this work we extend the linear alignment model with inclusion of halo bias and find that the halo alignment with its mass and formation time dependence can be explained by halo bias. However, the model can not account for the environment dependence, as it is found that halo bias is lower in cluster and higher in filament. Our results suggest that halo bias and environment are independent factors in determining halo alignment. We also study the halo alignment correlation function and find that halos are strongly clustered along their major axes and less clustered along the minor axes. The correlated halo alignment can extend to scale as large as $100h^{-1}$Mpc where its feature is mainly driven by the baryon acoustic oscillation effect.

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.

We propose an Analytical method of Blind Separation (ABS) of cosmic magnification from the intrinsic fluctuations of galaxy number density in the observed galaxy number density distribution. The ABS method utilizes the different dependences of the signal (cosmic magnification) and contamination (galaxy intrinsic clustering) on galaxy flux, to separate the two. It works directly on the measured cross galaxy angular power spectra between different flux bins. It determines/reconstructs the lensing power spectrum analytically, without assumptions of galaxy intrinsic clustering and cosmology. It is unbiased in the limit of infinite number of galaxies. In reality the lensing reconstruction accuracy depends on survey configurations, galaxy biases, and other complexities, due to finite number of galaxies and the resulting shot noise fluctuations in the cross galaxy power spectra. We estimate its performance (systematic and statistical errors) in various cases. We find that, stage IV dark energy surveys such as SKA and LSST are capable of reconstructing the lensing power spectrum at $z\simeq 1$ and $\ell\la 5000$ accurately. This lensing reconstruction only requires counting galaxies, and is therefore highly complementary to the cosmic shear measurement by the same surveys.

PKS 1155+251 is a radio-loud quasar source at z=0.203. Observations using very long baseline interferometry (VLBI) at ~2, 5, 8 and 15 GHz show that the structure of the radio source is quite complicated on parsec scales and that the outer hot spots are apparently undergoing a significant contraction. Because these results cannot be fully explained based on the compact symmetric object (CSO) scenario with a radio core located between the northern and southern complexes, we made observations with the Very Long Baseline Array (VLBA) at 24 and 43 GHz to search for compact substructures and alternative interpretations. The results show that the radio core revealed in the previous VLBI observations remains compact with a flat spectrum in our sub-milli-arcsecond--resolution images; the northern lobe emission becomes faint at 24 GHz and is mostly resolving out at 43 GHz; the southern complex is more bright but has been resolved into the brightest southern-end (S1) and jet or tail alike components westwards. Explaining the southern components aligned westward with a standard CSO scenario alone remains a challenge. As for the flatter spectral index of the southern-end component S1 between 24 and 43 GHz in our observations and the significant 15 GHz VLBA flux variability of S1, an alternative scenario is that the southern complex may be powered by a secondary black hole residing at S1. But more sensitive and high-resolution VLBI monitoring is required to discriminate the CSO and the binary black hole scenarios.

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.

We show that the ratio between the stellar mass of central galaxy and the mass of its host halo, $f_c \equiv M_{*,c}/M_{\rm h}$, can be used as an observable proxy of halo assembly time, in that galaxy groups with higher $f_c$ assembled their masses earlier. Using SDSS groups of Yang et al., we study how $f_c$ correlates with galaxy properties such as color, star formation rate, metallicity, bulge to disk ratio, and size. Central galaxies of a given stellar mass in groups with $f_c>0.02$ tend to be redder in color, more quenched in star formation, smaller in size, and more bulge dominated, as $f_c$ increases. The trends in color and star formation appear to reverse at $f_c<0.02$, reflecting a down-sizing effect that galaxies in massive halos formed their stars earlier although the host halos themselves assembled later (lower $f_c$). No such reversal is seen in the size of elliptical galaxies, suggesting that their assembly follows halo growth more closely than their star formation. Satellite galaxies of a given stellar mass in groups of a given halo mass tend to be redder in color, more quenched in star formation and smaller in size as $f_c$ increases. For a given stellar mass, satellites also tend to be smaller than centrals. The trends are stronger for lower mass groups. For groups more massive than $\sim 10^{13}{\rm M}_\odot$, a weak reversed trend is seen in color and star formation. The observed trends in star formation are qualitatively reproduced by an empirical model based on halo age abundance matching, but not by a semi-analytical model tested here.

The ELUCID project aims to build a series of realistic cosmological simulations that reproduce the spatial and mass distribution of the galaxies as observed in the Sloan Digital Sky Survey (SDSS). This requires powerful reconstruction techniques to create constrained initial conditions. We test the reconstruction method by applying it to several $N$-body simulations. We use 2 medium resolution simulations from each of which three additional constrained $N$-body simulations were produced. We compare the resulting friend of friend catalogs by using the particle indexes as tracers, and quantify the quality of the reconstruction by varying the main smoothing parameter. The cross identification method we use proves to be efficient, and the results suggest that the most massive reconstructed halos are effectively traced from the same Lagrangian regions in the initial conditions. Preliminary time dependence analysis indicates that high mass end halos converge only at a redshift close to the reconstruction redshift. This suggests that, for earlier snapshots, only collections of progenitors may be effectively cross-identified.

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.

The "unidentified" infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 $\mu$m are ubiquitously seen in various astrophysical regions. The UIE features are characteristic of the stretching and bending vibrations of aromatic hydrocarbons. The 3.3 $\mu$m feature resulting from aromatic C--H stretches is often accompanied by a weaker feature at 3.4 $\mu$m often attributed to aliphatic C--H stretches. The ratio of the observed intensity of the 3.3 $\mu$m aromatic C--H feature ($I_{3.3}$) to that of the 3.4 $\mu$m aliphatic C--H feature ($I_{3.4}$) allows one to estimate the aliphatic fraction (i.e. $N_{\rm C,aliph}/N_{\rm C,arom}$, the number of C atoms in aliphatic units to that in aromatic rings) of the UIE carriers, provided the intrinsic oscillator strengths of the 3.3 $\mu$m aromatic C--H stretch ($A_{3.3}$) and the 3.4 $\mu$m aliphatic C--H stretch ($A_{3.4}$) are known. In this article we summarize the computational results on $A_{3.3}$ and $A_{3.4}$ and their implications for the aromaticity and aliphaticity of the UIE carriers. We use density functional theory and second-order perturbation theory to derive $A_{3.3}$ and $A_{3.4}$ from the infrared vibrational spectra of seven PAHs with various aliphatic substituents (e.g., methyl-, dimethyl-, ethyl-, propyl-, butyl-PAHs, and PAHs with unsaturated alkyl-chains). The mean band strengths of the aromatic ($A_{3.3}$) and aliphatic ($A_{3.4}$) C--H stretches are derived and then employed to estimate the aliphatic fraction of the UIE carriers by comparing $A_{3.4}$/$A_{3.3}$ with $I_{3.4}$/$I_{3.3}$. We conclude that the UIE emitters are predominantly aromatic, as revealed by the observationally-derived ratio <$I_{3.4}$/$I_{3.3}$> ~ 0.12 and the computationally-derived ratio <$A_{3.4}$/$A_{3.3}$> ~ 1.76 which suggest an upper limit of $N_{\rm C,aliph}/N_{\rm C,arom}$ ~ 0.02 for the aliphatic fraction of the UIE carriers.

The so-called unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 $\mu$m ubiquitously seen in a wide variety of astrophysical regions are generally attributed to polycyclic aromatic hydrocarbon (PAH) molecules. Astronomical PAHs may have an aliphatic component as revealed by the detection in many UIE sources of the aliphatic C-H stretching feature at 3.4 $\mu$m. The ratio of the observed intensity of the 3.4 $\mu$m feature to that of the 3.3 $\mu$m aromatic C-H feature allows one to estimate the aliphatic fraction of the UIE carriers. This requires the knowledge of the intrinsic oscillator strengths of the 3.3 $\mu$m aromatic C-H stretch ($A_{3.3}$) and the 3.4 $\mu$m aliphatic C-H stretch ($A_{3.4}$). Lacking experimental data on $A_{3.3}$ and $A_{3.4}$ for the UIE candidate materials, one often has to rely on quantum-chemical computations. Although the second-order Moller-Plesset (MP2) perturbation theory with a large basis set is more accurate than the B3LYP density functional theory, MP2 is computationally very demanding and impractical for large molecules. Based on methylated PAHs, we show here that, by scaling the band strengths computed at an inexpensive level (e.g., B3LYP/6-31G*) we are able to obtain band strengths as accurate as that computed at far more expensive levels (e.g., MP2/6-311+G(3df,3pd)).

We analyze the photometric data and spectroscopic data that collect on the $\delta$ Scuti star AE UMa. The fundamental and the first overtone frequencies are confirmed as $f_{0} = 11.62560$\,c\,d$^{-1}$ and $f_{1} = 15.03124$\,c\,d$^{-1}$, respectively, from the frequency content by analyzing of the 40 nights light curve spanning from 2009 to 2012. Additionally, another 37 frequencies are identified as either the harmonics or the linear combinations of the fundamental and the first overtone frequencies, among which 25 are newly detected. The rate of period change of the fundamental mode is determined as $(1/P_{0})(dP_{0}/dt) = 5.4(\pm 1.9) \times 10^{-9}$\,yr$^{-1}$ as revealed from the $O-C$ diagram based on the 84 newly determined times of maximum light combined with those derived from the literature. The spectroscopic data suggests that AE UMa is a population\,I $\delta$ Scuti star. With these physical properties, we perform theoretical explorations based on the stellar evolution code MESA on this target, considering that the variation of pulsation period is caused by secular evolutionary effects. We finally constraint the AE UMa with the physical parameters as: the mass of $1.805 \pm 0.055\ M_{\odot}$, the radius of $1.647 \pm 0.032 \times 10^{11}$\,cm, the luminosity of $1.381 \pm 0.048~(\log L/L_{\odot})$ and the age of $1.055 \pm 0.095 \times 10^{9}$\,yr. AE UMa can be the (Pop. I) $\delta$ Scuti star that locates just after the second turn-off of its evolutional track leaving the main sequence, a star in the phase of the Hertzsprung Gap with a helium core and a hydrogen-burning shell.

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

The X-ray source CXO J133815.6+043255 has counterparts in the UV, optical, and radio bands. Based on the multi-band investigations, it has been recently proposed by Kim et al. (2015) as a rarely-seen off-nucleus ultraluminous X-ray (ULX) source with a black hole mass of >= 10^4 solar mass in the nearby Seyfert galaxy NGC 5252. To explore its radio properties at very high angular resolution, we performed very long-baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.7 GHz. We find that the radio counterpart is remarkably compact among the known ULXs. It does not show a resolved structure with a resolution of a few milliarcsecond (mas), and the total recovered flux density is comparable to that measured in earlier sub-arcsecond-resolution images. The compact radio structure, the relatively flat spectrum, and the high radio luminosity are consistent with a weakly accreting supermassive black hole in a low-luminosity active galactic nucleus. The nucleus of NGC 5252 itself has similar radio properties. We argue that the system represents a relatively rare pair of active galactic nuclei, where both components emit in the radio.

A method we developed recently for the reconstruction of the initial density field in the nearby Universe is applied to the Sloan Digital Sky Survey Data Release 7. A high-resolution N-body constrained simulation (CS) of the reconstructed initial condition, with $3072^3$ particles evolved in a 500 Mpc/h box, is carried out and analyzed in terms of the statistical properties of the final density field and its relation with the distribution of SDSS galaxies. We find that the statistical properties of the cosmic web and the halo populations are accurately reproduced in the CS. The galaxy density field is strongly correlated with the CS density field, with a bias that depend on both galaxy luminosity and color. Our further investigations show that the CS provides robust quantities describing the environments within which the observed galaxies and galaxy systems reside. Cosmic variance is greatly reduced in the CS so that the statistical uncertainties can be controlled effectively even for samples of small volumes.

The so-called unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer are ubiquitously seen in a wide variety of astrophysical regions. The UIE features are characteristic of the stretching and bending vibrations of aromatic hydrocarbon materials, e.g., polycyclic aromatic hydrocarbon (PAH) molecules. The 3.3 micrometer aromatic C--H stretching feature is often accompanied by a weaker feature at 3.4 micrometer. The latter is often thought to result from the C--H stretch of aliphatic groups attached to the aromatic systems. The ratio of the observed intensity of the 3.3 micrometer aromatic C--H feature to that of the 3.4 micrometer aliphatic C--H feature allows one to estimate the aliphatic fraction of the UIE carriers, provided that the intrinsic oscillator strengths of the 3.3 micrometer aromatic C--H stretch (A3.3) and the 3.4 micrometer aliphatic C--H stretch (A3.4) are known. While previous studies on the aliphatic fraction of the UIE carriers were mostly based on the A3.4/A3.3 ratios derived from the mono-methyl derivatives of small PAH molecules, in this work we employ density functional theory to compute the infrared vibrational spectra of several PAH molecules with a wide range of sidegroups including ethyl, propyl, butyl, and several unsaturated alkyl chains, as well as all the isomers of dimethyl-substituted pyrene. We find that, except PAHs with unsaturated alkyl chains, the corresponding A3.4/A3.3 ratios are close to that of mono-methyl PAHs. This confirms the predominantly-aromatic nature of the UIE carriers previously inferred from the A3.4/A3.3 ratio derived from mono-methyl PAHs.

Although it is generally accepted that the so-called "unidentified" infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer are characteristic of the stretching and bending vibrations of aromatic hydrocarbon materials, the exact nature of their carriers remains unknown: whether they are free-flying, predominantly aromatic gas-phase molecules, or amorphous solids with a mixed aromatic/aliphatic composition are being debated. Recently, the 3.3 and 3.4 micrometer features which are commonly respectively attributed to aromatic and aliphatic C-H stretches have been used to place an upper limit of ~2\% on the aliphatic fraction of the UIE carriers (i.e. the number of C atoms in aliphatic chains to that in aromatic rings). Here we further explore the aliphatic versus aromatic content of the UIE carriers by examining the ratio of the observed intensity of the 6.2 micrometer aromatic C-C feature (I6.2) to that of the 6.85 micrometer aliphatic C-H deformation feature (I6.85). To derive the intrinsic oscillator strengths of the 6.2 micrometer stretch (A6.2) and the 6.85 micrometer deformation (A6.85), we employ density functional theory to compute the vibrational spectra of seven methylated polycyclic aromatic hydrocarbon molecules and their cations. By comparing I6.85/I6.2 with A6.85/A6.2, we derive the fraction of C atoms in methyl(ene) aliphatic form to be at most ~10\%, confirming the earlier finding that the UIE emitters are predominantly aromatic. We have also computed the intrinsic strength of the 7.25 micrometer feature (A7.25), another aliphatic C-H deformation band. We find that A6.85 appreciably exceeds A7.25. This explains why the 6.85 micrometer feature is more frequently detected in space than the 7.25 micrometer feature.

We examine the consequences of the effective field theory (EFT) of dark matter-nucleon scattering for current and proposed direct detection experiments. Exclusion limits on EFT coupling constants computed using the optimum interval method are presented for SuperCDMS Soudan, CDMS II, and LUX, and the necessity of combining results from multiple experiments in order to determine dark matter parameters is discussed. We demonstrate that spectral differences between the standard dark matter model and a general EFT interaction can produce a bias when calculating exclusion limits and when developing signal models for likelihood and machine learning techniques. We also discuss the implications of the EFT for the next-generation (G2) direct detection experiments and point out regions of complementarity in the EFT parameter space.

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

As the first paper in a series on the study of the galaxy-galaxy lensing from Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our image processing pipeline that corrects the systematics primarily introduced by the Point Spread Function (PSF). Using this pipeline, we processed SDSS DR7 imaging data in $r$ band and generated a background galaxy catalog containing the shape information of each galaxy. Based on our own shape measurements of the galaxy images from SDSS DR7, we extract the galaxy-galaxy (GG) lensing signals around foreground spectroscopic galaxies binned in different luminosity and stellar mass. The overall signals are in good agreement with those obtained by \citet{Mandelbaum2005, Mandelbaum2006} from the SDSS DR4. The results in this paper with higher signal to noise ratio is due to the larger survey area than SDSS DR4, confirm that more luminous/massive galaxies bear stronger GG lensing signal. We also divide the foreground galaxies into red/blue and star forming/quenched subsamples and measured their GG lensing signals, respectively. We find that, at a specific stellar mass/luminosity, the red/quenched galaxies have relatively stronger GG lensing signals than their counterparts especially at large radii. These GG lensing signals can be used to probe the galaxy-halo mass relations and their environmental dependences in the halo occupation or conditional luminosity function framework.

All of the 14 subfields of the Kepler field have been observed at least once with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, Xinglong Observatory, China) during the 2012-2014 observation seasons. There are 88,628 reduced spectra with SNR$_g$ (signal-to-noise ratio in g band) $\geq$ 6 after the first round (2012-2014) of observations for the LAMOST-Kepler project (LK-project). By adopting the upgraded version of the LAMOST Stellar Parameter pipeline (LASP), we have determined the atmospheric parameters ($T_{\rm eff}$ , $\log g$, and $\rm [Fe/H]$) and heliocentric radial velocity $v_{\rm rad}$ for 51,406 stars with 61,226 spectra. Compared with atmospheric parameters derived from both high-resolution spectroscopy and asteroseismology method for common stars in Huber et al. (2014), an external calibration of LASP atmospheric parameters was made, leading to the determination of external errors for the giants and dwarfs, respectively. Multiple spectroscopic observations for the same objects of the LK-project were used to estimate the internal uncertainties of the atmospheric parameters as a function of SNR$_g$ with the unbiased estimation method. The LASP atmospheric parameters were calibrated based on both the external and internal uncertainties for the giants and dwarfs, respectively. A general statistical analysis of the stellar parameters leads to discovery of 106 candidate metal-poor stars, 9 candidate very metal-poor stars, and 18 candidate high-velocity stars. Fitting formulae were obtained segmentally for both the calibrated atmospheric parameters of the LK-project and the KIC parameters with the common stars. The calibrated atmospheric parameters and radial velocities of the LK-project will be useful for studying stars in the Kepler field.

Active galactic nuclei (AGN) is one of the main drivers for transition from star-forming disk to passive spheroidal galaxies. However, the role of large-scale environment versus one-on-one interactions in triggering different types of AGN is still uncertain. We present a statistical study of the prevalence of the nuclear activity in isolated galaxies and physically bound isolated pairs. For the purpose of this study we considered optically and radio selected nuclear activity types. We aim to assess the effect of one-on-one interaction on the fraction of AGN and the role of their large-scale environment. To study the effect of one-on-one interaction on the fraction of AGN in isolated galaxy pairs, we compare with a sample of isolated galaxies homogeneously selected under the same isolation criterion. We examine the effect of the large-scale environment by comparing with control samples of single galaxies and galaxy pairs. In general we found no difference in the prevalence of optical AGN for the considered samples. For massive galaxies, the fraction of optical AGN in isolated galaxies is slightly higher than that in control samples. Also the fraction of passives in high mass isolated galaxies is smaller than in any other sample. Generally, there is no dependence on optical nuclear activity with local environment. On the other hand, we found evidence that radio AGN are strongly affected by the local environment. Optical AGN phenomenon is related to cold gas accretion, while radio AGN is related to hot gas accretion. In this context, there is more cold gas, fueling the central optical AGN, in isolated systems. Our results are in agreement with a scenario where cold gas accretion by secular evolution is the main driver of optical AGN, while hot gas accretion and one-on-one interactions are the main drivers of radio AGN activity.

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

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.

A procedure is introduced to recognise sunspots automatically in solar full-disk photosphere images obtained from Huairou Solar Observing Station, National Astronomical Observatories of China. The images are first pre-processed through Gaussian algorithm. Sunspots are then recognised by the morphological Bot-hat operation and Otsu threshold. Wrong selection of sunspots is eliminated by a criterion of sunspot properties. Besides, in order to calculate the sunspots areas and the solar centre, the solar limb is extracted by a procedure using morphological closing and erosion operations and setting an adaptive threshold. Results of sunspot recognition reveal that the number of the sunspots detected by our procedure has a quite good agreement with the manual method. The sunspot recognition rate is 95% and error rate is 1.2%. The sunspot areas calculated by our method have high correlation (95%) with the area data from USAF/NOAA.

Numerical simulations of black hole hot accretion flows have shown the existence of strong wind. Those works focus only on the region close to black hole thus it is unknown whether or where the wind production stops at large radii. To address this question, Bu et al. (2016) have performed hydrodynamic (HD) simulations by taking into account the gravitational potential of both the black hole and the nuclear star clusters. The latter is assumed to be $\propto \sigma^2 \ln(r)$, with $\sigma$ being the velocity dispersion of stars and $r$ be the distance from the center of the galaxy. It was found that when the gravity is dominated by nuclear stars, i.e., outside of radius $R_A\equiv GM_{\rm BH}/\sigma^2$, winds can no longer be produced. That work, however, neglects the magnetic field, which is believed to play a crucial dynamical role in the accretion and thus must be taken into account. In this paper, we revisit this problem by performing magneto-hydrodynamical (MHD) simulations. We confirm the result of Bu et al. (2016), namely wind can't be produced at the region of $R>R_A$. Our result, combined with the results of Yuan et al. (2015), indicates that the formula describing the mass flux of wind $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$ can only be applied to the region where the black hole potential is dominant. Here $\dot{M}_{\rm BH}$ is the mass accretion rate at the black hole horizon and the value of $R_A$ is similar to the Bondi radius.

Previous works show strong winds exist in hot accretion flows around black holes. Those works focus only on the region close to the black hole thus it is unknown whether or where the wind production stops at large radii. In this paper, we investigate this problem by hydrodynamical simulations. We take into account gravities of both the black hole and the nuclear star clusters. For the latter, we assume that the velocity dispersion of stars is a constant and its gravitational potential $\propto \sigma^2 \ln (r)$, where $\sigma$ is the velocity dispersion of stars and $r$ is the distance from the center of the galaxy. We focus on the region where the gravitational potential is dominated by the star cluster. We find, same as the accretion flow at small radii, the mass inflow rate decreases inward and the flow is convectively unstable. However, trajectory analysis shows that there is very few wind launched from the flow. Our result, combined with the results of Yuan et al. (2015), indicates that the mass flux of wind launched from hot accretion flow $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$, with $r\la R_A\equiv GM_{\rm BH}/\sigma^2$. Here $\dot{M}_{\rm BH}$ is accretion rate at black hole horizon. $R_A$ is similar to Bondi radius. We argue that the inward decrease of inflow rate is not because of mass loss via wind, but because of convective motion. The disappearance of wind outside of $R_A$ must be because of the change of the gravitational potential, but the exact reason remains to be probed.

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.

CDMS II data from the 5-tower runs at the Soudan Underground Laboratory were reprocessed with an improved charge-pulse fitting algorithm. Two new analysis techniques to reject surface-event backgrounds were applied to the 612 kg days germanium-detector WIMP-search exposure. An extended analysis was also completed by decreasing the 10 keV analysis threshold to $\sim$5 keV, to increase sensitivity near a WIMP mass of 8 GeV/$c^2$. After unblinding, there were zero candidate events above a deposited energy of 10 keV and 6 events in the lower-threshold analysis. This yielded minimum WIMP-nucleon spin-independent scattering cross-section limits of $1.8 \times 10^{-44}$ and $1.18 \times 10 ^{-41}$ cm$^2$ at 90\% confidence for 60 and 8.6 GeV/$c^2$ WIMPs, respectively. This improves the previous CDMS II result by a factor of 2.4 (2.7) for 60 (8.6) GeV/$c^2$ WIMPs.

The nearly continuous light curves with micromagnitude precision provided by the space mission Kepler are revolutionising our view of pulsating stars. They have revealed a vast sea of low-amplitude pulsation modes that were undetectable from Earth. The long time base of Kepler light curves allows an accurate determination of frequencies and amplitudes of pulsation modes needed for in-depth asteroseismic modeling. However, for an asteroseismic study to be successful, the first estimates of stellar parameters need to be known and they can not be derived from the Kepler photometry itself. The Kepler Input Catalog (KIC) provides values for the effective temperature, the surface gravity and the metallicity, but not always with a sufficient accuracy. Moreover, information on the chemical composition and rotation rate is lacking. We are collecting low-resolution spectra for objects in the Kepler field of view with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, Xinglong observatory, China). All of the requested fields have now been observed at least once. In this paper we describe those observations and provide a database of use to the whole astronomical community.

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.

We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.

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 empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

By combining the data of the Two Micron All Sky Survey, the \textit{Wide Field Infrared Survey Explorer} and the \textit{AKARI} satellite, we study the infrared colour properties of a sample of 2712 nearby radio-luminous galaxies (RLGs). These RLGs are divided into radio-loud (RL) active galactic nuclei (AGNs), mainly occurring at redshifts of $0.05<z<0.3$ and star-forming-dominated RLGs (SFGs), mainly occurring at redshifts of $0.01<z<0.15$. RL AGNs and SFGs are separately distributed in the ([3.4]$-$[4.6])$-$([4.6]$-$[12]) two-colour diagram, in which the RL AGNs display a double-core distribution, and the SFGs display a single-core distribution. SFGs have a redder [4.6]$-$[12] colour than RL AGNs due to the significant contribution from the dust component of SFGs. We find simple criteria of MIR colour separation between RL AGNs and SFGs such that: 95$\%$ of RL AGNs have [4.6]$-$[12] $<$ 3.0 and 94$\%$ of SFGs have [4.6]$-$[12] $>$ 3.0. We also analyse the MIR colours of RL AGNs divided into low- and high-excitation radio galaxies (LERGs and HERGs, respectively). The ([3.4]$-$[4.6])$-$([4.6]$-$[12]) diagram clearly shows separate distributions of LERGs and HERGs and a region of overlap, which suggests that LERGs and HERGs have different MIR properties. LERGs are responsible for the double-core distribution of RL AGNs on the ([3.4]$-$[4.6])$-$([4.6]$-$[12]) diagram. In addition, we also suggest 90$-$140$\mu$m band spectral index $\alpha(90,140)<-1.4$ as a criterion of selecting nearby active galaxies with non-thermal emissions at far-infrared wavelengths.

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.

Reliable halo mass estimation for a given galaxy system plays an important role both in cosmology and galaxy formation studies. Here we set out to find the way that can improve the halo mass estimation for those galaxy systems with limited brightest member galaxies been observed. Using four mock galaxy samples constructed from semi-analytical formation models, the subhalo abundance matching method and the conditional luminosity functions, respectively, we find that the luminosity gap between the brightest and the subsequent brightest member galaxies in a halo (group) can be used to significantly reduce the scatter in the halo mass estimation based on the luminosity of the brightest galaxy alone. Tests show that these corrections can significantly reduce the scatter in the halo mass estimations by $\sim 50\%$ to $\sim 70\%$ in massive halos depending on which member galaxies are considered. Comparing to the traditional ranking method, we find that this method works better for groups with less than five members, or in observations with very bright magnitude cut.

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 derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the value of red centrals. 2) The intrinsic scatters of the SHMRs of red and blue centrals are ~0.14 and ~0.11dex, respectively. The intrinsic scatter of the average SHMR of all central galaxies changes from ~0.20dex to ~0.14dex in the 11.3<logMh<15 range. 3) The fraction of halos hosting blue centrals at Mh=1E11Msun is 87%, but at 2x1E12Msun decays to ~20%, approaching to a few per cents at higher masses. The characteristic mass at which this fraction is the same for blue and red galaxies is Mh~7x1E11Msun. Our results suggest that the SHMR of central galaxies at large masses is shaped by halo mass quenching (likely through shock virial heating and AGN feedback), but group richness also plays an important role: central galaxies living in less dense environments quenched their star formation later or did not quench it yet. At low masses, processes that delay star formation without invoking too strong supernova-driven outflows could explain the high Ms/Mh ratios of blue centrals as compared to those of the scarce red centrals.

Using a sample of spiral galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and Galaxy Zoo 2 (GZ2), we investigate the alignment of spin axes of spiral galaxies with their surrounding large scale structure, which is characterized by the large-scale tidal field reconstructed from the data using galaxy groups above a certain mass threshold. We find that the spin axes of only have weak tendency to be aligned with (or perpendicular to) the intermediate (or minor) axis of the local tidal tensor. The signal is the strongest in a \cluster environment where all the three eigenvalues of the local tidal tensor are positive. Compared to the alignments between halo spins and local tidal field obtained in N-body simulations, the above observational results are in best agreement with those for the spins of inner regions of halos, suggesting that the disk material traces the angular momentum of dark matter halos in the inner regions.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) at the Xinglong observatory in China is a new 4-m telescope equipped with 4,000 optical fibers. In 2010, we initiated the LAMOST-Kepler project. We requested to observe the full field-of-view of the nominal Kepler mission with the LAMOST to collect low-resolution spectra for as many objects from the KIC10 catalogue as possible. So far, 12 of the 14 requested LAMOST fields have been observed resulting in more than 68,000 low-resolution spectra. Our preliminary results show that the stellar parameters derived from the LAMOST spectra are in good agreement with those found in the literature based on high-resolution spectroscopy. The LAMOST data allows to distinguish dwarfs from giants and can provide the projected rotational velocity for very fast rotators.

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input non-linear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle Mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. The constrained simulations so obtained accurately reproduce both the amplitudes and phases of the input simulations at various $z$. Using a PM model with a grid cell size of 0.75 Mpc/h and 40 time-steps in the HMC can recover more than half of the phase information down to a scale k~0.85 h/Mpc at high z and to k~3.4 h/Mpc at z=0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.

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.

Major mergers and interactions between gas-rich galaxies with comparable masses are thought to be the main triggers of starburst. In this work, we study, for a large stellar mass range, the interaction rate of the starburst galaxies in the local universe. We focus independently on central and satellite star forming galaxies extracted from the Sloan Digital Sky Survey. Here the starburst galaxies are selected in the star formation rate (SFR) stellar mass plane with SFR five times larger than the median value found for "star forming" galaxies of the same stellar mass. Through visual inspection of their images together with close companions determined using spectroscopic redshifts, we find that ~50% of the "starburst" populations show evident merger features, i.e., tidal tails, bridges between galaxies, double cores and close companions. In contrast, in the control sample we selected from the normal star forming galaxies, only ~19% of galaxies are associated with evident mergers. The interaction rates may increase by ~5% for the starburst sample and 2% for the control sample if close companions determined using photometric redshifts are considered. The contrast of the merger rate between the two samples strengthens the hypothesis that mergers and interactions are indeed the main causes of starburst.

We use 12000 stars from Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic data to show that the metallicities of Kepler field stars as given in the Kepler Input Catalog (KIC) systematically underestimate both the true metallicity and the dynamic range of the Kepler sample. Specifically, to the first order approximation, we find [Fe/H]_KIC = -0.20 + 0.43 [Fe/H]_LAMOST, with a scatter of ~0.25 dex, due almost entirely to errors in KIC. This relation is most secure for -0.3<[Fe/H]_LAMOST<+0.4 where we have >200 comparison stars per 0.1 dex bin and good consistency is shown between metallicities determined by LAMOST and high-resolution spectra. It remains approximately valid in a slightly broader range. When the relation is inverted, the error in true metallicity as derived from KIC is (0.25 dex)/0.43~0.6 dex. We thereby quantitatively confirm the cautionary note by Brown et al. (2011) that KIC estimates of [Fe/H] should not be used by "anyone with a particular interest in stellar metallicities". Fortunately, many more LAMOST spectroscopic metallicities will be available in the near future.

We present high spatial resolution X-ray spectroscopy of supernova remnant G292.0+1.8 with the {\sl Chandra} observations. The X-ray emitting region of this remnant was divided into 25 $\times$ 25 pixels with a scale of 20$\arcsec$ $\times$ 20$\arcsec$ each. Spectra of 324 pixels were created and fitted with an absorbed one component non-equilibrium ionization model. With the spectral analysis results we obtained maps of absorbing column density, temperature, ionization age, and the abundances for O, Ne, Mg, Si, S, and Fe. The abundances of O, Ne and Mg show tight correlations between each other in the range of about two orders of magnitude, suggesting them all from explosive C/Ne burning. Meanwhile, the abundances of Si and S are also well correlated, indicating them to be the ashes of explosive explosive O-burning or incomplete Si-burbing. The Fe emission lines are not prominent among the whole remnant, and its abundance are significantly deduced, indicating that the reverse shock may have not propagated to the Fe-rich ejecta. Based on relative abundances of O, Ne, Mg, Si and Fe to Si, we suggest a progenitor mass of $25-30 M_{\odot}$ for this remnant.

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{\theta\phi}$ component of the viscous stress and consider various values of the viscous parameter $\alpha$. We find that when $T_{\theta\phi}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $\alpha$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $\rho(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $\alpha$ ranging from $\alpha=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $\alpha$, the flow is marginally stable (when $\alpha$ is small) or unstable (when $\alpha$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

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.

Context : We still do not know which mechanisms are responsible for the transport of angular momentum inside stars. The recent detection of mixed modes that contain the signature of rotation in the spectra of Kepler subgiants and red giants gives us the opportunity to make progress on this issue. Aims: Our aim is to probe the radial dependance of the rotation profiles for a sample of Kepler targets. For this purpose, subgiants and early red giants are particularly interesting targets because their rotational splittings are more sensitive to the rotation outside the deeper core than is the case for their more evolved counterparts. Methods: We first extract the rotational splittings and frequencies of the modes for six young Kepler red giants. We then perform a seismic modeling of these stars using the evolutionary codes CESAM2k and ASTEC. By using the observed splittings and the rotational kernels of the optimal models, we perform inversions of the internal rotation profiles of the six stars. Results: We obtain estimates of the mean rotation rate in the core and in the convective envelope of these stars. We show that the rotation contrast between the core and the envelope increases during the subgiant branch. Our results also suggest that the core of subgiants spins up with time, contrary to the RGB stars whose core has been shown to spin down. For two of the stars, we show that a discontinuous rotation profile with a deep discontinuity reproduces the observed splittings significantly better than a smooth rotation profile. Interestingly, the depths that are found most probable for the discontinuities roughly coincide with the location of the H-burning shell, which separates the layers that contract from those that expand. These results will bring observational constraints to the scenarios of angular momentum transport in stars.

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions.

We examine the statistical properties of the brightest group galaxies (BGGs) using a complete spectroscopic sample of groups/clusters of galaxies selected from the Data Release 7 of the Sloan Digital Sky Survey. We test whether BGGs and other bright members of groups are consistent with an ordered population among the total population of group galaxies. We find that the luminosity distributions of BGGs do not follow the predictions from the order statistics (OS). The average luminosities of BGGs are systematically brighter than OS predictions. On the other hand, by properly taking into account the brightening effect of the BGGs, the luminosity distributions of the second brightest galaxies are in excellent agreement with the expectations of OS. The brightening of BGGs relative to the OS expectation is consistent with a scenario that the BGGs on average have over-grown about 20 percent masses relative to the other member galaxies. The growth ($\Delta M$) is not stochastic but correlated with the magnitude gap ($G_{1,2}$) between the brightest and the second brightest galaxy. The growth ($\Delta M$) is larger for the groups having more prominent BGGs (larger $G_{1,2}$) and averagely contributes about 30 percent of the final $G_{1,2}$ of the groups of galaxies.

We use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}.

We select satellite galaxies from the galaxy group catalog constructed with the SDSS spectroscopic galaxies and measure the tangential shear around these galaxies with source catalog extracted from CFHT/MegaCam Stripe-82 Survey to constrain the mass of subhalos associated with them. The lensing signal is measured around satellites in groups with masses in the range [10^{13}, 5x10^{14}]h^{-1}M_{sun}, and is found to agree well with theoretical expectation. Fitting the data with a truncated NFW profile, we obtain an average subhalo mass of log M_{sub}= 11.68 \pm 0.67 for satellites whose projected distances to central galaxies are in the range [0.1, 0.3] h^{-1}Mpc, and log M_{sub}= 11.68 \pm 0.76 for satellites with projected halo-centric distance in [0.3, 0.5] h^{-1}Mpc. The best-fit subhalo masses are comparable to the truncated subhalo masses assigned to satellite galaxies using abundance matching and about 5 to 10 times higher than the average stellar mass of the lensing satellite galaxies.

In this paper we present a new public code, named $ynogkm$, for the fast calculation of time-like geodesics in the Kerr-Newmann (K-N) spacetime, which is a direct extension of $ynogk$ calculating null geodesics in a Kerr spacetime. Following the strategies used in $ynogk$, we also solve the equations of motion analytically and semi-analytically by using Weierstrass' and Jacobi's elliptic functions and integrals, in which the Boyer-Lidquist (B-L) coordinates $r$, $\theta$, $\phi$, $t$ and the proper time $\sigma$ are expressed as functions of an independent variable $p$ (Mino time). All of the elliptic integrals are computed by Carlson's elliptic integral method, which guarantees the fast speed of the code. Finally the code is applied to a couple of toy problems.

Using a sample of galaxy groups selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we examine the alignment between the orientation of galaxies and their surrounding large scale structure in the context of the cosmic web. The latter is quantified using the large-scale tidal field, reconstructed from the data using galaxy groups above a certain mass threshold. We find that the major axes of galaxies in filaments tend to be preferentially aligned with the directions of the filaments, while galaxies in sheets have their major axes preferentially aligned parallel to the plane of the sheets. The strength of this alignment signal is strongest for red, central galaxies, and in good agreement with that of dark matter halos in N-body simulations. This suggests that red, central galaxies are well aligned with their host halos, in quantitative agreement with previous studies based on the spatial distribution of satellite galaxies. There is a luminosity and mass dependence that brighter and more massive galaxies in filaments and sheets have stronger alignment signals. We also find that the orientation of galaxies is aligned with the eigenvector associated with the smallest eigenvalue of the tidal tensor. These observational results indicate that galaxy formation is affected by large-scale environments, and strongly suggests that galaxies are aligned with each other over scales comparable to those of sheets and filaments in the cosmic web.

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.

Many researchers have performed cosmological-model-independent tests for the distance duality (DD) relation. Theoretical work has been conducted based on the results of these tests. However, we find that almost all of these tests were perhaps not cosmological-model-independent after all, because the distance moduli taken from a given type Ia supernovae (SNe Ia) compilation are dependent on a given cosmological model and Hubble constant. In this Letter, we overcome these defects and by creating a new cosmological-model-independent test for the DD relation. We use the original data from the Union2 SNe Ia compilation and the angular diameter distances from two galaxy cluster samples compiled by De Filippis et al. and Bonamente et al. to test the DD relation. Our results suggest that the DD relation is compatible with observations, and the spherical model is slightly better than the elliptical model at describing the intrinsic shape of galaxy clusters if the DD relation is valid. However, these results are different from those of previous work.

A cosmological preferred direction was reported from the type Ia supernovae (SNe Ia) data in recent years. We use the Union2.1 data to give a simple classification of such studies for the first time. Because the maximum anisotropic direction is independent of isotropic dark energy models, we adopt two cosmological models ($\Lambda$CDM, $w$CDM) for the hemisphere comparison analysis and $\Lambda$CDM model for dipole fit approach. In hemisphere comparison method, the matter density and the equation of state of dark energy are adopted as the diagnostic qualities in the $\Lambda$CDM model and $w$CDM model, respectively. In dipole fit approach, we fit the fluctuation of distance modulus. We find that there is a null signal for the hemisphere comparison method, while a preferred direction ($b=-14.3^\circ \pm 10.1^\circ, l=307.1^\circ \pm 16.2^\circ$) for the dipole fit method. This result indicates that the dipole fit is more sensitive than the hemisphere comparison method.

The existence of galaxy intrinsic clustering severely hampers the weak lensing reconstruction from cosmic magnification. In paper I \citep{Yang2011}, we proposed a minimal variance estimator to overcome this problem. By utilizing the different dependences of cosmic magnification and galaxy intrinsic clustering on galaxy flux, we demonstrated that the otherwise overwhelming galaxy intrinsic clustering can be significantly suppressed such that lensing maps can be reconstructed with promising accuracy. This procedure relies heavily on the accuracy of determining the galaxy bias from the same data. Paper I adopts an iterative approach, which degrades toward high redshift. The current paper presents an alternative method, improving over paper I. We prove that the measured galaxy clustering between flux bins allows for simultaneous determination of the lensing power spectrum and the flux dependence of galaxy bias, at this redshift bin. Comparing to paper I, the new approach is not only more straightforward, but also more robust. It identifies an ambiguity in determining the galaxy bias and further discovers a mathematically robust way to suppress this ambiguity to non-negligible level ($\sim 0.1%$). The accurately determined galaxy bias can then be applied to the minimal variance estimator proposed in paper I to improve the lensing map-making. The gain at high redshift is significant. These maps can be used to measure other statistics, such as cluster finding and peak statistics. Furthermore, by including galaxy clustering measurement between different redshift bins, we can also determine the lensing cross power spectrum between these bins, up to a small and correctable multiplicative factor.

The unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer, commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules, have been recently ascribed to mixed aromatic/aliphatic organic nanoparticles. More recently, an upper limit of <9% on the aliphatic fraction (i.e., the fraction of carbon atoms in aliphatic form) of the UIE carriers based on the observed intensities of the 3.4 and 3.3 micrometer emission features by attributing them to aliphatic and aromatic C-H stretching modes, respectively, and assuming A_34./A_3.3~0.68 derived from a small set of aliphatic and aromatic compounds, where A_3.4 and A_3.3 are respectively the band strengths of the 3.4 micrometer aliphatic and 3.3 micrometer aromatic C-H bonds. To improve the estimate of the aliphatic fraction of the UIE carriers, here we analyze 35 UIE sources which exhibit both the 3.3 and 3.4 micrometer C-H features and determine I_3.4/I_3.3, the ratio of the power emitted from the 3.4 micrometer feature to that from the 3.3 micrometer feature. We derive the median ratio to be <I_3.4/I_3.3> ~ 0.12. We employ density functional theory and second-order perturbation theory to compute A_3.4/A_3.3 for a range of methyl-substituted PAHs. The resulting A_3.4/A_3.3 ratio well exceeds 1.4, with an average ratio of <A_3.4/A_3.3> ~1.76. By attributing the 3.4 micrometer feature exclusively to aliphatic C-H stretch (i.e., neglecting anharmonicity and superhydrogenation), we derive the fraction of C atoms in aliphatic form to be ~2%. We therefore conclude that the UIE emitters are predominantly aromatic.

Based on several magnetic nonpotentiality parameters obtained from the vector photospheric active region magnetograms obtained with the Solar Magnetic Field Telescope at the Huairou Solar Observing Station over two solar cycles, a machine learning model has been constructed to predict the occurrence of flares in the corresponding active region within a certain time window. The Support Vector Classifier, a widely used general classifier, is applied to build and test the prediction models. Several classical verification measures are adopted to assess the quality of the predictions. We investigate different flare levels within various time windows, and thus it is possible to estimate the rough classes and erupting times of flares for particular active regions. Several combinations of predictors have been tested in the experiments. The True Skill Statistics are higher than 0.36 in 97% of cases and the Heidke Skill Scores range from 0.23 to 0.48. The predictors derived from longitudinal magnetic fields do perform well, however they are less sensitive in predicting large flares. Employing the nonpotentiality predictors from vector fields improves the performance of predicting large flares of magnitude $\geq$M5.0 and $\geq$X1.0.

Assembly bias describes the finding that the clustering of dark matter haloes depends on halo formation time at fixed halo mass. In this paper, we analyse the influence of assembly bias on galaxy clustering using both semi-analytical models (SAMs) and observational data. At fixed stellar mass, SAMs predict that the clustering of {\it central} galaxies depends on the specific star formation rate (sSFR), with more passive galaxies having a higher clustering amplitude. We find similar trends using SDSS group catalogues, and verify that these are not affected by possible biases due to the group finding algorithm. Low mass central galaxies reside in narrow bins of halo mass, so the observed trends of higher clustering amplitude for galaxies with lower sSFR is not driven by variations of the parent halo mass. We argue that the clustering dependence on sSFR represent a direct detection of assembly bias. In addition, contrary to what expected based on clustering of dark matter haloes, we find that low-mass central galaxies in SAMs with larger host halo mass have a {\it lower} clustering amplitude than their counter-parts residing in lower mass haloes. This results from the fact that, at fixed stellar mass, assembly bias has a stronger influence on clustering than the dependence on the parent halo mass.

We perform time-dependent, 2DHD numerical simulations to study the dynamics of a slowly rotating accretion flow from sub-pc to pc scales under the irradiation from the central AGN. Compared to previous work, we improve the calculation of the radiative force due to X-rays. More importantly, in addition to radiative pressure and radiative heating/cooling directly from the central AGN, in the momentum equation we also include the force due to the scattered and reprocessed photons. We find that the accretion flow properties change significantly due to this "re-radiation" effect. The inflow rate at the inner boundary is reduced, while the outflow rate at the outer boundary is enhanced by about one order of magnitude. This effect is more significant when the density at the outer boundary is higher. The properties of outflows such as velocity, momentum and energy fluxes, and the ratio of outflow rate and the accretion rate, are calculated. We find that the efficiency of transferring the radiation power into the kinetic power of outflow is typically $10^{-3}$, far below the value of $\sim 0.05$ which is assumed in some cosmological simulations. The effect of the temperature of the gas at the outer boundary ($T_0$) is investigated. When $T_0$ is high, the emitted luminosity of the accretion flow oscillates. This is because in this case the gas around the Bondi radius can be more easily heated to be above the virial temperature due to its high internal energy. Another question we hope to address is the so-called "sub-Eddington" puzzle. Observationally, the luminosity of almost all AGNs are sub-Eddington, while theoretically the luminosity of an accretion flow can easily be super-Eddington. We find that even when the re-radiation effect is included and outflow does become much stronger, the luminosity, while reduced, can still be super-Eddington.

Our real universe is locally inhomogeneous. Dyer and Roeder introduced the smoothness parameter $\alpha$ to describe the influence of local inhomogeneity on angular diameter distance, and they obtained the angular diameter distance-redshift approximate relation (Dyer-Roeder equation) for locally inhomogeneous universe. Furthermore, the Distance-Duality (DD) relation, $D_L(z)(1+z)^{-2}/D_A(z)=1$, should be valid for all cosmological models that are described by Riemannian geometry, where $D_L$ and $D_A$ are, respectively, the luminosity and angular distance distances. Therefore, it is necessary to test whether if the Dyer-Roeder approximate equation can satisfy the Distance-Duality relation. In this paper, we use Union2.1 SNe Ia data to constrain the smoothness parameter $\alpha$ and test whether the Dyer-Roeder equation satisfies the DD relation. By using $\chi^2$ minimization, we get $\alpha=0.92_{-0.32}^{+0.08}$ at $1\sigma$ and $0.92_{-0.65}^{+0.08}$ at $2\sigma$, and our results show that the Dyer-Roeder equation is in good consistency with the DD relation at $1\sigma$.

Our research objective in this paper is to reconstruct an initial linear density field, which follows the multivariate Gaussian distribution with variances given by the linear power spectrum of the current CDM model and evolves through gravitational instability to the present-day density field in the local Universe. For this purpose, we develop a Hamiltonian Markov Chain Monte Carlo method to obtain the linear density field from a posterior probability function that consists of two components: a prior of a Gaussian density field with a given linear spectrum, and a likelihood term that is given by the current density field. The present-day density field can be reconstructed from galaxy groups using the method developed in Wang et al. (2009a). Using a realistic mock SDSS DR7, obtained by populating dark matter haloes in the Millennium simulation with galaxies, we show that our method can effectively and accurately recover both the amplitudes and phases of the initial, linear density field. To examine the accuracy of our method, we use $N$-body simulations to evolve these reconstructed initial conditions to the present day. The resimulated density field thus obtained accurately matches the original density field of the Millennium simulation in the density range 0.3 <= rho/rho_mean <= 20 without any significant bias. Especially, the Fourier phases of the resimulated density fields are tightly correlated with those of the original simulation down to a scale corresponding to a wavenumber of ~ 1 h/Mpc, much smaller than the translinear scale, which corresponds to a wavenumber of ~ 0.15 h\Mpc.

Using analytical method, we study the Singular Isothermal Quadrupole (SIQ) lens system, which is the simplest lens model that can produce four images. In this case, the radial mass distribution is in accord with the profile of the Singular Isothermal Sphere (SIS) lens, and the tangential distribution is given by adding a quadrupole on the monopole component. The basic properties of the SIQ lens have been studied in this paper, including deflection potential, deflection angle, magnification, critical curve, caustic, pseudo-caustic and transition locus. Analytical solutions of the image positions and magnifications for the source on axes are derived. As have been found, naked cusps will appear when the relative intensity $k$ of quadrupole to monopole is larger than 0.6. According to the magnification invariant theory of the SIQ lens, the sum of the signed magnifications of the four images should be equal to unity \citep{dal98}. However, if a source lies in the naked cusp, the summed magnification of the left three images is smaller than the invariant 1. With this simple lens system, we study the situations that a point source infinitely approaches a cusp or a fold. The sum of magnifications of cusp image triplet is usually not equal to 0, and it is usually positive for major cusp while negative for minor cusp. Similarly, the sum of magnifications of fold image pair is usually neither equal to 0. Nevertheless, the cusp and fold relations are still equal to 0, in that the sum values are divided by infinite absolute magnifications by definition.

Following \cite{dexagol2009} we present a new public code for the fast calculation of null geodesics in the Kerr spacetime. Using Weierstrass' and Jacobi's elliptic functions, we express all coordinates and affine parameters as analytical and numerical functions of a parameter $p$, which is an integral value along the geodesic. This is a main difference of our code compares with previous similar ones. The advantage of this treatment is that the information about the turning points do not need to be specified in advance by the user, and many applications such as imaging, the calculation of line profiles or the observer-emitter problem, etc become root finding problems. All elliptic integrations are computed by Carlson's elliptic integral method as \cite{dexagol2009} did, which guarantees the fast computational speed of our code. The formulae to compute the constants of motion given by \cite{cunnbard1973} have been extended, which allow one readily to handle the situations, in which the emitter or the observer has arbitrary distance and motion state with respect to the central compact object. The validation of the code has been extensively tested by its application to toy problems from the literature. The source FORTRAN code is freely available for download on the web.

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

We study the alignment signal between the distribution of brightest satellite galaxies (BSGs) and the major axis of their host groups using SDSS group catalog constructed by Yang et al. (2007). After correcting for the effect of group ellipticity, a statistically significant (~ 5\sigma) major-axis alignment is detected and the alignment angle is found to be 43.0 \pm 0.4 degrees. More massive and richer groups show stronger BSG alignment. The BSG alignment around blue BCGs is slightly stronger than that around red BCGs. And red BSGs have much stronger major-axis alignment than blue BSGs. Unlike BSGs, other satellites do not show very significant alignment with group major axis. We further explore the BSG alignment in semi-analytic model (SAM) constructed by Guo et al. (2011). We found general good agreement with observations: BSGs in SAM show strong major-axis alignment which depends on group mass and richness in the same way as observations; and none of other satellites exhibit prominent alignment. However, discrepancy also exists in that the SAM shows opposite BSG color dependence, which is most probably induced by the missing of large scale environment ingredient in SAM. The combination of two popular scenarios can explain the detected BSG alignment. The first one: satellites merged into the group preferentially along the surrounding filaments, which is strongly aligned with the major axis of the group. The second one: BSGs enter their host group more recently than other satellites, then will preserve more information about the assembling history and so the major-axis alignment. In SAM, we found positive evidence for the second scenario by the fact that BSGs merged into groups statistically more recently than other satellites. On the other hand, although is opposite in SAM, the BSG color dependence in observation might indicate the first scenario as well.

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.

With galaxy groups constructed from the Sloan Digital Sky Survey (SDSS), we analyze the expected galaxy-galaxy lensing signals around satellite galaxies residing in different host haloes and located at different halo-centric distances. We use Markov Chain Monte Carlo (MCMC) method to explore the potential constraints on the mass and density profile of subhaloes associated with satellite galaxies from SDSS-like surveys and surveys similar to the Large Synoptic Survey Telescope (LSST). Our results show that for SDSS-like surveys, we can only set a loose constraint on the mean mass of subhaloes. With LSST-like surveys, however, both the mean mass and the density profile of subhaloes can be well constrained.

We perform a survey of the Cr, Mn and Fe-K emission lines in young supernova remnants (SNRs) with the Japanese X-ray astronomy satellite {\sl Suzaku}. The Cr and/or Mn emission lines are detected in 3C\,397 and 0519-69.0 for the first time. We also confirm the detection of these lines in Kepler, W49B, N103B and Cas A. We derive the line parameters (i.e., the line centroid energy, flux and equivalent width [EW]) for these six sources and perform a correlation analysis for the line center energies of Cr, Mn and Fe. Also included in the correlation analysis are Tycho and G344.7-0.1 for which the Cr, Mn and Fe-K line parameters were available in the literature through {\sl Suzaku} observations. We find that the line center energies of Cr correlates very well with that of Fe and that of Mn. This confirms our previous findings that the Cr, Mn and Fe are spatially co-located, share a similar ionization state, and have a common origin in the supernova nucleo-synthesis. We find that the ratio of the EW of the Cr emission line to that of Fe ($\EWR\equiv \EWCr/\EWFe$) provides useful constraints on the SNR progenitors and on the SN explosion mechanisms: for SNRs with $\EWR > 2%$, a type Ia origin is favored (e.g., N103B, G344.7-0.1, 3C\,397 and 0519-69.0); for SNRs with $\EWR < 2%$, they could be of either core-collpase origin or carbon-deflagration Ia origin.

The statistics of dark matter halos is an essential component of understanding the nonlinear evolution in modified gravity cosmology. Based on a series of modified gravity N-body simulations, we investigate the halo mass function, concentration and bias. We model the impact of modified gravity by a single parameter \zeta, which determines the enhancement of particle acceleration with respect to GR, given the identical mass distribution (\zeta=1 in GR). We select snapshot redshifts such that the linear matter power spectra of different gravity models are identical, in order to isolate the impact of gravity beyond modifying the linear growth rate. At the baseline redshift corresponding to z_S=1.2 in the standard \Lambda CDM, for a 10% deviation from GR(|\zeta-1|=0.1), the measured halo mass function can differ by about 5-10%, the halo concentration by about 10-20%, while the halo bias differs significantly less. These results demonstrate that the halo mass function and/or the halo concentration are sensitive to the nature of gravity and may be used to make interesting constraints along this line.

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the \Lambda CDM model we used; some results deviate from expectation at a level of ~ 3\sigma but the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.

The pulsating DA white dwarfs (ZZ Ceti stars) are $g$-mode non-radial pulsators. Asteroseismology provides strong constraints on their global parameters and internal structure. Since all the DA white dwarfs falling in the ZZ Ceti instability strip do pulsate, the internal structure derived from asteroseismology brings knowledge for the DA white dwarfs as a whole group. HS 0507+0434B is one of the ZZ Ceti stars which lies approximately in the middle of the instability strip for which we have undertaken a detailed asteroseismological study. We carried out multisite observation campaigns in 2007 and from December 2009 to January 2010. In total, 206 hours of photometric time-series have been collected. They have been analysed by means of Fourier analysis and simultaneous multi-frequency sine-wave fitting. In total, 39 frequency values are resolved including 6 triplets and a number of linear combinations. We identify the triplets as $\ell$=1 $g$-modes split by rotation. We derived the period spacing, the rotational splitting and the rotation rate. From the comparison of the observed periods with the theoretical periods of a series of models we estimate the fundamental parameters of the star: its total mass M$_{*}$/M$_{\odot}$ = 0.675, its luminosity L/L$_{\odot}$=3.5$\times 10^{-3}$, and its hydrogen mass fraction M$_{H}$/M$_{*}$= 10$^{-8.5}$.

Peaks in two-dimensional weak lensing (WL) maps contain significant cosmological information, complementary to the WL power spectrum. This has recently been demonstrated using N-body simulations which neglect baryonic effects. Here we employ ray-tracing N-body simulations in which we manually steepen the density profile of each dark matter halo, mimicking the cooling and concentration of baryons into dark matter potential wells. We find, in agreement with previous works, that this causes a significant increase in the amplitude of the WL power spectrum on small scales (spherical harmonic index l>1,000). We then study the impact of the halo concentration increase on the peak counts, and find the following. (i) Low peaks (with convergence 0.02 < kappa_peak < 0.08), remain nearly unaffected. These peaks are created by a constellation of several halos with low masses (10^12-10^13 M_sun) and large angular offsets from the peak center (> 0.5 R_vir); as a result, they are insensitive to the central halo density profiles. These peaks contain most of the cosmological information, and thus provide an unusually sensitive and unbiased probe. (ii) The number of high peaks (with convergence kappa_peak > 0.08) is increased. However, when the baryon effects are neglected in cosmological parameter estimation, then the high peaks lead to a modest bias, comparable to that from the power spectrum on relatively large-scales (l<2000), and much smaller than the bias from the power spectrum on smaller scales (l>2,000). (iii) In the 3D parameter space (sigma_8, Omega_m, w), the biases from the high peaks and the power spectra are in different directions. This suggests the possibility of "self-calibration": the combination of peak counts and power spectrum can simultaneously constrain baryonic physics and cosmological parameters.

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.

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space cross-correlation function (CCF) with respect to a reference galaxy sample, xi(r_p, pi), the projected CCF, w_p(r_p), and the real-space CCF, xi(r). The VDP is then extracted from the redshift distortion in xi(r_p, pi), by comparing xi(r_p, pi) with xi(r). We find that the velocity dispersion (VD) within virial radius (R_200) shows a roughly flat profile, with a slight increase at radii below ~0.3 R_200 for high mass systems. The average VD within the virial radius, sigma_v, is a strongly increasing function of central galaxy mass. We apply the same methodology to N-body simulations with the concordance Lambda cold dark matter cosmology but different values of the density fluctuation parameter sigma_8, and we compare the results to the SDSS results. We show that the sigma_v-M_* relation from the data provides stringent constraints on both sigma_8 and sigma_ms, the dispersion in log M_* of central galaxies at fixed halo mass. Our best-fitting model suggests sigma_8 = 0.86 +/- 0.03 and sigma_ms = 0.16 +/- 0.03. The slightly higher value of sigma_8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for central galaxies of different luminosities and masses, in good agreement with the results obtained by Mandelbaum et al. (2006) from stacking the gravitational lensing signals of the SDSS galaxies.

We simultaneously constrain cosmology and galaxy bias using measurements of galaxy abundances, galaxy clustering and galaxy-galaxy lensing taken from the Sloan Digital Sky Survey. We use the conditional luminosity function (which describes the halo occupation statistics as function of galaxy luminosity) combined with the halo model (which describes the non-linear matter field in terms of its halo building blocks) to describe the galaxy-dark matter connection. We explicitly account for residual redshift space distortions in the projected galaxy-galaxy correlation functions, and marginalize over uncertainties in the scale dependence of the halo bias and the detailed structure of dark matter haloes. Under the assumption of a spatially flat, vanilla {\Lambda}CDM cosmology, we focus on constraining the matter density, {\Omega}m, and the normalization of the matter power spectrum, {\sigma}8, and we adopt WMAP7 priors for the spectral index, the Hubble parameter, and the baryon density. We obtain that \Omegam = 0.278_{-0.026}^{+0.023} and {\sigma}8 = 0.763_{-0.049}^{+0.064} (95% CL). These results are robust to uncertainties in the radial number density distribution of satellite galaxies, while allowing for non-Poisson satellite occupation distributions results in a slightly lower value for {\sigma}8 (0.744_{-0.047}^{+0.056}). These constraints are in excellent agreement (at the 1{\sigma} level) with the cosmic microwave background constraints from WMAP. This demonstrates that the use of a realistic and accurate model for galaxy bias, down to the smallest non-linear scales currently observed in galaxy surveys, leads to results perfectly consistent with the vanilla {\Lambda}CDM cosmology.

We quantify the accuracy with which the cosmological parameters characterizing the energy density of matter (\Omega_m), the amplitude of the power spectrum of matter fluctuations (\sigma_8), the energy density of neutrinos (\Omega_{\nu}) and the dark energy equation of state (w_0) can be constrained using data from large galaxy redshift surveys. We advocate a joint analysis of the abundance of galaxies, galaxy clustering, and the galaxy-galaxy weak lensing signal in order to simultaneously constrain the halo occupation statistics (i.e., galaxy bias) and the cosmological parameters of interest. We parameterize the halo occupation distribution of galaxies in terms of the conditional luminosity function and use the analytical framework of the halo model described in our companion paper (van den Bosch et al. 2012), to predict the relevant observables. By performing a Fisher matrix analysis, we show that a joint analysis of these observables, even with the precision with which they are currently measured from the Sloan Digital Sky Survey, can be used to obtain tight constraints on the cosmological parameters, fully marginalized over uncertainties in galaxy bias. We demonstrate that the cosmological constraints from such an analysis are nearly uncorrelated with the halo occupation distribution constraints, thus, minimizing the systematic impact of any imperfections in modeling the halo occupation statistics on the cosmological constraints. In fact, we demonstrate that the constraints from such an analysis are both complementary to and competitive with existing constraints on these parameters from a number of other techniques, such as cluster abundances, cosmic shear and/or baryon acoustic oscillations, thus paving the way to test the concordance cosmological model.

We present a new method that simultaneously solves for cosmology and galaxy bias on non-linear scales. The method uses the halo model to analytically describe the (non-linear) matter distribution, and the conditional luminosity function (CLF) to specify the halo occupation statistics. For a given choice of cosmological parameters, this model can be used to predict the galaxy luminosity function, as well as the two-point correlation functions of galaxies, and the galaxy-galaxy lensing signal, both as function of scale and luminosity. In this paper, the first in a series, we present the detailed, analytical model, which we test against mock galaxy redshift surveys constructed from high-resolution numerical $N$-body simulations. We demonstrate that our model, which includes scale-dependence of the halo bias and a proper treatment of halo exclusion, reproduces the 3-dimensional galaxy-galaxy correlation and the galaxy-matter cross-correlation (which can be projected to predict the observables) with an accuracy better than 10 (in most cases 5) percent. Ignoring either of these effects, as is often done, results in systematic errors that easily exceed 40 percent on scales of $\sim 1 h^{-1}\Mpc$, where the data is typically most accurate. Finally, since the projected correlation functions of galaxies are never obtained by integrating the redshift space correlation function along the line-of-sight out to infinity, simply because the data only cover a finite volume, they are still affected by residual redshift space distortions (RRSDs). Ignoring these, as done in numerous studies in the past, results in systematic errors that easily exceed 20 perent on large scales ($r_\rmp \gta 10 h^{-1}\Mpc$). We show that it is fairly straightforward to correct for these RRSDs, to an accuracy better than $\sim 2$ percent, using a mildly modified version of the linear Kaiser formalism.

It is an indisputable fact that solar magnetic fields are force-free in the corona, where force free fields means that current and magnetic fields are parallel and there is no Lorentz force in the fields. While the force-free extent of photospheric magnetic fields remains open. In this paper, the statistical results about it is given. The vector magnetograms (namely, $B_{x}$, $B_{y}$ and $B_{z}$ in heliocentric coordinates) are employed, which are deduced and calibrated from Stokes spectra, observed by Solar Magnetic Field Telescope (SMFT) at Huairou Solar Observing Station (HSOS) are used. We study and calibrated 925 magnetograms calibrated by two sets of calibration coefficients, that indicate the relations between magnetic fields and the strength of Stokes spectrum and can be calculated either theoretically or empirically. The statistical results show that the majority of active region magnetic fields are not consistent with the force-free model.

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

A statistical study is carried out on the photospheric magnetic nonpotentiality in solar active regions and its relationship with associated flares. We select 2173 photospheric vector magnetograms from 1106 active regions observed by the Solar Magnetic Field Telescope at Huairou Solar Observing Station, National Astronomical Observatories of China, in the period of 1988-2008, which covers most of the 22nd and 23rd solar cycles. We have computed the mean planar magnetic shear angle (\bar{\Delta\phi}), mean shear angle of the vector magnetic field (\bar{\Delta\psi}), mean absolute vertical current density (\bar{|J_{z}|}), mean absolute current helicity density (\bar{|h_{c}|}), absolute twist parameter (|\alpha_{av}|), mean free magnetic energy density (\bar{\rho_{free}}), effective distance of the longitudinal magnetic field (d_{E}), and modified effective distance (d_{Em}) of each photospheric vector magnetogram. Parameters \bar{|h_{c}|}, \bar{\rho_{free}}, and d_{Em} show higher correlation with the evolution of the solar cycle. The Pearson linear correlation coefficients between these three parameters and the yearly mean sunspot number are all larger than 0.59. Parameters \bar{\Delta\phi}, \bar{\Delta\psi}, \bar{|J_{z}|}, |\alpha_{av}|, and d_{E} show only weak correlations with the solar cycle, though the nonpotentiality and the complexity of active regions are greater in the activity maximum periods than in the minimum periods. All of the eight parameters show positive correlations with the flare productivity of active regions, and the combination of different nonpotentiality parameters may be effective in predicting the flaring probability of active regions.

We present a method that uses observations of galaxies to simultaneously constrain cosmological parameters and the galaxy-dark matter connection (aka halo occupation statistics). The latter describes how galaxies are distributed over dark matter haloes, and is an imprint of the poorly understood physics of galaxy formation. A generic problem of using galaxies to constrain cosmology is that galaxies are a biased tracer of the mass distribution, and this bias is generally unknown. The great advantage of simultaneously constraining cosmology and halo occupation statistics is that this effectively allows cosmological constraints marginalized over the uncertainties regarding galaxy bias. Not only that, it also yields constraints on the galaxy-dark matter connection, this time properly marginalized over cosmology, which is of great value to inform theoretical models of galaxy formation. We use a combination of the analytical halo model and the conditional luminosity function to describe the galaxy-dark matter connection, which we use to model the abundance, clustering and galaxy-galaxy lensing properties of the galaxy population. We use a Fisher matrix analysis to gauge the complementarity of these different observables, and present some preliminary results from an analysis based on data from the Sloan Digital Sky Survey. Our results are complementary to and perfectly consistent with the results from the 7 year data release of the WMAP mission, strengthening the case for a true 'concordance' cosmology.

We present a new model to describe the galaxy-dark matter connection across cosmic time, which unlike the popular subhalo abundance matching technique is self-consistent in that it takes account of the facts that (i) subhalos are accreted at different times, and (ii) the properties of satellite galaxies may evolve after accretion. Using observations of galaxy stellar mass functions out to $z \sim 4$, the conditional stellar mass function at $z\sim 0.1$ obtained from SDSS galaxy group catalogues, and the two-point correlation function (2PCF) of galaxies at $z \sim 0.1$ as function of stellar mass, we constrain the relation between galaxies and dark matter halos over the entire cosmic history from $z \sim 4$ to the present. This relation is then used to predict the median assembly histories of different stellar mass components within dark matter halos (central galaxies, satellite galaxies, and halo stars). We also make predictions for the 2PCFs of high-$z$ galaxies as function of stellar mass. Our main findings are the following: (i) Our model reasonably fits all data within the observational uncertainties, indicating that the $\Lambda$CDM concordance cosmology is consistent with a wide variety of data regarding the galaxy population across cosmic time. (ii) ... [abridged]

We present a concept study on weak lensing map reconstruction through the cosmic magnification effect in galaxy number density distribution. We propose a minimal variance linear estimator to minimize both the dominant systematical and statistical errors in the map reconstruction. It utilizes the distinctively different flux dependences to separate the cosmic magnification signal from the overwhelming galaxy intrinsic clustering noise. It also minimizes the shot noise error by an optimal weighting scheme on the galaxy number density in each flux bin. Our method is in principle applicable to all galaxy surveys with reasonable redshift information. We demonstrate its applicability against the planned Square Kilometer Array survey, under simplified conditions. Weak lensing maps reconstructed through our method are complementary to that from cosmic shear and CMB and 21cm lensing. They are useful for cross checking over systematical errors in weak lensing reconstruction and for improving cosmological constraints.

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 use the ROSAT all sky survey X-ray cluster catalogs and the optical SDSS DR7 galaxy and group catalogs to cross-identify X-ray clusters with their optical counterparts, resulting in a sample of 201 X-ray clusters in the sky coverage of SDSS DR7. We investigate various correlations between the optical and X-ray properties of these X-ray clusters, and find that the following optical properties are correlated with the X-ray luminosity: the central galaxy luminosity, the central galaxy mass, the characteristic group luminosity ($\propto \Lx^{0.43}$), the group stellar mass ($\propto \Lx^{0.46}$), with typical 1-$\sigma$ scatter of $\sim 0.67$ in $\log \Lx$. Using the observed number distribution of X-ray clusters, we obtain an unbiased scaling relation between the X-ray luminosity, the central galaxy stellar mass and the characteristic satellite stellar mass as ${\log L_X} = -0.26 + 2.90 [\log (M_{\ast, c} + 0.26 M_{\rm sat}) -12.0]$ (and in terms of luminosities, as ${\log L_X} = -0.15 + 2.38 [\log (L_{c} + 0.72 L_{\rm sat}) -12.0]$). We find that the systematic difference between different halo mass estimations, e.g., using the ranking of characteristic group stellar mass or using the X-ray luminosity scaling relation can be used to constrain cosmology. Comparing the properties of groups of similar stellar mass (or optical luminosities) and redshift that are X-ray luminous or under-luminous, we find that X-ray luminous groups have more faint satellite galaxies and higher red fraction in their satellites. The cross-identified X-ray clusters together with their optical properties are provided in Appendix B.

[abridge]Cosmic velocity and tidal fields are important for the understanding of the cosmic web and the environments of galaxies, and can also be used to constrain cosmology. In this paper, we reconstruct these two fields in SDSS volume from dark matter halos represented by galaxy groups. Detailed mock catalogues are used to test the reliability of our method against uncertainties arising from redshift distortions, survey boundaries, and false identifications of groups by our group finder. We find that both the velocity and tidal fields, smoothed on a scale of ~2Mpc/h, can be reliably reconstructed in the inner region (~66%) of the survey volume. The reconstructed tidal field is used to split the cosmic web into clusters, filaments, sheets, and voids, depending on the sign of the eigenvalues of tidal tensor. The reconstructed velocity field nicely shows how the flows are diverging from the centers of voids, and converging onto clusters, while sheets and filaments have flows that are convergent along one and two directions, respectively. We use the reconstructed velocity field and the Zel'dovich approximation to predict the mass density field in the SDSS volume as function of redshift, and find that the mass distribution closely follows the galaxy distribution even on small scales. We find a large-scale bulk flow of about 117km/s in a very large volume, equivalent to a sphere with a radius of ~170Mpc/h, which seems to be produced by the massive structures associated with the SDSS Great Wall. Finally, we discuss potential applications of our reconstruction to study the environmental effects of galaxy formation, to generate initial conditions for simulations of the local Universe, and to constrain cosmological models. The velocity, tidal and density fields in the SDSS volume, specified on a Cartesian grid with a spatial resolution of ~700kpc/h, are available from the authors upon request.

Recent studies have shown that the number counts of convergence peaks N(kappa) in weak lensing (WL) maps, expected from large forthcoming surveys, can be a useful probe of cosmology. We follow up on this finding, and use a suite of WL convergence maps, obtained from ray-tracing N-body simulations, to study (i) the physical origin of WL peaks with different heights, and (ii) whether the peaks contain information beyond the convergence power spectrum P_ell. In agreement with earlier work, we find that high peaks (with amplitudes >~ 3.5 sigma, where sigma is the r.m.s. of the convergence kappa) are typically dominated by a single massive halo. In contrast, medium-height peaks (~0.5-1.5 sigma) cannot be attributed to a single collapsed dark matter halo, and are instead created by the projection of multiple (typically, 4-8) halos along the line of sight, and by random galaxy shape noise. Nevertheless, these peaks dominate the sensitivity to the cosmological parameters w, sigma_8, and Omega_m. We find that the peak height distribution and its dependence on cosmology differ significantly from predictions in a Gaussian random field. We directly compute the marginalized errors on w, sigma_8, and Omega_m from the N(kappa) + P_ell combination, including redshift tomography with source galaxies at z_s=1 and z_s=2. We find that the N(kappa) + P_ell combination has approximately twice the cosmological sensitivity compared to P_ell alone. These results demonstrate that N(kappa) contains non-Gaussian information complementary to the power spectrum.

An analytical model is developed for the mass function of cold dark matter subhalos at the time of accretion and for the distribution of their accretion times. Our model is based on the model of Zhao et al. (2009) for the median assembly histories of dark matter halos, combined with a simple log-normal distribution to describe the scatter in the main-branch mass at a given time for halos of the same final mass. Our model is simple, and can be used to predict the un-evolved subhalo mass function, the mass function of subhalos accreted at a given time, the accretion-time distribution of subhalos of a given initial mass, and the frequency of major mergers as a function of time. We test our model using high-resolution cosmological $N$-body simulations, and find that our model predictions match the simulation results remarkably well. Finally, we discuss the implications of our model for the evolution of subhalos in their hosts and for the construction of a self-consistent model to link galaxies and dark matter halos at different cosmic times.

It is generally assumed that the central galaxy in a dark matter halo, that is, the galaxy with the lowest specific potential energy, is also the brightest halo galaxy (BHG), and that it resides at rest at the centre of the dark matter potential well. This central galaxy paradigm (CGP) is an essential assumption made in various fields of astronomical research. In this paper we test the validity of the CGP using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. For each group we compute two statistics, ${\cal R}$ and ${\cal S}$, which quantify the offsets of the line-of-sight velocities and projected positions of brightest group galaxies relative to the other group members. By comparing the cumulative distributions of $|{\cal R}|$ and $|{\cal S}|$ to those obtained from detailed mock group catalogues, we rule out the null-hypothesis that the CGP is correct. Rather, the data indicate that in a non-zero fraction $f_{\rm BNC}(M)$ of all haloes of mass $M$ the BHG is not the central galaxy, but instead, a satellite galaxy. In particular, we find that $f_{\rm BNC}$ increases from $\sim 0.25$ in low mass haloes ($10^{12} h^{-1} {\rm M_{\odot}} \leq M \lsim 2 \times 10^{13} h^{-1}{\rm M_{\odot}}$) to $\sim 0.4$ in massive haloes ($M \gsim 5 \times 10^{13} h^{-1} {\rm M_{\odot}}$). We show that these values of $f_{\rm BNC}$ are uncomfortably high compared to predictions from halo occupation statistics and from semi-analytical models of galaxy formation. We end by discussing various implications of a non-zero $f_{\rm BNC}(M)$, with an emphasis on the halo masses inferred from satellite kinematics.

It has been a long-standing question whether fossil groups are just sampling the tail of the distribution of ordinary groups, or whether they are a physically distinct class of objects, characterized by an unusual and special formation history. To study this question, we here investigate fossil groups identified in the hydrodynamical simulations of the GIMIC project, which consists of resimulations of five regions in the Millennium Simulation (MS) that are characterized by different large-scale densities, ranging from a deep void to a proto-cluster region. For comparison, we also consider semi-analytic models built on top of the MS, as well as a conditional luminosity function approach. We identify galaxies in the GIMIC simulations as groups of stars and use a spectral synthesis code to derive their optical properties. The X-ray luminosity of the groups is estimated in terms of the thermal bremsstrahlung emission of the gas in the host halos, neglecting metallicity effects. We focus on comparing the properties of fossil groups in the theoretical models and observational results, highlighting the differences between them, and trying to identify possible dependencies on environment for which our approach is particularly well set-up. We find that the optical fossil fraction in all of our theoretical models declines with increasing halo mass, and there is no clear environmental dependence. Combining the optical and X-ray selection criteria for fossil groups, the halo mass dependence of the fossil groups seen in optical vanishes. Over the GIMIC halo mass range we resolve best, 9.0\times1012 \sim 4.0\times1013 h-1 M, the central galaxies in the fossil groups show similar properties as those in ordinary groups, in terms of age, metallicity, color, concentration, and mass-to-light ratio. [abridged]

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.

Cosmic magnification has been detected through cross correlation between foreground and background populations (galaxies or quasars). It has been shown that weighing each background object by its $\alpha-1$ can significantly improve the cosmic magnification measurement \citep{Menard02,Scranton05}. Here, $\alpha$ is the logarithmic slope of the luminosity function of background populations. However, we find that this weighting function is optimal only for sparse background populations in which intrinsic clustering is negligible with respect to shot noise. We derive the optimal weighting function for general case including scale independent and scale dependent weights. The optimal weighting function improves the S/N (signal to noise ratio) by $\sim 20%$ for a BigBOSS-like survey and the improvement can reach a factor of $\sim 2$ for surveys with much denser background populations.

We propose that a non-thermal X-ray arc inside the remnant of Tycho's supernova (SN) represents the interaction between the SN ejecta and the companion star's envelope lost in the impact of the explosion. The X-ray emission of the remnant further shows an apparent shadow casted by the arc in the opposite direction of the explosion site, consistent with the blocking of the SN ejecta by the envelope. This scenario supports the single degenerate binary origin of Tycho's SN. The properties of the X-ray arc, together with the previous detection of the companion candidate and its space velocity by Ruiz-Lapuente et al. (2004) and Hern\'{a}ndez et al.(2009), enables us to further infer 1) the progenitor binary has a period of 4.9$^{+5.3}_{-3.0}$ days, 2) the companion gained a kick velocity of 42$\pm$30 km s$^{-1}$, and 3) the stripped envelope mass is about 0.0016($\leq0.0083$) $M_{\sun}$. However, we notice that the nature of the companion candidate is still under debate, and the above parameters need to be revised according to the actual properties of the companion candidate. Further work to measure the proper motion of the arc and to check the capability of the interaction to emit the amount of X-rays observed from the arc is also needed to validate the current scenario.

We use seven high-resolution $N$-body simulations to study the correlations among different halo properties (assembly time, spin, shape and substructure), and how these halo properties are correlated with the large-scale environment in which halos reside. The large-scale tidal field estimated from halos above a mass threshold is used as our primary quantity to characterize the large-scale environment, while other parameters, such as the local overdensity and the morphology of large-scale structure, are used for comparison. For halos at a fixed mass, all the halo properties depend significantly on environment, particularly the tidal field. The environmental dependence of halo assembly time is primarily driven by local tidal field. The mass of the unbound fraction in substructure is boosted in strong tidal force region, while the bound fraction is suppressed. Halos have a tendency to spin faster in stronger tidal field and the trend is stronger for more massive halos. The spin vectors show significant alignment with the intermediate axis of the tidal field, as expected from the tidal torque theory. Both the major and minor axes of halos are strongly aligned with the corresponding principal axes of the tidal field. In general, a halo that can accrete more material after the formation of its main halo on average is younger, is more elongated, spins faster, and contains a larger amount of substructure. Higher density environments not only provide more material for halo to accrete, but also are places of stronger tidal field that tends to suppress halo accretion. The environmental dependencies are the results of these two competing effects. The tidal field based on halos can be estimated from observation, and we discuss the implications of our results for the environmental dependence of galaxy properties.

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 investigate the potential of exploiting the Sunyaev-Zeldovich effect (SZE) to study the properties of hot gas in galaxy groups. It is shown that, with upcoming SZE surveys, one can stack SZE maps around galaxy groups of similar halo masses selected from large galaxy redshift surveys to study the hot gas in halos represented by galaxy groups. We use various models for the hot halo gas to study how the expected SZE signals are affected by gas fraction, equation of state, halo concentration, and cosmology. Comparing the model predictions with the sensitivities expected from the SPT, ACT and Planck surveys shows that a SPT-like survey can provide stringent constraints on the hot gas properties for halos with masses M ~> 10^{13} h^{-1}Msun. We also explore the idea of using the cross correlation between hot gas and galaxies of different luminosity to probe the hot gas in dark matter halos without identifying galaxy groups to represent dark halos. Our results show that, with a galaxy survey as large as the Sloan Digital Sky Survey and with the help of the conditional luminosity function (CLF) model, one can obtain stringent constraints on the hot gas properties in halos with masses down to 10^{13} h^{-1}Msun. Thus, the upcoming SZE surveys should provide a very promising avenue to probe the hot gas in relatively low-mass halos where the majority of L*-galaxies reside.

We present an {\sl XMM-Newton} observation of the eclipsing binary Algol which contains an X-ray dark B8V primary and an X-ray bright K2IV secondary. The observation covered the optical secondary eclipse and captured an X-ray flare that was eclipsed by the B star. The EPIC and RGS spectra of Algol in its quiescent state are described by a two-temperature plasma model. The cool component has a temperature around 6.4$\times 10^{6}$ K while that of the hot component ranges from 2 to 4.0$\times 10^{7}$ K. Coronal abundances of C, N, O, Ne, Mg, Si and Fe were obtained for each component for both the quiescent and the flare phases, with generally upper limits for S and Ar, and C, N, and O for the hot component. F-tests show that the abundances need not to be different between the cool and the hot component and between the quiescent and the flare phase with the exception of Fe. Whereas the Fe abundance of the cool component remains constant at $\sim$0.14, the hot component shows an Fe abundance of $\sim$0.28, which increases to $\sim$0.44 during the flare. This increase is expected from the chromospheric evaporation model. The absorbing column density $N_H$ of the quiescent emission is 2.5$\times10^{20}$ cm$^{-2}$, while that of the flare-only emission is significantly lower and consistent with the column density of the interstellar medium. This observation substantiates earlier suggestions of the presence of X-ray absorbing material in the Algol system.

We study the topology of cosmic large-scale structure through the genus statistics, using galaxy catalogues generated from the Millennium Simulation and observational data from the latest Sloan Digital Sky Survey Data Release (SDSS DR7). We introduce a new method for constructing galaxy density fields and for measuring the genus statistics of its isodensity surfaces. It is based on a Delaunay tessellation field estimation (DTFE) technique that allows the definition of a piece-wise continuous density field and the exact computation of the topology of its polygonal isodensity contours, without introducing any free numerical parameter. Besides this new approach, we also employ the traditional approaches of smoothing the galaxy distribution with a Gaussian of fixed width, or by adaptively smoothing with a kernel that encloses a constant number of neighboring galaxies. Our results show that the Delaunay-based method extracts the largest amount of topological information. Unlike the traditional approach for genus statistics, it is able to discriminate between the different theoretical galaxy catalogues analyzed here, both in real space and in redshift space, even though they are based on the same underlying simulation model. In particular, the DTFE approach detects with high confidence a discrepancy of one of the semi-analytic models studied here compared with the SDSS data, while the other models are found to be consistent.

A large fraction of cosmological information on dark energy and gravity is encoded in the nonlinear regime. Precision cosmology thus requires precision modeling of nonlinearities in general dark energy and modified gravity models. We modify the Gadget-2 code and run a series of N-body simulations on modified gravity cosmology to study the nonlinearities. The modified gravity model that we investigate in the present paper is characterized by a single parameter \zeta, which determines the enhancement of particle acceleration with respect to general relativity (GR), given the identical mass distribution (\zeta = 1 in GR). The first nonlinear statistics we investigate is the nonlinear matter power spectrum at k < 3h/Mpc, which is the relevant range for robust weak lensing power spectrum modeling at l < 2000. In this study, we focus on the relative difference in the nonlinear power spectra at corresponding redshifts where different gravity models have the same linear power spectra. This particular statistics highlights the imprint of modified gravity in the nonlinear regime and the importance to include the nonlinear regime in testing GR. By design, it is less susceptible to the sample variance and numerical artifacts. We adopt a mass assignment method based on wavelet to improve the power spectrum measurement. We run a series of tests to determine the suitable simulation specifications (particle number, box size and initial redshift). We find that, the nonlinear power spectra can differ by ~30% for 10% deviation from GR (|\zeta-1| = 0.1) where the rms density fluctuations reach 10. This large difference, on one hand, shows the richness of information on gravity in the corresponding scales, and on the other hand, invalidates simple extrapolations of some existing fitting formulae to modified gravity cosmology.

We present preliminary observational results of the multi-site STEPHI campaign on the cluster NGC 1817. The three observatories involved are San Pedro Martir (Mexico), Xing Long (China) and the Observatorio del Teide (Spain) - giving an ideal combination to maximise the duty cycle. The cluster has 12 known delta Scuti stars and at least two detached eclipsing binary systems. This combination of characteristics is ideal for extracting information about global parameters of the targets, which will in turn impose strict constraints on the stellar models. From an initial comparison with stellar models using the known fundamental parameters, and just the observed pulsation frequencies and measured effective temperatures, it appears that a lower value of initial helium mass fraction will most likely explain the observations of these stars.

The variations of gravity were measured with a high precision LaCoste-Romberg D gravimeter during a total solar eclipse to investigate the effect of solar eclipse on the gravitational field. The observed anomaly $(7.0 \pm 2.7) \times 10^{-8}$ m/s$^2$ during the eclipse implies that there may be a shielding property of gravitation.

We use the kinematics of satellite galaxies that orbit around the central galaxy in a dark matter halo to infer the scaling relations between halo mass and central galaxy properties. Using galaxies from the Sloan Digital Sky Survey, we investigate the halo mass-luminosity relation (MLR) and the halo mass-stellar mass relation (MSR) of central galaxies. In particular, we focus on the dependence of these scaling relations on the colour of the central galaxy. We find that red central galaxies on average occupy more massive haloes than blue central galaxies of the same luminosity. However, at fixed stellar mass there is no appreciable difference in the average halo mass of red and blue centrals, especially for M* $\lsim$ 10^{10.5} h^{-2} Msun. This indicates that stellar mass is a better indicator of halo mass than luminosity. Nevertheless, we find that the scatter in halo masses at fixed stellar mass is non-negligible for both red and blue centrals. It increases as a function of stellar mass for red centrals but shows a fairly constant behaviour for blue centrals. We compare the scaling relations obtained in this paper with results from other independent studies of satellite kinematics, with results from a SDSS galaxy group catalog, from galaxy-galaxy weak lensing measurements, and from subhalo abundance matching studies. Overall, these different techniques yield MLRs and MSRs in fairly good agreement with each other (typically within a factor of two), indicating that we are converging on an accurate and reliable description of the galaxy-dark matter connection. We briefly discuss some of the remaining discrepancies among the various methods.

We compare the stellar masses of central and satellite galaxies predicted by three independent semianalytical models with observational results obtained from a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. In particular, we compare the stellar mass functions of centrals and satellites, the relation between total stellar mass and halo mass, and the conditional stellar mass functions, which specify the average number of galaxies of stellar mass M_* that reside in a halo of mass M_h. The semi-analytical models only predict the correct stellar masses of central galaxies within a limited mass range and all models fail to reproduce the sharp decline of stellar mass with decreasing halo mass observed at the low mass end. In addition, all models over-predict the number of satellite galaxies by roughly a factor of two. The predicted stellar mass in satellite galaxies can be made to match the data by assuming that a significant fraction of satellite galaxies are tidally stripped and disrupted, giving rise to a population of intra-cluster stars in their host halos. However, the amount of intra-cluster stars thus predicted is too large compared to observation. This suggests that current galaxy formation models still have serious problems in modeling star formation in low-mass halos.

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.

We present new results of the HESS J1731-347/SNR G353.6-0.7 system from XMM-NEWTOM and Suzaku X-ray observations, and Delinha CO observations. We discover extended hard X-rays coincident with the bright, extended TeV source HESS J1731-347 and the shell of the radio SNR. We find that spatially-resolved X-ray spectra can generally be characterized by an absorbed power-law model, with photon-index of ~ 2, typical of non-thermal emission. A bright X-ray compact source, XMMS J173203-344518, is also detected near the center of the SNR. We find no evidence of a radio counterpart or an extended X-ray morphology for this source, making it unlikely to be a pulsar wind nebular (PWN). The spectrum of the source can be well fitted by an absorbed blackbody with a temperature of ~ 0.5 keV plus a power-law tail with a photon-index of ~ 5, reminiscent of the X-ray emission of a magnetar. CO observations toward the inner part of the HESS source reveal a bright cloud component at -20+/-4 km s^{-1}, which is likely located at the same distance of ~ 3.2 kpc as the SNR. Based on the probable association between the X-ray and $\gamma$-ray emissions and likely association between the CO cloud and the SNR, we argue that the extended TeV emission originates from the interaction between the SNR shock and the adjacent CO clouds rather than from a PWN.

Using a large SDSS galaxy group catalogue, we study how the stellar ages and metallicities of central and satellite galaxies depend on stellar mass and halo mass. We find that satellites are older and metal-richer than centrals of the same stellar mass. In addition, the slopes of the age-stellar mass and metallicity-stellar mass relations are found to become shallower in denser environments. This is due to the fact that the average age and metallicity of low mass satellite galaxies increase with the mass of the halo in which they reside. A comparison with the semi-analytical model of Wang et al. (2008) shows that it succesfully reproduces the fact that satellites are older than centrals of the same stellar mass and that the age difference increases with the halo mass of the satellite. This is a consequence of strangulation, which leaves the stellar populations of satellites to evolve passively, while the prolonged star formation activity of centrals keeps their average ages younger. The resulting age offset is larger in more massive environments because their satellites were accreted earlier. The model fails, however, in reproducing the halo mass dependence of the metallicities of low mass satellites, yields metallicity-stellar mass and age-stellar mass relations that are too shallow, and predicts that satellite galaxies have the same metallicities as centrals of the same stellar mass, in disagreement with the data. We argue that these discrepancies are likely to indicate the need to (i) modify the recipes of both supernova feedback and AGN feedback, (ii) use a more realistic description of strangulation, and (iii) include a proper treatment of the tidal stripping, heating and destruction of satellite galaxies. [Abridged]

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.

Clusters, filaments, sheets and voids are the building blocks of the cosmic web. In this study, we present and compare two distinct algorithms for finding cosmic filaments and sheets, a task which is far less well established than the identification of dark matter halos or voids. One method is based on the smoothed dark matter density field, the other uses the halo distributions directly. We apply both techniques to one high resolution N-body simulation and reconstruct the filamentary/sheet like network of the dark matter density field. We focus on investigating the properties of the dark matter halos inside these structures, in particular on the directions of their spins and the orientation of their shapes with respect to the directions of the filaments and sheets. We find that both the spin and the major axes of filament-halos with masses <= 10^{13} M_sun/h are preferentially aligned with the direction of the filaments. The spins and major axes of halos in sheets tend to lie parallel to the sheets. There is an opposite mass dependence of the alignment strengths for the spin (negative) and major (positive) axes, i.e. with increasing halo mass the major axis tends to be more strongly aligned with the direction of the filament whereas the alignment between halo spin and filament becomes weaker with increasing halo mass. The alignment strengths as a function of distance to the most massive node halo indicate that there is a transit large scale environment impact: from the 2-D collapse phase of the filament to the 3-D collapse phase of the cluster/node halo at small separation. Overall, the two algorithms for filament/sheet identification investigated here agree well with each other. The method based on halos alone can be easily adapted for use with observational data sets.

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.

Using a sample of galaxy groups found in the Sloan Digital Sky Survey Data Release 4, we measure the following four types of alignment signals: (1) the alignment between the distributions of the satellites of each group relative to the direction of the nearest neighbor group (NNG); (2) the alignment between the major axis direction of the central galaxy of the host group (HG) and the direction of the NNG; (3) the alignment between the major axes of the central galaxies of the HG and the NNG; and (4) the alignment between the major axes of the satellites of the HG and the direction of the NNG. We find strong signal of alignment between the satellite distribution and the orientation of central galaxy relative to the direction of the NNG, even when the NNG is located beyond $3r_{\rm vir}$ of the host group. The major axis of the central galaxy of the HG is aligned with the direction of the NNG. The alignment signals are more prominent for groups that are more massive and with early type central galaxies. We also find that there is a preference for the two major axis of the central galaxies of the HG and NNG to be parallel for the system with both early central galaxies, however not for the systems with both late type central galaxies. For the orientation of satellite galaxies, we do not find any significant alignment signals relative to the direction of the NNG. From these four types of alignment measurements, we conclude that the large scale environment traced by the nearby group affects primarily the shape of the host dark matter halo, and hence also affects the distribution of satellite galaxies and the orientation of central galaxies. (Shortened version)

Pictor A is a nearby Fanaroff-Riley class II (FR II) radio galaxy with a bright hot spot, the western hot spot. Observation of high polarization in the optical emission of the hot spot indicates that the optical emission could be synchrotron radiation of relativistic electrons in the hot spot. These electrons may be able to produce high energy gamma-ray photons through inverse Compton (IC) scattering. We use single-zone and multi-zone synchrotron + synchrotron-self-Compton (SSC) models to fit the observed spectral energy distribution (SED) from the radio to the X-ray band of the hot spot. Our esults show that in the case of a much weaker magnetic field strength than the equipartition magnetic field, both the single-zone and multi-zone models can fit the SED, but the multi-zone model significantly improves the fit. The two models predict the hot spot as a GeV-TeV source, which might be marginally detectable with Fermi/LAT and HESS. The inverse Compton scattering of cosmic microwave background (IC/CMB) is also considered, but its contribution to GeV-TeV emission is negligible. Note that under the equipartition condition, the SED can also be fit with the multi-zone model, but the predicted flux at 10^22 Hz is too weak to be detectable. The detection of TeV gamma-rays from this FR II radio galaxy, if confirmed, would establish a new subclass of extragalactic source in this energy regime since most of the AGNs detected to date at TeV energies are high-energy-peaked BL Lac objects.

The time dependence of the changes in the emission spectra of Comet 9P/Tempel 1 after Deep Impact are derived and discussed. This was a unique event because for the first time it gave astronomers the opportunity to follow the time history of the formation and decay of O(1S), OH, CN, C2, C3, NH, and NH2. Least squares fits of a modified Haser model with constraints using known rate constants were fit to the observed data. In the case of OH a simple two-step Haser model provides a reasonable fit to the observations. Fitting the emissions from O(1S), CN, C2, C3, NH, and NH2 requires the addition of a delayed component to a regular two or three step Haser model. From this information a picture of the Deep Impact encounter emerges where there is an initial formation of gas and dust, which is responsible for the prompt emission that occurs right after impact. A secondary source of gas starts later after impact when the initial dust has dissipated enough so that solar radiation can reach the surface of freshly exposed material. The implications of this and other results are discussed in terms of the implications on the structure and composition of the comet's nucleus.

Using dark matter halos traced by galaxy groups selected from the Sloan Digital Sky Survey Data Release 4, we find that about 1/4 of the faint galaxies ($\rmag >-17.05$, hereafter dwarfs) that are the central galaxies in their own halo are not blue and star forming, as expected in standard models of galaxy formation, but are red. In contrast, this fraction is about 1/2 for dwarf satellite galaxies. Many red dwarf galaxies are physically associated with more massive halos. In total, about $\sim 45$% of red dwarf galaxies reside in massive halos as satellites, while another $\sim 25$% have a spatial distribution that is much more concentrated towards their nearest massive haloes than other dwarf galaxies. We use mock catalogs to show that the reddest population of non-satellite dwarf galaxies are distributed within about 3 times the virial radii of their nearest massive halos. We suggest that this population of dwarf galaxies are hosted by low-mass halos that have passed through their massive neighbors, and that the same environmental effects that cause satellite galaxies to become red are also responsible for the red colors of this population of galaxies. We do not find any significant radial dependence of the population of dwarf galaxies with the highest concentrations, suggesting that the mechanisms operating on these galaxies affect color more than structure. However, over 30% of dwarf galaxies are red and isolated and their origin remains unknown.

We report the first results from the GammeV search for chameleon particles, which may be created via photon-photon interactions within a strong magnetic field. Chameleons are hypothesized scalar fields that could explain the dark energy problem. We implement a novel technique to create and trap the reflective particles within a jar and to detect them later via their afterglow as they slowly convert back into photons. These measurements provide the first experimental constraints on the couplings of chameleons to photons.

We use galaxy groups selected from the Sloan Digital Sky Survey (SDSS) together with mass models for individual groups to study the galaxy-galaxy lensing signals expected from galaxies of different luminosities and morphological types. We compare our model predictions with the observational results obtained from the SDSS by Mandelbaum et al. (2006) for the same samples of galaxies. The observational results are well reproduced in a $\Lambda$CDM model based on the WMAP 3-year data, but a $\Lambda$CDM model with higher $\sigma_8$, such as the one based on the WMAP 1-year data,significantly over-predicts the galaxy-galaxy lensing signal. We model, separately, the contributions to the galaxy-galaxy lensing signals from different galaxies: central versus satellite, early-type versus late-type, and galaxies in halos of different masses. We also examine how the predicted galaxy-galaxy lensing signal depends on the shape, density profile, and the location of the central galaxy with respect to its host halo.

Using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey Data Release 4 (SDSS DR4) with an adaptive halo-based group finder, we investigate the luminosity and stellar mass functions for different populations of galaxies (central versus satellite; red versus blue; and galaxies in groups of different masses) and for groups themselves. The conditional stellar mass function (CSMF), which describes the stellar distribution of galaxies in halos of a given mass for central and satellite galaxies can be well modeled with a log-normal distribution and a modified Schechter form, respectively. On average, there are about 3 times as many central galaxies as satellites. Among the satellite population, there are in general more red galaxies than blue ones. For the central population, the luminosity function is dominated by red galaxies at the massive end, and by blue galaxies at the low mass end. At the very low-mass end ($M_\ast \la 10^9 h^{-2}\Msun$), however, there is a marked increase in the number of red centrals. We speculate that these galaxies are located close to large halos so that their star formation is truncated by the large-scale environments. The stellar-mass function of galaxy groups is well described by a double power law, with a characteristic stellar mass at $\sim 4\times 10^{10}h^{-2}\Msun$. Finally, we use the observed stellar mass function of central galaxies to constrain the stellar mass - halo mass relation for low mass halos, and obtain $M_{\ast, c}\propto M_h^{4.9}$ for $M_h \ll 10^{11} \msunh$.

Using a representative sample of 911 central galaxies (CENs) from the SDSS DR4 group catalogue, we study how the structure of the most massive members in groups and clusters depend on (1) galaxy stellar mass (Mstar), (2) dark matter halo mass of the host group (Mhalo), and (3) their halo-centric position. We establish and thoroughly test a GALFIT-based pipeline to fit 2D Sersic models to SDSS data. We find that the fitting results are most sensitive to the background sky level determination and strongly recommend using the SDSS global value. We find that uncertainties in the background translate into a strong covariance between the total magnitude, half-light size (r50), and Sersic index (n), especially for bright/massive galaxies. We find that n depends strongly on Mstar for CENs, but only weakly or not at all on Mhalo. Less (more) massive CENs tend to be disk (spheroid)-like over the full Mhalo range. Likewise, there is a clear r50-Mstar relation for CENs, with separate slopes for disks and spheroids. When comparing CENs with satellite galaxies (SATs), we find that low mass (<10e10.75 Msun/h^2) SATs have larger median n than CENs of similar Mstar. Low mass, late-type SATs have moderately smaller r50 than late-type CENs of the same Mstar. However, we find no size differences between spheroid-like CENs and SATs, and no structural differences between CENs and SATs matched in both mass and colour. The similarity of massive SATs and CENs shows that this distinction has no significant impact on the structure of spheroids. We conclude that Mstar is the most fundamental property determining the basic structure of a galaxy. The lack of a clear n-Mhalo relation rules out a distinct group mass for producing spheroids, and the responsible morphological transformation processes must occur at the centres of groups spanning a wide range of masses. (abridged)

$f(R)$ gravity, capable of driving the late-time acceleration of the universe, is emerging as a promising alternative to dark energy. Various $f(R)$ gravity models have been intensively tested against probes of the expansion history, including type Ia supernovae (SNIa), the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO). In this paper we propose to use the statistical lens sample from Sloan Digital Sky Survey Quasar Lens Search Data Release 3 (SQLS DR3) to constrain $f(R)$ gravity models. This sample can probe the expansion history up to $z\sim2.2$, higher than what probed by current SNIa and BAO data. We adopt a typical parameterization of the form $f(R)=R-\alpha H^2_0(-\frac{R}{H^2_0})^\beta$ with $\alpha$ and $\beta$ constants. For $\beta=0$ ($\Lambda$CDM), we obtain the best-fit value of the parameter $\alpha=-4.193$, for which the 95% confidence interval that is [-4.633, -3.754]. This best-fit value of $\alpha$ corresponds to the matter density parameter $\Omega_{m0}=0.301$, consistent with constraints from other probes. Allowing $\beta$ to be free, the best-fit parameters are $(\alpha, \beta)=(-3.777, 0.06195)$. Consequently, we give $\Omega_{m0}=0.285$ and the deceleration parameter $q_0=-0.544$. At the 95% confidence level, $\alpha$ and $\beta$ are constrained to [-4.67, -2.89] and [-0.078, 0.202] respectively. Clearly, given the currently limited sample size, we can only constrain $\beta$ within the accuracy of $\Delta\beta\sim 0.1$ and thus can not distinguish between $\Lambda$CDM and $f(R)$ gravity with high significance, and actually, the former lies in the 68% confidence contour. We expect that the extension of the SQLS DR3 lens sample to the SDSS DR5 and SDSS-II will make constraints on the model more stringent.

We investigate the correlation of star formation quenching with internal galaxy properties and large scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based Group Catalog of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer (GALEX) are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of "central" and "satellite" galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with AGN feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This {\em satellite overquenching problem} is caused by the too-rapid stripping of the satellites' hot gas halos, which leads to rapid strangulation of star formation.

Using the SDSS DR4 group catalogue of Yang et al. (2007), we investigate sizes, concentrations, colour gradients and surface brightness profiles of central and satellite galaxies. We compare central and satellite galaxies at fixed stellar mass, in order to disentangle environmental from stellar mass dependencies. Early and late type galaxies are defined according to concentration. We find that at fixed stellar mass, late type satellite galaxies have smaller radii and larger concentrations than late type central galaxies. No such differences are found for early-type galaxies. We have also constructed surface brightness and colour profiles for the central and satellite galaxies in our sample. We find that late-type satellite galaxies have a lower surface brightness and redder colours than late-type central galaxies. We show that all observed differences between satellite and central galaxies can be explained by a simple fading model, in which the star formation in the disk decreases over timescales of 2-3 Gyr after a galaxy becomes a satellite. Processes that induce strong morphological changes (e.g. harassment) and processes that strip the galaxy of its entire ISM need not to be invoked in order to explain the environmental dependencies we find.

We develop a new method to reconstruct the cosmic density field from the distribution of dark matter haloes above a certain mass threshold. Our motivation is that well-defined samples of galaxy groups/clusters, which can be used to represent the dark halo population, can now be selected from large redshift surveys of galaxies, and our ultimate goal is to use such data to reconstruct the cosmic density field in the local universe. Our reconstruction method starts with a sample of dark matter haloes above a given mass threshold. Each volume element in space is assigned to the domain of the nearest halo according to a distance measure that is scaled by the virial radius of the halo. The distribution of the mass in and around dark matter haloes of a given mass is modelled using the cross-correlation function between dark matter haloes and the mass distribution within their domains. We use N-body cosmological simulations to show that the density profiles required in our reconstruction scheme can be determined reliably from large cosmological simulations, and that our method can reconstruct the density field accurately using haloes with masses down to $\sim 10^{12}\msun$ (above which samples of galaxy groups can be constructed from current large redshift surveys of galaxies). Working in redshift space, we demonstrate that the redshift distortions due to the peculiar velocities of haloes can be corrected in an iterative way. We also describe some applications of our method.

In the standard paradigm, satellite galaxies are believed to be associated with the population of dark matter subhalos. In this paper, we use the conditional stellar mass functions of {\it satellite galaxies} obtained from a large galaxy group catalogue together with models of the subhalo mass functions to explore the fraction and fate of stripped stars from satellites in galaxy groups and clusters of different masses. The majority of the stripped stars in massive halos are predicted to end up as intra-cluster stars, and the predicted amounts of the intra-cluster component as a function of the velocity dispersion of galaxy system match well the observational results obtained by Gonzalez et al. (2007). The fraction of the mass in the stripped stars to that remain bound in the central and satellite galaxies is the highest ($\sim 40%$ of the total stellar mass) in halos with masses $M_h\sim 10^{14}\msunh$. If all these stars end up in the intra-cluster component (Max), or maximum of them are accreted into the central galaxy (Min), then we can predict that a maximum $\sim 19%$ and a minimum $\sim 5%$ of the total stars in the whole universe are in terms of the diffused intra-cluster component. In the former case, in massive halos with $M_h \sim 10^{15} \msunh$, the stellar mass of the intra-cluster component is roughly 6 times as large as that of the central galaxy. This factor decreases to $\sim 2$, 1 and 0.1 in halos with $M_h \sim 10^{14}$, $10^{13}$, and $10^{12} \msunh$, respectively. The total amount of stars stripped from satellite galaxies is insufficient to build up the central galaxies in halos with masses $\la 10^{12.5}\msunh$, and so the quenching of star formation must occur in halos with higher masses. Abridged.

We performed a Cr-K emission line survey in young supernova remnants (SNRs) with the Chandra archival data. Our sample includes W49B, Cas A, Tycho and Kepler. We confirmed the existence of the Cr line in W49B and discovered this emission line in the other three SNRs. The line center energies, equivalent widths (EWs) and fluxes of the Cr lines are given. The Cr in Cas A is in a high ionization state while that in Tycho and Kepler is in a much lower one. We find a good positive correlation between Cr and Fe line center energies, suggesting a common origin of Cr and Fe in the nucleosynthesis, which is consistent with the theoretical predictions. We propose that the EW ratio between Cr and Fe can be used as a supplementary constraint on the progenitors' properties and the explosion mechanism.

The development of observational techniques has inriched our knowledge of supernova remnants. In this paper, we review the main progresses in the last decade, including new discoveries of supernova remnants and the associated (rare type of) pulsars, nucleosynthesis, the interaction between supernova remnants and molecular clouds, dust in the supernova remnants, shock physics, and cosmic ray accelerations.

We use a modified version of the halo-based group finder developed by Yang et al. to select galaxy groups from the Sloan Digital Sky Survey (SDSS DR4). In the first step, a combination of two methods is used to identify the centers of potential groups and to estimate their characteristic luminosity. Using an iterative approach, the adaptive group finder then uses the average mass-to-light ratios of groups, obtained from the previous iteration, to assign a tentative mass to each group. This mass is then used to estimate the size and velocity dispersion of the underlying halo that hosts the group, which in turn is used to determine group membership in redshift space. Finally, each individual group is assigned two different halo masses: one based on its characteristic luminosity, and the other based on its characteristic stellar mass. Applying the group finder to the SDSS DR4, we obtain 301237 groups in a broad dynamic range, including systems of isolated galaxies. We use detailed mock galaxy catalogues constructed for the SDSS DR4 to test the performance of our group finder in terms of completeness of true members, contamination by interlopers, and accuracy of the assigned masses. This paper is the first in a series and focuses on the selection procedure, tests of the reliability of the group finder, and the basic properties of the group catalogue (e.g. the mass-to-light ratios, the halo mass to stellar mass ratios, etc.). The group catalogues including the membership of the groups are available at http://gax.shao.ac.cn/data/Group.html and http://www.astro.umass.edu/~xhyang/Group.html

The kinematics of satellite galaxies reflect the masses of the extended dark matter haloes in which they orbit, and thus shed light on the mass-luminosity relation (MLR) of their corresponding central galaxies. In this paper we select a large sample of centrals and satellites from the Sloan Digital Sky Survey (SDSS) and measure the kinematics (velocity dispersions) of the satellite galaxies as a function of the $r$-band luminosity of the central galaxies. Using the analytical framework presented in Paper I, we use these data to infer {\it both} the mean and the scatter of the MLR of central galaxies, carefully taking account of selection effects and biases introduced by the stacking procedure. As expected, brighter centrals on average reside in more massive haloes. In addition, we find that the scatter in halo masses for centrals of a given luminosity, $\sigma_{\log M}$, also increases with increasing luminosity. As we demonstrate, this is consistent with $\sigma_{\log L}$, which reflects the scatter in the conditional probability function $P(L_c|M)$, being independent of halo mass. Our analysis of the satellite kinematics yields $\sigma_{\log L}=0.16\pm0.04$, in excellent agreement with constraints from clustering and group catalogues, and with predictions from a semi-analytical model of galaxy formation. We thus conclude that the amount of stochasticity in galaxy formation, which is characterized by $\sigma_{\log L}$, is well constrained, is independent of halo mass, and is in good agreement with current models of galaxy formation.

Galaxy clustering and galaxy-galaxy lensing probe the connection between galaxies and their dark matter haloes in complementary ways. On one hand, the halo occupation statistics inferred from the observed clustering properties of galaxies are degenerate with the adopted cosmology. Consequently, different cosmologies imply different mass-to-light ratios for dark matter haloes. On the other hand, galaxy-galaxy lensing yields direct constraints on the actual mass-to-light ratios and it can be used to break this degeneracy, and thus to constrain cosmological parameters. In this paper we establish the link between galaxy luminosity and dark matter halo mass using the conditional luminosity function (CLF). We constrain the CLF parameters using the galaxy luminosity function and the luminosity dependence of the correlation lengths of galaxies. The resulting CLF models are used to predict the galaxy-galaxy lensing signal. For a cosmology with $(\Omega_{\rm m},\sigma_8)=(0.238,0.734)$, the model accurately fits the galaxy-galaxy lensing data obtained from the SDSS. For a comparison cosmology with $(\Omega_{\rm m},\sigma_8)=(0.3,0.9)$, however, we can accurately fit the luminosity function and clustering properties of the galaxy population, but the model predicts mass-to-light ratios that are too high, resulting in a strong overprediction of the galaxy-galaxy lensing signal. We conclude that the combination of galaxy clustering and galaxy-galaxy lensing is a powerful probe of the galaxy-dark matter connection, with the potential to yield tight constraints on cosmological parameters. Since this method mainly probes the mass distribution on non-linear scales, it is complementary to constraints obtained from the galaxy power-spectrum, which mainly probes the large-scale (linear) matter distribution.

We use a SDSS galaxy group catalogue to study the dependence of galaxy activity on stellar mass, halo mass, and group hierarchy (centrals vs. satellites). We split our galaxy sample in star-forming galaxies, galaxies with optical AGN activity and radio sources. We find a smooth transition in halo mass as the activity of central galaxies changes from star formation to optical AGN activity to radio emission. Star-forming centrals preferentially reside in haloes with M<10^{12} Msun, central galaxies with optical-AGN activity typically inhabit haloes with M \sim 10^{13} Msun, and centrals emitting in the radio mainly reside in haloes more massive than 10^{14} Msun. Although this seems to suggest that the environment (halo mass) determines the type of activity of its central galaxy, we find a similar trend with stellar mass: central star formers typically have stellar masses below 10^{10} Msun, while optical-AGN hosts and central radio sources have characteristic stellar masses of 10^{10.8} Msun and 10^{11.6} Msun, respectively. Since more massive haloes typically host more massive centrals, it is unclear whether the activity of a central galaxy is causally connected to its stellar mass or to its halo mass. In general, satellite galaxies have their activity suppressed wrt central galaxies of the same stellar mass. At fixed stellar mass, we find that the activity of satellite galaxies depends only weakly on halo mass. In fact, for satellite galaxies the dependence of galaxy activity on halo mass is more than four times weaker than the dependence on stellar mass. As we discuss, all these results are consistent with a picture in which low mass haloes accrete cold gas, while massive haloes have coronae of hot gas that promote radio activity of their central galaxies. [Abridged]

In measuring the power spectrum of the distribution of large numbers of dark matter particles in simulations, or galaxies in observations, one has to use Fast Fourier Transforms (FFT) for calculational efficiency. However, because of the required mass assignment onto grid points in this method, the measured power spectrum $\la |\delta^f(k)|^2\ra$ obtained with an FFT is not the true power spectrum $P(k)$ but instead one that is convolved with a window function $|W(\vec k)|^2$ in Fourier space. In a recent paper, Jing (2005) proposed an elegant algorithm to deconvolve the sampling effects of the window function and to extract the true power spectrum, and tests using N-body simulations show that this algorithm works very well for the three most commonly used mass assignment functions, i.e., the Nearest Grid Point (NGP), the Cloud In Cell (CIC) and the Triangular Shaped Cloud (TSC) methods. In this paper, rather than trying to deconvolve the sampling effects of the window function, we propose to select a particular function in performing the mass assignment that can minimize these effects. An ideal window function should fulfill the following criteria: (i) compact top-hat like support in Fourier space to minimize the sampling effects; (ii) compact support in real space to allow a fast and computationally feasible mass assignment onto grids. We find that the scale functions of Daubechies wavelet transformations are good candidates for such a purpose. Our tests using data from the Millennium Simulation show that the true power spectrum of dark matter can be accurately measured at a level better than 2% up to $k=0.7k_N$, without applying any deconvolution processes. The new scheme is especially valuable for measurements of higher order statistics, e.g. the bi-spectrum,........

We use a sample of galaxy groups selected from the SDSS DR 4 with an adaptive halo-based group finder to probe how the clustering strength of groups depends on their masses and colors. In particular, we determine the relative biases of groups of different masses, as well as that of groups with the same mass but with different colors. In agreement with previous studies, we find that more massive groups are more strongly clustered, and the inferred mass dependence of the halo bias is in good agreement with predictions for the $\Lambda$CDM cosmology. Regarding the color dependence, we find that groups with red centrals are more strongly clustered than groups of the same mass but with blue centrals. Similar results are obtained when the color of a group is defined to be the total color of its member galaxies. The color dependence is more prominent in less massive groups and becomes insignificant in groups with masses $\gta 10^{14}\msunh$. We construct a mock galaxy redshift survey constructed from the large Millenium simulation that is populated with galaxies according to the semi-analytical model of Croton et al. Applying our group finder to this mock survey, and analyzing the mock data in exactly the same way as the true data, we are able to accurately recover the intrinsic mass and color dependencies of the halo bias in the model. This suggests that our group finding algorithm and our method of assigning group masses do not induce spurious mass and/or color dependencies in the group-galaxy correlation function. The semi-analytical model reveals the same color dependence of the halo bias as we find in our group catalogue. In halos with $M\sim 10^{12}\msunh$, though, the strength of the color dependence is much stronger in the model than in the data.

Using the Sloan Digital Sky Survey (SDSS) galaxy group catalogue of Yang et al. (2007), we study the average colour and average concentration of satellite galaxies as function of (i) their stellar mass, (ii) their group mass, and (iii) their group-centric radius. We find that the colours and concentrations of satellite galaxies are (almost) completely determined by their stellar mass. In particular, at fixed stellar mass, the average colours and concentrations of satellite galaxies are independent of either halo mass or halo-centric radius. We find clear evidence for mass segregation of satellite galaxies in haloes of all masses, and argue that this explains why satellites at smaller halo-centric radii are somewhat redder and somewhat more concentrated. In addition, the weak colour and concentration dependence of satellite galaxies on halo mass is simply a reflection of the fact that more massive haloes host, on average, more massive satellites. Combining these results with the fact that satellite galaxies are, on average, redder and somewhat more concentrated than central galaxies of the same stellar mass, the following picture emerges: galaxies become redder and somewhat more concentrated once they fall into a bigger halo (i.e., once they become a satellite galaxy). This is a clear manifestation of environment dependence. However, there is no indication that the magnitude of the transformation (or its timescale) depends on environment; a galaxy undergoes a transition when it becomes a satellite, but it does not matter whether it becomes a satellite of a small (Milky Way sized) halo, or of a massive cluster. We discuss the implication of this `dearth' of environment dependence for the physical processes responsible for transforming satellite 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.

We report the first results of the GammeV experiment, a search for milli-eV mass particles with axion-like couplings to two photons. The search is performed using a "light shining through a wall" technique where incident photons oscillate into new weakly interacting particles that are able to pass through the wall and subsequently regenerate back into detectable photons. The oscillation baseline of the apparatus is variable, thus allowing probes of different values of particle mass. We find no excess of events above background and are able to constrain the two-photon couplings of possible new scalar (pseudoscalar) particles to be less than 3.1x10^{-7} GeV^{-1} (3.5x10^{-7} GeV^{-1}) in the limit of massless particles.

We study the three-dimensional and projected shapes of galaxy groups in the Sloan Digital Sky Survey Data Release 4, and examine the alignment between the orientation of the central galaxy and the spatial distribution of satellite galaxies. The projected ellipticity of a group is measured using the moments of the discrete distribution of its member galaxies. We infer the three-dimensional and projected axis ratios of their dark matter haloes by comparing the measured ellipticity distributions with those obtained from Monte Carlo simulations of projected, triaxial dark matter haloes with different axis ratios. We find that the halo shape has a strong dependence on the halo mass. While the haloes of low-mass groups are nearly spherical, those of massive groups tend to be prolate. For groups containing at least four members, the statistical distribution of their measured ellipticities does not have a strong dependence on the colors of their central galaxies. Our analysis further shows that the average three-dimensional axis ratio for haloes with $12<{\rm log}[M/(h^{-1}M_{\odot})]\leq15$ is about $1:0.46:0.46$, resulting in a projected axis ratio of $\sim 0.77$. Our results for the alignment between the orientation of the central galaxy of a group and the distribution of their satellite galaxies are in broad agreement with those obtained by Yang et al. The distribution of satellite galaxies preferentially aligns with the major axis of the central galaxy, with a clear dependence on both halo mass and galaxy colors. (abridged)

We investigate various halo occupation statistics using a large galaxy group catalogue constructed from the SDSS DR4 with an adaptive halo-based group finder. The conditional luminosity function (CLF) is measured separately for all, red and blue galaxies, as well as in terms of central and satellite galaxies. The CLFs for central and satellite galaxies can be well modelled with a log-normal distribution and a modified Schechter form, respectively. About 85% of the central galaxies and about 80% of the satellite galaxies in halos with masses $M_h\ga 10^{14}\msunh$ are red galaxies. These numbers decrease to 50% and 40%, respectively, in halos with $M_h \sim 10^{12}\msunh$. For halos of a given mass, the distribution of the luminosities of central galaxies, $L_c$, has a dispersion of about 0.15 dex. The mean luminosity (stellar mass) of the central galaxies scales with halo mass as $L_c\propto M_h^{0.17}$ ($M_{*,c}\propto M_h^{0.22}$) for halos with masses $M\gg 10^{12.5}\msunh$, and both relations are significantly steeper for less massive halos. We also measure the luminosity (stellar mass) gap between the first and second brightest (most massive) member galaxies, $\log L_1 - \log L_2$ ($\log M_{*,1}-\log M_{*,2}$). These gap statistics, especially in halos with $M_h \la 10^{14.0}\msunh$, indicate that the luminosities of central galaxies are clearly distinct from those of their satellites. The fraction of fossil groups, defined as those groups with $\log L_1 - \log L_2\ge 0.8$, ranges from $\sim 2.5%$ for groups with $M_h\sim 10^{14}\msunh$ to 18-60% for groups with $M_h\sim 10^{13}\msunh$. Finally, we measure the fraction of satellites, which changes from $\sim 5.0%$ for galaxies with $\rmag\sim -22.0$ to $\sim40%$ for galaxies with $\rmag\sim -17.0$. (abridged)

We present high spatial resolution X-ray spectroscopy of supernova remnant Cassiopeia A with the {\sl Chandra} observations. The X-ray emitting region of this remnant was divided into 38 $\times$ 34 pixels with a scale of 10$\arcsec$ $\times$ 10$\arcsec$ each. Spectra of 960 pixels were created and fitted with an absorbed two component non-equilibrium ionization model. With the spectral analysis results we obtained maps of absorbing column density, temperatures, ionization ages, and the abundances for Ne, Mg, Si, S, Ca and Fe. The Si, S and possibly Ca abundance maps show obviously jet structures, while Fe doesn't follow the jet but seems to be distributed perpendicular to it. In the range of about two orders of magnitude, the abundances of Si, S and Ca show tight correlations between each other, suggesting them to be ejecta from explosive O-burning and incomplete Si-burning. Meanwhile, Ne abundance is well correlated with that of Mg, indicating them to be the ashes of explosive C/Ne burning. The Fe abundance is positively correlated with that of Si when Si abundance is lower than 3 solar abundances, but a negative correlation appears when the Si abundance is higher. We suggest that such a two phase correlation is the results of different ways in which Fe is synthesized.

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.

In the current paradigm, red sequence galaxies are believed to have formed as blue disk galaxies that subsequently had their star formation quenched. Since red-sequence galaxies typically have an early-type morphology, the transition from the blue to the red sequence also involves a morphological transformation. In this paper we study the impact of transformation mechanisms that operate only on satellite galaxies, such as strangulation, ram-pressure stripping and galaxy harassment. Using a large galaxy group catalogue constructed from the SDSS, we compare the colors and concentrations of satellites galaxies to those of central galaxies of the same stellar mass, adopting the hypothesis that the latter are the progenitors of the former. On average, satellites are redder and more concentrated than central galaxies of the same stellar mass. Central-satellite pairs that are matched in both stellar mass and color, however, show no average concentration difference, indicating that the transformation mechanisms affect color more than morphology. The color and concentration differences of matched central-satellite pairs are completely independent of the halo mass of the satellite galaxy, indicating that satellite-specific transformation mechanisms are equally efficient in haloes of all masses. This strongly favors strangulation as the main quenching mechanism for satellite galaxies. Finally, we determine the relative importance of satellite quenching for the build-up of the red sequence. We find that roughly 70 percent of red sequence satellite galaxies with a stellar mass of 10^9 Msun had their star formation quenched as satellites. This drops rapidly to zero with increasing stellar mass, indicating that a significant fraction of red satellites were already quenched before they became a satellite.

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.

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.

We investigate the incidence of major mergers creating >10e11 Msun galaxies in present-day groups and clusters more massive than 2.5e13 Msun. We identify 38 pairs of massive galaxies with mutual tidal interaction signatures selected from >5000 galaxies with >5e10 Msun that reside in 845 such groups. We fit the images of each galaxy pair as the line-of-sight projection of symmetric models and identify mergers by the presence of residual asymmetries around each progenitor, such as off-center isophotes, broad tidal tails, and dynamical friction wakes. At the resolution and sensitivity of the SDSS, such mergers are found in 16% of high-mass, galaxy-galaxy pairs with magnitude differences of <1.5 and <30 kpc projected separations. We find that 90% of these mergers have nearly equal-mass progenitors with red-sequence colors and centrally-concentrated morphologies, the hallmarks of dissipationless merger simulations. Mergers at group centers are more common than between 2 satellites, but both are morphologically indistinguishable and we tentatively conclude that the latter are likely located at the dynamical centers of recently accreted subhalos. The frequency of central and satellite merging diminishes with group mass consistent with dynamical friction expectations. Based on reasonable assumptions, the centers of these massive halos are growing in stellar mass by 1-9% per Gyr, on average. Compared to all LRG-LRG mergers, we find a 2-9 times higher rate for their merging when restricted to these dense environments. Our results imply that the massive end of the galaxy population continues to evolve hierarchically at a measurable level, and that the centers of massive groups are the preferred environment for merger-driven galaxy assembly. (abridged)

Using a large galaxy group catalogue based on the Sloan Digital Sky Survey Data Release 4 we measure three different types of intrinsic galaxy alignment within groups: halo alignment between the orientation of the brightest group galaxies (BGG) and the distribution of its satellite galaxies, radial alignment between the orientation of a satellite galaxy and the direction towards its BGG, and direct alignment between the orientation of the BGG and that of its satellites. In agreement with previous studies we find that satellite galaxies are preferentially located along the major axis. In addition, on scales r < 0.7 Rvir we find that red satellites are preferentially aligned radially with the direction to the BGG. The orientations of blue satellites, however, are perfectly consistent with being isotropic. Finally, on scales r < 0.1 \Rvir, we find a weak but significant indication for direct alignment between satellites and BGGs. We briefly discuss the implications for weak lensing measurements.

We propose a new concept for a very low background multi-ton liquid xenon Dark Matter experiment. The detector consists of two concentric spheres and a charge readout device in the centre. Xenon between the two spheres forms a self-shield and veto device. The inner surface of the central sphere is coated with CsI to form an internal photocathode with minimum of 2\pi coverage for any event in the active volume. Photoelectrons from the CsI photocathode drift toward the charge readout micro-structure in the centre of the detector. Both scintillation and ionisation are measured simultaneously for background rejection and 3-D event mapping. In addition to external shielding, the low background is achieved by eliminating PMTs and by using low radioactivity pure materials throughout the detector. We present detailed calculations of the charge readout system and design details. The detector is expected to probe the full SUSY parameter space.

We study the cross-correlation between galaxies of different luminosities and colors, using a sample selected from the SDSS Dr 4. Galaxies are divided into 6 samples according to luminosity, and each of these samples is divided into red and blue subsamples. Projected auto-correlation and cross-correlation is estimated for these subsample. At projected separations r_p > 1\mpch, all correlation functions are roughly parallel, although the correlation amplitude depends systematically on luminosity and color. On r_p < 1\mpch, the auto- and cross-correlation functions of red galaxies are significantly enhanced relative to the corresponding power laws obtained on larger scales. Such enhancement is absent for blue galaxies and in the cross-correlation between red and blue galaxies. We esimate the relative bias factor on scales r > 1\mpch for each subsample using its auto-correlation function and cross-correlation functions. The relative bias factors obtained from different methods are similar. For blue galaxies the luminosity-dependence of the relative bias is strong over the luminosity range probed (-23.0<M_r < -18.0),but for red galaxies the dependence is weaker and becomes insignificant for luminosities below L^*. To examine whether a significant stochastic/nonlinear component exists in the bias relation, we study the ratio R_ij= W_{ii}W_{jj}/W_{ij}^2, where W_{ij} is the projected correlation between subsample i and j. We find that the values of R_ij are all consistent with 1 for all-all, red-red and blue-blue samples, however significantly larger than 1 for red-blue samples. For faint red - faint blue samples the values of R_{ij} are as high as ~ 2 on small scales r_p < 1 \mpch and decrease with increasing r_p.

Recent studies have shown that the distribution of satellite galaxies is preferentially aligned with the major axis of their central galaxy. The strength of this alignment has been found to depend strongly on the colours of the satellite and central galaxies, and only weakly on the mass of the halo in which the galaxies reside. In this paper we study whether these alignment signals, and their dependence on galaxy and halo properties, can be reproduced in a hierarchical structure formation model of a $\Lambda$CDM concordance cosmology. To that extent we use a large $N$-body simulation which we populate with galaxies following a semi-analytical model for galaxy formation. We find that if the orientation of the central galaxy is perfectly aligned with that of its dark matter halo, then the predicted central-satellite alignment signal is much stronger than observed. If, however, the minor axis of a central galaxy is perfectly aligned with the angular momentum vector of its dark matter halo, we can accurately reproduce the observed alignment strength as function of halo mass and galaxy color. Although this suggests that the orientation of central galaxies is governed by the angular momentum of their dark matter haloes, we emphasize that any other scenario in which the minor axes of central galaxy and halo are misaligned by $\sim 40^{\circ}$ (on average) will match the data equally well. Finally, we show that dependence of the alignment strength on the color of the central galaxy is most likely an artefact due to interlopers in the group catalogue. The dependence on the color of the satellite galaxies, on the other hand, is real and owes to the fact that red satellites are associated with subhaloes that were more massive at their time of accretion.

We use the conditional luminosity function (CLF) and data from the 2dFGRS to constrain the average relation between light and mass in a LCDM cosmology with Omega_m=0.23 and sigma_8=0.74 (hereafter WMAP3 cosmology). Reproducing the observed luminosity dependence of the galaxy two-point correlation function results in average mass-to-light ratios that are about 35 percent lower than in a LCDM cosmology with Omega_m=0.3 and sigma_8=0.9 (hereafter WMAP1 cosmology). This removes an important problem with previous halo occupation models which had a tendency to predict cluster mass-to-light ratios that were too high. For the WMAP3 cosmology our model yields average mass-to-light ratios, central galaxy luminosities, halo occupation numbers, satellite fractions, and luminosity-gap statistics, that are all in excellent agreement with those obtained from a 2dFGRS group catalogue and from other independent studies. We also use our CLF model to compute the probability distribution P(M|L_cen), that a central galaxy of luminosity L_cen resides in a halo of mass M. We find this distribution to be much broader than what is typically assumed in HOD models, which has important implications for the interpretation of galaxy-galaxy lensing data. Finally, reproducing the luminosity dependence of the pairwise velocity dispersions in the 2dFGRS requires relatively low mass-to-light ratios for clusters and a satellite fraction that decreases strongly with increasing luminosity. This is only marginally consistent with our CLF constraints. We argue that a cosmology with parameters between those of the WMAP1 and WMAP3 cosmologies is likely to yield results with a higher level of consistency.

Galaxy-galaxy lensing is rapidly becoming one of the most promising means to accurately measure the average relation between galaxy properties and halo mass. In order to obtain a signal of sufficient signal-to-noise, one needs to stack many lens galaxies according to their property of interest, such as luminosity or stellar mass. Since such a stack consists of both central and satellite galaxies, which contribute very different lensing signals, the resulting shear measurements can be difficult to interpret. In the past, galaxy-galaxy lensing studies have either completely ignored this problem, have applied rough isolation criteria in an attempt to preferentially select `central' galaxies, or have tried to model the contribution of satellites explicitely. However, if one is able to {\it a priori} split the galaxy population in central and satellite galaxies, one can measure their lensing signals separately. This not only allows a much cleaner measurement of the relation between halo mass and their galaxy populations, but also allows a direct measurement of the sub-halo masses around satellite galaxies. In this paper, we use a realistic mock galaxy redshift survey to show that galaxy groups, properly selected from large galaxy surveys, can be used to accurately split the galaxy population in centrals and satellites. Stacking the resulting centrals according to their group mass, estimated from the total group luminosity, allows a remarkably accurate recovery of the masses and density profiles of their host haloes. In addition, stacking the corresponding satellite galaxies according to their projected distance from the group center yields a lensing signal that can be used to accurate measure the masses of both sub-haloes and host haloes. (Abridged)

We examine the relation between galaxy properties and environment in the SDSS DR2, quantifying environment in terms of the mass of the host halo, which is obtained with a new iterative group finder. We find that galaxy type fractions scale strongly and smoothly with halo mass, but, at fixed mass, not with luminosity. We compare these findings with the semi-analytical galaxy formation model of Croton et al. (2006). The discrepancies we find can be explained with an oversimplified implementation of strangulation, the neglect of tidal stripping, and shortcomings in the treatments of dust extinction and/or AGN feedback.

Successfully reproducing the galaxy luminosity function and the bimodality in the galaxy distribution requires a mechanism that can truncate star formation in massive haloes. Current models of galaxy formation consider two such truncation mechanisms: strangulation, which acts on satellite galaxies, and AGN feedback, which predominantly affects central galaxies. The efficiencies of these processes set the blue fraction of galaxies as function of galaxy luminosity and halo mass. In this paper we use a galaxy group catalogue extracted from the Sloan Digital Sky Survey (SDSS) to determine these fractions. To demonstrate the potential power of this data as a benchmark for galaxy formation models, we compare the results to the semi-analytical model for galaxy formation of Croton et al. (2006). Although this model accurately fits the global statistics of the galaxy population, as well as the shape of the conditional luminosity function, there are significant discrepancies when the blue fraction of galaxies as a function of mass and luminosity is compared between the observations and the model. In particular, the model predicts (i) too many faint satellite galaxies in massive haloes, (ii) a blue fraction of satellites that is much too low, and (iii) a blue fraction of centrals that is too high and with an inverted luminosity dependence. In the same order, we argue that these discrepancies owe to (i) the neglect of tidal stripping in the semi-analytical model, (ii) the oversimplified treatment of strangulation, and (iii) improper modeling of dust extinction and/or AGN feedback. The data presented here will prove useful to test and calibrate future models of galaxy formation and in particular to discriminate between various models for AGN feedback and other star formation truncation mechanisms.

We use galaxy groups selected from the Sloan Digital Sky Survey to examine the alignment between the orientation of the central galaxy (defined as the brightest group member) and the distribution of satellite galaxies. By construction, we therefore only address the alignment on scales smaller than the halo virial radius. We find a highly significant alignment of satellites with the major axis of their central galaxy. This is in qualitative agreement with the recent study of Brainerd (2005), but inconsistent with several previous studies who detected a preferential minor axis alignment. The alignment strength in our sample is strongest between red central galaxies and red satellites. On the contrary, the satellite distribution in systems with a blue central galaxy is consistent with isotropic. We also find that the alignment strength is stronger in more massive haloes and at smaller projected radii from the central galaxy. In addition, there is a weak indication that fainter (relative to the central galaxy) satellites are more strongly aligned. We present a detailed comparison with previous studies, and discuss the implications of our findings for galaxy formation.

We study the dependence of the cross-correlation between galaxies and galaxy groups on group properties. Confirming previous results, we find that the correlation strength is stronger for more massive groups, in good agreement with the expected mass dependence of halo bias. We also find, however, that for groups of the same mass, the correlation strength depends on the star formation rate (SFR) of the central galaxy: at fixed mass, the bias of galaxy groups decreases as the SFR of the central galaxy increases. We discuss these findings in light of the recent findings by Gao et al (2005) that halo bias depends on halo formation time, in that halos that assemble earlier are more strongly biased. We also discuss the implication for galaxy formation, and address a possible link to galaxy conformity, the observed correlation between the properties of satellite galaxies and those of their central galaxy.

Using a large galaxy group catalogue constructed from the SDSS, we investigate the correlation between various galaxy properties and halo mass. We split the population of galaxies in early types, late types, and intermediate types, based on their colour and specific star formation rate. At fixed luminosity, the early type fraction increases with increasing halo mass. Most importantly, this mass dependence is smooth and persists over the entire mass range probed, without any break or feature at any mass scale. We argue that the previous claim of a characteristic feature on galaxy group scales is an artefact of the environment estimators used. At fixed halo mass, the luminosity dependence of the type fractions is surprisingly weak: galaxy type depends more strongly on halo mass than on luminosity. We also find that the early type fraction decreases with increasing halo-centric radius. Contrary to previous studies, we find that this radial dependence is also present in low mass haloes. The properties of satellite galaxies are strongly correlated with those of their central galaxy. In particular, the early type fraction of satellites is significantly higher in a halo with an early type central galaxy than in a halo of the same mass but with a late type central galaxy. This phenomenon, which we call `galactic conformity', is present in haloes of all masses and for satellites of all luminosities. Finally, the fraction of intermediate type galaxies is always ~20 percent, independent of luminosity, independent of halo mass, independent of halo-centric radius, and independent of whether the galaxy is a central galaxy or a satellite galaxy. We discuss the implications of all these findings for galaxy formation and evolution.

We use recent observations of the HI-mass function to constrain galaxy formation. The data conflicts with the standard model where most of the gas in a low-mass dark matter halo is assumed to settle into a disk of cold gas that is depleted by star formation and supernova-driven outflows until the disk becomes gravitationally stable. A consistent model can be found if low-mass haloes are embedded in a preheated medium, with a specific gas entropy ~ 10Kev cm^2. Such a model simultaneously matches the faint-end slope of the galaxy luminosity function. We propose a preheating model where the medium around low-mass haloes is preheated by gravitational pancaking. Since gravitational tidal fields suppress the formation of low-mass haloes while promoting that of pancakes, the formation of massive pancakes precedes that of the low-mass haloes within them. We demonstrate that the progenitors of present-day dark matter haloes with M<10^{12}h^{-1}\msun were embedded in pancakes of masses $~5x10^{12}h^{-1}\msun$ at z~2. The formation of such pancakes heats the gas to a temperature of 5x10^5K and compresses it to an overdensity of ~10. Such gas has a cooling time that exceeds the age of the Universe at z~2, and has a specific entropy of ~15Kev cm^2, almost exactly the amount required to explain the stellar and HI mass functions. (Abridged)

We determine the cross-correlation function between galaxies and galaxy groups, using both the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey (SDSS). We study the cross-correlation as a function of group mass, and as a function of the luminosity, stellar mass, colour, spectral type and specific star formation rate of the galaxies. All these cross-correlation functions show a clear transition from the `1-halo' to the `2-halo' regimes on a scale comparable to the virial radius of the groups in consideration. On scales larger than the virial radius, all cross-correlation functions are roughly parallel, consistent with the linear bias model. In particular, the large scale correlation amplitudes are higher for more massive groups, and for brighter and redder galaxies. In the `1-halo' regime, the cross-correlation function depends strongly on the definition of the group center. We consider both a luminosity-weighted center (LWC) and a center defined by the location of the brightest group galaxy (BGC). With the first definition, the bright early-type galaxies in massive groups are found to be more centrally concentrated than the fainter, late-type galaxies. Using the BGC, and excluding the brightest galaxy from the cross correlation analysis, we only find significant segregation in massive groups ($M \gta 10^{13}h^{-1}\msun$) for galaxies of different spectral types (or colours or specific star formation rates). In haloes with masses $\la 10^{13}h^{-1}\msun$, there is a significant deficit of bright satellite galaxies. Comparing the results from the 2dFGRS with those obtained from realistic mock samples, we find that the distribution of galaxies in groups is much less concentrated than dark matter haloes predicted by the current $\Lambda$CDM model. (Abridged)

The brightest galaxy in a dark matter halo is expected to reside at rest at the center of the halo. In this paper we test this `Central Galaxy Paradigm' using group catalogues extracted from the 2dFGRS and the SDSS. For each group we compute a parameter, R, defined as the difference between the velocity of the brightest group galaxy and the average velocity of the other group members (hereafter satellites), normalized by the unbiased estimator of the velocity dispersion of the satellite galaxies. Using detailed mock galaxy redshift surveys, we show that the central galaxy paradigm is inconsistent with the data at high confidence, and that instead the brightest halo galaxies have a specific kinetic energy that is about 25 percent of that of the satellites. This indicates that either central galaxies reside at the minimum of the dark matter potential, but that the halo itself is not yet fully relaxed, or, that the halo is relaxed, but that the central galaxy oscillates in its potential well. The former is consistent with the fact that we find a weak hint that the velocity bias of brightest halo galaxies is larger in more massive haloes, while the latter may be indicative of cored, rather than cusped, dark matter haloes. We discuss several implications of these findings, including mass estimates based on satellite kinematics, strong gravitational lensing, halo occupation models, and the frequency and longevity of lopsidedness in disk galaxies.

We study the occupation statistics of galaxies in dark matter haloes using galaxy groups identified from the 2-degree Field Galaxy Redshift Survey with the halo-based group finder of Yang et al. The occupation distribution is considered separately for early and late type galaxies, as well as in terms of central and satellite galaxies. The mean luminosity of the central galaxies scales with halo mass approximately as $L_c\propto M^{2/3}$ for haloes with masses $M<10^{13}h^{-1}\msun$, and as $L_c\propto M^{1/4}$ for more massive haloes. The characteristic mass of $10^{13} h^{-1} \Msun$ is consistent with the mass scale where galaxy formation models suggest a transition from efficient to inefficient cooling. Another characteristic halo mass scale, $M\sim 10^{11} h^{-1}\msun$, which cannot be probed directly by our groups, is inferred from the conditional luminosity function (CLF) that matches the observed galaxy luminosity function and clustering. For a halo of given mass, the distribution of $L_c$ is rather narrow. The satellite galaxies are found to follow a Poissonian number distribution. The central galaxies in low-mass haloes are mostly late type galaxies, while those in massive haloes are almost all early types. We also measure the CLF of galaxies in haloes of given mass. Over the mass range that can be reliably probed with the present data ($13.3 \lta {\rm log}[M/(h^{-1}\Msun)] \lta 14.7$), the CLF is reasonably well fit by a Schechter function. Contrary to recent claims based on semi-analytical models of galaxy formation, the presence of central galaxies does not show up as a strong peak at the bright end of the CLF. (Abridged)

What galaxy lives in what halo? The answer to this simple question holds important information regarding galaxy formation and evolution. We describe a new statistical technique to link galaxies to their dark matter haloes, or light to mass, using the clustering properties of galaxies as function of their luminosity. The galaxy-dark matter connection thus established, and parameterized through the conditional luminosity function, indicates the presence of two characteristic scales in galaxy formation: one at ~10^{11} Msun, where galaxy formation is most efficient, and another at ~10^{13} Msun, where a transition occurs from systems dominated by one brightest, central galaxy to systems with several dominant galaxies of comparable luminosity. We show that the relation between light and mass established from galaxy clustering alone is in excellent agreement with the dynamical masses inferred from satellite kinematics. We also present a new (halo-based) galaxy-group finder, and show that the multiplicity function of galaxy groups identified in the 2dFGRS suggests a relatively high mass-to-light ratio on the scales of galaxy clusters, or, alternatively, a relatively low value of sigma_8. These findings are also supported by our studies of pairwise peculiar velocities and satellite abundances.

We analyze the 2-point correlation function (2PCF) of galaxy groups identified from the 2-degree Field Galaxy Redshift Survey with the halo based group finder recently developed by Yang et al (2004b). With this group catalogue we are able to estimate the 2PCFs for systems ranging from isolated galaxies to rich clusters of galaxies. The real-space correlation length obtained for these systems ranges from $\sim 4\mpch$ to $\sim 15\mpch$, respectively. The observed correlation amplitude (and the corresponding bias factor) as a function of group abundance is well reproduced by associating galaxy groups with dark matter haloes in the standard $\Lambda$CDM model. Redshift distortions are clearly detected in the redshift-space correlation function, the degree of which is consistent with the assumption of gravitational clustering and halo bias in the cosmic density field. In agreement with previous studies we find a strong increase of the correlation length with the mean inter-group separation. Although well determined observationally, we show that current theoretical predictions are not yet accurate enough to allow for stringent constraints on cosmological parameters. Finally, we use our results to explore the power-law nature of the 2PCF of galaxies. We split the 2PCF in 1-group and 2-group terms, equivalent to the 1-halo and 2-halo terms in halo occupation models, and show that the power-law form of the 2PCF is broken, when only including galaxies in the more massive systems.

We use the halo occupation model to calibrate galaxy group finders in magnitude-limited redshift surveys. Since, according to the current scenario of structure formation, galaxy groups are associated with cold dark matter halos, we make use of the properties of the halo population in the design of our group finder. The method starts with an assumed mass-to-light ratio to assign a tentative mass to each group. This mass is used to estimate the size and velocity dispersion of the underlying halo that hosts the group, which in turn is used to determine group membership (in redshift space). This procedure is repeated until no further changes occur in group memberships. We find that the final groups selected this way are insensitive to the mass-to-light ratio assumed. We use mock catalogues, constructed using the conditional luminosity function (CLF), to test the performance of our group finder in terms of completeness of true members and contamination by interlopers. Our group finder is more successful than the conventional Friends-of-Friends group finder in assigning galaxies in common dark matter halos to a single group. We apply our group finder to the 2-degree Field Galaxy Redshift Survey and compare the resulting group properties with model predictions based on the CLF. For the $\Lambda$CDM `concordance' cosmology we find a clear discrepancy between the model and data in the sense that the model predicts too many rich groups. In order to match the observational results we have to either increase the mass-to-light ratios of rich clusters to a level significantly higher than current observational estimates, or to assume $\sigma_8 \simeq 0.7$, compared to the `concordance' value of 0.9.

Using the data acquired in the Time To Spill (TTS) mode for long gamma-ray bursts (GRBs) collected by the Burst and Transient Source Experiment on board the Compton Gamma Ray Observatory (BATSE/CGRO), we have carefully measured spectral lags in time between the low (25-55 keV) and high (110-320 keV) energy bands of individual pulses contained in 64 multi-peak GRBs. We find that the temporal lead by higher-energy gamma-ray photons (i.e., positive lags) is the norm in this selected sample set of long GRBs. While relatively few in number, some pulses of several long GRBs do show negative lags. This distribution of spectral lags in long GRBs is in contrast to that in short GRBs. This apparent difference poses challenges and constraints on the physical mechanism(s) of producing long and short GRBs. The relation between the pulse peak count rates and the spectral lags is also examined. Observationally, there seems to be no clear evidence for systematic spectral lag-luminosity connection for pulses within a given long GRB.

Using detailed mock galaxy redshift surveys (MGRSs) we investigate the abundance and radial distribution of satellite galaxies. The mock surveys are constructed using the conditional luminosity function (CLF), and are compared against data from the 2dFGRS. We use Monte Carlo Markov Chains to explore the full posterior distribution of the CLF parameter space, and show that the average relation between light and mass is tightly constrained. The radial number density distribution of satellite galaxies in the 2dFGRS reveals a pronounced absence of satellites at small projected separations from their host galaxies. This is (at least partially) due to the overlap and merging of galaxy images in the APM catalogue. Because of the resulting close-pair incompleteness we are unfortunately unable to put meaningful constraints on the radial distribution of satellite galaxies. Marginalizing over the full CLF parameter space, we show that in a LCDM concordance cosmology the observed abundances of host and satellite galaxies in the 2dFGRS indicate a power spectrum normalization of sigma_8=0.7. The same cosmology but with sigma_8=0.9 is unable to simultanously match the abundances of host and satellite galaxies. This confirms our previous conclusions based on the pairwise peculiar velocity dispersions and the group multipliticy function.

We present models for the three-point correlation function (3PCF) of both dark matter and galaxies. We show that models based on the halo model can reasonably match the dark matter 3PCF obtained from high-resolution $N$-body simulations. On small scales ($r\la 0.5\mpch$) the 3PCF is sensitive to details regarding the density distributions of dark matter haloes. On larger scales ($r \ga 2.0\mpch$) the results are very sensitive to the abundance of the few most prominant haloes. Using the conditional luminosity function, we also construct models for the 3PCF of galaxies, which we test against large mock galaxy samples. The bias of the galaxy distrubution with respect to the dark matter, and the finite number of galaxies that can be hosted by individual haloes, significantly reduce the normalized three-point correlation function with respect to that of dark matter. Contrary to the 3PCF of the dark matter, the galaxy 3PCF is much less sensitive to details regarding the spatial number density distribution of galaxies in individual haloes or to the abundance of the few most massive systems. Finally, we show that our model based on the conditional luminosity function is in good agreement with results obtained from the 2-degree Field Galaxy Redshift Survey. In particular, the model nicely reproduces the observational finding that the 3PCF for early-type galaxies is slightly higher than that of late-type galaxies, and that there is no significant dependence of the 3PCF on galaxy luminosity.

Using detailed mock galaxy redshift surveys we investigate to what extent the kinematics of large samples of satellites galaxies extracted from flux-limited surveys can be used to constrain halo masses. Previous host-satellite selection criteria yield relatively large fractions of interlopers and with a velocity distribution that, contrary to what has typically been assumed, differs strongly from uniform. A new adaptive selection criterion is proposed which yields much larger host-satellite samples and with strongly reduced interloper fractions. We show that satellite weighting, which occurs naturally when stacking many host-satellite pairs, introduces a bias towards higher velocity dispersions compared to the true, host-averaged mean. A further bias, in the same direction, is introduced when using flux-limited, rather than volume-limited surveys. We apply our adaptive selection criterion to the 2dFGRS and obtain a sample of 12613 satellite galaxies and 8004 host galaxies. The satellite kinematics are in excellent agreement with predictions based on the conditional luminosity function (CLF), providing independent, dynamical confirmation of the average mass-to-light ratios predicted by the CLF formalism.

A basic assumption in current halo occupation model is that the properties of a galaxy depend only on the mass of its dark matter halo. An important consequence of this is that the segregation of the galaxy population by large-scale environment is entirely due to the environmental dependence of the halo population. In this paper we use such a model to predict how the galaxy luminosity function depends on large-scale environment. The latter is represented by the density contrast (delta) averaged over a spherical volume of radius R=8Mpc/h. The model predicts that the Schechter function is a good approximation to the luminosity functions of galaxies brighter than ~10^9 h^{-2}L_sun (b_j-band) in virtually all environments. The characteristic luminosity, L^{\star}, increases moderately with \delta. The faint-end slope, \alpha, on the other hand, is quite independent of $\delta$. However, when splitting the galaxy population into early and late types, it is found that for late-types \alpha is virtually constant, whereas for early-types \alpha increases from ~-0.3 in underdense regions (\delta~-0.5) to ~-0.8 in highly overdense regions with \delta ~ 10. The luminosity function at L_{b_j}<10^9 h^{-2} L_sun is significantly steeper than the extrapolation of the Schechter function that fits the brighter galaxies. This steepening is more significant for early-types and in low-density environments. The model also predicts that the luminosity density and mass density are closely correlated. The relation between the two is monotonic but highly non-linear. This suggests that one can use the luminosity density, averaged over a large volume, to rank the mass density. This, in turn, allows the environmental effects predicted here to be tested by observations.

In two recent papers, we developed a powerful technique to link the distribution of galaxies to that of dark matter haloes by considering halo occupation numbers as function of galaxy luminosity and type. In this paper we use these distribution functions to populate dark matter haloes in high-resolution N-body simulations of the standard LCDM cosmogony with Omega_m=0.3 and sigma_8=0.9. Stacking simulation boxes of 100 Mpc/h and 300 Mpc/h with 512^3 CDM particles each we construct Mock Galaxy Redshift Surveys out to a redshift of z=0.2. We use these mock surveys to investigate various clustering statistics. The projected correlation functions for galaxies with different luminosities and types match the observations well on scales larger than about 3 Mpc/h. On smaller scales, however, the model overpredicts the clustering power by about a factor two. Furthermore, the model predicts pairwise velocity dispersions (PVD) that are about 400 km/s too high at projected pair separations of ~1 Mpc/h. A strong velocity bias in massive haloes can reduce the predicted PVD to the observed level, but does not help to resolve the over-prediction of clustering power on small scales. Consistent results can be obtained within the standard LCDM model only when the mass-to-light ratio of clusters is almost twice the observed value. Alternatively, a LCDM model with sigma_8=0.75 can also reproduce the observational results.

We use the observed abundance and clustering of galaxies from the 2dFGRS to determine the matter density Omega_m and the linear amplitude of mass fluctuations sigma_8. We use a method based on the conditional luminosity function, which allows straightforward computation of the luminosity dependent bias of galaxies. In addition, it allows the inclusion of constraints on the redshift space distortion parameter beta=Omega_m^{0.6}/b, and yields average mass-to-light ratios as function of halo mass. Using only the luminosity function and the correlation lengths as function of luminosity we obtain constraints on Omega_m and sigma_8 that are in good agreement with COBE. Using additional constraints on the mass-to-light ratios of clusters and on the value of beta as determined from the 2dFGRS yields Omega_m = (0.27 +/- 0.06) and sigma_8 = (0.77 +/- 0.06). Adding further constraints from current CMB data, we obtain Omega_m = (0.25 +/- 0.04) and sigma_8 = (0.78 +/- 0.06). Thus, we find strong indications that both Omega_m and sigma_8 are significantly lower than their ``standard'' concordance values of 0.3 and 0.9, respectively. Finally, we show that the cosmologies favored here predict dark matter haloes that are significantly less centrally concentrated than for the standard LCDM concordance model. This may solve both the problem with the rotation curves of dwarf and LSB galaxies, as well as the problem of simultaneously matching the galaxy luminosity function and the Tully-Fisher zero-point. (Abridged)

(Abridged) Using data from the 2dFGRS we compute the conditional luminosity functions (CLFs) of early- and late-type galaxies. These functions give the average number of galaxies with luminosity L that reside in a halo of mass M, and are a powerful statistical tool to link the distribution of galaxies to that of dark matter haloes. We find that the average mass-to-light ratios of dark matter haloes have a minimum of about 100 Msun/Lsun around a halo mass of about 3x10^{11} Msun. In addition, haloes with M < 10^{10} Msun are virtually devoid of galaxies, while clusters have b_J band mass-to-light ratios in the range 500-1000 Msun/Lsun. Finally, the fact that early-type galaxies are more strongly clustered than late-type galaxies requires that the fraction of late-type galaxies is a strongly declining function of halo mass. We also compute two-point correlation functions as function of both luminosity and galaxy type, and find a weak (strong) luminosity dependence for the late (early) type galaxies. Finally we compare our CLFs with predictions from several semi-analytical models for galaxy formation. As long as these models accurately fit the 2dFGRS luminosity function the agreement with our predictions is remarkably good. This indicates that we have recovered a statistical description of how galaxies populate dark matter haloes which fits nicely within the standard framework for galaxy formation.

We use the conditional luminosity function (CLF), which gives the number of galaxies with luminosities in the range [L, L+dL] that reside in a halo of mass M, to link the distribution of galaxies to that of dark matter haloes. We seek the CLF that reproduces the galaxy luminosity function and the luminosity dependence of the galaxy clustering strength and test the models by comparing the resulting mass-to-light ratios (M/L) with constraints from the Tully-Fisher (TF) relation. We obtain a number of stringent constraints on both galaxy formation and cosmology. In particular, this method can break the degeneracy between Omega_0 and the power-spectrum normalization sigma_8, inherent in current weak-lensing and cluster-abundance studies. For flat LCDM cosmogonies with sigma_8 normalized by recent weak lensing observations, the best results are obtained for Omega_0~0.3; Omega_0 < 0.2 leads to too large galaxy correlation lengths, while Omega_0 > 0.4 gives too high M/L. The best-fit model for the LCDM concordance cosmology (Omega_0=0.3) predicts M/L that are slightly too high. We discuss a number of possible effects that might remedy this problem, including small modification of cosmological parameters, warm in stead of cold dark matter, systematic errors in current observational data, and the existence of dark galaxies. Finally we use the CLF to predict several statistics about the distribution of galaxy light in the local Universe.

McKay et al. (2002) have recently used measurements of weak galaxy-galaxy lensing in the Sloan Digital Sky Survey to estimate the cross correlation between galaxies and the projected dark matter density field. They derive a relation between aperture mass within a radius of $260\kpch$, $M_{260}$, and lens galaxy luminosity, that does not depend on galaxy luminosity, type or environment. In this paper, we study the cross-correlation between galaxies and dark matter using galaxy catalogs constructed from a high-resolution N-body simulation of a $\Lambda$CDM Universe. We show that our simulations reproduce the McKay et al. results reasonably well. In the simulation, $M_{260}$ is approximately equal to the halo virial mass for $L_*$ galaxies. $M_{260}$ overestimates the virial mass for fainter galaxies and underestimates it for brighter galaxies. We use the simulations to show that under certain circumstances the halo virial mass may be recovered by fitting an NFW profile to the projected galaxy-mass correlation function. If we apply our method to the observations, we find that $L_*$ galaxies typically reside in halos of $\sim 2\times 10^{12} h^{-1} \msun$, consistent with halo masses estimated from the observed Tully-Fisher relation. In the simulations, the halo virial mass scales with galaxy luminosity as $L^{1.5}$ for central galaxies in halos and for galaxies in low-density regions. For satellite galaxies, and for galaxies in high-density regions, there is no simple relation between galaxy luminosity and halo mass and care must be exercised when interpreting the lensing results.

In this paper, we develop a theory of redshift distortion of the galaxy power spectrum in the discrete wavelet transform (DWT) representation. Because the DWT power spectrum is dependent of both the scale and shape (configuration) of the decomposition modes, it is sensitive to distortion of shape of the field. On the other hand, the redshift distortion causes a shape distortion of distributions in real space with respect to redshift space. Therefore, the shape-dependent DWT power spectrum is useful to detect the effect of redshift distortion. We first established the mapping between the DWT power spectra in redshift and real space. The mapping depends on the redshift distortion effects of (1) bulk velocity, (2) selection function and (3) pairwise peculiar velocity. We then proposed $\beta$-estimators using the DWT off-diagonal power spectra. These $\beta$-estimators are model-free even when the non-linear redshift distortion effect is not negligible. Moreover, these estimators do not rely on the assumption of whether the pairwise velocity dispersion being scale-dependent. The tests with N-body simulation samples show that the proposed $\beta$-estimators can yield reliable measurements of $\beta$ with about 20% uncertainty for all popular dark matter models. We also develop an algorithm for reconstruction of the power spectrum in real space from the redshift distorted power spectrum. The numerical test also shows that the real power spectrum can be well recovered from the redshift distorted power spectrum.

In this paper, we develop the method of analyzing the velocity field of cosmic matter with a multiresolution decomposition. This is necessary in calculating the redshift distortion of power spectrum in the discrete wavelet transform (DWT) representation. We show that, in the DWT analysis, the velocity field can be described by discrete variables, which are given by assignment of the number density and velocity into the DWT modes. These DWT variables are complete and not redundant. In this scheme, the peculiar velocity and pairwise velocity of galaxies or particles are given by field variables. As a consequence, the velocity dispersion (VD) and pairwise velocity dispersion (PVD) are no longer measured by number-counting or pair-counting statistic, but with the ensemble of the field variables, and therefore, they are free from the bias due to the number-counting and pair-counting. We analyzed the VD and PVD of the velocity fields given by the N-body simulation for models of the SCDM, $\tau$CDM and $\Lambda$CDM. The spectrum (scale-dependence) of the VD and PVD show that the length scale of the two-point correlation of the velocity field is as large as few tens h$^{-1}$ Mpc. Although the VD and PVD show similar behavior in some aspects, they are substantially different from each other. The VD-to-PVD ratio shows the difference between the scale-dependencies of the VD and PVD. More prominent difference between the VD and PVD is shown by probability distribution function. The one-point distribution of peculiar velocity is approximately exponential, while the pairwise velocity's is lognormal, i.e. of long tail. This difference indicates that the cosmic velocity field is typically intermittent.

The power spectrum estimator based on the discrete wavelet transform (DWT) for 3-dimensional samples has been studied. The DWT estimator for multi-dimensional samples provides two types of spectra with respect to diagonal and off-diagonal modes, which are very flexible to deal with configuration-related problems in the power spectrum detection. With simulation samples and mock catalogues of the Las Campanas redshift survey (LCRS), we show (1) the slice-like geometry of the LCRS doesn't affect the off-diagonal power spectrum with ``slice-like'' mode; (2) the Poisson sampling with the LCRS selection function doesn't cause more than 1-$\sigma$ error in the DWT power spectrum; and (3) the powers of peculiar velocity fluctuations, which cause the redshift distortion, are approximately scale-independent. These results insure that the uncertainties of the power spectrum measurement are under control. The scatter of the DWT power spectra of the six strips of the LCRS survey is found to be rather small. It is less than 1-$\sigma$ of the cosmic variance of mock samples in the wavenumber range $0.1 < k < 2$ h Mpc$^{-1}$. To fit the detected LCRS diagonal DWT power spectrum with CDM models, we find that the best-fitting redshift distortion parameter $\beta$ is about the same as that obtained from the Fourier power spectrum. The velocity dispersions $\sigma_v$ for SCDM and $\Lambda$CDM models are also consistent with other $\sigma_v$ detections with the LCRS. A systematic difference between the best-fitting parameters of diagonal and off-diagonal power spectra has been significantly measured. This indicates that the off-diagonal power spectra are capable of providing information about the power spectrum of galaxy velocity field.