Convolutional Neutral Networks have been successfully applied in searching for strong lensing systems, leading to discoveries of new candidates from large surveys. On the other hand, systematic investigations about their robustness are still lacking. In this paper, we first construct a neutral network, and apply it to $r$-band images of Luminous Red Galaxies (LRGs) of the Kilo Degree Survey (KiDS) Data Release 3 to search for strong lensing systems. We build two sets of training samples, one fully from simulations, and the other one using the LRG stamps from KiDS observations as the foreground lens images. With the former training sample, we find 48 high probability candidates after human-inspection, and among them, 27 are newly identified. Using the latter training set, about 67\% of the aforementioned 48 candidates are also found, and there are 11 more new strong lensing candidates identified. We then carry out tests on the robustness of the network performance with respect to the variation of PSF. With the testing samples constructed using PSF in the range of 0.4 to 2 times of the median PSF of the training sample, we find that our network performs rather stable, and the degradation is small. We also investigate how the volume of the training set can affect our network performance by varying it from 0.1 millions to 0.8 millions. The output results are rather stable showing that within the considered range, our network performance is not very sensitive to the volume size.

We develop an analytic mass model for lensing galaxies, based on a broken power-law (BPL) density profile, which is a power-law profile with a mass deficit or surplus in the central region. Under the assumption of an elliptically symmetric surface mass distribution, the deflection angle and magnification can be evaluated analytically for this new model. We compute the theoretical prediction for various quantities, including the volume and surface mass density profiles of the galaxies, and the aperture and luminosity-weighted line-of-sight velocity dispersions, and compare them to those measured from the Illustris simulation. We find an excellent agreement between our model prediction and the simulation, which validates our modeling. The high efficiency and accuracy of our model manifests itself as a promising tool for studying properties of galaxies with strong lensing.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, in synergy with the SUDARE survey, is a deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). The observations for the CDFS field comprise about 4.38 deg$^2$ down to $r\sim26$ mag. The total on-sky time spans over four years in this field, distributed over four adjacent sub-fields. In this paper, we use the multi-epoch $r$-band imaging data to measure the variability of the detected objects and search for transients. We perform careful astrometric and photometric calibrations and point spread function (PSF) modeling. A new method, referring to as differential running-average photometry, is proposed to measure the light curves of the detected objects. With the method, the difference of PSFs between different epochs can be reduced, and the background fluctuations are also suppressed. Detailed uncertainty analysis and detrending corrections on the light curves are performed. We visually inspect the light curves to select variable objects, and present some objects with interesting light curves. Further investigation of these objects in combination with multi-band data will be presented in our forthcoming paper.

We apply the inverse-Gaussianization method proposed in \citealt{arXiv:1607.05007} to fast produce weak lensing convergence maps and investigate the peak statistics, including the peak height counts and peak steepness counts, in these mocks. We find that the distribution of peak height and steepness is in good agreement with the simulation. The difference is $\lesssim 20\%$ for these peak statistics in the maps at source redshift $z_s=1$. Besides, the loss of off-diagonal elements in peak covariance motivates us to consider the super sample variance in weak lensing peak statistics. We propose correction methods to effectively recover the (anti-)correlation among different bins by adding scatters in the mean value of these mocks. Finally, as an example of the application, we adopt the improved inverse-Gaussianization method with the above improvement to fast generate 40,000 mocks to calculate precision matrices between the power spectrum and peak statistics.

Weak lensing peak abundance analyses have been applied in different surveys and demonstrated to be a powerful statistics in extracting cosmological information complementary to cosmic shear two-point correlation studies. Future large surveys with high number densities of galaxies enable tomographic peak analyses. Focusing on high peaks, we investigate quantitatively how the tomographic redshift binning can enhance the cosmological gains. We also perform detailed studies about the degradation of cosmological information due to photometric redshift (photo-z) errors. We show that for surveys with the number density of galaxies $\sim40\,{\rm arcmin^{-2}}$, the median redshift $\sim1$, and the survey area of $\sim15000\,{\rm deg^{2}}$, the 4-bin tomographic peak analyses can reduce the error contours of $(\Omega_{{\rm m}},\sigma_{8})$ by a factor of $5$ comparing to 2-D peak analyses in the ideal case of photo-z error being absent. More redshift bins can hardly lead to significantly better constraints. The photo-z error model here is parametrized by $z_{{\rm bias}}$ and $\sigma_{{\rm ph}}$ and the fiducial values of $z_{{\rm bias}}=0.003$ and $\sigma_{{\rm ph}}=0.02$ is taken. We find that using tomographic peak analyses can constrain the photo-z errors simultaneously with cosmological parameters. For 4-bin analyses, we can obtain $\sigma(z_{{\rm bias}})/z_{{\rm bias}}\sim10\%$ and $\sigma(\sigma_{{\rm ph}})/\sigma_{{\rm ph}}\sim5\%$ without assuming priors on them. Accordingly, the cosmological constraints on $\Omega_{{\rm m}}$ and $\sigma_{8}$ degrade by a factor of $\sim2.2$ and $\sim1.8$, respectively, with respect to zero uncertainties on photo-z parameters. We find that the uncertainty of $z_{{\rm bias}}$ plays more significant roles in degrading the cosmological constraints than that of $\sigma_{{\rm ph}}$.

The Chinese Space Station Optical Survey (CSS-OS) is a planned full sky survey operated by the Chinese Space Station Telescope (CSST). It can simultaneously perform the photometric imaging and spectroscopic slitless surveys, and will probe weak and strong gravitational lensing, galaxy clustering, individual galaxies and galaxy clusters, active galactic nucleus (AGNs), and so on. It aims to explore the properties of dark matter and dark energy and other important cosmological problems. In this work, we focus on two main CSS-OS scientific goals, i.e. the weak gravitational lensing (WL) and galaxy clustering surveys. We generate the mock CSS-OS data based on the observational COSMOS and zCOSMOS catalogs. We investigate the constraints on the cosmological parameters from the CSS-OS using the Markov Chain Monte Carlo (MCMC) method. The intrinsic alignments, galaxy bias, velocity dispersion, and systematics from instrumental effects in the CSST WL and galaxy clustering surveys are also included, and their impacts on the constraint results are discussed. We find that the CSS-OS can improve the constraints on the cosmological parameters by a factor of a few (even one order of magnitude in the optimistic case), compared to the current WL and galaxy clustering surveys. The constraints can be further enhanced when performing joint analysis with the WL, galaxy clustering, and galaxy-galaxy lensing data. Therefore, the CSS-OS is expected to be a powerful survey for exploring the Universe. Since some assumptions may be still optimistic and simple, it is possible that the results from the real survey could be worse. We will study these issues in details with the help of simulations in the future.

The tilt, rotation, or offset of each CCD with respect to the focal plane, as well as the distortion of the focal plane itself, cause shape distortions to the observed objects, an effect typically known as field distortion (FD). We point out that FD provides a unique way of quantifying the accuracy of cosmic shear measurement. The idea is to stack the shear estimators from galaxies that share similar FD-induced shape distortions. Given that the latter can be calculated with parameters from astrometric calibrations, the accuracy of the shear estimator can be directly tested on real images. It provides a way to calibrate the multiplicative and additive shear recovery biases within the scientific data itself, without requiring simulations or any external data sets. We use the CFHTLenS images to demonstrate the accuracy of the Fourier_Quad shear recovery method. We highlight some details in our image processing pipeline, including background removal, source identification and deblending, astrometric calibration, star selection for PSF reconstruction, noise reduction, etc.. We show that in the shear ranges of -0.005 < g_1 < 0.005 and -0.008 < g_2 < 0.008, the multiplicative biases are at the level of < 0.04. Slight additive biases on the order of 5E-4 (6 sigma) are identified for sources provided by the official CFHTLenS catalog (not using its shear catalog), but are minor (4 sigma) for source catalog generated by our Fourier_Quad pipeline.

It is well-known that fitting the Color Magnitude Diagrams (CMDs) to the theoretical isochrones is the main method to determine star cluster ages. However, when the CMDs are not available, the Spectral Energy Distribution (SED)-fitting technique is the only other approach, although it suffers the age-metallicity-reddening degeneracy. In this work, we gather the ages, metallicities and masses of dozens of M31 star clusters from the CMD-fitting with HST images from the literatures for comparison. We check the reliability of the SED-fitting results with different models, i.e., Bruzual \& Charlot 2003 model (BC03), Galaxy Evolutionary Synthesis Models ({\sl GALEV}) and Advanced Stellar Population Synthesis (ASPS) for the simple stellar populations (SSPs) with single stars (ss)-SSP/binary star (bs)-SSPs models. The photometry bands includes the Galaxy Evolution Explorer {\sl GALEX} FUV/NUV bands, optical/near-infrared $UBVRIJHK$ bands, as well as the Wide-field Infrared Survey Explorer ({\sl WISE}) $W1$/$W2$ bands. The comparisons show that the SED-fitting ages agree well with the CMD-fitting ages, either with the fixed metallicity or with the free metallicity for both the BC03 and the {\sl GALEV} model. However, for the ASPS models, it seems that SED-fitting results are systematic older than the CMD ages, especially for the ages log t <9.0 (yr). The fitting also shows that the {\sl GALEX} FUV/NUV-band are more important than the {\sl WISE} $W1$/$W2$ for constraining the ages, which confirms the previous findings. We also derived the masses of our sample star clusters from the BC03 and {\sl GALEV} models and it is found that the values agree well with that in the literature.

We present our deep $Y$-band imaging data of a two square degree field within the F22 region of the VIMOS VLT Deep Survey. The observations were conducted using the WIRCam instrument mounted at the Canada--France--Hawaii Telescope (CFHT). The total on-sky time was 9 hours, distributed uniformly over 18 tiles. The scientific goals of the project are to select faint quasar candidates at redshift $z>2.2$, and constrain the photometric redshifts for quasars and galaxies. In this paper, we present the observation and the image reduction, as well as the photometric redshifts that we derived by combining our $Y$-band data with the CFHTLenS $u^*g'r'i'z'$ optical data and UKIDSS DXS $JHK$ near-infrared data. With $J$-band image as reference total $\sim$80,000 galaxies are detected in the final mosaic down to $Y$-band $5\sigma$ point source limiting depth of 22.86 mag. Compared with the $\sim$3500 spectroscopic redshifts, our photometric redshifts for galaxies with $z<1.5$ and $i'\lesssim24.0$ mag have a small systematic offset of $|\Delta{z}|\lesssim0.2$, 1$\sigma$ scatter $0.03<\sigma_{\Delta z} < 0.06$, and less than 4.0% of catastrophic failures. We also compare to the CFHTLenS photometric redshifts, and find that ours are more reliable at $z\gtrsim0.6$ because of the inclusion of the near-infrared bands. In particular, including the $Y$-band data can improve the accuracy at $z\sim 1.0-2.0$ because the location of the 4000\AA-break is better constrained. The $Y$-band images, the multi-band photometry catalog and the photometric redshifts are released at \url{http://astro.pku.edu.cn/astro/data/DYI.html}.

This study expands the coverage and improves the homogeneity of the distribution of MILES template stars in the parameter space, as well as extends the wavelength coverage of the template spectra to the far red beyond the Ca\,{\sc{ii}} triplet. To achieve this we have carried out a major observational campaign using the OMR long-slit spectrograph mounted on the NAOC 2.16\,m telescope and the YFOSC long-slit spectrograph mounted on the YNAO 2.4\,m telescope. The original sample is based on the MILES library, supplemented by 918 stars selected from PASTEL database. In total, 822 OMR and 1,324 YFOSC spectra, covering respectively the wavelength ranges $\lambda\lambda$3800--5180 and $\lambda\lambda$5150--9000, have been collected and reduced. The spectra have a mean resolution FWHM (full-width at half-maximum) of $\sim 3.3$\,{\AA} and are wavelength- and flux-calibrated to an accuracy of $\sim 20$\,km\,s$^{-1}$ and $\sim 5$\,per\,cent, respectively. The spectra are further corrected for systematic errors in the wavelength calibration to an accuracy of $\sim 4$\,km\,s$^{-1}$ by cross-correlating with the theoretical spectra. Almost all the spectra have an average signal to noise ratio (SNR) better than 100 per pixel. Combined with the MILES spectra, there are now 1,731, 1,542, 1,324 and 1,273 stars with spectra covering respectively $\lambda\lambda$3800--5180, $\lambda\lambda$3800--7500, $\lambda\lambda$5150--9000 and $\lambda\lambda$3800--9000. This paper describes our template star selection, the observation and data reduction, and presents the reduced spectra collected hitherto.

To investigate in more details of Stellar Abundance and Galactic Evolution (SAGE) and in a huge sample, we are performing a northern sky photometric survey named SAGES with the SAGE photometric system, which consists of 8 filters: Str\"omgren-$u$, SAGE-$v$, SDSS $g$, $r$, $i$, DDO-$51$, $H\alpha_{wide}$, and $H\alpha _{narrow}$, including three Sloan broadband filters, three intermediate-band filters and two narrow-band filters, and one newly-designed narrow-band filter. SAGES covers $\sim$12,000 square degrees of the northern sky with $\delta > -5 ^{\circ}$, excluding the Galactic disk ($|b|<10^{\circ}$) and the sky area of 12 hr $<$ R.A. $<$ 18\,hr. The photometric detection limit depth at signal-to-noise ratio $5\sigma$ can be as deep as $V\sim$20\,mag. The SAGES will produce a depth-uniformed photometric catalogue for $\sim$500 million stars with atmospheric parameters including effective temperature $T_{\rm eff}$, surface gravity log\,g, and metallicity [Fe/H], as well as interstellar extinction to each individual target. In this work, we will briefly introduce the SAGE photometric system, the SAGE survey, and a preliminary test field of the open cluster NGC\,6791 and around.

Sky survey is one of the most important motivations to improve the astrophysics development, especially when using new photometric bands. We are performing the SAGE (Stellar Abundance and Galactic Evolution) survey with a self-designed SAGE photometric system, which is composed of eight photometric bands. The project mainly aims to study the stellar atmospheric parameters of $\sim$0.5 billion stars in the $\sim12,000$ deg$^2$ of the northern sky, which mainly focuses on the Galactic sciences, as well as some extragalactic sciences. This work introduces the detailed data reduction process of the testing field NGC\,6791, including the data reduction of single-exposure image and stacking multi-exposure images, and properties of the final catalogue.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is a Guaranteed Time program carried out with the ESO/VST telescope to provide deep optical imaging over two 4 deg$^2$ patches of the sky centred on the CDFS and ES1 pointings. We present the cosmic shear measurement over the 4 deg$^2$ covering the CDFS region in the $r$-band using LensFit. Each of the four tiles of 1 deg$^2$ has more than one hundred exposures, of which more than 50 exposures passed a series of image quality selection criteria for weak lensing study. The $5\sigma$ limiting magnitude in $r$- band is 26.1 for point sources, which is $\sim$1 mag deeper than other weak lensing survey in the literature (e.g. the Kilo Degree Survey, KiDS, at VST). The photometric redshifts are estimated using the VOICE $u,g,r,i$ together with near-infrared VIDEO data $Y,J,H,K_s$. The mean redshift of the shear catalogue is 0.87, considering the shear weight. The effective galaxy number density is 16.35 gal/arcmin$^2$, which is nearly twice the one of KiDS. The performance of LensFit on such a deep dataset was calibrated using VOICE-like mock image simulations. Furthermore, we have analyzed the reliability of the shear catalogue by calculating the star-galaxy cross-correlations, the tomographic shear correlations of two redshift bins and the contaminations of the blended galaxies. As a further sanity check, we have constrained cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. For a flat $\Lambda$CDM model we have obtained $\Sigma_8$ = $\sigma_8(\Omega_m/0.3)^{0.5}$ = $0.68^{+0.11}_{-0.15}$.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is proposed to obtain deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). At present, the observations for the CDFS field have been completed, and comprise in total about 4.9 deg$^2$ down to $r_\mathrm{AB}$$\sim$26 mag. In the companion paper by Fu et al. (2018), we present the weak lensing shear measurements for $r$-band images with seeing $\le$ 0.9 arcsec. In this paper, we perform image simulations to calibrate possible biases of the measured shear signals. Statistically, the properties of the simulated point spread function (PSF) and galaxies show good agreements with those of observations. The multiplicative bias is calibrated to reach an accuracy of $\sim$3.0%. We study the bias sensitivities to the undetected faint galaxies and to the neighboring galaxies. We find that undetected galaxies contribute to the multiplicative bias at the level of $\sim$0.3%. Further analysis shows that galaxies with lower signal-to-noise ratio (SNR) are impacted more significantly because the undetected galaxies skew the background noise distribution. For the neighboring galaxies, we find that although most have been rejected in the shape measurement procedure, about one third of them still remain in the final shear sample. They show a larger ellipticity dispersion and contribute to $\sim$0.2% of the multiplicative bias. Such a bias can be removed by further eliminating these neighboring galaxies. But the effective number density of the galaxies can be reduced considerably. Therefore efficient methods should be developed for future weak lensing deep surveys.

