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

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

The Five-Hundred-Meter Aperture Spherical Radio Telescope(FAST) is the largest single-dish radio telescope in the world. In this paper, we make forecast on the FAST HI large scale structure survey by mock observations. We consider a drift scan survey with the L-band 19 beam receiver, which may be commensal with the pulsar search and Galactic HI survey. We also consider surveys at lower frequency, either using the current single feed wide band receiver, or a future multi-beam phased array feed (PAF) in the UHF band. We estimate the number density of detected HI galaxies and the measurement error in positions, the precision of the surveys are evaluated using both Fisher matrix and simulated observations. The measurement error in the HI galaxy power spectrum is estimated, and we find that the error is relatively large even at moderate redshifts, as the number of positively detected galaxies drops drastically with increasing redshift. However, good cosmological measurement could be obtained with the intensity mapping technique where the large scale HI distribution is measured without resolving individual galaxies. The figure of merit (FoM) for the dark energy equation of state with different observation times are estimated, we find that with the existing L-band multi-beam receiver, a good measurement of low redshift large scale structure can be obtained, which complements the existing optical surveys. With a PAF in the UHF band, the constraint can be much stronger, reaching the level of a dark energy task force (DETF) stage IV experiment.

We investigate the model of dark matter-dark energy (DM-DE) interaction with coupling strength proportional to the multiplication of dark sector densities with different power indices $Q = \gamma \rho_{\rm c}^{\alpha} \rho_{\rm d}^{\beta}$. We first investigate the modification of the cosmic expansion history, and then further develop the formalism to take into account the cosmological perturbations and dark matter temperature evolution. We then use the latest observational cosmology data, including cosmic microwave background (CMB) data, baryon acoustic oscillations (BAO) data, redshift-space distortion (RSD) data and Type Ia supernovae (SNe) data to constrain the model parameters. We find in the phantom region, a positive $\alpha$ is preferred by the data above $2\, \sigma$ statistic significance. If we choose the power indices to be integers or half-integers for {\it plausible} physics of particle interaction, the allowed values within $1\, \sigma$ confidence regions are $\alpha = 0.5$ and $\beta = 0, 0.5, 1$. The inclusion of BAO and RSD data from large-scale structure and SNe data improves the constraints significantly. Our model predicts lower values of $f(z) \sigma_8(z)$ at $z<1$ comparing to $\Lambda$CDM model, which alleviates the tension of $\Lambda$CDM with various RSD data from optical galaxy surveys. Overall, the DM-DE interaction model is consistent with the current observational data, especially providing a better fit to the RSD data.

M87 is one of the best available source for studying the AGN jet-launching region. To enrich our knowledge of this region, with quasi-simultaneous observations using VLBA at 22, 43 and 86 GHz, we capture the images of the radio jet in M87 on a scale within several thousand R s . Based on the images, we analyze the transverse jet structure and obtain the most accurate spectral-index maps of the jet in M87 so far, then for the first time, we compare the results of the two analyses and find a spatial association between the jet collimations and the local enhancement of the density of external medium in the jet-launching region. We also find the external medium is not uniform, and greatly contributes to the free-free absorption in this region. In addition, we find for the jet in M87, its temporal morphology in the launching region may be largely affected by the local, short-lived kink instability growing in itself.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

Space very long baseline interferometry (VLBI) has unique applications in high-resolution imaging of fine structure of astronomical objects and high-precision astrometry due to the key long space-Earth or space-space baselines beyond the Earth's diameter. China has been actively involved in the development of space VLBI in recent years. This review briefly summarizes China's research progress in space VLBI and the future development plan.

The digital correlator is a crucial element in a modern radio telescope. In this paper we describe a scalable design of the correlator system for the Tianlai pathfinder array, which is an experiment dedicated to test the key technologies for conducting 21cm intensity mapping survey. The correlator is of the FX design, which firstly performs Fast Fourier Transform (FFT) including Polyphase Filter Bank (PFB) computation using a Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) Reconfigurable Open Architecture Computing Hardware-2 (ROACH2) board, then computes cross-correlations using Graphical Processing Units (GPUs). The design has been tested both in laboratory and in actual observation.

The neutral hydrogen 21cm line is potentially a very powerful probe of the observable universe, and a number of on-going experiments are trying to detect it at cosmological distances. However, the presence of strong foreground radiations such as the galactic synchrotron radiation, galactic free-free emission and extragalactic radio sources make it a very challenging task. For the design of 21cm experiments and analysis of their data, simulation is an essential tool, and good sky foreground model is needed. With existing data the whole sky maps are available only in low angular resolutions or for limited patches of sky, which is inadequate in the simulation of these new 21cm experiments. In this paper, we present the method of constructing a high resolution self-consistent sky model at low frequencies, which incorporates both diffuse foreground and point sources. Our diffuse map is constructed by generating physical foreground components including the galactic synchrotron emission and galactic free-free emission. The point source sample is generated using the actual data from the NRAO VLA Sky Survey (NVSS) and the Sydney University Molonglo Sky Survey (SUMSS) where they are available and complete in flux limit, and mock point sources according to statistical distributions. The entire model is made self-consistent by removing the integrated flux of the point sources from the diffuse map so that this part of radiation is not double counted. We show that with the point sources added, a significant angular power is introduced in the mock sky map, which may be important for foreground subtraction simulations. Our sky maps and point source catalogues are available to download.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

