Gravitational wave (GW) data can be used to test the parity symmetry of gravity by investigating the difference between left-hand and right-hand circular polarization modes. In this article, we develop a method to decompose the circular polarizations of GWs produced during the inspiralling stage of compact binaries, with the help of stationary phase approximation. The foremost advantage is that this method is simple, clean, independent of GW waveform, and is applicable to the existing detector network. Applying it to the mock data, we test the parity symmetry of gravity by constraining the velocity birefringence of GWs. If a nearly edge-on binary neutron-stars with observed electromagnetic counterparts at 40 Mpc is detected by the second-generation detector network, one could derive the model-independent test on the parity symmetry in gravity: the lower limit of the energy scale of parity violation can be constrained within $\mathcal{O}(10^4{\rm eV})$.

In this paper, we systematically study spherically symmetric static spacetimes in the framework of Einstein-aether theory, and pay particular attention to the existence of black holes (BHs). In the theory, one scalar and one vector two gravitational modes additionally appear. To avoid the vacuum gravi-\v{C}erenkov radiation, they must all propagate with speeds greater than or at least equal to the speed of light. In the spherical case, only the scalar mode is relevant, so BH horizons are defined by this mode. In the present studies we first clarify several subtle issues. In particular, we find that, out of the five non-trivial field equations, only three are independent, so the problem is well-posed, as now there are only three unknown functions, two, $F(r)$ and $B(r)$, describe the spacetime and one, $A(r)$, describes the aether field. In addition, the two second-order differential equations for $A$ and $F$ are independent of $B$, and once they are found, $B$ is given simply by an algebraic expression of $F,\; A$ and their derivatives. To simplify the problem further, we explore the symmetry of field redefinitions, and work first with the redefined metric and aether field, and then obtain the physical ones by the inverse transformations. These clarifications significantly simplify the computational labor, which is important, as the problem is highly involved mathematically. In fact, it is exactly because of these, we find various numerical BH solutions with an accuracy that is at least two orders higher than previous ones. More important, these BH solutions are the only ones that satisfy the self-consistent conditions and meantime are consistent with all the observational constraints obtained so far. The locations of universal horizons are also identified, together with several other observationally interesting quantities.

The discovery of gravitational-wave (GW) signals, produced by the coalescence of stellar-mass binary black holes (SBBHs), opens a new window to study the astrophysical origins and dynamical evolutions of compact binaries. In addition, these GW events can be treated as the standard sirens to constrain various cosmological parameters. Both issues require the host identification for these GW events, with help of the spatial resolution of GW detector networks. In this paper, we investigate the capabilities of various detector networks for identifying the SBBHs' host galaxy groups, rather than their host galaxies, which can overcome the influence of galaxies' proper motions in dark matter halos for measuring the cosmological parameters. In our analysis, the group catalog of SDSS DR7 with redshift $z\in(0.01,0.1)$ is considered as an example of the application. We find that for the second-generation (2G) detector network, the host galaxy groups of around $17\%$ $30-30\ M_{\odot}$ SBBHs can be identified. For the 3G detector network, this fraction becomes $\sim 99\%$. We also investigate the potential constraint on the Hubble constant $H_0$ by these GW events, if their redshift information is extracted from the candidates of host galaxy groups. We find that, by five-year observations, 2G detector network is expected to give a constraint of $\Delta H_0/H_0\sim (1\%,4\%)$, which can be more than two order smaller if considering the 3G detector network.

We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

Testing gravity by binary pulsars nowadays becomes a key issue. Screened modified gravity is a kind of scalar-tensor theory with screening mechanism in order to satisfy the tight Solar System tests. In this paper, we investigate how the screening mechanism affects the orbital dynamics of binary pulsars, and calculate in detail the five post-Keplerian (PK) parameters in this theory. These parameters differ from those of general relativity (GR), and the differences are quantified by the scalar charges, which lead to the dipole radiation in this theory. We combine the observables of PK parameters for the ten binary pulsars, respectively, to place the constraints on the scalar charges and possible deviations from GR. The dipole radiation in the neutron star (NS) - white dwarf (WD) binaries leads to more stringent constraints on deviations from GR. The most constraining systems for the scalar charges of NS and WD are PSR~B1913$+$16 and PSR~J1738$+$0333, respectively. The results of all tests exclude significant strong-field deviations and show good agreement with GR.

Gravitational waves generated by the coalescence of compact binary open a new window to test the fundamental properties of gravity in the strong-field and dynamical regime. In this work, we focus on the parity symmetry of gravity which, if broken, can leave imprints on the waveform of gravitational wave. We construct generalized waveforms with amplitude and velocity birefringence due to parity violation in the effect field theory formalism, then analyze the open data of the ten binary black-hole merger events and the two binary neutron-star merger events detected by LIGO and Virgo collaboration. We do not find any signatures of violation of gravitational parity conservation, thereby setting the lower bound of the parity-violating energy scale to be $0.07$ GeV. This presents the first observational evidence of the parity conservation of gravity at high energy scale, about 17 orders of magnitude tighter than the constraints from the Solar system tests and binary pulsar observation. The third-generation gravitational-wave detector is capable of probing the parity-violating energy scale at $\mathcal{O}(10^2)$ GeV.

The tests of parity symmetry in the gravitational interaction is an attractive issue in the gravitational-wave astronomy. In the general theories of gravity with parity violation, one of the fundamental results is that the primordial gravitational waves (PGWs) produced during the slow-roll inflation is circularly polarized. In this article, we investigate the polarization of PGWs in the recently proposed ghost-free parity-violating gravity, which generalizes the Chern-Simons gravity by including the higher derivatives of the coupling scalar field. For this purpose, we first construct the approximate analytical solution to the mode function of the PGWs during the slow-roll inflation by using the uniform asymptotic approximation. With the approximate solution, we calculate explicitly the power spectrum and the corresponding circular polarization of the PGWs analytically, and find that the contributions of the higher derivatives of the coupling scalar field to the circular polarization is of the same order of magnitude as that of Chern-Simons gravity. The degree of circular polarization of PGWs is suppressed by the energy scale of parity violation in gravity, which is hardly detected on the basis solely of two-point statistics from future cosmic microwave background data.

As part of the survey component of the Megamaser Cosmology Project, we have discovered a disk megamaser system in the galaxy CGCG 074-064. Using the GBT and the VLA, we have obtained spectral monitoring observations of this maser system at a monthly cadence over the course of two years. We find that the systemic maser features display line-of-sight accelerations of ~4.4 km s$^{-1}$ yr$^{-1}$ that are nearly constant with velocity, while the high-velocity maser features show accelerations that are consistent with zero. We have also used the HSA to make a high-sensitivity VLBI map of the maser system in CGCG 074-064, which reveals that the masers reside in a thin, edge-on disk with a diameter of ~1.5 mas (0.6 pc). Fitting a three-dimensional warped disk model to the data, we measure a black hole mass of $2.42^{+0.22}_{-0.20} \times 10^7$ M$_{\odot}$ and a geometric distance to the system of $87.6^{+7.9}_{-7.2}$ Mpc. Assuming a CMB-frame recession velocity of $7308 \pm 150$ km s$^{-1}$, we constrain the Hubble constant to $H_0 = 81.0^{+7.4}_{-6.9}$ (stat.) $\pm 1.4$ (sys.) km s$^{-1}$ Mpc$^{-1}$.

As an extension of our previous work [J.Qiao, T.Zhu, W.Zhao & A.Wang, arXiv:1909.03815], in this article, we calculate the effects of parity violation on gravitational-wave (GW) waveforms during their propagation in the most general parity-violating gravities, including Chern-Simons modified gravity, ghost-free scalar-tensor gravity, symmetric teleparallel equivalence of GR theory, Ho\v{r}ava-Lifshitz gravity and so on. For this purpose, we consider the GWs generated by the coalescence of compact binaries and concentrate on the imprints of the parity violation in the propagation of GWs. With a unified description of GW in the theories of parity-violating gravity, we study the effects of velocity and amplitude birefringence on the GW waveforms. Decomposing the GWs into the circular polarization modes, the two birefringence effects exactly correspond to the modifications in phase and amplitude of GW waveforms respectively. We find that, for each circular polarization mode, the amplitude, phase and velocity of GW can be modified by both the parity-violating terms and parity-conserving terms in gravity. Therefore, in order to test the parity symmetry in gravity, we should compare the difference between two circular polarization modes, rather than measuring an individual mode. Combining two circular modes, we obtain the GW waveforms in the Fourier domain, and obtain the deviations from those in General Relativity. The GW waveforms derived in this paper are also applicable to the theories of parity-conserving gravity, which have the modified dispersion relations (e.g. massive gravity, double special relativity theory, extra-dimensional theories, etc), or/and have the modified friction terms (e.g. nonlocal gravity, gravitational theory with time-dependent Planck mass, etc).

Gravitational waves (GWs) provide an excellent opportunity to test the gravity in the strong gravitational fields. In this article, we calculate the waveform of GWs, produced by the coalescence of compact binaries, in an extension of the Chern-Simons gravity by including higher derivatives of the coupling scalar field. By comparing the two circular polarization modes, we find the effects of amplitude birefringence and velocity birefringence of GWs in their propagation caused by the parity violation in gravity, which are explicitly presented in the GW waveforms by the amplitude and phase modifications respectively. Combining the two modes, we obtain the GW waveforms in the Fourier domain, and find that the deviations from those in General Relativity are dominated by effects of velocity birefringence of GWs. In addition, we also map the effects of the parity violation on the waveform onto the parameterized post-Einsteinian (PPE) framework and identify explicitly the PPE parameters.

Supermassive binary black holes (SMBBHs) in galactic centers may radiate gravitational wave (GW) in the nano-Hertz frequency band, which are expected to be detected by pulsar timing arrays (PTAs) in the near future. GW signals from individual SMBBHs at cosmic distances, if detected by PTAs, are potentially powerful standard sirens that can be used to independently measure distances and thus put constraints on cosmological parameters. In this paper, we investigate the constraint that may be obtained on the equation of state ($w$) of dark energy by using those SMBBHs, expected to be detected by the PTAs in the Square Kilometre Array (SKA) era. By considering both the currently available SMBBH candidates and mock SMBBHs in the universe resulting from a simple galaxy major merger model, we find that $\sim 200$ to $3000$ SMBBHs with chirp mass $>10^9M_\odot$ are expected to be detected with signal-to-noise ratio $>10$ by SKA-PTA with conservative and optimistic settings and they can be used to put a constraint on $w$ to an uncertainty of $\Delta w\sim 0.02-0.1$. If further information on the mass and mass ratio of those SMBBHs can be provided by electromagnetic observations (e.g., chirp mass uncertainty $\lesssim 50\%$), the constraint may be further improved to $\lesssim 0.01$ level, as many more SMBBHs will be detected by SKA-PTA with relatively better distance measurements and can be used as the standard sirens.

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.

The similarity of the absolute luminosity profiles of Type Ia supernovae (SNIe), as one kind of distance indicator, has led their use in extragalactic astronomy as secondary standard candles. In general, the empirical relationship of SNIa on the absolute peak magnitude $M_{\rm B}$ is calibrated by Cepheid variables in the near distance scale and directly extrapolated to much farther distances. Two main problems arise. First, their calibration, in particular the determination of $M_{\rm B}$, depends on the empirical relationship of Cepheid variables, which suffers from various uncertainties. The second is related to the homogeneity of SNIa in their true $M_{\rm B}$, which is known to be poor in different environments. The observed GW signal of the coalescence of compact binary systems and their electromagnetic counterparts provide the novel and model-independent way to address these two problems. In the era of second-generation GW detectors, the low-redshift GW sources provide a novel method to calibrate the empirical relationship of SNIa, using their self-calibrated distances. Here, we use the event GW170817 to calibrate the empirical relationship in different low redshift ranges, and find that the calibration results are consistent with the ones derived from the Cepheid variables. Moreover, the uncertainties of $M_{\rm B}$ in both methods are also comparable. By the observations of third-generation GW detectors, GW sources can also be used to measure the values of $M_{\rm B}$ for the high-redshift SNIe, which provides a unique opportunity to study the dependence of $M_{\rm B}$ on the local environment, strength of gravity, and the intrinsic properties of the explosion, in addition to test the homogeneity of standard candles. We find that the uncertainties of $M_{\rm B}$ in both high and low redshifts are more than one order of magnitude smaller than the current accuracy.

