The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

In the era of big data astronomy, next generation telescopes and large sky surveys produce data sets at the TB or even PB level. Due to their large data volumes, these astronomical data sets are extremely difficult to transfer and analyze using personal computers or small clusters. In order to offer better access to data, data centers now generally provide online science platforms that enable analysis close to the data. The Chinese Virtual Observatory (China-VO) is one of the member projects in the International Virtual Observatory Alliance and it is dedicated to providing a research and education environment where globally distributed astronomy archives are simple to find, access, and interoperate. In this study, we summarize highlights of the work conducted at the China-VO, as well the experiences and lessons learned during the full life-cycle management of astronomical data. Finally, We discuss the challenges and future trends for astronomical science platforms.

Context: Instrumental profile (IP) is the basic property of a spectrograph. Accurate IP characterisation is the prerequisite of accurate wavelength solution. It also facilitates new spectral acquisition methods such as the forward modeling and deconvolution. Aims: We investigate an IP modeling method for the fibre-fed echelle spectrograph with the emission lines of the ThAr lamp, and explore the method to evaluate the accuracy of IP characterisation. Methods: The backbone-residual (BR) model is put forward and tested on the fibre-fed High Resolution Spectrograph (HRS) at the Chinese Xinglong 2.16-m Telescope, which is the sum of the backbone function and the residual function. The backbone function is a bell-shaped function to describe the main component and the spatial variation of IP. The residual function, which is expressed as the cubic spline function, accounts for the difference between the bell-shaped function and the actual IP. The method of evaluating the accuracy of IP characterisation is based on the spectral reconstruction and Monte Carlo simulation. Results: The IP of HRS is characterised with the BR model, and the accuracy of the characterised IP reaches 0.006 of the peak value of the backbone function. This result demonstrates that the accurate IP characterisation has been achieved on HRS with the BR model, and the BR model is an excellent choice for accurate IP characterisation of fibre-fed echelle spectrographs.

We perform an extensive review of the numerous studies and methods used to determine the total mass of the Milky Way. We group the various methods into seven broad classes, including: i) estimating Galactic escape velocity using high velocity objects; ii) measuring the rotation curve through terminal and circular velocities; iii) modeling halo stars, globular clusters and satellite galaxies with the Spherical Jeans equation and iv) with phase-space distribution functions; v) simulating and modeling the dynamics of stellar streams and their progenitors; vi) modeling the motion of the Milky Way, M31 and other distant satellites under the framework of Local Group timing argument; and vii) measurements made by linking the brightest Galactic satellites to their counterparts in simulations. For each class of methods, we introduce their theoretical and observational background, the method itself, the sample of available tracer objects, model assumptions, uncertainties, limits and the corresponding measurements that have been achieved in the past. Both the measured total masses within the radial range probed by tracer objects and the extrapolated virial masses are discussed and quoted. We also discuss the role of modern numerical simulations in terms of helping to validate model assumptions, understanding systematic uncertainties and calibrating the measurements. While measurements in the last two decades show a factor of two scatters, recent measurements using \textit{Gaia} DR2 data are approaching a higher precision. We end with a detailed discussion of future developments, especially as the size and quality of the observational data will increase tremendously with current and future surveys. In such cases, the systematic uncertainties will be dominant and thus will necessitate a much more rigorous testing and characterization of the various mass determination methods.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$\alpha$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$\alpha$ and high-$\alpha$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$\alpha$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$\alpha$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$\alpha$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$\alpha$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M=1.23_{-0.18}^{+0.21}\times 10^{12} M_\odot$ and a concentration of $c=9.4_{ -2.1}^{ +2.8}$ with the prior based on the $M$-$c$ relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. By varying the disruption rate of such satellites in this simulation, we estimate a $\sim 5\%$ scatter in the inferred MW halo mass among hydrodynamics-based simulations.

The WorldWide Telescope(WWT) is a scientific visualization platform which can browse deep space images, star catalogs, and planetary remote sensing data from different observation facilities in a three-dimensional virtual scene. First launched and then open-sourced by Microsoft Research, the WWT is now managed by the American Astronomical Society (AAS). Hierarchical Progressive Survey (HiPS) is an astronomical data release scheme proposed by Centre de Donn\'ees astronomiques de Strasbourg (CDS) and has been accepted as a recommendation by International Virtual Observatory Alliance (IVOA). The HiPS solution has been adopted widely by many astronomical institutions for data release. Since WWT selected Hierarchical Triangular Mesh (HTM) as the standard for data visualization in the early stage of development, data released by HiPS cannot be visualized in WWT, which significantly limits the application of WWT. This paper introduces the implementation method for HiPS dataset visualization in WWT, and introduces HiPS data projection, mesh rendering, and data index implementation in WWT. Taking Chang'E-2 lunar probe data as an example, this paper introduces how to convert planetary remote sensing data into a HiPS dataset and integrate it into WWT. This paper also compares the efficiency and memory consumption of WWT loading its native data and HiPS data, and illustrates the application of HiPS in scientific data visualization and science education in WWT.

The intergalactic medium (IGM) at $z\sim$ 5 to 6 is largely ionized, and yet the main source for the IGM ionization in the early universe is uncertain. Of the possible contributors are faint quasars with $-26 \lesssim M_{\rm 1450} \lesssim -23$, but their number density is poorly constrained at $z\sim5$. In this paper, we present our survey of faint quasars at $z\sim5$ in the European Large-Area {\it ISO} Survey-North 1 (ELAIS-N1) field over a survey area of 6.51 deg$^2$ and examine if such quasars can be the dominant source of the IGM ionization. We use the deep optical/near-infrared data of the ELAIS-N1 field as well as the additional medium-band observations to find $z \sim 5$ quasars through a two-step approach using the broadband color selection, and SED fitting with the medium-band information included. Adopting Bayesian information criterion, we identify ten promising quasar candidates. Spectra of three of the candidates are obtained, confirming all of them to be quasars at $z\sim5$ and supporting the reliability of the quasar selection. Using the promising candidates, we derive the $z\sim5$ quasar luminosity function at $-26 \lesssim M_{\rm 1450} \lesssim -23$. The number density of faint $z\sim5$ quasars in the ELAIS-N1 field is consistent with several previous results that quasars are not the main contributors to the IGM-ionizing photons at $z\sim5$

Sensitive radio continuum surveys of the Galactic plane are ideal for discovering new supernova remnants (SNRs). From the Sino-German {\lambda}6 cm polarisation survey of the Galactic plane, an extended shell-like structure has been found at l = 21.8 degree, b = -3.0 degree, which has a size of about 1 degree. New observations were made with the Effelsberg 100-m radio telescope at {\lambda}11 cm to estimate the source spectrum together with the Urumqi {\lambda}6 cm and the Effelsberg {\lambda}21 cm data. The spectral index of G21.8-3.0 was found to be {\alpha} = -0.72 {\pm} 0.16. Polarised emission was mostly detected in the eastern half of G21.8-3.0 at both {\lambda}6 cm and {\lambda}11 cm. These properties, together with the H{\alpha} filament along its northern periphery and the lack of infrared emission, indicate that the emission is non-thermal as is usual in shell-type SNRs.

We apply the Fourier Power Function Shapelets (FPFS) shear estimator to the first year data of the Hyper Suprime-Cam survey to construct a shape catalog. The FPFS shear estimator has been demonstrated to have multiplicative bias less than one percent in the absence of blending, regardless of complexities of galaxy shapes, smears of point spread functions (PSFs) and contamination from noise. The blending bias is calibrated with realistic image simulations, which include the impact of neighboring objects, using the COSMOS Hubble Space Telescope images. Here we carefully test the influence of PSF model residual on the FPFS shear estimation and the uncertainties in the shear calibration. Internal null tests are conducted to characterize potential systematics in the FPFS shape catalog and the results are compared with those measured using a catalog where the shapes were estimated using the re-Gaussianization algorithms. Furthermore, we compare various weak lensing measurements between the FPFS shape catalog and the re-Gaussianization shape catalog and conclude that the weak lensing measurements between these two shape catalogs are consistent with each other within the statistical uncertainty.

We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.

We propose a versatile and accurate method to estimate the halo mass and concentration from the kinematics of satellite galaxies. We construct the 6D phase-space distribution function of satellites from a cosmological simulation based on the similarity of internal dynamics for different halos. Within the Bayesian statistical framework, not only can we infer the halo mass and concentration efficiently, but also treat various observational effects, including the selection function, incomplete data, and measurement errors, in a rigorous and straightforward manner. Through tests with mock samples, we show that our method is valid and accurate, and more precise than pure steady-state methods. It can constrain the halo mass to within ~ 20% using only 20 tracers and has a small intrinsic uncertainty of ~ 10%. In addition to the clear application to the Milky Way and similar galaxies, our method can be extended to galaxy groups or clusters.

The inference of the Milky Way halo mass requires modelling the phase space structure of dynamical tracers, with different tracers following different models and having different levels of sensitivity to the halo mass. For steady-state models, phase correlations among tracer particles lead to an irreducible stochastic bias. This bias is small for satellite galaxies and dark matter particles, but as large as a factor of 2 for halo stars. This is consistent with the picture that satellite galaxies closely trace the underlying phase space distribution of dark matter particles, while halo stars are less phase-mixed. As a result, the use of only $\sim 100$ satellite galaxies can achieve a significantly higher accuracy than that achievable with a much larger sample of halo stars.

Much data on the Galactic polarized radio emission has been gathered in the last five decades. All-sky surveys have been made, but only in narrow, widely spaced frequency bands, and the data are inadequate for the characterization of Faraday rotation, the main determinant of the appearance of the polarized radio sky at decimetre wavelengths. We describe a survey of the polarized radio emission from the Southern sky, aiming to characterize the magneto-ionic medium, particularly the strength and configuration of the magnetic field. This work is part of the Global Magneto-Ionic Medium Survey (GMIMS). We have designed and built a feed and receiver covering the band 300 to 900 MHz for the CSIRO Parkes 64-m Telescope. We have surveyed the entire sky between declinations -90 and +20 degrees. We present data covering 300 to 480 MHz with angular resolution 81' to 45'. The survey intensity scale is absolutely calibrated, based on measurements of resistors at known temperatures and on an assumed flux density and spectral index for Taurus A. Data are presented as brightness temperatures. We have applied Rotation Measure Synthesis to the data to obtain a Faraday depth cube of resolution 5.9 radians per metre squared, sensitivity of 60 mK of polarized intensity, and angular resolution 1.35 degrees. The data presented in this paper are available at the Canadian Astronomy Data Centre.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

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

LAMOST DR5 released more than 200,000 low resolution spectra of early-type stars with S/N>50. Searching for metallic-line (Am) stars in such a large database and study of their statistical properties are presented in this paper. Six machine learning algorithms were experimented with using known Am spectra, and both the empirical criteria method(Hou et al. 2015) and the MKCLASS package(Gray et al. 2016) were also investigated. Comparing their performance, the random forest (RF) algorithm won, not only because RF has high successful rate but also it can derives and ranks features. Then the RF was applied to the early type stars of DR5, and 15,269 Am candidates were picked out. Manual identification was conducted based on the spectral features derived from the RF algorithm and verified by experts. After manual identification, 9,372 Am stars and 1,131 Ap candidates were compiled into a catalog. Statistical studies were conducted including temperature distribution, space distribution, and infrared photometry. The spectral types of Am stars are mainly between F0 and A4 with a peak around A7, which is similar to previous works. With the Gaia distances, we calculated the vertical height Z from the Galactic plane for each Am star. The distribution of Z suggests that the incidence rate of Am stars shows a descending gradient with increasing jZj. On the other hand, Am stars do not show a noteworthy pattern in the infrared band. As wavelength gets longer, the infrared excess of Am stars decreases, until little or no excess in W1 and W2 bands.

We study the cross-correlation between the Planck CMB lensing convergence map and the eBOSS quasar overdensity obtained from the Sloan Digital Sky Survey (SDSS) IV, in the redshift range $0.9 < z < 2.2$. We detect the CMB lensing convergence-quasar cross power spectrum at $5.4 \sigma$ significance. The cross power spectrum provides a quasar clustering bias measurement that is expected to be particularly robust against systematic effects. The redshift distribution of the quasar sample has a median redshift $z \approx 1.55$, and an effective redshift about $1.51$. The best fit bias of the quasar sample is $b_q = 2.43 \pm 0.45$, corresponding to a host halo mass of $\log_{10}\left( \frac{M}{h^{-1} M_\odot} \right) = 12.54^{+0.25}_{-0.36}$. This is broadly consistent with the previous literature on quasars with a similar redshift range and selection. Since our constraint on the bias comes from the cross-correlation between quasars and CMB lensing, we expect it to be robust to a wide range of possible systematic effects that may contaminate the auto correlation of quasars. We checked for a number of systematic effects from both CMB lensing and quasar overdensity, and found that all systematics are consistent with null within $2 \sigma$. The data is not sensitive to a possible scale dependence of the bias at present, but we expect that as the number of quasars increases (in future surveys such as DESI), it is likely that strong constraints on the scale dependence of the bias can be obtained.

Large-scale radio continuum surveys provide data to get insights into the physical properties of radio sources. HII regions are prominent radio sources produced by thermal emission of ionised gas around young massive stars. We identify and analyse HII regions in the Sino-German 6cm polarisation survey of the Galactic plane. Objects with flat radio continuum spectra together with infrared and/or Halpha emission were identified as HII regions. For HII regions with small apparent sizes, we cross-matched the 6cm small-diameter source catalogue with the radio HII region catalogue compiled by Paladini and the infrared HII region catalogue based on the WISE data. Extended HII regions were identified by eye by overlaying the Paladini and the WISE HII regions onto the 6cm survey images for coincidences. The TT-plot method was employed for spectral index verification. A total of 401 HII regions were identified and their flux densities were determined with the Sino-German 6cm survey data. In the surveyed area, 76 pairs of sources are found to be duplicated in the Paladini HII region catalogue, mainly due to the non-distinction of previous observations with different angular resolutions, and 78 objects in their catalogue are misclassified as HII regions, being actually planetary nebulae, supernova remnants or extragalactic sources that have steep spectra. More than 30 HII regions and HII region candidates from our 6cm survey data, especially extended ones, do not have counterparts in the WISE HII region catalogue, of which 9 are identified for the first time. Based on the newly derived radio continuum spectra and the evidence of infrared emission, the previously identified SNRs G11.1-1.0, G20.4+0.1 and G16.4-0.5 are believed to be HII regions.

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

Intensive Monitoring Survey of Nearby Galaxies (IMSNG) is a high cadence observation program monitoring nearby galaxies with high probabilities of hosting supernovae (SNe). IMSNG aims to constrain the SN explosion mechanism by inferring sizes of SN progenitor systems through the detection of the shock-heated emission that lasts less than a few days after the SN explosion. To catch the signal, IMSNG utilizes a network of 0.5-m to 1-m class telescopes around the world and monitors the images of 60 nearby galaxies at distances D < 50 Mpc to a cadence as short as a few hours. The target galaxies are bright in near-ultraviolet (NUV) with M_NUV < -18.4 AB mag and have high probabilities of hosting SNe (0.06 SN/yr per galaxy). With this strategy, we expect to detect the early light curves of 3.4 SNe per year to a depth of R ~ 19.5 mag, enabling us to detect the shock-heated emission from a progenitor star with a radius as small as 0.1 R_sun. The accumulated data will be also useful for studying faint features around the target galaxies and other science projects. So far, 18 SNe have occurred in our target fields (16 in IMSNG galaxies) over 5 years, confirming our SN rate estimate of 0.06 SN/yr per galaxy.

We spectroscopically identify 101 Galactic HII regions using spectra from the Large Sky Area Multi- Object Fiber Spectroscopic Telescope (LAMOST) survey, cross-matched with an HII region catalog derived from the all-sky Wide-Field Infrared Survey Explorer(WISE) data. Among all HII regions in our sample, 47 sources are newly confirmed. Spatially, most of our identified HII regions are located in the anti-center area of the Galaxy. For each of the HII regions, we accurately extract and measure the nebular emission lines of the spectra, and estimate the oxygen abundances using the strong-line method. We focus on the abundance distribution of HII regions in the Galactic anti-center area. Accordingly, we derive the oxygen abundance gradient with a slope of -0.036 +/-0.004 dex/kpc, covering a range of RG from 8.1 to 19.3 kpc. In particular, we also fit the outer disk objects with a slope of -0.039 +/- 0.012 dex /kpc, which indicates that there is no flattening of the radial oxygen gradient in the outer Galactic disk.

We develop a novel approach in exploring the joint dependence of halo bias on multiple halo properties using Gaussian process regression. Using a $\Lambda$CDM $N$-body simulation, we carry out a comprehensive study of the joint bias dependence on halo structure, formation history and environment. We show that the bias is a multivariate function of halo properties that falls into three regimes. For massive haloes, halo mass explains the majority of bias variation. For early-forming haloes, bias depends sensitively on the recent mass accretion history. For low-mass and late-forming haloes, bias depends more on the structure of a halo such as its shape and spin. Our framework enables us to convincingly prove that $V_\mathrm{max}/V_\mathrm{vir}$ is a lossy proxy of formation time for bias modelling, whereas the mass, spin, shape and formation time variables are non-redundant with respect to each other. Combining mass and formation time largely accounts for the mass accretion history dependence of bias. Combining all the internal halo properties fully accounts for the density profile dependence inside haloes, and predicts the clustering variation of individual haloes to a $20\%$ level at $\sim 10\mathrm{Mpc}h^{-1}$. When an environmental density is measured outside $1\mathrm{Mpc}h^{-1}$ from the halo centre, it outperforms and largely accounts for the bias dependence on the internal halo structure, explaining the bias variation above a level of $30\%$.

We identify a sample of 1959 massive clusters of galaxies in the redshift range of 0.7<z<1.0 from the survey data of Sloan Digital Sky Survey (SDSS) and Wide-field Infrared Survey Explorer (WISE). These clusters are recognized as the overdensity regions around the SDSS luminous red galaxies, having a richness greater than 15 or an equivalent mass M_{500} \ge 2.5*10^{14}M_{\odot}. Among them, 1505 clusters are identified for the first time, which significantly enlarge the number of high-redshift clusters of z>0.75. By comparing them with clusters at lower redshifts, we confirm that richer clusters host more luminous brightest cluster galaxies (BCGs) also at high redshifts, and that the fraction of blue galaxies is larger in clusters at higher redshifts. A small fraction of BCGs show ongoing star formation or active nuclei. The number density profile of member galaxies in stacked samples of clusters shows no significant redshift evolution.

Context. Pulsars scintillate. Dynamic spectra show brightness variation of pulsars in the time and frequency domain. Secondary spectra demonstrate the distribution of fluctuation power in the dynamic spectra. Aims. Dynamic spectra strongly depend on observational frequencies, but were often observed at frequencies lower than 1.5 GHz. Scintillation observations at higher frequencies help to constrain the turbulence feature of the interstellar medium over a wide frequency range and can detect the scintillations of more distant pulsars. Methods. Ten pulsars were observed at 2250 MHz (S-band) with the Jiamusi 66 m telescope to study their scintillations. Their dynamic spectra were first obtained, from which the decorrelation bandwidths and time scales of diffractive scintillation were then derived by autocorrelation. Secondary spectra were calculated by forming the Fourier power spectra of the dynamic spectra. Results. Most of the newly obtained dynamic spectra are at the highest frequency or have the longest time span of any published data for these pulsars. For PSRs B0540+23, B2324+60 and B2351+61, these were the first dynamic spectra ever reported. The frequencydependence of scintillation parameters indicates that the intervening medium can rarely be ideally turbulent with a Kolmogorov spectrum. The thin screen model worked well at S-band for the scintillation of PSR B1933+16. Parabolic arcs were detected in the secondary spectra of three pulsars, PSRs B0355+54, B0540+23 and B2154+40, all of which were asymmetrically distributed. The inverted arclets of PSR B0355+54 were seen to evolve along the main parabola within a continuous observing session of 12 hours, from which the angular velocity of the pulsar was estimated that was consistent with the measurement by very long baseline interferometry (VLBI).

Context. Pulsars scintillate. Dynamic spectra show brightness variation of pulsars in the time and frequency domain. Secondary spectra demonstrate the distribution of fluctuation power in the dynamic spectra. Aims. Dynamic spectra strongly depend on observational frequencies, but were often observed at frequencies lower than 1.5 GHz. Scintillation observations at higher frequencies help to constrain the turbulence feature of the interstellar medium over a wide frequency range and can detect the scintillations of more distant pulsars. Methods. Ten pulsars were observed at 2250 MHz (S-band) with the Jiamusi 66 m telescope to study their scintillations. Their dynamic spectra were first obtained, from which the decorrelation bandwidths and time scales of diffractive scintillation were then derived by autocorrelation. Secondary spectra were calculated by forming the Fourier power spectra of the dynamic spectra. Results. Most of the newly obtained dynamic spectra are at the highest frequency or have the longest time span of any published data for these pulsars. For PSRs B0540+23, B2324+60 and B2351+61, these were the first dynamic spectra ever reported. The frequencydependence of scintillation parameters indicates that the intervening medium can rarely be ideally turbulent with a Kolmogorov spectrum. The thin screen model worked well at S-band for the scintillation of PSR B1933+16. Parabolic arcs were detected in the secondary spectra of three pulsars, PSRs B0355+54, B0540+23 and B2154+40, all of which were asymmetrically distributed. The inverted arclets of PSR B0355+54 were seen to evolve along the main parabola within a continuous observing session of 12 hours, from which the angular velocity of the pulsar was estimated that was consistent with the measurement by very long baseline interferometry (VLBI).

Pulsar dispersion measure (DM) accounts for the total electron content between a pulsar and us. High-precision observations for projects of pulsar timing arrays show temporal DM variations of millisecond pulsars. The aim of this paper is to decompose the DM variations of 30 millisecond pulsars by using Hilbert-Huang Transform (HHT) method, so that we can determine the general DM trends from interstellar clouds and the annual DM variation curves from solar wind, interplanetary medium and/or ionosphere. We find that the decomposed annual variation curves of 22 pulsars exhibit quasi-sinusoidal, one component and double components features of different origins. The amplitudes and phases of the curve peaks are related to ecliptic latitude and longitude of pulsars, respectively.

We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a synthetic stellar halo from mergers of nine satellites in the Milky Way. We construct the mock catalogue by extracting a heliocentric volume of 10 kpc from our simulations and assigning expected observational uncertainties corresponding to bright stars from Gaia DR2 and LAMOST DR5. We compare the results from StarGO against that from a Friends-of-Friends (FoF) based method in the space of orbital energy and angular momentum. We show that StarGO is able to systematically identify more satellites and achieve higher number fraction of identified stars for most of the satellites within the extracted volume. When applied to data from Gaia DR2, StarGO will enable us to reveal the origins of the inner stellar halo in unprecedented detail.