We propose an eigenvector-based formalism for the calibration of radio interferometer arrays. In the presence of a strong dominant point source, the complex gains of the array can be obtained by taking the first eigenvector of the visibility matrix. We use the stable principle component analysis (SPCA) method to help separate outliers and noise from the calibrator signal to improve the performance of the method. This method can be applied with poorly known beam model of the antenna, and is insensitive to outliers or imperfections in the data, and has low computational complexity. It thus is particularly suitable for the initial calibration of the array, which can serve as the initial point for more accurate calibrations. We demonstrate this method by applying it to the cylinder pathfinder of the Tianlai experiment, which aims to measure the dark energy equation of state using the baryon acoustic oscillation (BAO) features in the large scale structure by making intensity mapping observation of the redshifted 21~cm emission of the neutral hydrogen (HI). The complex gain of the array elements and the beam profile in the East-West direction (short axis of the cylinder) are successfully obtained by applying this method to the transit data of bright radio sources.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

Radio astronomical observation below 30 MHz is hampered by the refraction and absorption of the ionosphere, and the radio frequency interference (RFI), so far high angular resolution sky intensity map is not available. An interferometer array on lunar orbit provides a perfect observatory in this frequency band: it is out of the ionosphere and the Moon helps to block the RFIs from the Earth. The satellites can make observations on the far side of the Moon and then send back the data on the near side part of the orbit. However, for such array the traditional imaging algorithm is not applicable: the field of view is very wide (almost whole sky), and for baselines distributed on a plane, there is a mirror symmetry between the two sides of the plane. A further complication is that for each baseline, the Moon blocks part of the sky, but as the satellites orbit the Moon, both the direction of the baseline and the blocked sky change, so even imaging algorithms which can deal with non-coplanar baseline may not work in this case. Here we present an imaging algorithm based on solving the linear mapping equations relating the sky intensity to the visibilities. We show that the mirror symmetry can be broken by the three dimensional baseline distribution generated naturally by the precession of the orbital plane of the satellites. The algorithm is applicable and good maps could be reconstructed, even though for each baseline the sky blocking by the Moon is different. We also investigate how the map-making is affected by inhomogeneous baseline distributions.

Extracting the neutral hydrogen (HI) signal is a great challenge for cosmological 21cm experiments, both the astrophysical foregrounds and the receiver noise are typically several orders of magnitude greater than the 21cm signal. However, the different properties of the 21cm signal, foreground, and noise can be exploited to separate these components. The foregrounds are generally smooth or correlated over the frequency space along a line of sight (l.o.s.), while both the 21cm signal and the noise varies stochastically along the same l.o.s. The foreground can be removed by filtering out the smooth component in the frequency space. The receiver noise is basically uncorrelated for observations at different times, hence for surveys they are also uncorrelated in the different directions, while the 21cm signal which traces the large scale structure are correlated up to certain scales. In this exercise, we apply Wiener filters in frequency and angular space, to extract the 21cm signals. We found that the method works well. The inaccurate knowledge about beam could degrade the reconstruction, but the overall result is still good, showing that the method is fairly robust.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

A wide bandwidth, dual polarized, modified four-square antenna is presented as a feed antenna for radio astronomical measurements. A linear array of these antennas is used as a line-feed for cylindrical reflectors for Tianlai, a radio interferometer designed for 21~cm intensity mapping. Simulations of the feed antenna beam patterns and scattering parameters are compared to experimental results at multiple frequencies across the 650 - 1420 MHz range. Simulations of the beam patterns of the combined feed array/reflector are presented as well.

In this work the relativistic mean field (RMF) FSUGold model extended to include hyperons is employed to study the properties of neutron stars with strong magnetic fields. The chaotic magnetic field approximation is utilized. The effect of the anomalous magnetic moments (AMMs) is also investigated. It is shown that the equation of state (EOS)of neutron star matter is stiffened by the presence of the magnetic field, which increases the maximum mass of neutron star by around 6%. The AMMs only have a small influence on the EOS of neutron star matter, and increase the maximum mass of neutron star by 0.02 M$_{sun}$. Neutral particles are spin polarized due to the presence of the AMMs.

We report on four radio-detected cosmic-ray (CR) or CR-like events observed with the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload. Two of the four were previously identified as stratospheric CR air showers during the ANITA-I flight. A third stratospheric CR was detected during the ANITA-II flight. Here we report on characteristics these three unusual CR events, which develop nearly horizontally, 20-30~km above the surface of the Earth. In addition, we report on a fourth steeply upward-pointing ANITA-I CR-like radio event which has characteristics consistent with a primary that emerged from the surface of the ice. This suggests a possible $\tau$-lepton decay as the origin of this event, but such an interpretation would require significant suppression of the Standard Model ${\tau}$-neutrino cross section.