In this paper, we study the gravitational waveforms, polarizations and radiation powers of the first relativistic triple systems PSR J0337 + 1715, observed in 2014, by using the post-Newtonian approximations to their lowest order. Although they cannot be observed either by current or next generation of the detectors, they do provide useful information to test different theories of gravity. In particular, we carry out the studies in three different theories, general relativity (GR), Brans-Dicke (BD) gravity, and Einstein-aether ($\ae$) theory. The tensor modes $h_{+}$ and $h_{\times}$ exist in all three theories and have almost equal amplitudes. Their frequencies are all peaked at two locations, $ f^{(+, \times)}_1 \simeq 0.0686656 \mu$Hz and $ f^{(+, \times)}_2 \simeq 14.2138 \mu$Hz, which are about twice of the outer and inner orbital frequencies of the triple system, as predicted in GR. In $\ae$-theory, all the six polarization modes are different from zero, but the breathing $h_b$ and longitudinal $h_L$ modes are not independent and also peaked at two frequencies, but at the frequencies, $ f^{(b, L, \ae)}_{1} \simeq 0.0457771 \mu$Hz and $f^{(b, L, \ae)}_{ 2} \simeq 7.09545 \mu$Hz. A similar phenomenon is also observed in BD gravity, in which only the three modes $h_{+},\; h_{\times}$ and $h_{b}$ exit. We also study the radiation powers, and find that the quadrupole emission in each of the three theories has almost the same amplitude, but the dipole emission can be as big as the quadrupole emission in $\ae$-theory. This can provide a very promising window to obtain severe constraints on $\ae$-theory by the multi-band gravitational wave astronomy.

The Analytical Blind Separation (ABS) method is dedicated to extract the CMB signal from its contaminants without any prior information about foregrounds. In this study, we apply the ABS method to simulated polarization sky maps in order to test its ability in the recovery of the CMB $E$- and $B$-mode power spectra. We simulate multi-frequency $Q$ and $U$ microwave sky maps in the presence of two polarized foreground components (synchrotron and thermal dust), together with an uncorrelated Gaussian distributed instrumental noise, measured with a hypothetical PICO-like experiment operating in 10-frequency bands in the range of 30 to 321 GHz. Both full and partial sky surveys are taken into account in our analysis. For the full sky case with the tensor-to-scalar ratios of $r=0$ and $r=0.05$, the ABS can accurately recover the $E$-mode and $B$-mode power spectra with a relative error below $20\%$ with respect to the input CMB power spectra for the full multipole range ($2 \leq \ell \leq 1050$), close to theoretical expectations. For the partial sky analysis with a E/B separation proposed by \cite{ebmixture10} to avoid "$E$-to-$B$" leakage, the results are comparable to the full sky case. We find that the ABS is able to estimate both CMB $E$-mode and $B$-mode power spectra within 1-$\sigma$ at all scales for full sky, and for most of the scales ($50 \leq \ell \leq 1050$) for partial sky. For low-$\ell$ multipoles ($\ell \leq 50$), a noticeable divergence occurs in the partial sky analysis, which might be due to uncertainties originated from both the ABS and the E/B separation. A more detailed analysis for the partial sky case is thus needed for estimating primordial $B$-modes unbiasedly and accurately.

Flat-spectrum radio quasar PKS~1229$-$02 with a knotty and asymmetric radio morphology was identified as the optical and radio counterpart of a $\gamma$-ray source. In this paper, we study the properties, e.g. morphology, opacity, polarization and kinematics of the jet in PKS~1229$-$02 using radio interferometry. With our results, we find that the knotty and asymmetric morphology of this source may probably shaped by the interaction between its anterograde jet and the nonuniform dense ambient medium. By reproducing a Spectral Energy Distribution of PKS~1229$-$02 with the obtained kinematic parameters, we find that the relativistic beaming effect in PKS~1229$-$02 is not strong enough to produce the reported $\gamma$-ray emission, i.e. PKS~1229$-$02 may not be a $\gamma$-ray AGN. The misidentification may probably due to the poor spatial resolution of the $\gamma$-ray detector of the previous generation.

Screened modified gravity (SMG) and Brans-Dicke (BD) gravity are typical examples of scalar-tensor theories with and without screening mechanisms, which can suppress the scalar field in dense regions. In this paper, we investigate the tests of time-varying gravitational constant $G$, gravitational dipole radiation, and Nordtvedt effect in BD and SMG theories, respectively. We place new constraints on these theories by combining Cassini experiment, lunar laser ranging (LLR) measurements, and pulsar observations from PSRs J1738$+$0333 and J0348$+$0432. We find that screening mechanism has important influence on theoretical constraints. The strongest, second, and weakest constraints on BD are from Cassini, pulsar, and LLR tests, respectively. The most stringent constraint on SMG comes from LLR measurements and improves the previous best constraint by more than seven orders of magnitude. We derive the bounds on the cosmological evolution of the scalar background in these theories using the time variation of $G$. The results of all tests agree well with general relativity (GR) and give more stringent constraints on the deviations from GR. Finally, as an example, we consider the chameleon model and derive the constraints on the model parameters.

Gravitationally bound hierarchies containing three or more components are very common in our Universe. In this paper we study {\em periodic} gravitational wave (GW) form, their polarizations, response function, its Fourier transform, and energy loss rate of a triple system through three different channels of radiation, the scalar, vector and tensor modes, in Einstein-aether theory of gravity. In the weak-field approximations and with the recently obtained constraints of the theory, we first analyze the energy loss rate of a binary system, and find that the dipole contributions from the scalar and vector modes could be of the order of ${\cal{O}}\left(c_{14}\right){\cal{O}}\left(G_Nm/d\right)^2$, where $c_{14} \; (\equiv c_{1} + c_{4})$ is constrained to $c_{14} \lesssim {\cal{O}}\left(10^{-5}\right)$ by current observations, where $c_i$'s are the four coupling constants of the theory. On the other hand, the "strong-field" effects for a binary system of neutron stars are about six orders lower than that of GR. So, in this paper we ignore these "strong-field" effects and first develop the general formulas to the lowest post-Newtonian order, by taking the coupling of the aether field with matter into account. Within this approximation, we find that the scalar breather mode and the scalar longitudinal mode are all suppressed by a factor of ${\cal{O}}\left(c_{14}\right)$ with respect to the transverse-traceless modes ($h_{+}$ and $h_{\times}$), while the vectorial modes $(h_{X}$ and $h_{Y}$) are suppressed by a factor of $c_{13} \lesssim {\cal{O}}\left(10^{-15}\right)$. Applying the general formulas to a triple system with periodic orbits, we find that the corresponding GW form, response function, and its Fourier transform depend sensitively on both the configuration of the triple system and their orientations with respect to the detectors.

The intrinsic peak luminosity of Type Ia supernovae (SNIa) depends on the value of Newton's gravitational constant $G$, through the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. If the luminosity distance can be independently determined, the SNIa can be treated as a tracker to constrain the possible time variation of $G$ in different redshift ranges. The gravitational-wave (GW) standard sirens, caused by the coalescence of binary neutron stars, provide a model-independent way to measure the distance of GW events, which can be used to determine the luminosity distances of SNIa by interpolation, provided the GW and SNIa samples have similar redshift ranges. We demonstrate that combining the GW observations of third-generation detectors with SNIa data provides a powerful and model-independent way to measure $G$ in a wide redshift range, which can constrain the ratio $G/G_0$, where $G$ and $G_0$ are respectively the values in the redshift ranges $z>0.1$ and $z<0.1$, at the level of $1.5\%$.

We use the Fisher information matrix to investigate the angular resolution and luminosity distance uncertainty for coalescing binary neutron stars (BNSs) and neutron star-black hole binaries (NSBHs) detected by the third-generation (3G) gravitational-wave (GW) detectors. Our study focuses on an individual 3G detector and a network of up to four 3G detectors at different locations including the US, Europe, China and Australia for the proposed Einstein Telescope (ET) and Cosmic Explorer (CE) detectors. We in particular examine the effect of the Earth's rotation, as GW signals from BNS and low mass NSBH systems could be hours long for 3G detectors. We find that, a time-dependent antenna beam-pattern function can help better localize BNS and NSBH sources, especially those edge-on ones. The medium angular resolution for one ET-D detector is around 150 deg$^2$ for BNSs at a redshift of $z=0.1$. The medium angular resolution for a network of two CE detectors in the US and Europe respectively is around 20 deg$^2$ at $z=0.2$ for the simulated BNS and NSBH samples. While for a network of two ET-D detectors, the similar angular resolution can be achieved at a much higher redshift of $z=0.5$. The angular resolution of a network of three detectors is mainly determined by the baselines between detectors regardless of the CE or ET detector type. We discuss the implications of our results to constrain the Hubble constant $H_0$, the deceleration parameter $q_0$ and the equation-of-state (EoS) of dark energy. We find that in general, if 10 BNSs or NSBHs at $z=0.1$ with known redshifts are detected, $H_0$ can be measured with an accuracy of $0.9\%$. If 1000 face-on BNSs at $z<2$ are detected with known redshifts, we are able to achieve $\Delta q_0=0.002$, or $\Delta w_0=0.03$ and $\Delta w_a=0.2$ for dark energy.(Abridged version).

There has been speculation of a class of relativistic explosions with an initial Lorentz factor smaller than that of classical Gamma-Ray Bursts (GRBs). These "dirty fireballs" would lack prompt GRB emission but could be pursued via their optical afterglow, appearing as transients that fade overnight. Here we report a search for such transients (transients that fade by 5-$\sigma$ in magnitude overnight) in four years of archival photometric data from the intermediate Palomar Transient Factory (iPTF). Our search criteria yielded 45 candidates. Of these, two were afterglows to GRBs that had been found in dedicated follow-up observations to triggers from the Fermi GRB Monitor (GBM). Another (iPTF14yb; Cenko et al. 2015) was a GRB afterglow discovered serendipitously. Two were spurious artifacts of reference image subtraction and one was an asteroid. The remaining 37 candidates have red stellar counterparts in external catalogs. The photometric and spectroscopic properties of the counterparts identify these transients as strong flares from M dwarfs of spectral type M3-M7 at distances of d ~ 0.15-2.1 kpc; two counterparts were already spectroscopically classified as late-type M stars. With iPTF14yb as the only confirmed relativistic outflow discovered independently of a high-energy trigger, we constrain the all-sky rate of transients that peak at m = 18 and fade by $\Delta$2 mag in $\Delta$3 hr to be 680 per year with a 68% confidence interval of 119-2236 per year. This implies that the rate of visible dirty fireballs is at most comparable to that of the known population of long-duration GRBs.

We present high resolution (sub-mas) VLBI maps of nuclear H2O megamasers for seven galaxies. In UGC6093, the well-aligned systemic masers and high-velocity masers originate in an edge-on, flat disk and we determine the mass of the central SMBH to be M_SMBH = 2.58*10^7Msun(+-7%). For J1346+5228, the distribution of masers is consistent with a disk, but the faint high-velocity masers are only marginally detected, and we constrain the mass of the SMBH to be in the range 1.5-2.0*10^7Msun. The origin of the masers in Mrk1210 is less clear, as the systemic and high-velocity masers are misaligned and show a disorganized velocity structure. We present one possible model in which the masers originate in a tilted, warped disk, but we do not rule out the possibility of other explanations including outflow masers. In NGC6926, we detect a set of redshifted masers, clustered within a pc of each other, and a single blueshifted maser about 4.4pc away, an offset that would be unusually large for a maser disk system. Nevertheless, if it is a disk system, we estimate the enclosed mass to be M_SMBH<4.8*10^7 Msun . For NGC5793, we detect redshifted masers spaced about 1.4pc from a clustered set of blueshifted features. The orientation of the structure supports a disk scenario as suggested by Hagiwara et al.(2001). We estimate the enclosed mass to be M SMBH<1.3*10^7 Msun. For NGC2824 and J0350-0127, the masers may be associated with pc or sub-pc scale jets or outflows.

The $\gamma$-ray bright blazar CTA 102 is studied using imaging (new 15 GHz and archival 43 GHz Very Long Baseline Array, VLBA data) and time variable optical flux density, polarization degree and electric vector position angle (EVPA) spanning between 2015 June 1 and 2016 October 1, covering a prominent $\gamma$-ray flare during 2016 January. The pc-scale jet indicates expansion with oscillatory features upto 17 mas. Component proper motions are in the range 0.04 - 0.33 mas/yr with acceleration upto 1.2 mas followed by a slowing down beyond 1.5 mas. A jet bulk Lorentz factor $\geq$ 17.5, position angle of 128.3 degrees, inclination angle $\leq$ 6.6 degrees and intrinsic half opening angle $\leq$ 1.8 degrees are derived from the VLBA data. These inferences are employed in a helical jet model to infer long term variability in flux density, polarization degree, EVPA and a rotation of the Stokes Q and U parameters. A core distance of $r_{\rm core,43 \ GHz}$ = 22.9 pc, and a magnetic field strength at 1 pc and the core location of 1.57 G and 0.07 G respectively are inferred using the core shift method. The study is useful in the context of estimating jet parameters and in offering clues to distinguish mechanisms responsible for variability over different timescales.