The convergence peaks, constructed from galaxy shape measurement in weak lensing, is a powerful probe of cosmology as the peaks can be connected with the underlined dark matter haloes. However the capability of convergence peak statistic is affected by the noise in galaxy shape measurement, signal to noise ratio as well as the contribution from the projected mass distribution from the large-scale structures along the line of sight (LOS). In this paper we use the ray-tracing simulation on a curved sky to investigate the correspondence between the convergence peak and the dark matter haloes at the LOS. We find that, in case of no noise and for source galaxies at $z_{\rm s}=1$, more than $65\%$ peaks with $\text{SNR} \geq 3$ (signal to noise ratio) are related to more than one massive haloes with mass larger than $10^{13} {\rm M}_{\odot}$. Those massive haloes contribute $87.2\%$ to high peaks ($\text{SNR} \geq 5$) with the remaining contributions are from the large-scale structures. On the other hand, the peaks distribution is skewed by the noise in galaxy shape measurement, especially for lower SNR peaks. In the noisy field where the shape noise is modelled as a Gaussian distribution, about $60\%$ high peaks ($\text{SNR} \geq 5$) are true peaks and the fraction decreases to $20\%$ for lower peaks ($ 3 \leq \text{SNR} < 5$). Furthermore, we find that high peaks ($\text{SNR} \geq 5$) are dominated by very massive haloes larger than $10^{14} {\rm M}_{\odot}$.

The Beijing Faint Object Spectrograph and Camera (BFOSC) is one of the most important instruments of the 2.16-m telescope of the Xinglong Observatory. Every year there are ~ 20 SCI-papers published based on the observational data of this telescope. In this work, we have systemically measured the total efficiency of the BFOSC of the 2.16-m reflector, based on the observations of two ESO flux standard stars. We have obtained the total efficiencies of the BFOSC instrument of different grisms with various slit widths in almost all ranges, and analysed the factors which effect the efficiency of telescope and spectrograph. For the astronomical observers, the result will be useful for them to select a suitable slit width, depending on their scientific goals and weather conditions during the observation; For the technicians, the result will help them systemically find out the real efficiency of telescope and spectrograph, and further to improve the total efficiency and observing capacity of the telescope technically.

The 85-cm telescope at the Xinglong station is a well-operated prime focus system with high science outputs. The telescope has been upgraded since 2014 with new corrector, filters and camera, which are provided by Beijing Normal University (BNU). The filter set is Johnson-Cousins UBVRI system. We report the test results of the new system including the bias, dark current, linearity, gain and readout noise of the CCD camera . Then we derive accurate instrumental calibration coefficients in UBVRI bands with Landolt standard stars in the photometric nights. Finally, we give the limiting magnitudes with various exposure time and signal-to-noise ratio for observers as references.

High peaks in weak lensing (WL) maps originate dominantly from the lensing effects of single massive halos. Their abundance is therefore closely related to the halo mass function and thus a powerful cosmological probe. On the other hand, however, besides individual massive halos, large-scale structures (LSS) along lines of sight also contribute to the peak signals. In this paper, with ray tracing simulations, we investigate the LSS projection effects. We show that for current surveys with a large shape noise, the stochastic LSS effects are subdominant. For future WL surveys with source galaxies having a median redshift $z_{\mathrm{med}}\sim1$ or higher, however, they are significant. For the cosmological constraints derived from observed WL high peak counts, severe biases can occur if the LSS effects are not taken into account properly. We extend the model of \citet{Fan2010} by incorporating the LSS projection effects into the theoretical considerations. By comparing with simulation results, we demonstrate the good performance of the improved model and its applicability in cosmological studies.

By using the Hectospec 6.5 m Multiple Mirror Telescope (MMT) and the 2.16 m telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC), we obtained 188 high signal-to-noise ratio (S/N) spectra of HII regions in the nearby galaxy M101, which are the largest spectroscopic sample of HII regions for this galaxy so far. These spectra cover a wide range of regions on M101, which enables us to analyze two dimensional distributions of its physical properties. The physical parameters are derived from emission lines or stellar continuum, including stellar population age, electron temperature, oxygen abundance and etc. The oxygen abundances are derived using two empirical methods based on O3N2 and R$_{23}$ indicators, as well as the direct Te method when OIII$\lambda4363$ is available. By applying the harmonic decomposition analysis to the velocity field, we obtained line-of-sight rotation velocity of 71 km s$^{-1}$ and a position angle of 36 degree. The stellar age profile shows an old stellar population in galaxy center and a relative young stellar population in outer regions, suggesting an old bulge and a young disk. Oxygen abundance profile exhibits a clear break at $\sim$18 kpc, with a gradient of $-$0.0364 dex kpc$^{-1}$ in the inner region and $-$0.00686 dex kpc$^{-1}$ in the outer region. Our results agree with the "inside-out" disk growth scenario of M101.

In this paper, we will give a general introduction to the project of Ali CMB Polarization Telescope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics (IHEP) and has involved many different institutes in China. It is the first ground-based Cosmic Microwave Background (CMB) polarization experiment in China and an integral part of China's Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primordial gravitational waves (PGWs) originated from the very early Universe. The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the worldwide highest ground-based CMB observatory and open a new window for probing PGWs in northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of detectors more than 20,000. Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar ratio $r$ by one order of magnitude with 3 years' observation. Besides the PGWs, the AliCPT will also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT symmetry. We show 3 years' observation will improve the current limit by two order of magnitude.

We report the result of a faint quasar survey in a one square degree field. The aim is to test the Y-K/g-z and J-K/i-Y color selection criteria for quasars at faint magnitude, to obtain a complete sample of quasars based on deep optical and near-infrared color-color selection, and to measure the faint end of quasar luminosity function (QLF) over a wide redshift range. We carried out a quasar survey based on the Y-K/g-z and J-K/i-Y quasar selection criteria, using the deep Y-band data obtained from our CFHT/WIRCam Y-band images in a two-degree field within the F22 field of the VIMOS VLT deep survey, optical co-added data from Sloan Digital Sky Survey Stripe 82 and deep near-infrared data from the UKIDSS Deep Extragalactic Survey in the same field. We discovered 25 new quasars at 0.5 < z < 4.5 and i < 22.5 mag within one square degree field. The survey significantly increases the number of faint quasars in this field, especially at z ~ 2-3. It confirms that our color selections are highly complete in a wide redshift range (z < 4.5), especially over the quasar number density peak at z ~ 2-3, even for faint quasars. Combining all previous known quasars and new discoveries, we construct a sample with 109 quasars, and measure the binned QLF and parametric QLF. Although the sample is small, our results agree with a pure luminosity evolution at lower redshift and luminosity evolution and density evolution model at redshift z > 2.5.

Accurate reconstruction of the spatial distributions of the Point Spread Function (PSF) is crucial for high precision cosmic shear measurements. Nevertheless, current methods are not good at recovering the PSF fluctuations of high spatial frequencies. In general, the residual PSF fluctuations are spatially correlated, therefore can significantly contaminate the correlation functions of the weak lensing signals. We propose a method to correct for this contamination statistically, without any assumptions on the PSF and galaxy morphologies or their spatial distribution. We demonstrate our idea with the data from the W2 field of CFHTLenS.

In this work, we present a novel centroiding method based on Fourier space Phase Fitting(FPF) for Point Spread Function(PSF) reconstruction. We generate two sets of simulations to test our method. The first set is generated by GalSim with elliptical Moffat profile and strong anisotropy which shifts the center of the PSF. The second set of simulation is drawn from CFHT i band stellar imaging data. We find non-negligible anisotropy from CFHT stellar images, which leads to $\sim$0.08 scatter in unit of pixels using polynomial fitting method Vakili and Hogg (2016). And we apply FPF method to estimate the centroid in real space, this scatter reduces to $\sim$0.04 in SNR=200 CFHT like sample. In low SNR (50 and 100) CFHT like samples, the background noise dominates the shifting of the centroid, therefore the scatter estimated from different methods are similar. We compare polynomial fitting and FPF using GalSim simulation with optical anisotropy. We find that in all SNR$\sim$50, 100 and 200) samples, FPF performs better than polynomial fitting by a factor of $\sim$3. In general, we suggest that in real observations there are anisotropy which shift the centroid, and FPF method is a better way to accurately locate it.

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

This paper is the first of a series of papers constraining cosmological parameters with weak lensing peak statistics using $\sim 450~\rm deg^2$ of imaging data from the Kilo Degree Survey (KiDS-450). We measure high signal-to-noise ratio (SNR: $\nu$) weak lensing convergence peaks in the range of $3<\nu<5$, and employ theoretical models to derive expected values. These models are validated using a suite of simulations. We take into account two major systematic effects, the boost factor and the effect of baryons on the mass-concentration relation of dark matter haloes. In addition, we investigate the impacts of other potential astrophysical systematics including the projection effects of large scale structures, intrinsic galaxy alignments, as well as residual measurement uncertainties in the shear and redshift calibration. Assuming a flat $\Lambda$CDM model, we find constraints for $S_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.5}=0.746^{+0.046}_{-0.107}$ according to the degeneracy direction of the cosmic shear analysis and $\Sigma_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.38}=0.696^{+0.048}_{-0.050}$ based on the derived degeneracy direction of our high-SNR peak statistics. The difference between the power index of $S_{\rm 8}$ and in $\Sigma_{\rm 8}$ indicates that combining the two probes has the potential to break the degeneracy in $\sigma_{\rm 8}$ and $\Omega_{\rm m}$. Our results are consistent with the cosmic shear tomographic correlation analysis of the same dataset and $\sim 2\sigma$ lower than the Planck 2016 results.

Spectral energy distribution (SED) fitting of stellar population synthesis models is an important and popular way to constrain the physical parameters ---e.g., the ages, metallicities, masses for stellar population analysis. The previous works suggest that both blue-bands and red-bands photometry works for the SED-fitting. Either blue-domained or red-domained SED-fitting usually lead to the unreliable or biased results. Meanwhile, it seems that extending the wavelength coverage could be helpful. Since the Galaxy Evolution Explorer ({\it GALEX}) and Wide-field Infrared Survey Explorer (WISE) provide the FUV/NUV and mid-infrared $W1$/$W2$ band data, we extend the SED-fitting to a wider wavelength coverage. In our work, we analyzed the effect of adding the FUV/NUV and $W1$/$W2$ band to the optical and near-infrared $UBVRIJHK$ bands for the fitting with the Bruzual \& Charlot 2003 (BC03) models and {\sc galev} models. It is found that the FUV/NUV bands data affect the fitting results of both ages and metallicities much more significantly than that of the WISE $W1$/$W2$ band with the BC03 models. While for the {\sc galev} models, the effect of the WISE $W1$/$W2$ band for the metallicity fitting seems comparable to that of {\it GALEX} FUV/NUV bands, but for age the effect of the $W1$/$W2$ band seems less crucial than that of the FUV/NUV bands. Thus we conclude that the {\it GALEX} FUV/NUV bands are more crucial for the SED-fitting of ages and metallicities, than the other bands, and the high-quality UV data (with high photometry precision) are required.

Reconstruction of the point spread function (PSF) is a critical process in weak lensing measurement. We develop a real-data based and galaxy-oriented pipeline to compare the performances of various PSF reconstruction schemes. Making use of a large amount of the CFHTLenS data, the performances of three classes of interpolating schemes - polynomial, Kriging, and Shepard - are evaluated. We find that polynomial interpolations with optimal orders and domains perform the best. We quantify the effect of the residual PSF reconstruction error on shear recovery in terms of the multiplicative and additive biases, and their spatial correlations using the shear measurement method of Zhang et al. (2015). We find that the impact of PSF reconstruction uncertainty on the shear-shear correlation can be significantly reduced by cross correlating the shear estimators from different exposures. It takes only 0.2 stars (SNR > 100) per square arcmin on each exposure to reach the best performance of PSF interpolation, a requirement that is generally satisfied in the CFHTlenS data.

We present a new tool for colour-magnitude diagram (CMD) studies, $Powerful~CMD$. This tool is built on the basis of the advanced stellar population synthesis (ASPS) model, in which single stars, binary stars, rotating stars, and star formation history have been taken into account. Via $Powerful~CMD$, the distance modulus, colour excess, metallicity, age, binary fraction, rotating star fraction, and star formation history of star clusters can be determined simultaneously from observed CMDs. The new tool is tested via both simulated and real star clusters. Five parameters of clusters NGC6362, NGC6652, NGC6838 and M67 are determined and compared to other works. It is shown that this tool is useful for CMD studies, in particular for those with the data of the Hubble Space Telescope (HST). Moreover, we find that the inclusion of binaries in theoretical stellar population models may lead to smaller colour excess compared to the case of single star population models.

We study the integral Galactic extinction and reddening based on the galaxy catalog of the South Galactic Cap U-band Sky Survey (SCUSS), where $u$ band galaxy number counts and $u-r$ color distribution are used to derive the Galactic extinction and reddening respectively. We compare these independent statistical measurements with the reddening map of \citet{Schlegel1998}(SFD) and find that both the extinction and reddening from the number counts and color distribution are in good agreement with the SFD results at low extinction regions ($E(B-V)^{SFD}<0.12$ mag). However, for high extinction regions ($E(B-V)^{SFD}>0.12$ mag), the SFD map overestimates the Galactic reddening systematically, which can be approximated by a linear relation $\Delta E(B-V)= 0.43[E(B-V)^{SFD}-0.12$]. By combing the results of galaxy number counts and color distribution together, we find that the shape of the Galactic extinction curve is in good agreement with the standard $R_V=3.1$ extinction law of \cite{ODonnell1994}.

We present and analyze the possibility of using optical ${\it u}$-band luminosities to estimate star-formation rates (SFRs) of galaxies based on the data from the South Galactic Cap ${\it u }$ band Sky Survey (SCUSS), which provides a deep ${\it u}$-band photometric survey covering about 5000 $deg^2$ of the South Galactic Cap. Based on two samples of normal star-forming galaxies selected by the BPT diagram, we explore the correlations between ${\it u}$-band, H$\alpha$, and IR luminosities by combing SCUSS data with the Sloan Digital Sky Survey (SDSS) and ${\it Wide}$-${\it field\ Infrared\ Survey\ Explorer}$ (${\it WISE}$). The attenuation-corrected ${\it u}$-band luminosities are tightly correlated with the Balmer decrement-corrected H$\alpha$ luminosities with an rms scatter of $\sim$ 0.17 dex. The IR-corrected ${\it u }$ luminosities are derived based on the correlations between the attenuation of ${\it u}$-band luminosities and ${\it WISE}$ 12 (or 22) $\mu$m luminosities, and then calibrated with the Balmer-corrected H$\alpha$ luminosities. The systematic residuals of these calibrations are tested against the physical properties over the ranges covered by our sample objects. We find that the best-fitting nonlinear relations are better than the linear ones and recommended to be applied in the measurement of SFRs. The systematic deviations mainly come from the pollution of old stellar population and the effect of dust extinction; therefore, a more detailed analysis is needed in the future work.