In this paper, we apply our sky map reconstruction method for transit type interferometers to the Tianlai cylinder array. The method is based on the spherical harmonic decomposition, and can be applied to cylindrical array as well as dish arrays and we can compute the instrument response, synthesised beam, transfer function and the noise power spectrum. We consider cylinder arrays with feed spacing larger than half wavelength, and as expected, we find that the arrays with regular spacing have grating lobes which produce spurious images in the reconstructed maps. We show that this problem can be overcome, using arrays with different feed spacing on each cylinder. We present the reconstructed maps, and study the performance in terms of noise power spectrum, transfer function and beams for both regular and irregular feed spacing configurations.

The spherical harmonics $m$-mode decomposition is a powerful sky map reconstruction method suitable for radio interferometers operating in transit mode. It can be applied to various configurations, including dish arrays and cylinders. We describe the computation of the instrument response function, the point spread function (PSF), transfer function, the noise covariance matrix and noise power spectrum. The analysis in this paper is focused on dish arrays operating in transit mode. We show that arrays with regular spacing have more pronounced side lobes as well as structures in their noise power spectrum, compared to arrays with irregular spacing, specially in the north-south direction. A good knowledge of the noise power spectrum $C^{\mathrm{noise}}(\ell)$ is essential for intensity mapping experiments as non uniform $C^{\mathrm{noise}}(\ell)$ is a potential problem for the measurement of the HI power spectrum. Different configurations have been studied to optimise the PAON-4 and Tianlai dish array layouts. We present their expected performance and their sensitivities to the 21-cm emission of the Milky Way and local extragalactic HI clumps

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. For 14 of these events, this radiation was reflected from the ice. The dominant contribution to the radiation from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. This radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of 36km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 EeV. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations. In addition, we find that the Monte Carlo simulation of the ANITA data set is in agreement with the total number of observed events and with the properties of those events.

We combine new Cosmic Microwave Background (CMB) data from Planck with Baryon Acoustic Oscillation (BAO) data to constrain the Brans-Dicke (BD) theory, in which the gravitational constant $G$ evolves with time. Observations of type Ia supernovae (SNeIa) provide another important set of cosmological data, as they may be regarded as standard candles after some empirical corrections. However, in theories that include modified gravity like the BD theory, there is some risk and complication when using the SNIa data because their luminosity may depend on $G$. In this paper, we assume a power law relation between the SNIa luminosity and $G$, but treat the power index as a free parameter. We then test whether the difference in distances measured with SNIa data and BAO data can be reduced in such a model. We also constrain the BD theory and cosmological parameters by making a global fit with the CMB, BAO and SNIa data set. For the CMB+BAO+SNIa data set, we find $0.08\times10^{-2} < \zeta <0.33\times10^{-2} $ at the 68\% confidence level (CL) and $-0.01\times10^{-2} <\zeta <0.43\times 10^{-2} $ at the 95\% CL, where $\zeta$ is related to the {BD} parameter $\omega$ by $\zeta=\ln(1+1/\omega)$.

Morphology of the complex HI gas distribution can be quantified by statistics like the Minkowski functionals, and can provide a way to statistically study the large scale structure in the HI maps both at low redshifts, and during the epoch of reionization (EoR). At low redshifts, the 21cm emission traces the underlying matter distribution. Topology of the HI gas distribution, as measured by the genus, could be used as a "standard ruler". This enables the determination of distance-redshift relation and also the discrimination of various models of dark energy and of modified gravity. The topological analysis is also sensitive to certain primordial non-Gaussian features. Compared with two-point statistics, the topological statistics are more robust against the nonlinear gravitational evolution, bias, and redshift-space distortion. The HI intensity map observation naturally avoids the sparse sampling distortion, which is an important systematic in optical galaxy survey. The large cosmic volume accessible to SKA would provide unprecedented accuracy using such a measurement... [abridged]

A large inflationary tensor-to-scalar ratio $r_\mathrm{0.002} = 0.20^{+0.07}_{-0.05}$ is reported by the BICEP2 team based on their B-mode polarization detection, which is outside of the $95\%$ confidence level of the Planck best fit model. We explore several possible ways to reduce the tension between the two by considering a model in which $\alpha_\mathrm{s}$, $n_\mathrm{t}$, $n_\mathrm{s}$ and the neutrino parameters $N_\mathrm{eff}$ and $\Sigma m_\mathrm{\nu}$ are set as free parameters. Using the Markov Chain Monte Carlo (MCMC) technique to survey the complete parameter space with and without the BICEP2 data, we find that the resulting constraints on $r_\mathrm{0.002}$ are consistent with each other and the apparent tension seems to be relaxed. Further detailed investigations on those fittings suggest that $N_\mathrm{eff}$ probably plays the most important role in reducing the tension. We also find that the results obtained from fitting without adopting the consistency relation do not deviate much from the consistency relation. With available Planck, WMAP, BICEP2 and BAO datasets all together, we obtain $r_{0.002} = 0.14_{-0.11}^{+0.05}$, $n_\mathrm{t} = 0.35_{-0.47}^{+0.28}$, $n_\mathrm{s}=0.98_{-0.02}^{+0.02}$, and $\alpha_\mathrm{s}=-0.0086_{-0.0189}^{+0.0148}$; if the consistency relation is adopted, we get $r_{0.002} = 0.22_{-0.06}^{+0.05}$.