The $f(R)$ gravity can be cast into the form of a scalar-tensor theory, and scalar degree of freedom can be suppressed in high-density regions by the chameleon mechanism. In this article, for the general $f(R)$ gravity, using a scalar-tensor representation with the chameleon mechanism, we calculate the parameterized post-Newtonian parameters $\gamma$ and $\beta$, the effective gravitational constant $G_{\rm eff}$, and the effective cosmological constant $\Lambda_{\rm eff}$. In addition, for the general $f(R)$ gravity, we also calculate the rate of orbital period decay of the binary system due to gravitational radiation. Then we apply these results to specific $f(R)$ models (Hu-Sawicki model, Tsujikawa model and Starobinsky model) and derive the constraints on the model parameters by combining the observations in solar system, cosmological scales and the binary systems.

Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.

In this paper, we perform a forecast analysis to test the capacity of future baryon acoustic oscillation (BAO) and cosmic microwave background (CMB) experiments to constrain phenomenological interacting dark energy models using the Fisher matrix formalism. We consider a Euclid-like experiment, in which BAO measurements is one of the main goals, to constrain the cosmological parameters of alternative cosmological models. Moreover, additional experimental probes can more efficiently provide information on the parameters forecast, justifying also the inclusion in the analysis of a future ground-based CMB experiment mainly designed to measure the polarization signal with high precision. In the interacting dark energy scenario, a coupling between dark matter and dark energy modifies the conservation equations such that the fluid equations for both constituents are conserved as the total energy density of the dark sector. In this context, we consider three phenomenological models which have been deeply investigated in literature over the past years. We find that the combination of both CMB and BAO information can break degeneracies among the dark sector parameters for all three models, although to different extents. We found powerful constraints on, for example, the coupling constant when comparing it with present limits for two of the models, and their future statistical 3-$\sigma$ bounds could potentially exclude the null interaction for the combination of probes that is considered. However, for one of the models, the constraint on the coupling parameter does not improve the present result (achieved using a large combination of surveys), and a larger combination of probes appears to be necessary to eventually claim whether or not interaction is favored in that context.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

In both WMAP and Planck observations on the temperature anisotropy of cosmic microwave background (CMB) radiation a number of large-scale anomalies were discovered in the past years, including the CMB parity asymmetry in the low multipoles. By defining a directional statistics, we find that the CMB parity asymmetry is directional dependent, and the preferred axis is stable, which means that it is independent of the chosen CMB map, the definition of the statistic, or the CMB masks. Meanwhile, we find that this preferred axis strongly aligns with those of the CMB quadrupole, octopole, as well as those of other large-scale observations. In addition, all of them aligns with the CMB kinematic dipole, which hints to the non-cosmological origin of these directional anomalies in cosmological observations.

Measuring the imprint of primordial gravitational waves in the cosmic microwave background (CMB) polarisation field is one of the main goals in modern cosmology. However, the so called $B$-mode polarisation can be generated by different sources besides the primary one predicted by inflationary theories, known as secondary $B$-mode signal. Among them, CMB lensing and astrophysical foregrounds play an important role. Moreover, a partial sky analysis leads to a leakage between $E$-modes and $B$-modes. In this article, we use the well known Minkowski functionals (MF) statistics to study the significance of this leakage in the CMB lensing $B$-mode signal. We find that the MF can detect the $E$-to-$B$ leakage contamination, thus it should not be neglected in future CMB data analysis.

Gravitational radiation is an excellent field for testing theories of gravity in strong gravitational fields. The current observations on the gravitational-wave (GW) bursts by LIGO have already placed various constraints on the alternative theories of gravity. In this paper, we investigate the possible bounds which could be placed on the Brans-Dicke gravity using GW detection from inspiralling compact binaries with the proposed Einstein Telescope, a third-generation GW detector. We first calculate in details the waveforms of gravitational radiation in the lowest post-Newtonian approximation, including the tensor and scalar fields, which can be divided into the three polarization modes, i.e. "plus mode", "cross mode" and "breathing mode". Applying the stationary phase approximation, we obtain their Fourier transforms, and derive the correction terms in amplitude, phase and polarization of GWs, relative to the corresponding results in General Relativity. Imposing the noise level of Einstein Telescope, we find that the GW detection from inspiralling compact binaries, composed of a neutron star and a black hole, can place stringent constraints on the Brans-Dicke gravity. The bound on the coupling constant $\omega_{\rm BD}$ depends on the mass, sky-position, orbital angle, polarization angle, luminosity distance, redshift distribution and total observed number $N_{\rm GW}$ of the binary systems. Taking into account all the burst events up to redshift $z=5$, we find that the bound could be $\omega_{\rm BD}\gtrsim 10^{6}\times(N_{\rm GW}/10^4)^{1/2}$. Even for the conservative estimation with $10^{4}$ observed events, the bound is still more than one order tighter than the current limit from Solar System experiments. So, we conclude that Einstein Telescope will provide a powerful platform to test alternative theories of gravity.

Screened modified gravity (SMG) is a kind of scalar-tensor theory with screening mechanisms, which can suppress the fifth force in dense regions and allow theories to evade the solar system and laboratory tests. In this paper, we investigate how the screening mechanisms in SMG affect the gravitational radiation damping effects, calculate in detail the rate of the energy loss due to the emission of tensor and scalar gravitational radiations, and derive their contributions to the change in the orbital period of the binary system. We find that the scalar radiation depends on the screened parameters and the propagation speed of scalar waves, and the scalar dipole radiation dominates the orbital decay of the binary system. For strongly self-gravitating bodies, all effects of scalar sector are strongly suppressed by the screening mechanisms in SMG. By comparing our results to observations of binary system PSR J1738+0333, we place the stringent constraints on the screening mechanisms in SMG. As an application of these results, we focus on three specific models of SMG (chameleon, symmetron, and dilaton), and derive the constraints on the model parameters, respectively.

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

Parity violation in the cosmic microwave background (CMB) radiation has been confirmed by recent Planck observation. In this paper, we extend our previous work [P. Naselsky et al. Astrophys. J. {\bf 749}, 31 (2012)] on the directional properties of CMB parity asymmetry by considering the Planck data of the CMB temperature anisotropy. We define six kinds of the directional statistics and find that they all indicate the odd-parity preference of CMB data. In addition, we find the preferred axes of all these statistics strongly correlate with the preferred axes of the CMB kinematic dipole, quadrupole, and octopole. The alignment between them is confirmed at more than $3\sigma$ confidence level, which implies that the CMB parity asymmetry, and the anomalies of the CMB quadrupole and octopole, may have the common physical, contaminated or systematic origin. In addition, they all should be connected with the possible contamination of the residual dipole component.

We investigate the primordial gravitational waves (PGWs) in the general scenario where the inflation is preceded by a pre-inflationary stage with the effective equation of state $w$. Comparing with the results in the usual inflationary models, the power spectrum of PGWs is modified in two aspects: One is the mixture of the perturbation modes caused by he presence of the pre-inflationary period, and the other is the thermal initial state formed at the Planck era of the early Universe. By investigating the observational imprints of these modifications on the B-mode polarization of cosmic microwave background (CMB) radiation, we obtain the constraints on the conformal temperature of the thermal gravitational-wave background $T<5.01\times 10^{-4}$Mpc$^{-1}$ and a tensor-to-scalar ratio $r<0.084$ ($95\%$ confident level), which follows the bounds on total number of e-folds $N>63.5$ for the model with $w=1/3$, and $N>65.7$ for that with $w=1$. By taking into account various noises and the foreground radiations, we forecast the detection possibility of the thermal gravitational-wave background by the future CMBPol mission, and find that if $r>0.01$, the detection is possible as long as $T>1.5\times 10^{-4}$Mpc$^{-1}$. However, the effect of different $w$ is quite small, and it seems impossible to determine its value from the potential observations of CMBPol mission.

One of the main goals of modern cosmology is to search for primordial gravitational waves by looking on their imprints in the B-type polarization in the cosmic microwave background radiation. However, this signal is contaminated by various sources, including cosmic weak lensing, foreground radiations, instrumental noises, as well as the E-to-B leakage caused by the partial sky surveys, which should be well understood to avoid the misinterpretation of the observed data. In this paper, we adopt the E/B decomposition method suggested by Smith in 2006, and study the imprints of E-to-B leakage residuals in the constructed B-type polarization maps, $\mathcal{B}(\hat{n})$, by employing various statistical tools. We find that the effects of E-to-B leakage are negligible for the $\mathcal{B}$-mode power spectrum, as well as the skewness and kurtosis analyses of $\mathcal{B}$-maps. However, if employing the morphological statistical tools, including Minkowski functionals and/or Betti numbers, we find the effect of leakage can be detected at very high confidence level, which shows that in the morphological analysis, the leakage can play a significant role as a contaminant for measuring the primordial B-mode signal and must be taken into account for a correct explanation of the data.

The quasar 3C~286 is one of two compact steep spectrum sources detected by the {\it Fermi}/LAT. Here, we investigate the radio properties of the parsec(pc)-scale jet and its (possible) association with the $\gamma$-ray emission in 3C~286. The Very Long Baseline Interferometry (VLBI) images at various frequencies reveal a one-sided core--jet structure extending to the southwest at a projected distance of $\sim$1 kpc. The component at the jet base showing an inverted spectrum is identified as the core, with a mean brightness temperature of $2.8\times 10^{9}$~K. The jet bends at about 600 pc (in projection) away from the core, from a position angle of $-135^\circ$ to $-115^\circ$. Based on the available VLBI data, we inferred the proper motion speed of the inner jet as $0.013 \pm 0.011$ mas yr$^{-1}$ ($\beta_{\rm app} = 0.6 \pm 0.5$), corresponding to a jet speed of about $0.5\,c$ at an inclination angle of $48^\circ$ between the jet and the line of sight of the observer. The brightness temperature, jet speed and Lorentz factor are much lower than those of $\gamma$-ray-emitting blazars, implying that the pc-scale jet in 3C~286 is mildly relativistic. Unlike blazars in which $\gamma$-ray emission is in general thought to originate from the beamed innermost jet, the location and mechanism of $\gamma$-ray emission in 3C~286 may be different as indicated by the current radio data. Multi-band spectrum fitting may offer a complementary diagnostic clue of the $\gamma$-ray production mechanism in this source.

Most of the gamma-ray emitting active galactic nuclei (AGN) are blazars, although there is still a small fraction of non-blazar AGN in the Fermi/LAT catalog. Among these misaligned gamma-ray-emitting AGN, a few can be classified as Compact Symmetric Objects (CSOs). In contrast to blazars in which gamma-ray emission is generally thought to originate from highly beamed relativistic jets, the source of gamma-ray emission in unbeamed CSOs remains an open question. The rarity of the gamma-ray emitting CSOs is a mystery as well. Here we present the radio properties of the gamma-ray CSO candidate 2234+282.

As part of the Megamaser Cosmology Project (MCP), we present VLBI maps of nuclear water masers toward five galaxies. The masers originate in sub-parsec circumnuclear disks. For three of the galaxies, we fit Keplerian rotation curves to estimate their supermassive black hole (SMBH) masses, and determine (2.9 $\pm$ 0.3) $\times~10^{6}M_\odot$ for J0437+2456, (1.7 $\pm$ 0.1) $\times~10^{7}M_\odot$ for ESO 558$-$G009, and (1.1 $\pm$ 0.2) $\times~10^{7}M_\odot$ for NGC 5495. In the other two galaxies, Mrk 1029 and NGC 1320, the geometry and dynamics are more complicated and preclude robust black hole mass estimates. Including our new results, we compiled a list of 15 VLBI-confirmed disk maser galaxies with robust SMBH mass measurements. With this sample, we confirm the empirical relation of $R_{out} \propto 0.3 M_{SMBH}$ reported in Wardle & Yusef-Zadeh (2012). We also find a tentative correlation between maser disk outer radii and WISE luminosity. We find no correlations of maser disk size with X-ray 2-10 keV luminosity or [O III] luminosity.

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