Application of fitting techniques to obtain physical parameters---such as ages, metallicities, and $\alpha$-element to iron ratios---of stellar populations is an important approach to understand the nature of both galaxies and globular clusters (GCs). In fact, fitting methods based on different underlying models may yield different results, and with varying precision. In this paper, we have selected 22 confirmed M31 GCs for which we do not have access to previously known spectroscopic metallicities. Most are located at approximately one degree (in projection) from the galactic center. We performed spectroscopic observations with the 6.5 m MMT telescope, equipped with its Red Channel Spectrograph. Lick/IDS absorption-line indices, radial velocities, ages, and metallicities were derived based on the $\rm EZ\_Ages$ stellar population parameter calculator. We also applied full spectral fitting with the ULySS code to constrain the parameters of our sample star clusters. In addition, we performed $\chi^2_{\rm min}$ fitting of the clusters' Lick/IDS indices with different models, including the Bruzual & Charlot models (adopting Chabrier or Salpeter stellar initial mass functions and 1994 or 2000 Padova stellar evolutionary tracks), the GALEV, and the Thomas et al. models. For comparison, we collected their $UVBRIJK$ photometry from the Revised Bologna Catalogue (v.5) to obtain and fit the GCs' spectral energy distributions (SEDs). Finally, we performed fits using a combination of Lick/IDS indices and SEDs. The latter results are more reliable and the associated error bars become significantly smaller than those resulting from either our Lick/IDS indices-only or our SED-only fits.

We present spectroscopic measurements on the night sky of Xinglong Observatory for a period of 12 years from 2004 to 2015. The spectra were obtained on moonless clear nights using the OMR spectrograph mounted on a 2.16-m reflector with a wavelength coverage of 4000-7000A. The night sky spectrum shows the presence of emission lines from Hg I and Na I due to local artificial sources, along with the atmospheric emission lines, i.e., O I and OH molecules, indicating the existence of light pollution. We have monitored the night sky brightness during the whole night and found some decrement in the sky brightness with time, but the change is not significant. Also, we monitored the light pollution level in different azimuthal directions and found that the influence of light pollution from the direction of Beijing is stronger compared with that from the direction of Tangshan and other areas. An analysis of night sky spectra for the entire data set suggested that the zenith sky brightness of Xinglong Observatory has brightened by about 0.5 mag arcsec -2 in the V and B bands from 2004 to 2015. We recommend consecutive spectroscopic measurements of the night sky brightness at Xinglong Observatory in the future, not only for monitoring but also for scientific reference.

In this Letter, we report the observational constraints on the Hu-Sawicki $f(R)$ theory derived from weak lensing peak abundances, which are closely related to the mass function of massive halos. In comparison with studies using optical or x-ray clusters of galaxies, weak lensing peak analyses have the advantages of not relying on mass-baryonic observable calibrations. With observations from the Canada-France-Hawaii-Telescope Lensing Survey, our peak analyses give rise to a tight constraint on the model parameter $|f_{R0}|$ for $n=1$. The $95\%$ CL limit is $\log_{10}|f_{R0}| < -4.82$ given WMAP9 priors on $(\Omega_{\rm m}, A_{\rm s})$. With Planck15 priors, the corresponding result is $\log_{10}|f_{R0}| < -5.16$.

We selected 82 u-band variable objects based on the u-band photometry data from SCUSS and SDSS, in the field of LAMOST Complete Spectroscopic Survey of Pointing Area at Southern Galactic Cap. The magnitude variation of the targets is restricted to larger than 0.2 mag and limiting magnitude down to 19.0 mag in u-band. According to the spectra from LAMOST, there are 11 quasars with red-shift between 0.4 and 1.8, 60 variable stars and 11 unidentified targets. The variable stars include one active M-dwarf with a series of Balmer emission lines, seven Horizontal Branch stars containing six RR Lyrae stars matching with SIMBAD, and one giant, one AGB star and two RR Lyrae candidates by different colour selections. All these variable stars mainly locate near the main sequence in the g-r verse u-g diagram. The quasars are well distinguished from stars by both u-g colour and variation in u-band.

The Xinglong 2.16-m reflector is the first 2-meter class astronomical telescope in China. It was jointly designed and built by the Nanjing Astronomical Instruments Factory (NAIF), Beijing Astronomical Observatory (now National Astronomical Observatories, Chinese Academy of Sciences, NAOC) and Institute of Automation, Chinese Academy of Sciences in 1989. It is Ritchey-Chr\'{e}tien (R-C) reflector on an English equatorial mount and the effective aperture is 2.16 meters. It had been the largest optical telescope in China for $\sim18$ years until the Guoshoujing Telescope (also called Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST) and the Lijiang 2.4-m telescope were built. At present, there are three main instruments on the Cassegrain focus available: the Beijing Faint Object Spectrograph and Camera (BFOSC) for direct imaging and low resolution ($R\sim500-2000$) spectroscopy, the spectrograph made by Optomechanics Research Inc. (OMR) for low resolution spectroscopy (the spectral resolutions are similar to those of BFOSC) and the fiber-fed High Resolution Spectrograph (HRS, $R\sim30000-65000$). The telescope is widely open to astronomers all over China as well as international astronomical observers. Each year there are more than 40 ongoing observing projects, including 6-8 key projects. Recently, some new techniques and instruments (e.g., astro-frequency comb calibration system, polarimeter and adaptive optics) have been or will be tested on the telescope to extend its observing abilities.

The SCUSS is a deep $u$-band imaging survey in the south Galactic cap using the 2.3m Bok telescope. The survey observations were completed in the end of 2013, covering an area of about 5000 square degrees. We release the data in the region with an area of about 4000 deg$^2$ that is mostly covered by the Sloan digital sky survey. The data products contain calibrated single-epoch images, stacked images, photometric catalogs, and a catalog of star proper motions derived by Peng et al, 2015. The median seeing and magnitude limit ($5\sigma$) are about 2".0 and 23.2 mag, respectively. There are about 8 million objects having measurements of absolute proper motions. All the data and related documentations can be accessed through the SCUSS data release website of \url{http://batc.bao.ac.cn/Uband/data.html}.

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

Based on weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), in this paper we study the mass-concentration ($M$-$c$) relation for $\sim 200$ redMaPPer clusters in the fields. We extract the $M$-$c$ relation by measuring the density profiles of individual clusters instead of using stacked weak lensing signals. By performing Monte Carlo simulations, we demonstrate that although the signal-to-noise ratio for each individual cluster is low, the unbiased $M$-$c$ relation can still be reliably derived from a large sample of clusters by carefully taking into account the impacts of shape noise, cluster center offset, dilution effect from member or foreground galaxies, and the projection effect. Our results show that within error bars the derived $M$-$c$ relation for redMaPPer clusters is in agreement with simulation predictions. There is a weak deviation in that the halo concentrations calibrated by Monte Carlo simulations are somewhat higher than that predicted from ${\it Planck}$ cosmology.

Whether or not the initial star cluster mass function is established through a universal, galactocentric-distance-independent stochastic process, on the scales of individual galaxies, remains an unsolved problem. This debate has recently gained new impetus through the publication of a study that concluded that the maximum cluster mass in a given population is not solely determined by size-of-sample effects. Here, we revisit the evidence in favor and against stochastic cluster formation by examining the young ($\lesssim$ a few $\times 10^8$ yr-old) star cluster mass--galactocentric radius relation in M33, M51, M83, and the Large Magellanic Cloud. To eliminate size-of-sample effects, we first adopt radial bin sizes containing constant numbers of clusters, which we use to quantify the radial distribution of the first- to fifth-ranked most massive clusters using ordinary least-squares fitting. We supplement this analysis with an application of quantile regression, a binless approach to rank-based regression taking an absolute-value-distance penalty. Both methods yield, within the $1\sigma$ to $3\sigma$ uncertainties, near-zero slopes in the diagnostic plane, largely irrespective of the maximum age or minimum mass imposed on our sample selection, or of the radial bin size adopted. We conclude that, at least in our four well-studied sample galaxies, star cluster formation does not necessarily require an environment-dependent cluster formation scenario, which thus supports the notion of stochastic star cluster formation as the dominant star cluster-formation process within a given galaxy.

The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic Cap, using the 90Prime wide-field imager on the 2.3m Bok telescope at Kitt Peak. The survey observations started in 2010 and ended in 2013. The final survey area is about 5000 deg2 with a median 5-sigma point source limiting magnitude of about 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS photometric catalog with those of the SDSS and CFHTLS.

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

Based on SDSS g, r and SCUSS (South Galactic Cap of u-band Sky Survey) $u$ photometry, we develop a photometric calibration for estimating the stellar metallicity from $u-g$ and $g-r$ colors by using the SDSS spectra of 32,542 F- and G-type main sequence stars, which cover almost $3700$ deg$^{2}$ in the south Galactic cap. The rms scatter of the photometric metallicity residuals relative to spectrum-based metallicity is $0.14$ dex when $g-r<0.4$, and $0.16$ dex when $g-r>0.4$. Due to the deeper and more accurate magnitude of SCUSS $u$ band, the estimate can be used up to the faint magnitude of $g=21$. This application range of photometric metallicity calibration is wide enough so that it can be used to study metallicity distribution of distant stars. In this study, we select the Sagittarius (Sgr) stream and its neighboring field halo stars in south Galactic cap to study their metallicity distribution. We find that the Sgr stream at the cylindrical Galactocentric coordinate of $R\sim 19$ kpc, $\left| z\right| \sim 14$ kpc exhibits a relative rich metallicity distribution, and the neighboring field halo stars in our studied fields can be modeled by two-Gaussian model, with peaks respectively at [Fe/H]$=-1.9$ and [Fe/H]$=-1.5$.

The Chinese Small Telescope ARray (CSTAR) has observed an area around the Celestial South Pole at Dome A since 2008. About $20,000$ light curves in the i band were obtained lasting from March to July, 2008. The photometric precision achieves about 4 mmag at i = 7.5 and 20 mmag at i = 12 within a 30 s exposure time. These light curves are analyzed using Lomb--Scargle, Phase Dispersion Minimization, and Box Least Squares methods to search for periodic signals. False positives may appear as a variable signature caused by contaminating stars and the observation mode of CSTAR. Therefore the period and position of each variable candidate are checked to eliminate false positives. Eclipsing binaries are removed by visual inspection, frequency spectrum analysis and locally linear embedding technique. We identify 53 eclipsing binaries in the field of view of CSTAR, containing 24 detached binaries, 8 semi-detached binaries, 18 contact binaries, and 3 ellipsoidal variables. To derive the parameters of these binaries, we use the Eclipsing Binaries via Artificial Intelligence (EBAI) method. The primary and the secondary eclipse timing variations (ETVs) for semi-detached and contact systems are analyzed. Correlated primary and secondary ETVs confirmed by false alarm tests may indicate an unseen perturbing companion. Through ETV analysis, we identify two triple systems (CSTAR J084612.64-883342.9 and CSTAR J220502.55-895206.7). The orbital parameters of the third body in CSTAR J220502.55-895206.7 are derived using a simple dynamical model.

We derived constraints on cosmological parameters using weak lensing peak statistics measured on the $\sim130~{\rm deg}^2$ of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). This analysis demonstrates the feasibility of using peak statistics in cosmological studies. For our measurements, we considered peaks with signal-to-noise ratio in the range of $\nu=[3,6]$. For a flat $\Lambda$CDM model with only $(\Omega_{\rm m}, \sigma_8)$ as free parameters, we constrained the parameters of the following relation $\Sigma_8=\sigma_8(\Omega_{\rm m}/0.27)^{\alpha}$ to be: $\Sigma_8=0.82 \pm 0.03 $ and $\alpha=0.43\pm 0.02$. The $\alpha$ value found is considerably smaller than the one measured in two-point and three-point cosmic shear correlation analyses, showing a significant complement of peak statistics to standard weak lensing cosmological studies. The derived constraints on $(\Omega_{\rm m}, \sigma_8)$ are fully consistent with the ones from either WMAP9 or Planck. From the weak lensing peak abundances alone, we obtained marginalised mean values of $\Omega_{\rm m}=0.38^{+0.27}_{-0.24}$ and $\sigma_8=0.81\pm 0.26$. Finally, we also explored the potential of using weak lensing peak statistics to constrain the mass-concentration relation of dark matter halos simultaneously with cosmological parameters.

We present the discovery of eight quasars at z~6 identified in the Sloan Digital Sky Survey (SDSS) overlap regions. Individual SDSS imaging runs have some overlap with each other, leading to repeat observations over an area spanning >4000 deg^2 (more than 1/4 of the total footprint). These overlap regions provide a unique dataset that allows us to select high-redshift quasars more than 0.5 mag fainter in the z band than those found with the SDSS single-epoch data. Our quasar candidates were first selected as i-band dropout objects in the SDSS imaging database. We then carried out a series of follow-up observations in the optical and near-IR to improve photometry, remove contaminants, and identify quasars. The eight quasars reported here were discovered in a pilot study utilizing the overlap regions at high galactic latitude (|b|>30 deg). These quasars span a redshift range of 5.86<z<6.06 and a flux range of 19.3<z_AB<20.6 mag. Five of them are fainter than z_AB=20 mag, the typical magnitude limit of z~6 quasars used for the SDSS single-epoch images. In addition, we recover eight previously known quasars at z~6 that are located in the overlap regions. These results validate our procedure for selecting quasar candidates from the overlap regions and confirming them with follow-up observations, and provide guidance to a future systematic survey over all SDSS imaging regions with repeat observations.

Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent.

The South Galactic Cap $u$-band Sky Survey (SCUSS) provides a deep $u$-band imaging of about 5000 deg$^2$ in south Galactic cap. It is about 1.5 mag deeper than the SDSS $u$-band. In this paper we evaluate the capability of quasar selection using both SCUSS and SDSS data, based on considerations of the deep SCUSS $u$-band imaging and two-epoch $u$-band variability. We find that the combination of the SCUSS $u$-band and the SDSS $griz$ band allows us to select more faint quasars and more quasars at redshift around 2.2 than the selection only with the SDSS $ugriz$ data. Quasars have significant $u$-band variabilities. The fraction of quasars with large two-epoch variability is much higher than that of stars. The selection by variability can select both low-redshift quasars with ultraviolet excess and mid-redshift ($2 < z <3.5$) quasars where quasar selection by optical colors is inefficient. The above two selections are complementary and make full use of the SCUSS u-band advantages.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

With the halo catalog from the {\it Millennium Simulation}, we analyze the weak-lensing measured density profiles for clusters of galaxies, paying attention to the determination of the concentration-mass ($c$-$M$) relation which can be biased by the center offset, selection effect, and shape noise from intrinsic ellipticities of background galaxies. Several different methods of locating the center of a cluster from weak-lensing effects alone are explored. We find that, for intermediate redshift clusters, the highest peak from our newly proposed two-scale smoothing method applied to the reconstructed convergence field, first with a smoothing scale of $2^\prime$ and then $0\overset{\prime}{.}5$, corresponds best to the true center. Assuming the parameterized Navarro-Frenk-White profile, we fit the reduced tangential shear signals around different centers identified by different methods. It is shown that, for the ensemble median values, a center offset larger than one scale radius $r_s$ can bias the derived mass and concentration significantly lower than the true values, especially for low-mass halos. However, the existence of noise can compensate for the offset effect and reduce the systematic bias, although the scatter of mass and concentration becomes considerably larger. Statistically, the bias effect of center offset on the $c$-$M$ relation is insignificant if an appropriate center finding method is adopted. On the other hand, noise from intrinsic ellipticities can bias the $c$-$M$ relation derived from a sample of weak-lensing analyzed clusters if a simple $\chi^2$ fitting method is used. To properly account for the scatter and covariance between $c$ and $M$, we apply a Bayesian method to improve the statistical analysis of the $c$-$M$ relation. It is shown that this new method allows us to derive the $c$-$M$ relation with significantly reduced biases.