Inspired by holographic principle, we suggest that the density of dark energy is proportional to the spatial Ricci scalar curvature (SRDE). Such model is phenomenologically viable. The best fit values of its parameters at 68% confidence level are found to be: $\Omega_{\rm m0}=0.259\pm0.016$ and $\alpha=0.261\pm0.0122$, constrained from the Union+CFA3 sample of 397 SNIa and the BAO measurement. We find the equation of state of SRDE crosses -1 at $z\simeq-0.14$. The present values of the deceleration parameter $q(z)$ for SRDE is found to be $q_{z=0}\sim -0.85$. The phase transition from deceleration to acceleration of the Universe for SRDE occurs at the redshift $z_{q=0}\sim 0.4$. After studying on the perturbation of each component of the Universe, we show that the matter power spectra and cosmic microwave background temperature anisotropy is slightly affected by SRDE, compared with $\Lambda$CDM.

Combining the latest Planck, Wilkinson Microwave Anisotropy Probe (WMAP), and baryon acoustic oscillation (BAO) data, we exploit the recent cosmic microwave background (CMB) B-mode power spectra data released by the BICEP2 collaboration to constrain the cosmological parameters of the $\Lambda$CDM model, esp. the primordial power spectra parameters of the scalar and the tensor modes, $n_s,\alpha_s, r, n_t$. We obtained constraints on the parameters for a lensed $\Lambda$CDM model using the Markov Chain Monte Carlo (MCMC) technique, the marginalized $68\%$ bounds are $r = 0.1043\pm_{0.0914}^{0.0307}$, $n_s = 0.9617\pm_{0.0061}^{ 0.0061}$, $\alpha_s = -0.0175\pm_{0.0097}^{0.0105}$, $n_t = 0.5198\pm_{0.4579}^{0.4515}$. We found that a blue tilt for $n_t$ is favored slightly, but it is still well consistent with flat or even red tilt. Our $r$ value is slightly smaller than the one obtained by the BICEP team, as we allow $n_t$ as a free parameter without imposing the single-field slow roll inflation consistency relation. If we impose this relation, $r=0.2130\pm_{0.0609}^{0.0446}$. For most other parameters, the best fit values and measurement errors are not altered much by the introduction of the BICEP2 data.

LAMOST has released more than two million spectra, which provide the opportunity to search for double-peaked narrow emission line (NEL) galaxies and AGNs. The double-peaked narrow-line profiles can be well modeled by two velocity components, respectively blueshifted and redshifted with respect to the systemic recession velocity. This paper presents 20 double-peaked NEL galaxies and AGNs found from LAMOST DR1 using a search method based on multi-gaussian fit of the narrow emission lines. Among them, 10 have already been published by other authors, either listed as genuine double-peaked NEL objects or as asymmetric NEL objects, the remaining 10 being first discoveries. We discuss some possible origins for double-peaked narrow-line features, as interaction between jet and narrow line regions, interaction with companion galaxies and black hole binaries. Spatially resolved optical imaging and/or follow-up observations in other spectral bands are needed to further discuss the physical mechanisms at work.

Carbon stars and DZ white dwarfs are two types of rare objects in the Galaxy. In this paper, we have applied the label propagation algorithm to search for these two types of stars from Data Release Eight (DR8) of the Sloan Digital Sky Survey (SDSS), which is verified to be efficient by calculating precision and recall. From nearly two million spectra including stars, galaxies and QSOs, we have found 260 new carbon stars in which 96 stars have been identified as dwarfs and 7 identified as giants, and 11 composition spectrum systems (each of them consists of a white dwarf and a carbon star). Similarly, using the label propagation method, we have obtained 29 new DZ white dwarfs from SDSS DR8. Compared with PCA reconstructed spectra, the 29 findings are typical DZ white dwarfs. We have also investigated their proper motions by comparing them with proper motion distribution of 9,374 white dwarfs, and found that they satisfy the current observed white dwarfs by SDSS generally have large proper motions. In addition, we have estimated their effective temperatures by fitting the polynomial relationship between effective temperature and g-r color of known DZ white dwarfs, and found 12 of the 29 new DZ white dwarfs are cool, in which nine are between 6000K and 6600K, and three are below 6000K.