We use new precision measurements of black hole masses from water megamaser disks to investigate scaling relations between macroscopic galaxy properties and supermassive black hole (BH) mass. The megamaser-derived BH masses span 10^6-10^8 M_sun, while all the galaxy properties that we examine (including stellar mass, central mass density, central velocity dispersion) lie within a narrow range. Thus, no galaxy property correlates tightly with M_BH in ~L* spiral galaxies. Of them all, stellar velocity dispersion provides the tightest relation, but at fixed sigma* the mean megamaser M_BH are offset by -0.6+/-0.1 dex relative to early-type galaxies. Spiral galaxies with non-maser dynamical BH masses do not show this offset. At low mass, we do not yet know the full distribution of BH mass at fixed galaxy property; the non-maser dynamical measurements may miss the low-mass end of the BH distribution due to inability to resolve the spheres of influence and/or megamasers may preferentially occur in lower-mass BHs.

The detection of the magnetic type $B$-mode polarization is the main goal of future cosmic microwave background (CMB) experiments. In the standard model, the $B$-mode map is a strongly non-gaussian field due to the lensed component. Besides the two-point correlation function, the other statistics are also very important to dig the information of the polarization map. In this paper, we employ the Minkowski functionals to study the morphological properties of the lensed $B$-mode maps. We find that the deviations from Gaussianity are very significant for both full and partial-sky surveys. As an application of the analysis, we investigate the morphological imprints of the foreground residuals in the $B$-mode map. We find that even for very tiny foreground residuals, the effects on the map can be detected by the Minkowski functional analysis. Therefore, it provides a complementary way to investigate the foreground contaminations in the CMB studies.

Screened modified gravity (SMG) is a kind of scalar-tensor theories with screening mechanisms, which can generate screening effect to suppress the fifth force in high density environments and pass the solar system tests. Meanwhile, the potential of scalar field in the theories can drive the acceleration of the late universe. In this paper, we calculate the parameterized post-Newtonian (PPN) parameters $\gamma$ and $\beta$, the effective gravitational constant $G_{\rm eff}$ and the effective cosmological constant $\Lambda$ for SMG with a general potential $V$ and coupling function $A$. The dependence of these parameters on the model parameters of SMG and/or the physical properties of the source object are clearly presented. As an application of these results, we focus on three specific theories of SMG (chameleon, symmetron and dilaton models). Using the formulae to calculate their PPN parameters and cosmological constant, we derive the constraints on the model parameters by combining the observations on solar system and cosmological scales.

The foundation of modern cosmology relies on the so-called cosmological principle which states an homogeneous and isotropic distribution of matter in the universe on large scales. However, recent observations, such as the temperature anisotropy of the cosmic microwave background (CMB) radiation, the motion of galaxies in the universe, the polarization of quasars and the acceleration of the cosmic expansion, indicate preferred directions in the sky. If these directions have a cosmological origin, the cosmological principle would be violated, and modern cosmology should be reconsidered. In this paper, by considering the preferred axis in the CMB parity violation, we find that it coincides with the preferred axes in CMB quadrupole and CMB octopole, and they all align with the direction of the CMB kinematic dipole. In addition, the preferred directions in the velocity flows, quasar alignment, anisotropy of the cosmic acceleration, the handedness of spiral galaxies, and the angular distribution of the fine-structure constant are also claimed to be aligned with the CMB kinematic dipole. Since CMB dipole was confirmed to be caused by the motion of our local group of galaxies relative to the reference frame of the CMB, the coincidence of all these preferred directions hints that these anomalies have a common origin, which is not cosmological or due to a gravitational effect. The systematical or contaminative errors in observation or in data analysis, which can be directly related to the motion of our local group of galaxies, can play an important role in explaining the anomalies.

Both WMAP and Planck data show a significant odd-multipole preference in the large scales of the cosmic microwave background (CMB) temperature anisotropies. If this pattern originates from cosmological effects, then it can be considered a crucial clue for a violation in the cosmological principle. By defining various direction dependent statistics in the full-sky Planck 2015 maps (see, for instance, Naselsky et al. (2012); W. Zhao (2014)), we found that the CMB parity asymmetry has a preferred direction, which is independent of the choices of the statistics. In particular, this preferred axis is strongly aligned with those in the CMB quadrupole and octopole, as well as that in the CMB kinematic dipole, which hints to their non-cosmological origin. In realistic observations, the foreground residuals are inevitable, and should be properly masked out in order to avoid possible misinterpretation of the results. In this paper, we extend our previous analyses to the masked Planck 2015 data. By defining a similar direction dependent statistic in the masked map, we find a preferred direction of the CMB parity asymmetry, in which the axis also coincides with that found in the full-sky analysis. Therefore, our conclusions on the CMB parity violation and its directional properties are confirmed.

We report on the discovery of broad Balmer absorption lines variability in the QSO SDSS J125942.80+121312.6, based on the optical and near-infrared spectra taken from the SDSS-I, SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and TripleSpec observations over a timescale of 5.8 years in the QSO's rest-frame. The blueshifted absorption profile of H$\beta$ shows a variation of more than 5$\sigma$ at a high velocity portion ($>3000\ \mathrm{km\ s}^{-1}$) of the trough. We perform a detailed analysis for the physical conditions of the absorber using Balmer lines as well as metastable He I and optical Fe II absorptions ($\lambda 4233$ from b$^4$P$_{5/2}$ level and $\lambda 5169$ from a$^6$S$_{5/2}$) at the same velocity. These Fe II lines are identified in the QSO spectra for the first time. According to the photoionization simulations, we estimate a gas density of $n(\mathrm{H})\approx 10^{9.1}\ \mathrm{cm}^{-3}$ and a column density of $N_{\mathrm{col}}(\mathrm{H})\approx 10^{23}\ \mathrm{cm}^{-2}$ for the BOSS data, but the model fails to predict the variations of ionic column densities between the SDSS and BOSS observations if changes in ionizing flux are assumed. We thus propose transverse motion of the absorbing gas being the cause of the observed broad Balmer absorption line variability. In fact, we find that the changes in covering factors of the absorber can well-reproduce all of the observed variations. The absorber is estimated $\sim 0.94$ pc away from the central engine, which is where the outflow likely experiences deceleration due to the collision with the surrounding medium. This scheme is consistent with the argument that LoBAL QSOs may represent the transition from obscured star-forming galaxies to classic QSOs.

We explore a plausible mechanism that the hemispherical power asymmetry in the CMB is produced by the spatial variation of the primordial sound speed parameter. We suggest that in a generalized approach of the $\delta N$ formalism the local e-folding number may depend on some other primordial parameters besides the initial values of inflaton. Here the $\delta N$ formalism is extended by considering the effects of a spatially varying sound speed parameter caused by a super-Hubble perturbation of a light field. Using this generalized $\delta N$ formalism, we systematically calculate the asymmetric primordial spectrum in the model of multi-speed inflation by taking into account the constraints of primordial non-Gaussianities. We further discuss specific model constraints, and the corresponding asymmetry amplitudes are found to be scale-dependent, which can accommodate current observations of the power asymmetry at different length scales.

Relic gravitational waves (RGWs) generated in the early Universe form a stochastic GW background, which can be directly probed by measuring the timing residuals of millisecond pulsars. In this paper, we investigate the constraints on the RGWs and on the inflationary parameters by the observations of current and potential future pulsar timing arrays. In particular, we focus on effects of various cosmic phase transitions (e.g. $e^{+}e^{-}$ annihilation, QCD transition and SUSY breaking) and relativistic free-streaming gases (neutrinos and dark fluids) in the general scenario of the early Universe, which have been neglected in the previous works. We find that the phase transitions can significantly damp the RGWs in the sensitive frequency range of pulsar timing arrays, and the upper limits of tensor-to-scalar ratio $r$ increase by a factor $\sim 2$ for both current and future observations. However, the effects of free-steaming neutrinos and dark fluids are all too small to be detected. Meanwhile, we find that, if the effective equation of state $w$ in the early Universe is larger than $1/3$, i.e. deviating from the standard hot big bang scenario, the detection of RGWs by pulsar timing arrays becomes much more promising.

As part of the Megamaser Cosmology Project (MCP), here we present a new geometric distance measurement to the megamaser galaxy NGC 5765b. Through a series of VLBI observations, we have confirmed the water masers trace a thin, sub-parsec Keplerian disk around the nucleus, implying an enclosed mass of 4.55 $\pm$ 0.40 $\times~10^{7}M_\odot$. Meanwhile, from single dish monitoring of the maser spectra over two years, we measured the secular drifts of maser features near the systemic velocity of the galaxy with rates between 0.5 and 1.2 km s$^{-1}$ yr$^{-1}$. Fitting a warped, thin disk model to these measurements, we determine a Hubble Constant $H_{0}$ of 66.0 $\pm$ 6.0 km s$^{-1}$ Mpc$^{-1}$ with the angular-diameter distance to NGC 5765b of 126.3 $\pm$ 11.6 Mpc. Apart from the distance measurement, we also investigate some physical properties related to the maser disk in NGC 5765b. The high-velocity features are spatially distributed into several clumps, which may indicate the existence of a spiral density wave associated with the accretion disk. For the red-shifted features, the envelope defined by the peak maser intensities increases with radius. The profile of the systemic masers in NGC 5765b is smooth and shows almost no structural changes over the two years of monitoring time, which differs from the more variable case of NGC 4258.

The anisotropies of the B-mode polarization in the cosmic microwave background radiation play a crucial role for the study of the very early Universe. However, in the real observation, the mixture of the E-mode and B-mode can be caused by the partial sky surveys, which must be separated before applied to the cosmological explanation. The separation method developed by Smith (\citealt{PhysRevD.74.083002}) has been widely adopted, where the edge of the top-hat mask should be smoothed to avoid the numerical errors. In this paper, we compare three different smoothing methods, and investigate the leakage residuals of the E-B mixture. We find that, if the less information loss is needed and the smaller region is smoothed in the analysis, the \textit{sin}- and \textit{cos}-smoothing methods are better. However, if we need a clean constructed B-mode map, the larger region around the mask edge should be smoothed. In this case, the \textit{Gaussian}-smoothing method becomes much better. In addition, we find that the leakage caused by the numerical errors in the \textit{Gaussian}-smoothing method mostly concentrates on two bands, which is quite easy to be reduced for the further E-B separations.

Both WMAP and PLANCK missions reported the extremely Cold Spot (CS) centered at Galactic coordinate ($l=209^{\circ}$, $b=-57^{\circ}$) in CMB map. In this paper, we study the local non-Gaussianity of CS by defining the local Minkowski functions. We find that the third Minkowski function $\nu_2$ is quite sensitive to the non-Gaussianity caused by CS. Compared with the random Gaussian simulations, WMAP CS deviates from Gaussianity at more than 99% confident level at the scale $R\sim10^{\circ}$. Meanwhile, we find that cosmic texture provides an excellent explanation for these anomalies related to WMAP CS, which could be further tested by the future polarization data.

The measured CMB angular distribution shows a great consistency with the LCDM model. However, isotropy violations were reported in CMB temperature maps of both WMAP and Planck data. We investigate the influence of different masks employed in the analysis of CMB angular distribution, in particular in the excess of power in the Southeastern quadrant (SEQ) and the lack of power in the Northeastern quadrant (NEQ). We compare the two-point correlation function (TPCF) computed for each quadrant of the CMB foreground-cleaned temperature maps to 1000 simulations generated assuming the LCDM best-fit power spectrum using four different masks. In addition to the quadrants, we computed the TPCF for circular regions in the map where the excess and lack of power are present. We also compare the effect of Galactic cuts in the TPCF calculations as compared to the simulations. We found consistent results for three masks, namely mask-rulerminimal, U73 and U66. The results indicate that the excess of power in the SEQ tends to vanish as the portion of the sky covered by the mask increases and the lack of power in the NEQ remains virtually unchanged. When UT78 mask is applied, the NEQ becomes no longer anomalous and the excess of power in the SEQ becomes the most significant one among the masks. Nevertheless, the asymmetry between the SEQ and NEQ is independent of the mask and it is in disagreement with the isotropic model with at least 95% C.L. We find that UT78 is in disagreement with the other analysed masks, specially considering the SEQ and the NEQ individual analysis. Most importantly, the use of UT78 washes out the anomaly in the NEQ. Furthermore, we found excess of kurtosis, compared with simulations, in the NEQ for the regions not masked by UT78 but masked by the other masks, indicating that the previous result could be due to non-removed residual foregrounds by UT78.