We present a comprehensive study of star formation toward the H II region S155. Star-formation activities therein were investigated based on multi-wavelength data from optical to the far-infrared. The surface density distribution of selected 2MASS sources toward S155 indicates the existence of a compact cluster, which is spatially consistent with the position of the exciting source of the H II region, HD 217086. A sample of more than 200 excessive emission sources in the infrared were selected based on their 2MASS color indices. The spatial distribution of the sample sources reveals the existence of three young sub-clusters in this region, among which, sub-cluster A is spatially coincident with the bright rim of the H II region. In addition, photometric data from the WISE survey were used to identify and classify young stellar objects (YSOs). To further explore the evolutionary stages of the candidate YSOs, we fit the spectral energy distribution (SEDs) of 44 sources, which led to the identification of 14 Class I, 27 Class II, and 3 Class III YSOs. The spatial distribution of the classified YSOs at different evolutionary stages presents a spatio-temporal gradient, which is consistent with a scenario of sequential star formation. On the other hand, Herschel PACS observations toward the interface between S155 and the ambient molecular cloud disclose an arc-shaped dust layer, the origin of which could be attributed to the UV dissipation from the early type stars e.g. HD 217061 in S155. Four dusty cores were revealed by the Herschel data, which hints for new generations of star formation.

Based on SDSS and South Galactic Cap of u-band Sky Survey (SCUSS) early data, we use star counts method to estimate the Galactic structure parameters in an intermediate latitude with 10,180 main-sequence (MS) stars in absolute magnitude interval of $4 \leq M_r \leq 13$. We divide the absolute magnitude into five intervals:$4 \leq M_r < 5$, $5 \leq M_r < 6$, $6 \leq M_r < 8$, $8 \leq M_r < 10$, $10 \leq M_r \leq 13$, and estimate the Galactic structure parameters in each absolute magnitude interval to explore their possible variation with the absolute magnitude. Our study shows the parameters depend on absolute magnitude. For the thin disk, the intrinsic faint MS stars have large local space density and they tend to stay close to the Galactic plane. A plausible explanation is that faint MS stars with long lifetime experience long gravitational interaction time result in a short scaleheight. However, for the thick disk, the parameters show a complex trend with absolute magnitude, which may imply the complicated original of the thick disk. For the halo, the intrinsic faint MS stars have large local density and small axial ratio, which indicate a flattened inner halo and a more spherical outer halo.

The Chinese Small Telescope ARray (CSTAR) is a group of four identical, fully automated, static 14.5 cm telescopes. CSTAR is located at Dome A, Antarctica and covers 20 square degree of sky around the South Celestial Pole. The installation is designed to provide high-cadence photometry for the purpose of monitoring the quality of the astronomical observing conditions at Dome A and detecting transiting exoplanets. CSTAR has been operational since 2008, and has taken a rich and high-precision photometric data set of 10,690 stars. In the first observing season, we obtained 291,911 qualified science frames with 20-second integrations in the i-band. Photometric precision reaches about 4 mmag at 20-second cadence at i=7.5, and is about 20 mmag at i=12. Using robust detection methods, ten promising exoplanet candidates were found. Four of these were found to be giants using spectroscopic follow-up. All of these transit candidates are presented here along with the discussion of their detailed properties as well as the follow-up observations.

The photometric characterization of M33 star clusters is far from complete. In this paper, we present homogeneous $UBVRI$ photometry of 708 star clusters and cluster candidates in M33 based on archival images from the Local Group Galaxies Survey, which covers 0.8 deg$^2$ along the galaxy's major axis. Our photometry includes 387, 563, 616, 580, and 478 objects in the $UBVRI$ bands, respectively, of which 276, 405, 430, 457, and 363 do not have previously published $UBVRI$ photometry. Our photometry is consistent with previous measurements (where available) in all filters. We adopted Sloan Digital Sky Survey $ugriz$ photometry for complementary purposes, as well as Two Micron All-Sky Survey near-infrared $JHK$ photometry where available. We fitted the spectral-energy distributions of 671 star clusters and candidates to derive their ages, metallicities, and masses based on the updated {\sc parsec} simple stellar populations synthesis models. The results of our $\chi^2$ minimization routines show that only 205 of the 671 clusters ($31\%$) are older than 2 Gyr, which represents a much smaller fraction of the cluster population than that in M31 ($56\%$), suggesting that M33 is dominated by young star clusters ($<1$ Gyr). We investigate the mass distributions of the star clusters---both open and globular clusters---in M33, M31, the Milky Way, and the Large Magellanic Cloud. Their mean values are $\log(M_{\rm cl}/M_{\odot})=4.25$, 5.43, 2.72, and 4.18, respectively. The fraction of open to globular clusters is highest in the Milky Way and lowest in M31. Our comparisons of the cluster ages, masses, and metallicities show that our results are basically in agreement with previous studies (where objects in common are available); differences can be traced back to differences in the models adopted, the fitting methods used, and stochastic sampling effects.

In this paper, we analyze in detail with numerical simulations how the mask effect can influence the weak lensing peak statistics reconstructed from the shear measurement of background galaxies. It is found that high peak fractions are systematically enhanced due to masks, the larger the masked area, the higher the enhancement. In the case with about $13\%$ of the total masked area, the fraction of peaks with SNR $\nu\ge 3$ is $\sim 11\%$ in comparison with $\sim 7\%$ of the mask-free case in our considered cosmological model. This can induce a large bias on cosmological studies with weak lensing peak statistics. Even for a survey area of $9\hbox{ deg}^2$, the bias in $(\Omega_m, \sigma_8)$ is already close to $3\sigma$. It is noted that most of the affected peaks are close to the masked regions. Therefore excluding peaks in those regions can reduce the bias but at the expense of loosing usable survey areas. Further investigations find that the enhancement of high peaks number can be largely attributed to higher noise led by the fewer number of galaxies usable in the reconstruction. Based on Fan et al. (2010), we develop a model in which we exclude only those large masks with radius larger than $3\arcmin. For the remained part, we treat the areas close to and away from the masked regions separately with different noise levels. It is shown that this two-noise-level model can account for the mask effect on peak statistics very well and the cosmological bias is significantly reduced.

In this paper we investigate the nature and distribution of large neutral regions during the late epoch of reionization. In the "bubble model" of reionization, the mass distribution of large ionized regions ("bubbles") during the early stage of reionization is obtained by using the excursion set model, where the ionization of a region corresponds to the first up-crossing of a barrier by random trajectories. We generalize this idea, and develop a method to predict the distribution of large scale neutral regions during the late stage of reionization, taking into account the ionizing background after the percolation of HII regions. The large scale neutral regions which we call "neutral islands" are not individual galaxies or minihalos, but larger regions where fewer galaxies formed and hence ionized later, and they are identified in the excursion set model with the first down-crossings of the island barrier. Assuming that the consumption rate of ionizing background photons is proportional to the surface area of the neutral islands, we obtained the size distribution of the neutral islands. We also take the "bubbles-in-island" effect into account by considering the conditional probability of up-crossing a bubble barrier after down-crossing the island barrier. We find that this effect is very important. An additional barrier is set to avoid islands being percolated through. We find that there is a characteristic scale for the neutral islands, while the small islands are rapidly swallowed up by the ionizing background, this characteristic scale does not change much as the reionization proceeds.

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

The quasi-simultaneous multi-wavelength emission of TeV blazar 3C 66A is studied by using a one-zone multi-component leptonic jet model. It is found that the quasi-simultaneous spectral energy distribution (SED) of 3C 66A can be well reproduced, especially its {\it Fermi}-LAT first 3 months average spectrum can be well reproduced by the synchrotron-self Compton (SSC) component plus external Compton (EC) component of the broad line region (BLR). Clues on its redshift and gamma-ray emission location are obtained. The results indicate the following. (i) On the redshift; The theoretical intrinsic TeV spectra can be predicted by extrapolating the reproduced GeV spectra. Through comparing this extrapolated TeV spectra with the extragalactic background light (EBL) corrected observed TeV spectra, it is suggested that the redshift of 3C 66A could be between 0.1 and 0.3, the most likely value is $\sim$ 0.2. (ii) On the gamma-ray emission location; To well reproduce the GeV emission of 3C 66A under different assumptions on BLR, the gamma-ray emission region is always required to be beyond the inner zone of BLR. The BLR absorption effect on gamma-ray emission confirms this point.

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

In this paper, we present photometry for young star clusters in M31, which are selected from Caldwell et al. These star clusters have been observed as part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey from 1995 February to 2008 March. The BATC images including these star clusters are taken with 15 intermediate-band filters covering 3000--10000 \AA. Combined with photometry in the {\sl GALEX} far- and near-ultraviolet, broad-band $UBVRI$, SDSS $ugriz$, and infrared $JHK_{\rm s}$ of Two Micron All Sky Survey, we obtain their accurate spectral energy distributions (SEDs) from 1538-20000 \AA. We derive these star clusters' ages and masses by comparing their SEDs with stellar population synthesis models. Our results are in good agreement with previous determinations. The mean value of age and mass of young clusters ($<2$ Gyr) is about 385 Myr and $2\times 10^4 {M_\odot}$, respectively. There are two distinct peaks in the age distribution, a highest peak at age $\sim$ 60 Myr and a secondary peak around 250 Myr, while the mass distribution shows a single peak around $10^4 {M_\odot}$. A few young star clusters have two-body relaxation times greater than their ages, indicating that those clusters have not been well dynamically relaxed and therefore have not established the thermal equilibrium. There are several regions showing aggregations of young star clusters around the 10 kpc ring and the outer ring, indicating that the distribution of the young star clusters is well correlated with M31's star-forming regions. The young massive star clusters (age $\leq 100$ Myr and mass $\geq 10^4 {M_\odot}$) show apparent concentration around the ring splitting region, suggesting a recent passage of a satellite galaxy (M32) through M31 disk.

We present extensive optical observations of a Type IIn supernova (SN) 2010jl for the first 1.5 years after the discovery. The UBVRI light curves demonstrated an interesting two-stage evolution during the nebular phase, which almost flatten out after about 90 days from the optical maximum. SN 2010jl has one of the highest intrinsic H_alpha luminosity ever recorded for a SN IIn, especially at late phase, suggesting a strong interaction of SN ejecta with the dense circumstellar material (CSM) ejected by the progenitor. This is also indicated by the remarkably strong Balmer lines persisting in the optical spectra. One interesting spectral evolution about SN 2010jl is the appearance of asymmetry of the Balmer lines. These lines can be well decomposed into a narrow component and an intermediate-width component. The intermediate-width component showed a steady increase in both strength and blueshift with time until t ~ 400 days after maximum, but it became less blueshifted at t ~ 500 days when the line profile appeared relatively symmetric again. Owing to that a pure reddening effect will lead to a sudden decline of the light curves and a progressive blueshift of the spectral lines, we therefore propose that the asymmetric profiles of H lines seen in SN 2010jl is unlikely due to the extinction by newly formed dust inside the ejecta, contrary to the explanation by some early studies. Based on a simple CSM-interaction model, we speculate that the progenitor of SN 2010jl may suffer a gigantic mass loss (~ 30-50 M_sun) in a few decades before explosion. Considering a slow moving stellar wind (e.g., ~ 28 km/s) inferred for the preexisting, dense CSM shell and the extremely high mass-loss rate (1-2 M_sun per yr), we suggest that the progenitor of SN 2010jl might have experienced a red supergiant stage and explode finally as a post-red supergiant star with an initial mass above 30-40 M_sun.

Simple stellar population (SSP) synthesis models are useful tools for studying the nature of unresolved star clusters in external galaxies. However, the plethora of currently available SSP models gives rise to significant and poorly documented systematic differences. Here we consider the outputs of the commonly used Bruzual & Charlot and GALEV models, as well as a recently updated SSP model suite which attempts to include the contributions of binary merger products in the form of blue straggler stars (BS-SSP). We rederive the ages, metallicities, extinction values and masses of 445 previously observed globular-like clusters in M31 based on chi-square minimisation of their spectral energy distributions with respect to these three different SSP models and adopting a Chabrier-like stellar initial mass function. A comparison between our new results and previous estimates of the same parameters shows that the Bruzual & Charlot models yield the youngest ages and lowest masses, while adoption of the BS-SSP models results in the oldest ages and highest mass estimates. Similarly, the GALEV SSP models produce the lowest metallicities, with the highest values resulting from the BS-SSP model suite. These trends are caused by intrinsic differences associated with the models, and are not significantly affected by the well-known age-metallicity degeneracy. Finally, we note that the mass function of the massive M31 star clusters is similar to that of the Milky Way's globular clusters, which implies that the two star cluster systems likely formed under similar environmental conditions.

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

Stochastic acceleration of charged particles due to their interactions with plasma waves may be responsible for producing superthermal particles in a variety of astrophysical systems. This process can be described as a diffusion process in the energy space with the Fokker-Planck equation. In this paper, a time-dependent numerical code is used to solve the reduced Fokker-Planck equation involving only time and energy variables with general forms of the diffusion coefficients. We also propose a self-similar model for particle acceleration in Sedov explosions and use the TeV SNR RX J1713.7-3946 as an example to demonstrate the model characteristics. Markov Chain Monte Carlo method is utilized to constrain model parameters with observations.

(Abridged) We performed the spectroscopic observations of 11 confirmed GCs in M31 with the Xinglong 2.16m telescope and we mainly focus on the fits method and the metallicity gradient for the M31 GC sample. We analyzed and discussed more about the dynamics, metallicity and age, and their distributions as well as the relationships between these parameters. Eight more confirmed GCs in the halo of M31 were observed, most of which lack the spectroscopic information before. These star clusters are located far from the galactic center at a projected radius of ~14 to ~117 kpc. The Lick absorption-line indices and the radial velocities have been measured and ages, metallicities [Fe/H] and alpha-element [alpha/Fe] have also been fitted by comparing the observed spectral feature indices and Thomas et al.SSP model. Our results show that most of the star clusters of our sample are older than 10 Gyr except B290 ~5.5 Gyr, and most of them are metal-poor with the metallicity [Fe/H]<-1, suggesting that these clusters were born at the early stage of the galaxy's formation. We find that the metallicity gradient for the outer halo clusters with r_p>25 kpc may not exist with a slope of -0.005+-0.005 dex kpc^-1. We also find that the metallicity is not a function of age for the GCs with age < 7 Gyr while for the old GCs with age >7 Gyr there seems to be a trend that the older ones have lower metallicity. Besides, We plot metallicity distributions with the largest sample of M31 GCs so far and it shows the bimodality is not significant and the number of the metal-poor and metal-rich groups becomes comparable. The spatial distributions shows that the metal-rich group is more centrally concentrated while the metal-poor group is occupy a more extended halo and the young population is centrally concentrated while the old population is more extended spatially to the outer halo.