Based on the new cosmic CMB temperature data from the Planck satellite, the 9 year polarization data from the WMAP, the BAO distance ratio data from the SDSS and 6dF surveys, we place a new constraint on the Brans-Dicke theory. We adopt a parametrization $\zeta=\ln(1+1/\omega})$, where the general relativity (GR) limit corresponds to $\zeta = 0$. We find no evidence of deviation from general relativity. At 95% probability, $-0.00246 < \zeta < 0.00567$, correspondingly, the region $-407.0 < \omega <175.87$ is excluded. If we restrict ourselves to the $\zeta>0$ (i.e. $\omega >0$) case, then the 95% probability interval is $\zeta<0.00549$, corresponding to $\omega> 181.65$. We can also translate this result to a constraint on the variation of gravitational constant, and find the variation rate today as $\dot{G}=-1.42^{+2.48}_{-2.27} \times 10^{-13} $yr$^{-1} $ ($1\sigma$ error bar), the integrated change since the epoch of recombination is $\delta G/G = 0.0104^{+0.0186}_{-0.0067} $ ($1\sigma$ error bar). These limits on the variation of gravitational constant are comparable with the precision of solar system experiments.

Aims. A long timeline kinematic study of the archetypal CSO OQ 208 sheds light on the physical properties of the most compact radio sources. Methods. Archival data from the VLBA at 15 GHz over a time span of 13.6 yr are used to investigate the kinematics of the radio source. The flux density monitoring data obtained at the Michigan 26-meter radio telescope are also used as supplementary information. Results. At 8.4 and 15 GHz, the two lobes are resolved into two sub-components, identified as hotspots. A knotty jet is linked with the NE hotspot and traces back toward the geometric center. The core is too weak to be detected. Significant flux density variation is found in the primary hotspots with the maximum level of 62% (NE1) and 19% (SW1). The peak in the flux density of NE1 leads that of SW1 by approximately 5.00 yr, suggesting that the northeast lobe is advancing and the southwest lobe is receding. This light travel difference indicates a radial distance difference between the two hotspots of 1.53 pc, which indicates an inclination angle of about 80.8 degree between the radio jet and the line of sight. The angular separation rate between NE1 and SW1 is 0.027 mas/yr (or 0.133 c). The inner jet knot moves at 0.047 mas/yr (or 0.230 c), about 3.5 times the hotspot advancing speed. Conclusions. The large viewing angle and the modest jet speed suggest a mildly relativistic jet. The jet axis is close to the plane of the sky. The separation rate and the distance between the two primary hotspots result in a kinematic age of 255$\pm$17 yr, confirming that OQ 208 is indeed a young radio source. In addition to the hotspot advancing motions, sideways motions provide evidence that the lobes are obstructed by the external interstellar medium.

Observations of ten Compact Symmetric Objects ({\rm CSO}) candidates have been made with the Very Long Baseline Array at 8.4 GHz in 2005 and with a combined Chinese and European VLBI array at 8.4 GHz in 2009. The 2009 observations incorporate for the first time the two new Chinese telescopes at Miyun and Kunming for international astrophysical observations. The observational data, in combination with archival VLBA data from previous epochs, have been used to derive the proper motions of the VLBI components. Because of the long time baseline of $\sim$16 years of the VLBI data sets, the expansion velocities of the hot spots can be measured at an accuracy as high as $\sim$1.3 $\mu$as yr$^{-1}$. Six of the ten sources are identified as CSOs with a typical double or triple morphology on the basis of both spectral index maps and their mirror-symmetry of proper motions of the terminal hot spots. The compact double source J1324+4048 is also identified as a CSO candidate. Among the three remaining sources, J1756+5748 and J2312+3847 are identified as core-jet sources with proper motions of their jet components relating to systemic source expansion. The third source J0017+5312 is likely also a core-jet source, but a robust detection of a core is needed for an unambiguous identification. The kinematic ages of the CSOs derived from proper motions range from 300 to 2500 years. The kinematic age distribution of the CSOs confirm an overabundance of compact young CSOs with ages less than 500 years. CSOs with known kinematic ages may be used to study the dynamical evolution of extragalactic radio sources at early stages.

The blue-shifted broad absorption lines (BAL) or troughs are observed in Active Galactic Nuclei (AGNs) when our line of sight is intercepted by a high speed outflow (wind), likely originating in the accretion disc. The outflow or wind can shed light on the internal structure obscured by the AGN torus. Recently, it has been shown that this outflow is rotating in the BAL quasar PG 1700+518, further supporting the accretion disc origin of the wind. With the purpose of giving independent constraints on the wind geometry, we performed high-resolution European VLBI Network (EVN) observations at 1.6 GHz in 2010. Combining the VLBI (Very Long Baseline Interferometry) results with the Very Large Array (VLA) archival data at 8.4 GHz, we present its jet structure on scales of parsec (pc) to kiloparsec (kpc) for the first time. The source shows two distinct jet features in East-West direction with a separation of around 4 kpc. The Eastern feature, which has so far been assumed to hide the core, is a kpc-scale hot spot, which is completely resolved out in the EVN image. In the western jet feature, we find a compact jet component, which pinpoints the position of the central black hole in the galaxy. Jet components on both sides of the core are additionally detected in the Northwest-Southeast direction, and they show a symmetric morphology on scale of <1 kpc. This two-sided jet feature is not common in the known BAL quasars and indicates that the jet axis is far away from the line of sight. Furthermore, it is nearly parallel to the scattering plane revealed earlier by optical polarimetry. By analogy to polar-scattered Seyfert 1 galaxies, we conclude that the jet likely has a viewing angle around 45 degree. The analogy is further supported by the recent report of significant cold absorption in the soft X-rays, a nearly unique feature to polar-scattered Seyfert galaxies.