The spherical Jeans equation (SJE) is widely used in dynamical modelling of the Milky Way (MW) halo potential. We use haloes and galaxies from the cosmological Millennium-II simulation and hydrodynamical APOSTLE simulations to investigate the performance of the SJE in recovering the underlying mass profiles of MW mass haloes. The best-fitting halo mass and concentration parameters scatter by 25% and 40% around their input values, respectively, when dark matter particles are used as tracers. This scatter becomes as large as a factor of 3 when using star particles instead. This is significantly larger than the estimated statistical uncertainty associated with the use of the SJE. The existence of correlated phase-space structures that violate the steady state assumption of the SJE as well as non-spherical geometries are the principal sources of the scatter. Binary haloes show larger scatter because they are more aspherical in shape and have a more perturbed dynamical state. Our results confirm the previous study of Wang et al. (2017) that the number of independent phase-space structures sets an intrinsic limiting precision on dynamical inferences based on the steady state assumption. Modelling with a radius-independent velocity anisotropy, or using tracers within a limited outer radius, result in significantly larger scatter, but the ensemble-averaged measurement over the whole halo sample is approximately unbiased.

We identify 47,600 clusters of galaxies from photometric data of Two Micron All Sky Survey (2MASS), Wide-field Infrared Survey Explorer (WISE), and SuperCOSMOS, among which 26,125 clusters are recognized for the first time and mostly in the sky outside the Sloan Digital Sky Survey (SDSS) area. About 90% of massive clusters of M_{500}>3*10^{14} M_{\odot} in the redshift range of 0.025<z<0.3 have been detected from such survey data, and the detection rate drops down to 50% for clusters with a mass of M_{500}~1*10^{14} M_{\odot}. Monte Carlo simulations show that the false detection rate for the whole cluster sample is less than 5%. By cross-matching with ROSAT and XMM-Newton sources, we get 779 new X-ray cluster candidates which have X-ray counterparts within a projected offset of 0.2 Mpc.

With many large science equipment constructing and putting into use, astronomy has stepped into the big data era. The new method and infrastructure of big data processing has become a new requirement of many astronomers. Cloud computing, Map/Reduce, Hadoop, Spark, etc. many new technology has sprung up in recent years. Comparing to the high performance computing(HPC), Data is the center of these new technology. So, a new computing architecture infrastructure is necessary, which can be shared by both HPC and big data processing. Based on Astronomy Cloud project of Chinese Virtual Observatory (China-VO), we have made much efforts to optimize the designation of the hybrid computing platform. which include the hardware architecture, cluster management, Job and Resource scheduling.

Public observatory project is playing more and more important role in science popularization education and scientific research, and many amateur astronomers also have began to build their own observatories in remote areas. As a result of the limitation of technical condition and construction funds for amateur astronomers, their system often breaks down, and then a stable remote unattended operation system becomes very critical. Hardware connection and control is the basic and core part in observatory design. Here we propose a conception of engineering hardware design for public observatory operation as a bridge between observatory equipment and observation software. It can not only satisfy multiple observation mode requirement, but also save cost.

We present the measurements of Faraday rotation for 477 pulsars observed by the Parkes 64-m radio telescope and the Green Bank 100-m radio telescope. Using these results along with previous measurements for pulsars and extra-galactic sources, we analyse the structure of the large-scale magnetic field in the Galactic disk. Comparison of rotation measures of pulsars in the disk at different distances as well as with rotation measures of background radio sources beyond the disk reveals large-scale reversals of the field directions between spiral arms and interarm regions. We develop a model for the disk magnetic field, which can reproduce not only these reversals but also the distribution of observed rotation measures of background sources.

We study the X-ray phases of the Crab pulsar utilizing the 11-year observations from the Rossi X-ray Timing Explorer, 6-year radio observations from the Nanshan Telescope, and the ephemeris from Jodrell Bank Observatory. It is found that the X-ray phases in different energy bands and the radio phases from Nanshan Telescope show similar behaviors, including long-time evolution and short-time variations. Such strong correlations between the X-ray and radio phases imply that the radio and X-ray timing noises are both generated from the pulsar spin that cannot be well described by the the monthly ephemeris from the Jodrell Bank observatory. When using the Nanshan phases as references to study the X-ray timing noise, it has a significantly smaller variation amplitude and shows no long-time evolution, with a change rate of $(-1.1\pm1.1)\times10^{-7}$ periods per day. These results show that the distance of the X-ray and radio emission regions on the Crab pulsar has no detectable secular change, and it is unlikely that timing-noises resulted from any unique physical processes in the radio or X-ray emitting regions. The similar behaviors of the X-ray and radio timing noises also imply that the variation of the interstellar medium is not the origin of the Crab pulsar's timing noises, which is consistent with the results obtained from the multi-frequency radio observations of PSR B1540-06.

Dark matter subhalos are the remnants of (incomplete) halo mergers. Identifying them and establishing their evolutionary links in the form of merger trees is one of the most important applications of cosmological simulations. The Hierachical Bound-Tracing (HBT) code identifies halos as they form and tracks their evolution as they merge, simultaneously detecting subhalos and building their merger trees. Here we present a new implementation of this approach, HBT+, that is much faster, more user friendly, and more physically complete than the original code. Applying HBT+ to cosmological simulations we show that both the subhalo mass function and the peak-mass function are well fit by similar double-Schechter functions.The ratio between the two is highest at the high mass end, reflecting the resilience of massive subhalos that experience substantial dynamical friction but limited tidal stripping. The radial distribution of the most massive subhalos is more concentrated than the universal radial distribution of lower mass subhalos. Subhalo finders that work in configuration space tend to underestimate the masses of massive subhalos, an effect that is stronger in the host centre. This may explain, at least in part, the excess of massive subhalos in galaxy cluster centres inferred from recent lensing observations. We demonstrate that the peak-mass function is a powerful diagnostic of merger tree defects, and the merger trees constructed using HBT+ do not suffer from the missing or switched links that tend to afflict merger trees constructed from more conventional halo finders. We make the HBT+ code publicly available.

Many dynamical models of the Milky Way halo require assumptions that the distribution function of a tracer population should be independent of time (i.e., a steady state distribution function) and that the underlying potential is spherical. We study the limitations of such modelling by applying a general dynamical model with minimal assumptions to a large sample of galactic haloes from cosmological $N$-body and hydrodynamical simulations. Using dark matter particles as dynamical tracers, we find that the systematic uncertainties in the measured mass and concentration parameters typically have an amplitude of 25% to 40%. When stars are used as tracers, however, the systematic uncertainties can be as large as a factor of $2-3$. The systematic uncertainties are not reduced by increasing the tracer sample size and vary stochastically from halo to halo. These systematic uncertainties are mostly driven by underestimated statistical noise caused by correlated phase-space structures that violate the steady state assumption. The number of independent phase-space structures inferred from the uncertainty level sets a limiting sample size beyond which a further increase no longer significantly improves the accuracy of dynamical inferences. The systematic uncertainty level is determined by the halo merger history, the shape and environment of the halo. Our conclusions apply generally to any spherical steady-state model.

In certain theories of modified gravity, solar system constraints on deviations from general relativity (GR) are satisfied by virtue of a so-called screening mechanism, which enables the theory to revert to GR in regions where the matter density is high or the gravitational potential is deep. In the case of chameleon theories, the screening has two contributions -- self-screening, which is due to the mass of an object itself, and environmental screening, which is caused by the surrounding matter -- which are often entangled, with the second contribution being more crucial for less massive objects. A quantitative understanding of the effect of the environment on the screening can prove critical in observational tests of such theories using systems such as the Local Group and dwarf galaxies, for which the environment may be inferred in various ways. We use the high-resolution {\sc liminality} simulation of Shi et al. (2015) to test the fidelity of different definitions of environment. We find that, although the different ways to define environment in practice do not agree with one another perfectly, they can provide useful guidance, and cross checks about how well a dark matter halo is screened. In addition, the screening of subhaloes in dark matter haloes is primarily determined by the environment, with the subhalo mass playing a minor role, which means that lower-resolution simulations where subhaloes are not well resolved can still be useful for understanding the modification of gravity inside subhaloes.

Primary cosmic-ray elemental spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment since 2004. The third CREAM payload (CREAM-III) flew for 29 days during the 2007-2008 Antarctic season. Energies of incident particles above 1 TeV are measured with a calorimeter. Individual elements are clearly separated with a charge resolution of ~0.12 e (in charge units) and ~0.14 e for protons and helium nuclei, respectively, using two layers of silicon charge detectors. The measured proton and helium energy spectra at the top of the atmosphere are harder than other existing measurements at a few tens of GeV. The relative abundance of protons to helium nuclei is 9.53+-0.03 for the range of 1 TeV/n to 63 TeV/n. The ratio is considerably smaller than other measurements at a few tens of GeV/n. The spectra become softer above ~20 TeV. However, our statistical uncertainties are large at these energies and more data are needed.

We determine spectral indices of 228 pulsars by using Parkes pulsar data observed at 1.4 GHz, among which 200 spectra are newly determined. The indices are distributed in the range from -4.84 to -0.46.Together with known pulsar spectra from literature, we tried to find clues to the pulsar emission process. The weak correlations between the spectral index, the spin-down energy loss rate $\dot{E}$ and the potential drop in the polar gap $\Delta\Psi$ hint that emission properties are related to the particle acceleration process in a pulsar's magnetosphere.

We determine spectral indices of 228 pulsars by using Parkes pulsar data observed at 1.4 GHz, among which 200 spectra are newly determined. The indices are distributed in the range from -4.84 to -0.46.Together with known pulsar spectra from literature, we tried to find clues to the pulsar emission process. The weak correlations between the spectral index, the spin-down energy loss rate $\dot{E}$ and the potential drop in the polar gap $\Delta\Psi$ hint that emission properties are related to the particle acceleration process in a pulsar's magnetosphere.

Highly polarized components of pulse profiles are investigated by analyzing observational data and simulating the emission processes. The highly polarized components appear at the leading or trailing part of a pulse profile, which preferably have a flat spectrum and a flat polarization angle curve compared with the low polarized components. By considering the emission processes and propagation effects, we simulate the distributions of wave modes and fractional linear polarization within the entire pulsar emission beam. We show that the highly polarized components can appear at the leading, central, and/or trailing parts of pulse profiles in the models, depending on pulsar geometry. The depolarization is caused by orthogonal modes or scattering. When a sight line cuts across pulsar emission beam with a small impact angle, the detected highly polarized component will be of the O mode, and have a flat polarization angle curve and/or a flat spectrum as observed. Otherwise, the highly polarized component will be of the X mode and have a steep polarization angle curve.

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

We calculate the correlation function of 79,091 galaxy clusters in the redshift region of $z \leq 0.5$ selected from the WH15 cluster catalog. With a weight of cluster mass, a significant baryon acoustic oscillation (BAO) peak is detected on the correlation function with a significance of $3.7 \sigma$. By fitting the correlation function with a $\Lambda$CDM model curve, we find $D_v(z = 0.331) r_d^{fid}/r_d = 1261.5 \pm 48$~Mpc which is consistent with the Planck 2015 cosmology. We find that the correlation function of the higher mass sub-sample shows a higher amplitude at small scales of $r < 80~h^{-1}{\rm Mpc}$, which is consistent with our previous result. The 2D correlation function of this large sample of galaxy clusters shows a faint BAO ring with a significance of $1.8\sigma$, from which we find that the distance scale parameters on directions across and along the line-of-sight are $\alpha_{\sigma} = 1.02 \pm 0.06$ and $\alpha_{\pi} = 0.94 \pm 0.10$, respectively.

By cross-matching the currently largest optical catalog of galaxy clusters and the NVSS radio survey database, we obtain the largest complete sample of brightest cluster galaxies (BCGs) in the redshift range of 0.05<z<0.45, which have radio emission and redshift information. We confirm that more powerful radio BCGs tend to be these optically very bight galaxies located in more relaxed clusters. We derived the radio luminosity functions of BCGs from the largest complete sample of BCGs, and find that the functions depend on the optical luminosity of BCGs and the dynamical state of galaxy clusters. However, the radio luminosity function does not show significant evolution with redshift.

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

We report the discovery of a radio halo in the massive merging cluster MACSJ2243.3-0935, as well as a new radio relic candidate, using the Giant Meterwave Radio Telescope and the KAT-7 telescope. The radio halo is coincident with the cluster X-ray emission and has a largest linear scale of approximately 0.9 Mpc. We measure a flux density of $10.0\pm 2.0$ mJy at 610 MHz for the radio halo. We discuss equipartition estimates of the cluster magnetic field and constrain the value to be of the order of 1 $\mu$G. The relic candidate is detected at the cluster virial radius where a filament meets the cluster. The relic candidate has a flux density of $5.2\pm 0.8$ mJy at 610 MHz. We discuss possible origins of the relic candidate emission and conclude that the candidate is consistent with an infall relic.

We observed total and polarized radio continuum emission from the spiral galaxy M 101 at 6.2 cm and 11.1 cm wavelengths with the Effelsberg telescope. We use these data to study various emission components in M 101 and properties of the magnetic field. Separation of thermal and non-thermal emission shows that the thermal emission is closely correlated with the spiral arms, while the non-thermal emission is more smoothly distributed indicating diffusion of cosmic ray electrons away from their places of origin. The radial distribution of both emissions has a break near R=16 kpc, where it steepens to an exponential scale length of about 5 kpc, which is about 2.5 times smaller than at R<16 kpc. The distribution of the polarized emission has a broad maximum near R=12 kpc and beyond R=16 kpc also decreases with about 5 kpc scalelength. It seems that near R=16 kpc a major change in the structure of M 101 takes place, which also affects the distributions of the strength of the random and ordered magnetic field. Beyond R=16 kpc the radial scale length of both fields is about 20 kpc, which implies that they decrease to about 0.3 \mu G at R=70 kpc, which is the largest optical extent. The equipartition strength of the total field ranges from nearly 10 \mu G at R<2 kpc to 4 \mu G at R=22-24 kpc. As the random field dominates in M 101, wavelength-independent polarization is the main polarization mechanism. We show that energetic events causing HI shells of mean diameter <625 pc could partly be responsible for this. At radii <24 kpc, the random magnetic field depends on the star formation rate per area with a power-law exponent of 0.28+-0.02. The ordered magnetic field is generally aligned with the spiral arms with pitch angles that are about 8{\deg} larger than those of HI filaments.

PSR B0919+06 generally radiates radio pulses in a normal phase range. It has been known for its occasional perplexing abnormal emission events wherein individual pulses come to an earlier phase range for a few tens of periods and then returns to its usual phase. Heretofore, only a few such events have been available for study. We observed PSR B0919+06 for about 30 hours using the Jiamusi 66-m telescope at Jiamusi Deep Space Station at S-band, and detected 92 abnormal emission events. We identify four types of events based on the abrupted or gradual phase-shifting of individual pulses. The abnormal emission events are seen to occur randomly some every 1000 to 3000 periods, and they affect the leading edge of the mean profile by up to 2\% in amplitude. The abnormal emission events are probably related to gradual changes of emission processing in the pulsar magnetosphere.

PSR B0919+06 generally radiates radio pulses in a normal phase range. It has been known for its occasional perplexing abnormal emission events wherein individual pulses come to an earlier phase range for a few tens of periods and then returns to its usual phase. Heretofore, only a few such events have been available for study. We observed PSR B0919+06 for about 30 hours using the Jiamusi 66-m telescope at Jiamusi Deep Space Station at S-band, and detected 92 abnormal emission events. We identify four types of events based on the abrupted or gradual phase-shifting of individual pulses. The abnormal emission events are seen to occur randomly some every 1000 to 3000 periods, and they affect the leading edge of the mean profile by up to 2\% in amplitude. The abnormal emission events are probably related to gradual changes of emission processing in the pulsar magnetosphere.

N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/