Using volume-limited samples drawn from The Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we measure the tidal environment of galaxies, which we characterize by the ellipticity e of the potential field calculated from the smoothed spatial number density 1 + {\delta} of galaxies. We analyze if galaxy properties, including color, Dn4000, concentration and size correlate with e, in addition to depending on 1 + {\delta}. We find that there exists a transition smoothing scale at which correlations/anti-correlations with e reverse. This transition scale is well represented by the distance to the 3rd- nearest-neighbor of a galaxy in a volume limited sample with Mr < -20 which has a distribution peaked at ~ 2 h-1Mpc. We further demonstrate that this scale corresponds to that where the correlation between the color of galaxies and environmental density 1+{\delta} is the strongest. For this optimal smoothing R0 no additional correlations with e are observed. The apparent dependence on tidal ellipticity e at other smoothing scales Rs can be viewed as a geometric effect, arising from the cross correlation be- tween (1+{\delta}o) and e(Rs). We perform the same analysis on numerical simulations with semi-analytical modeling (SAM) of galaxy formation. The e dependence of the galaxy properties shows similar behavior to that in the SDSS, although the color-density cor- relation is significantly stronger in the SAM. The 'optimal adaptive smoothing scale' in the SAM is also closely related to the distance to the 3rd-nearest-neighbor of a galaxy, and its characteristic value is consistent with, albeit slightly smaller than for SDSS.

We present the first weak gravitational lensing analysis of the completed Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We study the 64 square degrees W1 field, the largest of the CFHTLS-Wide survey fields, and present the largest contiguous weak lensing convergence "mass map" yet made. 2.66 million galaxy shapes are measured, using a KSB pipeline verified against high-resolution Hubble Space Telescope imaging that covers part of the CFHTLS. Our i'-band measurements are also consistent with an analysis of independent r'-band imaging. The reconstructed lensing convergence map contains 301 peaks with signal-to-noise ratio {\nu}>3.5, consistent with predictions of a {\Lambda}CDM model. Of these peaks, 126 lie within 3.0' of a BCG identified from multicolor optical imaging in an independent, red sequence survey. We also identify 7 counterparts for massive clusters previously seen in X-ray emission within 6 square degrees XMM-LSS survey. With photometric redshift estimates for the source galaxies, we use a tomographic lensing method to fit the redshift and mass of each convergence peak. Matching these to the optical observations, we confirm 85 groups/clusters with \chi^2 reduced < 3.0, at a mean redshift <z_c> = 0.36 and velocity dispersion <\sigma_c> = 658.8 km/s. Future surveys, such as DES, LSST, KDUST and EUCLID, will be able to apply these same techniques to map clusters in much larger volumes and thus tightly constrain cosmological models.

In this paper, we present a simulation method within the two-component spherical collapse model to investigate dark energy perturbations associated with the formation of dark matter halos. The realistic mass accretion history of a dark matter halo taking into account its fast and slow growth is considered by imposing suitable initial conditions and isotropized virializations for the spherical collapse process. The dark energy component is treated as a perfect fluid described by two important parameters, the equation of state parameter $w$ and the sound speed $c_s$. Quintessence models with $w>-1$ are analyzed. We adopt the Newtonian gauge to describe the spacetime which is perturbed mainly by the formation of a dark matter halo. It is found that the dark energy density perturbation $\delta_{DE}$ depends on $w$ and $c_s$, and its behavior follows closely the gravitational potential $\Phi$ of the dark matter halo with $\delta_{DE}\approx -(1+w)\Phi/c_s^2$. For $w>-1$, the dark energy perturbation presents a clustering behavior with $\delta_{DE}>0$ during the entire formation of the dark matter halo, from linear to nonlinear and virialized stages. The value of $\delta_{DE}$ increases with the increase of the halo mass. For a cluster of mass $M\sim 10^{15} M_{\odot}$, $\delta_{DE}\sim 10^{-5}$ within the virialized region for $c_s^2 \in [0.5, 1]$, and it can reach $\delta_{DE}=O(1)$ with $c_s^2=0.00001$. For a scalar-field dark energy model, we find that with suitably modeled $w$ and $c_s$, its perturbation behavior associated with the nonlinear formation of dark matter halos can well be analyzed using the fluid approach, demonstrating the validity of the fluid description for dark energy even considering its perturbation in the stage of nonlinear dark matter structure formation.

B514 is a remote M31 globular cluster which locating at a projected distance of R_p~55 kpc. Deep observations with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) are used to provide the accurate integrated light and star counts of B514. By coupling analysis of the distribution of the integrated light with star counts, we are able to reliably follow the profile of the cluster out to ~40". Based on the combined profile, we study in detail its surface brightness distribution in F606W and F814W filters, and determine its structural parameters by fitting a single-mass isotropic King model. The results showed that, the surface brightness distribution departs from the best-fit King model for r>10". B514 is quite flatted in the inner region, and has a larger half-light radius than majority of normal globular clusters of the same luminosity. It is interesting that, in the M_V versus log R_h plane, B514 lies nearly on the threshold for ordinary globular clusters as defined by Mackey & van den Bergh. In addition, B514 was observed as part of the Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey, using 13 intermediate-band filters covering a wavelength range of 3000--8500 \AA. Based on aperture photometry, we obtain its SEDs as defined by the 13 BATC filters. We determine the cluster's age and mass by comparing its SEDs (from 2267 to 20000{\AA}, comprising photometric data in the near-ultraviolet of GALEX, 5 SDSS bands, 13 BATC intermediate-band, and 2MASS near-infrared JHKs} filters) with theoretical stellar population synthesis models, resulting in age of $11.5\pm3.5$ Gyr. This age confirms the previous suggestion that B514 is an old GC in M31. B514 has a mass of $0.96-1.08 \times 10^6 \rm M_sun$, and is a medium-mass globular cluster in M31.

We present the spectroscopic observations for 11 confirmed globular clusters of M31 with the OMR spectrograph on 2.16m telescope at Xinglong site of National Astronomical Observatories, Chinese Academy of Sciences. Nine of our sample clusters are located in the halo of M31 and the most remote one is out to a projected radius of 78.75 kpc from the galactic center. For all our sample clusters, we measured the Lick absorption-line indices and the radial velocities. It is noted that most GCs of our sample are distinct from the HI rotation curve of M31 galaxy, especially for B514, MCGC5, H12 and B517, suggesting that most of our sample clusters do not have kinematic association with the star forming young disk of the galaxy. We fitted the absorption line indices with the updated stellar population model Thomas et al. (2010) with two different tracks of Cassisi and Padova, separately, by applying the $\chi^2-$minimization method. The fitting results show that all our sample clusters are older than 10 Gyr, and metal-poor ($-0.91 \le $ [Fe/H] $\le -2.38$ dex). After merging the spectroscopic metallicity of our work with the previously published ones, we extended the cluster sample out to a projected radius of 117 kpc from the galaxy's center. We found the metallicity gradient for all the confirmed clusters exists with a slope of $-0.028\pm0.001$ dex kpc$^{-1}$. However, the slope turns to be $-0.018\pm0.001$ dex kpc$^{-1}$ for all the halo clusters, which is much shallower. If we only consider the outer halo clusters with $r_{\rm p}>25$ kpc, the slope becomes $-0.010\pm0.002$ dex kpc$^{-1}$ and if one cluster G001 is excluded from the outer halo sample, the slope is $-0.004\pm0.002$ dex kpc$^{-1}$. Thus we conclude that metallicity gradient for M31 outer halo clusters is not significant, which agrees well with the previous findings.

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