Using the covariant formalism developed in a companion paper (paper I), we derive observational constraint on the Brans-Dicke model in a flat FLRW universe with cosmological constant and cold dark matter. The CMB observations we use include the WMAP five year data, the ACBAR 2007 data, the CBI polarization data, and the BOOMERanG 2003 flight data. For the large scale structure we use the matter power spectrum data measured with the LRG survey of the SDSS DR4. We parameterize the Brans-Dicke parameter $\omega$ with a new parameter $\zeta=\ln(1/\omega+1)$, and use the Markov-Chain Monte Carlo method to explore the parameter space. We find that using CMB data alone, one could place some constraint on positive $\zeta$ or $\omega$, but negative $\zeta$ or $\omega$ could not be constrained effectively. However, with additional large scale structure data, one could break the degeneracy at $\zeta<0$. The $2\sigma$ (95.5%) bound on $\zeta$ is $-0.00837<\zeta<0.01018$ (corresponding to $\omega < -120.0$ or $\omega > 97.8$). We also obtained constraint on $\dot{G}/G$, the rate of change of $G$ at present, as $-1.75 \times 10^{-12} \yr^{-1}<\dot{G}/G < 1.05 \times 10^{-12}\yr^{-1}$, and $\delta G/G$, the total variation of $G$ since the epoch of recombination, as $ -0.083 < \delta{G}/G < 0.095$ at $2\sigma$ confidence level.

In the covariant cosmological perturbation theory, a 1+3 decomposition ensures that all variables in the frame-independent equations are covariant, gauge-invariant and have clear physical interpretations. We develop this formalism in the case of Brans-Dicke gravity, and apply this method to the calculation of cosmic microwave background (CMB) anisotropy and large scale structures (LSS). We modify the publicly available Boltzmann code CAMB to calculate numerically the evolution of the background and adiabatic perturbations, and obtain the temperature and polarization spectra of the Brans-Dicke theory for both scalar and tensor mode, the tensor mode result for the Brans-Dicke gravity are obtained numerically for the first time. We first present our theoretical formalism in detail, then explicitly describe the techniques used in modifying the CAMB code. These techniques are also very useful to other gravity models. Next we compare the CMB and LSS spectra in Brans-Dicke theory with those in the standard general relativity theory. At last, we investigate the ISW effect and the CMB lensing effect in the Brans-Dicke theory. Constraints on Brans-Dicke model with current observational data is presented in a companion paper (paper II).

We report the observation of sixteen cosmic ray events of mean energy of 1.5 x 10^{19} eV, via radio pulses originating from the interaction of the cosmic ray air shower with the Antarctic geomagnetic field, a process known as geosynchrotron emission. We present the first ultra-wideband, far-field measurements of the radio spectral density of geosynchrotron emission in the range from 300-1000 MHz. The emission is 100% linearly polarized in the plane perpendicular to the projected geomagnetic field. Fourteen of our observed events are seen to have a phase-inversion due to reflection of the radio beam off the ice surface, and two additional events are seen directly from above the horizon.

The letter presents 25 discovered supernova candidates from SDSS-DR7 with our dedicated method, called Sample Decrease, and 10 of them were confirmed by other research groups, and listed in this letter. Another 15 are first discovered including 14 type Ia and one type II based on Supernova Identification (SNID) analysis. The results proved that our method is reliable, and the description of the method and some detailed spectra analysis procedures were also presented in this letter.

The Large Area Multi-Object Spectroscopic Telescope (LAMOST) is a dedicated spectroscopic survey telescope being built in China, with an effective aperture of 4 meters and equiped with 4000 fibers. Using the LAMOST telescope, one could make redshift survey of the large scale structure (LSS). The baryon acoustic oscillation (BAO) features in the LSS power spectrum provide standard rulers for measuring dark energy and other cosmological parameters. In this paper we investigate the meaurement precision achievable for a few possible surveys: (1) a magnitude limited survey of all galaxies, (2) a survey of color selected red luminous galaxies (LRG), and (3) a magnitude limited, high density survey of z<2 quasars. For each survey, we use the halo model to estimate the bias of the sample, and calculate the effective volume. We then use the Fisher matrix method to forecast the error on the dark energy equation of state and other cosmological parameters for different survey parameters. In a few cases we also use the Markov Chain Monte Carlo (MCMC) method to make the same forecast as a comparison. The fiber time required for each of these surveys is also estimated. These results would be useful in designing the surveys for LAMOST.