We develop two general methods to infer the gravitational potential of a system using steady-state tracers, i.e., tracers with a time-independent phase-space distribution. Combined with the phase-space continuity equation, the time independence implies a universal Orbital Probability Density Function (oPDF) $\mathrm{d} P(\lambda|{\rm orbit})\propto \mathrm{d} t$, where $\lambda$ is the coordinate of the particle along the orbit. The oPDF is equivalent to Jeans theorem, and is the key physical ingredient behind most dynamical modelling of steady-state tracers. In the case of a spherical potential, we develop a likelihood estimator that fits analytical potentials to the system, and a non-parametric method ("phase-mark") that reconstructs the potential profile, both assuming only the oPDF. The methods involve no extra assumptions about the tracer distribution function and can be applied to tracers with any arbitrary distribution of orbits, with possible extension to non-spherical potentials. The methods are tested on Monte Carlo samples of steady-state tracers in dark matter haloes to show that they are unbiased as well as efficient. A fully documented \textsc{C/Python} code implementing our method is freely available at a GitHub repository linked from \url{http://icc.dur.ac.uk/data/#oPDF}.

Using realistic cosmological simulations of Milky Way sized haloes, we study their dynamical state and the accuracy of inferring their mass profiles with steady-state models of dynamical tracers. We use a new method that describes the phase-space distribution of a steady-state tracer population in a spherical potential without any assumption regarding the distribution of their orbits. Applying the method to five haloes from the Aquarius $\Lambda$CDM N-body simulation, we find that dark matter particles are an accurate tracer that enables the halo mass and concentration parameters to be recovered with an accuracy of $5\%$. Assuming a potential profile of the NFW form does not significantly affect the fits in most cases, except for halo A whose density profile differs significantly from the NFW form, leading to a $30\%$ bias in the dynamically fitted parameters. The existence of substructures in the dark matter tracers only affects the fits by $\sim 1\%$. Applying the method to mock stellar haloes generated by a particle-tagging technique, we find the stars are farther from equilibrium than dark matter particles, yielding a systematic bias of $\sim 20\%$ in the inferred mass and concentration parameter. The level of systematic biases obtained from a conventional distribution function fit to stars is comparable to ours, while similar fits to DM tracers are significantly biased in contrast to our fits. In line with previous studies, the mass bias is much reduced near the tracer half-mass radius.

We use weak gravitational lensing to measure mean mass profiles around Locally Brightest Galaxies (LBGs). These are selected from the SDSS/DR7 spectroscopic and photometric catalogues to be brighter than any neighbour projected within 1.0 Mpc and differing in redshift by $<1000$ km/s. Most ($> 83\%$) are expected to be the central galaxies of their dark matter halos. Previous stacking analyses have used this LBG sample to measure mean Sunyaev-Zeldovich flux and mean X-ray luminosity as a function of LBG stellar mass. In both cases, a simulation of the formation of the galaxy population was used to estimate effective halo mass for LBGs of given stellar mass, allowing the derivation of scaling relations between the gas properties of halos and their mass. By comparing results from a variety of simulations to our lensing data, we show that this procedure has significant model dependence reflecting: (i) the failure of any given simulation to reproduce observed galaxy abundances exactly; (ii) a dependence on the cosmology underlying the simulation; and (iii) a dependence on the details of how galaxies populate halos. We use our lensing results to recalibrate the scaling relations, eliminating most of this model dependence and explicitly accounting both for residual modelling uncertainties and for observational uncertainties in the lensing results. The resulting scaling relations link the mean gas properties of dark halos to their mass over an unprecedentedly wide range, $10^{12.5}<M_{500}/ \mathrm{M_\odot}<10^{14.5}$, and should fairly and robustly represent the full halo population.

Diffuse radio emission in galaxy clusters is known to be related to cluster mass and cluster dynamical state. We collect the observed fluxes of radio halos, relics, and mini-halos for a sample of galaxy clusters from the literature, and calculate their radio powers. We then obtain the values of cluster mass or mass proxies from previous observations, and also obtain the various dynamical parameters of these galaxy clusters from optical and X-ray data. The radio powers of relics, halos, and mini-halos are correlated with the cluster masses or mass proxies, as found by previous authors, with the correlations concerning giant radio halos being, in general, the strongest ones. We found that the inclusion of dynamical parameters as the third dimension can significantly reduce the data scatter for the scaling relations, especially for radio halos. We therefore conclude that the substructures in X-ray images of galaxy clusters and the irregular distributions of optical brightness of member galaxies can be used to quantitatively characterize the shock waves and turbulence in the intracluster medium responsible for re-accelerating particles to generate the observed diffuse radio emission. The power of radio halos and relics is correlated with cluster mass proxies and dynamical parameters in the form of a fundamental plane.

The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $\Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40\% underestimate to a 5\% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $\beta$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.

Fast radio bursts show large dispersion measures, much larger than the Galactic dispersion measure foreground. Therefore,they evidently have an extragalactic origin. We investigate possible contributions to the dispersion measure from host galaxies. We simulate the spatial distribution of fast radio bursts and calculate the dispersion measures along the sightlines from fast radio bursts to the edge of host galaxies by using the scaled NE2001 model for thermal electron density distributions. We find that contributions to the dispersion measure of fast radio bursts from the host galaxy follow a skew Gaussian distribution. The peak and the width at half maximum of the dispersion measure distribution increase with the inclination angle of a spiral galaxy, to large values when the inclination angle is over $70\degr$. The largest dispersion measure produced by an edge-on spiral galaxy can reach a few thousand pc~cm$^{-3}$, while the dispersion measures from dwarf galaxies and elliptical galaxies have a maximum of only a few tens of pc~cm$^{-3}$. Notice, however, that additional dispersion measures of tens to hundreds of pc~cm$^{-3}$ can be produced by high density clumps in host galaxies. Simulations that include dispersion measure contributions from the Large Magellanic Cloud and the Andromeda Galaxy are shown as examples to demonstrate how to extract the dispersion measure from the intergalactic medium.

Spiral arms shown by different components may not be spatially coincident, which can constrain formation mechanisms of spiral structure in a galaxy. We reassess the spiral arm tangency directions in the Milky Way through identifying the bump features in the longitude plots of survey data for infrared stars, radio recombination lines (RRLs), star formation sites, CO, high density regions in clouds, and HI. The bump peaks are taken as indications for arm tangencies, which are close to the real density peaks near the spiral arm tangency point but often have $\sim$ 1$^\circ$ offset to the interior of spiral arms. The arm tangencies identified from the longitudes plots for RRLs, HII regions, methanol masers, CO, high density gas regions, and HI gas appear nearly the same Galactic longitude, and therefore there is no obvious offset for spiral arms traced by different gas components. However, we find obvious displacements of 1.3$^\circ-$ 5.8$^\circ$ between gaseous bump peaks from the directions of the maximum density of old stars near the tangencies of the Scutum-Centaurus Arm, the northern part of the Near 3 kpc Arm, and maybe also the Sagittarius Arm. The offsets between the density peaks of gas and old stars for spiral arms are comparable with the arm widths, which is consistent with expectations for quasi-stationary density wave in our Galaxy.

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($\Lambda$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $\Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $\Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $\Lambda$CDM predictions. In fact, ~10% of $\Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.

Accurately determining the mass of galaxy clusters is fundamental for many studies on cosmology and galaxy evolution. We collect and rescale the cluster masses of 1191 clusters of 0.05<z<0.75 estimated by X-ray or Sunyaev-Zeldovich measurements and use them to calibrate the optical mass proxy. The total r-band luminosity (in units of L^{\ast}) of these clusters are obtained by using spectroscopic and photometric data of the Sloan Digital Sky Survey (SDSS). We find that the correlation between the cluster mass M_{500} and total r-band luminosity L_{500} significantly evolves with redshift. After correcting for the evolution, we define a new cluster richness R_{L\ast,500}=L_{500}*E(z)^{1.40} as the optical mass proxy. By using this newly defined richness and the recently released SDSS DR12 spectroscopic data, we update the WHL12 cluster catalog and identify 25,419 new rich clusters at high redshift. In the SDSS spectroscopic survey region, about 89% of galaxy clusters have spectroscopic redshifts. The mass can be estimated with a scatter of 0.17 dex for the clusters in the updated catalog.

We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

Superbubbles and supershells are the channels for transferring mass and energy from the Galactic disk to the halo. Magnetic fields are believed to play a vital role in their evolution. We study the radio continuum and polarized emission properties of the W4 superbubble to determine its magnetic field strength. New sensitive radio continuum observations were made at 6 cm, 11 cm, and 21 cm. The total intensity measurements were used to derive the radio spectrum of the W4 superbubble. The linear polarization data were analysed to determine the magnetic field properties within the bubble shells. The observations show a multi-shell structure of the W4 superbubble. A flat radio continuum spectrum that stems from optically thin thermal emission is derived from 1.4 GHz to 4.8 GHz. By fitting a passive Faraday screen model and considering the filling factor fne , we obtain the thermal electron density ne = 1.0/\sqrt{fne} (\pm5%) cm^-3 and the strength of the line-of-sight component of the magnetic field B// = -5.0/\sqrt{fne} (\pm10%) {\mu}G (i.e. pointing away from us) within the western shell of the W4 superbubble. When the known tilted geometry of the W4 superbubble is considered, the total magnetic field Btot in its western shell is greater than 12 {\mu}G. The electron density and the magnetic field are lower and weaker in the high-latitude parts of the superbubble. The rotation measure is found to be positive in the eastern shell but negative in the western shell of the W4 superbubble, which is consistent with the case that the magnetic field in the Perseus arm is lifted up from the plane towards high latitudes. The magnetic field strength and the electron density we derived for the W4 superbubble are important parameters for evolution models of superbubbles breaking out of the Galactic plane.

The origin of the 6.67 hr period X-ray source, 1E161348-5055, in the young supernova remnant RCW 103 is puzzling. We propose that it may be the descendant of a Thorne-Zytkow Object (TZO). A TZO may at its formation have a rapidly spinning neutron star as a core, and a slowly rotating envelope. We found that the core could be braked quickly to an extremely long spin period by the coupling between its magnetic field and the envelope, and that the envelope could be disrupted by some powerful bursts or exhausted via stellar wind. If the envelope is disrupted after the core has spun down, the core will become an extremely long-period compact object, with a slow proper motion speed, surrounded by a supernova-remnant-like shell. These features all agree with the observations of 1E161348-5055. TZOs are expected to have produced extraordinary high abundances of lithium and rapid proton process elements that would remain in the remnants and could be used to test this scenario.

Frequency dependence of pulsar linear polarization is investigated by simulations of emission and propagation processes. Linearly polarized waves are generated through curvature radiation by relativistic particles streaming along curved magnetic field lines, which have ordinary mode (O-mode) and extra-ordinary mode (X-mode) components. As emitted waves propagate outwards, two mode components are separated due to re- fraction of the O mode, and their polarization states are also modified. According to the radius to frequency mapping, low frequency emission is generated from higher magnetosphere, where significant rotation effect leads the X and O modes to be sepa- rated. Hence, the low frequency radiation has a large fraction of linear polarization. As the frequency increases, emission is generated from lower heights, where the rotation effect becomes weaker and the distribution regions of two modes are more overlapped. Hence, more significant depolarization appears for emission at higher frequencies. In addition, refraction effect of the O mode becomes serious in very deep magnetosphere, which bends the O mode emission towards outer parts of a pulsar beam and also causes the separation of mode distribution regions and hence the fractional linear polariza- tion increasing with frequency. If emission of different frequencies is generated from a region of the same height, serious O mode refraction can result in the decrease of both profile width and fractional linear polarization. The observed frequency dependence of linear polarization for some pulsars can be naturally explained within the scope of our scenario.

Luminosity function of cluster galaxies provides a fundamental constraint on galaxy evolution in cluster environments. By using the bright member galaxies of a large sample of rich clusters identified from Sloan Digital Sky Survey, we obtain the bright end of composite luminosity functions of cluster galaxies, and study their dependence on cluster dynamical state. After a redshift-evolution correction of absolute magnitude, the luminosity function of member galaxies can be well fitted by a Schechter function when the brightest cluster galaxies (BCGs) are excluded. The absolute magnitudes of BCGs follow a Gaussian function with a characteristic width of about 0.36 mag. We find that the luminosity function of galaxies in more relaxed clusters has a fainter characteristic absolute magnitude (M_{\ast}), and these clusters have fewer bright non-BCG member galaxies but a brighter BCG. Our results suggest the co-evolution of galaxy population with cluster dynamical state and somewhat support the hierarchical formation scenario of the BCGs.

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $\Lambda$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $\sigma_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

The spiral structure of the Milky Way is not yet well determined. The keys to understanding this structure are to increase the number of reliable spiral tracers and to determine their distances as accurately as possible. HII regions, giant molecular clouds (GMCs), and 6.7-GHz methanol masers are closely related to high mass star formation, and hence they are excellent spiral tracers. We update the catalogs of Galactic HII regions, GMCs, and 6.7-GHz methanol masers, and then outline the spiral structure of the Milky Way. We collected data for more than 2500 known HII regions, 1300 GMCs, and 900 6.7-GHz methanol masers. If the photometric or trigonometric distance was not yet available, we determined the kinematic distance using a Galaxy rotation curve with the current IAU standard, $R_0$ = 8.5 kpc and $\Theta_0$ = 220 km s$^{-1}$, and the most recent updated values of $R_0$ = 8.3 kpc and $\Theta_0$ = 239 km s$^{-1}$, after we modified the velocities of tracers with the adopted solar motions. With the weight factors based on the excitation parameters of HII regions or the masses of GMCs, we get the distributions of these spiral tracers. The distribution of tracers shows at least four segments of arms in the first Galactic quadrant, and three segments in the fourth quadrant. The Perseus Arm and the Local Arm are also delineated by many bright HII regions. The arm segments traced by massive star forming regions and GMCs are able to match the HI arms in the outer Galaxy. We found that the models of three-arm and four-arm logarithmic spirals are able to connect most spiral tracers. A model of polynomial-logarithmic spirals is also proposed, which not only delineates the tracer distribution, but also matches the observed tangential directions.

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

We compiled a catalog of Faraday rotation measures (RMs) for 4553 extragalactic radio point sources ublished in literature. These RMs were derived from multi-frequency polarization observations. The RM data are compared to those in the NRAO VLA Sky Survey (NVSS) RM catalog. We reveal a systematic uncertainty of about $10.0 \pm 1.5$\,rad~m$^{-2}$ in the NVSS RM catalog. The Galactic foreground RM is calculated through a weighted averaging method by using the compiled RM catalog together with the NVSS RM catalog, with careful consideration of uncertainties in the RM data. The data from the catalog and the interface for the Galactic foreground RM calculations are publicly available on the webpage: http://zmtt.bao.ac.cn/RM/.

The survey data of Wide-field Infrared Survey Explorer (WISE) provide an opportunity for the identification of galaxy clusters. We present an efficient method for detecting galaxy clusters by combining the WISE data with SuperCOSMOS and 2MASS data. After performing star-galaxy separation, we calculate the number of companion galaxies around the galaxies with photometric redshifts previously estimated by the SuperCOSMOS, 2MASS and WISE data. A scaled richness is then defined to identify clusters. From a sky area of 275 deg^2 coincident with Sloan Digital Sky Survey Stripe 82 region, we identify 302 clusters in the redshift range of 0.1<z<0.35, 247 (82%) of which are previously known SDSS clusters. The results indicate that our method is efficient for identifying galaxy clusters by using the all sky data of the SuperCOSMOS, 2MASS and WISE.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

We present the first cross-correlation measurement between Sloan Digital Sky Survey (SDSS) Type 1 quasars and the cosmic infrared background (CIB) measured by Herschel. The distribution of the quasars at 0.15<z<3.5 covers the redshift range where we expect most of the CIB to originate. We detect the sub-mm emission of the quasars, which dominates on small scales, as well as correlated emission from dusty star-forming galaxies (DSFGs) dominant on larger scales. The mean sub-mm flux densities of the DR7 quasars (median redshift <z>=1.4) is $11.1^{+1.6}_{-1.4}$, $7.1^{+1.6}_{-1.3}$ and $3.6^{+1.4}_{-1.0}$ mJy at 250, 350 and 500 microns, respectively, while the mean sub-mm flux densities of the DR9 quasars (<z>=2.5) is $5.7^{+0.7}_{-0.6}$, $5.0^{+0.8}_{-0.7}$ and $1.8^{+0.5}_{-0.4}$ mJy. We find that the correlated sub-mm emission includes both the emission from satellite DSFGs in the same halo as the central quasar and the emission from DSFGs in separate halos (correlated with the quasar-hosting halo). The amplitude of the one-halo term is ~10 times smaller than the sub-mm emission of the quasars, implying the the satellites have a lower star-formation rate than the quasars. The satellite fraction for the DR7 quasars is $0.008^{+0.008}_{-0.005}$ and the host halo mass scale for the central and satellite quasars is $10^{12.36\pm0.87}$ M$_{\odot}$ and $10^{13.60\pm0.38}$ M$_{\odot}$, respectively. The satellite fraction of the DR9 quasars is $0.065^{+0.021}_{-0.031}$ and the host halo mass scale for the central and satellite quasars is $10^{12.29\pm0.62}$ M$_{\odot}$ and $10^{12.82\pm0.39}$ M$_{\odot}$, respectively. Thus, the typical halo environment of the SDSS Type 1 quasars is found to be similar to that of DSFGs, which supports the generally accepted view that dusty starburst and quasar activity are evolutionarily linked.

Large-scale radio continuum surveys have good potential for discovering new Galactic supernova remnants (SNRs). Surveys of the Galactic plane are often limited in the Galactic latitude of |b| ~ 5 degree. SNRs at high latitudes, such as the Cygnus Loop or CTA~1, cannot be detected by surveys in such limited latitudes. Using the available Urumqi 6 cm Galactic plane survey data, together with the maps from the extended ongoing 6 cm medium latitude survey, we wish to discover new SNRs in a large sky area. We searched for shell-like structures and calculated radio spectra using the Urumqi 6 cm, Effelsberg 11 cm, and 21 cm survey data. Radio polarized emission and evidence in other wavelengths are also examined for the characteristics of SNRs. We discover an enclosed oval-shaped object G150.3+4.5 in the 6 cm survey map. It is about 2.5 degree wide and 3 degree high. Parts of the shell structures can be identified well in the 11 cm, 21 cm, and 73.5 cm observations. The Effelsberg 21 cm total intensity image resembles most of the structures of G150.3+4.5 seen at 6 cm, but the loop is not closed in the northwest. High resolution images at 21 cm and 73.5 cm from the Canadian Galactic Plane Survey confirm the extended emission from the eastern and western shells of G150.3+4.5. We calculated the radio continuum spectral indices of the eastern and western shells, which are $\beta \sim -2.4$ and $\beta \sim -2.7$ between 6 cm and 21 cm, respectively. The shell-like structures and their non-thermal nature strongly suggest that G150.3+4.5 is a shell-type SNR. For other objects in the field of view, G151.4+3.0 and G151.2+2.6, we confirm that the shell-like structure G151.4+3.0 very likely has a SNR origin, while the circular-shaped G151.2+2.6 is an HII region with a flat radio spectrum, associated with optical filamentary structure, H$\alpha$, and infrared emission.

Context. The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is a variability survey of the Milky Way bulge and an adjacent section of the disk carried out from 2010 on ESO Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV will eventually deliver a deep near-IR atlas with photometry and positions in five passbands (ZYJHK_S) and a catalogue of 1-10 million variable point sources - mostly unknown - which require classifications. Aims. The main goal of the VVV Templates Project, that we introduce in this work, is to develop and test the machine-learning algorithms for the automated classification of the VVV light-curves. As VVV is the first massive, multi-epoch survey of stellar variability in the near-infrared, the template light-curves that are required for training the classification algorithms are not available. In the first paper of the series we describe the construction of this comprehensive database of infrared stellar variability. Methods. First we performed a systematic search in the literature and public data archives, second, we coordinated a worldwide observational campaign, and third we exploited the VVV variability database itself on (optically) well-known stars to gather high-quality infrared light-curves of several hundreds of variable stars. Results. We have now collected a significant (and still increasing) number of infrared template light-curves. This database will be used as a training-set for the machine-learning algorithms that will automatically classify the light-curves produced by VVV. The results of such an automated classification will be covered in forthcoming papers of the series.

We have observed 79 supernova remnants (SNRs) with the Urumqi 25m telescope at 6cm during the Sino-German 6cm polarization survey of the Galactic plane. We measured flux densities of SNRs at 6cm, some of which are the first ever measured or the measurements at the highest frequency, so that we can determine or improve spectra of SNRs. Our observations have ruled out spectral breaks or spectral flattening that were suggested for a few SNRs, and confirmed the spectral break of S147. By combining our 6cm maps with 11cm and 21cm maps from the Effelsberg 100 m telescope, we calculated the spectral index maps of several large SNRs. For many remnants we obtained for the first time polarization images, which show the intrinsic magnetic field structures at 6cm. We disapproved three objects as being SNRs, OA184, G192.8-1.1 and G16.8-1.1, which show a thermal spectrum and no polarization. We have discovered two large supernova remnants, G178.2-4.2 and G25.1-2.3., in the survey maps.

We study how to measure the galaxy merger rate from the observed close pair count. Using a high-resolution N-body/SPH cosmological simulation, we find an accurate scaling relation between galaxy pair counts and merger rates down to a stellar mass ratio of about 1:30. The relation explicitly accounts for the dependence on redshift (or time), on pair separation, and on mass of the two galaxies in a pair. With this relation, one can easily obtain the mean merger timescale for a close pair of galaxies. The use of virial masses, instead of stellar masses, is motivated by the fact that the dynamical friction time scale is mainly determined by the dark matter surrounding central and satellite galaxies. This fact can also minimize the error induced by uncertainties in modeling star formation in the simulation. Since the virial mass can be read from the well-established relation between the virial masses and the stellar masses in observation, our scaling relation can be easily applied to observations to obtain the merger rate and merger time scale. For major merger pairs (1:1-1:4) of galaxies above a stellar mass of 4*10^10 M_sun/h at z=0.1, it takes about 0.31 Gyr to merge for pairs within a projected distance of 20 kpc/h with stellar mass ratio of 1:1, while the time taken goes up to 1.6 Gyr for mergers with stellar mass ratio of 1:4. Our results indicate that a single timescale usually used in literature is not accurate to describe mergers with the stellar mass ratio spanning even a narrow range from 1:1 to 1:4.

We obtained rotation measures of 2642 quasars by cross-identification of the most updated quasar catalog and rotation measure catalog. After discounting the foreground Galactic Faraday rotation of the Milky Way, we get the residual rotation measure (RRM) of these quasars. We carefully discarded the effects from measurement and systematical uncertainties of RRMs as well as large RRMs from outliers, and get marginal evidence for the redshift evolution of real dispersion of RRMs which steady increases to 10 rad m$^{-2}$ from $z=0$ to $z\sim1$ and is saturated around the value at higher redshifts. The ionized clouds in the form of galaxy, galaxy clusters or cosmological filaments could produce the observed RRM evolutions with different dispersion width. However current data sets can not constrain the contributions from galaxy halos and cosmic webs. Future RM measurements for a large sample of quasars with high precision are desired to disentangle these different contributions.

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

The propagation of polarized emission in pulsar magnetosphere is investigated in this paper. The polarized waves are generated through curvature radiation from the relativistic particles streaming along curved magnetic field lines and co-rotating with the pulsar magnetosphere. Within the 1/{\deg} emission cone, the waves can be divided into two natural wave mode components, the ordinary (O) mode and the extraord nary (X) mode, with comparable intensities. Both components propagate separately in magnetosphere, and are aligned within the cone by adiabatic walking. The refraction of O-mode makes the two components separated and incoherent. The detectable emission at a given height and a given rotation phase consists of incoherent X-mode and O-mode components coming from discrete emission regions. For four particle-density models in the form of uniformity, cone, core and patches, we calculate the intensities for each mode numerically within the entire pulsar beam. If the co-rotation of relativistic particles with magnetosphere is not considered, the intensity distributions for the X-mode and O-mode components are quite similar within the pulsar beam, which causes serious depolarization. However, if the co-rotation of relativistic particles is considered, the intensity distributions of the two modes are very different, and the net polarization of out-coming emission should be significant. Our numerical results are compared with observations, and can naturally explain the orthogonal polarization modes of some pulsars. Strong linear polarizations of some parts of pulsar profile can be reproduced by curvature radiation and subsequent propagation effect.

The very small braking index of PSR J1734-3333, $n=0.9\pm0.2$, challenges the current theories of braking mechanisms in pulsars. We present a possible interpretation that this pulsar is surrounded by a fall-back disk and braked by it. A modified braking torque is proposed based on the competition between the magnetic energy density of a pulsar and the kinetic energy density of a fall-back disk. With this torque, a self-similar disk can fit all the observed parameters of PSR J1734-3333 with natural initial parameters. In this regime, the star will evolve to the region having anomalous X-ray pulsars and soft gamma repeaters in the $P-\dot{P}$ diagram in about 20000 years and stay there for a very long time. The mass of the disk around PSR J1734-3333 in our model is about $10M_{\oplus}$, similar to the observed mass of the disk around AXP 4U 0142+61.

We present a Revised IRAS-FSC Redshift Catalogue (RIFSCz) of 60,303 galaxies selected at 60 microns from the IRAS Faint Source Catalogue (FSC). This revision merges in data from the WISE All-Sky Data Release, the Tenth SDSS Data Release (DR10), the GALEX All-Sky Survey Source Catalog (GASC), the 2MASS Redshift Survey (2MRS) and the Planck Catalogue of Compact Sources (PCCS). The RIFSCz consists of accurate position, ultra-violet (UV), optical, near-, mid- and far-infrared, sub-millimetre (sub-mm) and/or radio identifications, spectroscopic redshift (if available) or photometric redshift (if possible), predicted far-infrared and sub-mm fluxes ranging from 12 to 1380 microns based upon the best-fit infrared template. We also provide stellar masses, star-formation rates and dust masses derived from the optical and infrared template fits, where possible. 56 of the galaxies in the RIFSCz have spectroscopic redshifts and a further 26 have photometric redshifts obtained through the template-fitting method. At S60 > 0.36 Jy, the 90% completeness limit of the FSC, 93 of the sources in the RIFSCz have either spectroscopic or photometric redshifts. An interesting subset of the catalogue is the sources detected by Planck at sub-mm wavelengths. 1200 sources have a detection at better than 5 sigma in at least one Planck band and a further 1186 sources have detections at 3-5 sigma in at least one Planck band.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

Information of small-diameter sources is extracted from the Sino-German 6cm polarisation survey of the Galactic plane carried out with the Urumqi 25-m telescope. We performed two-dimensional elliptical Gaussian fits to the 6cm maps to obtain a list of sources with total-intensity and polarised flux densities. The source list contains 3832 sources with a fitted diameter smaller than 16 arcmin and a peak flux density exceeding 30 mJy, so about 5 times the rms noise, of the total-intensity data. The cumulative source count indicates completeness for flux densities exceeding about 60 mJy. We identify 125 linearly polarised sources at 6cm with a peak polarisation flux density greater than 10 mJy, so about 3 times the rms noise, of the polarised-intensity data. Despite lacking compact steep spectrum sources, the 6cm catalogue lists about 20 percent more sources than the Effelsberg 21cm source catalogue at the same angular resolution and for the same area. Most of the faint 6cm sources must have a flat spectrum and are either HII regions or extragalactic. When compared with the Green Bank 6cm (GB6) catalogue, we obtain higher flux densities for a number of extended sources with complex structures. Polarised 6cm sources density are uniformly distributed in Galactic latitude. Their number density decreases towards the inner Galaxy. More than 80 percent of the polarised sources are most likely extragalactic. With a few exceptions, the sources have a higher percentage polarisation at 6cm than at 21cm. Depolarisation seems to occur mostly within the sources with a minor contribution from the Galactic foreground emission.

Pulsar nulling is a phenomenon of sudden cessation of pulse emission for a number of periods. The nulling fraction was often used to characterize the phenomenon. We propose a new method to analyse pulsar nulling phenomenon, by involving two key parameters, the nulling degree, $\chi$, which is defined as the angle in a rectangular coordinates for the numbers of nulling periods and bursting periods, and the nulling scale, $ N $, which is defined as the effective length of the consecutive nulling periods and bursting periods. The nulling degree $\chi$ can be calculated by $\tan \chi = N_{\rm nulling} / N_{\rm bursting} $ and the mean is related to the nulling fraction, while the nulling scale, $ N $, is also a newly defined fundamental parameter which indicates how often the nulling occurs. We determined the distributions of $\chi$ and $ N $ for 10 pulsars by using the data in literature. We found that the nulling degree $\chi$ indicates the relative length of nulling to that of bursting, and the nulling scale $ N $ is found to be related to the derivative of rotation frequency and hence the loss rate of rotational energy of pulsars. Their deviations reflect the randomness of the nulling process.

Clusters of galaxies have a huge mass which can act as gravitational lenses. Galaxies behind clusters can be distorted to form arcs in images by the lenses. Herein a search was done for giant lensed arcs by galaxy clusters using the SDSS data. By visually inspecting SDSS images of newly identified clusters in the SDSS DR8 and Stripe 82 data, we discover 8 strong lensing clusters together with additional 3 probable and 6 possible cases. The lensed arcs show bluer colors than the member galaxies of clusters. The masses and optical luminosities of galaxy clusters interior to the arcs are calculated, and the mass-to-light ratios are found to be in the range of a few tens of M_Solar/L_Solar, consistent with the distribution of previously known lensing clusters.

Dynamical state of galaxy clusters is closely related to their observational properties in X-ray, optical and radio wavelengths. We develop a method to diagnose the substructure and dynamical state of galaxy clusters by using photometric data of Sloan Digital Sky Survey (SDSS). To trace mass distribution, the brightness distribution of member galaxies is smoothed by using a Gaussian kernel with a weight of their optical luminosities. After deriving the asymmetry, the ridge flatness and the normalized deviation of the smoothed optical map, we define a relaxation parameter, Gamma, to quantify dynamical state of clusters. This method is applied to a test sample of 98 clusters of 0.05<z\lesssim0.42 collected from literature with known dynamical states and can recognize dynamical state for relaxed (Gamma\ge0) and unrelaxed (Gamma<0) clusters with a success rate of 94%. We then calculate relaxation parameters of 2092 rich clusters previously identified from the SDSS, of which 28% clusters are dynamically relaxed with Gamma\ge0. We find that the dominance and absolute magnitude of the brightest cluster galaxy closely correlate with dynamical state of clusters. The emission power of radio halos is quantitatively related to cluster dynamical state, beside the known dependence on the X-ray luminosity.

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $\mu$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $\mu$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016.

Radio continuum observations detect non-thermal synchrotron and thermal bremsstrahlung radiation. Separation of the two different emission components is crucial to study the properties of diffuse interstellar medium. The Cygnus X region is one of the most complex areas in the radio sky which contains a number of massive stars and HII regions on the diffuse thermal and non-thermal background. More supernova remnants are expected to be discovered. We aim to develop a method which can properly separate the non-thermal and thermal radio continuum emission and apply it to the Cygnus X region. The result can be used to study the properties of different emission components and search for new supernova remnants in the complex. Multi-frequency radio continuum data from large-scale surveys are used to develop a new component separation method. Spectral analysis is done pixel by pixel for the non-thermal synchrotron emission with a realistic spectral index distribution and a fixed spectral index of beta = -2.1 for the thermal bremsstrahlung emission. With the new method, we separate the non-thermal and thermal components of the Cygnus X region at an angular resolution of 9.5arcmin. The thermal emission component is found to comprise 75% of the total continuum emission at 6cm. Thermal diffuse emission, rather than the discrete HII regions, is found to be the major contributor to the entire thermal budget. A smooth non-thermal emission background of 100 mK Tb is found. We successfully make the large-extent known supernova remnants and the HII regions embedded in the complex standing out, but no new large SNRs brighter than Sigma_1GHz = 3.7 x 10^-21 W m^-2 Hz^-1 sr^-1 are found.

Magnetic fields in our Galaxy and nearby galaxies have been revealed by starlight polarization, polarized emission from dust grains and clouds at millimeter and submillimeter wavelength, the Zeeman effect of spectral lines or maser lines from clouds or clumps, diffuse radio synchrotron emission from relativistic electrons in interstellar magnetic fields, and the Faraday rotation of background radio sources as well as pulsars for our Milky Way. It is easy to get a global structure for magnetic fields in nearby galaxies, while we have observed many details of magnetic fields in our Milky Way, especially by using pulsar rotation measure data. In general, magnetic fields in spiral galaxies probably have a large-scale structure. The fields follow the spiral arms with or without the field direction reversals. In the halo of spiral galaxies magnetic fields exist and probably also have a large-scale structure as toroidal and poloidal fields, but seem to be slightly weaker than those in the disk. In the central region of some galaxies, poloidal fields have been detected as vertical components. Magnetic field directions in galaxies seem to have been preserved during cloud formation and star formation, from large-scale diffuse interstellar medium to molecular clouds and then to the cloud cores in star formation regions or clumps for the maser spots. Magnetic fields in galaxies are passive to dynamics.

The Five-hundred-metre Aperture Spherical Telescope (FAST) uses adaptive spherical panels to achieve a huge collecting area for radio waves. In this paper, we try to explore the optimal parameters for the curvature radius of spherical panels and the focal distance by comparison of the calculated beam patterns. We show that to get the best beam shape and maximum gain, the optimal curvature radius of panels is around 300 m, and a small shift in the focal distance of a few cm is needed. The aperture efficiency can be improved by ~10% at 3 GHz by this small shift. We also try to optimise the panel positioning for the best beam, and find that panel shifts of a few mm can improve the beam pattern by a similar extent. Our results indicate that accurate control of the feed and panel positions to the mm level is very crucial for the stability of FAST's observational performance.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20%$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40%$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18%$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.

The Sino-German 6cm polarization survey has mapped in total intensity and polarization intensity over an area of approximately 2200 square degrees in the Galactic disk. This survey provides an opportunity to search for Galactic supernova remnants (SNRs) that were previously unknown. We discovered the new SNRs G178.2-4.2 and G25.1-2.3 which have non-thermal spectra, using the 6cm data together with the observations with the Effelsberg telescope at 11 cm and 21 cm. Both G178.2-4.2 and G25.1-2.3 are faint and have an apparent diameter greater than 1deg. G178.2-4.2 shows a polarized shell. HI data suggest that G25.1-2.3 might have a distance of about 3 kpc. The 6cm survey data were also very important to identify two other new SNRs, G152.4-2.1 and G190.9-2.2.

The spiral structure of our Milky Way has not yet been well outlined. HII regions, giant molecular clouds (GMCs) and 6.7-GHz methanol masers are primary tracers for spiral arms. We collect and update the database of these tracers which has been used in Hou, Han & Shi (2009) for mapping the spiral structure.

Beam radii for cone-dominant pulsars follow a power-law relation with frequency, $\vartheta = (\nu/\nu_0)^k+\vartheta_0$, which has not yet well explained in previous works. We study this frequency dependence of beam radius (FDB) for cone-dominant pulsars by using the curvature radiation mechanism. Considering various density and energy distributions of particles in the pulsar open field line region, we numerically simulate the emission intensity distribution across emission height and rotation phase, and get integrated profiles at different frequencies and obtain the FDB curves. For the density model of a conal-like distribution, the simulated profiles always shrink to one component at high frequencies. In the density model with two separated density patches, the profiles always have two distinct components, and the power-law indices $k$ are found to be in the range from -0.1 to -2.5, consistent with observational results. Energy distributions of streaming particles have significant influence on the frequency-dependence behavior. Radial energy decay of particles are necessary to get proper $\vartheta_0$ in models. We conclude that by using the curvature radiation mechanism, the observed frequency dependence of beam radius for the cone-dominant pulsars can only be explained by the emission model of particles in two density patches with a Gaussian energy distribution and a radial energy loss.

We identified 132,684 clusters in the redshift range of 0.05<z<0.8 from SDSS DR8. The spectroscopic redshifts of 52,683 clusters have been included in the catalog using SDSS DR9 data. We found that BCGs are more luminous in richer clusters and at higher redshifts.

Cold dark matter models predict the existence of a large number of substructures within dark matter halos. If the cold dark matter consists of weakly interacting massive particles, their annihilation within these substructures could lead to diffuse GeV emission that would dominate over the annihilation signal of the host halo. In this work we search for GeV emission from three nearby galaxy clusters: Coma, Virgo and Fornax. We first remove known extragalactic and galactic diffuse gamma-ray backgrounds and point sources from the Fermi 2-year catalog and find a significant residual diffuse emission in all three clusters. We then investigate whether this emission is due to (i) unresolved point sources; (ii) dark matter annihilation; or (iii) cosmic rays (CR). Using 45 months of Fermi-LAT data we detect several new point sources (not present in the Fermi 2-year point source catalogue) which contaminate the signal previously analyzed by Han et al.(arxiv:1201.1003). Including these and accounting for the effects of undetected point sources, we find no significant detection of extended emission from the three clusters studied. Instead, we determine upper limits on emission due to dark matter annihilation and cosmic rays. For Fornax and Virgo the limits on CR emission are consistent with theoretical models, but for Coma the upper limit is a factor of 2 below the theoretical expectation. Allowing for systematic uncertainties associated with the treatment of CR, the upper limits on the cross section for dark matter annihilation from our clusters are more stringent than those from analyses of dwarf galaxies in the Milky Way. We rule out the thermal cross section for supersymmetric dark matter particles for masses as large as 100 GeV (depending on the annihilation channel).

The Origem Loop in the Galactic anticentre was discovered in 1970s and suggested to be a large supernova remnant. It was argued later to be a chance superposition of unrelated radio sources. We attempt to understand the properties of the Origem Loop. Available multi-frequency radio data were used for the determination of radio spectra of different parts of the Origem Loop and the polarization properties of the loop. Newly available sensitive observations show that the Origem Loop is a loop of more than 6 deg in diameter. It consists of a large non-thermal arc in the north, which we call the Origem Arc, and several known thermal H II regions in the south. Polarized radio emission associated with the arc was detected at 6 cm, revealing tangential magnetic fields. The arc has a brightness temperature spectral index of \beta = -2.70, indicating its non-thermal nature as a supernova remnant. We estimate the distance to the Origem Arc to be about 1.7 kpc, similar to those of some H II regions in the southern part of the loop. The Origem Loop is a visible loop in the sky, which consists of a supernova remnant arc in the north and H II regions in the south.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

The concept of a "magnetar" was proposed mainly because of two factors. First, the X-ray luminosity of Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) is larger than the rotational energy loss rate, and second, the magnetic field strength calculated from "normal method" is super strong. It is proposed that the radiation energy of magnetar comes from its magnetic fields. Here it is argued that the magnetic field strength calculated through the normal method is incorrect when X-ray luminosity is larger than rotational energy loss rate, because the wind braking is not taken into account. Besides, the "anti-magnetar" and some other X-ray and radio observations are difficult to understand with a magnetar model. Instead of the magnetar, we propose a "quarctar", which is a crusted quark star in an accretion disk, to explain the observations. In this model, the persistent X-ray emission, burst luminosity, spectrum of AXPs and SGRs can be understood naturally. The radio-emitting AXPs, which are challenging the magnetar, can also be explained by the quarctar model.

In this invited talk, I first discuss the advantages and disadvantages of many probes for the magnetic fields of the Milky Way. I conclude that pulsars are the best probes for the magnetic structure in our Galaxy, because magnetic field strength and directions can be derived from their dispersion measures (DMs) and rotation measures (RMs). Using the pulsars as probes, magnetic field structures in the Galactic disk, especially the field reversals between the arms and interarm regions, can be well revealed from the distribution of RM data. The field strengths on large scales and small scales can be derived from RM and DM data. RMs of extragalactic radio sources can be used as the indication of magnetic field directions in the spiral tangential regions, and can be used as probes for the magnetic fields in the regions farther away than pulsars when their median RMs are compared with pulsar RMs.

We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dark matter subhaloes in simulations based on the merger hierarchy of dark matter haloes. Application of this code to a recent benchmark test of finding subhaloes demonstrates that HBT stands as one of the best codes to trace the evolutionary history of subhaloes. The success of the code lies in its careful treatment of the complex physical processes associated with the evolution of subhaloes and in its robust unbinding algorithm with an adaptive source subhalo management. We keep a full record of the merger hierarchy of haloes and subhaloes, and allow growth of satellite subhaloes through accretion from its "satellite-of-satellites", hence allowing mergers among satellites. Local accretion of background mass is omitted, while rebinding of stripped mass is allowed. The justification of these treatments is provided by case studies of the lives of individual subhaloes and by the success in finding the complete subhalo catalogue. We compare our result to other popular subhalo finders and show that HBT is able to well resolve subhaloes in high density environment and keep strict physical track of subhaloes' merger history. This code is fully parallelized and freely available upon request to the authors.

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

We report evidence for extended gamma-ray emission from the Virgo, Fornax and Coma clusters based on a maximum-likelihood analysis of the 3-year Fermi-LAT data. For all three clusters, excess emission is observed within three degrees of the center, peaking at the GeV scale. This emission cannot be accounted for by known Fermi sources or by the galactic and extragalactic backgrounds. If interpreted as annihilation emission from supersymmetric dark matter (DM) particles, the data prefer models with a particle mass in the range 20-60 GeV annihilating into the b-bbar channel, or 2-10 GeV and >1 TeV annihilating into mu-mu final states. Our results are consistent with those obtained by Hooper and Linden from a recent analysis of Fermi-LAT data in the region of the Galactic Centre. An extended DM annihilation profile dominated by emission from substructures is preferred over a simple point source model. The significance of DM detection is 4.4 sigma in Virgo and lower in the other two clusters. We also consider the possibility that the excess emission arises from cosmic ray (CR) induced gamma-rays, and infer a CR level within a factor of three of that expected from analytical models. However, the significance of a CR component is lower than the significance of a DM component, and there is no need for such a CR component in the presence of a DM component in the preferred DM mass range. We also set flux and cross-section upper limits for DM annihilation into the b-bbar and mu-mu channels in all three clusters.

Using the photometric redshifts of galaxies from the Sloan Digital Sky Survey III (SDSS-III), we identify 132,684 clusters in the redshift range of 0.05<z<0.8. Monte Carlo simulations show that the false detection rate is less than 6% for the whole sample. The completeness is more than 95% for clusters with a mass of M_{200}>1.0*10^{14} M_{\odot} in the redshift range of 0.05<z<0.42, while clusters of z>0.42 are less complete and have a biased smaller richness than the real one due to incompleteness of member galaxies. We compare our sample with other cluster samples, and find that more than 90% of previously known rich clusters of 0.05<z<0.42 are matched with clusters in our sample. Richer clusters tend to have more luminous brightest cluster galaxies (BCGs). Correlating with X-ray and the Planck data, we show that the cluster richness is closely related to the X-ray luminosity, temperature and Sunyaev-Zel'dovich measurements. Comparison of the BCGs with the SDSS luminous red galaxy (LRG) sample shows that 25% of LRGs are BCGs of our clusters and 36% of LRGs are cluster member galaxies. In our cluster sample, 63% of BCGs of r_{petro}<19.5 satisfy the SDSS LRG selection criteria.

We investigate the influence of major mergers on the radio emission of elliptical galaxies. We use a complete sample of close pairs, which contains 475 merging and 1828 non-merging paired elliptical galaxies of M_r<-21.5 selected from the Sloan Digital Sky Survey. In addition, a control sample of 2000 isolated field galaxies is used for comparison. We cross-identify the optical galaxies with the radio surveys of FIRST and NVSS. We find that the radio fraction of merging paired galaxies is about 6%, which is slightly higher than the 5% obtained for non-merging paired galaxies, although these values are consistent with each other owing to the large uncertainty caused by the limited sample. The radio fraction is twice as that of isolated galaxies, which is less than 3%. Radio emission of elliptical galaxies is only slightly affected by major mergers, but predominantly depends on their optical luminosities. Therefore, merging is not important in triggering the radio emission of elliptical galaxies.

To investigate the missing compact star of Supernova 1987A, we analyzed both the cooling and the heating processes of a possible compact star based on the upper limit of observational X-ray luminosity. From the cooling process we found that a solid quark-cluster star, which has a stiffer equation of state than that of conventional liquid quark star, has a heat capacity much smaller than a neutron star. It can cool down quickly, which can naturally explain the non-detection of a point source (neutron star or quark star) in X-ray band. On the other hand, we consider the heating process from magnetospheric activity and possible accretion, and obtain some constraints to the parameters of a possible pulsar. We conclude that a solid quark-cluster star can be fine with the observational limit in a large and acceptable parameter space. A pulsar with a short period and a strong magnetic field (or with a long period and a weak field) would has luminosity higher than the luminosity limit if the optical depth is not large enough to hide the compact star. The constraints of the pulsar parameters can be tested if the central compact object in 1987A is discovered by advanced facilities in the future.

We consider the curvature emission properties from relativistic particles streaming along magnetic field lines and co-rotating with pulsar magnetosphere. The co-rotation affects the trajectories of the particles and hence the emission properties, especially the polarization. We consider the modification of the particle velocity and acceleration due to the co-rotation. Curvature radiation from a single particle is calculated using the approximation of a circular path to the particle trajectory. Curvature radiation from particles at a given height actually contains the contributions from particles streaming along all the nearby field lines around the tangential point, forming the emission cone of 1/{\gamma}. The polarization patterns from the emission cone are distorted by the additional rotation, more serious for emission from a larger height. Net circular polarization can be generated by the density gradient in the emission cone. For three typical density models in the form of core, cone and patches, we calculate the polarization profiles for emission generated at a given height. We find that the circular polarization could have a single sign or sign reversal, depending on the density gradient along the rotation phase. The polarization profiles of the total curvature radiation from the whole open field line region, calculated by adding the emission from all possible heights, are similar to that from a dominating emission height. The circular polarization of curvature radiation has sign reversals in the patchy emission, while it has a single sign for the core emission, and is negligible for the cone emission.

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

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

We have finished the 6cm polarization survey of the Galactic plane using the Urumqi 25m radio telescope. It covers 10deg<l<230deg in Galactic longitude and |b| <5deg in Galactic latitude. The new polarization maps not only reveal new properties of the diffuse magnetized interstellar medium, but also are very useful for studying individual objects such as Hii regions, which may act as Faraday screens with strong regular magnetic fields inside, and supernova remnants for their polarization properties and spectra. The high sensitivity of the survey enables us to discover two new SNRs G178.2-4.2 and G25.3-2.1 and a number of Hii regions.

We calculate the correlation function of 13,904 galaxy clusters of z \leq 0.4 selected from the cluster catalog of Wen, Han & Liu. The correlation function can be fitted with a power-law model \xi(r)=(r/R_0)^{\gamma} on the scales of 10h^(-1)Mpc \leq r \leq 50h^(-1)Mpc, with a larger correlation length of R_0=18.84\pm0.27 h^(-1)Mpc for clusters with a richness of R \geq 15 and a smaller length of R_0=16.15 \pm 0.13 h^(-1)Mpc for clusters with a richness of R \geq 5. The power law index of \gamma=2.1 is found to be almost the same for all cluster subsamples. A pronounced baryon acoustic oscillations (BAO) peak is detected at r ~ 110 h^(-1)Mpc with a significance of ~ 1.9\sigma. By analyzing the correlation function in the range of 20 h^(-1)Mpc \leq r \leq 200 h^(-1)Mpc, we find the constraints on distance parameters are D_v(0.276)=1077\pm55(1\sigma) Mpc and h=0.73 \pm 0.039 (1\sigma), which are consistent with the WMAP 7-year cosmology. However, the BAO signal from the cluster sample is stronger than expected and leads to a rather low matter density \Omega_m h^2=0.093\pm0.0077 (1\sigma), which deviates from the WMAP 7-year result by more than 3 \sigma. The correlation function of the GMBCG cluster sample is also calculated and our detection of the BAO feature is confirmed.

Identification of high redshift clusters is important for studies of cosmology and cluster evolution. Using photometric redshifts of galaxies, we identify 631 clusters from the Canada-France-Hawaii Telescope (CFHT) Wide field, 202 clusters from the CHFT Deep field, 187 clusters from the Cosmic Evolution Survey (COSMOS) and 737 clusters from the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) field. The redshifts of these clusters are in the range of 0.1<z<1.6. Merging these cluster samples gives 1644 clusters in the four survey fields, of which 1088 are newly identified and more than half are from the large SWIRE field. Among 228 clusters of z>1, 191 clusters are newly identified, and most of them from the SWIRE field. With this large sample of high redshift clusters, we study the color evolution of the brightest cluster galaxies (BCGs). The colors r'-z' and r^+-m_{3.6\mu m} of the BCGs are consistent with a stellar population synthesis model in which the BCGs are formed at redshift z_f>2 and evolved passively. The colors g'-z' and B-m_{3.6\mu m} of the BCGs at redshifts z>0.8 are systematically bluer than the passive evolution model for galaxy formed at z_f~2, indicating star formation in high redshift BCGs.