Panchromatic spectral energy distribution (SED) from the ultraviolet (UV), optical to infrared (IR) photometry of NGC 628, combined with the evolutionary stellar population synthesis, is used to derive the spatially resolved age, metallicity and reddening maps. These parameter distributions show that the bulge of this galaxy is a disk-like pseudobulge, which has the S{\'e}rsic index close to the exponential law, rich gas, and a young circumnuclear ring structure. We also discover the disk has two distinct regions with different radial age and metallicity gradients. The inner region is older and has a much steeper age gradient than the outer region of the disk. Both these two regions and the central young structure can be seen in the radial profile of the optical color. Based on the age and reddening distributions, we consider that the pseudobulge and disk are likely to have grown via the secular evolution, which is the redistribution of mass and energy through the angular momentum transport caused by the non-axisymmetric potential of the spirals. However, possible gas accretion events could affect the outer region of the disk, due to abundant H{\sc i} gas accumulating in the outer disk.

We study large-scale structures from numerical simulations, paying particular attention to supercluster-like structures. A grid-density-contour based algorithm is adopted to locate connected groups. With the increase of the linking density threshold, the foam- like cosmic web is subsequently broken into individual supercluster-like groups and further halos which are in accordance to groups with the linking density threshold {\rho}/{\rho}= 1 + {\delta} = 80. By analyzing sets of cosmological simulations with varying cosmological parameters, we find that an universal mass function exists not only for halos but also for low-density supercluster-like groups until the linking density threshold decreases to a density where the global percolation of large-scale structures occurs. We further show that the mass functions of different groups can be well described by the Jenkins form with the parameters being dependent on the linking density threshold. On the other hand, these low- density supercluster-like groups cannot be directly associated with the predictions from the excursion set theory with effective barriers obtained from dynamical collapse models, and the peak exclusion effect must be taken into account. Including such an effect, the mass function of groups with the linking density threshold 1 + {\delta} = 16 is in good agreements with that from the excursion set theory with a nearly flat effective barrier. A simplified analysis of the ellipsoidal collapse model indicates that the barrier for collapses along two axes to form filaments is approximately flat in scales, thus we define groups identified with 1 + {\delta} = 16 as filaments. We then further study the halo-filament conditional mass function and the filament-halo conditional mass function, and compare them with the predictions from the two-barrier excursion set theory. The shape statistics for filaments are also presented.

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.

VDB0-B195D is a massive, blue star cluster in M31. It was observed as part of the Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey using 15 intermediate-band filters covering a wavelength range of 3000--10,000 \AA. Based on aperture photometry, we obtain its spectral-energy distribution (SED) as defined by the 15 BATC filters. We apply previously established relations between the BATC intermediate-band and the Johnson-Cousins $UBVRI$ broad-band systems to convert our BATC photometry to the standard system. A detailed comparison shows that our newly derived $VRI$ magnitudes are fully consistent with previous results, while our new $B$ magnitude agrees to within $2\sigma$. In addition, we determine the cluster's age and mass by comparing its SED (from 3000 to 20,000{\AA}, comprising photometric data in the 15 BATC intermediate bands, optical broad-band $BVRI$, and 2MASS near-infrared $JHK_{\rm s}$ data) with theoretical stellar population synthesis models, resulting in age and mass determinations of $60.0\pm 8.0$~Myr and $(1.1-1.6) \times 10^5 M_\odot$, respectively. This age and mass confirms previous suggestions that VDB0-B195D is a young massive cluster in M31.

Taking into account the noise from intrinsic ellipticities of source galaxies, we study the efficiency and completeness of halo detections from weak lensing convergence maps. Particularly, with numerical simulations, we compare the Gaussian filter with the so called MRLens treatment based on the modification of the Maximum Entropy Method. For a pure noise field without lensing signals, a Gaussian smoothing results a residual noise field that is approximately Gaussian in statistics if a large enough number of galaxies are included in the smoothing window. On the other hand, the noise field after the MRLens treatment is significantly non-Gaussian, resulting complications in characterizing the noise effects. Considering weak-lensing cluster detections, although the MRLens treatment effectively deletes false peaks arising from noise, it removes the real peaks heavily due to its inability to distinguish real signals with relatively low amplitudes from noise in its restoration process. The higher the noise level is, the larger the removal effects are for the real peaks. For a survey with a source density n_g~30 arcmin^(2), the number of peaks found in an area of 3x3 sq.deg after MRLens filtering is only ~50 for the detection threshold kappa=0.02, while the number of halos with M>5x10^{13} M_{\odot} and with redshift z<=2 in the same area is expected to be ~530. For the Gaussian smoothing treatment, the number of detections is ~260, much larger than that of the MRLens. The Gaussianity of the noise statistics in the Gaussian smoothing case adds further advantages for this method to circumvent the problem of the relatively low efficiency in weak-lensing cluster detections. Therefore, in studies aiming to construct large cluster samples from weak-lensing surveys, the Gaussian smoothing method performs significantly better than the MRLens.

We present an updated UBVRI photometric catalog containing 970 objects in the field of M31, selected from the Revised Bologna Catalog (RBC v.4.0), including 965, 967, 965, 953, and 827 sources in the individual UBVRI bands, respectively, of which 205, 123, 14, 126, and 109 objects do not have previously published photometry. Photometry is performed using archival images from the Local Group Galaxies Survey, which covers 2.2 deg^2 along the major axis of M31. We focus on 445 confirmed `globular-like' clusters and candidates, comprising typical globular and young massive clusters. The ages and masses of these objects are derived by comparison of their observed spectral-energy distributions with simple stellar population synthesis. Approximately half of the clusters are younger than 2 Gyr, suggesting that there has been significant recent active star formation in M31, which is consistent with previous results. We note that clusters in the halo (r_ projected>30kpc) are composed of two different components, older clusters with ages >10 Gyr and younger clusters with ages around 1 Gyr. The spatial distributions show that the young clusters (<2 Gyr) are spatially coincident with the galaxy's disk, including the `10 kpc ring,' the `outer ring,' and the halo of M31, while the old clusters (> 2 Gyr) are spatially correlated with the bulge and halo. We also estimate the masses of the 445 confirmed clusters and candidates in M31 and find that our estimates agree well with previously published values. We find that none of the young disk clusters can survive the inevitable encounters with giant molecular clouds in the galaxy's disk and that they will eventually disrupt on timescales of a few Gyr. Specifically, young disk clusters with a mass of 10^4 M_\odot are expected to dissolve within 3.0 Gyr and will, thus, not evolve to become globular clusters.

M31 globular cluster B379 is the first extragalactic cluster, the age of which was determined by main-sequence photometry. In this method, the age of a cluster is obtained by fitting its CMD with stellar evolutionary models. However, different stellar evolutionary models use different parameters of stellar evolution, such as range of stellar masses, different opacities and equations of state, and different recipes, and so on. So, it is interesting to check whether different stellar evolutionary models can give consistent results for the same cluster. Brown et al. (2004a) constrained the age of B379 by comparing its CMD with isochrones of the 2006 VandenBerg models. Using SSP models of BC03 and its multi-photometry, Ma et al. (2007) independently determined the age of B379, which is in good agreement with the determination of Brown et al. (2004a). The BC03 models are calculated based on the Padova evolutionary tracks. It is necessary to check whether the age of B379 which, being determined based on the Padova evolutionary tracks, is in agreement with the determination of Brown et al. (2004a). So, in this paper, we re-determine its age using isochrones of the Padova stellar evolutionary models. In addition, the metal abundance, the distance modulus, and the reddening value for B379 are also determined in this paper. The results obtained in this paper are consistent with the previous determinations, which including the age obtained by Brown et al. (2004a). So, this paper confirms the consistence of the age scale of B379 between the Padova isochrones and the 2006 VandenBerg isochrones, i.e. the results' comparison between Brown et al. (2004a) and Ma et al. (2007) is meaningful. The results obtained in this paper are: the metallicity [M/H]=-0.325, the age $\tau=11.0\pm1.5$ Gyr, the reddening value E(B-V)=0.08, and the distance modulus $(m-M)_{0}=24.44\pm0.10$.

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

Taking into account noise from intrinsic ellipticities of source galaxies, in this paper, we study the peak statistics in weak-lensing convergence maps around clusters of galaxies and beyond. We emphasize how the noise peak statistics is affected by the density distribution of nearby clusters, and also how cluster-peak signals are changed by the existence of noise. These are the important aspects to be understood thoroughly in weak-lensing analyses for individual clusters as well as in cosmological applications of weak-lensing cluster statistics. We adopt Gaussian smoothing with the smoothing scale $\theta_G=0.5\hbox{ arcmin}$ in our analyses. It is found that the noise peak distribution near a cluster of galaxies depends sensitively on the density profile of the cluster. For a cored isothermal cluster with the core radius $R_c$, the inner region with $R\le R_c$ appears noisy containing on average $\sim 2.4$ peaks with $\nu\ge 5$ for $R_c= 1.7\hbox{ arcmin}$ and the true peak height of the cluster $\nu=5.6$, where $\nu$ denotes the convergence signal to noise ratio. For a NFW cluster of the same mass and the same central $\nu$, the average number of peaks with $\nu\ge 5$ within $R\le R_c$ is $\sim 1.6$. Thus a high peak corresponding to the main cluster can be identified more cleanly in the NFW case. In the outer region with $R_c<R\le 5R_c$, the number of high noise peaks is considerably enhanced in comparison with that of the pure noise case without the nearby cluster. (abridged)

In the leptonic scenario for TeV emission from a few well-observed shell-type TeV supernova remnants (STTSNRs), very weak magnetic fields are inferred. If fast-mode waves are produced efficiently in the shock downstream, we show that they are viable agents for acceleration of relativistic electrons inferred from the observed spectra even in the subsonic phase, in spite that these waves are subject to strong damping by thermal background ions at small dissipation scales. Strong collisionless non-relativistic astrophysical shocks are studied with the assumption of a constant Aflven speed in the downstream. The turbulence evolution is modeled with both the Kolmogorov and Kraichnan phenomenology. Processes determining the high-energy cutoff of nonthermal electron distributions are examined. The Kraichnan models lead to a shallower high-energy cutoff of the electron distribution and require a lower downstream density than the Kolmogorov models to fit a given emission spectrum. With reasonable parameters, the model explains observations of STTSNRs, including recent data obtained with the Fermi gamma-ray telescope. More detailed studies of the turbulence generation and dissipation processes, supernova explosions and progenitors are warranted for better understanding the nature of supernova shocks.

Dark energy is an important science driver of many upcoming large-scale surveys. With small, stable seeing and low thermal infrared background, Dome A, Antarctica, offers a unique opportunity for shedding light on fundamental questions about the universe. We show that a deep, high-resolution imaging survey of 10,000 square degrees in \emph{ugrizyJH} bands can provide competitive constraints on dark energy equation of state parameters using type Ia supernovae, baryon acoustic oscillations, and weak lensing techniques. Such a survey may be partially achieved with a coordinated effort of the Kunlun Dark Universe Survey Telescope (KDUST) in \emph{yJH} bands over 5000--10,000 deg$^2$ and the Large Synoptic Survey Telescope in \emph{ugrizy} bands over the same area. Moreover, the joint survey can take advantage of the high-resolution imaging at Dome A to further tighten the constraints on dark energy and to measure dark matter properties with strong lensing as well as galaxy--galaxy weak lensing.

We present the correlations between the spectroscopic metallicities and ninety-three different intrinsic colors of M31 globular clusters, including seventy-eight BATC colors and fifteen SDSS and near infrared ugrizK colors. The BATC colors were derived from the archival images of thirteen filters (from c to p), which were taken by Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey with a 60/90 cm f/3 Schmidt telescope. The spectroscopic metallicities adopted in our work were from literature. We fitted the correlations of seventy-eight different BATC colors and the metallicities for 123 old confirmed globular clusters, and the result implies that correlation coefficients of twenty-three colors r>0.7. Especially, for the colors $(f-k)_0$, $(f-o)_0$, and $(h-k)_0$, the correlation coefficients are r>0.8. Meanwhile, we also note that the correlation coefficients (r) approach zero for $(g-h)_0$, $(k-m)_0$, $(k-n)_0$, and $(m-n)_0$, which are likely to be independent of metallicity. Similarity, we fitted the correlations of metallicity and ugrizK colors for 127 old confirmed GCs. The result indicates that all these colors are metal-sensitive (r>0.7), of which $(u-K)_0$ is the most metal-sensitive color. Our work provides an easy way to simply estimate the metallicity from colors.

In 2008 January the 24th Chinese expedition team successfully deployed the Chinese Small Telescope ARray (CSTAR) to DomeA, the highest point on the Antarctic plateau. CSTAR consists of four 14.5cm optical telescopes, each with a different filter (g, r, i and open) and has a 4.5degree x 4.5degree field of view (FOV). It operates robotically as part of the Plateau Observatory, PLATO, with each telescope taking an image every 30 seconds throughout the year whenever it is dark. During 2008, CSTAR #1 performed almost flawlessly, acquiring more than 0.3 million i-band images for a total integration time of 1728 hours during 158 days of observations. For each image taken under good sky conditions, more than 10,000 sources down to 16 mag could be detected. We performed aperture photometry on all the sources in the field to create the catalog described herein. Since CSTAR has a fixed pointing centered on the South Celestial Pole (Dec =-90 degree), all the sources within the FOV of CSTAR were monitored continuously for several months. The photometric catalog can be used for studying any variability in these sources, and for the discovery of transient sources such as supernovae, gamma-ray bursts and minor planets.

We present photometry of 104 M31 globular clusters (GCs) and GC candidates in 15 intermediate-band filters of the Beijing-Arizona-Taiwan-Connecticut (BATC) photometric system. The GCs and GC candidates were selected from the Revised Bologna Catalog (v.3.5). We obtain the cluster ages by comparing the photometric data with up-to-date theoretical synthesis models. The photometric data used are {\sl GALEX} far- and near-ultraviolet and 2MASS near-infrared $JHK_{\rm s}$ magnitudes, combined with optical photometry. The ages of our sample clusters cover a large range, although most clusters are younger than 10 Gyr. Combined with the ages obtained in our series of previous papers focusing on the M31 GC system, we present the full M31 GC age distribution. The M31 GC system contains populations of young and intermediate-age GCs, as well as the `usual' complement of well-known old GCs, i.e., GCs of similar age as the majority of the Galactic GCs. In addition, young GCs (and GC candidates) are distributed nearly uniformly in radial distance from the center of M31, while most old GCs (and GC candidates) are more strongly concentrated.

Using high redshift radio sources as background, the 21cm forest observations probe the neutral hydrogen absorption signatures of early structures along the lines of sight. Directly sensitive to the spin temperature, it complements the 21cm tomography observations, and provides information on the temperature as well as the ionization state of the intergalactic medium (IGM). We use a radiative transfer simulation to investigate the 21cm forest signals during the epoch of reionization. We first confirmed that the optical depth and equivalent width (EW) are good representations of the ionization and thermal state of the IGM. The features selected by their relative optical depth are excellent tracers of ionization fields, while those selected by absolute optical depth are very sensitive to the IGM temperature, so the IGM temperature information could potentially be extracted, thus breaking a degeneracy in 21cm tomographic observation. From the distributions of EWs and the number evolution of absorbers and leakers with different EWs, we see clearly the cosmological evolution of ionization state of the IGM. The number density of potentially observable features decreases rapidly with increasing gas temperature. The sensitivity of the proposed EW statistic to the IGM temperature makes it a unique and potentially powerful probe of reionization. Missing small-scale structures, such as small filaments and minihalos that are unresolved in our current simulation, and lack of an accurate calculation of the IGM temperature, however, likely have rendered the presented signals quantitatively inaccurate. Finally we discuss the requirements of the background radio sources for such observations, and find that signals with equivalent widths larger than 1kHz are hopeful to be detected.(Abridged)

We have collected all available spectra and photometric data from SDSS catalogue for bright AGNs complied from the first three months of the Fermi large area telescope all-sky survey. Based on the 106 high-confidence and 11 low-confidence association bright AGN list, the photometry data are collected from SDSS DR7 for 28 sources (12 BL Lacs and 16 FSRQs), two of which are low-confidence association bright AGNs. Among these 28 SDSS photometric sources, SDSS spectra are available for 20 sources (6 BL Lacs and 14 FSRQs). The black hole mass M_BH and the broad line region (BLR) luminosity were obtained for 14 FSRQs by measuring the line-width and strength of broad emission lines from SDSS spectra. The broad emission lines measurements of five FSRQs are presented for the first time in this work. The optical continuum emission of these 14 FSRQs is found to be likely dominated by the nonthermal jet emission through comparing the relationship between the broad Mg II line and continuum luminosity to that of radio quiet AGNs. The black hole mass of 14 FSRQs ranges from 10^8.2 to 10^9.9 solar mass, with most of sources larger than 10^9 solar mass. The Eddington ratio L_bol/L_Edd ranges from 10^-1.5 to ~ 1. This implies that the optically thin, geometrically thick accretion disk may exist in these FSRQs.

In this paper, we analyze the dynamical evolution of quintessence dark energy induced by the collapse of dark matter halos. Different from other previous studies, we develop a numerical strategy which allows us to calculate the dark energy evolution for the entire history of the spherical collapse of dark matter halos, without the need of separate treatments for linear, quasi-linear and nonlinear stages of the halo formation. It is found that the dark energy perturbations evolve with redshifts, and their specific behaviors depend on the quintessence potential as well as the collapsing process. The overall energy density perturbation is at the level of $10^{-6}$ for cluster-sized halos. The perturbation amplitude decreases with the decrease of the halo mass. At a given redshift, the dark energy perturbation changes with the radius to the halo center, and can be either positive or negative depending on the contrast of $\partial_t \phi$, $\partial_r \phi$ and $\phi$ with respect to the background, where $\phi$ is the quintessence field. For shells where the contrast of $\partial_r \phi$ is dominant, the dark energy perturbation is positive and can be as high as about $10^{-5}$.

We present the photometry of 30 globular clusters (GCs) and GC candidates in 15 intermediate-band filters covering from ~3000 to ~10000 \AA using the archival CCD images of M31 observed as part of the Beijing - Arizona - Taiwan - Connecticut (BATC) Multicolour Sky Survey. We transform these intermediate-band photometric data to the photometry in the standard UBVRI broad-bands. These M31 GC candidates are selected from the Revised Bologna Catalogue (RBC V.3.5), and most of these candidates do not have any photometric data. Therefore the present photometric data are supplement to RBC V.3.5. We find that 4 out of 61 GCs and GC candidates in RBC V.3.5 do not show any signal on the BATC images at their locations. By linear fit of the distribution in colour-magnitude diagram of blue GCs and GC candidates using the data from RBC V.3.5, and in this study we find the ``blue-tilt'' of blue M31 GCs with a high confidence at 99.95% or 3.47 sigma for the confirmed GCs, and >99.99% or 4.87 sigma for GCs and GC candidates.

We study the impact of catastrophic errors occurring in the photometric redshifts of galaxies on cosmological parameter estimates with cosmic shear tomography. We consider a fiducial survey with 9-filter set and perform photo-z measurement simulations. It is found that a fraction of 1% galaxies at z_{spec}~0.4 is misidentified to be at z_{phot}~3.5. We then employ both chi^2 fitting method and the extension of Fisher matrix formalism to evaluate the bias on the equation of state parameters of dark energy, w_0 and w_a, induced by those catastrophic outliers. By comparing the results from both methods, we verify that the estimation of w_0 and w_a from the fiducial 5-bin tomographic analyses can be significantly biased. To minimize the impact of this bias, two strategies can be followed: (A) the cosmic shear analysis is restricted to 0.5<z<2.5 where catastrophic redshift errors are expected to be insignificant; (B) a spectroscopic survey is conducted for galaxies with 3<z_{phot}<4. We find that the number of spectroscopic redshifts needed scales as N_{spec} \propto f_{cata}\times A where f_{cata}=1% is the fraction of catastrophic redshift errors (assuming a 9-filter photometric survey) and A is the survey area. For A=1000 {deg}^2, we find that N_{spec}>320 and 860 respectively in order to reduce the joint bias in (w_0,w_a) to be smaller than 2\sigma and 1\sigma. This spectroscopic survey (option B) will improve the Figure of Merit of option A by a factor \times 1.5 thus making such a survey strongly desirable.

Mayall II = G1 is one of the most luminous globular clusters (GCs) in M31. Here, we determine its age and mass by comparing multicolor photometry with theoretical stellar population synthesis models. Based on far- and near-ultraviolet GALEX photometry, broad-band UBVRI, and infrared JHK_s 2MASS data, we construct the most extensive spectral energy distribution of G1 to date, spanning the wavelength range from 1538 to 20,000 A. A quantitative comparison with a variety of simple stellar population (SSP) models yields a mean age that is consistent with G1 being among the oldest building blocks of M31 and having formed within ~1.7 Gyr after the Big Bang. Irrespective of the SSP model or stellar initial mass function adopted, the resulting mass estimates (of order $10^7 M_\odot$) indicate that G1 is one of the most massive GCs in the Local Group. However, we speculate that the cluster's exceptionally high mass suggests that it may not be a genuine GC. We also derive that G1 may contain, on average, $(1.65\pm0.63)\times10^2 L_\odot$ far-ultraviolet-bright, hot, extreme horizontal-branch stars, depending on the SSP model adopted. On a generic level, we demonstrate that extensive multi-passband photometry coupled with SSP analysis enables one to obtain age estimates for old SSPs to a similar accuracy as from integrated spectroscopy or resolved stellar photometry, provided that some of the free parameters can be constrained independently.

This paper supplements Jiang et al. (2003), who studied 172 M31 globular clusters (GCs) and globular cluster candidates from Battistini et al. (1987) on the basis of integrated photometric measurements in the Beijing-Arizona-Taiwan-Connecticut (BATC) photometric system. Here, we present multicolor photometric CCD data (in the BATC system) for the remaining 39 M31 GCs and candidates. In addition, the ages of 35 GCs are constrained by comparing our accurate photometry with updated theoretical stellar synthesis models. We use photometric measurements from GALEX in the far- and near-ultraviolet and 2MASS infrared $JHK_s$ data, in combination with optical photometry. Except for two clusters, the ages of the other sample GCs are all older than 1 Gyr. Their age distribution shows that most sample clusters are younger than 6 Gyr, with a peak at ~3 Gyr, although the `usual' complement of well-known old GCs (i.e., GCs of similar age as the majority of the Galactic GCs) is present as well.

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.

In this paper, we study the cosmological implications of the 100 square degree Weak Lensing survey (the CFHTLS-Wide, RCS, VIRMOS-DESCART and GaBoDS surveys). We combine these weak lensing data with the cosmic microwave background (CMB) measurements from the WMAP5, BOOMERanG, CBI, VSA, ACBAR, the SDSS LRG matter power spectrum and the Type Ia Supernoave (SNIa) data with the "Union" compilation (307 sample), using the Markov Chain Monte Carlo method to determine the cosmological parameters. Our results show that the \Lambda CDM model remains a good fit to all of these data. For the dynamical dark energy model with time evolving EoS parameterized as w_{\DE}(a) = w_0 + w_a (1-a), we find that the best-fit model implying the mildly preference of Quintom model whose EoS gets across the cosmological constant boundary during evolution. Regarding the total neutrino mass limit, we obtain the upper limit, \sum m_{\nu}< 0.471 eV (95% C.L.) within the framework of the flat \Lambda CDM model. Due to the obvious degeneracies between the neutrino mass and the EoS of dark energy model, this upper limit will be relaxed by a factor of 2 in the framework of dynamical dark energy models. For the constraints on the inflation parameters, we find that the upper limit on the ratio of the tensor to scalar is r<0.35 (95% C.L.) and the inflationary models with the slope n_s\geq1 are excluded at more than 2 \sigma confidence level. In this paper we pay particular attention to the contribution from the weak lensing data and find that the current weak lensing data do improve the constraints on matter density \Omega_m, \sigma_8, \sum{m_{\nu}}, and the EoS of dark energy.

In this paper we study the sensitivity of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) project to the determination of cosmological parameters, employing the Monte Carlo Markov Chains (MCMC) method. For comparison, we first analyze the constraints on cosmological parameters from current observational data, including WMAP, SDSS and SN Ia. We then simulate the 3D matter power spectrum data expected from LAMOST, together with the simulated CMB data for PLANCK and the SN Ia from 5-year Supernovae Legacy Survey (SNLS). With the simulated data, we investigate the future improvement on cosmological parameter constraints, emphasizing the role of LAMOST. Our results show the potential of LAMOST in probing for the cosmological parameters, especially in constraining the equation-of-state (EoS) of the dark energy and the neutrino mass.