Motivated by the holographic principle, it has been suggested that the dark energy density may be inversely proportional to the area of the event horizon of the Universe. However, such a model would have a causality problem. In this paper, we propose to replace the future event horizon area with the inverse of the Ricci scalar curvature. We show that this model does not only avoid the causality problem and is phenomenologically viable, but also naturally solves the coincidence problem of dark energy. Our analysis of the evolution of density perturbations show that the matter power spectra and CMB temperature anisotropy is only slightly affected by such modification.

We report initial results of the first flight of the Antarctic Impulsive Transient Antenna (ANITA-1) 2006-2007 Long Duration Balloon flight, which searched for evidence of a diffuse flux of cosmic neutrinos above energies of 3 EeV. ANITA-1 flew for 35 days looking for radio impulses due to the Askaryan effect in neutrino-induced electromagnetic showers within the Antarctic ice sheets. We report here on our initial analysis, which was performed as a blind search of the data. No neutrino candidates are seen, with no detected physics background. We set model-independent limits based on this result. Upper limits derived from our analysis rule out the highest cosmogenic neutrino models. In a background horizontal-polarization channel, we also detect six events consistent with radio impulses from ultra-high energy extensive air showers.

We present a detailed report on the experimental details of the Antarctic Impulsive Transient Antenna (ANITA) long duration balloon payload, including the design philosophy and realization, physics simulations, performance of the instrument during its first Antarctic flight completed in January of 2007, and expectations for the limiting neutrino detection sensitivity. Neutrino physics results will be reported separately.

Radiowave detection of the Cherenkov radiation produced by neutrino-ice collisions requires an understanding of the radiofrequency (RF) response of cold polar ice. We herein report on a series of radioglaciological measurements performed approximately 10 km north of Taylor Dome Station, Antarctica from Dec. 6, 2006 - Dec. 16, 2006. Using RF signals broadcast from: a) an englacial discone, submerged to a depth of 100 meters and broadcasting to a surface dual polarization horn receiver, and b) a dual-polarization horn antenna on the surface transmitting signals which reflect off the underlying bed and back up to the surface receiver, we have made time-domain estimates of both the real (index-of-refraction) and imaginary (attenuation length) components of the complex ice dielectric constant. We have also measured the uniformity of ice response along two orthogonal axes in the horizontal plane. We observe a wavespeed asymmetry of order 0.1%, projected onto the vertical propagation axis, consistent with some previous measurements, but somewhat lower than others.

The holographic dark energy model is proposed by Li as an attempt for probing the nature of dark energy within the framework of quantum gravity. The main characteristic of holographic dark energy is governed by a numerical parameter $c$ in the model. The parameter $c$ can only be determined by observations. Thus, in order to characterize the evolving feature of dark energy and to predict the fate of the universe, it is of extraordinary importance to constrain the parameter $c$ by using the currently available observational data. In this paper, we derive constraints on the holographic dark energy model from the latest observational data including the gold sample of 182 Type Ia supernovae (SNIa), the shift parameter of the cosmic microwave background (CMB) given by the three-year {\it Wilkinson Microwave Anisotropy Probe} ({\it WMAP}) observations, and the baryon acoustic oscillation (BAO) measurement from the Sloan Digital Sky Survey (SDSS). The joint analysis gives the fit results in 1-$\sigma$: $c=0.91^{+0.26}_{-0.18}$ and $\Omega_{\rm m0}=0.29\pm 0.03$. That is to say, though the possibility of $c<1$ is more favored, the possibility of $c>1$ can not be excluded in one-sigma error range, which is somewhat different from the result derived from previous investigations using earlier data. So, according to the new data, the evidence for the quintom feature in the holographic dark energy model is not as strong as before.

A new model for the description of the behaviour of strangelets in the Earth's atmosphere is presented. Strangelet fission induced by collision with air nuclei is included. It is shown that strangelets with certain parameters of initial mass and energy may reach depths near sea level, which can be examined by ground-based experiments.

We report on the first observations of the Askaryan effect in ice: coherent impulsive radio Cherenkov radiation from the charge asymmetry in an electromagnetic (EM) shower. Such radiation has been observed in silica sand and rock salt, but this is the first direct observation from an EM shower in ice. These measurements are important since the majority of experiments to date that rely on the effect for ultra-high energy neutrino detection are being performed using ice as the target medium. As part of the complete validation process for the Antarctic Impulsive Transient Antenna (ANITA) experiment, we performed an experiment at the Stanford Linear Accelerator Center (SLAC) in June 2006 using a 7.5 metric ton ice target, yielding results fully consistent with theoretical expectations.