We study the spectral and polarization properties of supernova remnants (SNRs) based on our 6cm survey data. The observations were taken from the Sino-German 6cm polarization survey of the Galactic plane. By using the integrated flux densities at 6cm together with measurements at other wavelengths from the literature we derive the global spectra of 50 SNRs. In addition, we use the observations at 6cm to present the polarization images of 24 SNRs. We derived integrated flux densities at 6cm for 51 small SNRs with angular sizes less than 1 degree. Global radio spectral indices were obtained in all the cases except for Cas A. For SNRs G15.1-1.6, G16.2-2.7, G16.4-0.5, G17.4-2.3, G17.8-2.6, G20.4+0.1, G36.6+2.6, G43.9+1.6, G53.6-2.2, G55.7+3.4, G59.8+1.2, G68.6-1.2, and G113.0+0.2, the spectra have been significantly improved. From our analysis we argue that the object G16.8-1.1 is probably an HII region instead of a SNR. Cas A shows a secular decrease in total intensity, and we measured a flux density of 688+/-35 Jy at 6cm between 2004 and 2008. Polarized emission from 25 SNRs were detected. For G16.2-2.7, G69.7+1.0, G84.2-0.8 and G85.9-0.6, the polarized emission is detected for the first time confirming them as SNRs. High frequency observations of SNRs are rare but important to establish their spectra and trace them in polarization in particular towards the inner Galaxy where Faraday effects are important.

We identify new strong lensing clusters of galaxies from the Sloan Digital Sky Survey III (SDSS DR8) by visually inspecting color images of a large sample of clusters of galaxies. We find 68 new clusters showing giant arcs in addition to 30 known lensing systems. Among 68 cases, 13 clusters are "almost certain" lensing systems with tangential giant arcs, 22 clusters are "probable" and 31 clusters are "possible" lensing systems. We also find two exotic systems with blue rings. The giant arcs have angular separations of 2.0"-25.7" from the bright central galaxies. We note that the rich clusters are more likely to be lensing systems, and the separations between arcs and the central galaxies increase with cluster richness.

(Abridged) We conducted new 6cm and 11cm observations of CTA 1 using the Urumqi 25-m and Effelsberg 100-m telescopes. Data at other wavelengths were included to investigate the spectrum and polarisation properties. We obtained new total intensity and polarisation maps at 6cm and 11cm with angular resolutions of 9.5 arcmin and 4.4 arcmin, respectively. We derived a spectral index of alpha=-0.63+/-0.05 based on the integrated flux densities at 408 MHz, 1420 MHz, 2639 MHz, and 4800 MHz. The spectral index map calculated from data at the four frequencies shows a clear steepening of the spectrum from the strong shell emission towards the north-western breakout region with weak diffuse emission. The decrease of the spectral index is up to about 0.3. The RM map derived from polarisation data at 6cm and 11cm shows a sharp transition between positive RMs in the north-eastern and negative RMs in the south-western part of the SNR. We note a corresponding RM pattern of extragalactic sources and propose the existence of a large-diameter Faraday screen in front of CTA 1, which covers the north-eastern part of the SNR. The RM of the Faraday screen is estimated to be about +45 rad/m2. A RM structure function of CTA 1 indicates a very regular magnetic field within the Faraday screen, which is larger than about 2.7 microG in case of 500 pc distance.