We discuss the generic characteristics of stochastic particle acceleration by a fully developed turbulence spectrum and show that resonant interactions of particles with high speed waves dominate the acceleration process. To produce the relativistic electrons inferred from the broadband spectrum of a few well-observed shell-type supernova remnants in the leptonic scenario for the TeV emission, fast mode waves must be excited effectively in the downstream and dominate the turbulence in the subsonic phase. Strong collisionless non-relativistic astrophysical shocks are studied with the assumption of a constant Aflven speed. The energy density of non-thermal electrons is found to be comparable to that of the magnetic field. With reasonable parameters, the model explains observations of shell-type supernova remnants. More detailed studies are warranted to better understand the nature of supernova shocks.

Recently, the WMAP group has published their five-year data and considered the constraints on the time evolving equation of state of dark energy for the first time from the WMAP distance information. In this paper, we study the effectiveness of the usage of these distance information and find that these compressed CMB information can give similar constraints on dark energy parameters compared with the full CMB power spectrum if dark energy perturbations are included, however, once incorrectly neglecting the dark energy perturbations, the difference of the results are sizable.

We study the stochastic electron acceleration by fast mode waves in the turbulent downstream of weakly magnetized collisionless astrophysical shocks. The acceleration is most efficient in a dissipative layer, and the model characteristics are determined by the shock speed, density, magnetic field, and turbulence decay length. The model explains observations of shell-type supernova remnants RX J1713.7-3946 and J0852.0-4622 and can be tested by observations in hard X-rays with the HXMT and NuSTAR or gamma-rays with the GLAST.

Due to the lack of low redshift long Gamma-Ray Bursts (GRBs), the circular problem has been a severe obstacle for using GRBs as cosmological candles. In this paper, we present a new method to deal with such a problem in MCMC global fitting analysis. Assuming that a certain type of correlations between different observables exists in a subsample of GRBs, for the parameters involved in the correlation relation, we treat them as free parameters and determine them simultaneously with cosmological parameters through MCMC analysis on GRB data together with other observational data. Then the circular problem is naturally eliminated in this procedure. We take the Ghirlanda relation as an example while keeping in mind the debate about its physical validity. Together with SNe Ia, WMAP and SDSS data, we include 27 GRBs with the reported Ghirlanda relation in our study, and perform MCMC global fitting. We consider the $\Lambda$CDM model and dynamical dark energy models. In each case, in addition to the constraints on the relevant cosmological parameters, we obtain the best fit values as well as the distributions of the correlation parameters $A$ and $C$. We find that the observational data sets other than GRBs can affect $A$ and $C$ considerably through their degeneracies with the cosmological parameters. The results on $A$ and $C$ for different cosmological models are in well agreement within $1\sigma$ range. The best fit value of $A$ in all models being analyzed is $A\sim 1.53$ with $\sigma \sim 0.08$. For $C$, we have the best value in the range of $0.94-0.98$ with $\sigma\sim 0.1$. It is also noted that the distributions of $A$ and $C$ are generally broader than the priors used in many studies in literature. (Abriged)

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)

Chandra's high resolution observations of radio galaxies require a revisit of the relevant electron acceleration processes. Although the diffusive shock particle acceleration model may explain spectra of spatially unresolved sources, it encounters difficulties in explaining the structure and spectral properties of recently discovered Chandra X-ray features in several low-power radio sources. We argue that these observations strongly suggest stochastic electron acceleration by magnetized turbulence, and show that the simplest stochastic particle acceleration model with energy independent acceleration and escape timescales can overcome most of these difficulties. We use the bright core of the western hotspot of Pictor A as an example to demonstrate the model characteristics, which may be tested with high energy observations.

Using metallicities from the literature, combined with the Revised Bologna Catalogue of photometric data for M31 clusters and cluster candidates (the latter of which is the most comprehensive catalogue of M31 clusters currently available, including 337 confirmed globular clusters -- GCs -- and 688 GC candidates), we determine 443 reddening values and intrinsic colours, and 209 metallicities for individual clusters without spectroscopic observations. This, the largest sample of M31 GCs presently available, is then used to analyse the metallicity distribution of M31 GCs, which is bimodal with peaks at $\rm {[Fe/H]}\approx -1.7$ and -0.7 dex. An exploration of metallicities as a function of radius from the M31 centre shows a metallicity gradient for the metal-poor GCs, but no such gradient for the metal-rich GCs. Our results show that the metal-rich clusters appear as a centrally concentrated spatial distribution; however, the metal-poor clusters tend to be less spatially concentrated. There is no correlation between luminosity and metallicity among the M31 sample clusters, which indicates that self-enrichment is indeed unimportant for cluster formation in M31. The reddening distribution shows that slightly more than half of the GCs are affected by a reddening of $E(B-V) \la 0.2$ mag; the mean reddening value is $E(B-V) = 0.28_{-0.14}^{+0.23}$ mag. The spatial distribution of the reddening values indicates that the reddening on the northwestern side of the M31 disc is more significant than that on the southeastern side, which is consistent with the conclusion that the northwestern side in nearer to us.

In this paper we analyze the constraints on the property of dark energy from cosmological observations. Together with SNe Ia Gold sample, WMAP, SDSS and 2dFGRS data, we include 69 long Gamma-Ray Bursts (GRBs) data in our study and perform global fitting using Markov Chain Monte Carlo (MCMC) technique. Dark energy perturbations are explicitly considered. We pay particular attention to the time evolution of the equation of state of dark energy parameterized as $w_{DE}=w_0+w_a(1-a)$ with $a$ the scale factor of the universe, emphasizing the complementarity of high redshift GRBs to other cosmological probes. It is found that the constraints on dark energy become stringent by taking into account high redshift GRBs, especially for $w_a$, which delineates the evolution of dark energy.

This paper is the third in the series of papers on M81 globular clusters. In this paper, we present spatial and metal abundance properties of 95 M81 globular clusters, which comprise nearly half of all the M81 globular cluster system. These globular clusters are divided into two M81 metallicity groups by a KMM test. Our results show that, the metal-rich clusters did not demonstrate a centrally concentrated spatial distribution as ones in M31, and metal-poor clusters tend to be less spatially concentrated. In other words, the distribution of the metal-rich clusters in M81 is not very similar to that of M31. Most of the metal-rich clusters distribute at projected radii of 4-8 kpc. It is also noted that the metal-rich clusters distribute within the inner 20 kpc, and the metal-poor ones do out to radii of ~40 kpc. Like our Galaxy and M31, the metallicity distribution of globular clusters in M81 along galactocentric radius suggests that some dissipation occurred during the formation of the globular cluster system, i.e. smooth, pressure-supported collapse models of galaxies are unlikely to produce such radial distribution of metallicity presented in this paper. There is not evident correlation between globular cluster luminosity and metallicity in M81 globular clusters. The overwhelming conclusion of this paper seems to be that a more complete and thorough cluster search is needed in M81.

In this paper we investigate the influence of the intrinsic alignment of background galaxies on weak lensing detections of mass concentrations. Specifically, we analyze the number counts of false peaks resulting from intrinsic ellipticities in lensing convergence maps. Including the alignment of source galaxies, the full noise variance from intrinsic ellipticites in convergence $\kappa$-maps can be written as $\sigma^2_0=\sigma^2_{0ran}+\sigma^2_{0corr}$, where $\sigma^2_{0ran}$ is the noise contributed from randomly oriented source galaxies and $\sigma^2_{0corr}$ denotes the additional noise from intrinsic alignments. However, it is observationally difficult to measure $\sigma^2_{0corr}$ and usually only $\sigma^2_{0ran}$ can be estimated in weak lensing observations. Thus the observational signal-to-noise ratio is often defined with respect to $\sigma_{0ran}$, which is denoted as $\nu_{ran}$ in this paper. The true signal-to-noise ratio $\nu$ in terms of $\sigma_0$ is then $\nu=\nu_{ran}/(1+\sigma^2_{0corr}/\sigma^2_{0ran})^{1/2}$. Given a detection threshold on $\nu_{ran}$, a larger value of $\sigma^2_{0corr}/\sigma^2_{0ran}$ leads to a lower threshold on $\nu$ and therefore a larger expected number of false peaks. With $\sigma^2_{0corr}/\sigma^2_{0ran}\sim 10%$, the average number of false peaks with $\nu_{ran}\ge 3.5$ nearly doubles compared to that without considering the alignment, and for $\nu_{ran}\ge 5$, the number is tripled. As a result, the efficiency of weak lensing cluster detection degrades significantly. (Abridged)

Mayall II = G1 is one of the most luminous globular clusters (GCs) known in M31. New deep, high-resolution observations with the Advanced Camera for Surveys on the {\sl Hubble Space Telescope} are used to provide accurate photometric data to the smallest radii yet. In particular, we present the precise variation of ellipticity and position angle, and of surface brightness for the core of the object. Based on these accurate photometric data, we redetermine the structural parameters of G1 by fitting a single-mass isotropic King model. We derive a core radius, $r_c=0.21\pm0.01\arcsec (=0.78\pm0.04 \rm{pc})$, a tidal radius, $r_t=21.8\pm1.1\arcsec (=80.7\pm3.9 \rm{pc})$, and a concentration index $c=\log (r_t/r_c)=2.01\pm0.02$. The central surface brightness is 13.510 mag arcsec$^{-2}$. We also calculate the half-light radius, at $r_h=1.73\pm0.07\arcsec(=6.5\pm0.3 \rm{pc})$. The results show that, within 10 core radii, a King model fits the surface brightness distribution well. We find that this object falls in the same region of the $M_V$ vs. $\log R_h$ diagram as $\omega$ Centauri, M54 and NGC 2419 in the Milky Way. All three of these objects have been claimed to be the stripped cores of now defunct dwarf galaxies. We discuss in detail whether GCs, stripped cores of dwarf spheroidals and normal dwarf galaxies form a continuous distribution in the $M_V$ versus $\log R_h$ plane, or if GCs and dwarf spheroidals constitute distinct classes of objects; we present arguments in favour of this latter view.

We have constrained the age of the globular cluster S312 in the Andromeda galaxy (M31) by comparing its multicolor photometry with theoretical stellar population synthesis models. This is both a check on the age of this globular cluster, as well a check on our methodology. Main-sequence photometry has been the most direct method for determining the age of a star cluster. S312 was observed as part of the Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey from 1995 February to 2003 December. The photometry of BATC images for S312 was taken with 9 intermediate-band filters covering 5000-10000\AA. Combined with photometry in the near-ultraviolet (NUV) of {\sl GALEX}, broad-band $UBVR$ and infrared $JHK_s$ of 2MASS, we obtained the accurate spectral energy distributions (SEDs) of S312 from 2267-20000\AA. A quantitative comparison to simple stellar population models yields an age of $9.5_{-0.99}^{+1.15}$ Gyr, which is in very good agreement with the previous determination by main-sequence photometry. S312 has a mass of $9.8\pm{1.85}\times 10^5 \rm M_\odot$, and is a medium-mass globular cluster in M31. By analysis of errors of ages determined based on the SED fitting method of this paper, secure age constraints are derived with errors of $< 3$ Gyr for ages younger than 9 Gyr. In fact, the theoretical SEDs are not sensitive to the variation of age for ages greater than $\sim 10$ Gyr. Therefore, for globular clusters as old as the majority of the Galactic GCs, our method do not distinguish them accurately. We emphasize that our results show that even with multiband photometry spanning NUV to $K_s$, our age constraints from SED fitting are distressingly uncertain, which has implications for age derivations in extragalactic globular cluster systems.

We combine recent long Gamma Ray Bursts (GRBs) sample including 52 objects out to z$=$6.3 compiled from Swift Gamma Ray Bursts by Schaefer (2006) with Type Ia Supernova (SNIa), Cosmic Microwave Background (CMB) and Baryon Oscillation (BAO) to constrain cosmologies. We study the constraints arising from GRBs alone and the complementarity of GRBs to other cosmological probes. To analyze the cosmological role of GRBs, we adopt the Hubble constant as a free parameter with a prior in the range of $H_0=72\pm8\enspace\kmsmpc$ instead of a fixed $H_0$ used in previous studies. By jointly using SNIa gold/SNLS samples and GRBs, the constraints on \Om -\Ol parameter space are dramatically improved in comparing with those from SNIa data alone. The complementarity of GRBs is mostly from the sub-sample with z$>$1.5 due to the different parameter degeneracies involved in luminosity distances at different redshifts. Including GRB data in addition to SNIa, BAO and CMB in our analysis, we find that the concordance model with \Om $\sim$0.3 and \Ol$\sim$0.7 is still well within 1$\sigma$ confidence range.

The large majority of extragalactic star cluster studies performed to date essentially use multi-colour photometry, combined with theoretical stellar synthesis models, to derive ages, masses, extinction estimates, and metallicities. M31 offers a unique laboratory for studies of globular cluster (GC) systems. In this paper, we obtain new age estimates for 91 M31 globular clusters, based on improved photometric data, updated theoretical stellar synthesis models and sophisticated new fitting methods. In particular, we used photometric measurements from the Two Micron All Sky Survey (2MASS), which, in combination with optical photometry, can partially break the well-known age-metallicity degeneracy operating at ages in excess of a few Gyr. We show robustly that previous age determinations based on photometric data were affected significantly by this age-metallicity degeneracy. Except for one cluster, the ages of our other sample GCs are all older than 1 Gyr. Their age distribution shows populations of young and intermediate-age GCs, peaking at $\sim3$ and 8 Gyr respectively, as well as the "usual" complement of well-known old GCs, i.e., GCs of similar age as the majority of the Galactic GCs. Our results also show that although there is significant scatter in metallicity at any age, there is a noticeable lack of young metal-poor and old metal-rich GCs, which might be indicative of an underlying age-metallicity relationship among the M31 GC population.

There is clear evidence that the $\gamma$-ray emission from active galactic nuclei (AGNs) is attributed to the inverse Compton scatterings in the relativistic blobs near the massive black holes. If the soft seed photons are from the regions outside the blobs, a linear relation between $(\nu F_{\nu,\gamma}/\nu F_{\nu, \rm syn} u^{*})^{1/(1+\alpha)}$ and Doppler factor $\delta$ is expected, where $\nu F_{\nu,\gamma}$ and $\nu F_{\nu, \rm syn}$ are monochromatic $\gamma$-ray and synchrotron fluxes, respectively, and $u^{*}$ is the energy density of soft seed photons \citep{D97}. We estimate the soft photon energy density in the relativistic blobs contributed by the broad line region (BLRs) in these $\gamma$-ray AGNs using their broad-line emission data. The Doppler factors $\delta$ are derived from their radio core and X-ray emission data, based on the assumption that the X-ray emission is produced through synchrotron self-Compton (SSC) scatterings. We find two nearly linear correlations: $(\nu F_{\nu,\gamma}/\nu F_{\rm opt} u^{*})^{1/(1+\alpha)} \propto \delta^{1.09}$, and $(\nu F_{\nu,\gamma}/\nu F_{\rm IR} u^{*})^{1/(1+\alpha)} \propto \delta^{1.20}$, which are roughly consistent with the linear correlation predicted by the theoretical model for external Compton scatterings. Our results imply that the soft seed photons are dominantly from the BLRs in these $\gamma$-ray AGNs.