By taking into account the contamination of foreground radiations, we employ the Fisher matrix to forecast the future sensitivity on the tilt of power spectrum of primordial tensor perturbations for several ground-based (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloon-borne (EBEX, Spider) and satellite (CMBPol, COrE, LiteBIRD) experiments of B-mode polarizations. For the fiducial model $n_t=0$, our results show that the satellite experiments give good sensitivity on the tensor tilt $n_t$ to the level $\sigma_{n_t}\lesssim0.1$ for $r\gtrsim2\times10^{-3}$, while the ground-based and balloon-borne experiments give worse sensitivity. By considering the BICEP2/Keck Array and Planck (BKP) constraint on the tensor-to-scalar ratio $r$, we see that it is impossible for these experiments to test the consistency relation $n_t=-r/8$ in the canonical single-field slow-roll inflation models.

Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.

The B-mode polarization of the cosmic microwave background (CMB) radiation is an excellent information channel for the detection of relic gravitational waves. However, the detection is contaminated by the B-mode polarization generated by some other effects. In this paper, we discuss the contaminations caused by the cosmological birefringence, which converts the CMB E-mode to the B-mode, and forms the effective noise for the detection of gravitational waves. We find that this contamination is significant, if the rotation angle is large. However, this kind of B-mode can be properly de-rotated, and the effective noises can be greatly reduced. We find that, comparing with the contaminations caused by cosmic weak lensing, the residual polarization generated by the cosmological birefringence is negligible for the detection of relic gravitational waves in the CMB.

The inorganic CsI(Tl) crystal scintillator is a candidate served as an anti-compton detector for China Dark matter EXperiment (CDEX). Studying the intrinsic radiopurity of CsI(Tl) crystal is an issue of major importance. The timing, energy and spatial correlations, as well as the capability of pulse shape discrimination provide powerful methods for the measurement of intrinsic radiopurities. The experimental design, detector performance and event-selection algorithms are described. A total of 359*3 kg-days data from three prototypes of CsI(Tl) crystals were taken at China Jinping underground laboratory where offers a good shielding environment. The contamination levels of internal isotopes from Cs-137, Th-232 and U-238 series, as well as the upper bounds of U-235 series are reported. Identification of the whole alpha8 peaks from U/Th decay chains and derivation of those corresponding quenching factors are achieved.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

Recently Background Imaging of Cosmic Extragalactic Polarization (B2) discovered the relic gravitational waves at $7.0\sigma$ confidence level. However, the other cosmic microwave background (CMB) data, for example Planck data released in 2013 (P13), prefer a much smaller amplitude of the primordial gravitational waves spectrum if a power-law spectrum of adiabatic scalar perturbations is assumed in the six-parameter $\Lambda$CDM cosmology. In this paper, we explore whether the $w$CDM model and the running spectral index can relax the tension between B2 and other CMB data. In particular, we find that a positive running of running of spectral index is preferred at $1.7\sigma$ level from the combination of B2, P13 and WMAP Polarization data.

The cosmological birefringence caused by the coupling between the cosmic scalar field and the cosmic microwave background (CMB) photons through the Chern-Simons term can rotate the polarization planes of the photons, and mix the CMB E-mode and B-mode polarizations. The rotation angle induced by the dynamical scalar field can be separated into the isotropic background part and the anisotropic fluctuations. The effect of the background part has been be studied in the previous work (Zhao & Li, arXiv:1402.4324). In this paper, we focus on the influence of the anisotropies of the rotation angle. We first assume that the cosmic scalar field is massless, consistent with other works, we find that the rotation spectrum can be quite large, which may be detected by the potential CMB observations. However, if the scalar field is identified as the quintessence field, by detailed discussing both the entropy and adiabatic perturbation modes for the first time, we find that the anisotropies of the rotation angle are always too small to be detectable.In addition, as the main goal of this paper, we investigate the effect of rotated polarization power spectrum on the detection of relic gravitational waves. we find that, the rotated B-mode polarization could be fairly large, and comparable with those generated by the gravitational waves. This forms a new contamination for the detection of relic gravitational waves in the CMB. In particular, we also propose the method to reconstruct and subtract the rotated B-mode polarization, by which the residuals become negligible for the gravitational-wave detection.

Relic gravitational waves (RGWs) leave well-understood imprints on the anisotropies in the temperature and polarization of cosmic microwave background (CMB) radiation. In the TT and TE information channels, which have been well observed by WMAP and Planck missions, RGWs compete with density perturbations mainly at low multipoles. It is dangerous to include high-multipole CMB data in the search for gravitational waves, as the spectral indices may not be constants. In this paper, we repeat our previous work [W.Zhao & L.P.Grishchuk, Phys.Rev.D {\bf 82}, 123008 (2010)] by utilizing the Planck TT and WMAP TE data in the low-multipole range $\ell\le100$. We find that our previous result is confirmed {with higher confidence}. The constraint on the tensor-to-scalar ratio from Planck TT and WMAP TE data is $r\in [0.06,~0.60]$ (68% C.L.) with the maximum likelihood at around $r\sim 0.2$. Correspondingly, the spectral index at the pivot wavenumber $k_*=0.002$Mpc$^{-1}$ is $n_s=1.13^{+0.07}_{-0.08}$, which is larger than 1 at more than $1\sigma$ level. So, we conclude that the new released CMB data indicate a stronger hint for the RGWs with the amplitude $r\sim 0.2$, which is hopeful to be confirmed by the imminent BICEP and Planck polarization data.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

In the non-standard model of relic gravitational waves (RGWs) generated in the early universe, the theoretical spectrum of is mainly described by an amplitude $r$ and a spectral index $\beta$, the latter usually being determined by the slope of the inflation potential. Pulsar timing arrays (PTAs) data have imposed constraints on the amplitude of strain spectrum for a power-law form as a phenomenological model. Applying these constraints to a generic, theoretical spectrum with $r$ and $\beta$ as independent parameters, we convert the PTAs constraint into an upper bound on the index $\beta$, which turns out to be less stringent than those upper bounds from BBN, CMB, and LIGO/VIRGO, respectively. Moreover, it is found that PTAs constrain the non-standard RGWs more stringent than the standard RGWs. If the condition of the quantum normalization is imposed upon a theoretical spectrum of RGWs, $r$ and $\beta$ become related. With this condition, a minimum requirement of the horizon size during inflation is greater than the Planck length results in an upper bound on $\beta$, which is comparable in magnitude to that by PTAs. When both PTAs and the quantum normalization are applied to a theoretical spectrum of RGWs, constraints can be obtained for other cosmic processes of the early universe, such as the reheating, a process less understood observationally so far. The resulting constraint is consistent with the slow-roll, massive scalar inflation model. The future SKA will be able to constrain RGWs further and might even detect RGWs, rendering an important probe to the very early universe.

In this paper, we study the effects of parity violation on non-gaussianities of primordial gravitational waves in the framework of Ho\v{r}ava-Lifshitz theory of gravity, in which high-order spatial derivative operators, including the ones violating parity, generically appear. By calculating the three point function, we find that the leading-order contributions to the non-gaussianities come from the usual second-order derivative terms, which produce the same bispectrum as that found in general relativity. The contributions from high-order spatial n-th derivative terms are always suppressed by a factor $(H/M_*)^{n-2} \; (n \ge 3)$, where $H$ denotes the inflationary energy and $M_*$ the suppression mass scale of the high-order spatial derivative operators of the theory. Therefore, the next leading-order contributions come from the 3-dimensional gravitational Chern-Simons term. With some reasonable arguments, it is shown that this 3-dimensional operator is the only one that violates the parity and in the meantime has non-vanishing contributions to non-gaussianities.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The hemispherical power asymmetry in the cosmic microwave background may be explained by the modulation of some primordial cosmological parameter, such as the sound speed of the inflaton. This modulation can be achieved by the perturbation of a light field during inflation. We numerically examine this mechanism and show it can be consistent with current observations. Further, this model predicts that the power asymmetry also exists in the temperature-polarization correlation and polarization autocorrelation with the same shape, and in primordial non-gaussianity of equilateral type with a particular shape. Therefore, our mechanism is observationally detectable in forthcoming experiments.

We study primordial gravitational waves (PGWs) in the Horava-Lifshitz (HL) theory of quantum gravity, in which high-order spatial derivative operators, including the ones violating parity, generically appear in order for the theory to be power-counting renormalizable and ultraviolet (UV) complete. Because of both parity violation and non-adiabatic evolution of the modes due to a modified dispersion relationship, a large polarization of PGWs becomes possible, and it could be well within the range of detection of the BB, TB and EB power spectra of the forthcoming cosmic microwave background (CMB) observations.

We investigate the local properties of WMAP Cold Spot (CS) by defining the local statistics: mean temperature, variance, skewness and kurtosis. We find that, compared with the \emph{coldest spots} in random Gaussian simulations, WMAP CS deviates from Gaussianity at $\sim 99%$ significant level. In the meanwhile, when compared with the spots at the same position in the simulated maps, the values of local variance and skewness around CS are all systematically larger in the scale of $R>5^{\circ}$, which implies that WMAP CS is a large-scale non-Gaussian structure, rather than a combination of some small structures. This is consistent with the finding that the non-Gaussianity of CS is totally encoded in the WMAP low multipoles. Furthermore, we find that the cosmic texture can excellently explain all the anomalies in these statistics.

Measurement of the pulsar timing residuals provides a direct way to detect relic gravitational waves at the frequency $f\sim 1/{\rm yr}$. In this paper, we investigate the constraints on the inflationary parameters, the tensor-to-scalar ratio $r$ and the tensor spectral index $n_t$, by the current and future Pulsar Timing Arrays (PTAs). We find that Five-hundred-meter Aperture Spherical radio Telescope (FAST) in China and the planned Square Kilometer Array (SKA) projects have the fairly strong abilities to test the phantom-like inflationary models. If $r=0.1$, FAST could give the constraint on the spectral index $n_t<0.56$, and SKA gives $n_t<0.32$. While an observation with the total time T=20 yr, the pulsar noise level $\sigma_w=30$ns and the monitored pulsar number $n=200$, could even constrain $n_t<0.07$. These are much tighter than those inferred from the current results of Parkers Pulsar Timing Array (PPTA), European Pulsar Timing Array (EPTA) and North American Nanohertz Observatory for Gravitational waves (NANOGrav). Especially, by studying the effects of various observational factors on the sensitivities of PTAs, we found that compared with $\sigma_w$ and $n$, the total observation time $T$ has the most significant effect.

In this paper, we study the anisotropy of cosmic acceleration by dividing the Union2 Type Ia supernova dataset into 12 subsets according to their positions in Galactic coordinate system. In each region, we derive the deceleration parameter $q_0$ as the diagnostic to quantify the anisotropy level in the corresponding direction, and construct $q_0$ anisotropic maps by combining these $q_0$ values. In addition to the monopole component, we find the significant dipole effect in the $q_0$-maps with the amplitude $A_1=0.466^{+0.255}_{-0.205}$, which deviates from zero at more than 2-$\sigma$ level. The direction of the best-fit dipole is ($\theta=108.8^{\circ}$, $\phi=187.0^{\circ}$) in Galactic system. Interesting enough, we find the direction of this dipole is nearly perpendicular to the CMB kinematic dipole, and the angle between them is $95.7^{\circ}$. The perpendicular relation is anomalous at the 1-in-10 level.

We compute the gravitational radiation of the supermassive binary black hole OJ287. By fitting the data of its recent seven outbursts, we obtain the orbital motion up to 3.5PN (Post-Newtonian) order, the energy and angular momentum fluxes and the waveform up to the 2PN. It is found that the 1PN term of the energy flux, has an opposite sign to the Newtonian flux, and the 2PN term has an opposite sign to the 1PN one. The same pattern is also found for the angular momentum fluxes. The total flux is reduced by 30.8% on average from the Newtonian flux for OJ287, resulting in a smaller orbital decay rate than that of 2.5PN calculation. Consisting with this, it is checked that, in the sequence of non-dissipative PN forces, each term of $i^{th}$ order (dissipative or non-dissipative) has a sign opposite to the $(i-1)^{th}$ order. The origin of this characteristic is traced to the appearance of the non-diagonal metric component $g_{0i}$ in PN approximation. This feature will have profound impact on estimation of gravitational waves from binary systems.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

The low multipole anomalies of the Cosmic Microwave Background has received much attention during the last few years. It is still not ascertained whether these anomalies are indeed primordial or the result of systematics or foregrounds. An example of a foreground, which could generate some non-Gaussian and statistically anisotropic features at low multipole range, is the very symmetric Kuiper Belt in the outer solar system. In this paper, expanding upon the methods presented by Maris et al. (2011), we investigate the contributions from the Kuiper Belt objects (KBO) to the WMAP ILC 7 map, whereby we can minimize the contrast in power between even and odd multipoles in the CMB, discussed discussed by Kim & Naselsky (2010). We submit our KBO de-correlated CMB signal to several tests, to analyze its validity, and find that incorporation of the KBO emission can decrease the quadrupole-octupole alignment and parity asymmetry problems, provided that the KBO signals has a non-cosmological dipole modulation, associated with the statistical anisotropy of the ILC 7 map. Additionally, we show that the amplitude of the dipole modulation, within a 2 sigma interval, is in agreement with the corresponding amplitudes, discussed by Lew (2008).

The contribution from the residuals of the foreground can have a significant impact on the temperature maps of the Cosmic Microwave Background (CMB). Mostly, the focus has been on the galactic plane, when foreground cleaning has taken place. However, in this paper, we will investigate the possible foreground contamination, from sources outside the galactic plane in the CMB maps. We will analyze the correlation between the Faraday rotation map and the CMB temperature map. The Faraday rotation map is dependent on the galactic magnetic field, as well as the thermal electron density, and both may contribute to the CMB temperature. We find that the standard deviation for the mean cross correlation deviate from that of simulations at the 99.9% level. Additionally, a comparison between the CMB temperature extrema and the extremum points of the Faraday rotation is also performed, showing a general overlap between the two. Also we find that the CMB Cold Spot is located at an area of strong negative cross correlation, meaning that it may be explained by a galactic origin. Further, we investigate nearby supernova remnants in the galaxy, traced by the galactic radio loops. These super nova remnants are located at high and low galactic latitude, and thus well outside the galactic plane. We find some correlation between the Faraday Rotation and the CMB temperature, at select radio loops. This indicate, that the galactic foregrounds may affect the CMB, at high galactic latitudes

Parity violation found in the Cosmic Microwave Background (CMB) radiation is a crucial clue for the non-standard cosmological model or the possible contamination of various foreground residuals and/or calibration of the CMB data sets. In this paper, we study the directional properties of the CMB parity asymmetry by excluding the $m=0$ modes in the definition of parity parameters. We find that the preferred directions of the parity parameters coincide with the CMB kinematic dipole, which implies that the CMB parity asymmetry may be connected with the possible contamination of the residual dipole component. We also find that such tendency is not only localized at $l=2,3$, but in the extended multipole ranges up to $l\sim 22$.

Testing the so-called consistency relations plays an important role for distinguishing the different classes of inflation models. In this paper, we investigate the possible testing based on the potential observations of the cosmic microwave background (CMB) radiation, including the planned CMBPol mission and the ideal CMB experiment where only the reduced cosmic weak lensing contamination for the B-mode polarization is considered. We find that for the canonical single-field inflation and the phantom inflation, the consistency relations are quite hard to be tested: the testing is possible only if $r>0.14$ for CMBPol mission, and $r>0.06$ for the ideal experiment. However, the situation could become much better for the general Lorentz invariant single-field inflation with large non-gaussian signal and the two-field inflation with strong correlation between the adiabatic and the isocurvature perturbations. We find that for these two classes of inflation the testing is possible if $r\gtrsim 10^{-2}$ or even smaller for both CMB experiments.

Analyzing the imprint of relic gravitational waves (RGWs) on the cosmic microwave background (CMB) power spectra provides a way to determine the signal of RGWs. In this Letter, we discuss a statistical bias, which could exist in the data analysis and has the tendency to overlook the RGWs. We also explain why this bias exists, and how to avoid it.

Recent pulsar timing observations by the Parkers Pulsar Timing Array and European Pulsar Timing Array teams obtained the constraint on the relic gravitational waves at the frequency $f_*=1/{\rm yr}$, which provides the opportunity to constrain $H_*$, the Hubble parameter when these waves crossed the horizon during inflation. In this paper, we investigate this constraint by considering the general scenario for the early Universe: we assume that the effective (average) equation-of-state $w$ before the big bang nucleosynthesis stage is a free parameter. In the standard hot big-bang scenario with $w=1/3$, we find that the current PPTA result follows a bound $H_*\leq 1.15\times10^{-1}\mpl$, and the EPTA result follows $H_*\leq 6.92\times10^{-2}\mpl$. We also find that these bounds become much tighter in the nonstandard scenarios with $w>1/3$. When $w=1$, the bounds become $H_*\leq5.89\times10^{-3}\mpl$ for the current PPTA and $H_*\leq3.39\times10^{-3}\mpl$ for the current EPTA. In contrast, in the nonstandard scenario with $w=0$, the bound becomes $H_*\leq7.76\mpl$ for the current PPTA.

1156+295 is a flat-spectrum quasar which is loud at radio and gamma-ray. Previous observations of the source revealed a radio morphology on pc to kpc scales consistent with a helical jet model. In our present research, this source was observed with the VLBA at 86, 43 and 15 GHz on four epochs from 10 May 2003 to 13 March 2005 aiming at studying the structure of the innermost jet in order to understand the relation between the helical structure and the astrophysical processes in the central engine. A core-jet structure with six jet components is identified. The apparent transverse velocities of the six jet components derived from proper motion measurements are in the range between 3.6 c and 11.6 c. The overall jet shape shows oscillatory morphology with multiple curvatures on pc scales which might be indicative of a helical pattern. Models of helical jet are discussed on the basis of both Kelvin-Helmholtz (K-H) instability and jet precession. The K-H instability model shows better agreement with the observed data. The overall radio structure on the scale from sub-pc to kpc appears to be fitted with a hydrodynamic model with the fundamental helical mode in Kelvin-Helmholtz (K-H) instability. This helical mode with an initial characteristic wavelength of 0.2 pc is excited at the base of the jet on the scale of 0.005 pc (or 1000R_s, the typical size of the broad line region for a super massive black hole of $4.3\times10^8M_{\odot}$). A presessing jet model can also fit the observed jet structure on the scale between 10 pc and 300 pc. However, additional astrophysical processes may be required for the presessing jet model in order to explain the bendings on the inner jet structure (1 to 10 pc) and re-collimation of the large scale jet outflow (>300 pc).

Detection of relic gravitational waves, through their imprint in the cosmic microwave background radiation, is one of the most important tasks for the planned CMBPol mission. In the simplest viable theoretical models the gravitational wave background is characterized by two parameters, the tensor-to-scalar ratio $r$ and the tensor spectral index $n_t$. In this paper, we analyze the potential joint constraints on these two parameters, $r$ and $n_t$, using the potential observations of the CMBPol mission, which is expected to detect the relic gravitational waves if $r\gtrsim0.001$. The influence of the contaminations, including cosmic weak lensing, various foreground emissions, and systematical errors, is discussed.

A design study is currently in progress for a third generation gravitational-wave (GW) detector called Einstein Telescope (ET). An important kind of source for ET will be the inspiral and merger of binary neutron stars (BNS) up to $z \sim 2$. If BNS mergers are the progenitors of short-hard $\gamma$-ray bursts, then some fraction of them will be seen both electromagnetically and through GW, so that the luminosity distance and the redshift of the source can be determined separately. An important property of these `standard sirens' is that they are \emph{self-calibrating}: the luminosity distance can be inferred directly from the GW signal, with no need for a cosmic distance ladder. Thus, standard sirens will provide a powerful independent check of the $\Lambda$CDM model. In previous work, estimates were made of how well ET would be able to measure a subset of the cosmological parameters (such as the dark energy parameter $w_0$) it will have access to, assuming that the others had been determined to great accuracy by alternative means. Here we perform a more careful analysis by explicitly using the potential Planck CMB data as prior information for these other parameters. We find that ET will be able to constrain $w_0$ and $w_a$ with accuracies $\Delta w_0 = 0.099$ and $\Delta w_a = 0.302$, respectively. These results are compared with projected accuracies for the JDEM Baryon Acoustic Oscillations project and the SNAP Type Ia supernovae observations.

The forthcoming release of data from the Planck mission, and possibly from the next round of Wilkinson Microwave Anisotropy Probe (WMAP) observations, make it necessary to revise the evaluations of relic gravitational waves in the existing data and, at the same time, to refine the assumptions and data analysis techniques in preparation for the arrival of new data. We reconsider with the help of the commonly used CosmoMC numerical package the previously found indications of relic gravitational waves in the 7-year (WMAP7) data. The CosmoMC approach reduces the confidence of these indications from approximately 2$\sigma$ level to approximately 1$\sigma$ level, but the indications do not disappear altogether. We critically analyze the assumptions that are currently used in the Cosmic Microwave Background (CMB) data analyzes and outline the strategy that should help avoid the oversight of relic gravitational waves in the future CMB data. In particular, it is important to keep away from the unwarranted assumptions about density perturbations. The prospects of confident detection of relic gravitational waves by the Planck satellite have worsened, but they are still good. It appears that more effort will be required in order to mitigate the foreground contamination.

We investigate the observational signatures of three models of the early Universe in the $B$-mode polarization of the Cosmic Microwave Background (CMB) radiation. In addition to the standard single field inflationary model, we also consider the constraints obtainable on the loop quantum cosmology model (from Loop Quantum Gravity) and on cosmic strings, expected to be copiously produced during the latter stages of Brane inflation. We first examine the observational features of the three models, and then use current $B$-mode polarization data from the BICEP and QUaD experiments to constrain their parameters. We also examine the detectability of the primordial $B$-mode signal predicted by these models and forecast the parameter constraints achievable with future CMB polarization experiments. We find that: (a) these three models of the early Universe predict different features in the CMB B-mode polarization power spectrum, which are potentially distinguishable from the CMB experiments; (b) since $B$-mode polarization measurements are mostly unaffected by parameter degeneracies, they provide the cleanest probe of these early Universe models; (c) using the BICEP and QUaD data we obtain the following parameter constraints: $r=0.02^{+0.31}_{-0.26}$ ($1\sigma$ for the tensor-to-scalar ratio in the single field inflationary model); $m < 1.36\times 10^{-8} \text{M}_{\text{pl}}$ and $k_{*} < 2.43 \times 10^{-4} \text{Mpc}^{-1}$ ($1\sigma$ for the mass and scale parameters in the loop quantum cosmology model); and $G\mu < 5.77 \times 10^{-7}$ ($1\sigma$ for the cosmic string tension); (d) future CMB observations (both satellite missions and forthcoming sub-orbital experiments) will provide much more rigorous tests of these early Universe models.

The new release of data from Wilkinson Microwave Anisotropy Probe improves the observational status of relic gravitational waves. The 7-year results enhance the indications of relic gravitational waves in the existing data and change to the better the prospects of confident detection of relic gravitational waves by the currently operating Planck satellite. We apply to WMAP7 data the same methods of analysis that we used earlier [W. Zhao, D. Baskaran, and L.P. Grishchuk, Phys. Rev. D 80, 083005 (2009)] with WMAP5 data. We also revised by the same methods our previous analysis of WMAP3 data. It follows from the examination of consecutive WMAP data releases that the maximum likelihood value of the quadrupole ratio $R$, which characterizes the amount of relic gravitational waves, increases up to $R=0.264$, and the interval separating this value from the point $R=0$ (the hypothesis of no gravitational waves) increases up to a $2\sigma$ level. The primordial spectra of density perturbations and gravitational waves remain blue in the relevant interval of wavelengths, but the spectral indices increase up to $n_s =1.111$ and $n_t=0.111$. Assuming that the maximum likelihood estimates of the perturbation parameters that we found from WMAP7 data are the true values of the parameters, we find that the signal-to-noise ratio $S/N$ for the detection of relic gravitational waves by the Planck experiment increases up to $S/N=4.04$, even under pessimistic assumptions with regard to residual foreground contamination and instrumental noises. We comment on theoretical frameworks that, in the case of success, will be accepted or decisively rejected by the Planck observations.

Detection of magnetic-type ($B$-type) polarization in the Cosmic Microwave Background (CMB) radiation plays a crucial role in probing the relic gravitational wave (RGW) background. In this paper, we propose a new method to deconstruct a polarization map on an incomplete sky in real space into purely electric and magnetic polarization type maps, ${\mathcal{E}}(\hat{\gamma})$ and ${\mathcal{B}}(\hat{\gamma})$, respectively. The main properties of our approach are as follows: Firstly, the fields ${\mathcal{E}}(\hat{\gamma})$ and ${\mathcal{B}}(\hat{\gamma})$ are constructed in real space with a minimal loss of information. This loss of information arises due to the removal of a narrow edge of the constructed map in order to remove various numerical errors, including those arising from finite pixel size. Secondly, this method is fast and can be efficiently applied to high resolution maps due to the use of the fast spherical harmonics transformation. Thirdly, the constructed fields, ${\mathcal{E}}(\hat{\gamma})$ and ${\mathcal{B}}(\hat{\gamma})$, are scalar fields. For this reason various techniques developed to deal with temperature anisotropy maps can be directly applied to analyze these fields. As a concrete example, we construct and analyze an unbiased estimator for the power spectrum of the $B$-mode of polarization $C_{\ell}^{BB}$. Basing our results on the performance of this estimator, we discuss the RGW detection ability of two future ground-based CMB experiments, QUIET and POLARBEAR.

In this paper, we improve the calculation of the relic gravitational waves (RGW) in two aspects: First, we investigate the transfer function after considering the redshift-suppression effect, the accelerating expansion effect, the damping effect of free-streaming relativistic particles, and the damping effect of cosmic phase transition, and give a simple approximate analytic expression, which clearly illustrates the dependent relations on the cosmological parameters. Second, we develop a numerical method to calculate the primordial power spectrum of RGW at a very wide frequency range, where the observed constraints on $n_s$ (the scalar spectral index) and $P_S(k_0)$ (the amplitude of primordial scalar spectrum) and the Hamilton-Jacobi equation are used. This method is applied to two kinds of inflationary models, which all satisfy the current constraints on $n_s$, $\alpha$ (the running of $n_s$) and $r$ (the tensor-scalar ratio). We plot them in the $r-\Omega_g$ diagram, where $\Omega_g$ is the strength of RGW, and study their detection by the CMB experiments and laser interferometers.

In this paper, we investigate the quintessence models with an oscillating equation of state (EoS) and their potentials. From the constructed potentials, which have the EoS of $\omega_{\phi}=\omega_0+\omega_1\sin z$, we find they are all the oscillating functions of the field $\phi$, and the oscillating amplitudes are decreasing (or increasing) with $\phi$. From the evolutive equation of the field $\phi$, we find this is caused by the expansion of the universe. This also makes that it is very difficult to build a model whose EoS oscillates forever. However one can build a model with EoS oscillating for a period. Then we discuss three quintessence models, which are the combinations of the invert power law functions and the oscillating functions of the field $\phi$. We find they all follow the oscillating EoS.

The accelerating expansion of the Universe at the present stage is a process that will change the spectrum of relic gravitational waves. Here we present a numerical calculation for the power spectrum of relic gravitational waves in the accelerating Universe. The results show that although the overall features of the power spectrum is similar to that in the non-accelerating models, the amplitude is smaller by an order of $10^{-1}$. We also find that the spectrum is very sensitive to the index $\beta$ of the inflationary expansion with the scale factor $a(\tau) \propto |\tau|^{1+\beta}$. With increase of $\beta$, the resulting spectrum is tilted to be flatter with more power on high frequencies, and the sensitivity of the second science run of the LIGO detectors puts a restriction on the parameter $\beta < -1.8$. The influence of reheating following the inflation has been examined.

This paper has been withdrawn.

This paper has been withdrawn.

This is a summary of presentations delivered at the OC1 parallel session "Primordial Gravitational Waves and the CMB" of the 12th Marcel Grossmann meeting in Paris, July 2009. The reports and discussions demonstrated significant progress that was achieved in theory and observations. It appears that the existing data provide some indications of the presence of gravitational wave contribution to the CMB anisotropies, while ongoing and planned observational efforts are likely to convert these indications into more confident statements about the actual detection.

The relic gravitational waves (gw) are the cleanest probe of the violent times in the very early history of the Universe. They are expected to leave signatures in the observed cosmic microwave background anisotropies. We significantly improved our previous analysis [1] of the 5-year WMAP $TT$ and $TE$ data at lower multipoles $\ell$. This more general analysis returned essentially the same maximum likelihood (ML) result (unfortunately, surrounded by large remaining uncertainties): the relic gw are present and they are responsible for approximately 20% of the temperature quadrupole. We identify and discuss the reasons by which the contribution of gw can be overlooked in a data analysis. One of the reasons is a misleading reliance on data from very high multipoles $\ell$, another - a too narrow understanding of the problem as the search for $B$-modes of polarization, rather than the detection of relic gw with the help of all correlation functions. Our analysis of WMAP5 data has led to the identification of a whole family of models characterized by relatively high values of the likelihood function. Using the Fisher matrix formalism we formulated forecasts for {\it Planck} mission in the context of this family of models. We explore in details various `optimistic', `pessimistic' and `dream case' scenarios. We show that in some circumstances the $B$-mode detection may be very inconclusive, at the level of signal-to-noise ratio $S/N =1.75$, whereas a smarter data analysis can reveal the same gw signal at $S/N= 6.48$. The final result is encouraging. Even under unfavourable conditions in terms of instrumental noises and foregrounds, the relic gw, if they are characterized by the ML parameters that we found from WMAP5 data, will be detected by {\it Planck} at the level $S/N = 3.65$.

We present results on the compact steep-spectrum quasar 3C 48 from observations with the VLBA, MERLIN and EVN at multiple radio frequencies. In the 1.5-GHz VLBI images, the radio jet is characterized by a series of bright knots. The active nucleus is embedded in the southernmost VLBI component A, which is further resolved into two sub-components A1 and A2 at 4.8 and 8.3 GHz. A1 shows a flat spectrum and A2 shows a steep spectrum. The most strongly polarized VLBI components are located at component C $\sim$0.25 arcsec north of the core. The polarization angles at C show gradual changes across the jet width at all observed frequencies, indicative of a gradient in the emission-weighted intrinsic polarization angle across the jet and possibly a systematic gradient in the rotation measure; moreover, the percentage of polarization increases near the curvature at C, likely consistent with the presence of a local jet-ISM interaction and/or changing magnetic-field directions. The hot spot B shows a higher rotation measure, and has no detected proper motion. These facts provide some evidence for a stationary shock in the vicinity of B. Comparison of the present VLBI observations with those made 8.43 years ago suggests a proper motion of $\beta_{app}=3.7\pm0.4 c$ for A2 to the north. The apparent superluminal motion suggests that the relativistic jet plasma moves at a velocity of $\gtrsim0.96 c$ if the jet is viewed at an inclination angle less than $20\degr$. A simple precessing jet model and a hydrodynamical isothermal jet model with helical-mode Kelvin-Helmholtz instabilities are used to fit the oscillatory jet trajectory of 3C 48 defined by the bright knots.

We review the quantum Yang-Mills condensate (YMC) dark energy models. As the effective Yang-Mills Lagrangian is completely determined by the quantum field theory, there is no adjustable parameter in the model except the energy scale. In this model, the equation-of-state (EOS) of the YMC dark energy, $w_y > -1$ and $w_y < -1$, can both be naturally realized. By studying the evolution of various components in the model, we find that, in the early stage of the universe, dark energy tracked the evolution of the radiation, i.e. $w_y \to 1/3$. However, in the late stage, $w_y$ naturally runs to the critical state with $w_y = -1$, and the universe transits from matter-dominated into dark energy dominated stage only at recently $z \sim 0.3$. These characters are independent of the choice of the initial condition, and the cosmic coincidence problem is avoided in the models. We also find that, if the possible interaction between YMC and dust matter is considered, the late time attractor solution may exist. In this case, the EOS of YMC must evolve from $w_y>0$ into $w_y < -1$, which is slightly suggested by the observations. At the same time, the total EOS in the attractor solution is $w_{tot} = -1$, the universe being the de Sitter expansion in the late stage, and the cosmic big rip is naturally avoided. These features are all independent of the interacting forms.

A thermal gravitational wave background can be produced in the early Universe if a radiation dominated epoch precedes the usual inflationary stage. This background provides a unique way to study the initial state of the Universe. We discuss the imprint of this thermal spectra of gravitons on the cosmic microwave background (CMB) power spectra, and its possible detection by CMB observations. Assuming the inflationary stage is a pure de Sitter expansion we find that, if the number of e-folds of inflation is smaller than 65, the signal of this thermal spectrum can be detected by the observations of Planck and PolarBear experiments, or the planned EPIC experiments. This bound can be even looser if inflation-like stage is the sub-exponential.

We investigate the attractor solution in the coupled Yang-Mills field dark energy models with the general interaction term, and obtain the constraint equations for the interaction if the attractor solution exists. The research also shows that, if the attractor solution exists, the equation-of-state of the dark energy must evolve from $w_y>0$ to $w_y\le-1$, which is slightly suggested by the observation. At the same time, the total equation-of-state in the attractor solution is $w_{tot}=-1$, the universe is a de Sitter expansion, and the cosmic big rip is naturally avoided. These features are all independent of the interacting forms.

It has been suggested that the Yang-Mills (YM) field can be a kind of candidate for the inflationary field at high energy scales or the dark energy at very low energy scales, which can naturally give the equation of state $-1<\omega<0$ or $\omega<-1$. We discuss the zero order and first order Einstein equations and YM field kinetic energy equations of the free YM field models. From the zero order equations, we find that $\omega+1\propto a^{-2}$, from which it follows that the equation of state of YM field always goes to -1, independent of the initial conditions. By solving the first order Einstein equations and YM field equations, we find that in the YM field inflationary models, the scale-invariant primordial perturbation power spectrum cannot be generated. Therefore, only this kind of YM field is not enough to account for inflationary sources. However, as a kind of candidate of dark energy, the YM field has the `sound speed' $c_s^2=-1/3<0$, which makes the perturbation $\phi$ have a damping behavior at the large scale. This provides a way to distinguish the YM field dark energy models from other kinds of models.

We compare the detection abilities for the relic gravitational waves by two kinds of forthcoming cosmic microwave background radiation (CMB) experiments, space-based Planck satellite and the various ground-based experiments. Comparing with the ground-based experiments, Planck satellite can observe all the CMB power spectra in all the multipole range, but having much larger instrumental noises. We find that, for the uncertainty of the tensor-to-scalar ratio $\Delta r$, PolarBear (II) as a typical ground-based experiment can give much smaller value than Planck satellite. However, for the uncertainty of the spectral index $\Delta n_t$, Planck can give the similar result with PolarBear (II). If combining these two experiments, the value of $\Delta n_t$ can be reduced by a factor 2. For the model with $r=0.1$, the constraint $\Delta n_t=0.10$ is expected to be achieved, which provides an excellent opportunity to study the physics in the very early universe. We also find the observation in the largest scale ($\ell<20$) is very important for constraining the spectral index $n_t$. So it is necessary to combine the observations of the future space-based and ground-based CMB experiments to determine the relic gravitational waves.

We discuss a method to constrain the fraction density $f$ of the relativistic gas in the radiation dominant stage, by their impacts on a relic gravitational waves and the cosmic microwave background (CMB) $B$ polarization power spectrum. We find that the uncertainty of $f$ strongly depends on the noise power spectra of the CMB experiments and the amplitude of the gravitational waves. Taking into account of the CMBPol instrumental noises, an uncertainty $\Delta f=0.046$ is obtained for the model with tensor-to scalar ratio $r=0.1$. For an ideal experiment with only the reduced cosmic lensing as the contamination of $B$-polarization, $\Delta f=0.008$ is obtained for the model with $r=0.1$. So the precise observation of the CMB $B$-polarization provides a great opportunity to study the relativistic components in the early Universe.

In this paper, we discuss the constraint on the relic gravitational waves by both temperature and polarization anisotropies power spectra of cosmic microwave background radiation. Taking into account the instrumental noises of Planck satellite, we calculate the signal-to-noise ratio $S/N$ by the simulation and the analytic approximation methods. We find that, comparing with the $BB$ channel, the value of $S/N$ is much improved in the case where all the power spectra, $TT$, $TE$, $EE$ and $BB$, are considered. If the noise power spectra of Planck satellite increase for some reasons, the value of $S/N$ in $BB$ channel is much reduced. However, in the latter case where all the power spectra of cosmic microwave background radiation are considered, the value of $S/N$ is less influenced. We also find that the free parameters $A_s$, $n_s$ and $n_t$ have little influence on the value of $S/N$ in both cases.

The prospect of detecting relic gravitational waves (RGWs), through their imprint in the cosmic microwave background radiation, provides an excellent opportunity to study the very early Universe. In simplest viable theoretical models the RGW background is characterized by two parameters, the tensor-to-scalar ratio $r$ and the tensor spectral index $n_t$. In this paper, we analyze the potential joint constraints on these two parameters, $r$ and $n_t$, using the data from the upcoming cosmic microwave background radiation experiments. Introducing the notion of the best pivot multipole $\ell_t^*$, we find that at this pivot multipole the parameters $r$ and $n_t$ are uncorrelated, and have the smallest variances. We derive the analytical formulae for the best pivot multipole number $\ell_t^*$, and the variances of the parameters $r$ and $n_t$. We verify these analytical calculations using numerical simulation methods, and find agreement to within 20%. The analytical results provides a simple way to estimate the detection ability for the relic gravitational waves by the future observations of the cosmic microwave background radiation.

The detection of primordial gravitational waves is one of the biggest challenges of the present time. The existing (Wilkinson Microwave Anisotropy Probe) observations are helpful on the road to this goal, and the forthcoming experiments (Planck) are likely to complete this mission. We show that the 5-year Wilkinson Microwave Anisotropy Probe $TE$ data contains a hint of the presence of gravitational wave contribution. In terms of the parameter $R$, which gives the ratio of contributions from gravitational waves and density perturbations to the temperature quadrupole, the best-fit model produced $R=0.24$. Because of large residual noises, the uncertainty of this determination is still large, and it easily includes the R=0 hypothesis. However, the uncertainty will be strongly reduced in the forthcoming observations which are more sensitive. We numerically simulated the Planck data and concluded that the relic gravitational waves with $R=0.24$ will be present at a better than 3$\sigma$ level in the $TE$ observational channel, and at a better than 2$\sigma$ level in the `realistic' $BB$ channel. The balloon-borne and ground-based observations may provide a healthy competition to Planck in some parts of the lower-$\ell$ spectrum.

An exact analytic calculation is presented for the spectrum of relic gravitational waves in the scenario of accelerating Universe $\Omega_{\Lambda}+\Omega_m = 1$. The spectrum formula contains explicitly the parameters of acceleration, inflation, reheating, and the (tensor/scalar) ratio, so that it can be employed for a variety of cosmological models. We find that the spectrum depends on the behavior of the present accelerating expansion. The amplitude of gravitational waves for the model $\Omega_{\Lambda}=0.65$ is about $\sim 50 %$ greater than that of the model $\Omega_{\Lambda}=0.7$, an effect accessible to the designed sensitivities of LIGO and LISA. The spectrum sensitively depends on the inflationary models with $a(\tau) \propto |\tau|^{1+\beta}$, and a larger $\beta$ yields a flatter spectrum, producing more power. The current LIGO results rule out the inflationary models of $\beta \geq -1.8$. The LIGO with its design sensitivity and the LISA will also be able to test the model of $\beta=-1.9$. We also examine the constraints on the spectral energy density of relic gravitational waves. Both the LIGO bound and the nucleosynthesis bound point out that the model $\beta=-1.8$ is ruled out, but the model $\beta=-2.0$ is still alive. The exact analytic results also confirm the approximate spectrum and the numerical one in our previous work.

The evolution of the electric and magnetic components in an effective Yang-Mills condensate dark energy model is investigated. If the electric field is dominant, the magnetic component disappears with the expansion of the Universe. The total YM condensate tracks the radiation in the earlier Universe, and later it becomes $w_y\sim-1$ thus is similar to the cosmological constant. So the cosmic coincidence problem can be avoided in this model. However, if the magnetic field is dominant, $w_y>1/3$ holds for all time, suggesting that it cannot be a candidate for the dark energy in this case.

In this paper, we study the possibility of building Yang-Mills(YM) field dark energy models with equation of state (EoS) crossing -1, and find that it can not be realized by the single YM field models, no matter what kind of lagrangian or initial condition. But the states of $-1<\omega<0$ and $\omega<-1$ all can be naturally got in this kind of models. The former is like a quintessence field, and the latter is like a phantom field. This makes that one can build a model with two YM fields, in which one with the initial state of $-1<\omega<0$, and the other with $\omega<-1$. We give an example model of this kind, and find that its EoS is larger than -1 in the past and less than -1 at the present time. We also find that this change must be from $\omega>-1$ to $<-1$, and it will go to the critical state of $\omega=-1$ with the expansion of the Universe, which character is same with the single YM field models, and the Big Rip is naturally avoided.

We study the statefinder parameters in the Yang-Mills condensate dark energy models, and find that the evolving trajectories of these models are different from those of other dark energy models. We also define two eigenfunctions of the Yang-Mills condensate dark energy models. The values of these eigenfunctions are quite close to zero if the equation-of-state of the Yang-Mills condensate is not far from -1, which can be used to simply differentiate between the Yang-Mills condensate models and other dark energy models.

By Polnarev's method we analytically calculate the polarization spectra of the cosmic microwave background radiation (CMB) generated by cosmic relic gravitational waves (RGW). In this analytic approach the physics involved in this generating process is more transparent. Consequently, the effects due to various elements of physics can be isolated easily. In solving the equation for evolution of RGW in the expanding universe, both the sudden transition and the WKB approximation for the scale factor during the radiation-matter transition have been taken. To describe more precisely the decoupling process, we have introduced an analytic expression for the visibility function, consisting of two pieces of half-Gaussian curves. We also include the damping on polarizations due to the photon diffusion up to the second order of the tight coupling. Analytic polarization spectra $C_l^{XX}$ have been obtained with the following several improvements over the previous results. 1. The approximate analytic result is quite close to the numerical one evaluated from the cmbfast code, especially, for the first three peaks of the spectrum that are observable (see the text).

As the strong evidence for inflation, the relic gravitational waves (RGW) have been extensively studied. Although, it has not been detected, yet some constraints have been achieved by the observations. Future experiments for the RGW detection are mainly two kinds: the CMB experiments and the laser interferometers. In this paper, we study these current constraints and the detective abilities of future experiments. We calculate the strength of RGW $\Omega_g(k)$ in two methods: the analytic method and the numerical method by solving the inflationary flow equations. By the first method we obtain a bound $\Omega_g<3.89\times10^{-16}$ at $\nu=0.1$Hz, where we have used the current constraints on the scalar spectral index, the tensor-scalar ratio, furthermore, we have taken into account of the redshift-suppression effect, the accelerating expansion effect, the neutrino damping effect on the RGW. But the analytic expression of $\Omega_g(k)$ depends on the specific inflationary models and not well applies for the waves with very high frequencies. The numerical method is more precise for the waves with high frequencies and applies to any single-field inflationary model. It gives a bound $\Omega_g<8.62\times10^{-14}$, which is independent of the inflationary parameters, and applies to any single-field slow-roll inflationary model. After considering the current constraints on the inflationary parameters, this bound becomes $\Omega_g<2\times10^{-17}$ (see text).

In this paper, we investigate a kind of special quintom model, which is made of a quintessence field $\phi_1$ and a phantom field $\phi_2$, and the potential function has the form of $V(\phi_1^2-\phi_2^2)$. This kind of quintom fields can be separated into two kinds: the hessence model, which has the state of $\phi_1^2>\phi_2^2$, and the hantom model with the state $\phi_1^2<\phi_2^2$. We discuss the evolution of these models in the $\omega$-$\omega'$plane ($\omega$ is the state equation of the dark energy, and $\omega'$ is its time derivative in unites of Hubble time), and find that according to $\omega>-1$ or $<-1$, and the potential of the quintom being climbed up or rolled down, the $\omega$-$\omega'$ plane can be divided into four parts. The late time attractor solution, if existing, is always quintessence-like or $\Lambda$-like for hessence field, so the Big Rip doesn't exist. But for hantom field, its late time attractor solution can be phantom-like or $\Lambda$-like, and sometimes, the Big Rip is unavoidable. Then we consider two special cases: one is the hessence field with an exponential potential, and the other is with a power law potential. We investigate their evolution in the $\omega$-$\omega'$ plane. We also develop a theoretical method of constructing the hessence potential function directly from the effective equation of state function $\omega(z)$. We apply our method to five kinds of parametrizations of equation of state parameter, where $\omega$ crossing -1 can exist, and find they all can be realized. At last, we discuss the evolution of the perturbations of the quintom field, and find the perturbations of the quintom $\delta_Q$ and the metric $\Phi$ are all finite even if at the state of $\omega=-1$ and $\omega'\neq0$.

The coincidence problem is studied in the effective Yang-Mills condensate dark energy model. As the effective YM Lagrangian is completely determined by quantum field theory, there is no adjustable parameter in this model except the energy scale, and the cosmic evolution only depends on the initial conditions. For generic initial conditions with the YM condensate subdominant to the radiation and matter, the model always has a tracking solution, the Universe transits from matter-dominated into the dark energy dominated stage only recently $z\sim 0. 3$, and evolve to the present state with $\Omega_{y}\sim 0.73$ and $\Omega_m\sim 0.27$.

We discuss the evolution of holographic hessence model, which satisfies the holographic principle and can naturally realizes the equation of state crossing -1. By discussing the evolution of the models in the $w w'$ plane, we find that, if $c\geq1$, $w_{he}\geq-1$ and $\dot{V}<0$ keep for all time, which are quintessence-like. However, if $c<-1$, which mildly favors the current observations, $w_{he}$ evolves from $w_{he}>-1$ to $w_{he}<-1$, and the potential is a nonmonotonic function. In the earlier time, the potential must be rolled down, and then be climbed up. Considered the current constraint on the parameter $c$, we reconstruct the potential of the holographic hessence model.

In this paper we briefly present our works on the relic gravitational waves (RGW) and the CMB polarization in the accelerating universe. The spectrum of RGW has been obtained, showing the influence of the dark energy. Compared with those from non-accelerating models, the shape of the spectrum is approximately similar, nevertheless, the amplitude of RGW now acquires a suppressing factor of the ratio of matter over dark energy $\propto \Omega_m/\Omega_{\Lambda}\sim 0.4$ over almost the whole range of frequencies. The RGW spectrum is then used as the source to calculate the spectra of CMB polarization. By a two half Gaussian function as an approximation to the visibility function during the photon decoupling, both the "electric" and "magnetic" spectra have been analytically derived, which are quite close to the numerical ones. Several physical elements that affect the spectra have been examined, such as the decoupling process, the inflation, the dark energy, the baryons, etc.

In the inflation models, the relic gravitational waves (RGW) generated in the inflation stage, and evolved in the Universe until now. In the different cosmological evolution models, one can get different gravitational waves power spectrum. In this paper, we give a simple formula to estimate this spectrum in a general cosmological model. From this formula, one can easily find the relation between this power spectrum and the cosmological evolution models. The spectrum includes all the information about the evolution of the scale factor $a$ from inflation to now. So RGW is a more clear fossil, which records the cosmological evolution information, and be a useful complementarity for another fossil$-$Cosmic Microwave Background Radiation (CMB).

From the damping of the Cosmic Microwave Background Radiation (CMB) anisotropy power spectrum at large scale and the recent accelerating expansion of the Universe, we find that, there may be a largest scale which we can detect in the Universe. From this, we can get the inflation parameters as spectrum index $n_s$, e-fold $N$, Hubble parameter $H$, the ratio of tensor and scalar $r$, the lasting time of reheating time $\alpha$ for special inflation models. We do them in three inflation models, and find that all the results fit very well with the observations and the inflation theory.

Observations have recently indicated that the Universe at the present stage is in an accelerating expansion, a process that has great implications. We evaluate the spectrum of relic gravitational waves in the current accelerating Universe and find that there are new features appearing in the resulting spectrum as compared to the decelerating models. In the low frequency range the peak of spectrum is now located at a frequency $\nu_E \simeq (\frac{\Omega_m}{\Omega_{\Lambda}})^{1/3} \nu_H$, where $\nu_H$ is the Hubble frequency, and there appears a new segment of spectrum between $\nu_E$ and $\nu_H$. In all other intervals of frequencies $\geq \nu_H$, the spectral amplitude acquires an extra factor $\frac{\Omega_m}{\Omega_{\Lambda}}$, due to the current acceleration, otherwise the shape of spectrum is similar to that in the decelerating models. The recent WMAP result of CMB anisotropies is used to normalize the amplitude for gravitational waves. The slope of the power spectrum depends sensitively on the scale factor $a(\tau) \propto |\tau|^{1+\beta}$ during the inflationary stage with $\beta=-2$ for the exact de Sitter space. With the increasing of $\beta$, the resulting spectrum is tilted to be flatter with more power on high frequencies, and the sensitivity of the second science run of the LIGO detectors puts a restriction on the parameter $\beta \leq -1.8$. We also give a numerical solution which confirms these features.