Supernova remnants (SNRs) were often discovered in radio surveys of the Galactic plane. Because of the surface-brightness limit of previous surveys, more faint or confused SNRs await discovery. The Sino-German $\lambda$6\ cm Galactic plane survey is a sensitive survey with the potential to detect new low surface-brightness SNRs. We want to identify new SNRs from the $\lambda$6\ cm survey map of the Galactic plane. We searched for new shell-like objects in the $\lambda$6\ cm survey maps, and studied their radio emission, polarization, and spectra using the $\lambda$6\ cm maps together with the $\lambda$11\ cm and $\lambda$21\ cm Effelsberg observations. Extended polarized objects with non-thermal spectra were identified as SNRs. We have discovered two new, large, faint SNRs, G178.2-4.2 and G25.1-2.3, both of which show shell structure. G178.2-4.2 has a size of 72 arcmin x 62 arcmin with strongly polarized emission being detected along its northern shell. The spectrum of G178.2-4.2 is non-thermal, with an integrated spectral index of $\alpha = -0.48\pm0.13$. Its surface brightness is $\Sigma_{1 GHz} = 7.2 x 10^{-23}{Wm^{-2} Hz^{-1} sr^{-1}}$, which makes G178.2-4.2 the second faintest known Galactic SNR. G25.1-2.3 is revealed by its strong southern shell which has a size of 80 arcmin x 30\arcmin. It has a non-thermal radio spectrum with a spectral index of $\alpha = -0.49\pm0.13$. Two new large shell-type SNRs have been detected at $\lambda$6\ cm in an area of 2200 deg^2 along the the Galactic plane. This demonstrates that more large and faint SNRs exist, but are very difficult to detect.

We develop a new method which measures the projected density distribution w_p(r_p)n of photometric galaxies surrounding a set of spectroscopically-identified galaxies, and simultaneously the projected correlation function w_p(r_p) between the two populations. In this method we are able to divide the photometric galaxies into subsamples in luminosity intervals when redshift information is unavailable, enabling us to measure w_p(r_p)n and w_p(r_p) as a function of not only the luminosity of the spectroscopic galaxy, but also that of the photometric galaxy. Extensive tests show that our method can measure w_p(r_p) in a statistically unbiased way. The accuracy of the measurement depends on the validity of the assumption in the method that the foreground/background galaxies are randomly distributed and thus uncorrelated with those galaxies of interest. Therefore, our method can be applied to the cases where foreground/background galaxies are distributed in large volumes, which is usually valid in real observations. We applied our method to data from SDSS including a sample of 10^5 LRGs at z~0.4 and a sample of about half a million galaxies at z~0.1, both of which are cross-correlated with a deep photometric sample drawn from the SDSS. On large scales, the relative bias factor of galaxies measured from w_p(r_p) at z~0.4 depends on luminosity in a manner similar to what is found at z~0.1, which are usually probed by autocorrelations of spectroscopic samples. On scales smaller than a few Mpc and at both z~0.4 and z~0.1, the photometric galaxies of different luminosities exhibit similar density profiles around spectroscopic galaxies at fixed luminosity and redshift. This provides clear support for the assumption commonly-adopted in HOD models that satellite galaxies of different luminosities are distributed in a similar way, following the dark matter distribution within their host halos.

The magnetosphere of a pulsar is composed of relativistic plasmas streaming along the magnetic field lines and corotating with the pulsar. We study the intrinsic Faraday rotation in the pulsar magnetosphere by critically examining the wave modes and the variations of polarization properties for the circularly polarized natural modes under various assumptions about the magnetosphere plasma properties. Since it is difficult to describe analytically the Faraday rotation effect in such a plasma, we use numerical integrations to study the wave propagation effects in the corotating magnetosphere. Faraday rotation effect is identified among other propagation effects, such as wave mode coupling and the cyclotron absorption. In a highly symmetrical electron-positron pair plasma, the Faraday rotation effect is found to be negligible. Only for asymmetrical plasmas, such as the electron-ion streaming plasma, can the Faraday rotation effect become significant, and the Faraday rotation angle is found to be approximately proportional to $\lambda^{0.5}$ instead of the usual $\lambda^2$-law. For such electrons-ion plasma of pulsar magnetosphere, the induced rotation measure becomes larger at higher frequencies, and should have opposite signs for the emissions from opposite magnetic poles.

Observations of large supernova remnants (SNRs) at high frequencies are rare, but provide valuable information about their physical properties. The total intensity and polarization properties of 16 large SNRs in the Galactic plane were investigated based on observations of the Urumqi $\lambda$6\ cm polarization survey of the Galactic plane with an angular resolution of 9$\farcm$5. We extracted total intensity and linear polarization maps of large SNRs from the Urumqi $\lambda$6\ cm survey, obtained their integrated flux densities, and derived the radio spectra in context with previously published flux densities at various frequencies. In particular, Effelsberg $\lambda$11\ cm and $\lambda$21\ cm survey data were used for calculating integrated flux densities. The $\lambda$6\ cm polarization data also delineate the magnetic field structures of the SNRs. We present the first total intensity maps at $\lambda$6\ cm for SNRs G106.3+2.7, G114.3+0.3, G116.5+1.1, G166.0+4.3 (VRO 42.05.01), G205.5+0.5 (Monoceros Nebula) and G206.9+2.3 (PKS 0646+06) and the first polarization measurements at $\lambda$6\ cm for SNRs G82.2+5.3 (W63), G106.3+2.7, G114.3+0.3, G116.5+1.1, G166.0+4.3 (VRO 42.05.01), G205.5+0.5 (Monoceros Nebula) and G206.9+2.3 (PKS 0646+06). Most of the newly derived integrated radio spectra are consistent with previous results. The new flux densities obtained from the Urumqi $\lambda$6\ cm, Effelsberg $\lambda$11\ cm and $\lambda$21\ cm surveys are crucial to determine the spectra of SNR G65.1+0.6, G69.0+2.7 (CTB 80), G93.7-0.2 and G114.3+0.3. We find that G192.8$-$1.1 (PKS 0607+17) consists of background sources, \ion{H}{II} regions and the extended diffuse emission of thermal nature, and conclude that G192.8$-$1.1 is not a SNR.

Ultraluminous infrared galaxies (ULIRGs) have several types according to dominance of starburst or AGN component. We made stellar population analysis for a sample of 160 ULIRGs to study the evolution of ULIRGs. We found that the dominance of intermediate-age and old stellar populations increases along the sequence of HII-like ULIRGs, Seyfert-HII composite ULIRGs, and Seyfert 2 ULIRGs. Consequently the typical mean stellar age and the stellar mass increase along the sequence. Comparing the gas mass estimated from the CO measurements with the stellar mass estimated from the optical spectra, we found that gas fraction is anti-correlated with the stellar mass. HII-like ULIRGs with small stellar masses do not possess enough gas and the total mass, and therefore have no evolution connections with massive Seyfert 2 ULIRGs. Only massive ULIRGs may follow the evolution sequence toward AGNs, and massive HII-like ULIRGs are probably in an earlier stage of the sequence.

Linear polarization of diffuse Galactic emission is a signature of magnetic fields in the interstellar medium of our Galaxy. Observations at high frequencies are less affected by Faraday depolarization than those at lower frequencies and are able to detect polarized emission from more distant Galactic regions. We attempt to perform a sensitive survey of the polarized emission from the Galactic disk at 6 cm wavelength. We made polarization observations of the Galactic plane using the Urumqi 25-m telescope at 6 cm covering the area of 60deg< l <129deg and |b|<5deg. Missing large-scale structures in polarization were restored by extrapolation of the WMAP polarization data. We present the 6cm total intensity and linear polarization maps of the surveyed region. We identify two new extended HII regions G98.3-1.6 and G119.6+0.4 in this region. Numerous polarized patches and depolarization structures are visible in the polarization maps. Depolarization along the periphery of a few HII complexes was detected and can be explained by a Faraday screen model. We discuss some prominent depolarization HII regions, which have regular magnetic fields of several $\mu$G. Structure functions of U, Q and PI images of the entire 6cm survey region of 10deg< l < 230deg exhibit much larger fluctuation power towards the inner Galaxy, which suggests a higher turbulence in the arm regions of the inner Galaxy. The Sino-German 6cm survey reveals new properties of the diffuse magnetized interstellar medium. The survey is also very useful for studying individual objects such as HII regions, which act as Faraday screens and have high rotation measures and therefore strong regular magnetic fields inside the regions.

Cosmic-ray proton and helium spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass experiment flown for 42 days in Antarctica in the 2004-2005 austral summer season. High-energy cosmic-ray data were collected at an average altitude of ~38.5 km with an average atmospheric overburden of ~3.9 g cm$^{-2}$. Individual elements are clearly separated with a charge resolution of ~0.15 e (in charge units) and ~0.2 e for protons and helium nuclei, respectively. The measured spectra at the top of the atmosphere are represented by power laws with a spectral index of -2.66 $\pm$ 0.02 for protons from 2.5 TeV to 250 TeV and -2.58 $\pm$ 0.02 for helium nuclei from 630 GeV/nucleon to 63 TeV/nucleon. They are harder than previous measurements at a few tens of GeV/nucleon. The helium flux is higher than that expected from the extrapolation of the power law fitted to the lower-energy data. The relative abundance of protons to helium nuclei is 9.1 $\pm$ 0.5 for the range from 2.5 TeV/nucleon to 63 TeV/nucleon. This ratio is considerably smaller than the previous measurements at a few tens of GeV/nucleon.

The Vela pulsar represents a distinct group of {\gamma}-ray pulsars. Fermi {\gamma}-ray observations reveal that it has two sharp peaks (P1 and P2) in the light curve with a phase separation of 0.42 and a third peak (P3) in the bridge. The location and intensity of P3 are energy-dependent. We use the 3D magnetospheric model for the annular gap and core gap to simulate the {\gamma}-ray light curves, phase-averaged and phase-resolved spectra. We found that the acceleration electric field along a field line in the annular gap region decreases with heights. The emission at high energy GeV band is originated from the curvature radiation of accelerated primary particles, while the synchrotron radiation from secondary particles have some contributions to low energy {\gamma}-ray band (0.1 - 0.3 GeV). The {\gamma}-ray light curve peaks P1 and P2 are generated in the annular gap region near the altitude of null charge surface, whereas P3 and the bridge emission is generated in the core gap region. The intensity and location of P3 at different energy bands depend on the emission altitudes. The radio emission from the Vela pulsar should be generated in a high-altitude narrow regions of the annular gap, which leads to a radio phase lag of ~ 0.13 prior to the first {\gamma}-ray peak.

Brightest Cluster Galaxies (BCGs) are mostly elliptical galaxies and very rarely have prominent star formation. We found that five out of 8,812 BCGs are E+A (i.e. post-starburst) galaxies, having the H$\delta$~absorption line with an equivalent width $>2.5\AA$ and no distinct emission lines in [O II] and H$\alpha$. The E+A features we identified from the BCGs for the first time are not as significant as those in general galaxies, indicating that historically the star formation were not very violent.

Estimating black hole masses of blazars is still a big challenge. Because of the contamination of jets, using the previously suggested size -- continuum luminosity relation can overestimate the broad line region (BLR) size and black hole mass for radio-loud AGNs, including blazars. We propose a new relation between the BLR size and $H_{\beta}$ emission line luminosity and present evidences for using it to get more accurate black hole masses of radio-loud AGNs. For extremely radio-loud AGNs such as blazars with weak/absent emission lines, we suggest to use the fundamental plane relation of their elliptical host galaxies to estimate the central velocity dispersions and black hole masses, if their velocity dispersions are not known but the host galaxies can be mapped. The black hole masses of some well-known blazars, such as OJ 287, AO 0235+164 and 3C 66B, are obtained using these two methods and the M - $\sigma$ relation. The implications of their black hole masses on other related studies are also discussed.

(Abridged) We have conducted a total intensity and polarization survey of the Galactic plane at 6 cm using the Urumqi 25 m telescope for the Galactic longitude range of 10 deg<l<60 deg and the Galactic latitude range of |b|<5 deg. Missing absolute zero levels of Stokes U and Q maps were restored by extrapolating the WMAP five-year K-band polarization data. For total intensities we recovered missing large-scale components by referring to the Effelsberg 11 cm survey. Total intensity and polarization maps are presented with an angular resolution of 9.5 arcmin and a sensitivity of 1 mK and 0.5 mK in total and polarized intensity, respectively. The 6 cm polarized emission in the Galactic plane originates within about 4 kpc distance, which increases for polarized emission out of the plane. The polarization map shows "patches", "canals" and "voids" with no correspondence in total intensity. We attribute the patches to turbulent magnetic field cells. Canals are caused by abrupt variation of polarization angles at the boundaries of patches rather than by foreground Faraday Screens. The superposition of foreground and Faraday Screen rotated background emission almost cancels polarized emission locally, so that polarization voids appear. By modelling the voids, we estimate the Faraday Screen's regular magnetic field along the line-of-sight to be larger than about 8 microG. We separated thermal (free-free) and non-thermal (synchrotron) emission according to their different spectral indices. The spectral index for the synchrotron emission was based on WMAP polarization data. The fraction of thermal emission at 6 cm is about 60% in the plane.

A sample of 129 unresolved radio sources with ultrahigh linear polarization (>30 per cent) has been selected from the NRAO VLA Sky Survey. Such high average linear polarization is unusual in extragalactic sources. Higher resolution Australia Telescope Compact Array and Very Large Array observations confirm the high average polarization but find that most of these sources are extended. The Sloan Digital Sky Survey spectroscopy, where available, shows that the optical counterparts are elliptical galaxies with no detectable emission lines. The optical spectra, radio luminosity, linear size and spectral index of these sources are typical of radio-loud active galactic nuclei. Galaxy counts within a 1 Mpc radius of the radio sources show that these highly polarized sources are in environments similar to their low polarization (<2 per cent) counterparts. Similarly, the line-of-sight environments of the ultrahigh polarization sources are on average indistinguishable from those of the low-polarization sources. We conclude that the extraordinarily high average polarization must be due to intrinsic properties of the sources, such as an extremely ordered source magnetic field, low internal thermal plasma density or a preferential orientation of the source magnetic field perpendicular to the line of sight.

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $\gamma-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $\gamma$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $\gamma$-B optical depth for $\gamma$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.

We study possible correlations between ultrahigh energy cosmic rays (UHECRs), observed by Auger, AGASA and Yakutsk, and nearby active galactic nuclei (AGNs) and $Fermi$ sources. We consider the deflection effects by a Galactic magnetic field (GMF) model constrained by the most updated measurements. We found that the average deflection angles of UHECRs by the Galactic magnetic fields are less than $4^\circ$. A correlation between the Auger cosmic-ray events and nearby AGNs with a significance level of $\sim 4\sigma$ was found for the Auger UHECR data sets with or without deflection correction. No correlation was found between the AGASA/Yakutsk events with nearby AGNs. Marginal correlations between the Auger events and the $Fermi$ sources, and between AGASA events and $Fermi$ AGNs were found when the deflections calculated by the GMF model were considered. However, no correlation was found between the Yakutsk data and $Fermi$ sources. Some $Fermi$ sources are close to the arrival directions of UHECR events detected by Auger, AGASA, and Yakutsk, most of which are probably chance coincidence rather than objects producing UHECRs in the nearby Universe. Four $Fermi$ sources, NGC 4945, ESO 323-G77, NGC 6951, and Cen A, within 100~Mpc have UHECR events within $3.1^{\circ}$ from their positions, which could potentially be cosmic ray accelerators. However, the association can only be confirmed if more UHECRs are preferably detected in these directions.

We present high resolution images of the bipolar outflow from W51e2, which are produced from the Submillimeter Array archival data observed for CO(3-2) and HCN(4-3) lines with angular resolutions of 0.8" x 0.6" and 0.3" x 0.2", respectively. The images show that the powerful outflow originates from the protostellar core W51e2-E rather than from the ultracompact HII region W51e2-W. The kinematic timescale of the outflow from W51e2-E is about 1000 yr, younger than the age (~5000 yr) of the ultracompact HII region W51e2-W. A large mass loss rate of ~1 x 10^{-3} M_sun yr^{-1} and a high mechanical power of 120 L_sun are inferred, suggesting that an O star or a cluster of B stars are forming in W51e2-E. The observed outflow activity along with the inferred large accretion rate indicates that at present W51e2-E is in a rapid phase of star formation.

The mass function of galaxy clusters is a powerful tool to constrain cosmological parameters, e.g., the mass fluctuation on the scale of 8 h^{-1} Mpc, sigma_8, and the abundance of total matter, Omega_m. We first determine the scaling relations between cluster mass and cluster richness, summed r-band luminosity and the global galaxy number within a cluster radius. These relations are then used to two complete volume-limited rich cluster samples which we obtained from the Sloan Digital Sky Survey (SDSS). We estimate the masses of these clusters and determine the cluster mass function. Fitting the data with a theoretical expression, we get the cosmological parameter constraints in the form of sigma_8(Omega_m/0.3)^{alpha}=beta and find out the parameters of alpha=0.40-0.50 and beta=0.8-0.9, so that sigma_8=0.8-0.9 if Omega_m=0.3. Our sigma_8 value is slightly higher than recent estimates from the mass function of X-ray clusters and the Wilkinson Microwave Anisotropy Probe (WMAP) data, but consistent with the weak lensing statistics.

Linearly polarized Galactic synchrotron emission provides valuable information about the properties of the Galactic magnetic field and the interstellar magneto-ionic medium, when Faraday rotation along the line of sight is properly taken into account. We aim to survey the Galactic plane at 6 cm including linear polarization. At such a short wavelength Faraday rotation effects are in general small and the Galactic magnetic field properties can be probed to larger distances than at long wavelengths. The Urumqi 25-m telescope is used for a sensitive 6 cm survey in total and polarized intensities. WMAP K-band (22.8 GHz) polarization data are used to restore the absolute zero-level of the Urumqi U and Q maps by extrapolation. Total intensity and polarization maps are presented for a Galactic plane region of 129 degree < l < 230 degree and |b| < 5 degree in the anti-centre with an angular resolution of 9'5 and an average sensitivity of 0.6 mK and 0.4 mK Tb in total and polarized intensity, respectively. We briefly discuss the properties of some extended Faraday Screens detected in the 6 cm polarization maps. The Sino-German 6 cm polarization survey provides new information about the properties of the magnetic ISM. The survey also adds valuable information for discrete Galactic objects and is in particular suited to detect extended Faraday Screens with large rotation measures hosting strong regular magnetic fields.

Pulsed high energy radiation from pulsars is not yet completely understood. In this paper, we use the 3D self-consistent annular gap model to study light curves for both young and millisecond pulsars observed by the Fermi Gamma-ray Space Telescope. The annular gap can generate high energy emission for short-period pulsars. The annular gap regions are so large that they have enough electric potential drop to accelerate charged particles to produce gamma-ray photons. For young pulsars, the emission region is from the neutron star surface to about half of the light cylinder radius, and the peak emissivity is in the vicinity of the null charge surface. The emission region for the millisecond pulsars is located much lower than that of the young pulsars. The higher energy gamma-ray emission comes from higher altitudes in the magnetosphere. We present the simulated light curves for three young pulsars (the Crab, the Vela, the Geminga) and three millisecond pulsars (PSR J0030+0451, PSR J0218+4232, PSR J0437-3715) using the annular gap model. Our simulations can reproduce the main properties of observed light curves.

The balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment launched five times from Antarctica has achieved a cumulative flight duration of about 156 days above 99.5% of the atmosphere. The instrument is configured with complementary and redundant particle detectors designed to extend direct measurements of cosmic-ray composition to the highest energies practical with balloon flights. All elements from protons to iron nuclei are separated with excellent charge resolution. Here we report results from the first two flights of ~70 days, which indicate hardening of the elemental spectra above ~200 GeV/nucleon and a spectral difference between the two most abundant species, protons and helium nuclei. These results challenge the view that cosmic-ray spectra are simple power laws below the so-called knee at ~1015 eV. This discrepant hardening may result from a relatively nearby source, or it could represent spectral concavity caused by interactions of cosmic rays with the accelerating shock. Other possible explanations should also be investigated.

The Cosmic Ray Energetics And Mass (CREAM) is a balloon-borne experiment designed to measure the composition and energy spectra of cosmic rays of charge Z = 1 to 26 up to an energy of ~ 10^15 eV. CREAM had two successful flights on long-duration balloons (LDB) launched from Mc- Murdo Station, Antarctica, in December 2004 and December 2005. CREAM-I achieves a substantial measurement redundancy by employing multiple detector systems, namely a Timing Charge Detector and a Silicon Charge Detector (SCD) for particle identification, and a Transition Radiation Detector and a sampling tungsten/scintillating-fiber ionization calorimeter (CAL) for energy measurement. In this paper, preliminary energy spectra of various elements measured with CAL/SCD during the first 42-day flight are presented.

During its second Antarctic flight, the CREAM (Cosmic Ray Energetics And Mass) balloon experiment collected data for 28 days, measuring the charge and the energy of cosmic rays (CR) with a redundant system of particle identification and an imaging thin ionization calorimeter. Preliminary direct measurements of the absolute intensities of individual CR nuclei are reported in the elemental range from carbon to iron at very high energy.

We present new measurements of the energy spectra of cosmic-ray (CR) nuclei from the second flight of the balloon-borne experiment CREAM (Cosmic Ray Energetics And Mass). The instrument (CREAM-II) was comprised of detectors based on different techniques (Cherenkov light, specific ionization in scintillators and silicon sensors) to provide a redundant charge identification and a thin ionization calorimeter capable of measuring the energy of cosmic rays up to several hundreds of TeV. The data analysis is described and the individual energy spectra of C, O, Ne, Mg, Si and Fe are reported up to ~ 10^14 eV. The spectral shape looks nearly the same for all the primary elements and can be expressed as a power law in energy E^{-2.66+/-0.04}. The nitrogen absolute intensity in the energy range 100-800 GeV/n is also measured.