Joint analyses on X-ray and Sunyaev-Zel'dovich (SZ) effect of a cluster of galaxies can give rise to an estimate on the angular diameter distance to the cluster. With the redshift information of the cluster, the Hubble constant $H_0$ can then be derived. Furthermore, such measurements on a sample of clusters with a range of redshift can potentially be used to discriminate different cosmological models. In this paper, we present statistical studies on the systematic errors in the determination of $H_0$ due to the triaxiality and non-isothermality of clusters of galaxies. Different from many other studies that assume artificially a specific distribution for the intracluster gas, we start from the triaxial model of dark matter halos obtained from numerical simulations. The distribution of the intracluster gas is then derived under the assumption of the hydrodynamic equilibrium. For the equation of state of the intracluster gas, both the isothermal and the polytropic cases are investigated. We run Monte Carlo simulations to generate samples of clusters according to the distributions of their masses, axial ratios, concentration parameters, as well as line-of-sight directions. To mimic observations, the estimation of the Hubble constant is done by fitting X-ray and SZ profiles of a triaxial cluster with the isothermal and spherical $\beta$-model. We find that for a sample of clusters with $M=10^{14}h^{-1}\hbox{M}_{\odot}$ and $z=0.1$, the value of the estimated $H_0$ is positively biased with $H_0^{peak}(estimated)\approx 1.05H_0(true)$ and $H_0^{ave}(estimated)\approx 1.05H_0(true)$ for the isothermal case. For the polytropic case with $\gamma=1.15$, the bias is rather large with $H_0^{peak}(estimated)\approx 1.35H_0(true)$ and $H_0^{ave}(estimated)\approx 3H_0(true)$. (abridged)

This paper presents accurate spectral energy distributions (SEDs) of 16 M31 globular clusters (GCs) confirmed by spectroscopy and/or high spatial-resolution imaging, as well as 30 M31 globular cluster candidates detected by Mochejska et al. Most of these candidates have m_V > 18, deeper than previous searches, and these candidates have not yet been confirmed to be globular clusters. The SEDs of these clusters and candidates are obtained as part of the BATC Multicolor Survey of the Sky, in which the spectrophotometrically-calibrated CCD images of M31 in 13 intermediate-band filters from 4000 to 10000 A were observed. These filters are specifically designed to exclude most of the bright and variable night-sky emission lines including the OH forest. In comparison to the SEDs of true GCs, we find that some of the candidate objects are not GCs in M31. SED fits show that theoretical simple stellar population (SSP) models can fit the true GCs very well. We estimate the ages of these GCs by comparing with SSP models. We find that, the M31 clusters range in age from a few ten Myr to a few Gyr old, as well as old GCs, confirming the conclusion that has been found by Barmby et a, Williams & Hodge, Beasley et al., Burstein et al. and Puzia et al. in their investigations of the SEDs of M31 globular clusters.

Gravitational lensing effects arise from the light ray deflection by all of the mass distribution along the line of sight. It is then expected that weak lensing cluster surveys can provide us true mass-selected cluster samples. With numerical simulations, we analyze the correspondence between peaks in the lensing convergence $\kappa$-map and dark matter halos. Particularly we emphasize the difference between the peak $\kappa$ value expected from a dark matter halo modeled as an isolated and spherical one, which exhibits a one-to-one correspondence with the halo mass at a given redshift, and that of the associated $\kappa$-peak from simulations. For halos with the same expected $\kappa$, their corresponding peak signals in the $\kappa$-map present a wide dispersion. At an angular smoothing scale of $\theta_G=1\hbox{arcmin}$, our study shows that for relatively large clusters, the complex mass distribution of individual clusters is the main reason for the dispersion. The projection effect of uncorrelated structures does not play significant roles. The triaxiality of dark matter halos accounts for a large part of the dispersion, especially for the tail at high $\kappa$ side. Thus lensing-selected clusters are not really mass-selected. (abridged)

This paper includes two parts. The first is to present the spectral energy distributions (SEDs) of 49 globular cluster (GC) X-ray sources in the BATC 13 intermediate-band filters from 3800 to 10000 A, and identify 8 unidentified X-ray sources in M31. Using the X-ray data of Einstein observation from 1979 to 1980, ROSAT HRI observation in 1990, Chandra HRC and ACIS-I observations from 1999 to 2001, and the BATC optical survey from 1995 to 1999, we find 49 GC X-ray sources and 8 new unidentified X-ray sources in the BATC M31 field. By analyzing SEDs and FWHMs, 4 of the 8 X-ray sources may be GC candidates. The second is to present some statistical relationships about 62 GC X-ray sources, of which 58 are already known, and 4 are identified in this paper. The distribution of M31 GC X-ray sources' V mags is bimodal, with peaks at m_v = 15.65 and m_v = 17.89, which is different from the distribution of GC candidates. The distribution of B-V color shows that,the GC X-ray sources seem to be associated preferentially with the redder GCs, in agreement with the previous results. Kolmogorov-Smirnov test shows that the maximum value of the absolute difference of B-V distributions of GC X-ray sources and GCs is D_{max}=0.181, and the probability P=0.068 which means we can reject the hypothesis that the two distributions are the same at the 90.0% confidence level. In the end, we study the correlation between X-ray luminosity (0.3-10 keV) and the optical luminosity (in V band) of the GC X-ray sources in M31, and find that there exits a weak relationship with the linear correlation coefficient r = 0.36 at the confidence level of 98.0%.

With the triaxial density profile of dark matter halos and the corresponding equilibrium gas distribution, we derive two-dimensional Sunyaev-Zeldovich (SZ) effect and X-ray surface brightness profiles for clusters of galaxies. It is found that the contour map of these observables can be well approximated by a series of concentric ellipses with scale-dependent eccentricities. The statistical distribution of their eccentricities (or equivalently axial ratios) is analyzed by taking into account the orientation of clusters with respect to the line of sight and the distribution of the axial ratios and the concentration parameters of dark matter halos. For clusters of mass $10^{13}h^{-1}{M}_{\odot}$ at redshift $z=0$, the axial ratio is peaked at $\eta \sim 0.9$ for both SZ and X-ray profiles. For larger clusters, the deviation from circular distributions is more apparent, with $\eta$ peaked at $\eta \sim 0.85$ for $M=10^{15}h^{-1}{M}_{\odot}$. To be more close to observations, we further study the axial-ratio distribution for mass-limited cluster samples with the number distribution of clusters at different redshifts described by a modified Press-Schechter model. For a mass limit of value $M_{lim}=10^{14}h^{-1}{M}_{\odot}$, the average axial ratio is $<\eta > \sim 0.84$ with a tail extended to $\eta \sim 0.6$. With fast advance of high quality imaging observations of both SZ effect and X-ray emissions, our analyses provide a useful way to probe cluster halo profiles and therefore to test theoretical halo-formation models.

The sizes of the broad-line regions (BLRs) in $\gamma$-ray active galactic nuclei (AGNs) are estimated from their optical continuum luminosity by using the empirical relation between $R_{\rm BLR}$ and $L_{\lambda,\rm opt}$. Using the broad emission line data, we derive the photon energy density in the relativistic blobs near the massive black holes in AGNs. We calculate the power of the broad-line photons Compton up-scattered by the relativistic electrons in the blobs. Compared with observed $\gamma$-ray emission data, the Doppler factors $\delta$ of the blobs for a sample of 36 $\gamma$-ray AGNs are derived, which are in the range of $\sim 3-17$. We estimate the central black hole masses of these $\gamma$-ray AGNs from the sizes of the BLR and the broad line widths. It is found that the black hole masses are in the range of $\sim 10^{8}-10^{10}M_{\odot}$. A significant correlation is found between the Doppler factor $\delta$ and the core dominance parameter $R$. The results are consistent with the external radiation Compton (ERC) models for $\gamma$-ray emission from AGNs. The soft seed photons are probably from the broad-line regions.

Since gravitational lensing effects directly probe inhomogeneities of dark matter, lensing-galaxy cross-correlations can provide us important information on the relation between dark matter and galaxy distributions, i.e., the bias. In this paper, we propose a method to measure the stochasticity/nonlinearity of the galaxy bias through correlation studies of the cosmic shear and galaxy number fluctuations. Specifically, we employ the aperture mass statistics $M_{ap}$ to describe the cosmic shear. We divide the foreground galaxy redshift $z_f<z_s$ into several bins, where $z_s$ is the redshift of the source galaxies, and calculate the quantity $<M_{ap}N_g(z_f)>^2/< N_g^2(z_f)>$ for each redshift bin. Then the ratio of the summation of $<M_{ap}N_g(z_f)>^2/< N_g^2(z_f)>$ over the bins to $<M_{ap}^2>$ gives a measure of the nonlinear/stochastic bias. Here $N_g(z_f)$ is the projected surface number density fluctuation of foreground galaxies at redshift $z_f$, and $M_{ap}$ is the aperture mass from the cosmic-shear analysis. We estimate that for a moderately deep weak-lensing survey with $z_s=1$, source galaxy surface number density $n_b=30 \hbox {gal}/\hbox {arcmin}^2$ and a survey area of $25 \hbox {deg}^2$, the effective $r$-parameter that represents the deviation from the linear and deterministic bias is detectable in the angular range of 1'-10' if $|r-1|\gsim 10%$. For shallow, wide surveys such as the Sloan Digital Sky Survey with $z_s=0.5$, $n_b=5 \hbox {gal}/\hbox {arcmin}^2$, and a survey area of $10^4 \hbox {deg}^2$, a 10% detection of $r$ is possible over the angular range $1'-100'$.

Sunyaev-Zel'dovich Effect (SZE) cluster surveys are anticipated to yield tight constraints on cosmological parameters such as the equation of state of dark energy. In this paper, we study the impact of relativistic corrections of the thermal SZE on the cluster number counts expected from a cosmological model and thus, assuming that other cosmological parameters are known to high accuracies, on the determination of the $w$ parameter and $\sigma_8$ from a SZE cluster survey, where $w=p/\rho$ with $p$ the pressure and $\rho$ the density of dark energy, and $\sigma_8$ is the rms of the extrapolated linear density fluctuation smoothed over $8\hbox{Mpc}h^{-1}$. For the purpose of illustrating the effects of relativistic corrections, our analyses mainly focus on $\nu=353 \hbox{GHz}$ and $S_{lim}=30\hbox{mJy}$, where $\nu$ and $S_{lim}$ are the observing frequency and the flux limit of a survey, respectively. These observing parameters are relevant to the {\it Planck} survey. It is found that from two measurable quantities, the total number of SZE clusters and the number of clusters with redshift $z\ge 0.5$, $\sigma_8$ and $w$ can be determined to a level of $\pm 1%$ and $\pm 8%$, respectively, with $1\sigma$ uncertainties from a survey of $10000\hbox{deg}^2$. Relativistic effects are important in determining the central values of $\sigma_8$ and $w$. If we choose the two quantities calculated relativistically from the flat cosmological model with $\sigma_8=0.8284$ and $w=-0.75$ as input, the derived $\sigma_8$ and $w$ would be 0.819 and -0.81, respectively, if relativistic effects are wrongly neglected. The location of the resulting $\sigma_8$ and $w$ in the $\sigma_8-w$ plane is outside the $3\sigma$ region around the real central $\sigma_8$ and $w$.

Multi-wavelength observations from radio to soft X-ray bands of large-scale galactic spiral structures offer synthesized and comprehensive views of nearby disk galaxies. In the presence of a massive dark-matter halo, the density-wave dynamics on galactic scales involves the stellar disk, the gas disk of interstellar medium (ISM), the magnetic field, and the cosmic-ray gas (CRG). In this paper, we explore the dynamic and electromagnetic interplay between the magnetized ISM disk and CRG disk so that structural and diagnostic features of optical, infrared, and synchrotron radio-continuum emissions from a spiral galaxy can be physically understood. On timescales of galactic density waves, cosmic rays collectively may be treated as a relativistically hot tenuous gas fluid that is tied to the large-scale mean magnetic field in transverse bulk motions but moves otherwise differently along the magnetic field relative to the ISM. For both fast and slow magnetohydrodynamic (MHD) density waves in a composite disk system of magnetized ISM and CRG, the minute CRG mass density enhancement is phase shifted relative to the enhancement of parallel magnetic field. Owing to an extremely small number of cosmic rays, the large-scale magnetic field enhancement dominates in synchrotron radio-continuum emissions (as if the CRG is almost unperturbed) for spiral structural manifestations. In addition to the fast and slow MHD density waves, there also exists a suprathermal MHD wave mode by which CRG adjusts itself with an effective suprathermal sound speed close to the speed of light c.

Mid-term quasi-periodicities in solar coronal mass ejections (CMEs) during the most recent solar maximum cycle 23 are reported here for the first time using the four-year data (February 5, 1999 to February 10, 2003) of the Large Angle Spectrometric Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO). In parallel, mid-term quasi-periodicities in solar X-ray flares (class >M5.0) from the Geosynchronous Operational Environment Satellites (GOES) and in daily averages of Ap index for geomagnetic disturbances from the World Data Center (WDC) at the International Association for Geomagnetism and Aeronomy (IAGA) are also examined for the same four-year time span. Several conceptual aspects of possible equatorially trapped Rossby-type waves at and beneath the solar photosphere are discussed.

The Seyfert galaxy NGC 1097 has an extended neutral hydrogen disk, a companion, a prominent bar and a luminous circumnuclear starburst "ring". Magnetic fields as revealed by nonthermal radio-continuum emissions correlate well with the optical barred spiral structure on large scales, have a gross enhancement overlapping with the optical/infrared "ring", and show a trailing swirl around and within the "ring". We propose a scenario of bar-excited long-trailing fast magnetohydrodynamic (MHD) density waves at the modified inner Lindblad resonance (mILR), physically identified with the outer rim of the "ring". These sustained outgoing long-waves are bounced back by the $Q_M$-barrier in the form of incoming short-trailing waves. The damping of these waves deposits a {\it negative} angular momentum into the magnetized circumnuclear gas disk. Thus, gas materials spiral inward, bring in frozen-in magnetic flux, and accumulate inside the mILR to create a circular zone of high density and magnetic flux vulnerable to massive star formation. Depending on the wave damping efficiency, this process may simultaneously sustain a net mass inflow across the "ring" and toward the nucleus. A wavelet analysis on a Hubble Space Telescope image of central NGC 1097 shows a distinct two-arm spiral structure extended down to the nucleus as a strong evidence for circumnuclear MHD density waves. We predict that magnetic-field observations with improved sensitivity and resolution would reveal a specific correspondence between circumnuclear optical and magnetic field spirals much as those known to exist on large scales in nearby spiral galaxies, including NGC 1097.

We study Sunyaev-Zeldovich Effect (SZE) cluster counts in different cosmologies. It is found that even without the full knowledge of the redshift distribution of SZE clusters, one can still readily distinguish a flat universe with a cosmological constant from an open universe. We divide clusters into a low redshift group (with redshift z<0.5) and a high redshift group (with z>1), and compute the relation r=N(z<0.5)/N(z>1), where N is the number of flux-limited (S_lim) SZE clusters. With about the same total number of SZE clusters N(z>0), the r-value for a flat universe with a consmological constant and that for an open universe occupy different regions in the S_lim -- r plot for the most likely cosmological parameters 0.25<\Omega_0<0.35 and 0.2<\Gamma<0.3, where \Gamma is the shape parameter of the initial power spectrum of density fluctuations. Thus, with a deep SZE cluster survey, the ratio r can reveal, independent of the normalization of the power spectrum, whether we are living in a low-density flat universe or an open universe. Within the flat universe scenario, the SZE cluster-normalized \sigma_8 has also been investigated in this work, where sigma_8 is the r.m.s. density fluctuations within the top-hat scale of 8 h^{-1} Mpc. A functional relation sigma_8 \sim \Omega_0^{-0.13} is found. Combining it with the X-ray cluster-normalized \sigma_8\sim\Omega_0^{-0.52+0.13\Omega_0}, one can put tight constraints on both \Omega_0 and \sigma_8 simulataneously.

We discuss the probability distributions of Fourier modes of cosmological density fields using the central limit theorem as it applies to weighted integrals of random fields. It is shown that if the cosmological principle holds in a certain sense, i.e., the universe approaches homogeneity and isotropy sufficiently rapidly on very large scales, the one-point distribution of each Fourier mode of the density field is Gaussian whether or not the density field itself is Gaussian. Therefore, one-point distributions of the power spectrum obtained from observational data or from simulations are not a good test of whether the density field is Gaussian.