Noncommutative inflation is based upon the consideration of some effects of the space-time uncertainty principle motivated by ideas from string/M theory. The CMB anisotropies may carry a signature of this very early Universe correction from the space-time uncertainty principle and can be used to place constraints on the parameters of the noncommutative inflation model. In this paper we analyze the noncommutative chaotic inflation model by means of the WMAP three year results. We show that the noncommutative chaotic inflation model can produce a large negative running of spectral index within a reasonable range of $e$-folding number $N$, provided that the value of $p$ in the potential $V(\phi)\propto \phi^p$ is enhanced roughly to $p\sim 12-18$. Applying a likelihood analysis to the WMAP results, we find the best-fit values of parameters for the noncommutative chaotic inflation model. In addition, this model predicts a rather big gravitational wave which may be tested by the upcoming experiments.

We report new limits on cosmic neutrino fluxes from the test flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment, which completed an 18.4 day flight of a prototype long-duration balloon payload, called ANITA-lite, in early 2004. We search for impulsive events that could be associated with ultra-high energy neutrino interactions in the ice, and derive limits that constrain several models for ultra-high energy neutrino fluxes. We rule out the long-standing Z-burst model as the source for the ultra-high energy cosmic rays.

We propose in this paper a new model for describing the unification of dark energy and dark matter. This new model is a further generalization of the generalized Chaplygin gas (GCG) model, thus dubbed new generalized Chaplygin gas (NGCG) model. The equation of state of the NGCG is given by $p = - \tilde{A}(a)/\rho^{\alpha}$, where $a$ is the scale factor and $\tilde{A}(a)=-w_XAa^{-3(1+w_X)(1+\alpha)}$. We show that the NGCG model is totally dual to an interacting XCDM parametrization scenario, in which the interaction between dark energy and dark matter is characterized by the constant $\alpha$. We discuss the cosmological consequences led by such an unified dark sectors model. Furthermore, we perform a statefinder analysis on this scenario and show the discrimination between this scenario and other dark energy models. Finally, a combined analysis of the data of Type Ia supernovae, cosmic microwave background, and large scale structure provides a fairly tight constraint on the parameters of the NGCG model.

We use the Chandra measurements of the X-ray gas mass fraction of 26 rich clusters released by Allen et al. to perform constraints on the holographic dark energy model. The constraints are consistent with those from other cosmological tests, especially with the results of a joint analysis of supernovae, cosmic microwave background, and large scale structure data. From this test, the holographic dark energy also tends to behave as a quintom-type dark energy.

In this paper, we use the type Ia supernovae data to constrain the holographic dark energy model proposed by Li. We also apply a cosmic age test to this analysis. We consider in this paper a spatially flat Friedmann-Robertson-Walker Universe with matter component and holographic dark energy component. The fit result shows that the case $c<1$ ($c=0.21$) is favored, which implies that the holographic dark energy behaves as a quintom-type dark energy. Furthermore, we also perform a joint analysis of SNe+CMB+LSS to this model; the result is well improved, and still upholds the quintom dark energy conclusion. The best fit results in our analysis are $c=0.81$, $\Omega_m^0=0.28$, and $h=0.65$, which lead to the present equation of state of dark energy $w_0=-1.03$ and the deceleration/acceleration transition redshift $z_T=0.63$. Finally, an expected SNAP simulation using $\Lambda$CDM as a fiducial model is performed on this model, and the result shows that the holographic dark energy model takes on $c<1$ ($c=0.92$) even though the dark energy is indeed a cosmological constant.

We describe a new experiment to search for neutrinos with energies above 3 x 10^18 eV based on the observation of short duration radio pulses that are emitted from neutrino-initiated cascades. The primary objective of the ANtarctic Impulse Transient Antenna (ANITA) mission is to measure the flux of Greisen-Zatsepin-Kuzmin (GZK) neutrinos and search for neutrinos from Active Galactic Nuclei (AGN). We present first results obtained from the successful launch of a 2-antenna prototype instrument (called ANITA-lite) that circled Antarctica for 18 days during the 03/04 Antarctic campaign and show preliminary results from attenuation length studies of electromagnetic waves at radio frequencies in Antarctic ice. The ANITA detector is funded by NASA, and the first flight is scheduled for December 2006.

Using the method proposed in a previous paper, we calculate pulsar braking indices in the models with torque contributions from both inner and outer accelerating regions, assuming that the interaction between them is negligible. We suggest that it is likely that the inverse Compton scattering induced polar vacuum gap and the outer gap coexist in the pulsar magnetosphere. We include the new near threshold vacuum gap models with curvature-radiation and inverse Compton scattering induced cascades, respectively; and find that these models can well reproduce the measured values of the braking indices.

The conjecture that ultra high energy cosmic rays (UHECRs) are actually strangelets is discussed. Besides the reason that strangelets can do as cosmic rays beyond the GZK-cutoff, another argument to support the conjecture is addressed in this letter via the study of formation of TeV-scale microscopic black holes when UHECRs bombarding bare strange stars. It is proposed that the exotic quark surface of a bare strange star could be an effective astro-laboratory in the investigations of the extra dimensions and of the detection of ultra-high energy neutrino fluxes. The flux of neutrinos (and other point-like particles) with energy >2.3 x 10^{20} eV could be expected to be smaller than 10^{-26} cm^{-2}$ s^{-1} if there are two extra spatial dimensions.