The Cosmic Ray Energetics And Mass (CREAM) balloon experiment had two successful flights in 2004/05 and 2005/06. It was designed to perform energy measurements from a few GeV up to 1000 TeV, taking advantage of different detection techniques. The first instrument, CREAM-1, combined a transition radiation detector with a calorimeter to provide independent energy measurements of cosmicraynuclei. Each detector was calibrated with particle beams in a limited range of energies. In order to assess the absolute energy scale of the instrument and to investigate the systematic effects of each technique, a cross-calibration was performed by comparing the two independent energy estimates on selected samples of oxygen and carbon nuclei.

The Cosmic Ray Energetics And Mass (CREAM) calorimeter is designed to measure the spectra of cosmic-ray particles over the energy range from ~10^11 eV to ~10^15 eV. Its first flight as part of the CREAM-I balloon-borne payload in Antarctica during the 2004/05 season resulted in a recordbreaking 42 days of exposure. Calorimeter calibration using various beam test data will be discussed in an attempt to assess the uncertainties of the energy measurements.

The CREAM calorimeter, designed to measure the spectra of cosmic-ray nuclei from under 1 TeV to 1000 TeV, is a 20 radiation length (X0) deep sampling calorimeter. The calorimeter is comprised of 20 layers of tungsten interleaved with 20 layers of scintillating fiber ribbons, and is preceded by a pair of graphite interaction targets providing about 0.42 proton interaction lengths (\lambda int). The calorimeter was placed in one of CERN's SPS accelerator beams for calibration and testing. Beams of 150 GeV electrons were used for calibration, and a variety of electron, proton, and nuclear fragment beams were used to test the simulation model of the detector. In this paper we discuss the performance of the calorimeter in the electron beam and compare electron beam data with simulation results.

We present high-resolution continuum images of the W51e2 complex processed from archival data of the Submillimeter Array (SMA) at 0.85 and 1.3 mm and the Very Large Array (VLA) at 7 and 13 mm. We also made line images and profiles of W51e2 for three hydrogen radio recombination lines (H26\alpha, H53\alpha, and H66\alpha) and absorption of two molecular lines of HCN(4-3) and CO(2-1). At least four distinct continuum components have been detected in the 3" region of W51e2 from the SMA continuum images at 0.85 and 1.3 mm with resolutions of 0.3"x0.2" and 1.4"x0.7", respectively. The west component, W51e2-W, coincides with the UC HII region reported from previous radio observations. The H26\alpha line observation reveals an unresolved hyper-compact ionized core (<0.06" or <310 AU) with a high electron temperature of 1.2x10^4 K, with corresponding emission measure EM>7x10^{10} pc cm^{-6} and electron density N_e>7x10^6 cm^{-3}. The inferred Lyman continuum flux implies that the HII region W51e2-W requires a newly formed massive star, an O8 star or a cluster of B-type stars, to maintain the ionization. The east component, W51e2-E, has a total mass of ~140 M_{\sun} according to our SED analysis and a large infall rate of > 1.3x10^{-3} M_{\sun}yr^{-1} inferred from the absorption of HCN. W51e2-E appears to be the accretion center in W51e2 and to host one or more growing massive proto-stars. Located 2" northwest from W51e2-E, W51e2-NW is not detected in the continuum emission at \lambda>=7 mm. Along with the maser activities previously observed, our analysis suggests that W51e2-NW is at an earlier phase of star formation. W51e2-N is located 2" north of W51e2-E and has only been detected at 1.3 mm with a lower angular resolution (~1"), suggesting that it is a primordial, massive gas clump in the W51e2 complex.

We present new measurements of the energy spectra of cosmic-ray (CR) nuclei from the second flight of the balloon-borne experiment Cosmic Ray Energetics And Mass (CREAM). The instrument included different particle detectors to provide redundant charge identification and measure the energy of CRs up to several hundred TeV. The measured individual energy spectra of C, O, Ne, Mg, Si, and Fe are presented up to $\sim 10^{14}$ eV. The spectral shape looks nearly the same for these primary elements and it can be fitted to an $E^{-2.66 \pm 0.04}$ power law in energy. Moreover, a new measurement of the absolute intensity of nitrogen in the 100-800 GeV/$n$ energy range with smaller errors than previous observations, clearly indicates a hardening of the spectrum at high energy. The relative abundance of N/O at the top of the atmosphere is measured to be $0.080 \pm 0.025 $(stat.)$ \pm 0.025 $(sys.) at $\sim $800 GeV/$n$, in good agreement with a recent result from the first CREAM flight.

Clusters of galaxies in most previous catalogs have redshifts z<0.3. Using the photometric redshifts of galaxies from the Sloan Digital Sky Survey Data Release 6 (SDSS DR6), we identify 39,668 clusters in the redshift range 0.05< z <0.6 with more than eight luminous (M_r<-21) member galaxies. Cluster redshifts are estimated accurately with an uncertainty less than 0.022. The contamination rate of member galaxies is found to be roughly 20%, and the completeness of member galaxy detection reaches to ~90%. Monte Carlo simulations show that the cluster detection rate is more than 90% for massive (M_{200}>2\times10^{14} M_{\odot}) clusters of z<0.42. The false detection rate is ~5%. We obtain the richness, the summed luminosity, and the gross galaxy number within the determined radius for identified clusters. They are tightly related to the X-ray luminosity and temperature of clusters. Cluster mass is related to the richness and summed luminosity with M_{200}\propto R^{1.90\pm0.04} and M_{200}\propto L_r^{1.64\pm0.03}, respectively. In addition, 685 new candidates of X-ray clusters are found by cross-identification of our clusters with the source list of the ROSAT X-ray survey.

We study the propagation effects of radio waves in a pulsar magnetosphere, composed of relativistic electron-positron pair plasmas streaming along the magnetic field lines and corotating with the pulsar. We critically examine the various physical effects that can potentially influence the observed wave intensity and polarization, including resonant cyclotron absorption, wave mode coupling due to pulsar rotation, wave propagation through quasi-tangential regions (where the photon ray is nearly parallel to the magnetic field) and mode circularization due to the difference in the electron/positron density/velocity distributions. We numerically integrate the transfer equations for wave polarization in the rotating magnetosphere, taking account of all the propagation effects in a self-consistent manner. For typical magnetospheric plasma parameters produced by pair cascade, we find that the observed radio intensity and polarization profiles can be strongly modified by the propagation effects. For relatively large impact parameter (the minimum angle between the magnetic dipole axis and the line of sight), the polarization angle profile is similar to the prediction from the Rotating Vector Model, except for a phase shift and an appreciable circular polarization. For smaller impact parameter, the linear polarization position angle may exhibit a sudden $90^o$ jump due to the quasi-tangential propagation effect, accompanied by complex circular polarization profile. Some applications of our results are discussed, including the origin of non-gaussion pulse profiles, the relationship between the position angle profile and circular polarization in conal-double pulsars, and the orthogonal polarization modes.

We have made a careful selection of a large complete volume-limited sample (1209) of projected close pairs (7<r_p<50 kpc) of luminous early-type galaxies (M_r<-21.5) in the local universe (z<0.12) from the SDSS data. 249 (21%) of them show interaction features, which suggests that about 0.8% of the galaxies are merging. We derived a comoving volume merger rate of ~(1.0+/-0.4)\times 10^{-5} Mpc^{-3} Gyr^{-1} for luminous early-type galaxies. This is a direct observational determination of the merger rate of luminous galaxies in the local universe. We also obtained the chirp mass distribution of supermassive black hole (SMBH) binary following log[Phi(log M /M_{odot})]=(21.7+/-4.2)-(3.0+/-0.5)log M/M_{odot}. With less assumptions than previous works, we estimated the strain amplitude of the gravitational wave (GW) background from coalescence of SMBH binaries at frequency 10^{-9}-10^{-7} Hz.

Ultra-luminous infrared galaxies (ULIRGs) are interesting objects with dramatic properties. Many efforts have been made to understand the physics of their luminous infrared emission and evolutionary stages. However, a large ULIRG sample is still needed to study the properties of their central black holes (BHs), the BH-host galaxy relation, and their evolution. We identified 308 ULIRGs from the Sloan Digital Sky Survey Data Release 6, and classified them into the NL ULIRGs (with only narrow emission lines) and the Type I ULIRGs (with broad emission lines). About 56% of ULIRGs in our total sample show interaction features, and this percentage is 79% for redshift z < 0.2. Optical identifications of these ULIRGs show that the AGN percentage is at least 49%, and the percentage increases with the infrared luminosity. We found 62 Type I ULIRGs, and estimated their BH masses and velocity dispersions from their optical spectra. Together with known Type I ULIRGs in the literature, a sample of 90 Type I ULIRGs enables us to make a statistical study. We found that the BH masses of Type I ULIRGs are typically smaller than those of PG QSOs, and most Type I ULIRGs follow the M_{BH}-sigma relation. However, some ULIRGs with larger Eddington ratio deviate from this relation, even the line width of the [OIII] narrow line (NL) core or the [SII] line was used as the surrogate of velocity dispersion. This implies that at least some ULIRGs are probably still in the early evolution stage toward QSOs. The anti-correlation between the mass deviation from the M_{BH}-sigma relation and the Eddington ratio supports that the evolution of Type I ULIRGs is probably followed by the building up of the M_{BH}-sigma relation and the evolution to the QSO phase.

Laboratory spectroscopy of non-thermal equilibrium plasmas photoionized by intense radiation is a key to understanding compact objects, such as black holes, based on astronomical observations. This paper describes an experiment to study photoionizing plasmas in laboratory under well-defined and genuine conditions. Photoionized plasma is here generated using a 0.5-keV Planckian x-ray source created by means of a laser-driven implosion. The measured x-ray spectrum from the photoionized silicon plasma resembles those observed from the binary stars Cygnus X-3 and Vela X-1 with the Chandra x-ray satellite. This demonstrates that an extreme radiation field was produced in the laboratory, however, the theoretical interpretation of the laboratory spectrum significantly contradicts the generally accepted explanations in x-ray astronomy. This model experiment offers a novel test bed for validation and verification of computational codes used in x-ray astronomy.

Pulsars have been recognized as normal neutron stars, but sometimes argued as quark stars. {\it Sub-millisecond pulsars, if detected, would play an essential and important role in distinguishing quark stars from neutron stars.} We focus on the formation of such sub-millisecond pulsars in this paper. A new approach to form a sub-millisecond pulsar (quark star) via accretion induced collapse (AIC) of a white dwarf is investigated here. Under this AIC process, we found that: (1) almost all the newborn quark stars could have an initial spin period of $\sim 0.1$ ms; (2) the nascent quark stars (even with a low mass) have sufficiently high spin-down luminosity and satisfy the conditions for pair production and sparking process to be as sub-millisecond radio pulsars; (3) in most cases, the timescales of newborn quark stars in the phase of spin period $< 1$ (or $<0.5$) ms can be long enough to be detected. As a comparison, an accretion spin-up process (for both neutron and quark stars) is also investigated. It is found that, quark stars formed through AIC process can have shorter periods ($\leq$ 0.5 ms); while the periods of neutron stars formed in accretion spin-up process must be longer than 0.5ms. Thus if a pulsar with a period less than 0.5 ms can be identified in the future, it should be a quark star.

Many physical properties of this SNR such as spectrum and polarization can only be investigated by radio observations. The $\lambda$11 cm and $\lambda$6\ cm continuum and polarization observations of SNR G65.2+5.7 were made with the Effelsberg 100-m and the Urumqi 25-m telescopes, respectively, to investigate the integrated spectrum, the spectral index distribution, and the magnetic field properties. $\lambda$21 cm archival data from the Effelsberg 100-m telescope have been also used. The integrated flux densities of G65.2+5.7 at $\lambda 11$ cm and $\lambda 6$ cm are $21.9\pm3.1$ Jy and 16.8$\pm$1.8 Jy, respectively. The power-law spectrum ($S\sim\nu^{\alpha}$) is well fitted by $\alpha = -0.58\pm0.07$ from 83 MHz to 4.8 GHz. Spatial spectral variations are small. Along the northern shell strong depolarizion is observed at both wavelengths. The southern filamentary shell of SNR G65.2+5.7 is polarized up to 54% at $\lambda 6$ cm. There is significant depolarization at $\lambda 11$ cm and confusion with diffuse polarized Galactic emission. Using equipartition principle, we estimated the magnetic field strength for the southern filamentary shell about 20 $\mu$G (filling factor 1) to 50 $\mu$G (filling factor 0.1). A faint HI shell may be associated with the SNR. Despite its unusual strong X-ray and optical emission and its very low surface brightness, the radio properties of SNR G65.2+5.7 are found to be typical for evolved shell type SNRs. SNR G65.2+5.7 may be expanding in a preblown cavity as indicated by a deficit of HI gas and a possible HI-shell.

The spiral structure of our Milky Way Galaxy is not yet known. HII regions and giant molecular clouds are the most prominent spiral tracers. We collected the spiral tracer data of our Milky Way from the literature, namely, HII regions and giant molecular clouds (GMCs). With weighting factors based on the excitation parameters of HII regions or the masses of GMCs, we fitted the distribution of these tracers with models of two, three, four spiral-arms or polynomial spiral arms. The distances of tracers, if not available from stellar or direct measurements, were estimated kinetically from the standard rotation curve of Brand & Blitz (1993) with $R_0$=8.5 kpc, and $\Theta_0$=220 km s$^{-1}$ or the newly fitted rotation curves with $R_0$=8.0 kpc and $\Theta_0$=220 km s$^{-1}$ or $R_0$=8.4 kpc and $\Theta_0$=254 km s$^{-1}$. We found that the two-arm logarithmic model cannot fit the data in many regions. The three- and the four-arm logarithmic models are able to connect most tracers. However, at least two observed tangential directions cannot be matched by the three- or four-arm model. We composed a polynomial spiral arm model, which can not only fit the tracer distribution but also match observed tangential directions. Using new rotation curves with $R_0$=8.0 kpc and $\Theta_0$=220 km s$^{-1}$ and $R_0$=8.4 kpc and $\Theta_0$=254 km s$^{-1}$ for the estimation of kinematic distances, we found that the distribution of HII regions and GMCs can fit the models well, although the results do not change significantly compared to the parameters with the standard $R_0$ and $\Theta_0$.

We report on polarimetric observations of 100 pulsars centered on 774 MHz, made using the Green Bank Telescope, presenting their polarization profiles and polarized flux densities and comparing them with previous observations when possible. For 67 pulsars, these are the first such measurements made. Polarization profiles of 8 millisecond pulsars in our sample show wide profiles and flat position-angle curves. Strong linear polarization, sometimes approaching 100% of the total intensity, has been detected in all or a part of the average pulse profiles of some pulsars. In general, circular polarization is very weak, although it is observed to be extremely strong in the leading component of PSR J1920+2650. Sense reversal of circular polarization as a function of pulse phase has been detected from both core and other components of more than 20 pulsars. Any relationship between the spin-down luminosity and the percentage of linear polarization is not evident in our data at this frequency.

Supermassive black hole (SMBH) coalescence in galaxy mergers is believed to be one of the primary sources of very low frequency gravitational waves (GWs). Significant contribution of the GWs comes from mergers of massive galaxies with redshifts z<2. Very few previous studies gave the merger rate of massive galaxies. % We selected a large sample (1209) of close pairs of galaxies with projected separations 7<r_p<50 kpc from 87,889 luminous early-type galaxies (M_r<-21.5) from the Sloan Digital Sky Survey Data Release 6. These pairs constitute a complete volume-limited sample in the local universe (z<0.12). Using our newly developed technique, 249 mergers have been identified by searching for interaction features. From them, we found that the merger fraction of luminous early-type galaxies is 0.8%, and the merger rate in the local universe is % R_g=(1.0+/-0.4)*10^{-5} Mpc^{-3} Gyr^{-1}} % with an uncertainty mainly depending on the merging timescale. % We estimated the masses of SMBHs in the centers of merging galaxies based on their luminosities. We found that the chirp mass distribution of the SMBH binaries follows a power law with an index of -3.0+/-0.5 in the range 5*10^8--5*10^{9} M_{\odot}. % Using the SMBH population in the mergers and assuming that the SMBHs can be efficiently driven into the GW regime, we investigated the stochastic GW background in the frequency range 10^{-9}--10^{-7} Hz. We obtained the spectrum of the GW background of h_c(f)=10^{-15}(f/yr^{-1})^{-2/3}, which is one magnitude higher than that obtained by Jaffe & Backer in 2003, but consistent with those calculated from galaxy-formation models.

We report the discovery of 4 strong gravitational lensing systems by visual inspections of the Sloan Digital Sky Survey images of galaxy clusters in Data Release 6 (SDSS DR6). Two of the four systems show Einstein rings while the others show tangential giant arcs. These arcs or rings have large angular separations (>8") from the bright central galaxies and show bluer color compared with the red cluster galaxies. In addition, we found 5 probable and 4 possible lenses by galaxy clusters.

The magnetic structure in the Galactic disk, the Galactic center and the Galactic halo can be delineated more clearly than ever before. In the Galactic disk, the magnetic structure has been revealed by starlight polarization within 2 or 3 kpc of the Solar vicinity, by the distribution of the Zeeman splitting of OH masers in two or three nearby spiral arms, and by pulsar dispersion measures and rotation measures in nearly half of the disk. The polarized thermal dust emission of clouds at infrared, mm and submm wavelengths and the diffuse synchrotron emission are also related to the large-scale magnetic field in the disk. The rotation measures of extragalactic radio sources at low Galactic latitudes can be modeled by electron distributions and large-scale magnetic fields. The statistical properties of the magnetized interstellar medium at various scales have been studied using rotation measure data and polarization data. In the Galactic center, the non-thermal filaments indicate poloidal fields. There is no consensus on the field strength, maybe mG, maybe tens of uG. The polarized dust emission and much enhanced rotation measures of background radio sources are probably related to toroidal fields. In the Galactic halo, the antisymmetric RM sky reveals large-scale toroidal fields with reversed directions above and below the Galactic plane. Magnetic fields from all parts of our Galaxy are connected to form a global field structure. More observations are needed to explore the untouched regions and delineate how fields in different parts are connected.

Pulsars are fantastic objects, which show the extreme states of matters and plasma physics not understood yet. Pulsars can be used as probes for the detection of interstellar medium and even the gravitational waves. Here I review the basic facts of pulsars which should attract students to choose pulsar studies as their future projects.

The magnetic fields of our Milky Way galaxy are the main agent for cosmic rays to transport. In the last decade, much new knowledge has been gained from measurements of the Galactic magnetic fields. In the Galactic disk, from the RMs of a large number of newly discovered pulsars, the large-scale magnetic fields along the spiral arms have been delineated in a much larger region than ever before, with alternating directions in the arm and interarm regions. The toroidal fields in the Galactic halo were revealed to have opposite directions below and above the Galactic plane, which is an indication of an A0 mode dynamo operating in the halo. The strength of large-scale fields obtained from pulsar RM data has been found to increase exponentially towards the Galactic center. Compared to the steep Kolmogorov spectrum of magnetic energy at small scales, the large-scale magnetic fields show a shallow broken spatial magnetic energy spectrum.

The Global Magneto-Ionic Medium Survey (GMIMS) is a project to map the diffuse polarized emission over the entire sky, Northern and Southern hemispheres, from 300 MHz to 1.8 GHz. With an angular resolution of 30 - 60 arcmin and a frequency resolution of 1 MHz or better, GMIMS will provide the first spectro-polarimetric data set of the large-scale polarized emission over the entire sky, observed with single-dish telescopes. GMIMS will provide an invaluable resource for studies of the magneto-ionic medium of the Galaxy in the local disk, halo, and its transition.

Pulsars have been recognized as normal neutron stars or quark stars. Sub-millisecond pulsars, if detected, would play an essential and important role in distinguishing quark stars from neutron stars. A key question is how sub-millisecond pulsars could form. Both sub-Keplerian (for neutron and quark stars) and super-Keplerian cases (only for quark stars, which are bound additionally by strong interaction) have been discussed in this paper in order to investigate possible ways of forming sub-millisecond pulsars. In the sub-Keplerian case, the equilibrium periods of both neutron and quark stars could be as low as ~0.4 ms when they are spun up through accretion in binary systems. In the super-Keplerian case, pulsars could very likely have an initial period of ~0.1 ms if quark stars with different mass could be formed from accretion-induced collapse (AIC) of white dwarfs. The timescale for a sub-millisecond pulsar to keep its period <1 ms is restricted due to gravitational wave radiation. We found that the timescales of neutron stars within sub-millisecond periods are approximately ~10 yr, but the timescales estimated (>10 Myr) for low-mass quark stars could be long enough for us to detect.

{Magnetic fields exist on all scales in our Galaxy. There is a controversy about whether the magnetic fields in molecular clouds are preserved from the permeated magnetic fields in the interstellar medium (ISM) during cloud formation. We investigate this controversy using available data in the light of the newly revealed magnetic field structure of the Galactic disk obtained from pulsar rotation measures (RMs).} % {We collected measurements of the magnetic fields in molecular clouds, including Zeeman splitting data of OH masers in clouds and OH or HI absorption or emission lines of clouds themselves.} % {The Zeeman data show structures in the sign distribution of the line-of-sight component of the magnetic field. Compared to the large-scale Galactic magnetic fields derived from pulsar RMs, the sign distribution of the Zeeman data shows similar large-scale field reversals. Previous such examinations were flawed in the over-simplified global model used for the large-scale magnetic fields in the Galactic disk.} % {We conclude that the magnetic fields in the clouds may still ``remember'' the directions of magnetic fields in the Galactic ISM to some extent, and could be used as complementary tracers of the large-scale magnetic structure. More Zeeman data of OH masers in widely distributed clouds are required.}

We present new measurements of heavy cosmic-ray nuclei at high energies per- formed during the first flight of the balloon-borne cosmic-ray experiment CREAM (Cosmic-Ray Energetics And Mass). This instrument uses multiple charge detectors and a transition radiation detector to provide the first high accuracy measurements of the relative abundances of elements from boron to oxygen up to energies around 1 TeV/n. The data agree with previous measurements at lower energies and show a relatively steep decline (~E$^-0.6$ to E$^-0.5$) at high energies. They further show the source abundance of nitrogen relative to oxygen is ~10% in the TeV/n region.

Evidence for a spectral flattening of the supernova remnant (SNR) HB3 (G132.7+1.3) was recently claimed in literature based on previously published total flux density data, and the flattening was further interpreted as the discovery of thermal bremsstrahlung emission in the shell of HB3. A spectral flattening has never been observed from any SNR before. Reliable observations of HB3 at frequencies above 3000 MHz are crucial to confirm such a spectral behaviour. We extracted 4800 MHz total intensity and polarisation data of HB3 from the Sino-German 6 cm polarisation survey of the Galactic plane made with the Urumqi 25 m telescope, and analysed the spectrum of HB3, together with Effelsberg data at 1408 MHz and 2675 MHz. We found an overall spectral index of HB3 of alpha=-0.61+-0.06 between 1408 MHz and 4800 MHz, similar to the index at lower frequencies. There is no spectral flattening at high frequencies. We detected strong polarised emission from HB3 at 4800 MHz. Our 4800 MHz data show a tangential field orientation in the HB3 shell.

Our purpose is to place firm observational constraints on the three most widely used theoretical models for the spatial configuration of the large-scale interstellar magnetic field in the Galactic disk, namely, the ring, the axisymmetric and the bisymmetric field models. We use the rotation measures (RMs) of low-latitude Galactic pulsars and combine them with their dispersion measures and estimated distances to map out the line-of-sight component of the interstellar magnetic field in the near half of the Galactic disk. We then fit our map of the line-of-sight field to the three aforementioned theoretical field models and discuss the acceptability of each fit, in order to determine whether the considered field model is allowed by the pulsar data or not. Strictly speaking, we find that all three field models are ruled out by the pulsar data. Furthermore, none of them appears to perform significantly better than the others. From this we conclude that the large-scale interstellar magnetic field in the Galactic disk has a more complex pattern than just circular, axisymmetric or bisymmetric.

Binary pulsar systems are superb probes of stellar and binary evolution and the physics of extreme environments. In a survey with the Arecibo telescope, we have found PSR J1903+0327, a radio pulsar with a rotational period of 2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass companion. Infrared observations identify a possible main-sequence companion star. Conventional binary stellar evolution models predict neither large orbital eccentricities nor main-sequence companions around millisecond pulsars. Alternative formation scenarios involve recycling a neutron star in a globular cluster then ejecting it into the Galactic disk or membership in a hierarchical triple system. A relativistic analysis of timing observations of the pulsar finds its mass to be 1.74+/-0.04 Msun, an unusually high value.

A sample of classical IntraDay Variable (IDV) and IDV candidate sources has been monitored with the Urumqi 25m telescope and the Effelsberg 100m telescope. Aim of the project is to investigate the origin of IntraDay Variability, a phenomenon which has been observed in about 30% of flat spectrum radio quasars. Simultaneous Effelsberg-Urumqi observations demonstrated that the Urumqi antenna, although relatively small in diameter, is well suitable for IDV experiments. A few Urumqi datasets, however, turned out to be affected by a spurious $\sim 24$ hours modulation, an effect which has been removed by means of a new procedure for data reduction. In about 14 months, 12 epochs of observation have been collected, for a total observing time of more than 45 days. The epochs are regularly distributed over the whole year, in order to check for the presence of systematic annual changes in the variability time scales - a crucial test for verifying the consistency of source-extrinsic models of the variability. Preliminary time-analysis of the monitored sources revealed some hint for a slowing down of the characteristic time scales of S5~0716+714, a result that, if confirmed, would be compatible with a source-extrinsic origin of the variability, in contrast to previous IDV studies. No significant modulation of the time scales has been detected for S4~0954+658.

(Abridged) S147 is a large faint shell-type supernova remnant (SNR). Its remarkable spectral break at cm-wavelengths is an important physical property to characterize the SNR evolution. However, the spectral break is based on radio observations with limited precision. We made new radio continuum and polarization observations of S147 at 11cm and at 6cm with the Effelsberg 100-m telescope and the Urumqi 25-m telescope, respectively. These new data were combined with published lower frequency data from the Effelsberg 100-m telescope and very high frequency data from WMAP to investigate the spectral turnover and polarization properties of S147. S147 consists of numerous filaments embedded in diffuse emission. We found that the integrated flux densities of S147 are 34.8+/-4.0 Jy at 11cm and 15.4+/-3.0Jy at 6cm. These new measurements confirm the known spectral turnover at ~1.5GHz, which can be entirely attributed to the diffuse emission component. The spectral index above the turnover is -1.35+/-0.20. The filamentary emission component has a constant spectral index over the entire wavelength range up to 40.7GHz of -0.35+/-0.15. The weak polarized emission of S147 is at the same level as the ambient diffuse Galactic polarization. The rotation measure of the eastern filamentary shell is about -70 rad/m2. The filamentary and diffuse emission components of S147 have different physical properties, which make S147 outstanding among shell type SNRs.The weak polarization of S147 at 11cm and at 6cm can be attributed to a section of the S147 shell showing a tangential magnetic field direction.

Variations in the solar wind density introduce variable delays into pulsar timing observations. Current pulsar timing analysis programs only implement simple models of the solar wind, which not only limit the timing accuracy, but can also affect measurements of pulsar rotational, astrometric and orbital parameters. We describe a new model of the solar wind electron density content which uses observations from the Wilcox Solar Observatory of the solar magnetic field. We have implemented this model into the tempo2 pulsar timing package. We show that this model is more accurate than previous models and that these corrections are necessary for high precision pulsar timing applications.

Magnetic fields have been observed on all scales in our Galaxy, from AU to kpc. With pulsar dispersion measures and rotation measures, we can directly measure the magnetic fields in a very large region of the Galactic disk. The results show that the large-scale magnetic fields are aligned with the spiral arms but reverse their directions many times from the inner-most arm (Norma) to the outer arm (Perseus). The Zeeman splitting measurements of masers in HII regions or star-formation regions not only show the structured fields inside clouds, but also have a clear pattern in the global Galactic distribution of all measured clouds which indicates the possible connection of the large-scale and small-scale magnetic fields.

G156.2+5.7 is a large supernova remnant (SNR) first discovered in the ROSAT X-ray survey but radio observations are rare because of its low surface brightness and large size. Radio continuum and polarization images of the SNR G156.2+5.7 at lambda 6 cm are presented for the first time to be discussed in the context with available data to investigate its physical properties. We obtained an integrated flux density for G156.2+5.7 of S6cm= 2.5+-0.5Jy. The spectral index of the integrated emission is alpha=-0.48+-0.08 between lambda 74 cm and 6 cm excluding compact radio sources. We also obtained a spectral index map, which shows little variation between the shells and the central area of the SNR. This is consistent with the constant integrated flux density spectrum. Highly polarized radio emission has been detected from the SNR shell, but also from a central patch, which probably originates in the front part of the spherical SNR shell. We derived the distribution of rotation measure from polarization data at 11 cm and 6 cm and found RM gradients of opposite direction in the SNR shell. The SNR G156.2+5.7 is unusual by its exceptionally high X-ray brightness and very low surface-brightness in the radio range. The magnetic field is very well ordered along the shell periphery as expected for a compressed ambient magnetic field. A toroidal magnetic field component is indicated by the RM distribution. (abridged)

Short time-scale radio variations of compact extragalactic radio quasars and blazars known as IntraDay Variability (IDV) can be explained in at least some sources as a propagation effect; the variations are interpreted as scintillation of radio waves in the turbulent interstellar medium of the Milky Way. One of the most convincing observational arguments in favor of a propagation-induced variability scenario is the observed annual modulation in the characteristic time scale of the variation due to the Earth's orbital motion. So far there are only two sources known with a well-constrained seasonal cycle. Annual modulation has been proposed for a few other less well-documented objects. However, for some other IDV sources source-intrinsic structural variations which cause drastic changes in the variability time scale were also suggested. J1128+592 is a recently discovered, highly variable IDV source. Previous, densely time-sampled flux-density measurements with the Effelsberg 100-m radio telescope (Germany) and the Urumqi 25-m radio telescope (China), strongly indicate an annual modulation of the time scale. The most recent 4 observations in 2006/7, however, do not fit well to the annual modulation model proposed before. In this paper, we investigate a possible explanation of this discrepancy.

Short time-scale radio variations of compact extragalactic radio sources, known as IntraDay Variability, can be explained in at least some sources by a source-extrinsic effect, in which the variations are interpreted as scintillation of radio waves caused by the turbulent ISM of the Milky Way. One of the most convincing observational arguments in favour of propagation-induced variability is the so called annual modulation of the characteristic variability time-scale, which is due to the orbital motion of the Earth. Data for the recently discovered and highly variable IDV source J1128+5925 are presented. We study the frequency and time dependence of the IDV in this compact quasar. We measure the characteristic variability time-scale of the IDV throughout the year, and analyze whether the observed changes in the variability time-scale are consistent with annual modulation. We monitored the flux density variability of J1128+5925 with dense time sampling between 2.7 and 10.45GHz with the 100m Effelsberg radio telescope of the MPIfR and with the 25m Urumqi radio telescope. From ten observing sessions, we determine the variability characteristics and time-scales. The observed pronounced changes of the variability time-scale of J1128+5925 are modelled with an anisotropic annual modulation model. The observed frequency dependence of the variation is in good agreement with the prediction from interstellar scintillation. Adopting a simple model for the annual modulation model and using also the frequency dependence of the IDV, we derive a lower limit to the distance of the scattering screen and an upper limit to the scintillating source size. The latter is found to be consistent with the measured core size from VLBI.

Because of the contamination of jets, using the size -- continuum luminosity relation can overestimate the broad line region (BLR) size and black hole mass for radio-loud AGNs. We propose a new relation between the BLR size and $H_{\beta}$ emission line luminosity and present evidences for using it to get more accurate black hole masses of radio-loud AGNs. For extremely radio-loud AGNs such as blazars with weak/absent emission lines, we suggest to use the fundamental plane relation of their elliptical host galaxies to estimate the central velocity dispersions and black hole masses, if their host galaxies can be mapped.

(Abridged) Polarization measurements of the Galactic plane at lambda6 cm probe the interstellar medium (ISM) to larger distances compared to measurements at longer wavelengths, hence enable us to investigate properties of the Galactic magnetic fields and electron density. We are conducting a new lambda6 cm continuum and polarization survey of the Galactic plane covering 10degr < l < 230degr and |b|<5degr. Missing large-scale structures in the U and Q maps are restored based on extrapolated polarization K-band maps from the WMAP satellite. The lambda6 cm data are analyzed together with maps at other bands. We discuss some results for the first survey region, 7degr X 10degr in size, centered at (l,b)=(125.5degr, 0degr). Two new passive Faraday screens, G125.6-1.8 and G124.9+0.1, were detected. They cause significant rotation of background polarization angles but little depolarization. G124.9+0.1 was identified as a new faint HII region at a distance of 2.8 kpc. G125.6-1.8, with a size of about 46 pc, has neither correspondence in enhanced Halpha emission nor a counterpart in total intensity. A model combining foreground and background polarization modulated by the Faraday screen was developed. Using this model, we estimated the strength of the ordered magnetic field along the line of sight to be 3.9 microGauss for G124.9+0.1, and exceeding 6.4 microGauss for G125.6-1.8. We obtained an estimate of 2.5 and 6.3 mK/kpc for the average polarized and total synchrotron emissivity towards G124.9+0.1. The synchrotron emission beyond the Perseus arm is quite weak.

Recent observations of pulsar wind nebulae and radio polarization profiles revealed a tendency of the alignment between the spin and velocity directions in neutron stars. We study the condition for spin-kick alignment using a toy model, in which the kick consists of many off-centered, randomly-oriented thrusts. Both analytical considerations and numerical simulations indicate that spin-kick alignment cannot be easily achieved if the proto-neutron star does not possess some initial angular momentum, contrary to some previous claims. To obtain the observed spin-kick misalignment angle distribution, the initial spin period of the neutron star must be smaller than the kick timescale. Typically, an initial period of a hundred milliseconds or less is required.

The relationship between the emission-line equivalent width and the continuum luminosity, so called the Baldwin effect, plays an important role in studying the physics of the broad line region of AGNs. Using the archived ultraviolet spectra obtained by IUE, HST, and HUT in 1978-2002, we investigated the intrinsic CIV Baldwin effect of a well-studied Seyfert 1 galaxy NGC 4151. Both its continuum flux and CIV emission-line flux varied about two orders of magnitude in more than two decades, making it one of the best targets for studying the slope variation of the Baldwin effect. We fitted the CIV line profile of the 490 archived UV spectra of NGC 4151 with the same model consisting of a few Gaussian components, and derived the slope in the log-log plot for the total flux of CIV emission line against the UV continuum flux in different observation epochs.We found that the slope is not constant for NGC 4151, as it varies from 0.58 in the highest flux epoch to 0.83 in the lowest flux epoch. The slope evidently decreases as the continuum flux increases, which reinforces the previous findings of the non-constant slope in the H_beta Baldwin effect of NGC 5548 and the CIV Baldwin effect of Fairall 9. Our result suggests that such a non-constant slope may not be unusual for AGNs. Its physical origin is probably related to the different non-linear responses of the emission line to the variable ionizing continuum caused by the different accretion modes at different luminosity levels. We briefly discuss the effects of various absorption components in the CIV line profile of NGC 4151 and argue that the slope variation is not driven mainly by the absorption effect. Intensive, long-term and high resolution spectral observations of strongly variable AGNs are needed in the future to confirm our result.

The determination of the central black hole mass is crucial to the understanding of AGN physics. In this paper we briefly review some methods that are currently used to estimate the black hole mass of AGNs. Particularly we demonstrate the importance of two correlations: one between the black hole mass and the stellar velocity dispersion and another one between the broad line region (BLR) size and the optical continuum luminosity. Besides applying these relations in deriving black hole masses of various types of AGNs, we also employed the fundamental plane of elliptical galaxies to estimate the central velocity dispersions of AGN host galaxies, and then the black hole masses of AGNs including BL Lac objects. In addition, we derived another empirical relation between the BLR size and H_beta emission line luminosity from AGNs with the BLR size measured by reverberation mapping studies, and argued that more accurate black hole masses of extremely radio-loud AGNs could be obtained with it than using the usual R-L_{5100\AA} relation because of the jet contribution to the optical continuum. Finally we pointed out that black hole mass estimation is very much helpful to determine the accretion rate and understand the accretion process in AGNs.

The recently started Arecibo L-band Feed Array (ALFA) pulsar survey aims to find ~1000 new pulsars. Due to its high time and frequency resolution the survey is especially sensitive to millisecond pulsars, which have the potential to test gravitational theories, detect gravitational waves and probe the neutron-star equation of state. Here we report the results of our preliminary analysis: in the first months we have discovered 21 new pulsars. One of these, PSR J1906+0746, is a young 144-ms pulsar in a highly relativistic 3.98-hr low-eccentricity orbit. The 2.61 +- 0.02 solar-mass system is expected to coalesce in ~300 Myr and contributes significantly to the computed cosmic inspiral rate of compact binary systems.

A decade ago, there was very limited knowledge of magnetic fields of our Galaxy. The local fields in the Solar vicinity were known to be directed towards a Galactic longitude $l\sim~90^{\circ}$ with reversed directions at smaller Galacto-radii. The regular field strength was found to be about 2~$\mu$G. The filaments near the Galactic Center show the possible poloidal fields there. There was no information about the magnetic fields in the Galactic halo. In last decade, there has been significant progress on measurements of the Galactic magnetic fields. In the Galactic disk, from the RMs of a large number of newly observed pulsars, large-scale magnetic fields along the spiral arms have been delineated in a much larger region then ever before, with alternating directions in the arm and interarm regions. The toroidal fields in the Galactic halo were revealed to have opposite directions below and above the Galactic plane, which is an indication of A0 dynamo operating in the halo. The strength of large-scale fields also has been found from pulsar RM data to exponentially increase at smaller Galacto-radii. Compared to the steep Kolmogorov spectrum of magnetic energy at small scales, the large-scale magnetic fields show a shallow broken spatial magnetic energy spectrum.

We update the systematic studies of circular polarization in integrated pulse profiles by Han et al (1998). Data of circular polarization profiles are compiled. Sense reversals can occur in core or cone components, or near the intersection between components. The correlation between the sense of circular polarization and the sense of position angle variation for conal-double pulsars is confirmed with a much large database. Circular polarization of some pulsars has clear changes with frequency. Circular polarization of millisecond pulsars is marginally different from that of normal pulsars.

The large-scale magnetic field of our Galaxy can be probed in three dimensions using Faraday rotation of pulsar signals. We report on the determination of 223 rotation measures from polarization observations of relatively distant southern pulsars made using the Parkes radio telescope. Combined with previously published observations these data give clear evidence for large-scale counterclockwise fields (viewed from the north Galactic pole) in the spiral arms interior to the Sun and weaker evidence for a counterclockwise field in the Perseus arm. However, in interarm regions, including the Solar neighbourhood, we present evidence that suggests that large-scale fields are clockwise. We propose that the large-scale Galactic magnetic field has a bisymmetric structure with reversals on the boundaries of the spiral arms. Streaming motions associated with spiral density waves can directly generate such a structure from an initial inwardly directed radial field. Large-scale fields increase toward the Galactic Center, with a mean value of about 2~$\mu$G in the Solar neighbourhood and 4~$\mu$G at a Galactocentric radius of 3 kpc.

Based on the measured sizes of broad line region of the reverberation-mapping AGN sample, two new empirical relations are introduced to estimate the central black hole masses of radio-loud high-redshift ($z > 0.5$) AGNs. First, using the archival $IUE/HST$ spectroscopy data at UV band for the reverberation-mapping objects, we obtained two new empirical relations between the BLR size and \Mg/\C emission line luminosity. Secondly, using the newly determined black hole masses of the reverberation-mapping sample for calibration, two new relationships for determination of black hole mass with the full width of half maximum and the luminosity of \Mg/\C line are also found. We then apply the relations to estimate the black hole masses of AGNs in Large Bright Quasar Surveyq and a sample of radio-loud quasars. For the objects with small radio-loudness, the black hole mass estimated using the $R_{\rm BLR} - L_{\eMg/\eC}$ relation is consistent with that from the $R_{BLR} - L_{3000\AA/1350 \AA}$ relation. But for radio-loud AGNs, the mass estimated from the $R_{BLR} - L_{\eMg/\eC}$ relation is systematically lower than that from the continuum luminosity $L_{3000\AA/1350\AA}$. Because jets could have significant contributions to the UV/optical continuum luminosity of radio-loud AGNs, we emphasized again that for radio-loud AGNs, the emission line luminosity may be a better tracer of the ionizing luminosity than the continuum luminosity, so that the relations between the BLR size and UV emission line luminosity should be used to estimate the black hole masses of high redshift radio-loud AGNs.

We report the discovery of PSR J1906+0746, a young 144-ms pulsar in a highly relativistic 3.98-hr orbit with an eccentricity of 0.085 and expected gravitational wave coalescence time of 300 Myr. The new pulsar was found during precursor survey observations with the Arecibo 1.4-GHz feed array system and retrospectively detected in the Parkes Multibeam plane pulsar survey data. From radio follow-up observations with Arecibo, Jodrell Bank, Green Bank, and Parkes, we have measured the spin-down and binary parameters of the pulsar and its basic spectral and polarization properties. We also present evidence for pulse profile evolution, which is likely due to geodetic precession, a relativistic effect caused by the misalignment of the pulsar spin and total angular momentum vectors. Our measurements show that PSR J1906+0746 is a young object with a characteristic age of 112 kyr. From the measured rate of orbital periastron advance 7.57+/-0.03 deg/yr, we infer a total system mass of 2.61+/-0.02 Msun. While these parameters suggest that the PSR J1906+0746 binary system might be a younger version of the double pulsar system, intensive searches for radio pulses from the companion have so far been unsuccessful. It is therefore not known whether the companion is another neutron star or a massive white dwarf. Regardless of the nature of the companion, a simple calculation suggests that the Galactic birth rate of binaries similar to PSR J1906+0746 is 60/Myr. This implies that PSR J1906+0746 will make a significant contribution to the computed cosmic inspiral rate of compact binary systems.

We study observational constraints on neutron star (NS) kicks for isolated pulsars and for neutron stars in binary systems. We are particularly interested in the evidence of kick-spin alignment/misalignment and its dependence on the neutron star initial spin period. For several young pulsars, X-ray observations of compact nebulae showed that pulsar proper motion is aligned with the spin direction as defined by the symmetry axis of the nebula. We also critically examine the measurements of the proper motion and the projected spin axis from a large sample of pulsars with well-calibrated polarization data. We find that among the two dozen pulsars for which reliable measurements are available, there is a significant correlation between the spin axis and the proper motion. For various NS binaries, including double NS systems, binaries with massive main-sequence star companion and binaries with massive white-dwarf companion, we obtain constraints on the kick magnitudes and directions from the observed orbital characteristics of the system. We find that the kick velocity is misaligned with the the NS spin axis in a number of systems, and the NS spin period (when available) in these systems is generally longer than several hundreds milliseconds. These constraints, together with spin-kick alignment observed in many isolated pulsars, suggest that the kick timescale is hundreds of milliseconds to 1 s, so that spin-kick alignment or misalignment can be obtained depending on the initial spin period of the NS. We discuss the implication of our result for various NS kick mechanisms.

Radio continuum and polarization observations of the entire Cygnus Loop at 6cm wavelength were made with the Urumqi 25m telescope. The 6cm map is analysed together with recently published maps from the Effelsberg 100m telescope at 21cm and 11cm. The integrated flux density of the Cygnus Loop at 6cm is 90+/-9Jy, which implies a spectral index of -0.40+/-0.06 being consistent with that of Uyaniker et al. (2004) in the wavelength range up to 11cm. This rules out any global spectral steepening up to 6cm. However, small spectral index variations in some regions of the source are possible, but there are no indications for any spectral curvature. The linear polarization data at 6cm show percentage polarizations up to 35% in some areas of the Cygnus Loop, exceeding those observed at 11cm. The Rotation Measure is around -21rad/m^2 in the southern area, which agrees with previous observations. However, the distribution of Rotation Measures is rather complex in the northern part of the Cygnus Loop, where the 21cm emission is totally depolarized. Rotation Measures based on 11cm and 6cm data are significantly larger than in the southern part. The difference in the polarization characteristic between the northern and southern part supports previous ideas that the Cygnus Loop consists of two supernova remnants.

We report results from the initial stage of a long-term pulsar survey of the Galactic plane using the Arecibo L-band Feed Array (ALFA), a seven-beam receiver operating at 1.4 GHz with 0.3 GHz bandwidth. The search targets Galactic latitudes |b| < 5 deg in the longitude ranges 32 deg < l < 77 deg and 168 deg < l < 77 deg. Data discussed here were collected over a 100 MHz passband centered on 1.42 GHz using a spectrometer that recorded 256 channels every 64 microsec. In a preliminary, standard period-DM analysis, we have detected 29 previously known pulsars and discovered 11 new ones. One of these, with a period of 69 ms and a low characteristic age of 82 kyr, is a plausible candidate for association with the unidentified EGRET source 3EG J1928+1733. Another is a non-recycled pulsar in a relativistic binary with orbital period of 3.98 hr. We also search the data for isolated dispersed pulses, a technique that yielded discovery of an extremely sporadic radio emitter with a spin period of 1.2 s. Simulations we have carried out indicate that about 1000 new pulsars will be found in the ALFA survey. In addition to providing a large sample for use in population analyses and for probing the magnetoionic interstellar medium, the survey maximizes the chances of finding rapidly spinning millisecond pulsars and pulsars in compact binary systems. Our search algorithms will exploit the multiple data streams from ALFA to discriminate between radio frequency interference and celestial signals, including pulsars and possibly new classes of transient radio sources.

Interstellar magnetic fields exist over a broad range of spatial scales, extending from the large Galactic scales ($\sim 10$ kpc) down to the very small dissipative scales ($\ll 1$ pc). In this paper, we use a set of 490 pulsars distributed over roughly one third of the Galactic disk out to a radius $R \simeq 10$ kpc (assuming $R_\odot = 8.5$ kpc) and combine their observed rotation and dispersion measures with their estimated distances to derive the spatial energy spectrum of the Galactic interstellar magnetic field over the scale range $0.5 - 15$ kpc. We obtain a nearly flat spectrum, with a 1D power-law index $\alpha=-0.37\pm0.10$ for $E_{\rm B}(k)=C k^{\alpha}$ and an rms field strength of approximately $6 \mu$G over the relevant scales. Our study complements the derivation of the magnetic energy spectrum over the scale range $0.03 - 100$ pc by \citet{ms96b}, showing that the magnetic spectrum becomes flatter at larger scales. This observational result is discussed in the framework of current theoretical and numerical models.

The inner annular gap (IAG), a new type of inner gap whose magnetic field lines intersect the null charge surface (NCS), is proposed to explain $\gamma$-ray and radio emission from pulsars. The IAG can be an important source for high-energy particles. The particles can radiate between the NCS and the IAG. Some observational characteristics in both $\gamma$-ray and radio bands, such as the $\gamma$-ray emission beams of Crab-like, Vela-like and Geminga-like, can be reproduced by numerical method. It is predicted that the view angle $\zeta$ should be larger than the inclination angle ($\zeta>\alpha$), otherwise the $\gamma$-ray radiation will have little possibility to be observed. Whether the inner annular gap (or cap) is sparking (or free flow) depends on the surface binding energy of the pulsar. In stead of neutron star models, the scenario of the IAG is favorable for bare strange star models, because bare strange stars can easily satisfy the requisite condition to form an IAG for both pulsars ($\vec{\Omega}\cdot \vec{B}<0$) and anti-pulsars ($\vec\Omega \cdot \vec{B}>0$).

An empirical relation between the broad line region (BLR) size and optical continuum luminosity is often adopted to estimate the BLR size and then the black hole mass of AGNs. However, optical luminosity may not be a good indicator of photoionizing luminosity for extremely radio-loud AGNs because the jets usually contribute significantly to the optical continuum. Therefore, the black hole masses derived for blazar-type AGNs with this method are probably overestimated. Here we first derived a tight empirical relation between the BLR size and the H$_\beta$ emission line luminosity, $R(\rm{light-days})= 24.05(L_{H_\beta}/10^{42} ergs s^{-1})^{0.68}$, from a sample of 34 AGNs with the BLR size estimated with the reverberation mapping technique. Then we applied this relation to estimate the black hole masses of some AGNs and found that for many extremely radio-loud AGNs the black hole masses obtained with the $R-L_{H_\beta}$ relation are systematically lower than those derived previously with the $R-L_{5100\AA}$ relation, while for radio-quiet and slightly radio-loud AGNs the results obtained with these two methods are almost the same. The difference of black hole masses estimated with these two relations increases with the radio-loudness for extremely radio-loud AGNs, which is consistent with the fact that their equivalent widths of H$_{\beta}$ emission line become smaller at higher radio-loudness. If the small H$_{\beta}$ equivalent widths of extremely radio-loud AGNs are indeed caused by the beaming effect, we argue that the optical emission line luminosity may be a better tracer of ionizing luminosity for blazar-type AGNs and the black hole mass derived with the $R-L_{H_\beta}$ relation are probably more accurate.

Mean pulse profiles and polarization properties are presented for nine southern pulsars. The observations were made using the Parkes radio telescope at frequencies near 1330 MHz; three of the nine pulsars were also observed at 660 MHz. A very high degree of circular polarization was observed in PSR J1603-7202. Complex position angle variations which are not well described by the rotating-vector model were observed in PSRs J2124-3358 and J2145-0750, both of which have very wide profiles. Rotation measures were obtained for all nine pulsars, with two implying relatively strong interstellar magnetic fields.

Pulsar motions in our Galaxy from their birth until 2 Gyr are studied statistically via Monte-Carlo simulation of $2\times 10^5$ pulsars with the best available representation of the Galactic potential. We find that the distribution of height above the Galactic plane for pulsars with characteristic ages less than about 8 Myr could be well fitted by a Gaussian function. For older pulsars, an extra exponential function is necessary to fit the distribution. The scale-height of the Gaussian component increases linearly with time until about 40 Myr. The height distribution becomes stabilized after about 200 Myr. These results are not sensitive to initial height or radial distributions. Taking the relationship between the initial velocity and height distribution, we found from the latest pulsar catalog that the height distribution of pulsars younger than 1 Myr directly implies the mean initial velocity of $280\pm96$ km s$^{-1}$. Comparison of simulated sample of pulsars with the current available millisecond pulsars shows that their 1D initial velocity dispersion should be most probably $60\pm10$ km s$^{-1}$.

We study structure functions of rotation measures in the Canadian Galactic Plane Survey (CGPS) region and the North Galactic Pole (NGP) to extract the interstellar medium (ISM) fluctuation information. The CGPS data are divided into three longitude intervals: $82^\circ<l<96^\circ$ (CGPS1), $115^\circ<l<130^\circ$ (CGPS2) and $130^\circ<l<146^\circ$ (CGPS3). The structure functions of all three regions have large uncertainties when the angular separation is smaller than $\delta\theta\approx 1^\circ$. A power law can fit the structure function well for $\delta\theta>1^\circ$. The power law indices get smaller from CGPS1 to CGPS3 and the amplitudes decrease. The variations of the large-scale field and the electron density have only negligible effects on the structure function and thus cannot account for the changes, indicating that the turbulent properties of the Galactic ISM are intrinsically longitude-dependent. The Kolmogorov-like fluctuation spectrum of the electron density or the magnetic field should produce a power law structure function with an index of 5/3 or 2/3, neither of which is consistent with our results of small indices in the three sub-CGPS regions. For the NGP region, the structure function is flat, showing that the rotation measures are mostly intrinsic to the extragalactic sources, and the ISM is very random in that part of our Galaxy. It is obvious that the ISM fluctuation is latitude-dependent when comparing the results in the NGP region and the CGPS regions.

In this review, I will first introduce possible methods to probe the large-scale magnetic fields in our Galaxy and discuss their limitations. The magnetic fields in the Galactic halo, mainly revealed by the sky distribution of rotation measures of extragalactic radio sources, probably have a global structure of a twisted dipole field. The large-scale field structure in the Galactic disk has been most efficiently deduced from pulsar rotation measures (RMs). There has been a lot of progress since the 1980s when the magnetic field in the local area of our Galaxy was first traced by a very small sample of local pulsars. Now we have pulsars distributed in about one third of the whole Galactic disk of the interstellar medium, which shows that the large-scale magnetic fields go along the spiral arms and that the field directions reverse from arm to arm. The RMs of newly discovered pulsars in the very inner Galaxy have been used to show the coherent magnetic field in the Norma arm. The magnetic fields in the Galactic disk most probably have a bisymmetric spiral structure.

Pulsars are among the most highly polarized sources in the universe. The NVSS has catalogued 2 million radio sources with linear polarization measurements, from which we have selected 253 sources, with polarization percentage greater than 25%, as targets for pulsar searches. We believe that such a sample is not biased by selection effects against ultra-short spin or orbit periods. Using the Parkes 64m telescope, we conducted searches with sample intervals of 0.05 ms and 0.08 ms, sensitive to submillisecond pulsars. Unfortunately we did not find any new pulsars.

The low-frequency profiles of some pulsars manifest temporal broadening due to scattering, usually accompanied by flat polarization position angle (PA) curves. Assuming that the scattering works on the 4 Stokes parameters in the same way, we have simulated the effect of scattering on polarization profiles and find that the scattering can indeed flatten the PA curves. Since the higher-frequency profiles suffer less from scattering, they are convolved with scattering models to fit the observed low-frequency profiles. The calculated flat PA curves exactly reproduce the corresponding observations.

Pulsars provide unique probes for the {\it large-scale} interstellar magnetic field in the Galactic disk. The Parkes multibeam pulsar survey has discovered many distant pulsars which enables us for the first time to explore the magnetic field in most of the nearby half of the Galactic disk. The fields are found to be {\it coherent in direction} over a linear scale of $\sim 10$ kpc between the Carina-Sagittarius and Crux-Scutum arms from $l\sim45\degr$ to $l\sim305\degr$. The coherent spiral structures and field direction reversals, including the newly determined counterclockwise field near the Norma arm, are consistent with bi-symmetric spiral model for the disk field. However, the antisymmetric rotation measure sky from the Galactic halo and the dipole field in the Galactic center suggest that the A0 dynamo is operating there.

Magnetic fields are observed everywhere in the universe. In this review, we concentrate on the observational aspects of the magnetic fields of Galactic and extragalactic objects. Readers can follow the milestones in the observations of cosmic magnetic fields obtained from the most important tracers of magnetic fields, namely, the star-light polarization, the Zeeman effect, the rotation measures (RMs, hereafter) of extragalactic radio sources, the pulsar RMs, radio polarization observations, as well as the newly implemented sub-mm and mm polarization capabilities. (Another long paragraph is omitted due to the limited space here)

Pulsars provide unique probes of the large-scale interstellar magnetic field in the Galactic disk. Up to now, the limited Galactic distribution of the known pulsar population has restricted these investigations to within a few kiloparsec of the Sun. The Parkes multibeam pulsar survey has discovered many more-distant pulsars which enables us for the first time to explore the magnetic field in most of the nearby half of the Galactic disk. Here we report the detection of counterclockwise magnetic fields in the Norma spiral arm using pulsar rotation measures. The fields are coherent in direction over a linear scale of $\sim 5$ kpc along the arm and have a strength of $-4.4\pm0.9$ $\mu$G. The magnetic field between the Carina-Sagittarius and Crux-Scutum arms is confirmed to be coherent from $l\sim45\degr$ to $l\sim305\degr$ over a length of $\sim 10$ kpc. These results strengthen arguments for a bisymmetric spiral model for the field configuration in the Galactic disk.

We identified 16 pulsars from the Faint Images of the Radio Sky at Twenty-cm (FIRST) at 1.4 GHz. Their positions and total flux densities are extracted from the FIRST catalog. Comparing the source positions with those in the PSRcatalog, we obtained better determined positions of PSRs J1022+1001, J1518+4904, J1652+2651, and proper motion upper limits of another three pulsars PSRs J0751+1807, J1012+5307, J1640+2224. Proper motions of the other 10 pulsars are consistent with the values in the catalog.

I review the large scale global magnetic field structure of our Galaxy, using all information available for disk fields, halo fields and magnetic fields near the Galactic center (GC). In the local disk of our Galaxy, RM and dispersion measure (DM) data of nearby pulsars yield the strength of regular field as 1.8$\mu$G, with a pitch angle of about $8\degr$, and a bisymmetric spiral structure. There are at least four, maybe five, field reversals from the Norma arm to the outskirts of our Galaxy. In the thick disk or Galactic halo, large scale toroidal magnetic fields, with opposite field directions in the Southern and Northern Galaxy, have been revealed by the antisymmetric RM sky towards the inner Galaxy. This signature of the A0 dynamo-mode field structure is strengthened by the indication of a poloidal field of dipole form, that is the transition of the RM signs probably shifted from $l\sim0\degr$ to $l\sim+10\degr$. The local vertical field is probably a part of this dipole field. The field structure of the A0 dynamo-mode strikingly continues towards the region near the GC. In short, the magnetic fields in the Galactic disk have a bisymmetric spiral structure of primordial nature, while in the halo and near the GC the A0 dynamo seems to dominate, so that the fields consist of toroidal fields with opposite directions below and above the Galactic plane and poloidal fields of dipole form.

We estimated black hole masses for 9 Seyfert 1 and 13 Seyfert 2 galaxies in the Palomar and CfA bright Seyfert samples using the tight correlation between black hole mass and bulge velocity dispersion. Combining other 13 Seyfert 1s and 2 Seyfert 2s in these samples but with black hole masses measured recently by reverberation mapping and stellar/gas dynamics, we studied the correlations of black hole masses with radio loudness and bulge luminosities for a sample of 37 Seyfert galaxies. We found that if radio-loudness is measured using the optical and radio luminosities of the nuclear components, the black hole masses of radio-loud Seyfert 1s tend to increase with the radio-loudness. The black hole masses of all Seyfert galaxies increase with the radio power, but Seyfert galaxies have larger radio powers than nearby galaxies with the same black hole masses. In addition, the correlation between black hole masses and bulge V-band luminosities for Seyfert galaxies is consistent with that found for quasars and normal galaxies. The combined sample of 37 Seyfert galaxies, 15 quasars and 30 normal galaxies suggests a possible universal nonlinear relation between black hole and bulge masses, $M(BH) \propto M(bulge)^{1.74\pm0.14}$, which is slightly steeper than that found recently by Laor (2001) for a smaller sample. This nonlinear relation is supported by a larger sample including 65 Seyfert galaxies. The different M(BH)/M(bulge) ratio for galaxies with different bulge luminosities or different black hole masses may be explained by this relation. These results are consistent with some theoretical implications and are important for understanding the nature of radio emissions and the formation and evolution of supermassive black holes and galaxies.

A tight correlation of black hole mass and central velocity dispersion has been found recently for both active and quiescent galaxies. By applying this correlation, we develop a simple method to derive the inclination angles for a sample of 11 Seyfert 1 galaxies that have both measured central velocity dispersions and black hole masses estimated by reverberation mapping. These angles, with a mean value of 36 degree that agrees well with the result obtained by fitting the iron K$\alpha$ lines of Seyfert 1s observed with ASCA, provide further support to the orientation-dependent unification scheme of AGN. A positive correlation of the inclinations with observed FWHMs of H$\beta$ line and a possible anti-correlation with the nuclear radio-loudness have been found. We conclude that more accurate knowledge on inclinations and broad line region dynamics is needed to improve the black hole mass determination of AGN with the reverberation mapping technique.

The shapes of pulse profiles, especially their variations with respect to observing frequencies, are very important to understand emission mechanisms of pulsars, while no previous attempt has been made in interpreting the complicated phenomenology. In this paper, we present theoretical simulations for the integrated pulse profiles and their frequency evolution within the framework of the inverse Compton scattering (ICS) model proposed by Qiao (1988) and Qiao & Lin (1998). Using the phase positions of the pulse components predicted by the ``beam-frequency figure'' of the ICS model, we present Gaussian fits to the multi-frequency pulse profiles for some pulsars. It is shown that the model can reproduce various types of the frequency evolution behaviors of pulse profiles observed.

The fundamental plane of spiral galaxies was searched from observational data, which can be approximately represented in observational parameters by $ L \propto V^2 R $, where $L$, $R$, and $V$ are the luminosity, the linear size of galactic disk, and the rotation velocity. This plane exists at all optical bands in our samples of more than 500 spiral galaxies in total. It is more fundamental than the relationship between any two of the parameters, and can reduce the residual of the Tully--Fisher $L$--$V$ relations by about 50\%. Noticed that the power index of $V$ is doubled from that of $R$, which implies that the total mass and the mass distribution of galaxies plays an important role in forming the fundamental plane. Involving the third parameter of galactic size has a strong physical implication on galaxy formation and the dark--mass distribution.

The large scale magnetic fields of our Galaxy have been mostly revealed by rotation measures ($RM$s) of pulsars and extragalactic radio sources. In the disk of our Galaxy, the average field strength over a few kpc scale is about 1.8 $\mu$G, while the total field, including the random fields on smaller scales, has a strength of about 5 $\mu$G. The local regular field, if it is part of the large scale field of a bisymmetric form, has a pitch angle of about $-8\degr$. There are at least three, and perhaps five, field reversals from the Norma arm to the outer skirt of our Galaxy.

Using all available multi-component radio pulse profiles for pulsars with medium to long periods and good polarisation data, we have constructed a two-dimensional image of the mean radio beamshape. This shows a peak near the centre of the beam but is otherwise relatively uniform with only mild enhancements in a few regions. This result supports the patchy-beam model for emission beams in which the mean beam shape represents the properties of the emission mechanism and observed pulse components result from emission sources distributed randomly across the beam.

Some characteristics of the inverse Compton scattering (ICS) model are reviewed. At least the following properties of radio pulsars can be reproduced in the model: core or central emission beam, one or two hollow emission cones, different emission heights of these components, diverse pulse profiles at various frequencies, linear and circular polarization features of core and cones.

Polarization profiles of several hundred pulsars have been published recently (eg. Gould & Lyne 1998, GL98; Weisberg et al. 1999). In this report, we summarize the characteristics of circular and high linear polarizations of pulsar profiles, based on all previously published data.

We carried out the Gaussian fitting to the profile of PSR B1237+25 and found that six components rather than five are necessary to make a good fit. In the central part, we found that the core emission is not filled pencil beam but is a small hollow cone. This implies that the impact angle could be $\beta<0.5^\circ$. The ``hollow core'' is in agreement with Inverse Compton Scattering model of radio pulsars.

We propose that the violent breakdown of the Ruderman-Sutherland type vacuum gap above the pulsar polar cap region may produce copious primary particles and make the charge density above the gap in excess of the Goldreich-Julian density. Such an excess residue of charge density may arise some re-acceleration regions above the gap due to the space-charge-limited flow mechanism. Such re-acceleration regions, if exist, may accelerate particles to carry away more spin-down energy, and may have some implications for pulsar $\gamma$-ray emission.

Qiao and his collaborators recently proposed an inverse Compton Scattering (ICS) model to explain radio emission of pulsars. In this paper, we investigate the polarization properties of pulsar emission in the model. First of all, using the lower frequency approximation, we derived the analytical amplitude of inverse Compton scattered wave of a single electron in strong magnetic field. We found that the out-going radio emission of a single relativistic electron scattering off the ``low frequency waves'' produced by gap-sparking should be linearly polarized and have no circular polarization at all. However, considering the coherency of the emission from a bunch of electrons, we found that the out-going radiation from the inner part of emission beam, i.e., that from the lower emission altitudes, prefers to have circular polarization. Computer simulations show that the polarization properties, such as the sense reversal of circular polarization near the pulse center, S-shape of position angle swing of the linear polarization, strong linear polarization in conal components, can be reproduced in the ICS model.

We identified 97 strong pulsars from the NRAO VLA Sky Survey (NVSS) at 1.4 GHz north of Dec(J2000) $>-40\degr$. The total flux density, linear polarization intensity and polarization angle (PA) of all pulsars are extracted from the NVSS catalog. The well-calibrated PA measurement of 5 pulsars can be used for absolute PA calibrations in other observations. Comparing the source positions with those in the pulsar catalog, we got the first measurement of the proper motion upper limit of PSR B0031-07, which is $\mu_{\alpha}\cos\delta = -102 \pm 74 mas/yr$ and $\mu_{\delta} = -105 \pm 78 mas/yr$.

We have obtained 63 rotation measures (RMs) from polarization observations of southern pulsars, of which 54 are new measurements and 3 are varied from previous values. The new pulsar RM data at high Galactic latitudes are mostly consistent with the antisymmetric RM distribution found previously. For the Galactic disc, evidence for a field reversal near the Perseus arm, and possibly another beyond it, is presented. Inside the Solar Circle, in addition to the two known field reversals in or near the Carina-Sagittartus arm and the Crux-Scutum arm, a further reversal in the Norma arm is tentatively identified. These reversals, together with the pitch angle derived from pulsar RM and stellar polarization distributions, are consistent with bisymmetric spiral (BSS) models for the large-scale magnetic field structure in the disc of our Galaxy. However, discrimination between models is complicated by the presence of smaller-scale irregularities in the magnetic field, as well as uncertainties in the theoretical modelling.

We present a systematic study of the circular polarization in pulsar integrated profiles, based on published polarization data. For core components, we find no significant correlation between the sense-change of circular polarization and the sense of linear position angle variation. Circular polarization is not restricted to core components and, in some cases, reversals of circular polarization sense are observed across the conal emission. In conal double profiles, the sense of circular polarization is found to be correlated with the sense of position-angle variation. Pulsars with a high degree of linear polarization often have one hand of circular polarization across the whole profile. For most pulsars, the sign of circular polarization is the same at 50-cm and 20-cm wavelength, and the degree of polarization is similar, albeit with a wide scatter. However, at least two cases of frequency-dependent sign reversals are known. This diverse behaviour may require more than one mechanism to generate circular polarization.

We observed 21 polarized background radio sources in the field of M31 at 1.365 GHz and 1.652 GHz, and determined their rotation measures (RMs). The RM data show that the regular magnetic field of the disk probably extends from about 5 kpc to 25 kpc from the center with similar structure. The RMs obtained from the polarized emission from M31 at $\lambda$6 cm and $\lambda$11 cm indicate that M31 might have a weak poloidal field in its tenuous halo. Observational features of the odd and even dynamo modes in a galaxy are discussed. An even mode (S0) dynamo may operate in M31.

The theory predicts that the effect of gravitational lensing by the matter associated with clusters of galaxies can magnify background sources, leading to an enhancement of source number density around foreground clusters of galaxies. We conduct a search for the associations of distant radio-bright quasars with Abell clusters using the 1-Jy and 2-Jy all-sky catalogs. Statistics turns to be very poor on the basis of the 1-Jy sample, which shows no correlations between the distant radio quasars with the foreground Abell clusters above $1\sigma$ level. However, an apparent association ($>1\sigma$) of the 2-Jy radio sources with the foreground Abell clusters has been detected on scale of $\sim20$ arcminutes. We point out that this enhancement is unlikely to be produced by the statistical gravitational lensing hypothesis utilizing the matter associated with a population of isolated clusters unless (1)their velocity dispersion is a few times larger than the presently adopted value and/or (2)the intrinsic counts of the radio-bright sources with flux have a steeper slope than the presently observed ones. This indicates that (1)the observed associations, if real, can be the integrated result of all the matter along the line of sight to the distant quasars, namely, the weak lensing effect of clusters of galaxies that trace large-scale structures of the universe, and/or (2)the number counts of the radio-bright sources have been seriously contaminated by lensing.